Constant pressure molecular dynamics algorithms
01 Sep 1994-Journal of Chemical Physics (American Institute of Physics)-Vol. 101, Iss: 5, pp 4177-4189
TL;DR: In this paper, a modularly invariant equations of motion are derived that generate the isothermal-isobaric ensemble as their phase space averages, and the resulting methods are tested on two problems, a particle in a one-dimensional periodic potential and a spherical model of C60 in the solid/fluid phase.
Abstract: Modularly invariant equations of motion are derived that generate the isothermal–isobaric ensemble as their phase space averages. Isotropic volume fluctuations and fully flexible simulation cells as well as a hybrid scheme that naturally combines the two motions are considered. The resulting methods are tested on two problems, a particle in a one‐dimensional periodic potential and a spherical model of C60 in the solid/fluid phase.
Citations
More filters
TL;DR: In this paper, a new method for performing molecular dynamics simulations under constant pressure is presented, which is based on the extended system formalism introduced by Andersen, the deterministic equations of motion for the piston degree of freedom are replaced by a Langevin equation; a suitable choice of collision frequency then eliminates the unphysical "ringing" of the volume associated with the piston mass.
Abstract: A new method for performing molecular dynamics simulations under constant pressure is presented. In the method, which is based on the extended system formalism introduced by Andersen, the deterministic equations of motion for the piston degree of freedom are replaced by a Langevin equation; a suitable choice of collision frequency then eliminates the unphysical ‘‘ringing’’ of the volume associated with the piston mass. In this way it is similar to the ‘‘weak coupling algorithm’’ developed by Berendsen and co‐workers to perform molecular dynamics simulation without piston mass effects. It is shown, however, that the weak coupling algorithm induces artifacts into the simulation which can be quite severe for inhomogeneous systems such as aqueous biopolymers or liquid/liquid interfaces.
3,799 citations
TL;DR: The presented lipid FF is developed and applied to phospholipid bilayers with both choline and ethanolamine containing head groups and with both saturated and unsaturated aliphatic chains and is anticipated to be of utility for simulations of pure lipid systems as well as heterogeneous systems including membrane proteins.
Abstract: A significant modification to the additive all-atom CHARMM lipid force field (FF) is developed and applied to phospholipid bilayers with both choline and ethanolamine containing head groups and with both saturated and unsaturated aliphatic chains. Motivated by the current CHARMM lipid FF (C27 and C27r) systematically yielding values of the surface area per lipid that are smaller than experimental estimates and gel-like structures of bilayers well above the gel transition temperature, selected torsional, Lennard-Jones and partial atomic charge parameters were modified by targeting both quantum mechanical (QM) and experimental data. QM calculations ranging from high-level ab initio calculations on small molecules to semiempirical QM studies on a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer in combination with experimental thermodynamic data were used as target data for parameter optimization. These changes were tested with simulations of pure bilayers at high hydration of the following six lipids: ...
3,489 citations
TL;DR: The optimal simulation protocol for each program has been implemented in CHARMM-GUI and is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.
Abstract: Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order para...
2,182 citations
TL;DR: Simulation of the folding of a WW domain showed a well-defined folding pathway and simulation of the dynamics of bovine pancreatic trypsin inhibitor showed interconversion between distinct conformational states.
Abstract: Molecular dynamics (MD) simulations are widely used to study protein motions at an atomic level of detail, but they have been limited to time scales shorter than those of many biologically critical conformational changes. We examined two fundamental processes in protein dynamics—protein folding and conformational change within the folded state—by means of extremely long all-atom MD simulations conducted on a special-purpose machine. Equilibrium simulations of a WW protein domain captured multiple folding and unfolding events that consistently follow a well-defined folding pathway; separate simulations of the protein’s constituent substructures shed light on possible determinants of this pathway. A 1-millisecond simulation of the folded protein BPTI reveals a small number of structurally distinct conformational states whose reversible interconversion is slower than local relaxations within those states by a factor of more than 1000.
1,650 citations
Cites methods from "Constant pressure molecular dynamic..."
...We removed the restraints and performed a 6-ns simulation at constant isotropic pressure of 1 atm and temperature of 300 K, using the MTK algorithm (34) with 10-ps relaxation time as implemented in the software Desmond (35)....
[...]
TL;DR: The new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems are described, including addition of new lipid types, including phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol.
Abstract: CHARMM-GUI Membrane Builder, http://www.charmm-gui.org/input/membrane, is a web-based user interface designed to interactively build all-atom protein/membrane or membrane-only systems for molecular dynamics simulations through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types, such as phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance-based algorithm for faster initial ion displacement, (5) CHARMM inputs for P21 image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: CL lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM-GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments. © 2014 Wiley Periodicals, Inc.
1,604 citations
References
More filters
Book•
11 Feb 1988TL;DR: In this paper, the gear predictor -corrector is used to calculate forces and torques in a non-equilibrium molecular dynamics simulation using Monte Carlo methods. But it is not suitable for the gear prediction problem.
Abstract: Introduction Statistical mechanics Molecular dynamics Monte Carlo methods Some tricks of the trade How to analyse the results Advanced simulation techniques Non-equilibrium molecular dynamics Brownian dynamics Quantum simulations Some applications Appendix A: Computers and computer simulation Appendix B: Reduced units Appendix C: Calculation of forces and torques Appendix D: Fourier transforms Appendix E: The gear predictor - corrector Appendix F: Programs on microfiche Appendix G: Random numbers References Index.
21,073 citations
TL;DR: In this paper, a numerical algorithm integrating the 3N Cartesian equations of motion of a system of N points subject to holonomic constraints is formulated, and the relations of constraint remain perfectly fulfilled at each step of the trajectory despite the approximate character of numerical integration.
18,394 citations
TL;DR: The dynamical steady-state probability density is found in an extended phase space with variables x, p/sub x/, V, epsilon-dot, and zeta, where the x are reduced distances and the two variables epsilus-dot andZeta act as thermodynamic friction coefficients.
Abstract: Nos\'e has modified Newtonian dynamics so as to reproduce both the canonical and the isothermal-isobaric probability densities in the phase space of an N-body system. He did this by scaling time (with s) and distance (with ${V}^{1/D}$ in D dimensions) through Lagrangian equations of motion. The dynamical equations describe the evolution of these two scaling variables and their two conjugate momenta ${p}_{s}$ and ${p}_{v}$. Here we develop a slightly different set of equations, free of time scaling. We find the dynamical steady-state probability density in an extended phase space with variables x, ${p}_{x}$, V, \ensuremath{\epsilon}\ifmmode \dot{}\else \.{}\fi{}, and \ensuremath{\zeta}, where the x are reduced distances and the two variables \ensuremath{\epsilon}\ifmmode \dot{}\else \.{}\fi{} and \ensuremath{\zeta} act as thermodynamic friction coefficients. We find that these friction coefficients have Gaussian distributions. From the distributions the extent of small-system non-Newtonian behavior can be estimated. We illustrate the dynamical equations by considering their application to the simplest possible case, a one-dimensional classical harmonic oscillator.
17,939 citations
TL;DR: In this article, the authors compared the canonical distribution in both momentum and coordinate space with three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 77, 63 (1983); and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983).
Abstract: Three recently proposed constant temperature molecular dynamics methods by: (i) Nose (Mol. Phys., to be published); (ii) Hoover et al. [Phys. Rev. Lett. 48, 1818 (1982)], and Evans and Morriss [Chem. Phys. 77, 63 (1983)]; and (iii) Haile and Gupta [J. Chem. Phys. 79, 3067 (1983)] are examined analytically via calculating the equilibrium distribution functions and comparing them with that of the canonical ensemble. Except for effects due to momentum and angular momentum conservation, method (1) yields the rigorous canonical distribution in both momentum and coordinate space. Method (2) can be made rigorous in coordinate space, and can be derived from method (1) by imposing a specific constraint. Method (3) is not rigorous and gives a deviation of order N−1/2 from the canonical distribution (N the number of particles). The results for the constant temperature–constant pressure ensemble are similar to the canonical ensemble case.
13,921 citations
Book•
01 Jan 1974
TL;DR: In this paper, Newtonian mechanics: experimental facts investigation of the equations of motion, variational principles Lagrangian mechanics on manifolds oscillations rigid bodies, differential forms symplectic manifolds canonical formalism introduction to pertubation theory.
Abstract: Part 1 Newtonian mechanics: experimental facts investigation of the equations of motion. Part 2 Lagrangian mechanics: variational principles Lagrangian mechanics on manifolds oscillations rigid bodies. Part 3 Hamiltonian mechanics: differential forms symplectic manifolds canonical formalism introduction to pertubation theory.
11,008 citations