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.

SourceResearch connectionFrom yearThrough year

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

At least 1,153 records · Page 64Linked to original sources

PubMed1999

Molecular dynamics simulations of the Ras:Raf and Rap:Raf complexes.

The protein Raf is an immediate downstream target of Ras in the MAP kinase signalling pathway. The complex of Ras with the Ras-binding domain (RBD) of Raf has been modelled by homology to the (E30D,K31E)-Rap1A:RBD complex, and both have been subjected to multiple molecular dynamics simulations in solution. While both complexes are stable, several rearrangements occur in the Ras:RBD simulations: the RBD loop 100-109 moves closer to Ras, Arg73 in the RBD moves towards Ras to form a salt bridge with Ras-Asp33, and Loop 4 of the Ras switch II region shifts upwards toward the RBD. The Ras:RBD interactions (including the RBD-Arg73 interaction) are consistent with available NMR and mutagenesis data on the Ras: RBD complex in solution. The Ras switch II region does not interact directly with the RBD, although indirect interactions exist through the effector domain and bridging water molecules. No large-scale RBD motion is seen in the Ras:RBD complex, compared to the Rap:RBD complex, to suggest an allosteric activation of Raf by Ras. This may be because the Raf kinase domain (whose structure is unknown) is not included in the model.

Amino Acid Sequence↗

PubMed2002

Molecular dynamics simulations on discoidal HDL particles suggest a mechanism for rotation in the apo A-I belt model.

Apolipoprotein A-I (apo A-I) is the major protein component of high-density lipoprotein (HDL) particles. Elevated levels of HDL in the bloodstream have been shown to correlate strongly with a reduced risk factor for atherosclerosis. Molecular dynamics simulations have been carried out on three separate model discoidal high-density lipoprotein particles (HDL) containing two monomers of apo A-I and 160 molecules of palmitoyloleoylphosphatidylcholine (POPC), to a time-scale of 1ns. The starting structures were on the basis of previously published molecular belt models of HDL consisting of the lipid-binding C-terminal domain (residues 44-243) wrapped around the circumference of a discoidal HDL particle. Subtle changes between two of the starting structures resulted in significantly different behavior during the course of the simulation. The results provide support for the hypothesis of Segrest et al. that helical registration in the molecular belt model of apo A-I is modulated by intermolecular salt bridges. In addition, we propose an explanation for the presence of proline punctuation in the molecular belt model, and for the presence of two 11-mer helical repeats interrupting the otherwise regular pattern of 22-mer helical repeats in the lipid-binding domain of apo A-I.

1,2-Dipalmitoylphosphatidylcholine↗

PubMed2005

Structures of a peptide fragment of beta2-microglobulin studied by replica-exchange molecular dynamics simulations - towards the understanding of the mechanism of amyloid formation.

We investigate in detail the structural properties of the monomeric peptide fragment that corresponds to residues 21-31 of beta(2)-microglobulin. As a first step towards the understanding of the mechanism of the amyloid formation, we have performed a replica-exchange molecular dynamics simulation of this peptide with explicit water molecules. We analyze various structural properties as functions of temperature. Although the corresponding part of the native protein is a fully extended beta-strand, our results show that beta-hairpin structures are formed with high frequency around 310 K. We conjecture that this beta-hairpin formation is closely related to the amyloid fibrillogenesis.

Amyloid↗

PubMed2006

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface.

Hydrophobins are small ( approximately 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily beta-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in alpha-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Amino Acid Sequence↗

PubMed2006

Molecular dynamics simulations of model trans-membrane peptides in lipid bilayers: a systematic investigation of hydrophobic mismatch.

Hydrophobic mismatch, which is the difference between the hydrophobic length of trans-membrane segments of a protein and the hydrophobic width of the surrounding lipid bilayer, is known to play a role in membrane protein function. We have performed molecular dynamics simulations of trans-membrane KALP peptides (sequence: GKK(LA)nLKKA) in phospholipid bilayers to investigate hydrophobic mismatch alleviation mechanisms. By varying systematically the length of the peptide (KALP15, KALP19, KALP23, KALP27, and KALP31) and the lipid hydrophobic length (DLPC, DMPC, and DPPC), a wide range of mismatch conditions were studied. Simulations of durations of 50-200 ns show that under positive mismatch, the system alleviates the mismatch predominantly by tilting the peptide and to a smaller extent by increased lipid ordering in the immediate vicinity of the peptide. Under negative mismatch, alleviation takes place by a combination of local bilayer bending and the snorkeling of the lysine residues of the peptide. Simulations performed at a higher peptide/lipid molar ratio (1:25) reveal slower dynamics of both the peptide and lipid relative to those at a lower peptide/lipid ratio (1:128). The lysine residues have favorable interactions with specific oxygen atoms of the phospholipid headgroups, indicating the preferred localization of these residues at the lipid/water interface.

Algorithms↗

PubMed2004

Molecular dynamics simulations of forced conformational transitions in 1,6-linked polysaccharides.

Recent atomic force microscopy stretching measurements of single polysaccharide molecules suggest that their elasticity is governed by force-induced conformational transitions of the pyranose ring. However, the mechanism of these transitions and the mechanics of the pyranose ring are not fully understood. Here we use steered molecular dynamics simulations of the stretching process to unravel the mechanism of forced conformational transitions in 1,6 linked polysaccharides. In contrast to most sugars, 1,6 linked polysaccharides have an extra bond in their inter-residue linkage, C5-C6, around which restricted rotations occur and this additional degree of freedom increases the mechanical complexity of these polymers. By comparing the computational results with the atomic force microscopy data we determine that forced rotations around the C5-C6 bond have a significant and different impact on the elasticity of alpha- and beta-linked polysaccharides. Beta-linkages of a polysaccharide pustulan force the rotation around the C5-C6 bonds and produce a Hookean-like elasticity but do not affect the conformation of the pyranose rings. However, alpha-linkages of dextran induce compound conformational transitions that include simultaneous rotations around the C5-C6 bonds and chair-boat transitions of the pyranose rings. These previously not-recognized transitions are responsible for the characteristic plateau in the force-extension relationship of dextran.

Carbohydrate Conformation↗

PubMed2004

Molecular dynamics simulations of peptides from the central domain of smooth muscle caldesmon.

The central domain of smooth muscle caldesmon contains a highly charged region consisting of ten 13-residue repeats. Experimental evidence obtained from the intact protein and fragments thereof suggests that this entire region forms a single stretch of stable alpha-helix. We have carried out molecular dynamics simulations on peptides consisting of one, two and three repeats to examine the mechanism of alpha-helical stability of the central domain at the atomic level. All three peptides show high helical stability on the timescale of the MD simulations. Deviations from alpha-helical structure in all the simulations arise mainly from the formation of long stretches of pi-helix. Interconversion between alpha-helical and pi-helical conformations occurs through insertion of water molecules into alpha-helical hydrogen bonds and subsequent formation of reverse turns. The alpha-helical structure is stabilized by electrostatic interactions (salt bridges) between oppositely charged sidechains with i,i+4 spacings, while the pi-helix is stabilized by i,i+5 salt bridge interactions. Possible i,i+3 salt bridges are of minor importance. There is a strong preference for salt bridges with a Glu residue N-terminal to a basic sidechain as compared to the opposite orientation. In the double and triple repeat peptides, strong i,i+4 salt bridges exist between the last Glu residue of one repeat and the first Lys residue of the next. This demonstrates a relationship between the repetitive nature of the central domain sequence and its ability to form very long stretches of alpha-helical structure.

Calmodulin-Binding Proteins↗

PubMed1999

Molecular dynamics simulations of polypeptide conformations in water: A comparison of alpha, beta, and poly(pro)II conformations.

A significant fraction of the so-called "random coil" residues in globular proteins exists in the left-handed poly(Pro)II conformation. In order to compare the behavior of this secondary structure with that of the other regular secondary structures, molecular dynamics simulations, with the GROMOS suite of programs, of an alanine octapeptide in water, in alpha-helix, beta-strand, and left-handed poly(Pro)II conformations, have been performed. Our results indicate a limited flexibility for the alpha-helix conformation and a relatively larger flexibility for the beta-strand and poly(Pro)II conformations. The behavior of oligopeptides with a starting configuration of beta-strand and poly(Pro)II conformations, both lacking interchain hydrogen bonds, were similar. The (phi, psi) angles reflect a continuum of structures including both beta and P(II) conformations, but with a preference for local P(II) regions. Differences in the network of water molecules involved in hydrogen bonding with the backbone of the polypeptide were observed in local regions of beta and P(II) conformations. Such water bridges help stabilize the P(II) conformation relative to the beta conformation. Proteins 1999;36:400-406.

Alanine↗

PubMed2006

Moving solvated electrons with light: nonadiabatic mixed quantum/classical molecular dynamics simulations of the relocalization of photoexcited solvated electrons in tetrahydrofuran (THF).

Motivated by recent ultrafast spectroscopic experiments [Martini et al., Science 293, 462 (2001)], which suggest that photoexcited solvated electrons in tetrahydrofuran (THF) can relocalize (that is, return to equilibrium in solvent cavities far from where they started), we performed a series of nonequilibrium, nonadiabatic, mixed quantum/classical molecular dynamics simulations that mimic one-photon excitation of the THF-solvated electron. We find that as photoexcited THF-solvated electrons relax to their ground states either by continuous mixing from the excited state or via nonadiabatic transitions, approximately 30% of them relocalize into cavities that can be over 1 nm away from where they originated, in close agreement with the experiments. A detailed investigation shows that the ability of excited THF-solvated electrons to undergo photoinduced relocalization stems from the existence of preexisting cavity traps that are an intrinsic part of the structure of liquid THF. This explains why solvated electrons can undergo photoinduced relocalization in solvents like THF but not in solvents like water, which lack the preexisting traps necessary to stabilize the excited electron in other places in the fluid. We also find that even when they do not ultimately relocalize, photoexcited solvated electrons in THF temporarily visit other sites in the fluid, explaining why the photoexcitation of THF-solvated electrons is so efficient at promoting recombination with nearby scavengers. Overall, our study shows that the defining characteristic of a liquid that permits the photoassisted relocalization of solvated electrons is the existence of nascent cavities that are attractive to an excess electron; we propose that other such liquids can be found from classical computer simulations or neutron diffraction experiments.

Computer Simulation↗

PubMed2004

Opening mechanism of G.T/U pairs in DNA and RNA duplexes: a combined study of imino proton exchange and molecular dynamics simulation.

The opening pathway of wobble pairs dG.T and rG.U has been investigated in four DNA and two RNA duplexes. Using NMR spectroscopy, we measured the imino proton exchange of both G(H1) and T/U(H3), catalyzed by ammonia, tris, and OH(-), and we calculated the free energy surface related to G.T/U opening by molecular dynamics simulations. Taken together the experimental and theoretical results, we suggest that wobble pairs open through a coupled rotation of the bases toward the major groove where exchange of both imino protons takes place with the surrounding water.

Ammonia↗

PubMed2004

Beta-hairpin conformation of fibrillogenic peptides: structure and alpha-beta transition mechanism revealed by molecular dynamics simulations.

Understanding the conformational transitions that trigger the aggregation and amyloidogenesis of otherwise soluble peptides at atomic resolution is of fundamental relevance for the design of effective therapeutic agents against amyloid-related disorders. In the present study the transition from ideal alpha-helical to beta-hairpin conformations is revealed by long timescale molecular dynamics simulations in explicit water solvent, for two well-known amyloidogenic peptides: the H1 peptide from prion protein and the Abeta(12-28) fragment from the Abeta(1-42) peptide responsible for Alzheimer's disease. The simulations highlight the unfolding of alpha-helices, followed by the formation of bent conformations and a final convergence to ordered in register beta-hairpin conformations. The beta-hairpins observed, despite different sequences, exhibit a common dynamic behavior and the presence of a peculiar pattern of the hydrophobic side-chains, in particular in the region of the turns. These observations hint at a possible common aggregation mechanism for the onset of different amyloid diseases and a common mechanism in the transition to the beta-hairpin structures. Furthermore the simulations presented herein evidence the stabilization of the alpha-helical conformations induced by the presence of an organic fluorinated cosolvent. The results of MD simulation in 2,2,2-trifluoroethanol (TFE)/water mixture provide further evidence that the peptide coating effect of TFE molecules is responsible for the stabilization of the soluble helical conformation.

Animals↗

PubMed2002

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations.

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.

Computer Simulation↗

PubMed2003

Insight into the mechanism of the IMP-1 metallo-beta-lactamase by molecular dynamics simulations.

Two models, a purely nonbonded model and a cationic dummy atom approach, were examined for the modeling of the binuclear zinc-containing IMP-1 metallo-beta-lactamase in complex with a mercaptocarboxylate inhibitor. The cationic dummy atom approach had substantial advantages as it maintained the initial, experimentally determined geometry of the metal-containing active site during molecular dynamics simulations in water. The method was extended to the modeling of the free enzyme and the enzyme in complex with a cephalosporin substrate docked in an intermediate structure. For all three systems, the modeled complexes and the tetrahedral coordination of the zinc ions were stable. The average zinc-zinc distance increased by approximately 1 A in the substrate complex compared with the inhibitor complex and the free enzyme in which a hydroxide ion acts as a bridging ligand. Thus, the zinc ions are predicted to undergo a back and forth movement upon the cycle of hydrolysis. In contrast to previous assumptions, no interaction of the Asn167 side chain with the bound cephalosporin substrate was observed. Our observations are in agreement with quantum-mechanical calculations and experimental data and indicate that the cationic dummy atom approach is useful to model zinc-containing metallo-beta-lactamases as free proteins, in complex with inhibitors and in complex with substrates.

Binding Sites↗

PubMed1991

Molecular dynamics simulation of the hydration shell of a B-DNA decamer reveals two main types of minor-groove hydration depending on groove width.

The conformation of the self-complementary B-DNA decamer C-C-A-A-C-G-T-T-G-G is known from a high-resolution x-ray crystal structure analysis. Molecular dynamics simulation of the hydration shell of the decamer has revealed two main types of minor-groove hydration, depending on groove width. The narrow part of the minor groove has a spine of hydration analogous to that described for the A + T-rich center of the minor groove in the dodecamer C-G-C-G-A-A-T-T-C-G-C-G [Drew, H. R. & Dickerson, R. E. (1981) J. Mol. Biol. 151, 535-556], the first hydration layer of which contains one water molecular per base pair. In contrast, in the wide part of the minor groove, each base is hydrated individually, water molecules lying predominantly in the base plane. In intermediate-width regions, preferred water-molecule sites are shifted away from the base plane in a 3'-to-5' direction. This shift becomes more pronounced as the minor groove narrows, until the two water molecules lie approximately midway between base pairs. If the minor groove is narrowed still further, it accommodates only one water molecule, and the hydration transforms to the well-known water spine. The observed pattern agrees with available crystallographic data and with our earlier calculations. The results confirm the assumption that preferred positions of water oxygens in the minor groove depend predominantly on groove width rather than on base sequence. However, the location of water hydrogens, and the network of hydrogen bonding, can depend on base sequence. We suggest a simple explanation of water-spine formation in the narrow minor groove of a random DNA sequence. The spine of hydration may be a property of the minor groove of overwound variants of B-DNA, the C and D forms, for which the middle part of the decamer C-C-A-A-C-G-T-T-G-G can serve as a model.

Base Sequence↗

PubMed1993

Molecular dynamics simulation of the stability of a 22-residue alpha-helix in water and 30% trifluoroethanol.

A molecular dynamics (MD) simulation was performed on the alpha-helix H8-HC5, the C-terminal part of myoglobin (residue 132-153), under periodic boundary conditions in two different solutions, water and water with 30% (v/v) 2,2,2-trifluoroethanol (TFE), at 300 K to investigate the stability of the helix. In both simulations, the initial configuration was a canonical right-handed alpha-helix. In the course of the MD trajectory in water (200 ps), the helix clearly destabilized and began to unfold after 100 ps. In the TFE solution, two stable parts of helical regions were observed after 70 ps of a 200-ps MD simulation, supporting the notion that TFE acts as a structure-forming solvent.

Amino Acid Sequence↗

PubMed2004

Structure and dynamics of the CrIII ion in aqueous solution: Ab initio QM/MM molecular dynamics simulation.

Structural and dynamical properties of the Cr(III) ion in aqueous solution have been investigated using a combined ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation. The hydration structure of Cr(III) was determined in terms of radial distribution functions, coordination numbers, and angular distributions. The QM/MM simulation gives coordination numbers of 6 and 15.4 for the first and second hydration shell, respectively. The first hydration shell is kinetically very inert but by no means rigid and variations of the first hydration shell geometry lead to distinct splitting in the vibrational spectra of Cr(H(2)O)(6) (3+). A mean residence time of 22 ps was obtained for water ligands residing in the second hydration shell, which is remarkably shorter than the experimentally estimated value. The hydration energy of -1108 +/- 7 kcal/mol, obtained from the QM/MM simulation, corresponds well to the experimental hydration enthalpy value.

Journal Article↗

PubMed2005

Effect of cholesterol on DMPC phospholipid membranes and QSAR model construction in membrane-interaction QSAR study through molecular dynamics simulation.

In this study, both pure DMPC and DMPC/cholesterol mixed membrane monolayer were built to compare the physical-chemical properties and dynamics properties through molecular dynamics simulation and normal-mode analysis. The results show that the addition of cholesterol decreases the area of per molecule of membrane, increases the lipid amplitude motion, and changes the solute diffusion coefficient. It is also found that the addition of cholesterol greatly changes the solute-membrane 1,4-nonbonded interaction energy (deltaE14). MI-QSAR models were constructed based on solute-membrane interaction energy descriptors and other intramolecular descriptors. The results show that deltaE14 substitutes deltaE(HB) as the second important descriptor compared with the previous study. Final results suggest that short range solute-membrane interaction energy changes due to the uptake of the solute may play an important decision on permeability in DMPC/cholesterol membrane. A test set was applied to evaluate the predictivity of MI-QSAR models. The result suggests that the combination of F(H2O) and deltaE14 not only improves r2 and q2, but also greatly improves the model predictivity. Based on the combination of q2 and r(pre)2 values, a two-term model is better used to predict the solute permeability in this study (r2 = 0.859, q2 = 0.803, and r(pre)2 = 0.540). Due to the small sample both in training set and test set, more datasets are necessary to make a final decision about the model construction and prediction.

Chemical Phenomena↗

PubMed2006

Neutron diffraction data and molecular dynamics simulations of the molten mixture Ag(Br0.7I0.3).

The structure factors of the ionic liquid mixture Ag(Br(0.7)I(0.3)) at three temperatures, 723, 923, and 1023 K, as well as of the pure molten AgI at 923 K and the pure molten AgBr at 773 and 923 K, were studied experimentally and by means of molecular dynamics simulations. The experiments were carried out using the high intensity total scattering time-of-flight spectrometer, HIT-II, at the KENS spallation neutron source in Japan. The experimental data are very reliable, with the possible exception of the small momentum transfer region, whose accessibility is limited by neutron energy and detector positions. The simulations made use of the semiempirical rigid ion potentials of the Vashishta-Rahman [Phys. Rev. Lett. 40, 1337 (1978)] type using a new set of parameters appropriate for the mixture. Within the known constraints of the pairwise rigid ion potentials, the simulated structure factors are in fair agreement with experiment. The results for the pair distribution functions suggest that the molten mixture retains the superionic character found in previous calculations of both the AgI and AgBr melts. This suggestion is confirmed by the results for the self-diffusion coefficients. Values obtained for the ionic conductivities are also presented.

Journal Article↗