Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer.

An extremely scalable computational strategy is described for calculations of the potential of mean force (PMF) in multidimensions on massively distributed supercomputers. The approach involves coupling thousands of umbrella sampling (US) simulation windows distributed to cover the space of order parameters with a Hamiltonian molecular dynamics replica-exchange (H-REMD) algorithm to enhance the sampling of each simulation. In the present application, US/H-REMD is carried out in a two-dimensional (2D) space and exchanges are attempted alternatively along the two axes corresponding to the two order parameters. The US/H-REMD strategy is implemented on the basis of parallel/parallel multiple copy protocol at the MPI level, and therefore can fully exploit computing power of large-scale supercomputers. Here the novel technique is illustrated using the leadership supercomputer IBM Blue Gene/P with an application to a typical biomolecular calculation of general interest, namely the binding of calcium ions to the small protein Calbindin D9k. The free energy landscape associated with two order parameters, the distance between the ion and its binding pocket and the root-mean-square deviation (rmsd) of the binding pocket relative the crystal structure, was calculated using the US/H-REMD method. The results are then used to estimate the absolute binding free energy of calcium ion to Calbindin D9k. The tests demonstrate that the 2D US/H-REMD scheme greatly accelerates the configurational sampling of the binding pocket, thereby improving the convergence of the potential of mean force calculation.

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

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

[3]  Alexander D. MacKerell,et al.  Development and current status of the CHARMM force field for nucleic acids , 2000, Biopolymers.

[4]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

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

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

[7]  S Linse,et al.  Symmetrical stabilization of bound Ca2+ ions in a cooperative pair of EF-hands through hydrogen bonding of coordinating water molecules in calbindin D(9k). , 2001, Biochemistry.

[8]  Wonpil Im,et al.  Transmembrane helix assembly by window exchange umbrella sampling. , 2012, Physical review letters.

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

[10]  Eric Darve,et al.  Calculating Free Energies Using a Scaled-Force Molecular Dynamics Algorithm , 2002 .

[11]  B. Roux,et al.  Calculation of absolute protein-ligand binding free energy from computer simulations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Daniel M Zuckerman,et al.  Equilibrium sampling in biomolecular simulations. , 2011, Annual review of biophysics.

[13]  Sara Linse,et al.  Coupling of ligand binding and dimerization of helix‐loop‐helix peptides: Spectroscopic and sedimentation analyses of calbindin D9k EF‐hands , 2002, Proteins.

[14]  Eric Darve,et al.  Assessing the efficiency of free energy calculation methods. , 2004, The Journal of chemical physics.

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

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

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

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

[19]  S. Rick,et al.  Increasing the Efficiency of Free Energy Calculations Using Parallel Tempering and Histogram Reweighting. , 2006, Journal of chemical theory and computation.

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

[21]  Christopher B. Barnett,et al.  Free Energies from Adaptive Reaction Coordinate Forces (FEARCF): an application to ring puckering , 2009 .

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

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

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

[25]  S. Forsén,et al.  Proline cis-trans isomers in calbindin D9k observed by X-ray crystallography. , 1992, Journal of Molecular Biology.

[26]  John Blankenship,et al.  Site–site communication in the EF-hand Ca2+-binding protein calbindin D9k , 2000, Nature Structural Biology.

[27]  Benoît Roux,et al.  Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations , 2001 .

[28]  B. Roux,et al.  Molecular dynamics study of calbindin D9k in the apo and singly and doubly calcium‐loaded states , 1998, Proteins.

[29]  E. Vanden-Eijnden,et al.  Single-sweep methods for free energy calculations. , 2007, The Journal of chemical physics.

[30]  T. Grundström,et al.  Electrostatic contributions to the binding of Ca2+ in calbindin D9k. , 1991, Biochemistry.

[31]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

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

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