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

How does huperzine A enter and leave the binding gorge of acetylcholinesterase? Steered molecular dynamics simulations.

The entering and leaving processes of Huperzine A (HupA) binding with the long active-site gorge of Torpedo californica acetylcholinesterase (TcAChE) have been investigated by using steered molecular dynamics simulations. The analysis of the force required along the pathway shows that it is easier for HupA to bind to the active site of AChE than to disassociate from it, which for the first time interprets at the atomic level the previous experimental result that unbinding process of HupA is much slower than its binding process to AChE. The direct hydrogen bonds, water bridges, and hydrophobic interactions were analyzed during two steered molecular dynamics (SMD) simulations. Break of the direct hydrogen bond needs a great pulling force. The steric hindrance of bottleneck might be the most important factor to produce the maximal rupture force for HupA to leave the binding site but it has a little effect on the binding process of HupA with AChE. Residue Asp72 forms a lot of water bridges with HupA leaving and entering the AChE binding gorge, acting as a clamp to take out HupA from or put HupA into the active site. The flip of the peptide bond between Gly117 and Gly118 has been detected during both the conventional MD and SMD simulations. The simulation results indicate that this flip phenomenon could be an intrinsic property of AChE and the Gly117-Gly118 peptide bond in both HupA bound and unbound AChE structures tends to adopt the native enzyme structure. At last, in a vacuum the rupture force is increased up to 1500 pN while in water solution the greatest rupture force is about 800 pN, which means water molecules in the binding gorge act as lubricant to facilitate HupA entering or leaving the binding gorge.

Acetylcholinesterase↗
PubMed1991

Molecular dynamics simulations of the unfolding of an alpha-helical analogue of ribonuclease A S-peptide in water.

Molecular dynamics simulations of the S-peptide analogue AETAAAKFLREHMDS have been conducted in aqueous solution for 300 ps at 278 K and for 500 ps in two different runs at 358 K. The results show agreement with experimental observations in that at low temperature, 5 degrees C, the helix is stable, while unfolding is observed at 85 degrees C. In the low-temperature simulation a solvent-separated ion pair was formed between Glu-2 and Arg-10, and the side chain of His-12 reoriented toward the C-terminal end of the alpha-helix. Detailed analyses of the unfolding pathways at high temperature have also revealed that the formation or disappearance of main-chain helical hydrogen bonds occurs frequently through an alpha in equilibrium with 3(10) in equilibrium with no hydrogen bond sequence.

Amino Acid Sequence↗
PubMed2002

Steered molecular dynamics simulations on the "tail helix latch" hypothesis in the gelsolin activation process.

The molecular basis of the "tail helix latch" hypothesis in the gelsolin activation process has been studied by using the steered molecular dynamics simulations. In the present nanosecond scale simulations, the tail helix of gelsolin was pulled away from the S2 binding surface, and the required forces were calculated, from which the properties of binding between the tail helix and S2 domain and their dynamic unbinding processes were obtained. The force profile provides a detailed rupture mechanism that includes six major unbinding steps. In particular, the hydrogen bonds formed between Arg-207 and Asp-744 and between Arg-221 and Leu-753 are of the most important interaction pairs. The two hydrogen bond "clamps" stabilize the complex. The subsequent simulation on Arg-207-Ala (R207A) mutation of gelsolin indicated that this mutation facilitates the unbinding of the tail helix and that the contribution of the hydrogen bond between Arg-207 and Asp-744 to the binding is more than 50%, which offers a new clue for further mutagenesis study on the activation mechanism of gelsolin. Surrounding water molecules enhance the stability of the tail helix and facilitate the rupture process. Additionally, temperature also has a significant effect on the conformation of the arginine and arginine-related interactions, which revealed the molecular basis of the temperature dependence in gelsolin activation.

Arginine↗
PubMed2005

Molecular dynamics simulation of dense carbon dioxide fluid on amorphous silica surfaces.

Molecular dynamics (MD) simulations of dense carbon dioxide on the amorphous dehydroxylated silica surfaces have been carried out. The adsorption potential surfaces of the silica solids have been obtained in order to evaluate the characteristics of the amorphous surfaces. The atom density profiles, adsorption free energy profiles, surface orientation order parameters, and radial distribution functions for the CO2 molecules have been presented in order to study the effect of the amorphous surfaces on the microscopic interfacial structure properties of the CO2 molecules. The translational diffusion and orientation rotation at silica surfaces have also been investigated. It was observed that there is marked hindrance of the translational diffusion and orientation rotation of CO2 molecules near amorphous silica surfaces.

Journal Article↗
PubMed2006

Molecular dynamics simulation of liquid sulfur dioxide.

A previously proposed model for molecular dynamics (MD) simulation of liquid sulfur dioxide, SO(2), has been reviewed. Thermodynamic, structural, and dynamical properties were calculated for a large range of thermodynamic states. Predicted (P,V,T) of simulated system agrees with an elaborated equation of state recently proposed for liquid SO(2). Calculated heat capacity, expansion coefficient, and isothermal compressibility are also in good agreement with experimental data. Calculated equilibrium structure agrees with X-ray and neutron scattering measurements on liquid SO(2). The model also predicts the same (SO(2))(2) dimer structure as previously determined by ab initio calculations. Detailed analysis of equilibrium structure of liquid SO(2) is provided, indicating that, despite the rather large dipole moment of the SO(2) molecule, the structure is mainly determined by the Lennard-Jones interactions. Both single-particle and collective dynamics are investigated. Temperature dependency of dynamical properties is given. The MD results are compared with previous findings obtained from the analysis of inelastic neutron scattering spectra of liquid SO(2), including wave-vector dependent structural relaxation, tau(k), and viscosity, eta(k).

Computer Simulation↗
PubMed2005

Molecular dynamics simulations of the surface tension of ionic liquids.

We report molecular dynamics computer simulations of the surface tension and interfacial thickness of ionic liquid-vapor interfaces modeled with a soft core primitive model potential. We find that the surface tension shows an anomalous oscillatory behavior with interfacial area. This observation is discussed in terms of finite size effects introduced by the periodic boundary conditions employed in computer simulations. Otherwise we show that the thickness of the liquid-vapor interface increases with surface area as predicted by the capillary wave theory. Data on the surface tension of size-asymmetric ionic liquids are reported and compared with experimental data of molten salts. Our data suggest that the surface tensions of size-asymmetric ionic liquids do not follow a corresponding states law.

Journal Article↗
PubMed2004

"Learn on the fly": a hybrid classical and quantum-mechanical molecular dynamics simulation.

We describe and test a novel molecular dynamics method which combines quantum-mechanical embedding and classical force model optimization into a unified scheme free of the boundary region, and the transferability problems which these techniques, taken separately, involve. The scheme is based on the idea of augmenting a unique, simple parametrized force model by incorporating in it, at run time, the quantum-mechanical information necessary to ensure accurate trajectories. The scheme is tested on a number of silicon systems composed of up to approximately 200 000 atoms.

Journal Article↗
PubMed1999

Molecular dynamics simulations of human glutathione transferase P1-1: conformational fluctuations of the apo-structure.

We have investigated by molecular dynamics simulations the conformational fluctuations of the monomer of human apo-glutathione transferase P1-1. After attainment of steady-state dynamics, the structural fluctuations involve mainly the protein segments that participate also in the holo-apo transition discussed in the accompanying article (Stella et al., 1999:37:1-9.). The most mobile region is the C-terminal segment of helix 2. In contrast, helices 1, 6, 7, and 8 constitute a relatively rigid protein core. An "essential dynamics" analysis of the simulation shows that the largest fluctuations involve specific regions of glutathione transferases. In such regions, atomic motions are correlated. Motions of helix 2 are accounted for by the second most prominent principal component, which reveals a fluctuation between two distinct conformations. The residues that constitute the H-site undergo a breathing motion, possibly relevant during the binding of hydrophobic cosubstrates. Based on our simulation, several experimental findings can be rationalized, including the viscosity-dependent reactivity of Cys 47 and Cys 101 as well as the selective proteolysis of the peptide bond between Lys 44 and Ala 45. We have also modeled the structural changes that lead to the formation of an intrachain disulfide bridge between cysteines 47 and 101 and to the inactivation of the enzyme. The resulting structure maintains essentially the native fold except for helix 2, which closes the G-site. Proteins 1999;37:10-19.

Apoproteins↗
PubMed2003

Steered molecular dynamics simulation on the binding of NNRTI to HIV-1 RT.

HIV-1 reverse transcriptase (RT) is the primary target for anti-AIDS chemotherapy. Nonnucleoside RT inhibitors (NNRTIs) are very potent and most promising anti-AIDS drugs that specifically inhibit HIV-1 RT. The binding and unbinding processes of alpha-APA, an NNRTI, have been studied using nanosecond conventional molecular dynamics and steered molecular dynamics simulations. The simulation results show that the unbinding process of alpha-APA consists of three phases based on the position of alpha-APA in relation to the entrance of the binding pocket. When alpha-APA is bound in the binding pocket, the hydrophobic interactions between HIV-1 RT and alpha-APA dominate the binding; however, the hydrophilic interactions (both direct and water-bridged hydrogen bonds) also contribute to the stabilizing forces. Whereas Tyr-181 makes significant hydrophobic interactions with alpha-APA, Tyr-188 forms a strong hydrogen bond with the acylamino group (N14) of alpha-APA. These two residues have very flexible side chains and appear to act as two "flexible clamps" discouraging alpha-APA to dissociate from the binding pocket. At the pocket entrance, two relatively inflexible residues, Val-179 and Leu-100, gauge the openness of the entrance and form the bottleneck of the inhibitor-unbinding pathway. Two special water molecules at the pocket entrance appear to play important roles in inhibitor recognition of binding and unbinding. These water molecules form water bridges between the polar groups of the inhibitor and the residues around the entrance, and between the polar groups of the inhibitor themselves. The water-bridged interactions not only induce the inhibitor to adopt an energetically favorable conformation so the inhibitor can pass through the pocket entrance, but also stabilize the binding of the inhibitor in the pocket to prevent the inhibitor's dissociation. The complementary steered molecular dynamics and conventional molecular dynamics simulation results strongly support the hypothesis that NNRTIs inhibit HIV-1 RT polymerization activity by enlarging the DNA-binding cleft and restricting the flexibility and mobility of the p66 thumb subdomain that are believed to be essential during DNA translocation and polymerization.

Acetamides↗
PubMed2001

Liquid-liquid phase transitions of phosphorus via constant-pressure first-principles molecular dynamics simulations.

Pressure-induced phase transitions in liquid phosphorus have been studied by constant-pressure first-principles molecular dynamics simulations. By compressing a low-pressure liquid which consists of the tetrahedral P4 molecules, a structural phase transition from the molecular to polymeric liquid (a high-pressure phase) observed in the recent experiment by Katayama et al. [Nature (London) 403, 170 (2000)] was successfully realized. It is found that this transition is caused by a breakup of the tetrahedral molecules with large volume contraction. The same transition is also realized by heating. This indicates that only the polymeric liquid can stably exist at high temperature.

Journal Article↗
PubMed1992

Differential helix propensity of small apolar side chains studied by molecular dynamics simulations.

A series of oligoalanine molecules with single amino acid replacements in the middle of the chain has been studied by molecular dynamics simulations. Differences in stability of the alpha-helix (as free energies delta delta G degrees) were estimated for the following series of residues: alpha-aminoisobutyric acid, alanine, alpha-amino-n-butyric acid, valine, glycine, D-alanine, t-leucine (= alpha-amino-beta,beta-dimethyl-n- butyric acid), and proline, arranged here in decreasing order of helix-forming potential. (The results for proline and valine had been reported earlier.) No experimental results were available for alpha-amino-n-butyric acid, D-alanine, and t-leucine at the time these calculations were done. The values of delta delta G degrees, including the three predictions, are in striking agreement with recent experimental results. A combination of free dynamics, dynamics with forced conformational change, and dynamics with forced molecular replacement was used. Conformational distributions were calculated for the peptide backbone of the dipeptides and, where appropriate, for the side chains of the dipeptide and the alpha-helix. The results demonstrate an unexpected level of accuracy for the all-atom model used to represent atomic interactions in the simulations. The simulations permit a detailed analysis of different factors responsible for conformational preferences and differences in stability. These conclusions drawn from this analysis agree with accepted qualitative explanations and allow these explanations to be quantitated to an extent not heretofore possible.

Alanine↗
PubMed2005

E. coli RNase HI and the phosphonate-DNA/RNA hybrid: molecular dynamics simulations.

A model for the complex between E. coli RNase HI and the DNA/RNA hybrid (previously refined by molecular dynamics simulations) was used to determine the impact of the internucleotide linkage modifications (either 3-O-CH2-P-O-5' or 3-O-P-CH2-O-5) on the ability of the modified-DNA/RNA hybrid to create a complex with the protein. Modified internucleotide linkages were incorporated systematically at different positions close to the 3-end of the DNA strand to interfere with the DNA binding site of RNase H. Altogether, six trajectories were produced (length 1.5ns). Mutual hydrogen bonds connecting both strands of the nucleic acids hybrid, DNA with RNase H, RNA with RNase H, and the scissile bond with the Mg++. 4H2O chelate complex (bound in the active site) were analyzed in detaiL Many residues were involved in binding of the DNA (Arg88, Asn84, Trp85, Trp104, Tyr73, Lys99, Asn100, Thr43, and Asn 16) and RNA (Gln76, Gln72, Tyr73, Lys122, Glu48, Asn44, and Cys13) strand to the substrate-binding site of the RNase H enzyme. The most remarkable disturbance of the hydrogen bonding net was observed for structures with modified internucleotide linkages positioned in a way to interact with the Trp104, Tyr73, Lys99, and Asn100 residues (situated in the middle of the DNA binding site, where a cluster of Trp residues forms a rigid core of the protein structure).

Binding Sites↗
PubMed1999

Molecular dynamics simulation of a synthetic four-alpha-helix bundle that binds the anesthetic halothane.

The structural features of binding sites for volatile anesthetics are examined by performing a molecular dynamics simulation study of the synthetic four-alpha-helix bundles (Aalpha2)2, which are formed by association of two 62-residue di-alpha-helical peptides. The peptide bundle (Aalpha2)2 was designed by Johansson et al. [Biochemistry 37 (1998) 1421-1429] and was shown experimentally to have a high affinity for the binding of the anesthetic halothane (CF3CBrCIH) in a hydrophobic cavity. Since (Aalpha2)2 can exhibit either the anti or syn topologies, the two distinct bundles are simulated both in the presence and in the absence of halothane. Nanosecond length molecular dynamics trajectories were generated for each system at room temperature (T = 298 K). The structural and dynamic effects of the inclusion of halothane are compared, illustrating that the structures are stable over the course of the simulation; that the (Aalpha2)2 bundles have suitable pockets that can accommodate halothane; that the halothane remains in the designed hydrophobic cavity in close proximity to the Trp residues with a preferred orientation; and that the dimensions of the peptide are perturbed by the inclusion of an anesthetic molecule.

Amino Acid Sequence↗
PubMed2004

Homology modeling and molecular dynamics simulations of transmembrane domain structure of human neuronal nicotinic acetylcholine receptor.

A three-dimensional model of the transmembrane domain of a neuronal-type nicotinic acetylcholine receptor (nAChR), (alpha4)2(beta2)3, was constructed from a homology structure of the muscle-type nAChR recently determined by cryo-electron microscopy. The neuronal channel model was embedded in a fully hydrated DMPC lipid bilayer, and molecular-dynamics simulations were performed for 5 ns. A comparative analysis of the neuronal- versus muscle-type nAChR models revealed many conserved pore-lining residues, but an important difference was found near the periplasmic mouth of the pore. A flickering salt-bridge of alpha4-E266 with its adjacent beta2-K260 was observed in the neuronal-type channel during the course of the molecular-dynamics simulations. The narrowest region, with a pore radius of approximately 2 A formed by the salt-bridges, does not seem to be the restriction site for a continuous water passage. Instead, two hydrophobic rings, formed by alpha4-V259, alpha4-L263, and the homologous residues in the beta2-subunits, act as the gates for water flow, even though the region has a slightly larger pore radius. The model offers new insight into the water transport across the (alpha4)2(beta2)3 nAChR channel, and may lead to a better understanding of the structures, dynamics, and functions of this family of ion channels.

Amino Acid Sequence↗
PubMed2006

Accelerated molecular dynamics simulation of the thermal desorption of n-alkanes from the basal plane of graphite.

We utilize accelerated molecular dynamics to simulate alkane desorption from the basal plane of graphite. Eight different molecules, ranging from n-pentane to n-hexadecane, are studied in the low coverage limit. Acceleration of the molecular dynamics simulations is achieved using two different methods: temperature acceleration and a compensating potential scheme. We find that the activation energy for desorption increases with increasing chain length. The desorption prefactor increases with chain length for molecules ranging from pentane to decane. This increase subsides and the value of the preexponential factor fluctuates about an apparently constant value for decane, dodecane, tetradecane, and hexadecane. These trends are consistent with data obtained in experimental temperature-programed desorption (TPD) studies. We explain the dependence of the preexponential factor on alkane chain length by examining conformational changes within the alkane molecules. For the shorter molecules, torsional motion is not activated over experimental temperature ranges. These molecules can be treated as rigid rods and their partial loss in translational and rotational entropies upon adsorption increases as chain length increases, leading to an increasing preexponential factor. At their typical TPD peak temperatures, torsions are activated in the longer adsorbed chain molecules to a significant extent which increases with increasing chain length, increasing the entropy of the adsorbed molecule. This increase counteracts the decrease in entropy due to a loss of translation and rotation, leading to a virtually constant prefactor.

Journal Article↗
PubMed1995

Acid and thermal denaturation of barnase investigated by molecular dynamics simulations.

The transition in barnase from the native state to a partially unfolded conformation has been studied by molecular dynamics simulations with explicit water molecules at 360 K and low pH(450 ps), and at 600 K and neutral pH (three simulations of 120, 250 and 200 ps each). The use of several simulations provides evidence that the results are not sensitive to initial conditions. To mimic low pH conditions, the acidic sidechains in barnase were neutralized and the two histidine residues were doubly protonated. Runs at 300 K showed that the solvated structures at low pH (300 ps) and neutral pH (310 ps) are very similar. The main structural differences involved the acidic residues, histidine residues, and the beta-turn connecting strands 4 and 5. When the temperature is raised to 360 K at low pH and to 600 K at neutral pH the barnase molecule begins to unfold. The molecule rapidly expands (Rg changes from 13.9 A to 15.3 A in 450 ps at 360 K and from 13.7 A to between 15.1 and 15.5 A in 120 ps at 600 K). However, the expansion is not uniform. In all the simulations, the chain termini, loops and the N-terminal parts of the main alpha-helix (helix 1) show a continuous and progressive unfolding. An essential step in the denaturation process is that the major alpha-helix (helix 1) separates from the beta-sheet; this is coupled to the exposure of the principal hydrophobic core, many of whose non-polar side chains become solvated by hydrogen-bonded water molecules. The barnase-water interaction energy improves during unfolding at the expense of the barnase self-energy. The deterioration of the intramolecular van der Waals energy suggests that the rupture of the tight packing during the initial unfolding phase contributes to the energy barrier of the denaturation process. The mutationally well-analyzed Asp8-Arg110-Asp12 double salt-bridge on the barnase surface is found to be marginally stable in the folded form in the simulations. A Poisson-Boltzmann calculation indicates that the salt-bridge is unstable; this is probably due to an overestimate of the solvation energy. A detailed analysis of the main hydrophobic core reveals that increase in solvent-accessible surface area and penetration of water molecules are simultaneous in the high-temperature simulation; at lower temperatures there is significant cavity formation and the entrance of the water molecules is somewhat delayed. The cavities occur in the neighborhood of the hydrophobic sidechains; the region formed by the sidechains of Val10, Leu14, Leu20, Tyr24, Ala74, Ile76 and Tyr90 is involved. The loosening of the core packing is coupled to an increase in the number of dihedral transitions.(ABSTRACT TRUNCATED AT 400 WORDS)

Bacterial Proteins↗
PubMed2006

Mimicking the action of GroEL in molecular dynamics simulations: application to the refinement of protein structures.

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results support the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Algorithms↗