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,711 records · Page 95Linked to original sources
PubMed2005

Sequence dependence of amyloid fibril formation: insights from molecular dynamics simulations.

The clarification of the physico-chemical determinants underlying amyloid deposition is critical for our understanding of misfolding diseases. With this purpose we have performed a systematic all-atom molecular dynamics (MD) study of a series of single point mutants of the de novo designed amyloidogenic peptide STVIIE. Sixteen different 50ns long simulations using explicit solvent have been carried out starting from four different conformations of a polymeric six-stranded beta-sheet. The simulations have provided evidence for the influence of a small number of site-specific hydrophobic interactions on the packing and stabilization of nascent aggregates, as well as the interplay between side-chain interactions and the net charge of the molecule on the strand arrangement of polymeric beta-sheets. This MD analysis has also shed light into the origin of the position dependence on mutation of beta-sheet polymerization that was found experimentally for this model system. Our results suggest that MD can be applied to detect critical positions for beta-sheet aggregation within a given amyloidogenic stretch. Studies similar to the one presented here can guide site-directed mutations or the design of drugs that specifically disrupt the key stabilizing interactions of beta-sheet aggregates.

Amino Acid Sequence↗
PubMed2005

Molecular dynamics simulations of alkylsilane stationary-phase order and disorder. 2. Effects of temperature and chain length.

In an effort to elucidate the molecular-level structural features that control shape-selective separations, we have investigated the molecular dynamics of chromatographic models that represent both monomeric and polymeric stationary phases with alkylsilane length and temperature conditions analogous to actual materials of low to high shape selectivity. The structural characterization of these models is consistent with previous experimental observations of alkyl chain order and disorder: alkyl chain order increases both with alkyl chain length and with reduced temperature. Models that represent shape-selective reversed-phase liquid chromatography (RPLC) phases possess a significant region of distal end chain order with primarily trans dihedral angle conformations; the extension of these ordered regions into the phase increases with an increase in chain length. Models with extended chain length (C30) possess a higher degree of conformational order and are relatively insensitive to changes in surface coverage, bonding chemistry, and temperature. Chromatography models of various chain lengths and over a temperature range that represents highly shape-selective RPLC stationary phases all contain a series of well-defined and rigid cavities; the size and depth of these "slots" increase for the C30 models, which may promote the enhanced separations of larger size shape-constrained solutes, such as carotenoids.

Journal Article↗
PubMed2004

Model-free analysis of protein dynamics: assessment of accuracy and model selection protocols based on molecular dynamics simulation.

The popular model-free approach to analyze NMR relaxation measurements has been examined using artificial amide (15)N relaxation data sets generated from a 10 nanosecond molecular dynamics trajectory of a dihydrofolate reductase ternary complex in explicit water. With access to a detailed picture of the underlying internal motions, the efficacy of model-free analysis and impact of model selection protocols on the interpretation of NMR data can be studied. In the limit of uncorrelated global tumbling and internal motions, fitting the relaxation data to the model-free models can recover a significant amount of quantitative information on the internal dynamics. Despite a slight overestimation, the generalized order parameter is quite accurately determined. However, the model-free analysis appears to be insensitive to the presence of nanosecond time scale motions with relatively small magnitude. For such cases, the effective correlation time can be significantly underestimated. As a result, proteins appear to be more rigid than they really are. The model selection protocols have a major impact on the information one can reliably obtain. The commonly employed protocol based on step-up hypothesis testing has severe drawbacks of oversimplification and underfitting. The consequences are that the order parameter is more severely overestimated and the correlation time more severely underestimated. Instead, model selection based on Bayesian Information Criteria (BIC), recently introduced to the model-free analysis by d'Auvergne and Gooley (2003), provides a better balance between bias and variance. More appropriate models can be selected, leading to improved estimate of both the order parameter and correlation time. In addition, the computational cost is significantly reduced and subjective parameters such as the significance level are unnecessary.

Anisotropy↗
PubMed1991

Experimentally verifying molecular dynamics simulations through fluorescence anisotropy measurements.

The fluorescence anisotropy decay of the single tryptophan residue in phospholipase A2 was studied by use of differential polarized phase fluorometry and computer simulations of protein dynamics. The results enable the verification of a simulated dynamic event by direct experimental measurement on the same time scale. When all hydrogen atoms are modeled explicitly, the simulations agree well with the experimental measurements. However, the measurements contradict simulations in which nonpolar hydrogens are incorporated into "extended" or "united" atoms. These simulations predict an anisotropy decay in excess of measured values and appear to seriously underestimate the electrostatic interactions occurring between water and aromatic side chains. The results support the general validity of studying protein dynamics with the molecular-mechanics approach and illustrate a potentially serious deficiency of simulations which do not explicitly model all hydrogen atoms.

Animals↗
PubMed1995

Internal mobility of the basic pancreatic trypsin inhibitor in solution: a comparison of NMR spin relaxation measurements and molecular dynamics simulations.

Order parameters as well as longitudinal and transverse relaxation rates are calculated for the backbone 15N and 13C alpha nuclei of the basic pancreatic trypsin inhibitor (BPTI) from a 1000 ps molecular dynamics trajectory in explicit water at 277 K using the "model free" approach of Lipari and Szabo. New NMR relaxation data at 277 K are presented, and a comparison is made between NMR relaxation measurements and molecular dynamics relaxation data. It is found that the relaxation processes determining the longitudinal (T1) relaxation rates are inadequately sampled even during this length of simulation. In effect, the calculated relaxation rates are determined almost solely by the order parameters and the overall rotational correlation time of the protein, which appears to be in clear contrast to experimental relaxation rates.

Amino Acid Sequence↗
PubMed1992

Trefoil knotting revealed by molecular dynamics simulations of supercoiled DNA.

Computer simulations of the supercoiling of DNA, largely limited to stochastic search techniques, can offer important information to complement analytical models and experimental data. Through association of an energy function, minimum-energy supercoiled conformations, fluctuations about these states, and interconversions among forms may be sought. In theory, the observation of such large-scale conformational changes is possible, but modeling and numerical considerations limit the picture obtained in practice. A new computational approach is reported that combines an idealized elastic energy model, a compact B-spline representation of circular duplex DNA, and deterministic minimization and molecular dynamics algorithms. A trefoil knotting result, made possible by a large time-step dynamics scheme, is described. The simulated strand passage supports and details a supercoiled-directed knotting mechanism. This process may be associated with collective bending and twisting motions involved in supercoiling propagation and interwound branching. The results also demonstrate the potential effectiveness of the Langevin/implicit-Euler dynamics scheme for studying biomolecular folding and reactions over biologically interesting time scales.

Chemical Phenomena↗
PubMed2006

Molecular dynamics simulation of hydrogen fluoride mixtures with 1-ethyl-3-methylimidazolium fluoride: a simple model for the study of structural features.

Mixtures of hydrogen fluoride with ionic liquids show unique physicochemical properties, including their ability to form polyfluoride species (pointed out for the first time in this media by von Rosenvinge et al. J. Chem. Phys. 1997, 107, 8012). Among those systems the acidic 1-ethyl-3-methylimidazolium fluoride (EMIF.2.3HF) has been widely studied experimentally since it is the more promising for electrochemical applications. Recent studies (Hagiwara et al. J. Electrochem. Soc. 2002, 149, D1), while yielding many results, raised some questions about structural features of the liquid: absence of hydrogen bonds between the EMI+ ring hydrogen atoms and the fluoride anions, persistence of stacks and layers of cations similar to those existing in the crystal, and interpretation of the X-ray diffraction spectra. To address these questions, we have developed a simple molecular dynamics model. Our simulations are very consistent with experimental results and complete them, providing an atomic scale interpretation.

Journal Article↗
PubMed1999

X-ray structure and conformational dynamics of the HIV-1 protease in complex with the inhibitor SDZ283-910: agreement of time-resolved spectroscopy and molecular dynamics simulations.

Based on the X-ray structure of the human immunodeficiency virus type-1 (HIV-1) protease in complex with the statine-derived inhibitor SDZ283-910, a 542 ps molecular dynamics trajectory was computed. For comparison with the 805 ps trajectory obtained for the uncomplexed enzyme, the theoretical fluorescence anisotropy decay of the unliganded protease and the inhibitor complex was calculated from the trajectories of the Trp6A/Trp6B and Trp42A/Trp42B transition dipole moments. This enabled us to directly compare the simulated data with the experimental picosecond time-resolved fluorescence data. Fitting both experimental and simulated data to the Kohlrausch-Williams-Watts (KWW) function exp(-t/tauk)beta revealed a very good agreement for the uncomplexed protease as well as for the SDZ283-910 complex. Binding of the inhibitor induced a faster decay of both the experimental and the computed protease fluorescence anisotropy decay. By this integrative approach, the atomic detail of inhibitor-induced changes in the conformational dynamics of the HIV-1 protease was experimentally verified and will be used for further inhibitor optimisation.

Anti-HIV Agents↗
PubMed2006

HD isotope effect on the dihydrogen bond of NH+4...BeH2 by ab initio path integral molecular dynamics simulation.

In order to investigate the HD isotope effect on a dihydrogen bonded cation system, we have studied NH+4...BeH2 and its isotopomers by ab initio path integral molecular dynamics. It is found that the dihydrogen bond can be exchanged by NH+(4) rotation. The deuterated isotopomer (ND+(4)...BeD(2); DD) can exchange the dihydrogen bond more easily than other isotopomers such as (NH+4...BeH2; HH). This unusual isotope effect is ascribed to the "quantum localization" which occurs when the effective energy barrier for the rotational mode becomes higher by the zero point energy of other modes. We also found that the binding energy of dihydrogen bonds for DD species is the smallest among the isotopomers.

Journal Article↗
PubMed2004

Quantum molecular dynamics simulation of proton transfer in cytochrome c oxidase.

Proton transfer/translocation is studied in cytochrome c oxidase (CcO) by a combination of quantum mechanics (QM) for the transferring protons and classical molecular dynamics (MD) for the protein and solvent. The possibility of a glutamate, Glu286 in the Rhodobacter sphaeroides numbering scheme, acting as a rely point for proton translocation is investigated. The MD finds a hydrogen-bonded cycle of two waters and the carboxylate oxygens of Glu286. The possibility of protonating Glu286 to form neutral GluH is studied and we find that, as experimentally inferred, this glutamate can spend most of its time as GluH. Since translocation relies on the presence of water chains within CcO channels, MD is used to assess their formation. Glu286 and Mg(2+) can be connected by continuous hydrogen-bonded chains that are robust, though transient, and the protein appears spongy above (toward the outer membrane) the Mg(2+). In contrast, the D-channel spanning Asp132, close to the inner membrane surface, to Glu286, forms water chains that are much sparser and do not continuously connect these residues. Rather, there are chains spanning Glu286 to the vicinity of Asn140, and other more robust and ramified water structures that connect Asp132 with waters close to Asn140.

Asparagine↗
PubMed1999

Molecular dynamics simulations of unfolding and refolding of a beta-hairpin fragment of protein G.

We have studied the unfolding and refolding pathway of a beta-hairpin fragment of protein G by using molecular dynamics. Although this fragment is small, it possesses several of the qualities ascribed to small proteins: cooperatively formed beta-sheet secondary structure and a hydrophobic "core" of packed side chains. At high temperatures, we find that the beta-hairpin unfolds through a series of sudden, discrete conformational changes. These changes occur between states that are identified with the folded state, a pair of partially unfolded kinetic intermediates, and the unfolded state. To study refolding at low temperatures, we perform a series of short simulations starting from the transition states of the discrete transitions determined by the unfolding simulations.

Amino Acid Sequence↗
PubMed2004

Severe acute respiratory syndrome coronavirus 3C-like proteinase N terminus is indispensable for proteolytic activity but not for enzyme dimerization. Biochemical and thermodynamic investigation in conjunction with molecular dynamics simulations.

Severe acute respiratory syndrome (SARS) coronavirus is a novel human coronavirus and is responsible for SARS infection. SARS coronavirus 3C-like proteinase (SARS 3CL(pro)) plays key roles in viral replication and transcription and is an attractive target for anti-SARS drug discovery. In this report, we quantitatively characterized the dimerization features of the full-length and N-terminal residues 1-7 deleted SARS 3CL(pro)s by using glutaraldehyde cross-linking SDS-PAGE, size-exclusion chromatography, and isothermal titration calorimeter techniques. Glutaraldehyde cross-linking SDS-PAGE and size-exclusion chromatography results show that, similar to the full-length SARS 3CL(pro), the N-terminal deleted SARS 3CL(pro) still remains a dimer/monomer mixture within a wide range of protein concentrations. Isothermal titration calorimeter determinations indicate that the equilibrium dissociation constant (K(d)) of the N-terminal deleted proteinase dimer (262 microm) is very similar to that of the full-length proteinase dimer (227 microm). Enzymatic activity assay using the fluorescence resonance energy transfer method reveals that N-terminal deletion results in almost complete loss of enzymatic activity for SARS 3CL(pro). Molecular dynamics and docking simulations demonstrate the N-terminal deleted proteinase dimer adopts a state different from that of the full-length proteinase dimer, which increases the angle between the two protomers and reduces the binding pocket that is not beneficial to the substrate binding. This conclusion is verified by the surface plasmon resonance biosensor determination, indicating that the model substrate cannot bind to the N-terminal deleted proteinase. These results suggest the N terminus is not indispensable for the proteinase dimerization but may fix the dimer at the active state and is therefore vital to enzymatic activity.

Binding Sites↗
PubMed2004

Generalized-ensemble algorithms: enhanced sampling techniques for Monte Carlo and molecular dynamics simulations.

In complex systems with many degrees of freedom such as spin glass and biomolecular systems, conventional simulations in canonical ensemble suffer from the quasi-ergodicity problem. A simulation in generalized ensemble performs a random walk in potential energy space and overcomes this difficulty. From only one simulation run, one can obtain canonical ensemble averages of physical quantities as functions of temperature by the single-histogram and/or multiple-histogram reweighting techniques. In this article we review the generalized ensemble algorithms. Three well-known methods, namely, multicanonical algorithm (MUCA), simulated tempering (ST), and replica-exchange method (REM), are described first. Both Monte Carlo (MC) and molecular dynamics (MD) versions of the algorithms are given. We then present five new generalized-ensemble algorithms which are extensions of the above methods.

Algorithms↗
PubMed2003

The role of side-chain interactions in the early steps of aggregation: Molecular dynamics simulations of an amyloid-forming peptide from the yeast prion Sup35.

Understanding the early steps of aggregation at atomic detail might be crucial for the rational design of therapeutics preventing diseases associated with amyloid deposits. In this paper, aggregation of the heptapeptide GNNQQNY, from the N-terminal prion-determining domain of the yeast protein Sup35, was studied by 20 molecular dynamics runs for a total simulation time of 20 micros. The simulations generate in-register parallel packing of GNNQQNY beta-strands that is consistent with x-ray diffraction and Fourier transform infrared data. The statistically preferred aggregation pathway does not correspond to a purely downhill profile of the energy surface because of the presence of enthalpic barriers that originate from out-of-register interactions. The parallel beta-sheet arrangement is favored over the antiparallel because of side-chain contacts; in particular, stacking interactions of the tyrosine rings and hydrogen bonds between amide groups. No ordered aggregation was found in control simulations with the mutant sequence SQNGNQQRG in accord with experimental data and the strong sequence dependence of aggregation.

Alanine↗
PubMed2006

The transition between the B and Z conformations of DNA investigated by targeted molecular dynamics simulations with explicit solvation.

The transition between the B and Z conformations of double-helical deoxyribonucleic acid (DNA) belongs to the most complex and elusive conformational changes occurring in biomolecules. Since the accidental discovery of the left-handed Z-DNA form in the late 1970s, research on this DNA morphology has been engaged in resolving questions relative to its stability, occurrence, and function in biological processes. While the occurrence of Z-DNA in vivo is now widely recognized and the major factors influencing its thermodynamical stability are largely understood, the intricate conformational changes that take place during the B-to-Z transition are still unknown at the atomic level. In this article, we report simulations of this transition for the 3'-(CGCGCG)-5' hexamer duplex using targeted molecular dynamics with the GROMOS96 force field in explicit water under different ionic-strength conditions. The results suggest that for this oligomer length and sequence, the transition mechanism involves: 1), a stretched intermediate conformation, which provides a simple solution to the important sterical constraints involved in this transition; 2), the transient disruption of Watson-Crick hydrogen-bond pairing, partly compensated energetically by an increase in the number of solute-solvent hydrogen bonds; and 3), an asynchronous flipping of the bases compatible with a zipperlike progression mechanism.

Computer Simulation↗
PubMed2004

Molecular dynamics simulation of the adsorption of a fibronectin module on a graphite surface.

We report atomistic simulations of the adsorption of a fibronectin type I module on a hydrophobic graphite surface. This module comprises only beta-sheets, unlike the albumin fragments previously investigated by us which contained only alpha-helices (Raffaini, G.; Ganazzoli, F. Langmuir 2003, 19, 3403-3412). As done in the latter case, most simulations are carried out in an effective dielectric medium by energy minimizations and molecular dynamics (MD). Further optimizations and MD runs in the explicit presence of water are also performed to assess the stability of the geometries found and to describe the solvation of the adsorbed fibronectin module. The initial adsorption is accompanied by local rearrangements of the strands in contact with the surface, but the overall molecular structure is largely preserved. Much larger rearrangements take place at longer times as found through the MD runs, with the molecule spreading as much as possible so as to maximize the surface coverage, hence the interaction energy, despite a significant strain energy. Energetic aspects of adsorption together with the concomitant size change are discussed in comparison with our previous results for two albumin fragments.

Adsorption↗
PubMed2005

pH-dependent conformational flexibility of the SARS-CoV main proteinase (M(pro)) dimer: molecular dynamics simulations and multiple X-ray structure analyses.

The SARS coronavirus main proteinase (M(pro)) is a key enzyme in the processing of the viral polyproteins and thus an attractive target for the discovery of drugs directed against SARS. The enzyme has been shown by X-ray crystallography to undergo significant pH-dependent conformational changes. Here, we assess the conformational flexibility of the M(pro) by analysis of multiple crystal structures (including two new crystal forms) and by molecular dynamics (MD) calculations. The MD simulations take into account the different protonation states of two histidine residues in the substrate-binding site and explain the pH-activity profile of the enzyme. The low enzymatic activity of the M(pro) monomer and the need for dimerization are also discussed.

Binding Sites↗