Dataset Information

Comparing alchemical and physical pathway methods for computing the absolute binding free energy of charged ligands.

ABSTRACT: Accurately predicting absolute binding free energies of protein-ligand complexes is important as a fundamental problem in both computational biophysics and pharmaceutical discovery. Calculating binding free energies for charged ligands is generally considered to be challenging because of the strong electrostatic interactions between the ligand and its environment in aqueous solution. In this work, we compare the performance of the potential of mean force (PMF) method and the double decoupling method (DDM) for computing absolute binding free energies for charged ligands. We first clarify an unresolved issue concerning the explicit use of the binding site volume to define the complexed state in DDM together with the use of harmonic restraints. We also provide an alternative derivation for the formula for absolute binding free energy using the PMF approach. We use these formulas to compute the binding free energy of charged ligands at an allosteric site of HIV-1 integrase, which has emerged in recent years as a promising target for developing antiviral therapy. As compared with the experimental results, the absolute binding free energies obtained by using the PMF approach show unsigned errors of 1.5-3.4 kcal mol-1, which are somewhat better than the results from DDM (unsigned errors of 1.6-4.3 kcal mol-1) using the same amount of CPU time. According to the DDM decomposition of the binding free energy, the ligand binding appears to be dominated by nonpolar interactions despite the presence of very large and favorable intermolecular ligand-receptor electrostatic interactions, which are almost completely cancelled out by the equally large free energy cost of desolvation of the charged moiety of the ligands in solution. We discuss the relative strengths of computing absolute binding free energies using the alchemical and physical pathway methods.

SUBMITTER: Deng N

PROVIDER: S-EPMC6061996 | biostudies-literature | 2018 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Comparing alchemical and physical pathway methods for computing the absolute binding free energy of charged ligands.

Deng Nanjie N Cui Di D Zhang Bin W BW Xia Junchao J Cruz Jeffrey J Levy Ronald R

Physical chemistry chemical physics : PCCP 20180601 25

Accurately predicting absolute binding free energies of protein-ligand complexes is important as a fundamental problem in both computational biophysics and pharmaceutical discovery. Calculating binding free energies for charged ligands is generally considered to be challenging because of the strong electrostatic interactions between the ligand and its environment in aqueous solution. In this work, we compare the performance of the potential of mean force (PMF) method and the double decoupling me ...[more]

PMID: 29896599

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

Comparing alchemical and physical pathway methods for computing the absolute binding free energy of charged ligands.

Publications

Comparing alchemical and physical pathway methods for computing the absolute binding free energy of charged ligands.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets