Dataset Information

Effects of Biomolecular Flexibility on Alchemical Calculations of Absolute Binding Free Energies.

ABSTRACT: The independent trajectory thermodynamic integration (IT-TI) approach (Lawrenz et. al J. Chem. Theory. Comput. 2009, 5:1106-1116(1)) for free energy calculations with distributed computing is employed to study two distinct cases of protein-ligand binding: first, the influenza surface protein N1 neuraminidase bound to the inhibitor oseltamivir, and second, the M. tuberculosis enzyme RmlC complexed with the molecule CID 77074. For both systems, finite molecular dynamics (MD) sampling and varied molecular flexibility give rise to IT-TI free energy distributions that are remarkably centered on the target experimental values, with a spread directly related to protein, ligand, and solvent dynamics. Using over 2 ?s of total MD simulation, alternative protocols for the practical, general implementation of IT-TI are investigated, including the optimal use of distributed computing, the total number of alchemical intermediates, and the procedure to perturb electrostatics and van der Waals interactions. A protocol that maximizes predictive power and computational efficiency is proposed. IT-TI outperforms traditional TI predictions and allows a straightforward evaluation of the reliability of free energy estimates. Our study has broad implications for the use of distributed computing in free energy calculations of macromolecular systems.

SUBMITTER: Lawrenz M

PROVIDER: S-EPMC3146769 | biostudies-literature | 2011 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Effects of Biomolecular Flexibility on Alchemical Calculations of Absolute Binding Free Energies.

Lawrenz Morgan M Baron Riccardo R Wang Yi Y McCammon J Andrew JA

Journal of chemical theory and computation 20110601 7

The independent trajectory thermodynamic integration (IT-TI) approach (Lawrenz et. al J. Chem. Theory. Comput. 2009, 5:1106-1116(1)) for free energy calculations with distributed computing is employed to study two distinct cases of protein-ligand binding: first, the influenza surface protein N1 neuraminidase bound to the inhibitor oseltamivir, and second, the M. tuberculosis enzyme RmlC complexed with the molecule CID 77074. For both systems, finite molecular dynamics (MD) sampling and varied mo ...[more]

PMID: 21811708

Similar Datasets

Project description:We present a new approach to more accurately and efficiently compute the absolute binding free energy for receptor-ligand complexes. Currently, the double decoupling method (DDM) and the potential of mean force method (PMF) are widely used to compute the absolute binding free energy of biomolecular complexes. DDM relies on alchemically decoupling the ligand from its environments, which can be computationally challenging for large ligands and charged ligands because of the large magnitude of the decoupling free energies involved. In contrast, the PMF method uses a physical pathway to directly transfer the ligand from solution to the receptor binding pocket and thus avoids some of the aforementioned problems in DDM. However, the PMF method has its own drawbacks: because of its reliance on a ligand binding/unbinding pathway that is free of steric obstructions from the receptor atoms, the method has difficulty treating ligands with buried atoms. To overcome the limitation in the standard PMF approach and enable buried ligands to be treated, here we develop a new method called AlchemPMF in which steric obstructions along the physical pathway for binding are alchemically removed. We have tested the new approach on two important drug targets involving charged ligands. One is HIV-1 integrase bound to an allosteric inhibitor; the other is the human telomeric DNA G-quadruplex in complex with a natural product protoberberine buried in the binding pocket. For both systems, the new approach leads to more reliable estimates of absolute binding free energies with smaller error bars and closer agreements with experiments compared with those obtained from the existing methods, demonstrating the effectiveness of the new method in overcoming the hysteresis often encountered in PMF binding free energy calculations of such systems. The new approach could also be used to improve the sampling of water equilibration and resolvation of the binding pocket as the ligand is extracted.

Dataset Information

Effects of Biomolecular Flexibility on Alchemical Calculations of Absolute Binding Free Energies.

Publications

Effects of Biomolecular Flexibility on Alchemical Calculations of Absolute Binding Free Energies.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets