Generalized scalable multiple copy algorithms for molecular dynamics simulations in NAMD

Computational methodologies that couple the dynamical evolution of a set of replicated copies of a system of interest offer powerful and flexible approaches to characterize complex molecular processes. Such multiple copy algorithms (MCAs) can be used to enhance sampling, compute reversible work and free energies, as well as refine transition pathways. Widely used examples of MCAs include temperature and Hamiltonian-tempering replica-exchange molecular dynamics (T-REMD and H-REMD), alchemical free energy perturbation with lambda replica-exchange (FEP/λ-REMD), umbrella sampling with Hamiltonian replica exchange (US/H-REMD), and string method with swarms-of-trajectories conformational transition pathways. Here, we report a robust and general implementation of MCAs for molecular dynamics (MD) simulations in the highly scalable program NAMD built upon the parallel programming system Charm++. Multiple concurrent NAMD instances are launched with internal partitions of Charm++ and located continuously within a single communication world. Messages between NAMD instances are passed by low-level point-to-point communication functions, which are accessible through NAMD's Tcl scripting interface. The communication-enabled Tcl scripting provides a sustainable application interface for end users to realize generalized MCAs without modifying the source code. Illustrative applications of MCAs with fine-grained inter-copy communication structure, including global lambda exchange in FEP/λ-REMD, window swapping US/H-REMD in multidimensional order parameter space, and string method with swarms-of-trajectories were carried out on IBM Blue Gene/Q to demonstrate the versatility and massive scalability of the present implementation.

[1]  K. Sanbonmatsu,et al.  Structure of Met‐enkephalin in explicit aqueous solution using replica exchange molecular dynamics , 2002, Proteins.

[2]  B. Roux,et al.  Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials. , 2006, Biophysical journal.

[3]  Ross C. Walker,et al.  An overview of the Amber biomolecular simulation package , 2013 .

[4]  B. Berne,et al.  Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[5]  A Mitsutake,et al.  Generalized-ensemble algorithms for molecular simulations of biopolymers. , 2000, Biopolymers.

[6]  Chao Mei,et al.  Optimizing fine-grained communication in a biomolecular simulation application on Cray XK6 , 2012, 2012 International Conference for High Performance Computing, Networking, Storage and Analysis.

[7]  Martin Zacharias,et al.  Advanced replica-exchange sampling to study the flexibility and plasticity of peptides and proteins. , 2013, Biochimica et biophysica acta.

[8]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[9]  E Weinan,et al.  Transition pathways in complex systems: Reaction coordinates, isocommittor surfaces, and transition tubes , 2005 .

[10]  Benoît Roux,et al.  Calculation of Standard Binding Free Energies:  Aromatic Molecules in the T4 Lysozyme L99A Mutant. , 2006, Journal of chemical theory and computation.

[11]  S. Harrison,et al.  Crystal structures of c-Src reveal features of its autoinhibitory mechanism. , 1999, Molecular cell.

[12]  Y. Sugita,et al.  Free energy calculations for DNA base stacking by replica-exchange umbrella sampling , 2004 .

[13]  Aaron R. Dinner,et al.  Erratum: “Nonequilibrium umbrella sampling in spaces of many order parameters” [J. Chem. Phys. 130, 074104 (2009)] , 2012 .

[14]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[15]  H. C. Andersen,et al.  Role of Repulsive Forces in Determining the Equilibrium Structure of Simple Liquids , 1971 .

[16]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[17]  Angel E. Garcia,et al.  Faculty Opinions recommendation of Can a continuum solvent model reproduce the free energy landscape of a beta -hairpin folding in water? , 2002 .

[18]  Benoît Roux,et al.  Publisher’s Note: “Computation of binding free energy with molecular dynamics and grand canonical Monte Carlo simulations” [J. Chem. Phys. 128, 115103 (2008)] , 2008 .

[19]  J. Andrew McCammon,et al.  Replica-Exchange Accelerated Molecular Dynamics (REXAMD) Applied to Thermodynamic Integration , 2008, Journal of chemical theory and computation.

[20]  Yilin Meng,et al.  Self-Learning Adaptive Umbrella Sampling Method for the Determination of Free Energy Landscapes in Multiple Dimensions. , 2013, Journal of chemical theory and computation.

[21]  Benoît Roux,et al.  Computation of Absolute Hydration and Binding Free Energy with Free Energy Perturbation Distributed Replica-Exchange Molecular Dynamics (FEP/REMD). , 2009, Journal of chemical theory and computation.

[22]  P. Cochat,et al.  Et al , 2008, Archives de pediatrie : organe officiel de la Societe francaise de pediatrie.

[23]  Alex Dickson,et al.  Flow-Dependent Unfolding and Refolding of an RNA by Nonequilibrium Umbrella Sampling. , 2011, Journal of chemical theory and computation.

[24]  Martin Zacharias,et al.  Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing‐potential , 2006, Proteins.

[25]  A. Dinner,et al.  Enhanced sampling of nonequilibrium steady states. , 2010, Annual review of physical chemistry.

[26]  E. Stellwagen,et al.  Distribution of Helicity within the Model Peptide Acetyl(AAQAA)3amide , 1994 .

[27]  Benoît Roux,et al.  Hydration of Amino Acid Side Chains: Nonpolar and Electrostatic Contributions Calculated from Staged Molecular Dynamics Free Energy Simulations with Explicit Water Molecules , 2004 .

[28]  P. R. ten Wolde,et al.  Sampling rare switching events in biochemical networks. , 2004, Physical review letters.

[29]  Toby W Allen,et al.  Molecular dynamics - potential of mean force calculations as a tool for understanding ion permeation and selectivity in narrow channels. , 2006, Biophysical chemistry.

[30]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[31]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[32]  Christophe Chipot,et al.  Standard binding free energies from computer simulations: What is the best strategy? , 2013, Journal of chemical theory and computation.

[33]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[34]  Yuko Okamoto,et al.  From multidimensional replica-exchange method to multidimensional multicanonical algorithm and simulated tempering. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[35]  Christopher J. Woods,et al.  Enhanced configurational sampling in binding free-energy calculations , 2003 .

[36]  Benoît Roux,et al.  Atomistic view of the conformational activation of Src kinase using the string method with swarms-of-trajectories. , 2009, Biophysical journal.

[37]  Benoît Roux,et al.  Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer. , 2012, Journal of chemical theory and computation.

[38]  Eric Vanden-Eijnden,et al.  Transition pathways in complex systems: Application of the finite-temperature string method to the alanine dipeptide. , 2005, The Journal of chemical physics.

[39]  J. Kästner Umbrella sampling , 2011 .

[40]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[41]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

[42]  Daniel J. Sindhikara,et al.  Exchange Often and Properly in Replica Exchange Molecular Dynamics. , 2010, Journal of chemical theory and computation.

[43]  Servaas Michielssens,et al.  Improved Replica Exchange Method for Native-State Protein Sampling. , 2011, Journal of chemical theory and computation.

[44]  R. Zhou Free energy landscape of protein folding in water: Explicit vs. implicit solvent , 2003, Proteins.

[45]  Christophe Chipot,et al.  Efficient determination of protein-protein standard binding free energies from first principles. , 2013, Journal of chemical theory and computation.

[46]  Danial S. Dashti,et al.  Computing Alchemical Free Energy Differences with Hamiltonian Replica Exchange Molecular Dynamics (H-REMD) Simulations. , 2011, Journal of chemical theory and computation.

[47]  J. Kirkwood Statistical Mechanics of Fluid Mixtures , 1935 .

[48]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[49]  Sameer Kumar,et al.  Acceleration of an Asynchronous Message Driven Programming Paradigm on IBM Blue Gene/Q , 2013, 2013 IEEE 27th International Symposium on Parallel and Distributed Processing.

[50]  Benoît Roux,et al.  Free Energy Perturbation Hamiltonian Replica-Exchange Molecular Dynamics (FEP/H-REMD) for Absolute Ligand Binding Free Energy Calculations. , 2010, Journal of chemical theory and computation.

[51]  Matthew A. Wyczalkowski,et al.  Satisfying the fluctuation theorem in free-energy calculations with Hamiltonian replica exchange. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[52]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[53]  W. E,et al.  Finite temperature string method for the study of rare events. , 2002, Journal of Physical Chemistry B.

[54]  A. Dinner,et al.  Nonequilibrium umbrella sampling in spaces of many order parameters. , 2009, The Journal of chemical physics.

[55]  Daniel J. Sindhikara,et al.  Exchange frequency in replica exchange molecular dynamics. , 2008, The Journal of chemical physics.

[56]  M. Zacharias,et al.  Simulation of DNA double-strand dissociation and formation during replica-exchange molecular dynamics simulations. , 2009, Physical chemistry chemical physics : PCCP.

[57]  Albert C. Pan,et al.  Finding transition pathways using the string method with swarms of trajectories. , 2008, The journal of physical chemistry. B.