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 1,603 records · Page 89Linked to original sources
PubMed2005

The role of electrostatic interaction in triggering the unraveling of stable helix 1 in normal prion protein. A molecular dynamics simulation investigation.

The conversion of normal prion protein (PrPC) into scrapie isoform (PrPSc) is a key event in the pathogenesis of prion diseases. However, the conversion mechanism has given rise to much controversy. For instance, there is much debate on the behavior of helix 1 (H1) in the conversion. A series of experiments demonstrated that H1 in isolated state was very stable under a variety of conditions. But, other experiments indicated that helices 2 and 3 rather than H1 were retained in PrPSc. In this paper, molecular dynamics (MD) simulation is employed to investigate the dynamic behavior of H1. It is revealed that although the helix 1 of Human PrPC (HuPrPC) is very stable in the isolated state, it becomes unstable when incorporated into native HuPrPC, which likely results from the long-range electrostatic interaction between Asp147 and Arg208 located in the helices 1 and 3, respectively. This explanation is supported by experimental evaluation and MD simulation on D147N mutant of HuPrPC that the mutant becomes a little more stable than the wild type HuPrPC. This finding not only help to reconcile the existing debate on the role of helix 1 in the PrPC-->PrPSc transition, but also reveals a possible mechanism for triggering the PrPC-->PrPSc conversion.

Aspartic Acid↗
PubMed2003

Structure and dynamics of metal ions in solution: QM/MM molecular dynamics simulations of Mn(2+) and V(2+).

Structural and dynamical properties of the transition metal ions V(2+) and Mn(2+) in aqueous solution, resulting from combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations have been compared. The necessity of polarization functions on the ligand's oxygen for a satisfactory description of such ions in aqueous solution is shown using V(2+) as test case. Radial distribution functions, coordination number distributions, and several angle distributions were pursued for a detailed structural comparison of the first hydration shells. Dynamical properties, such as the librational and vibrational motions of water molecules were evaluated by means of velocity autocorrelation functions. Approximative normal coordinate analyses were employed to calculate the rotational frequencies and vibrational motions around the three principal axes. The very low exchange rates for the first shell water exchanges only allow an investigation of the water exchange processes in the second shell, which take place within the picosecond range.

Cations, Divalent↗
PubMed1993

Molecular dynamics simulation of 1,2-dilauroyl-L-phosphatidylethanolamine binding to phospholipase A2: an attempt to explain the selective hydrolysis of substrate fatty acid ester at position 2.

To improve our understanding of why phospholipase A2 (PLA2) specifically catalyzes the hydrolysis of the fatty acid ester bond at position 2, not at position 1, of 1,2-diacyl-3-sn-phosphoglycerides, the binding of each fatty acid chain of 1,2-dilauroyl-L-phosphatidyl-ethanolamine (DLPE), a natural substrate, to bovine pancreas PLA2 was examined by molecular dynamics (MD) simulations. Two different binding modes were considered, i.e., the respective hydrocarbon chains of 1- and 2-lauroyl fatty acid esters were located at the PLA2 binding sites usually observed in the complex crystals (Form A2) and at the reverse sites (Form A1). Although the total energies of both forms fluctuated within nearly the same range during the 80 ps MD simulations, the binding mode of DLPE to the PLA2 catalytic site through the coordination to Ca2+ was much more advantageous in Form A2 than that in Form A1; significant deviation of the Ca2+ position from its starting structure was observed in the MD simulation of Form A1. The result suggests the importance of Ca2+ in the selective recognition and catalytic function of PLA2 toward the 2-positioned fatty acid ester of phosphoglyceride substrates.

Animals↗
PubMed2000

Successful molecular dynamics simulation of the zinc-bound farnesyltransferase using the cationic dummy atom approach.

Farnesyltransferase (FT) inhibitors can suppress tumor cell proliferation without substantially interfering with normal cell growth, thus holding promise for cancer treatment. A structure-based approach to the design of improved FT inhibitors relies on knowledge of the conformational flexibility of the zinc-containing active site of FT. Although several X-ray structures of FT have been reported, detailed information regarding the active site conformational flexibility of the enzyme is still not available. Molecular dynamics (MD) simulations of FT can offer the requisite information, but have not been applied due to a lack of effective methods for simulating the four-ligand coordination of zinc in proteins. Here, we report in detail the problems that occurred in the conventional MD simulations of the zinc-bound FT and a solution to these problems by employing a simple method that uses cationic dummy atoms to impose orientational requirement for zinc ligands. A successful 1.0 ns (1.0 fs time step) MD simulation of zinc-bound FT suggests that nine conserved residues (Asn127alpha, Gln162alpha, Asn165alpha, Gln195alpha, His248beta, Lys294beta, Leu295beta, Lys353beta, and Ser357beta) in the active site of mammalian FT are relatively mobile. Some of these residues might be involved in the ligand-induced active site conformational rearrangement upon binding and deserve attention in screening and design of improved FT inhibitors for cancer chemotherapy.

Alkyl and Aryl Transferases↗
PubMed2003

Determination of the active site protonation state of beta-secretase from molecular dynamics simulation and docking experiment: implications for structure-based inhibitor design.

Memapsin 2 (BACE) is an aspartyl protease known as beta-secretase that acts on the production of the beta-amyloid peptide in the human brain, a key event in the pathogenesis of Alzheimer's disease. Although it is expected that the net charge of the catalytic Asp diad would be -1 as in other kinds of aspartyl proteases, the exact protonation states of Asp32 and Asp228 have not been known without ambiguity. Two independent molecular dynamics (MD) simulations of BACE in complex with the potent inhibitor OM99-2 are carried out to determine the preferred protonation state of the Asp diad in the context that is consistent with the previous X-ray crystal structure. The results show that a strong hydrogen bond between the inhibitor hydroxyl group and Asp228 can be maintained only when Asp32 is neutral and Asp228 is ionized. The preference of this protonation state is further supported from the energetic and structural features found in the docking experiment of a novel potent inhibitor with the BACE active site. Thus, both MD and docking studies suggest that the role of hydrogen bond acceptor for the hydroxyl and piperazine groups of the inhibitors should be played by Asp228 instead of Asp32. This may be a key piece of information for the structure-based design/discovery of new inhibitor drugs.

Amyloid Precursor Protein Secretases↗
PubMed2006

Molecular dynamics simulations show that bound Mg2+ contributes to amino acid and aminoacyl adenylate binding specificity in aspartyl-tRNA synthetase through long range electrostatic interactions.

Molecular recognition between the aminoacyl-tRNA synthetase enzymes and their cognate amino acid ligands is essential for the faithful translation of the genetic code. In aspartyl-tRNA synthetase (AspRS), the co-substrate ATP binds preferentially with three associated Mg2+ cations in an unusual, bent geometry. The Mg2+ cations play a structural role and are thought to also participate catalytically in the enzyme reaction. Co-binding of the ATP x Mg3(2+) complex was shown recently to increase the Asp/Asn binding free energy difference, indicating that amino acid discrimination is substrate-assisted. Here, we used molecular dynamics free energy simulations and continuum electrostatic calculations to resolve two related questions. First, we showed that if one of the Mg2+ cations is removed, the Asp/Asn binding specificity is strongly reduced. Second, we computed the relative stabilities of the three-cation complex and the 2-cation complexes. We found that the 3-cation complex is overwhelmingly favored at ordinary magnesium concentrations, so that the protein is protected against the 2-cation state. In the homologous LysRS, the 3-cation complex was also strongly favored, but the third cation did not affect Lys binding. In tRNA-bound AspRS, the single remaining Mg2+ cation strongly favored the Asp-adenylate substrate relative to Asn-adenylate. Thus, in addition to their structural and catalytic roles, the Mg2+ cations contribute to specificity in AspRS through long range electrostatic interactions with the Asp side chain in both the pre- and post-adenylation states.

Adenosine Monophosphate↗
PubMedDate not supplied

Molecular dynamics simulation of POPC at low hydration near the liquid crystal phase transition.

We report results of a preliminary molecular dynamics study of a 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) bilayer at low hydration (5% in weight). Our results suggest a gel phase structure where the oleic chain is bent recalling the crystalline oleic acid structure. We have also found a discontinuity in the volume/temperature curve which might be related to the gel to liquid crystal phase transition in POPC.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1998

Alamethicin channels in a membrane: molecular dynamics simulations.

Alamethicin (Alm) is a 20 residue peptide which forms a kinked alpha-helix in membrane and membrane-mimetic environments. Ion channels formed by intramembraneous aggregates of Alm are thought to be formed by bundles of approximately parallel Alm helices surrounding a central bilayer pore. Different channel conductance levels correspond to different numbers of helices per bundle, ranging from N = 5 to N > 8. Calculation of the predicted pKA values of the ring of Glu18 sidechains at the C-terminal mouth of the pore suggests that at neutral pH most or all of these sidechains will remain protonated. Nanosecond molecular dynamics (MD) simulations of N = 5, 6, 7 and 8 bundles of Alm helices in a POPC bilayer have been run, corresponding to a total simulation time of 4 ns. These simulations explore the stability and conformational dynamics of these helix bundle channels when embedded in a full phospholipid bilayer in an aqueous environment. The structural and dynamic properties of water in these model channels are examined. As in earlier in vacuo simulations (J. Breed, R. Sankararamakrishnan, I. D. Kerr and M. S. P. Sansom, Biophys. J., 1996, 70, 1643) the dipole moments of water molecules within the pores are aligned antiparallel to the helix dipoles. This helps to contribute to the stability of the helix bundles.

Alamethicin↗
PubMed2005

Ultrafast dynamics of myoglobin without the distal histidine: stimulated vibrational echo experiments and molecular dynamics simulations.

Ultrafast protein dynamics of the CO adduct of a myoglobin mutant with the polar distal histidine replaced by a nonpolar valine (H64V) have been investigated by spectrally resolved infrared stimulated vibrational echo experiments and molecular dynamics (MD) simulations. In aqueous solution at room temperature, the vibrational dephasing rate of CO in the mutant is reduced by approximately 50% relative to the native protein. This finding confirms that the dephasing of the CO vibration in the native protein is sensitive to the interaction between the ligand and the distal histidine. The stimulated vibrational echo observable is calculated from MD simulations of H64V within a model in which vibrational dephasing is driven by electrostatic forces. In agreement with experiment, calculated vibrational echoes show slower dephasing for the mutant than for the native protein. However, vibrational echoes calculated for H64V do not show the quantitative agreement with measurements demonstrated previously for the native protein.

Histidine↗
PubMed1998

The geometry of intercalation complex of antitumor mitoxantrone and ametantrone with DNA: molecular dynamics simulations.

Intercalative binding of the antitumor drugs ametantrone and mitoxantrone to the dodecamer duplex d(CGCGAGCTCGCG)2 was studied by applying molecular dynamics in water with the GROMOS 87 force field. A number of reasonable binding orientations were tested by short pre-simulations. It was shown that in energetically favourable orientation the anthraquinone chromophore is perpendicular to the direction of inter-base hydrogen bonds. Helically shaped side-chains of the drugs fit to the minor groove. The best orientation obtained in pre-simulations was applied in the main simulations. Small but significant differences were found between structures of intercalation complexes of the two drugs with the dodecamer duplex, the mitoxantrone complex possessing more favourable energy. The molecular nature of interactions responsible for those differences has been discussed.

Antineoplastic Agents↗
PubMed2004

Molecular dynamics simulations of HPr under hydrostatic pressure.

The histidine-containing protein (HPr) plays an important role in the phosphotransferase system (PTS). The deformations induced on the protein structure at high hydrostatic pressure values (4, 50, 100, 150, and 200 MPa) were previously (H. Kalbitzer, A. Görler, H. Li, P. Dubovskii, A. Hengstenberg, C. Kowolik, H. Yamada, and K. Akasaka, Protein Science 2000, Vol. 9, pp. 693-703) analyzed by NMR experiments: the nonlinear variations of the amide chemical shifts at high pressure values were supposed to arise from induced shifts in the protein conformational equilibrium. Molecular dynamics (MD) simulations are here performed, to analyze the protein internal mobility at 0.1 MPa, and to relate the nonlinear variations of chemical shifts observed at high pressure, to variations in conformational equilibrium. The global features of the protein structure are only slightly modified along the pressure. Nevertheless, the values of the Voronoi residues volumes show that the residues of alpha-helices are more compressed that those belonging to the beta-sheet. The alpha-helices are also displaying the largest internal mobility and deformation in the simulations. The nonlinearity of the 1H chemical shifts, computed from the MD simulation snapshots, is in qualitative agreement with the nonlinearity of the experimentally observed chemical shifts.

Algorithms↗
PubMed2006

Distance between a native cofactor and a spin label in the reaction centre of Rhodobacter sphaeroides by a two-frequency pulsed electron paramagnetic resonance method and molecular dynamics simulations.

The distance between the paramagnetic state of a native cofactor and a spin label is measured in the photosynthetic reaction centre from the bacterium Rhodobacter sphaeroides R26. A two-frequency pulsed electron paramagnetic resonance method [double-electron-electron spin resonance (DEER)] is used. A distance of 3.05 nm between the semiquinone anion state of the primary acceptor (Q(A)) and the spin label at the native cysteine at position 156 in the H-subunit is found. Molecular-dynamics (MD) simulations are performed to interpret the distance. A 6 ns run comprising the entire RC protein yields a distance distribution that is close to the experimental one. The average distance found by the MD simulation is smaller than the distance obtained by DEER by at least 0.2 nm. To better represent the experiments performed at low temperature (60K), a MD method to mimic the freezing-in of the room-temperature conformations is introduced. Both MD methods yield similar distances, but the second method has a trend towards a wider distance distribution.

Anions↗
PubMed2006

Molecular dynamics simulations of the aqueous interface with the [BMI][PF6] ionic liquid: Comparison of different solvent models.

We report a Molecular Dynamics (MD) study of the interface between water and the hygroscopic room temperature Ionic Liquid "IL" [BMI][PF6] (1-butyl-3-methyl-imidazolium hexafluorophosphate), comparing the TIP3P, SPC/E and TIP5P models for water and two IL models where the ions are +/-1 or +/-0.9 charged. A recent MD study (A. Chaumont, R. Schurhammer and G. Wipff, J. Phys. Chem. B, 2005, 109, 18964) showed that using TIP3P water in conjunction with the IL(+/-1) model led to water-IL mixing without forming an interface, whereas a biphasic system could be obtained with the IL(+/-0.9) model. With the TIP5P and SPC/E models, the juxtaposed aqueous and IL phases are found to remain distinct for at least 20 ns. The resulting IL humidity, exaggerated with the IL(+/-1) model, is in better agreement with experiment using the IL(+/-0.9) model. We also report demixing simulations on the "randomly mixed" liquids, using the IL(+/-0.9) model for the ionic liquid. With the three tested water models, the phases separate very slowly ( approximately 20 ns or more) compared to "classical" chloroform-water mixtures (less than 1 ns), leading to biphasic systems similar to those obtained after equilibration of the juxtaposed liquids. The characteristics of the interface (size, polarity, ion orientation, electrostatic potential) are compared with the different models. Possible reasons why, among the three tested water models, the widely-used TIP3P model exaggerates the inter-solvent mixing, are analyzed. The difficulty in computationally and experimentally equilibrating water-IL mixtures is attributed to the slow dynamics and micro-heterogeneity of the IL and to the different states of water in the IL phase.

Computer Simulation↗
PubMed1995

Characterization of the internal motions of Escherichia coli ribonuclease HI by a combination of 15N-NMR relaxation analysis and molecular dynamics simulation: examination of dynamic models.

The backbone dynamics of Escherichia coli ribonuclease HI (RNase HI) in the picosecond to nanosecond time scale were characterized by a combination of measurements of 15N-NMR relaxation (T1, T2, and NOE), analyzed by a model-free approach, and molecular dynamics (MD) simulation in water. The MD simulations in water were carried out with long-range Coulomb interactions to avoid the artificial fluctuation caused by the cutoff approximation. The model-free analysis of the 15N-NMR relaxation indicated that RNase HI has a rotational correlation time of 10.9 ns at 27 degrees C. The generalized order parameter (S2) for the internal motions varied from 0.15 to 1.0, with an average value of 0.85, which is much larger than that of the RNase H domain of HIV-1 reverse transcriptase (0.78). Large internal motions (small order parameters) were observed in the N-terminal region (Leu2-Lys3), the loop between beta-strands A and B (Cys13-Gly15), the turn between alpha-helix I and beta-strand D (Glu61, His62), the loop between beta-strand D and alpha-helix II (Asp70-Tyr71), the loop between alpha-helices III and IV (Ala93-Lys96), the loop between beta-strand E and alpha-helix V (Gly123-His127), and the C-terminal region (Gln152-Val155). The effective correlation time observed in these regions varied from 0.45 ns (Glu61, Lys96) to 2.2 ns (Leu14). The order parameters calculated from the MD agreed well with those from the NMR experiment, with a few exceptions. The distributions of most of the backbone N-H vectors obtained by MD are approximately consistent with the diffusion-in-a-cone model. These distributions, however, were elliptic, with a long axis perpendicular to the plane defined by the N-H and N-C alpha vectors. Distributions supporting the axial fluctuation model or the jump-between-two-cones model were also observed in the MD simulation.

Binding Sites↗
PubMed1993

Molecular dynamics simulations of poly(dA).poly(dT): comparisons between implicit and explicit solvent representations.

The program AMBER 3.0 has been used to generate molecular dynamics trajectories of a poly(dA).poly(dT) decamer. The simulations were performed using different methods to treat solvent effects. Results of a simulation including 18 counterions NH4+ and 4109 water molecules under (N, P, T) conditions were compared to simulation runs with implicit solvent representation in which solvent screening effects were represented by the use of a sigmoidal distance-dependent dielectric function. In the latter case, the system was simulated under microcanonical (N, V, E) and canonical (N, V, T) conditions. For the fully hydrated system simulation, a preequilibration protocol was developed since it was observed that long and progressive periods of heating and equilibration on the overall system were necessary in order to avoid energetic collisions between the solute and the solvent molecules, leading to severe irreversible deformation of the solute. A detailed analysis of DNA conformations, sugar puckers, and stability of the hydrogen bonds, Watson-Crick and three-center H bonds, is reported. The results show that DNA remains essentially in the B conformer with a tendency in the hydrated model to adopt a slightly distorted, unwound, and stretched conformation in comparison to standard B-DNA. Concerning sugar puckers, the mean pseudorotation phases of the adenine residues are systematically higher than those of the thymine residues, except in the case of the hydrated model for which a articular behavior is observed for the adenine strand. In this case, the terminal bases oscillate between C2'-endo and O4'-endo and the central ones stay in the C3'-endo domain. The mean lifetimes of the internal Watson-Crick H-bond (A) HN6...O4(T) are also dependent on the base pairs included in the calculation, excepted for the implicit solvent simulation at constant temperature. The three-center H bonds have very small mean lifetimes in all three cases of MD simulation. In the minor groove of the hydrated model, a spine of hydration is found as observed by x-ray crystallography and other theoretical simulations. On the basis of the rms deviations, it appears that the fully hydrated simulation has not reached a plateau at the end of the run, while the implicit simulation at constant energy seems to have converged. At constant temperature, very large oscillations in rms deviations are observed.

Base Sequence↗
PubMed2005

Barnase thermal titration via molecular dynamics simulations: detection of early denaturation sites.

The thermal stability of barnase has been studied using constant pressure and temperature (CPT) molecular dynamics at different temperatures. Barnase X-ray coordinates were obtained from the Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank (PDB code:1rnb). Simulations were performed at 285, 295, 300, 335, 345, and 395 K in explicit water under periodic boundary conditions for 280 ps. For each simulation, conformations were saved every 0.2 ps. Root mean square deviation (RMSD) values were calculated relative to the starting structure at 300 K and at time t = 0. Root mean square fluctuation (RMSF) values were calculated relative to the average structure obtained from the 300K simulation. Both root mean square deviation and fluctuation analysis indicated the presence of discrete regions of hyper-sensitivity along the barnase polypeptide chain. These regions exhibited spikes in flexibility prior to any global structural changes. The specific changes in barnase backbone flexibility are accompanied by increased phi/psi angle fluctuations. These results suggest the presence of early denaturation sites or denaturation nuclei whose local structure is disrupted prior to global structure disruption. Identification of denaturation nuclei suggests that appropriate amino acid replacements at these sites may lead to the design and development of more stable barnase mutants. This strategy of identifying denaturation nuclei in protein structures may represent a first step in the design of more stable protein structures.

Amino Acid Sequence↗
PubMed2003

Investigation of the induced-fit mechanism and catalytic activity of the human cytomegalovirus protease homodimer via molecular dynamics simulations.

Human cytomegalovirus (HCMV) is a highly species-specific DNA virus infecting up to 80% of the general population. The viral genome contains the open reading frame UL80, which encodes the full-length 80 kDa HCMV serine protease and its substrate. Full-length HCMV protease is composed of an N-terminal 256-amino-acid proteolytic domain, called assemblin, a linker region, and a C-terminal structural domain, the assembly protein precursor. Biochemical studies have shown that dimerization activates assemblin because of an induced stabilization of the oxyanion hole (Arg166). Thus, we performed here molecular dynamics (MD) simulations on HCMV protease models to study the induced-fit mechanism of the enzyme upon the binding of substrates and peptidyl inhibitors, and structural and energetic factors that are responsible for the catalytic activity of the enzyme dimer. Long and stable trajectories were obtained for the models of the monomeric and dimeric states, free in solution and bound to a peptidyl-activated carbonyl inhibitor, with very good agreement between theoretical and experimental results. Our results suggest that HCMV protease is indeed a novel example of serine protease that operates by an induced-fit mechanism. Also, in agreement with mutagenesis studies, our MD simulations suggest that the dimeric form is necessary to activate the enzyme because of an induced stabilization of the oxyanion hole.

Arginine↗
PubMed2005

Ion solvation in water from molecular dynamics simulation with the ABEEM/MM force field.

A systematic study on monovalent ions in water clusters and in aqueous solution is presented for providing insight into their solvation structures, charge distributions, binding energies, as well as dynamic and thermodynamic properties in terms of the atom-bond electronegativity equalization method fused into molecular mechanics (ABEEM/MM) that is to take ABEEM charges into the Coulomb term in MM. For hydrated systems of M+(H2O)n, M+ being Li+, Na+, and K+, as well as X-(H2O)n, X- being F-, Cl-, and Br-, with n = 1-6, parameters for the effective interaction between the ion and the water molecules were determined, so as to reproduce the experimental or ab initio results. The corresponding parameters were tested with molecular dynamics (MD) simulations of these ions in liquid water and with solvation free energy calculations using the perturbation technique. The results of aqueous ionic solution simulations with the ABEEM/MM force field provide a reasonable description of many important properties, which are in good agreement with the experimental measurements. This work demonstrates that the combination of ABEEM/MM-MD provides a powerful tool in analyzing solvation processes of monovalent ions in water.

Chemistry, Physical↗