Calculation of absolute protein-ligand binding free energy using distributed replica sampling.

Distributed replica sampling [T. Rodinger et al., J. Chem. Theory Comput. 2, 725 (2006)] is a simple and general scheme for Boltzmann sampling of conformational space by computer simulation in which multiple replicas of the system undergo a random walk in reaction coordinate or temperature space. Individual replicas are linked through a generalized Hamiltonian containing an extra potential energy term or bias which depends on the distribution of all replicas, thus enforcing the desired sampling distribution along the coordinate or parameter of interest regardless of free energy barriers. In contrast to replica exchange methods, efficient implementation of the algorithm does not require synchronicity of the individual simulations. The algorithm is inherently suited for large-scale simulations using shared or heterogeneous computing platforms such as a distributed network. In this work, we build on our original algorithm by introducing Boltzmann-weighted jumping, which allows moves of a larger magnitude and thus enhances sampling efficiency along the reaction coordinate. The approach is demonstrated using a realistic and biologically relevant application; we calculate the standard binding free energy of benzene to the L99A mutant of T4 lysozyme. Distributed replica sampling is used in conjunction with thermodynamic integration to compute the potential of mean force for extracting the ligand from protein and solvent along a nonphysical spatial coordinate. Dynamic treatment of the reaction coordinate leads to faster statistical convergence of the potential of mean force than a conventional static coordinate, which suffers from slow transitions on a rugged potential energy surface.

[1]  Régis Pomès,et al.  Enhancing the accuracy, the efficiency and the scope of free energy simulations. , 2005, Current opinion in structural biology.

[2]  Justin L. MacCallum,et al.  Calculation of the water–cyclohexane transfer free energies of neutral amino acid side‐chain analogs using the OPLS all‐atom force field , 2003, J. Comput. Chem..

[3]  Laxmikant V. Kale,et al.  NAMD2: Greater Scalability for Parallel Molecular Dynamics , 1999 .

[4]  I. Andricioaei Specialized Methods for Improving Ergodic Sampling Using Molecular Dynamics and Monte Carlo Simulations , 2007 .

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

[6]  S Gnanakaran,et al.  Peptide folding simulations. , 2003, Current opinion in structural biology.

[7]  Michael R. Shirts,et al.  Extremely precise free energy calculations of amino acid side chain analogs: Comparison of common molecular mechanics force fields for proteins , 2003 .

[8]  Anna-Pitschna E. Kunz,et al.  A comparison of methods to compute the potential of mean force. , 2007, Chemphyschem : a European journal of chemical physics and physical chemistry.

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

[10]  C. Chipot,et al.  Overcoming free energy barriers using unconstrained molecular dynamics simulations. , 2004, The Journal of chemical physics.

[11]  M. Gilson,et al.  The statistical-thermodynamic basis for computation of binding affinities: a critical review. , 1997, Biophysical journal.

[12]  J. Hine,et al.  Structural effects on rates and equilibriums. XIX. Intrinsic hydrophilic character of organic compounds. Correlations in terms of structural contributions , 1975 .

[13]  Christophe Chipot,et al.  Free Energy Calculations. The Long and Winding Gilded Road , 2002 .

[14]  Thomas Simonson,et al.  Free energy simulations come of age: protein-ligand recognition. , 2002, Accounts of chemical research.

[15]  M. Shirts,et al.  The solvation interface is a determining factor in peptide conformational preferences. , 2006, Journal of molecular biology.

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

[17]  S. Whittington,et al.  Monte carlo study of the interacting self-avoiding walk model in three dimensions , 1996 .

[18]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

[19]  S. Joseph,et al.  Simulating movement of tRNA into the ribosome during decoding. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Régis Pomès,et al.  Absolute free energy calculations by thermodynamic integration in four spatial dimensions. , 2005, The Journal of chemical physics.

[21]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[22]  Lu Wang,et al.  Inclusion of Loss of Translational and Rotational Freedom in Theoretical Estimates of Free Energies of Binding. Application to a Complex of Benzene and Mutant T4 Lysozyme , 1997 .

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

[24]  A. Garcia,et al.  Reversible temperature and pressure denaturation of a protein fragment: a replica exchange molecular dynamics simulation study. , 2004, Physical review letters.

[25]  R. Pomès,et al.  Equilibrium exchange enhances the convergence rate of umbrella sampling , 2008 .

[26]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[27]  Christian Bartels,et al.  Multidimensional adaptive umbrella sampling: Applications to main chain and side chain peptide conformations , 1997 .

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

[29]  A. Pohorille,et al.  Free energy calculations : theory and applications in chemistry and biology , 2007 .

[30]  A. Lyubartsev,et al.  New approach to Monte Carlo calculation of the free energy: Method of expanded ensembles , 1992 .

[31]  Yuji Sugita,et al.  Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape , 2000, cond-mat/0009119.

[32]  K. Dill,et al.  On the use of orientational restraints and symmetry corrections in alchemical free energy calculations. , 2006, The Journal of chemical physics.

[33]  Yuko Okamoto,et al.  Replica-exchange extensions of simulated tempering method. , 2004, The Journal of chemical physics.

[34]  M. Mezei Adaptive umbrella sampling: Self-consistent determination of the non-Boltzmann bias , 1987 .

[35]  Gerhard Hummer,et al.  Cooperative water filling of a nonpolar protein cavity observed by high-pressure crystallography and simulation. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[36]  J. Onuchic,et al.  Folding a protein in a computer: An atomic description of the folding/unfolding of protein A , 2003, Proceedings of the National Academy of Sciences of the United States of America.

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

[38]  Wei Yang,et al.  Free energy simulations: use of reverse cumulative averaging to determine the equilibrated region and the time required for convergence. , 2004, The Journal of chemical physics.

[39]  Régis Pomès,et al.  Distributed Replica Sampling. , 2006, Journal of chemical theory and computation.

[40]  R. Zwanzig High‐Temperature Equation of State by a Perturbation Method. I. Nonpolar Gases , 1954 .

[41]  H. Nymeyer,et al.  Simulation of the folding equilibrium of α-helical peptides: A comparison of the generalized Born approximation with explicit solvent , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[42]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

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

[44]  Peter L. Freddolino,et al.  Molecular dynamics simulations of the complete satellite tobacco mosaic virus. , 2006, Structure.

[45]  Michael R Shirts,et al.  Parallelized-over-parts computation of absolute binding free energy with docking and molecular dynamics. , 2006, The Journal of chemical physics.

[46]  J A McCammon,et al.  Theory of biomolecular recognition. , 1998, Current opinion in structural biology.

[47]  Benoît Roux,et al.  Calculating excess chemical potentials using dynamic simulations in the fourth dimension , 1999 .

[48]  Wilfred F. van Gunsteren,et al.  Computation of free energy , 2002 .

[49]  Yuko Okamoto,et al.  Replica-exchange multicanonical and multicanonical replica-exchange Monte Carlo simulations of peptides. II. Application to a more complex system , 2003 .

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

[51]  Eric F Darve,et al.  Calculating free energies using average force , 2001 .

[52]  Y. Sugita,et al.  Replica-exchange multicanonical and multicanonical replica-exchange Monte Carlo simulations of peptides. I. Formulation and benchmark test , 2003 .

[53]  B. Brandsdal,et al.  Evaluation of protein-protein association energies by free energy perturbation calculations. , 2000, Protein engineering.

[54]  B. Matthews,et al.  Energetic origins of specificity of ligand binding in an interior nonpolar cavity of T4 lysozyme. , 1995, Biochemistry.

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

[56]  Stefan Boresch,et al.  Absolute Binding Free Energies: A Quantitative Approach for Their Calculation , 2003 .

[57]  Y. Okamoto Institute for Molecular Science,et al.  Replica-exchange simulated tempering method for simulations of frustrated systems , 2000 .