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 577 records · Page 32Linked to original sources
PubMed2003

Molecular dynamics simulations of a helicase.

Helicases are ubiquitous enzymes involved in nucleic acid metabolism. The PcrA DNA helicase is an essential bacterial protein involved in rolling circle plasmid replication and DNA repair. Recent crystal structures of PcrA bound to DNA indicate that a flexible loop mediates a functionally important rigid-body-domain rotation. In this study, we report stochastic boundary molecular dynamics simulations focused on this region for wild-type and mutants designed to increase the rigidity of the region. Residues in the region that were helix-disfavoring, such as glycine, threonine, and others, were mutated to alanine. The simulated dynamics, analyzed with a variety of measures of structure and mobility, indicate that a few point mutations will substantially increase helix formation in this region. Subnanosecond stochastic boundary molecular dynamics simulations at several temperatures offer a rapid protocol for assessing large numbers of mutants and provides a novel strategy for the design of experiments to test the role of this flexible loop region in the function of PcrA.

Amino Acid Sequence↗
PubMed2006

Insights into correlated motions and long-range interactions in CheY derived from molecular dynamics simulations.

CheY is a response regulator protein involved in bacterial chemotaxis. Much is known about its active and inactive conformations, but little is known about the mechanisms underlying long-range interactions or correlated motions. To investigate these events, molecular dynamics simulations were performed on the unphosphorylated, inactive structure from Salmonella typhimurium and the CheY-BeF(3)(-) active mimic structure (with BeF(3)(-) removed) from Escherichia coli. Simulations utilized both sequences in each conformation to discriminate sequence- and structure-specific behavior. The previously identified conformational differences between the inactive and active conformations of the strand-4-helix-4 loop, which are present in these simulations, arise from the structural, and not the sequence, differences. The simulations identify previously unreported structure-specific flexibility features in this loop and sequence-specific flexibility features in other regions of the protein. Both structure- and sequence-specific long-range interactions are observed in the active and inactive ensembles. In the inactive ensemble, two distinct mechanisms based on Thr-87 or Ile-95 rotameric forms, are observed for the previously identified g+ and g- rotamer sampling by Tyr-106. These molecular dynamics simulations have thus identified both sequence- and structure-specific differences in flexibility, long-range interactions, and rotameric form of key residues. Potential biological consequences of differential flexibility and long-range correlated motion are discussed.

Amino Acid Sequence↗
PubMed2002

Peptide loop-closure kinetics from microsecond molecular dynamics simulations in explicit solvent.

End-to-end contact formation rates of several peptides were recently measured by tryptophan triplet quenching (Lapidus et al. Proc. Natl. Acad. Sci. U.S.A. 2000, 97, 7220). Motivated by these experiments, we study loop-closure kinetics for two peptides of different lengths, Cys-(Ala-Gly-Gln)n-Trp (n = 1, 2), in multiple all-atom explicit-solvent molecular dynamics simulations with different initial conditions and force fields. In 150 simulations of approximately 20 ns each, we collect data covering 1.0 and 0.8 micros for the penta-peptide simulated with the AMBER and CHARMM force fields, respectively, and about 0.5 micros each with the two force fields for the octa-peptide. These extensive simulations allow us to analyze the dynamics of peptides in the unfolded state with atomic resolution, thus probing early events in protein folding, and to compare molecular dynamics simulations directly with experiment. The calculated lifetimes of the tryptophan triplet state are in the range of 50-100 ns, in agreement with experimental measurements. However, end-to-end contacts form more rapidly, with characteristic times less than 10 ns. The contact formation rates for the two force fields are similar despite differences in the respective ensembles of peptide conformations.

Computer Simulation↗
PubMed2001

Molecular dynamics simulation of a lipid diamond cubic phase.

This paper presents the first atomistic simulation of a cubic membrane phase. Using the molecular dynamics simulation technique both the global and the local organization of glycerolmonoolein molecules inside the diamond cubic phase are studied. Multinanosecond simulations reveal that the center of the cubic bilayer remains close to the infinite periodic minimal surface that describes the diamond geometry. We further show that the equilibrium structure of the surfactant molecules inside the cubic phase is very similar to their structure inside a simulated lamellar bilayer. The small differences arise from the packing constraints of the surfactants within the cubic phase which has an area per surfactant that increases toward the bilayer center.

Computer Simulation↗
PubMed2006

Prediction of the binding mode between BMS-378806 and HIV-1 gp120 by docking and molecular dynamics simulation.

BMS-378806 is a newly discovered small molecule that effectively blocks the binding of CD4 with gp120. The binding mode of this kind of inhibitor remains unknown. In this paper, AutoDock 3.0 in conjunction with molecular dynamics simulation, accommodating the receptor's flexibility, was used to explore the binding mode between BMS-378806 and gp120. Two structures, Mode I and Mode II, with the lowest docking energy were selected as different representative binding modes. The analysis of the results from the molecular dynamics simulation indicated that the binding of BMS-348806 in Mode II is more stable. The average structure of Mode II was analyzed and compared with the experimental data. The conclusion was that BMS-378806 inserts the azaindole ring deeply into the PHE43 cavity and makes contact with a number of residues in the cavity, on the cavity and near the cavity. This study benefits the understanding of the mechanism of this kind of inhibitor and may provide useful information for rational drug design.

Binding Sites↗
PubMed2005

Ion transport through membrane-spanning nanopores studied by molecular dynamics simulations and continuum electrostatics calculations.

Narrow hydrophobic regions are a common feature of biological channels, with possible roles in ion-channel gating. We study the principles that govern ion transport through narrow hydrophobic membrane pores by molecular dynamics simulation of model membranes formed of hexagonally packed carbon nanotubes. We focus on the factors that determine the energetics of ion translocation through such nonpolar nanopores and compare the resulting free-energy barriers for pores with different diameters corresponding to the gating regions in closed and open forms of potassium channels. Our model system also allows us to compare the results from molecular dynamics simulations directly to continuum electrostatics calculations. Both simulations and continuum calculations show that subnanometer wide pores pose a huge free-energy barrier for ions, but a small increase in the pore diameter to approximately 1 nm nearly eliminates that barrier. We also find that in those wider channels the ion mobility is comparable to that in the bulk phase. By calculating local electrostatic potentials, we show that the long range Coulomb interactions of ions are strongly screened in the wide water-filled channels. Whereas continuum calculations capture the overall energetics reasonably well, the local water structure, which is not accounted for in this model, leads to interesting effects such as the preference of hydrated ions to move along the pore wall rather than through the center of the pore.

Computer Simulation↗
PubMed2006

Molecular dynamics simulations of the complete satellite tobacco mosaic virus.

This work presents an all-atom molecular dynamics simulation of a complete virus, the satellite tobacco mosaic virus. Simulations with up to 1 million atoms for over 50 ns demonstrate the stability of the entire virion and of the RNA core alone, while the capsid without RNA exhibits a pronounced instability. Physical properties of the simulated virus particle including electrostatic potential, radial distribution of viral components, and patterns of correlated motion are analyzed, and the implications for the assembly and infection mechanism of the virus are discussed.

Capsid↗
PubMed2004

Molecular dynamics simulations of a guaiacyl beta-O-4 lignin model compound: examination of intramolecular hydrogen bonding and conformational flexibility.

The dynamical conformational behavior of a guaiacyl beta-O-4 lignin model compound has been investigated by molecular simulations. The potential energy surface of the molecule in vacuum has been examined by means of an adiabatic map, showing a large accessible conformational space with multiple energy minima separated by low barriers. Molecular dynamics simulations have been performed in vacuum and with explicit solvent molecules for 10 and 2.1 ns, respectively. Molecular dynamics trajectories recorded in vacuum have shown the molecule to be flexible and to visit a large number of conformations. Many intramolecular H-bonds have been observed, existing for more than 90% of the total simulation time. The presence of explicit solvent molecules induces a significant broadening of some regions of the accessible conformational space and also largely reduces the statistical significance of intramolecular H-bonding. Intramolecular H-bonds observed in vacuum do not persist significantly and are preferentially exchanged with intermolecular H-bonds to the surrounding solvent molecules. The theoretical results are in good agreement with experimental NMR data that do not support the existence of strong and persistent intramolecular H-bonds in solution but instead indicate that H-bonds to solvent predominate. Finally, both molecular modeling and NMR approaches predict the guaiacyl beta-O-4 structure to be flexible and indicate that intramolecular H-bonds are not strong and persistent enough to confer rigidity to the molecule in solution.

Guaiacol↗
PubMed2002

Molecular dynamics simulation of the response of a gas to a spherical piston: implications for sonoluminescence.

Sonoluminescence is the phenomena of light emission from a collapsing gas bubble in a liquid. Theoretical explanations of this extreme energy focusing are controversial and difficult to validate experimentally. We propose to use molecular dynamics simulations of the collapsing gas bubble to clarify the energy focusing mechanism, and determine physical parameters that restrict theories of the light emitting mechanism. In this paper, we model the interior of a collapsing noble gas bubble as a hard sphere gas driven by a spherical piston boundary moving according to the Rayleigh-Plesset equation. We also include a simplified treatment of ionization effects in the gas at high temperatures. The effects of water vapor are neglected in the model. By using fast, tree-based algorithms, we can exactly follow the dynamics of 10(6) particle systems during the collapse. Our preliminary model shows strong energy focusing within the bubble, including the formation of shocks, strong ionization, and temperatures in the range of 50 000-500 000 K. Our calculations show that the gas-liquid boundary interaction has a strong effect on the internal gas dynamics, and that the gas passes through states where the mean free path is greater than the characteristic distance over which the temperature varies. We also estimate the duration of the light pulse from our model, which predicts that it scales linearly with the ambient bubble radius. As the number of particles in a physical sonoluminescing bubble is within the foreseeable capability of molecular dynamics simulations, we also propose that fine scale sonoluminescence experiments can be viewed as excellent test problems for advancing the art of molecular dynamics.

Journal Article↗
PubMed1995

Estimation of binding free energies for HIV proteinase inhibitors by molecular dynamics simulations.

Absolute binding free energies for three inhibitors of HIV-1 proteinase were estimated from molecular dynamics simulations by a recently reported linear approximation procedure. The results were in fairly good agreement with experimental binding data. Two of the inhibitors were very similar and, for comparison, their relative free energies of binding were also calculated by free energy perturbation methods, giving virtually the same result. Effects of cut-off radii and charge states of the protein model were examined. The effects of pH on binding of one of the inhibitors were predicted.

Computer Simulation↗
PubMed1995

Molecular dynamics simulations indicate that F87W,T185F-cytochrome P450cam may reductively dehalogenate 1,1,1-trichloroethane.

Cytochrome P450cam is capable of reductively dehalogenating several chlorinated alkanes at low, but measurable, rates. In previous investigations of structure-function relationships in this enzyme using molecular dynamics simulations, we noticed that 1,1,1-trichloroethane (TCA) exhibits a very high degree of mobility in the active site due to its smaller molecular volume relative to the native substrate, camphor(1,2). Several amino acid sidechains lining the active site also exhibit significant dynamic fluctuations, possibly as a result of poor steric complementary to TCA. Guided by these results, we modeled double (F87W, T185F) and triple (F87W, T185F, V295I) mutants of P450cam, which provide additional bulk in the active site and increase the frequency of heme-substrate collision. Molecular dynamics simulations (300 ps on each protein) indicate that these mutants do not significantly perturb the three-dimensional fold of the enzyme, or local structure in the region of the active site. Both mutants bind the substrate more stably near the heme than the wild-type. Interestingly, however, the bulkier triple mutant seems to actually inhibit heme-substrate interactions relative to the double mutant. Over the final 200 ps of simulation, TCA is within 1 A of nonbonded contact with the heme 25% more often in the double mutant versus the wild-type. The triple mutant, on the other hand, binds TCA within 1 A of the heme only 15% as often as the wild-type. These results indicate that the double mutant may reductively dehalogenate TCA, a property not observed for the native protein. Implications for other experimentally measurable parameters are discussed.

Binding Sites↗
PubMed2005

Molecular dynamics simulations of human butyrylcholinesterase.

Herein, we present results from molecular dynamics (MD) simulations of the human butyrylcholinesterase (BuChE) enzyme in aqueous solution. Two configurations of the unbound form of BuChE differing in the presence or absence of a sodium ion inside the protein gorge were simulated for 10 and 5 ns, respectively. Besides complementing the structural information provided by X-ray data, the MD simulations give insight into the structure of the native BuChE enzyme. For example, it is shown that: the nucleophilic Ser(198) residue and the various binding subsites in the BuChE catalytic cavity are readily accessible from the exterior of the protein; the presence of the sodium ion dynamically explores two different binding sites in the gorge leading to the active site and stabilizes the productive conformation of the Glu(325)/His(438)/Ser(198) catalytic triad; several long-lived water bridges are fully integrated into the architecture of the active site; the positions of the residues at the rim of the gorge region display large deviations with respect to the crystal structure; and two side doors, constituted by residues situated at the tip of the acyl- and Omega-loops, respectively, open wide enough to allow the passage of water molecules. In conclusion, we compare our theoretical results with those from previous work on mouse acetylcholinesterase and discuss their implications for substrate binding and catalysis in BuChE.

Animals↗
PubMed1995

Thermodynamic parameters for the helix-coil transition of oligopeptides: molecular dynamics simulation with the peptide growth method.

The helix-coil transition equilibrium of polypeptides in aqueous solution was studied by molecular dynamics simulation. The peptide growth simulation method was introduced to generate dynamic models of polypeptide chains in a statistical (random) coil or an alpha-helical conformation. The key element of this method is to build up a polypeptide chain during the course of a molecular transformation simulation, successively adding whole amino acid residues to the chain in a predefined conformation state (e.g., alpha-helical or statistical coil). Thus, oligopeptides of the same length and composition, but having different conformations, can be incrementally grown from a common precursor, and their relative conformational free energies can be calculated as the difference between the free energies for growing the individual peptides. This affords a straightforward calculation of the Zimm-Bragg sigma and s parameters for helix initiation and helix growth. The calculated sigma and s parameters for the polyalanine alpha-helix are in reasonable agreement with the experimental measurements. The peptide growth simulation method is an effective way to study quantitatively the thermodynamics of local protein folding.

Alanine↗
PubMed1991

Molecular dynamics simulation in vacuo and in solution of cyclolinopeptide A: a conformational study.

The conformation of cyclolinopeptide A [c-(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val)], a naturally occurring peptide with remarkable cytoprotective activity, has been investigated by means of molecular dynamics simulations in various molecular environments. Structural and dynamical properties have been analyzed and compared with those experimentally determined. A detailed analysis of hydrogen bonds is reported.

Amino Acid Sequence↗
PubMed2005

Examination of dynamic facilitation in molecular dynamics simulations of glass-forming liquids.

Using data from molecular dynamics computer simulations of the one-component Dzugutov liquid and of BKS silica in metastable equilibrium supercooled states, we examine ideas introduced by Garrahan and Chandler (GC) in their dynamic facilitation (DF) model of the glass transition. Utilizing a recently introduced measure of DF, we find that DF is important for particle motion in both the supercooled Dzugutov liquid and in the BKS silica melt, that mobility propagates continuously, and that this effect becomes increasingly pronounced with decreasing T. We show that, in both systems, dynamic facilitation is strongest on the time scale of the late-beta relaxation, where clusters of highly mobile neighboring particles escaping from their cages are largest and, except for the silicon atoms in BKS silica, stringlike motion is most prominent. By comparing the two systems, we show that the temperature dependence of one measure of DF as the mode-coupling temperature is approached from high temperature is similar, once the temperature dependence of the structural relaxation time in each system is scaled out.

Journal Article↗
PubMed2006

B-DNA under stress: over- and untwisting of DNA during molecular dynamics simulations.

The twist flexibility of DNA is central to its many biological functions. Explicit solvent molecular dynamics simulations in combination with an umbrella sampling restraining potential have been employed to study induced twist deformations in DNA. Simulations allowed us to extract free energy profiles for twist deformations and were performed on six DNA dodecamer duplexes to cover all 10 possible DNA basepair steps. The shape of the free energy curves was similar for all duplexes. The calculated twist deformability was in good agreement with experiment and showed only modest variation for the complete duplexes. However, the response of the various basepair steps on twist stress was highly nonuniform. In particular, pyrimidine/purine steps were much more flexible than purine/purine steps followed by purine/pyrimidine steps. It was also possible to extract correlations of twist changes and other helical as well as global parameters of the DNA molecules. Twist deformations were found to significantly alter the local as well as global shape of the DNA modulating the accessibility for proteins and other ligands. Severe untwisting of DNA below an average of 25 degrees per basepair step resulted in the onset of a global structural transition with a significantly smaller twist at one end of the DNA compared to the other.

Base Pairing↗
PubMed2005

Water molecules and hydrogen-bonded networks in bacteriorhodopsin--molecular dynamics simulations of the ground state and the M-intermediate.

Protein crystallography provides the structure of a protein, averaged over all elementary cells during data collection time. Thus, it has only a limited access to diffusive processes. This article demonstrates how molecular dynamics simulations can elucidate structure-function relationships in bacteriorhodopsin (bR) involving water molecules. The spatial distribution of water molecules and their corresponding hydrogen-bonded networks inside bR in its ground state (G) and late M intermediate conformations were investigated by molecular dynamics simulations. The simulations reveal a much higher average number of internal water molecules per monomer (28 in the G and 36 in the M) than observed in crystal structures (18 and 22, respectively). We found nine water molecules trapped and 19 diffusive inside the G-monomer, and 13 trapped and 23 diffusive inside the M-monomer. The exchange of a set of diffusive internal water molecules follows an exponential decay with a 1/e time in the order of 340 ps for the G state and 460 ps for the M state. The average residence time of a diffusive water molecule inside the protein is approximately 95 ps for the G state and 110 ps for the M state. We have used the Grotthuss model to describe the possible proton transport through the hydrogen-bonded networks inside the protein, which is built up in the picosecond-to-nanosecond time domains. Comparing the water distribution and hydrogen-bonded networks of the two different states, we suggest possible pathways for proton hopping and water movement inside bR.

Bacteriorhodopsins↗
PubMed2003

Molecular dynamics simulations of the transmembrane domain of the oncogenic ErbB2 receptor dimer in a DMPC bilayer.

Molecular dynamics simulations of an atomic model of the transmembrane domain of the oncogenic ErbB2 receptor dimer embedded in an explicit dimyristoylphosphatidylcholine (DMPC) bilayer were performed for more than 4 ns. The oncogenic Glu mutation in the membrane spanning segment plays a major role in tyrosine kinase activity and receptor dimerization, and is thought to be partly responsible for the structure of the transmembrane domain of the active receptor. MD results show that the interactions between the two transmembrane helices are characteristic of a left-handed packing as previously demonstrated from in vacuo simulations. Moreover, MD results reveal the absence of persistent hydrogen bonds between the Glu side chains in a membrane environment, which raise the question of the ability for Glu alone to stabilize the TM domain of the ErbB2 receptor. Interestingly the formation of the alpha-pi motif in the two ErbB2 transmembrane helices confirms the concept of intrinsic sequence-induced conformational flexibility. From a careful analysis of our MD results, we suggest that the left-handed helix-helix packing could be the key to correctly orient the intracellular domain of the activated receptor dimer. The prediction of such interactions from computer simulations represents a new step towards the understanding of signaling mechanisms.

Computer Simulation↗