RepEx: A Flexible Framework for Scalable Replica Exchange Molecular Dynamics Simulations

Replica Exchange (RE) simulations have emerged as an important algorithmic tool for the molecular sciences. Typically RE functionality is integrated into the molecular simulation software package. A primary motivation of the tight integration of RE functionality with simulation codes has been performance. This is limiting at multiple levels. First, advances in the RE methodology are tied to the molecular simulation code for which they were developed. Second, it is difficult to extend or experiment with novel RE algorithms, since expertise in the molecular simulation code is required. The tight integration results in difficulty to gracefully handle failures, and other runtime fragilities. We propose the RepEx framework which is addressing aforementioned shortcomings, while striking the balance between flexibility (any RE scheme) and scalability (several thousand replicas) over a diverse range of HPC platforms. The primary contributions of the RepEx framework are: (i) its ability to support different Replica Exchange schemes independent of molecular simulation codes, (ii) provide the ability to execute different exchange schemes and replica counts independent of the specific availability of resources, (iii) provide a runtime system that has first-class support for task-level parallelism, and (iv) provide a required scalability along multiple dimensions.

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

[2]  Preston Moore,et al.  Metadynamics as a tool for exploring free energy landscapes of chemical reactions. , 2006, Accounts of chemical research.

[3]  Adrian E. Roitberg,et al.  Multidimensional Replica Exchange Molecular Dynamics Yields a Converged Ensemble of an RNA Tetranucleotide , 2013, Journal of chemical theory and computation.

[4]  Wei Dai,et al.  A framework for flexible and scalable replica-exchange on production distributed CI , 2013, XSEDE.

[5]  Phil Andrews,et al.  Project Summary: XSEDE: eXtreme Science and Engineering Discovery Environment , 2010 .

[6]  Bin W. Zhang,et al.  Large‐Scale Asynchronous and Distributed Multidimensional Replica Exchange Molecular Simulations and Efficiency Analysis , 2016, Journal of computational chemistry.

[7]  Maria T. Panteva,et al.  Multiscale methods for computational RNA enzymology. , 2015, Methods in enzymology.

[8]  Shantenu Jha,et al.  RADICAL-Pilot: Scalable Execution of Heterogeneous and Dynamic Workloads on Supercomputers , 2015, ArXiv.

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

[10]  Peter M. Kasson,et al.  GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit , 2013, Bioinform..

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

[12]  J. L. Suter,et al.  Structure, dynamics, and function of the hammerhead ribozyme in bulk water and at a clay mineral surface from replica exchange molecular dynamics. , 2015, Langmuir : the ACS journal of surfaces and colloids.

[13]  Takashi Imai,et al.  Temperature and pressure effects on conformational equilibria of alanine dipeptide in aqueous solution. , 2004, Biopolymers.

[14]  J. Mccammon,et al.  Population Based Reweighting of Scaled Molecular Dynamics , 2013, The journal of physical chemistry. B.

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

[16]  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.

[17]  Guohui Li,et al.  Free energy simulations with the AMOEBA polarizable force field and metadynamics on GPU platform , 2016, J. Comput. Chem..

[18]  Lei Huang,et al.  Generalized scalable multiple copy algorithms for molecular dynamics simulations in NAMD , 2014, Comput. Phys. Commun..

[19]  Darrin M York,et al.  Interpretation of pH-activity profiles for acid-base catalysis from molecular simulations. , 2015, Biochemistry.

[20]  Shantenu Jha,et al.  Characterization of the three-dimensional free energy manifold for the uracil ribonucleoside from asynchronous replica exchange simulations. , 2015, Journal of chemical theory and computation.

[21]  Tai-Sung Lee,et al.  A New Maximum Likelihood Approach for Free Energy Profile Construction from Molecular Simulations. , 2013, Journal of chemical theory and computation.

[22]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[23]  Tai-Sung Lee,et al.  Roadmaps through free energy landscapes calculated using the multi-dimensional vFEP approach. , 2014, Journal of chemical theory and computation.

[24]  Eric Vanden-Eijnden,et al.  Some recent techniques for free energy calculations , 2009, J. Comput. Chem..

[25]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[26]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[27]  Adrian E Roitberg,et al.  Constant pH replica exchange molecular dynamics in biomolecules using a discrete protonation model. , 2010, Journal of chemical theory and computation.

[28]  Richard A. Friesner,et al.  Integrated Modeling Program, Applied Chemical Theory (IMPACT) , 2005, J. Comput. Chem..

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