A framework for flexible and scalable replica-exchange on production distributed CI

Replica exchange represents a powerful class of algorithms used for enhanced configurational and energetic sampling in a range of physical systems. Computationally it represents a type of application with multiple scales of communication. At a fine-grained level there is often communication with a replica, typically an MPI process. At a coarse-grained level, the replicas communicate with other replicas -- both temporally as well as in amount of data exchanged. This paper outlines a novel framework developed to support the flexible execution of large-scale replica exchange. The framework is flexible in the sense that it supports different coupling schemes between replicas and is agnostic to the specific underlying simulation -- classical or quantum, serial or parallel simulation. The scalability of the framework is assessed using standard simulation benchmarks. In spite of the increasing communication and coordination requirements as a function of the number of replicas, our framework supports the execution of hundreds replicas without significant overhead. Although there are several specific aspects that will benefit from further optimization, a first working prototype has the ability to fundamentally change the scale of replica exchange simulations possible on production distributed cyberinfrastructure such as XSEDE, as well as support novel usage modes. This paper also represents the release of the framework to the broader biophysical simulation community and provides details on its usage.

[1]  Shantenu Jha,et al.  P∗: A model of pilot-abstractions , 2012, 2012 IEEE 8th International Conference on E-Science.

[2]  Anthony K. Felts,et al.  Free Energy Surfaces of -Hairpin and -Helical Peptides Generated by Replica Exchange Molecular Dynamics with the AGBNP Implicit Solvent Model , 2004 .

[3]  Michael R. Shirts,et al.  Replica exchange and expanded ensemble simulations as Gibbs sampling: simple improvements for enhanced mixing. , 2011, The Journal of chemical physics.

[4]  Ruhong Zhou,et al.  Hydrophobic aided replica exchange: an efficient algorithm for protein folding in explicit solvent. , 2006, The journal of physical chemistry. B.

[5]  A. Laio,et al.  Free-energy landscape for beta hairpin folding from combined parallel tempering and metadynamics. , 2006, Journal of the American Chemical Society.

[6]  R. Levy,et al.  Protein folding pathways from replica exchange simulations and a kinetic network model. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[7]  Emilio Gallicchio,et al.  The Binding Energy Distribution Analysis Method (BEDAM) for the Estimation of Protein-Ligand Binding Affinities. , 2010, Journal of chemical theory and computation.

[8]  Zaida Ann Luthey-Schulten,et al.  Variations in the fast folding rates of the λ-repressor: A hybrid molecular dynamics study , 2004 .

[9]  John Shalf,et al.  SAGA: A Simple API for Grid Applications. High-level application programming on the Grid , 2006 .

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

[11]  Manish Parashar,et al.  Asynchronous replica exchange for molecular simulations , 2008, J. Comput. Chem..

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

[13]  Anthony K. Felts,et al.  Free energy surfaces of beta-hairpin and alpha-helical peptides generated by replica exchange molecular dynamics with the AGBNP implicit solvent model. , 2004, Proteins.

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

[15]  Emilio Gallicchio,et al.  Advances in all atom sampling methods for modeling protein-ligand binding affinities. , 2011, Current opinion in structural biology.

[16]  Ronald M. Levy,et al.  Prediction of SAMPL3 host-guest affinities with the binding energy distribution analysis method (BEDAM) , 2012, Journal of Computer-Aided Molecular Design.

[17]  Anders Wallqvist,et al.  Free-energy profiles of membrane insertion of the M2 transmembrane peptide from influenza A virus. , 2008, Biophysical journal.

[18]  Yuko Okamoto,et al.  Multicanonical Algorithm, Simulated Tempering, Replica-Exchange Method, and All That , 2008, 0810.1117.

[19]  Jonathan W. Essex,et al.  The development of replica-exchange-based free-energy methods , 2003 .

[20]  U H Hansmann,et al.  New Monte Carlo algorithms for protein folding. , 1999, Current opinion in structural biology.

[21]  Emilio Gallicchio,et al.  Conformational Transitions and Convergence of Absolute Binding Free Energy Calculations. , 2012, Journal of chemical theory and computation.

[22]  Shantenu Jha,et al.  SAGA BigJob: An Extensible and Interoperable Pilot-Job Abstraction for Distributed Applications and Systems , 2010, 2010 10th IEEE/ACM International Conference on Cluster, Cloud and Grid Computing.