Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD 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,423 records · Page 79Linked to original sources
PubMed2006

UMMS: constrained harmonic and anharmonic analyses of macromolecules based on elastic network models.

UMass Morph Server (UMMS) has been developed for the broad impact on the study of molecular dynamics (MD). The elastic network model (ENM) of a given macromolecule has been proven as a useful tool for analyzing thermal behaviors locally and predicting folding pathways globally. UMMS utilizes coarse-grained ENMs at various levels. These simplifications remarkably save computation time compared with all-atom MD simulations so that one can bring down massive computational problems from a supercomputer to a PC. To improve computational efficiency and physical reality of ENMs, the symmetry-constrained, rigid-cluster, hybrid and chemical-bond ENMs have been developed and implemented at UMMS. One can request both harmonic normal mode analysis of a single macromolecule and anharmonic pathway generation between two conformations of a same molecule using elastic network interpolation at http://biomechanics.ecs.umass.edu/umms.html.

Computational Biology↗
PubMed2006

The molecular mechanism of DNA damage recognition by MutS homologs and its consequences for cell death response.

We determined the molecular mechanism of cell death response by MutS homologs in distinction to the repair event. Key protein-DNA contacts differ in the interaction of MutS homologs with cisplatinated versus mismatched DNA. Mutational analyses of protein-DNA contacts, which were predicted by molecular dynamics (MD) simulations, were performed. Mutations in suggested interaction sites can affect repair and cell death response independently, and to different extents. A glutamate residue is identified as the key contact with cisplatin-DNA. Mutation of the residue increases cisplatin resistance due to increased non-specific DNA binding. In contrast, the conserved phenylalanine that is instrumental and indispensable for mismatch recognition during repair is not required for cisplatin cytotoxicity. These differences in protein-DNA interactions are translated into localized conformational changes that affect nucleotide requirements and inter-subunit interactions. Specifically, the ability for ATP binding/hydrolysis has little consequence for the MMR-dependent damage response. As a consequence, intersubunit contacts are altered that most likely affect the interaction with downstream proteins. We here describe the interaction of MutS homologs with DNA damage, as it differs from the interaction with a mismatch, and its structural translation into all other functional regions of the protein as a mechanism to initiate cell death response and concomitantly inhibit repair.

Adenosine Triphosphatases↗
PubMed2006

Mutagenic nucleotide incorporation and hindered translocation by a food carcinogen C8-dG adduct in Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4): modeling and dynamics studies.

Bulky carcinogen-DNA adducts commonly cause replicative polymerases to stall, leading to a switch to bypass polymerases. We have investigated nucleotide incorporation opposite the major adduct of 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) in the DinB family polymerase, Dpo4, using molecular modeling and molecular dynamics (MD) simulations. PhIP, the most prevalent heterocyclic aromatic amine formed by cooking of proteinaceous food, is mutagenic in mammalian cells and is implicated in mammary and colon tumors. Our results show that the dG-C8-PhIP adduct can be accommodated in the spacious major groove Dpo4 open pocket, with Dpo4 capable of incorporating dCTP, dTTP or dATP opposite the adduct reasonably well. However, the PhIP ring system on the minor groove side would seriously disturb the active site, regardless of the presence and identity of dNTP. Furthermore, the simulations indicate that dATP and dTTP are better incorporated in the damaged system than in their respective mismatched but unmodified controls, suggesting that the PhIP adduct enhances incorporation of these mismatches. Finally, bulky C8-dG adducts, situated in the major groove, are likely to impede translocation in this polymerase (Rechkoblit et al. (2006), PLoS Biol., 4, e11). However, N2-dG adducts, which can reside on the minor groove side, appear to cause less hindrance when in this position.

Binding Sites↗
PubMed2002

The molecular dynamics thermocycler. A new approach to sample conformational space, as exemplified by the RNA hairpin.

The approach called the molecular dynamics thermocycler employs repeated cycles of heating and cooling of the model during the course of unrestrained in aqua MD simulation. It is validated by confronting results of three 1.6 ns simulations using different heating temperatures, starting from a de novo constructed 12-mer RNA hairpin model containing a UGAA loop vs. the set of NMR derived coordinates (PDB code 1AFX). The lowest average RMS deviation and best coherence with the experimentally delivered structures was identified in case of the molecular dynamics thermocycler using Tmax = 400 K.

Nuclear Magnetic Resonance, Biomolecular↗
PubMed1991

Relative binding free energy calculations of inhibitors to two mutants (Glu46----Ala/Gln) of ribonuclease T1 using molecular dynamics/free energy perturbation approaches.

We present free energy perturbation calculations on the complexes of Glu46----Ala46 (E46A) and Glu46----Gln46 (E46Q) mutants of ribonuclease T1 (RNaseT1) with inhibitors 2'-guanosine monophosphate (GMP) and 2'-adenosine monophosphate (AMP) by a thermodynamic perturbation method implemented with molecular dynamics (MD). Using the available crystal structure of the RNaseT1-GMP complex, the structures of E46A-GMP and E46Q-GMP were model built and equilibrated with MD simulations. The structures of E46A-AMP and E46Q-AMP were obtained as a final structure of the GMP----AMP perturbation calculation respectively. The calculated difference in the free energy of binding (delta delta Gbind) was 0.31 kcal/mol for the E46A system and -1.04 kcal/mol for the E46Q system. The resultant free energies are much smaller than the experimental and calculated value of approximately 3 kcal/mol for the native RNaseT1, which suggests that both mutants have greater relative adenine affinities than native RNaseT1. Especially E46Q is calculated to have a larger affinity for adenine than guanine, as we suggested previously from the calculation on the native RNaseT1. Thus, the molecular dynamics/free energy perturbation method may be helpful in protein engineering, directed toward increasing or changing the substrate specificity of enzymes.

Adenosine Monophosphate↗
PubMed1996

Construction and evaluation of a three-dimensional structure of cytochrome P450choP enzyme (CYP105C1).

The purpose of this work was to develop and carefully evaluate improved strategies for constructing reliable 3-D models of P450 isozymes. To this end, a unique combination of steps for building and evaluating a model structure was used to build a homology model of the P450choP isozyme, based on knowledge of the X-ray structures of P450cam, P450terp, P450BM-3 and P450eryF. Specifically, the reliability of this model was examined by systematic comparisons of its conformational, energetic, environmental and packing properties and those of the four reference proteins with corresponding properties from the database of proteins with known structures. The results showed that the examined properties of this model structure are well within the criteria established for reliable structures and are of nearly as good quality as those of the reference proteins. In addition, the result from a 120 ps unconstrained MD simulation of the model with structural waters provided evidence that the model is stable at room temperature. This 3-D model can now be reliably used for explicit characterization of substrate and inhibitor complexes. Most importantly, although it is envisioned that building models for mammalian P450s will be even more challenging, the steps described here should be very useful in future construction of 3-D models of mammalian P450 isozymes.

Amino Acid Sequence↗
PubMed2000

Test of molecular mode coupling theory for general rigid molecules

We report recent progress on the test of mode coupling theory for molecular liquids (MMCT) for molecules of arbitrary shape. The MMCT equations in the long time limit are solved for supercooled water including all molecular degrees of freedom. In contrast to our earlier treatment of water as a linear molecule, we find that the glass-transition temperature T(c) is overestimated by the theory as was found in the case of simple liquids. The nonergodicity parameters are calculated from the "full" set of MMCT equations truncated at l(co)=2. These results are compared (i) with the nonergodicity parameters from MMCT with l(co)=2 in the "dipole" approximation n=n(')=0 and the diagonalization approximation n=n(')=0, l=l(') and (ii) with the corresponding results from a MD simulation. This work supports the possibility that a reduction to the most prominent correlators may constitute a valid approximation for solving the MMCT equations for rigid molecules.

Journal Article↗
PubMed2000

Effects of dispersion forces on the structure and thermodynamics of fluid krypton

Semianalytical and numerical calculations are performed to predict the structural and thermodynamic properties of low-density Kr fluid. Assuming that the interatomic forces can be modelled by a pairwise potential plus the three-body Axilrod-Teller potential, two different routes are explored. The first one is based on the hybridized mean spherical approximation integral equation of the theory of liquids and the second one uses large-scale molecular dynamics (MD). Algorithms for MD simulation are constructed on parallel machines to reduce the amount of computer time induced by the calculations of the three-body forces and the pair-correlation function. Our results obtained with the two methods mentioned above are in quite good agreement with the recent small-angle neutron-scattering experiments [Formisano et al., Phys. Rev. Lett. 79, 221 (1997); Benmore et al., J. Phys.: Condens. Matter 11, 3091 (1999)]. Moreover, the reliability of the asymptotic form of the integral equation is assessed for the specific case of dispersion forces including the three-body contributions, by an analysis at low wave vector and low density. It is seen that the effects of the Axilrod-Teller triple-dipole potential cannot be ignored to describe the structure and the thermodynamic properties of fluid krypton even at low density.

Journal Article↗
PubMed2000

Microscopic expressions for the thermodynamic temperature

We show that arbitrary phase space vector fields can be used to generate phase functions whose ensemble averages give the thermodynamic temperature. We describe conditions for the validity of these functions in periodic boundary systems and the molecular dynamics (MD) ensemble, and test them with a short-ranged potential MD simulation.

Journal Article↗
PubMed2003

Freezing transition and correlated motion in a quasi-two-dimensional colloid suspension.

Recent experiments have demonstrated that the deviation of the single-particle displacement distribution from Gaussian form in a dense quasi-two-dimensional colloid suspension is a result of heterogenous dynamics that involves cooperative motions of neighboring colloid particles [J. Chem. Phys. 47, 9142 (2001)]. In this paper, we report the results of molecular dynamics (MD) simulations of a quasi-two-dimensional assembly of nearly hard-sphere colloid particles. The colloid-colloid interaction we use is short ranged and everywhere repulsive; it is related to the Marcus-Rice (MR) and modified MR interactions used in a previous study [Phys. Rev. E 58, 7529 (1998)]. As is the case for those systems, the one we study supports liquid, hexatic, and solid phases. Our calculations show that the deviation of the single-particle displacement distribution from Gaussian form is present in the liquid phase, and that a sharp increase in its magnitude occurs at the liquidus density and extends into the crystalline phase. For densities greater than the liquidus density we find three dynamical relaxation processes that include, at intermediate times, a slowing down in the rate of growth of the diffusive displacement of a particle due to the cage effect. As the density increases toward the solidus density, the dependence of the mean squared displacement on time, at intermediate times, changes from sublinear to zero. The onset of the long-time relaxation mode corresponds to the time at which the deviation of the particle displacement distribution from Gaussian form is a maximum. At this time, which increases exponentially with the density, the self-part of the van Hove function exhibits multiple maxima with respect to r while the distinct part of the van Hove function is a maximum at the origin, thereby signaling jump dynamics. At long times the particle mean square displacement has diffusive character at all densities including solid phase densities. A remarkable feature of our findings is the continuity of character of the particle displacement from the liquid phase through the hexatic phase and into the solid phase. Cooperative jumps that lead to diffusive process in crystals can be explained by a mechanism that involves many such correlated hops in random locations and random directions (but along the crystallographic axes) thereby generating effective random walk behavior. We argue that the collective motion we have found is generated by superpositions of instantaneous normal mode vibrations along diffusive paths. The diffusive paths are along the directions with strong bond orientation correlation, and start to grow in amplitude rapidly on entry into the hexatic phase.

Biophysical Phenomena↗
PubMed2006

Critical behavior of entropic shear rigidity.

We report on extensive molecular dynamics (MD) simulations of a model for gels in both two and three dimensions. The model consists of randomly cross-linked monomers with a concentration p of cross-links above the percolation concentration so that the system is in the amorphous solid phase. As the concentration of cross-links approaches the percolation concentration, the entropic shear modulus vanishes as G approximately(p-p(c))(t) with t approximately 1.9 in three dimensions and t approximately 1.3 in two dimensions. These results hold whether or not the background fluid consisting of finite clusters is retained in the system. These results are also consistent with our previous calculations and with a conjecture of de Gennes but not with recent analytical results and another body of simulations.

Journal Article↗
PubMed2002

Influence of carbon curvature on molecular adsorptions in carbon-based materials: a force field approach.

A general force field methodology is developed for description of molecular interactions in carbon-based materials. The method makes use of existing parameters of potential functions developed for sp(2) and sp(3) carbons and allows accurate representation of molecular forces in curved carbon environment. The potential parameters are explicitly curvature and site dependent. The proposed force field approach was used in molecular dynamics (MD) simulations for hydrogen adsorption in single-walled carbon nanotubes (SWNTs). The results reveal significant nanotube deformations and the calculated energies of adsorption are comparable to the reported experimental heat of adsorption for H2 in SWNTs.

Journal Article↗
PubMed2002

Patching and tearing single-wall carbon-nanotube ropes into multiwall carbon nanotubes.

Bundles of single-wall carbon nanotubes (SWCNTs) coalesce forming multiwall carbon nanotubes (MWCNTs), containing from two to six nested tubes, under thermal treatment at high temperatures [(2200-2400) degrees C]. This structural transformation is confirmed by extensive molecular dynamics (MD) simulations. The simulations suggest a "patching-and-tearing" mechanism for the single-wall-to-multiwall transformation underlying the "concerted" coalescence of the tubes that begins with their polymerization. Tubes of different sizes and chiralities are considered.

Journal Article↗
PubMed2003

Shock waves in high-energy materials: the initial chemical events in nitramine RDX.

We extend the reactive force field ReaxFF to describe the high energy nitramine RDX and use it with molecular dynamics (MD) to study its shock-induced chemistry. We studied shock propagation via nonequilibrium MD simulations at various collision velocities. We find that for high impact velocities (>6 km/s) the RDX molecules decompose and react to form a variety of small molecules in very short time scales (<3 ps). These products are consistent with those found experimentally at longer times. For lower velocities only NO2 is formed, also in agreement with experiments.

Aniline Compounds↗
PubMed2005

Ultrafast core-hole-induced dynamics in water probed by x-ray emission spectroscopy.

The isotope effect and excitation-energy dependence have been measured in the oxygen K-edge x-ray emission spectrum (XES). The use of XES to monitor core decay processes provides information about molecular dynamics (MD) on an ultrafast time scale through the O1s lifetime of a few femtoseconds. Different nuclear masses give rise to differences in the dynamics and the observed isotope effect in XES is direct evidence of the importance of such processes. MD simulations show that even the excitation-energy dependence in the XES is mainly related to differences in core-excited-state dynamics.

Journal Article↗
PubMed2005

Coarse nonlinear dynamics and metastability of filling-emptying transitions: water in carbon nanotubes.

Using a coarse-grained molecular dynamics (CMD) approach we study the apparent nonlinear dynamics of water molecules filling or emptying carbon nanotubes as a function of system parameters. Different levels of the pore hydrophobicity give rise to tubes that are empty, water-filled, or fluctuate between these two long-lived metastable states. The corresponding coarse-grained free-energy surfaces and their hysteretic parameter dependence are explored by linking MD to continuum fixed point and bifurcation algorithms. The results are validated through equilibrium MD simulations.

Algorithms↗
PubMed2001

Hydration of some large and highly charged metal ions.

EXAFS studies of metal ions with hydration numbers higher than six in aqueous solution, often show asymmetric distribution of the metal-oxygen bond distances. The hydration number can be determined from a correlation with the bond distance. The mean Ca-O distance 2.46(1) A shows the calcium(II) ion to be eight-hydrated in a wide asymmetric distribution. Theoretically calculated EXAFS oscillations for individual snapshots from an MD simulation show large variations. The scandium(III) ion is surrounded by two groups of about eight water molecules, with the mean Sc-O distance 2.185(6) A. The yttrium(III) ion coordinates eight waters in an asymmetric distribution at 2.368(5) A, and the lanthanum(III) ion 6 + 3 water molecules at 2.52(2) and 2.65(3) A, respectively. For the the uranium(IV) and thorium(IV) ions, the M-O distances 2.42(1) and 2.45(1) A, respectively, indicate hydration numbers close to 10.

Cations↗
PubMed2006

Engineered lanthanide-binding metallohomeodomains: designing folded chimeras by modular turn substitution.

A series of chimeric metallohomeodomains are described, engineered by rational design of a flexible Ca/Ln binding site into a DNA-binding scaffold. A modular turn-substitution approach was used to create proteins that both bind DNA and lanthanide ions, while retaining the secondary structure of the full homeodomain (determined by circular dichroism [CD]). Four similar metallohomeodomains were designed (C1-C4), their structural stability predicted by molecular dynamics (MD) simulation of loop-mutations into the known homeodomain structure, and each designed protein cloned, expressed, and purified using standard molecular biology techniques. Two of the four loop insertions resulted in folded, metal- and DNA-binding proteins (EuC2 Kd = 2.1 +/- 0.4 microM; EuC4 Kd = 3.2 +/- 1.0 microM). These results show the successful incorporation of a metal site into a full protein domain, without compromising long-range structure. This is an important achievement in biomolecular design, as it provides a critical starting point for exploring metallonuclease function and substrate accessibility in a well-organized chimeric protein domain (rather than only in small HTH peptide systems).

Amino Acid Motifs↗