Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “molecular dynamics”

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 91 records · Page 5Linked to original sources
PubMed1991

Cooperative ligand reorientations in cytochrome c3: a molecular dynamics simulation.

Molecular dynamics simulations of a tetraheme cytochrome c3 were performed to investigate dynamic aspects of the motion of the axial heme iron ligands. It was found that persistent transitions between alternate axial imidazole orientations of the histidine incorporated in the CXXCH heme binding sequence occurred via correlated motions. The correlated motions involved virtually all of the atoms comprising the polypeptide backbone of the heme binding sequence as well as the histidine imidazole side-chain.

Amino Acid Sequence↗
PubMed2003

Effects of pressure on the fragile nature of fluorozirconates studied by molecular dynamics simulations.

Molecular dynamics simulations for a multicomponent zirconium fluoride, ZBEN (55ZrF(4)-17BaF(2)-5EuF(3)-23NaF) have been carried out at various temperatures and pressures including those of the molten and glassy states. The effect of pressure on the fragile nature of ZBEN was investigated. The relaxation times were calculated from the decay curves of the intermediate scattering function. The fragile nature of ZBEN was reflected by the obvious non-Arrhenius behavior in the temperature dependence of the relaxation time. The so-called "Angell plot" of the relaxation time was represented by a single curve irrespective of applied pressure; the fragile nature changes negligibly with pressure. This seems to contradict the enhancement of the connectivity of ZrF(n) polyhedra under high pressure. The origin of the pressure effects is discussed in comparison with those for molecular liquids and covalent-network glasses.

Journal Article↗
PubMed2006

Study of the inhibition of cyclin-dependent kinases with roscovitine and indirubin-3'-oxime from molecular dynamics simulations.

Molecular dynamics simulations were performed to elucidate the interactions of CDK2 and CDK5 complexes with three inhibitors: R-roscovitine, S-roscovitine, and indirubin-3'-oxime. The preference of the two complexes for R-roscovitine over the S enantiomer, as reported by the experiment, was also found by the simulations. More importantly, the simulations showed that the cause of the stronger affinity for the R enantiomer is the presence of an important hydrogen bond between R-roscovitine and the kinases not found with S-roscovitine. The simulations also showed two amino acid mutations in the active site of CDK5/R-roscovitine that favor binding-enhanced electrostatic contributions, making the inhibitor more effective for CDK5 than for CDK2. This suggests that the effectiveness of roscovitine-like inhibitors can be improved by enhancing their electrostatic interaction with the kinases. Finally, molecular mechanics-Possion-Boltzmann/surface area calculations of the CDK5/indirubin-3'-oxime system in both water-excluded and water-included environments gave significantly different electrostatic contributions to the binding. The simulations detected the displacement of a water molecule in the active site of the water-included CDK/indirubin-3'-oxime system. This resulted in a more conserved binding pattern than the water-excluded structure. Hence, in the design of new indirubin-like inhibitors, it is important to include the water molecule in the analysis.

Binding Sites↗
PubMed1991

Conformational space of alpha 1-4 glycosidic linkage: a molecular dynamics study.

Molecular Dynamics simulations have been carried out for 100 ps on crystal structure of beta-cyclodextrin in vacuo and with explicit inclusion of solvent at constant pressure and constant temperature using the GROMOS MD algorithm, with a time step of 0.005 ps. The conformational space of the glycosidic linkage was studied by calculating two virtual dihedrals connecting the successive glucose units for the 2000 structures saved during the two simulations. Three preferred regions for alpha 1-4 glycosidic linkage were found in both the simulations. The use of these virtual dihedral angles in representing the glycosidic linkage is also brought out from these studies.

Algorithms↗
PubMed1998

The M2 channel of influenza A virus: a molecular dynamics study.

Molecular dynamics simulations have been performed on a tetramer of the 25-residue (SSDPLVVAASIIGILHLILWILDRL) synthetic peptide [1] which contains the transmembrane domain of the influenza A virus M2 coat protein. The peptide bundle was initially assembled as a parallel alpha-helix bundle in the octane portion of a phase separated water/octane system, which provided a membrane-mimetic environment. A 4-ns dynamics trajectory identified a left-handed coiled coil state of the neutral bundle, with a water filled funnel-like structural motif at the N-terminus involving the long hydrophobic sequence. The neck of the funnel begins at V27 and terminates at H37, which blocks the channel. The C-terminus is held together by inter-helix hydrogen bonds and contains water below H37. Solvation of the S23 and D24 residues, located at the rim of the funnel, appears to be important for stability of the structure. The calculated average tilt of the helices in the neutral bundle is 27 +/- 5 degrees, which agrees well with recent NMR data.

Amino Acid Sequence↗
PubMed2004

Dependence of DNA polymerase replication rate on external forces: a model based on molecular dynamics simulations.

Molecular dynamics simulations are presented for a Thermus aquaticus (Taq) DNA polymerase I complex (consisting of the protein, the primer-template DNA strands, and the incoming nucleotide) subjected to external forces. The results obtained with a force applied to the DNA template strand provide insights into the effect of the tension on the activity of the enzyme. At forces below 30 pN a local model based on the parameters determined from the simulations, including the restricted motion of the DNA bases at the active site, yields a replication rate dependence on force in agreement with experiment. Simulations above 40 pN reveal large conformational changes in the enzyme-bound DNA that may have a role in the force-induced exonucleolysis observed experimentally.

Binding Sites↗
PubMed2005

Accommodation of a 1S-(-)-benzo[c]phenanthrenyl-N6-dA adduct in the Y-family Dpo4 DNA polymerase active site: structural insights through molecular dynamics simulations.

Molecular modeling and molecular dynamics simulations have been performed to elucidate feasible structures in the Y-family Dpo4 DNA polymerase for the 1S-(-)-trans-anti-B[c]Ph-N6-dA adduct, derived from the fjord region polycyclic aromatic hydrocarbon (PAH) benzo[c]phenanthrene. Three types of models were delineated as follows: an intercalation model, a model with the aromatic ring system in the polymerase major groove open pocket, and a -1 deletion major groove model. All four 2'-deoxyribonucleoside 5'-triphosphates (dNTPs) were considered in the first two cases, and a normal Watson-Crick partner positioned to have skipped the modified template was employed as the incoming dNTP in the -1 deletion case. The trajectories derived from the dynamics simulations were analyzed in detail to evaluate the extents of distortion for each system. Overall, our results suggest that the major groove model is the least distorted, followed by the -1 deletion model, while the intercalation model is perturbed the most. The syn-dGTP and syn-dATP mismatches opposite the lesion are well-accommodated in the major groove model, as is the normal Watson-Crick partner dTTP. The intercalation model appears most likely to impede the polymerase. More broadly, these models look reasonable for other PAH metabolite-derived adducts to adenine with similar 1S stereochemistry. Furthermore, these models suggest how error-prone translesion synthesis by Y-family polymerases might produce mutations that may play a role in the initiation of cancer.

Carcinogens↗
PubMed2004

Characterization of the conformational space of a triple-stranded beta-sheet forming peptide with molecular dynamics simulations.

Molecular dynamics (MD) simulations have been performed on a series of mutants of the 20 amino acid peptide Betanova in order to critically assess the ability of MD simulations to reproduce the folding and stability of small beta-sheet-forming peptides on currently accessible timescales. Simulations were performed in both water and in 40% methanol solution, using an explicit solvent model. The simulations suggest that all mutants adopt a wide range of conformations in solution, that the structures are highly flexible, and that stabilization of compact structures is due to a delicate balance of hydrophobic and polar side-chain interactions. Simulations longer than 100 ns, although not sufficient for a complete thermodynamic and kinetic description of the system, sample an ensemble of compact conformations characterized by the loss of ordered beta-sheet secondary structure. This suggests that no significant free energy barrier separates the different conformations available.

Algorithms↗
PubMed2003

The Mg2+ binding sites of the 5S rRNA loop E motif as investigated by molecular dynamics simulations.

Molecular dynamics simulations have been used to investigate the binding of Mg(2+) ions to the deep groove of the eubacterial 5S rRNA loop E. The simulations suggest that long-lived and specific water-mediated interactions established between the hydrated ions and the RNA atoms lining up the binding sites contribute to the stabilization of this motif. The Mg(2+) binding specificity is modulated by two factors: (i) a required electrostatic complementarity and (ii) a structural correspondence between the hydrated ion and its binding pocket that can be estimated by its degree of dehydration and the resulting number and lifetime of the intervening water-mediated contacts. Two distinct binding modes for pentahydrated Mg(2+) ions that result in a significant freezing of the tumbling motions of the ions are described, and mechanistic details related to the stabilization of nucleic acids by divalent ions are provided.

Binding Sites↗
PubMed1993

Mechanism for the rotamase activity of FK506 binding protein from molecular dynamics simulations.

Molecular dynamics (MD) and free energy perturbation (FEP) methods are used to study the binding and mechanism of isomerization of a tetrapeptide (AcAAPFNMe) by FK506 binding protein (FKBP). Detailed structures are predicted for the complexes of FKBP with the peptide in both ground-state and transition-state forms. The results support a mechanism of catalysis by distortion, where a large number of nonbonded interactions act together to stabilize preferentially the twisted transition state. The two most important groups for the catalysis are suggested to be Trp59 and Asp37, but several other groups are identified as directly or indirectly involved in the binding and catalysis. However, the structural results do not support the notion that the keto oxygen of the immunosuppressive agents FK506 and rapamycin mimics the oxygen for the twisted peptide bond in the FKBP-transition-state complex.

Amino Acid Sequence↗
PubMed1993

Protein hydration elucidated by molecular dynamics simulation.

Molecular dynamics (MD) simulation covering a wide range of hydration indicate that myoglobin is fully hydrated by 350 water molecules, in agreement with experiment. These waters, originally placed uniformly about the protein, form clusters that hydrate every charged group throughout the entire simulation. Some atoms in charged groups are hydrated by two water layers while 37% of the protein surface remains uncovered. The locations of the 350 waters are consistent with those of crystallographic waters resolved by x-ray and neutron diffraction. Hydration by 350 waters at 300 K stabilizes the conformation of carboxymyoglobin measured by x-ray diffraction throughout the entire protein, halves the rate of torsional transitions, and promotes alternative conformations for surface atoms. The glass transition observed experimentally in hydrated myoglobin near 220 K is also seen in the simulations and correlates with an increase in the number of dihedral angles undergoing transitions. The anharmonic protein motion above 220 K is enhanced by protein hydration.

Computer Simulation↗
PubMed2002

Mechanism by which 2,2,2-trifluoroethanol/water mixtures stabilize secondary-structure formation in peptides: a molecular dynamics study.

Molecular dynamics simulation techniques have been used to investigate the effect of 2,2,2-trifluoroethanol (TFE) as a cosolvent on the stability of three different secondary structure-forming peptides: the alpha-helix from Melittin, the three-stranded beta-sheet peptide Betanova, and the beta-hairpin 41-56 from the B1 domain of protein G. The peptides were studied in pure water and 30% (vol/vol) TFE/water mixtures at 300 K. The simulations suggest that the stabilizing effect of TFE is induced by the preferential aggregation of TFE molecules around the peptides. This coating displaces water, thereby removing alternative hydrogen-bonding partners and providing a low dielectric environment that favors the formation of intrapeptide hydrogen bonds. Because TFE interacts only weakly with nonpolar residues, hydrophobic interactions within the peptides are not disrupted. As a consequence, TFE promotes stability rather than inducing denaturation.

Amino Acid Sequence↗
PubMed1992

The conformational behaviour of complexes of alpha-cyclodextrin with p-chlorophenol and p-hydroxybenzoic acid in water as studied by molecular dynamics simulations.

Molecular dynamics simulations were performed to obtain information about the conformational behaviour and stabilization of alpha-cyclodextrin (alpha CD) complexes in water. Simulations of p-chlorophenol and p-hydroxybenzoic acid in alpha CD showed that the complex is a very flexible system. The guest compound rotates inside the cavity and partly moves in and out. alpha CD continuously adapts its conformation to the orientation of the guest compound (or vice versa): the hexagon of the glycosidic oxygen atoms is stretched parallel to the plane of the aromatic ring of the guest compound during 80% of the simulation. This suggests that Van der Waals interactions play an important role in the stabilization of the complex. Each intramolecular hydrogen bond between neighbouring glucose units is formed during 30-80% of the simulation. Hydrogen bonds between alpha CD and the guest compound, on the other hand, are rarely formed. Thus, intermolecular hydrogen bonding seems to play a minor role in the stabilization of alpha CD complexes.

Chlorophenols↗
PubMed2003

Gating of MscL studied by steered molecular dynamics.

Steered molecular dynamics simulations of the mechanosensitive channel of large conductance, MscL, were used to investigate how forces arising from membrane tension induce gating of the channel. A homology model of the closed form of MscL from Escherichia coli was subjected to external forces of 35-70 pN applied to residues near the membrane-water interface. The magnitude and location of these forces corresponded to those determined from the lateral pressure profile computed from a lipid bilayer simulation. A fully expanded state was obtained on the 10-ns timescale that revealed the mechanism for transducing membrane forces into channel opening. The expanded state agrees well with proposed models of MscL gating, in that it entails an irislike expansion of the pore accompanied by tilting of the transmembrane helices. The channel was most easily opened when force was applied predominantly on the cytoplasmic side of MscL. Comparison of simulations in which gating progressed to varying degrees identified residues that pose steric hindrance to channel opening.

Binding Sites↗
PubMed1993

A model for the dynemicin-A cleavage of DNA using molecular dynamics simulation.

Molecular dynamics (MD) simulations to model possible reaction pathways of the dynemicin-A-DNA cleavage mechanism were performed. Two base-pairs sequences, ApCpTpApCpTp-3'/TpGpApTpGpAp-5' and CpApCpGpGpGp-3'/GpTpGpCpCpCp-5', were considered in the calculations. A model based on a prior study of intercalation of dynemicin-A and posterior activation of the drug was assumed in this study. The minimum energy minor groove intercalation complexes for dynemicin-A were used as starting structures in the MD simulations for the reactive intermediate species involved in the postulated action mechanism. The dynemicin-A diol derivative product of the opening of the epoxy ring was used as a "steric mimic" ligand for the DNA-reactive diaryl intermediate. The calculated changes in the geometry of the intercalation complex, due to the opening of the epoxy ring, correspond to the approach of the postulated intermolecular reaction centers in the intercalation states that are responsible for the highest observed DNA cleavage frequency observed. Conversely, unfavorable reaction geometries were found for the intercalation modes corresponding to the lowest observed DNA cutting frequencies.

Anthraquinones↗
PubMed1998

Tertiary contacts in alpha-lactalbumin at pH 7 and pH 2: a molecular dynamics study.

Molecular dynamics simulations of alpha-lactalbumin were performed under conditions of neutral pH and low pH in order to study the acid-induced molten globule state. Through the use of experimental techniques such as NMR and CD spectroscopy, molten globules have been characterized as being compact intermediates with secondary structure similar to that of the native protein but with tertiary structure that is disordered. The detailed structure of the molten globule state is unknown, however. Through the use of computer simulations we can study the structural changes which occur upon lowering pH. The simulations presented here differ from previous unfolding simulations in two important ways: the electrostatic interactions are treated more accurately than ever before, and artificially high temperatures are not used to force the protein to unfold. Simulations of 880 psec each were run at pH 7 (control simulation) and pH 2. We concentrate on the interesting changes in the tertiary interactions within the protein with lowering of pH. In particular, there is a loss of native tertiary contacts in the beta domain and interdomain region, and a large decrease in interdomain hydrogen bonds.

Computer Simulation↗
PubMed2001

Role of the helix capping in the stability of the mouse prion (180-213) segment: investigation through molecular dynamics simulations.

Molecular dynamics simulation of the 180-213 segment, forming the B and C helices in the mouse prion protein, and of three mutants, where the capping box residues or the hydrophobic staple motif residues were selectively mutated, have been carried out. The results indicate that the wild type segment is stable over all the trajectory, whilst the mutants display different degrees of destabilization. In detail mutation of Asp202 brings to a rapid unfolding of helix C likely because of the concomitant loss of a hydrogen bond and of a negative charge able to stabilize the dipole in the first turn of the helix. A lower destabilizing effect is observed upon mutation Thr199. On the other hand mutation of Phe198 and Val203, the hydrophobic staple residues, brings to an incorrect orientation of the first helix relative to the second one due to a weakening of the hydrophobic interaction. The results confirm the importance of the presence of both motifs for the structural integrity of the isolated fragment and suggest that these residues may have a main role in the structural transition observed in the inherited human prion diseases.

Amino Acid Sequence↗
PubMed2005

Adsorption and diffusion of benzene in the nanoporous catalysts FAU, ZSM-5 and MCM-22: a molecular dynamics study.

Molecular dynamics (MD) simulations of benzene in siliceous zeolites (FAU, ZSM-5, and MCM-22) were performed at loadings of 1, 2, 4, 8, and 16 molecules per supercell. The potential energy functions for these simulations were constructed in a semi-empirical way from existing potentials and experimental energetic data. The MD simulations were employed to analyze the dynamic properties of the benzene-zeolite systems. The adsorption energies of benzene/siliceous zeolite complexes increase with increasing loading number, due to the intermolecular attraction between benzene molecules. The self-diffusion coefficient of benzene in siliceous zeolites decreases with increasing loading due to the steric hindrance between the sorbates passing each other. From the zeolite-benzene radial distribution functions it was found that the benzene molecules are relatively far from each other, about 5.2A for siliceous FAU, 5.2A for siliceous ZSM-5, and 4.8A for siliceous MCM-22. In the case of FAU, the benzene molecules prefer to be adsorbed parallel to the surface of the sodalite cage above the six-membered-ring. In ZSM-5, we found a T-structure of the benzene molecules at loadings 2, 4, and 8 molecules per supercell. At loadings of 16 molecules per supercell, the molecules are lined up along the straight channel and their movement is highly correlated. For MCM-22 we found adjacent benzene molecules at a loading of 4 molecules with an orientation similar to the stacked conformation of benzene dimer in the gas phase.

Adsorption↗