SEARCH · Search PubMed

Results for “Molecular Dynamics 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 19 recordsLinked to original sources

PubMed2002

Reaction mechanism of soluble epoxide hydrolase: insights from molecular dynamics simulations.

Molecular dynamics simulations have been performed to gain insights into the catalytic mechanism of the hydrolysis of epoxides to vicinal diols by soluble epoxide hydrolase (sEH). The binding of a substrate, 1S,2S-trans-methylstyrene oxide, was studied in two conformations in the active site of the enzyme. It was found that only one is likely to be found in the active enzyme. In the preferred conformation the phenyl group of the substrate is pi-sandwiched between two aromatic residues, Tyr381 and His523, whereas the other conformation is pi-stacked with only one aromatic residue, Trp334. Two simulations were carried out to 1 ns for each conformation to evaluate the protonation state of active site residue His523. It was found that a protonated histidine is essential for keeping the active site from being disrupted. Long time scale, 4 ns, molecular dynamics simulation was done for the structure with the most likely combination of binding conformation and protonation state of His523. Near Attack Conformers (NACs) are present 5.3% of the time and nucleophilic attack on either epoxide carbon atom, approximately 75% on C(1) and approximately 25% on C(2), is found. A maximum of one hydrogen bond between the epoxide oxygen and either of the active site tyrosines, Tyr465 and Tyr381, is present, in agreement with experimental mutagenesis results that reveal a slight loss in activity if one tyrosine is mutated and essential loss of all activity upon double mutation of the two tyrosines in question. It was found that a hydrogen bond from Tyr465 to the substrate oxygen is essential for controlling the regioselectivity of the reaction. Furthermore, a relationship between the presence of this hydrogen bond and the separation of reactants was found. Two groups of amino acid segments were identified each as moving collectively. Furthermore, an overall anti-correlation was found between the movements of these two individually collectively moving groups, made up by parts of the cap-region, including the two tyrosines, and the site of the catalytic triad, respectively. This overall anti-correlated collective domain motion is, perhaps, involved in the conversion of E.NAC to E.TS.

Animals↗

PubMed2006

Vibrational spectra of methane clathrate hydrates from molecular dynamics simulation.

Molecular dynamics simulations were performed on methane clathrate hydrates at ambient conditions. Thermal expansion results over the temperature range 60-300 K show that the unit cell volume increases with temperature in agreement with experiment. Power spectra were obtained at 273 K from velocity autocorrelation functions for selected atoms, and normal modes were assigned. The spectra were further classified according to individual atom types, allowing the assignment of contributions from methane molecules located in small and large cages within the structure I unit cell. The symmetric C-H stretch of methane in the small cages occurs at a higher frequency than for methane located in the large cages, with a peak separation of 14 cm(-1). Additionally, we determined that the symmetric C-H stretch in methane gas occurs at the same frequency as methane in the large cages. Results of molecular dynamics simulations indicate the use of power spectra obtained from the velocity autocorrelation function is a reliable method to investigate the vibrational behavior of guest molecules in clathrate hydrates.

Letter↗

PubMed2007

Stereoselectivity by enantiomeric inhibitors of matrix metalloproteinase-8: new insights from molecular dynamics simulations.

Molecular Dynamics simulations in aqueous solution were performed for the matrix metalloproteinase-8 (MMP-8) free catalytic domain and for its complexes with the (R)- and (S)-[1-(4'-methoxybiphenyl-4-sulfonylamino)-2-methylpropyl] phosphonate. The 144-155 loop of the enzyme undergoes a drastic decrease of mobility once complexed with both enantiomers. The two enantiomers induce a different decrease of conformational entropy upon complexation. The higher affinity of the R-enantiomer can be related to the lower loss of conformational entropy accompanying its binding. The differences in the dynamical behavior of the protein induced by the two enantiomers are discussed at molecular level and the mode of binding of the simulated complexes is compared with that previously determined by X-ray crystallography.

Biphenyl Compounds↗

PubMed2000

Dynamics of the Hck-SH3 domain: comparison of experiment with multiple molecular dynamics simulations.

Molecular dynamics calculations provide a method by which the dynamic properties of molecules can be explored over timescales and at a level of detail that cannot be obtained experimentally from NMR or X-ray analyses. Recent work (Philippopoulos M, Mandel AM, Palmer AG III, Lim C, 1997, Proteins 28:481-493) has indicated that the accuracy of these simulations is high, as measured by the correspondence of parameters extracted from these calculations to those determined through experimental means. Here, we investigate the dynamic behavior of the Src homology 3 (SH3) domain of hematopoietic cell kinase (Hck) via 5N backbone relaxation NMR studies and a set of four independent 4 ns solvated molecular dynamics calculations. We also find that molecular dynamics simulations accurately reproduce fast motion dynamics as estimated from generalized order parameter (S2) analysis for regions of the protein that have experimentally well-defined coordinates (i.e., stable secondary structural elements). However, for regions where the coordinates are not well defined, as indicated by high local root-mean-square deviations among NMR-determined structural family members or high B-factors/low electron density in X-ray crystallography determined structures, the parameters calculated from a short to moderate length (less than 5-10 ns) molecular dynamics trajectory are dependent on the particular coordinates chosen as a starting point for the simulation.

Crystallography, X-Ray↗

PubMed2006

Understanding the population, coordination, and orientation of water species contributing to the nonlinear optical spectroscopy of the vapor-water interface through molecular dynamics simulations.

Molecular dynamics simulations are used to deconvolve the vibrational spectral features of the vapor-water interface based on molecular environment. A simple geometric description of hydrogen bonding is deployed to identify the OH stretch modes that comprise the vibrational sum-frequency spectrum of the vapor-water interface with direct comparison to our experimental results. The population densities of different species of water molecules are presented as functions of interfacial depth and orientation. It is found that surface water molecules that possess one proton donor bond and one proton acceptor bond make the dominant contribution to both the SSP- and SPS-polarized spectral responses and are located within an angstrom of the Gibbs dividing surface.

Journal Article↗

PubMed2004

Intramolecular diffusive motion in alkane monolayers studied by high-resolution quasielastic neutron scattering and molecular dynamics simulations.

Molecular dynamics simulations of a tetracosane (n-C24H50) monolayer adsorbed on a graphite basal-plane surface show that there are diffusive motions associated with the creation and annihilation of gauche defects occurring on a time scale of approximately 0.1-4 ns. We present evidence that these relatively slow motions are observable by high-energy-resolution quasielastic neutron scattering (QNS) thus demonstrating QNS as a technique, complementary to nuclear magnetic resonance, for studying conformational dynamics on a nanosecond time scale in molecular monolayers.

Journal Article↗

PubMed2006

Spin-labeled alkylphospholipids in a dipalmitoylphosphatidylcholine bilayer: molecular dynamics simulations.

Molecular dynamics simulation has been performed to investigate the structural properties of perifosine and its synthetic spin-labeled alkylphospholipid analogues. The conformations adopted by these compounds in water and in a dipalmitoylphosphatidylcholine bilayer as a function of the presence and position of the N-oxyl-4',4'-dimethyloxazolidine ring (doxyl group) have been investigated by all-atom molecular dynamics. No predominant conformation was observed in water, but the molecules adopt specific orientations and conformations in the lipid bilayer. As is expected, alkyl chains tend to insert into the hydrophobic core, while charged groups stay at the lipid-water interface. A doxyl group in the middle of the alkyl chain moves up to the interface region, thus preventing adoption of the extended conformation. Compounds with a doxyl group close to the polar head group adopt conformations similar to that of unlabeled perifosine within the first nanoseconds of simulation. When the doxyl group is at the end of alkyl chain, the spin-labeled molecule needs more time to reach equilibrium. These results indicate a considerable effect of the doxyl position within the alkyl chain on its localization in the lipid bilayer and can be extended further to other similar spin probes used in the electron paramagnetic resonance spectroscopy of biological membranes.

1,2-Dipalmitoylphosphatidylcholine↗

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↗

PubMed2001

Conformations of an adenine bulge in a DNA octamer and its influence on DNA structure from molecular dynamics simulations.

Molecular dynamics simulations have been applied to the DNA octamer d(GCGCA-GAAC). d(GTTCGCGC), which has an adenine bulge at the center to determine the pathway for interconversion between the stacked and extended forms. These forms are known to be important in the molecular recognition of bulges. From a total of ~35 ns of simulation time with the most recent CHARMM27 force field a variety of distinct conformations and subconformations are found. Stacked and fully looped-out forms are in excellent agreement with experimental data from NMR and x-ray crystallography. Furthermore, in a number of conformations the bulge base associates with the minor groove to varying degrees. Transitions between many of the conformations are observed in the simulations and used to propose a complete transition pathway between the stacked and fully extended conformations. The effect on the surrounding DNA sequence is investigated and biological implications of the accessible conformational space and the suggested transition pathway are discussed, in particular for the interaction of the MS2 replicase operator RNA with its coat protein.

Adenine↗

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↗

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↗

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↗

PubMed2002

Molecular dynamics simulations.

Molecular dynamics simulations have become a standard tool for the investigation of biomolecules. Simulations are performed of ever bigger systems using more realistic boundary conditions and better sampling due to longer sampling times. Recently, realistic simulations of systems as complex as transmembrane channels have become feasible. Simulations aid our understanding of biochemical processes and give a dynamic dimension to structural data; for example, the transformation of harmless prion protein into the disease-causing agent has been modeled.

Cell Membrane↗

PubMed2005

Pressure-dependent transition in protein dynamics at about revealed by molecular dynamics simulation.

Molecular dynamics simulations of a crystalline protein, Staphylococcal nuclease, over the pressure range 1 bar to 15 kbar reveal a qualitative change in the internal protein motions at approximately 4 kbar. This change involves the existence of two linear regimes in the mean-square displacement for internal protein motion, (P) with a twofold decrease in the slope for P>4 kbar. The major effect of pressure on the dynamics is a loss, with increasing pressure of large amplitude, collective protein modes below 2 THz effective frequency, accompanied by restriction of large-scale solvent translational motion.

Computer Simulation↗

PubMed2002

Conformational states of the glucocorticoid receptor DNA-binding domain from molecular dynamics simulations.

Molecular dynamics simulations (MD) have been performed on variant crystal and NMR-derived structures of the glucocorticoid receptor DNA-binding domain (GR DBD). A loop region five residues long, the so-called D-box, exhibits significant flexibility, and transient perturbations of the tetrahedral geometry of two structurally important Cys4 zinc finger are seen, coupled to conformational changes in the D-box. In some cases, one of the Cys ligands to zinc exchanges with water, although no global distortion of the protein structure is observed. Thus, from MD simulation, dynamics of the D-box could partly be explained by solvent effects in conjunction with structural reformation of the zinc finger.

Amino Acid Sequence↗

PubMed2005

Conformational folding of xyloglucan side chains in aqueous solution from molecular dynamics simulation.

Molecular dynamics simulation was carried out on xyloglucan with explicit water molecules to investigate the folding mechanism of side chains onto a main chain in aqueous solution. The model xyloglucan was composed of 12 beta-D-glucopyranoses as a main chain substituted with six galactoses and three xyloses as side chains. Two conditions were set for the ribbon-like main chain; one is restricted to be 'flat' and the other is without restriction. The free main chain of xyloglucan has a 'twisted' conformation as the major one. Conformational folding of side chains onto the main chain was analyzed with dihedral angles at each glycosidic linkage. In a 5-ns calculation, the xyloglucan has a tendency to contract in both the restricted and the free systems, but the mode of contraction is different. Side chains tend to stick onto the flat surface of the main chain in the restricted system, while they do not tightly do so in the free one; instead the main chain takes a twisted and sometimes embowed conformation. This result indicates that the main chain has greater attractive forces to bind side chains when it is flat, while it loses the ability as it is twisted.

Carbohydrate Conformation↗

PubMed2006

Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations.

Molecular dynamics simulations were performed to unfold a homologous pair of thermophilic and mesophilic cold shock proteins at high temperatures. The two proteins differ in just 11 of 66 residues and have very similar structures with a closed five-stranded antiparallel beta-barrel. A long flexible loop connects the N-terminal side of the barrel, formed by three strands (beta1-beta3), with the C-terminal side, formed by two strands (beta4-beta5). The two proteins were found to follow the same unfolding pathway, but with the thermophilic protein showing much slower unfolding. Unfolding started with the melting of C-terminal strands, leading to exposure of the hydrophobic core. Subsequent melting of beta3 and the beta-hairpin formed by the first two strands then resulted in unfolding of the whole protein. The slower unfolding of the thermophilic protein could be attributed to ion pair formation of Arg-3 with Glu-46, Glu-21, and the C-terminal. These ion pairs were also found to be important for the difference in folding stability between the pair of proteins. Thus electrostatic interactions appear to play similar roles in the difference in folding stability and kinetics between the pair of proteins.

Amino Acid Sequence↗

PubMed1992

Cholesterol in model membranes. A molecular dynamics simulation.

Molecular dynamics simulations of a model membrane with inserted cholesterol molecules have been performed to study the perturbing influence of cholesterol. In the fluid phase of a lipid bilayer at 13 mol% concentration of cholesterol, local ordering of the hydrocarbon chains is induced. This perturbation decays with the distance from the cholesterol, and the effect extends 1.25 nm. It can be monitored in several ways, e.g., by an order parameter corresponding to deuterium nuclear magnetic resonance quadrupolar splittings, by the fraction of gauche bonds, or by the local bilayer thickness. At constant surface density, the local ordering is accompanied by disordering of the bulk phase, and, consequently, the net ordering effect is small. After compressing the system laterally in accordance with experimentally known surface areas, the bulk order parameters agree with those of a pure system, and the average order parameters are in accordance with experimental data. The necessity for this lateral compression is supported by calculated lateral pressures. At lower cholesterol concentration (3%), no direct perturbing effect is observed. A smaller lateral pressure than in a pure system indicates that the system with cholesterol is expected to have a smaller surface area, which would result in an increase of the order parameters, thus accounting for the experimental observations. The lack of spatial variation is, however, puzzling and may indicate a cooperative ordering effect.

Cholesterol↗