Free Energy Perturbation Hamiltonian Replica-Exchange Molecular Dynamics (FEP\H-REMD) for absolute ligand binding free energy calculations.
- Biosciences Division
Free Energy Perturbation with Replica Exchange Molecular Dynamics (FEP/REMD) offers a powerful strategy to improve the convergence of free energy computations. In particular, it has been shown previously that a FEP/REMD scheme allowing random moves within an extended replica ensemble of thermodynamic coupling parameters '{lambda}' can improve the statistical convergence in calculations of absolute binding free energy of ligands to proteins [J. Chem. Theory Comput. 2009, 5, 2583]. In the present study, FEP/REMD is extended and combined with an accelerated MD simulations method based on Hamiltonian replica-exchange MD (H-REMD) to overcome the additional problems arising from the existence of kinetically trapped conformations within the protein receptor. In the combined strategy, each system with a given thermodynamic coupling factor {lambda} in the extended ensemble is further coupled with a set of replicas evolving on a biased energy surface with boosting potentials used to accelerate the interconversion among different rotameric states of the side chains in the neighborhood of the binding site. Exchanges are allowed to occur alternatively along the axes corresponding to the thermodynamic coupling parameter {lambda} and the boosting potential, in an extended dual array of coupled {lambda}- and H-REMD simulations. The method is implemented on the basis of new extensions to the REPDSTR module of the biomolecular simulation program CHARMM. As an illustrative example, the absolute binding free energy of p-xylene to the nonpolar cavity of the L99A mutant of the T4 lysozyme was calculated. The tests demonstrate that the dual {lambda}-REMD and H-REMD simulation scheme greatly accelerates the configurational sampling of the rotameric states of the side chains around the binding pocket, thereby improving the convergence of the FEP computations.
- Research Organization:
- Argonne National Laboratory (ANL)
- Sponsoring Organization:
- NSF
- DOE Contract Number:
- AC02-06CH11357
- OSTI ID:
- 1018895
- Report Number(s):
- ANL/BIO/JA-67686
- Journal Information:
- J. Chem. Theory Comput., Journal Name: J. Chem. Theory Comput. Journal Issue: 9 ; Sep. 2010 Vol. 6; ISSN 1549-9618
- Country of Publication:
- United States
- Language:
- ENGLISH
Similar Records
Computation of absolute hydration and binding free energy with free energy perturbation distributed replica-exchange molecular dynamics.
Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis
Journal Article
·
Sat Aug 01 00:00:00 EDT 2009
· Journal of Chemical Theory and Computation
·
OSTI ID:1034250
Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis
Journal Article
·
Mon Sep 24 20:00:00 EDT 2018
· Journal of Physical Chemistry. B, Condensed Matter, Materials, Surfaces, Interfaces and Biophysical Chemistry
·
OSTI ID:1488538