Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Molecular Dynamic simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41Linked to original sources
PubMed2001

Molecular dynamics simulations of pressure effects on hydrophobic interactions.

We report results on the pressure effects on hydrophobic interactions obtained from molecular dynamics simulations of aqueous solutions of methanes in water. A wide range of pressures that is relevant to pressure denaturation of proteins is investigated. The characteristic features of water-mediated interactions between hydrophobic solutes are found to be pressure-dependent. In particular, with increasing pressure we find that (1) the solvent-separated configurations in the solute-solute potential of mean force (PMF) are stabilized with respect to the contact configurations; (2) the desolvation barrier increases monotonically with respect to both contact and solvent-separated configurations; (3) the locations of the minima and the barrier move toward shorter separations; and (4) pressure effects are considerably amplified for larger hydrophobic solutes. Together, these observations lend strong support to the picture of the pressure denaturation process proposed previously by Hummer et al. (Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 1552): with increasing pressure, the transfer of water into protein interior becomes key to the pressure denaturation process, leading to the dissociation of close hydrophobic contacts and subsequent swelling of the hydrophobic protein interior through insertions of water molecules. The pressure dependence of the PMF between larger hydrophobic solutes shows that pressure effects on the interaction between hydrophobic amino acids may be considerably amplified compared to those on the methane-methane PMF.

Computer Simulation↗
PubMed2005

Molecular dynamics simulations on the effects of diameter and chirality on hydrogen adsorption in single walled carbon nanotubes.

We present systematic molecular dynamics simulation studies of hydrogen storage in single walled carbon nanotubes of various diameters and chiralities using a recently developed curvature-dependent force field. Our main objective is to address the following fundamental issues: 1. For a given H2 loading and nanotube type, what is the H2 distribution in the nanotube bundle? 2. For a given nanotube type, what is the maximal loading (H2 coverage)? 3. What is the diameter range and chirality for which H2 adsorption is most energetically favorable? Our simulation results suggest strong dependence of H2 adsorption energies on the nanotube diameter but less dependence on the chirality. Substantial lattice expansion upon H2 adsorption was found. The average adsorption energy increases with the lowering of nanotube diameter (higher curvature) and decreases with higher H2 loading. The calculated H2 vibrational power spectra and radial distribution functions indicate a strong attractive interaction between H2 and nanotube walls. The calculated diffusion coefficients are much higher than what has been reported for H2 in microporous materials such as zeolites, indicating that diffusivity does not present a problem for hydrogen storage in carbon nanotubes.

Journal Article↗
PubMed2006

Molecular dynamics simulation of laser desorption of a fragment of protein kinase A from two MALDI matrices.

We have carried out molecular dynamics simulations to study the desorption of a dephosphorylated fragment of protein kinase A from two matrices, sinapic acid (SA) and 2,5-dihydroxybenzoic acid (DHB), after laser excitation. We have examined the results as a function of the laser fluence and of the burial depth of the guest peptide in the matrices. In most cases, we found that the energy transferred from the matrix to the guest peptide was not sufficiently large to fragment the peptide. Exceptions occurred when the peptide was more buried. This finding suggested that protein analytes might be less likely to break into smaller fragments if they were placed closer to the surface of the matrix. We have also examined how likely the guest peptide could form small clusters with the matrix molecules and found that the results depended on the degree of burial of the peptide, on the laser fluence, and on which matrix was used. Generally, stable clusters were more likely to be formed for guest peptides that were more buried, at a lower laser fluence, and in the SA rather than the DHB matrix. In addition, we found that the DHB matrix was broken down more easily by the laser than the SA matrix.

Adsorption↗
PubMed2005

Molecular dynamics simulations of optical conductivity of dense plasmas.

The optical conductivity sigma (omega) for dense Coulomb systems is investigated using molecular dynamics simulations on the basis of pseudopotentials to mimic quantum effects. Starting from linear response theory, the response in the long-wavelength limit k=0 can be expressed by different types of autocorrelation functions (ACF's) such as the current ACF, the force ACF, or the charge density ACF. Consistent simulation data for transverse as well as longitudinal ACF's are shown which are based on calculations with high numerical accuracy. Results are compared with perturbation expansions which are restricted to small values of the plasma parameter. The relevance with respect to a quantum Coulomb plasma is discussed. Finally, results are presented showing a consistent description of these model plasmas in comparison to quantum statistical approaches and to experimental data.

Journal Article↗
PubMed2003

Heat-induced transformation of nanodiamond into a tube-shaped fullerene: a molecular dynamics simulation.

Heat-induced structural transformation in nanodiamond of diameter approximately 1.4 nm is investigated by tight-binding molecular dynamics simulations using the environment-dependent tight-binding carbon potential. The nanodiamond is found to transform into a tube-shaped fullerene via annealing. Three interesting mechanisms for promoting inner carbon atoms of the nanodiamond into the surface carbon atoms of the tubular structure are observed. The "flow-out" mechanism prevails at temperatures lower than 2500 K and the "direct adsorption" and "push-out" mechanisms are observed at higher temperatures.

Journal Article↗
PubMed1997

Conformational analysis and molecular dynamics simulation of alpha-(1-->2) and alpha-(1-->3) linked rhamnose oligosaccharides: reconciliation with optical rotation and NMR experiments.

Molecular mechanics and dynamics calculations were carried out on the disaccharides alpha-L-Rhap-(1-->2)-alpha-L-Rhap-(1-->OMe) (1) and alpha-L-Rhap-(1-->3)-alpha-L-Rhap-(1-->OMe) (2), and the trisaccharide alpha-L-Rhap-(1-->2)-alpha-L-Rhap-(1--> 3)-alpha-L-Rhap-(1-->OMe) (3). The semiflexible conformational behavior of these molecules was characterized by the occupation of a combination of different glycosidic linkage and side-chain conformational positions whose relative occupations were sensitive to dielectric screening. Molecular dynamics simulations of the trisaccharide 3 showed little difference between the linkage conformations in the trisaccharide and the component disaccharides 1 and 2 Experimental optical rotation data of 1 and 2 were obtained as a function of temperature in varying solvents. The molecular models were combined with the semiempirical theory of Stevens and Sathyanarayana to yield calculated optical rotations. Interpretation of the data of both 1 and 2 implied that a combination of conformations, both in glycosidic and side-chain positions, could explain the experimental data. Solvents effects were important in influencing the conformational mix and averaged optical rotation. Three-bond heteronuclear coupling constants 3JC H were obtained for the glycosidic linkages of 1 and 2 in D2O and DMSO. Analysis of the coupling constants with a Karplus curve showed that small reductions in the glycosidic torsion angles of the conformations of the models used here of ca. 10 degrees-15 degrees in phi and 5 degrees-10 degrees in psi were required to give better agreement with experiment; a combination of conformations for both 1 and 2 was consistent with the data. There was a negligible influence on the coupling constants of 1 on changing the solvent from D2O to DMSO.

Carbohydrate Conformation↗
PubMed2003

Molecular dynamics simulations of hydrated unsaturated lipid bilayers in the liquid-crystal phase and comparison to self-consistent field modeling.

Molecular dynamics simulations, using the collision dynamics method, were carried out for hydrated bilayers of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (18:0/18:1 omega 9cis PC, SOPC) and 1-stearoyl-2-docosahexaenoyl-sn-glycero-3-phosphatidylcholine (18:0/22:6 omega 3cis PC, SDPC). The simulation cells of the two bilayers consisted of 96 SOPC (or SDPC) molecules and 2304 water molecules: 48 lipid molecules per layer and 24 H2O molecules per lipid. The water was modeled by explicit TIP3P water molecules. The C-H bond-order-parameter -S(CH) profiles of the hydrocarbon tails, the bond orientation distribution functions and the root-mean-square values of the positional fluctuations of the lipid chain carbons were calculated. Simulation results are compared to the available experimental data and to other computer investigations of these lipid molecules. Several results of molecular-level self-consistent field calculations of these bilayers are also presented. Both theoretical methods reveal the same main characteristic features of the order-parameter profiles for the given bilayers. Some aspects of the physical properties of unsaturated lipids and their biological significance are discussed.

Journal Article↗
PubMed1989

An NMR-based molecular dynamics simulation of the interaction of the lac repressor headpiece and its operator in aqueous solution.

The results of a 125 psec molecular dynamics simulation of a lac headpiece-operator complex in aqueous solution are reported. The complex satisfies essentially all experimental distance information derived from two-dimensional nuclear magnetic resonance (2-D-NMR) studies. The interaction between lac repressor headpiece and its operator is based on many direct- and water-mediated hydrogen bonds and nonpolar contacts which allow the formation of a tight complex. No stable hydrogen bonds between side chains and bases are found, while specific contacts occur between both nonpolar groups and, to a lesser extent, through water-mediated hydrogen bonds. The simulated complex structure in water is intrinsically stable without application of nuclear Overhauser effect (NOE) distance restraints, while being compatible with most of the available biochemical, genetic, and chemically induced dynamic nuclear polarization (CIDNP) data.

Binding Sites↗
PubMed2004

Structure of the His44 --> Ala single point mutant of the distal finger motif of HIV-1 nucleocapsid protein: a combined NMR, molecular dynamics simulation, and fluorescence study.

The nucleocapsid protein (NCp7) of human immunodeficiency virus type 1 (HIV-1) contains two highly conserved CCHC zinc fingers that strongly bind Zn(2+) through coordination of one His and three Cys residues. It has been suggested that NCp7 function is conformation specific since substitution of any of the zinc coordinating residues in the zinc finger motifs leads to subsequent loss of viral infectivity. To further determine the structural requirements necessary for this specific conformation, we investigated by (1)H 2D NMR and molecular dynamics simulations the structure of the distal finger motif of NCp7 in which the zinc coordinating amino acid, His 44, was substituted by a noncoordinating Ala residue. While the fold of the N-terminal part of this mutated peptide was similar to that of the native peptide, an increased lability and significant conformational changes were observed in the vicinity of the His-to-Ala mutation. Moreover, molecular dynamics simulations suggested a mechanism by which the variant peptide can bind zinc ion even though one zinc-coordinating amino acid was lacking. Using the fluorescence of the naturally occurring Trp37 residue, the binding affinity of the variant peptide to the (TG)(3) model oligonucleotide was found to be decreased by about 2 orders of magnitude with respect with the native peptide. Modeling of the DNA:NCp7 complex using structures of the variant peptide suggests that the residues forming a hydrophobic cleft in the native protein are improperly oriented for efficient DNA binding by the variant peptide.

Alanine↗
PubMed2001

Molecular dynamics simulations of alpha-D-Manp-(1-->3)-beta-D-Glcp-OMe in methanol and in dimethyl sulfoxide solutions.

Molecular dynamics simulations have been performed of the disaccharide alpha-D-Manp-(1-->3)-beta-D-Glcp-OMe in two different solvents, namely in methanol and in dimethyl sulfoxide. The conformation of the disaccharide is similar to that previously determined in water. The three-dimensional structure around the solute was investigated by geometric hydrogen bonding criteria, radial distribution functions, coordination number analysis, residence times for hydrogen bonds, and spatial distribution functions. Differences and similarities between methanol and the aprotic dimethyl sulfoxide as solvent are analyzed.

Carbohydrate Conformation↗
PubMed2001

Homology modeling using multiple molecular dynamics simulations and docking studies of the human androgen receptor ligand binding domain bound to testosterone and nonsteroidal ligands.

To facilitate the rational design of novel and more potent androgen receptor ligands, three-dimensional models for the human androgen receptor ligand binding domain bound to testosterone have been developed. These models of the androgen receptor were based on the crystal structure of the highly homologous human progesterone receptor ligand binding domain. The homology modeled androgen receptor was refined using unrestrained multiple molecular dynamics simulations in explicit solvent. Key H-bonding partners with the 17-hydroxy group and 3-keto group of testosterone are Asn705 and Thr877, and Gln711 and Arg752, respectively. These models show the presence of a unique unoccupied cavity within the androgen receptor binding pocket which may be valuable in the development of novel selective androgen receptor ligands. A qualitative analysis of amino acid mutations within the hAR binding pocket that affect ligand binding are consistent with these androgen receptor models. In addition to testosterone, the binding modes of several hydroxyflutamide-like nonsteroidal ligands for the androgen receptor are investigated using flexible docking with FlexX followed by refinement of the initial complexes with molecular dynamics simulations. These docking studies indicate that Asn705 is an important determinant in binding hydroxyflutamide and its derivatives by participating in H-bond interactions with the alpha-hydroxy moiety of these ligands. In addition, the nitro functionality mimics the 3-keto group of the natural ligand testosterone and is involved in H-bonding interactions with Gln711 and Arg752. From these docking studies, we suggest a mechanism for the enantioselective binding of chiral hydroxyflutamide derivatives and expand upon the previously reported structure-activity relationship for hydroxyflutamide and its derivatives.

Amino Acid Substitution↗
PubMed1999

Long-term molecular dynamics simulation of copper azurin: structure, dynamics and functionality.

A long-term molecular dynamics simulation (1.1 ns), at 300 K, of fully hydrated azurin has been performed to put into relationship the protein dynamics to functional properties with particular attention to those structural elements involved in the electron transfer process. A detailed analysis of the root mean square deviations and fluctuations and of the intraprotein H-bonding pattern has allowed us to demonstrate that a rigid arrangement of the beta-stranded protein skeleton is maintained during the simulation run, while a large mobility is registered in the solvent-exposed connecting regions (turns) and in the alpha-helix. Moreover, the structural elements, likely involved in the electron transfer path, show a stable H-bonding arrangement and low fluctuations. Analysis of the dynamical cross-correlation map has revealed the existence of correlated motions among residues connected by hydrogen bonds and of correlated and anti-correlated motions between regions which are supposed to be involved in the functional process, namely the hydrophobic patch and the regions close to the copper reaction center. The results are briefly discussed also in connection to the current through-bond tunneling model for the electron transfer process. Finally, a comparison with the structural and the dynamical behaviour of plastocyanin, whose structure and functional role are very similar to those of azurin, has been performed.

Journal Article↗
PubMed2003

Similarity of force-induced unfolding of apomyoglobin to its chemical-induced unfolding: an atomistic molecular dynamics simulation approach.

We have compared force-induced unfolding with traditional unfolding methods using apomyoglobin as a model protein. Using molecular dynamics simulation, we have investigated the structural stability as a function of the degree of mechanical perturbation. Both anisotropic perturbation by stretching two terminal atoms and isotropic perturbation by increasing the radius of gyration of the protein show the same key event of force-induced unfolding. Our primary results show that the native structure of apomyoglobin becomes destabilized against the mechanical perturbation as soon as the interhelical packing between the G and H helices is broken, suggesting that our simulation results share a common feature with the experimental observation that the interhelical contact is more important for the folding of apomyoglobin than the stability of individual helices. This finding is further confirmed by simulating both helix destabilizing and interhelical packing destabilizing mutants.

Animals↗
PubMed2005

Influence of homologous disaccharides on the hydrogen-bond network of water: complementary Raman scattering experiments and molecular dynamics simulations.

A comparative investigation of trehalose, sucrose, and maltose in water solution has been performed using Raman scattering experiments and Molecular Dynamics simulations. From the analysis of the O-H stretching region in the [2500,4000] cm(-1) Raman spectral range, which includes for the first time the contribution of 'free' water, and the statistical distribution of water HB probabilities from MD simulations, this study confirms the privileged interaction of trehalose with water above a peculiar threshold weight concentration of about 30%. The role of the hydration number of sugars--found higher for trehalose--on the destructuring effect of the water hydrogen bond network is also addressed. The analysis of the water O-H-O bending spectral range [1500,1800] cm(-1) reveals a change of the homogeneity of water molecules influenced by sugars, but the three investigated sugars are found to behave similarly.

Computer Simulation↗
PubMed1998

Structural characteristics of 2'-O-(2-methoxyethyl)-modified nucleic acids from molecular dynamics simulations.

The structure and physical properties of 2'-sugar substituted O -(2-methoxyethyl) (MOE) nucleic acids have been studied using molecular dynamics simulations. Nanosecond simulations on the duplex MOE[CCAACGTTGG]-r[CCAACGUUGG] in aqueous solution have been carried out using the particle mesh Ewald method. Parameters for the simulation have been developed from ab initio calculations on dimethoxyethyl fragments in a manner consistent with the AMBER 4.1 force field database. The simulated duplex is compared with the crystal structure of the self-complementary duplex d[GCGTATMOEACGC]2, which contains a single modification in each strand. Structural details from each sequence have been analyzed to rationalize the stability imparted by substitution with 2'- O -(2-methoxyethyl) side chains. Both duplexes have an A-form structure, as indicated by several parameters, most notably a C3' endo sugar pucker in all residues. The simulated structure maintains a stable A-form geometry throughout the duration of the simulation with an average RMS deviation of 2.0 A from the starting A-form structure. The presence of the 2' substitution appears to lock the sugars in the C3' endo conformation, causing the duplex to adopt a stable A-form geometry. The side chains themselves have a fairly rigid geometry with trans , trans , gauche +/- and trans rotations about the C2'-O2', O2'-CA', CA'-CB' and CB'-OC' bonds respectively.

Base Sequence↗
PubMed1991

Motion of aromatic side chains, picosecond fluorescence, and internal energy transfer in Escherichia coli thioredoxin studied by site-directed mutagenesis, time-resolved fluorescence spectroscopy, and molecular dynamics simulations.

We have determined the picosecond fluorescence of the four aromatic amino acid residues (W28, W31, Y49, and Y70) in wild-type Escherichia coli thioredoxin (wt Trx) and a mutant Trx with W31 replaced by phenylalanine, Trx-W28-W31F. The internal motions of the four aromatic side chains were also analyzed. We examined the possibility of using internal energy transfer from tyrosine to tryptophan as a measure of long-range distances. The major features of the lifetime distribution of tryptophan fluorescence were unchanged in the W31F mutation, indicating that the environment of W28 is similar in both wt Trx and Trx-W28-W31F. However, the mutation of W31F changed the mobility of W28, situated close to the active-site disulfide/dithiol, but not the mobility of two tyrosines, Y49 and Y70, situated on the other side of the molecule. The mobility of the two tyrosine residues increased upon reduction of the active-site disulfide, indicating a looser structure with reduction. This increased motion could also be seen from molecular dynamics simulations. The change in energy transfer rates, as judged by tyrosine fluorescence lifetimes, was in agreement with energy transfer rates calculated from the molecular dynamics simulations. The anisotropy of tryptophan and tyrosine fluorescence could be separated in three parts: (I) overall rotation of the protein (10(-9)s), (II) internal mobility of side chains (10(-10)s), and (III) a very fast relaxation (10(-12)s). We can only experimentally detect this very fast relaxation when the internal motion is not present.

Amino Acids↗
PubMed2001

Variation in aspects of cysteine proteinase catalytic mechanism deduced by spectroscopic observation of dithioester intermediates, kinetic analysis and molecular dynamics simulations.

The possibility of a slow post-acylation conformational change during catalysis by cysteine proteinases was investigated by using a new chromogenic substrate, N-acetyl-Phe-Gly methyl thionoester, four natural variants (papain, caricain, actinidin and ficin), and stopped-flow spectral analysis to monitor the pre-steady state formation of the dithioacylenzyme intermediates and their steady state hydrolysis. The predicted reversibility of acylation was demonstrated kinetically for actinidin and ficin, but not for papain or caricain. This difference between actinidin and papain was investigated by modelling using QUANTA and CHARMM. The weaker binding of hydrophobic substrates, including the new thionoester, by actinidin than by papain may not be due to the well-known difference in their S2-subsites, whereby that of actinidin in the free enzyme is shorter due to the presence of Met211. Molecular dynamics simulation suggests that during substrate binding the sidechain of Met211 moves to allow full access of a Phe sidechain to the S2-subsite. The highly anionic surface of actinidin may contribute to the specificity difference between papain and actinidin. During subsequent molecular dynamics simulations the P1 product, methanol, diffuses rapidly (over<8 ps) out of papain and caricain but 'lingers' around the active centre of actinidin. Uniquely in actinidin, an Asp142-Lys145 salt bridge allows formation of a cavity which appears to constrain diffusion of the methanol away from the catalytic site. The cavity then undergoes large scale movements (over 4.8 A) in a highly correlated manner, thus controlling the motions of the methanol molecule. The changes in this cavity that release the methanol might be those deduced kinetically.

Catalysis↗
PubMed1996

How long does DNA keep the memory of its conformation? A time-dependent canonical correlation analysis of molecular dynamics simulation.

The time dependence of the correlation between motions of different parts of DNA is analyzed from a 200 ps molecular dynamics simulation of the double-stranded self-complementary d(CTGATCAG) in the B form. Each nucleotide is decomposed into three subunits corresponding to the furanose ring (SU), the base (BA), and the backbone (SK). The motion of each subunit is considered as the superimposition of rigid body translation, rigid body rotation, and internal deformation. Canonical time-dependent correlation functions calculated with coordinates describing the different components of the subunits motion are defined and computed. This allows us to probe how long a particular type of motion of one subunit influences the other types of motions of other subunits (cross correlation functions) or how long a particular subunit keeps the memory of its own conformation or location (autocorrelation functions). From autocorrelation analysis it is found that deformation decorrelates within a few tenths of picoseconds, rotational correlation times are on the order of 8 ps, while translational motions are long-time correlated. The deformation of a subunit is not correlated to the deformation of another one (at the 200 ps time scale of our simulation), but influences slightly their translation and orientation as time increases.

Base Sequence↗