skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis

Journal Article · · Journal of Physical Chemistry. B, Condensed Matter, Materials, Surfaces, Interfaces and Biophysical Chemistry

We report that replica-exchange molecular dynamics (REMD) has been proven to efficiently improve the convergence of free energy perturbation (FEP) calculations involving considerable reorganization of their surrounding. We previously introduced the FEP/(λ,H)-REMD algorithm for ligand binding, in which replicas along the alchemical thermodynamic coupling axis λ were expanded as a series of Hamiltonian boosted replicas along a second axis to form a two-dimensional (2D) replica-exchange exchange map [Jiang, W.; Roux, B., J. Chem. Theory Comput. 2010, 6 (9), 2559-2565]. Aiming to achieve a similar performance at a lower computational cost, we propose here a modified version of this algorithm in which only the end-states along the alchemical axis are augmented by boosted replicas. The reduced FEP/(λ,H)-REMD method with one-dimensional (1D) unbiased alchemical thermodynamic coupling axis λ is implemented on the basis of generic multiple copy algorithm (MCA) module of the biomolecular simulation program NAMD. The flexible MCA framework of NAMD enables a user to design customized replica-exchange patterns through Tcl scripting in the context of a highly parallelized simulation program without touching the source code. Two Hamiltonian tempering boosting scheme were examined with the new algorithm: a first one based on potential energy rescaling of a pre-identified “solute”, and a second one via the introduction of flattening torsional free energy barriers. As two illustrative examples with reliable experiment data, the absolute binding free energies of pxylene and n-butylbenzene to the nonpolar cavity of the L99A mutant of T4 lysozyme were calculated. Lastly, the tests demonstrate that the new protocol efficiently enhances the sampling of torsional motions for backbone and side chains around the binding pocket and accelerates the convergence of the free energy computations.

Research Organization:
Argonne National Lab. (ANL), Argonne, IL (United States)
Sponsoring Organization:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR); National Institutes of Health (NIH)
Grant/Contract Number:
AC02-06CH11357
OSTI ID:
1488538
Journal Information:
Journal of Physical Chemistry. B, Condensed Matter, Materials, Surfaces, Interfaces and Biophysical Chemistry, Vol. 122, Issue 41; ISSN 1520-6106
Publisher:
American Chemical SocietyCopyright Statement
Country of Publication:
United States
Language:
English
Citation Metrics:
Cited by: 29 works
Citation information provided by
Web of Science

References (45)

Efficient Drug Lead Discovery and Optimization journal June 2009
The Many Roles of Computation in Drug Discovery journal March 2004
The statistical-thermodynamic basis for computation of binding affinities: a critical review journal March 1997
Calculation of Protein-Ligand Binding Affinities journal June 2007
Binding of Small-Molecule Ligands to Proteins: “What You See” Is Not Always “What You Get” journal April 2009
Alchemical free energy methods for drug discovery: progress and challenges journal April 2011
Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations journal March 2001
Calculation of Standard Binding Free Energies:  Aromatic Molecules in the T4 Lysozyme L99A Mutant journal July 2006
Hydration of Amino Acid Side Chains:  Nonpolar and Electrostatic Contributions Calculated from Staged Molecular Dynamics Free Energy Simulations with Explicit Water Molecules journal October 2004
Predicting Absolute Ligand Binding Free Energies to a Simple Model Site journal August 2007
Confine-and-Release Method:  Obtaining Correct Binding Free Energies in the Presence of Protein Conformational Change journal May 2007
Accurate and Reliable Prediction of Relative Ligand Binding Potency in Prospective Drug Discovery by Way of a Modern Free-Energy Calculation Protocol and Force Field journal February 2015
On achieving high accuracy and reliability in the calculation of relative protein-ligand binding affinities journal January 2012
Advancing Drug Discovery through Enhanced Free Energy Calculations journal July 2017
Glycoside Hydrolase Processivity Is Directly Related to Oligosaccharide Binding Free Energy journal December 2013
Nanoparticle–Protein Interactions: A Thermodynamic and Kinetic Study of the Adsorption of Bovine Serum Albumin to Gold Nanoparticle Surfaces journal November 2013
Predicting Adsorption Affinities of Small Molecules on Carbon Nanotubes Using Molecular Dynamics Simulation journal October 2015
The lag between the Hamiltonian and the system configuration in free energy perturbation calculations journal December 1989
Calculation of absolute protein-ligand binding free energy from computer simulations journal May 2005
Absolute Binding Free Energy Calculations Using Molecular Dynamics Simulations with Restraining Potentials journal October 2006
Free Energy Perturbation Hamiltonian Replica-Exchange Molecular Dynamics (FEP/H-REMD) for Absolute Ligand Binding Free Energy Calculations journal July 2010
Binding Energy Distribution Analysis Method: Hamiltonian Replica Exchange with Torsional Flattening for Binding Mode Prediction and Binding Free Energy Estimation journal April 2016
Generalized scalable multiple copy algorithms for molecular dynamics simulations in NAMD journal March 2014
Sensitivity in Binding Free Energies Due to Protein Reorganization journal August 2016
Specificity of ligand binding in a buried nonpolar cavity of T4 lysozyme: Linkage of dynamics and structural plasticity journal July 1995
On the use of orientational restraints and symmetry corrections in alchemical free energy calculations journal August 2006
Computation of Absolute Hydration and Binding Free Energy with Free Energy Perturbation Distributed Replica-Exchange Molecular Dynamics journal August 2009
Free energy calculations for DNA base stacking by replica-exchange umbrella sampling journal February 2004
Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential journal November 2006
Accelerated molecular dynamics: A promising and efficient simulation method for biomolecules journal June 2004
Accelerated molecular dynamics simulations of protein folding journal June 2015
Replica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2) journal August 2011
A generic implementation of replica exchange with solute tempering (REST2) algorithm in NAMD for complex biophysical simulations journal December 2015
Replica-Exchange Accelerated Molecular Dynamics (REXAMD) Applied to Thermodynamic Integration journal September 2008
Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation journal May 2014
Scalable molecular dynamics with NAMD journal January 2005
CHARMM-GUI Ligand Binder for Absolute Binding Free Energy Calculations and Its Application journal December 2012
CHARMM-GUI 10 years for biomolecular modeling and simulation: Biomolecular Modeling and Simulation journal November 2016
Constant pressure molecular dynamics simulation: The Langevin piston method journal September 1995
CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields journal January 2009
Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations journal June 1994
Separation‐shifted scaling, a new scaling method for Lennard‐Jones interactions in thermodynamic integration journal June 1994
Improving the efficiency and reliability of free energy perturbation calculations using overlap sampling methods journal January 2003
THE weighted histogram analysis method for free-energy calculations on biomolecules. I. The method journal October 1992
Generalized solvent boundary potential for computer simulations journal February 2001

Cited By (9)

Ligand binding free energy and kinetics calculation in 2020 journal January 2020
Oversampling Free Energy Perturbation Simulation in Determination of the Ligand‐Binding Free Energy journal December 2019
Deciphering molecular mechanism behind conformational change of the São Paolo metallo-β-lactamase 1 by using enhanced sampling text January 2019
Conformational flexibility correlates with glucose tolerance for point mutations in β-glucosidases – a computational study journal March 2020
On Restraints in End‐Point Protein–Ligand Binding Free Energy Calculations journal December 2019
Deciphering molecular mechanism behind conformational change of the São Paolo metallo-β-lactamase 1 by using enhanced sampling journal December 2019
Conformational flexibility correlates with glucose tolerance for point mutations in β-glucosidases – a computational study text January 2020
On Restraints in End-Point Protein-Ligand Binding Free Energy Calculations journal February 2019
Deciphering molecular mechanism behind conformational change of the São Paolo metallo-β-lactamase 1 by using enhanced sampling text January 2019

Figures / Tables (6)