Path-Integral Method for Predicting Relative Binding Affinities of Protein−Ligand Complexes
- 10 March 2009
- journal article
- research article
- Published by American Chemical Society (ACS) in Journal of the American Chemical Society
- Vol. 131 (12), 4521-4528
- https://doi.org/10.1021/ja807460s
Abstract
We present a novel approach for computing biomolecular interaction binding affinities based on a simple path integral solution of the Fokker−Planck equation. Computing the free energy of protein−ligand interactions can expedite structure-based drug design. Traditionally, the problem is seen through the lens of statistical thermodynamics. The computations can become, however, prohibitively long for the change in the free energy upon binding to be determined accurately. In this work, we present a different approach based on a stochastic kinetic formalism. Inspired by Feynman’s path integral formulation, we extend the theory to classical interacting systems. The ligand is modeled as a Brownian particle subjected to the effective nonbonding interaction potential of the receptor. This allows the calculation of the relative binding affinities of interacting biomolecules in water to be computed as a function of the ligand’s diffusivity and the curvature of the potential surface in the vicinity of the binding minimum. The calculation is thus exceedingly rapid. In test cases, the correlation coefficient between actual and computed free energies is >0.93 for accurate data sets.Keywords
This publication has 52 references indexed in Scilit:
- Optimization of Azoles as Anti-Human Immunodeficiency Virus Agents Guided by Free-Energy CalculationsJournal of the American Chemical Society, 2008
- Assessment of programs for ligand binding affinity predictionJournal of Computational Chemistry, 2008
- Predicting Absolute Ligand Binding Free Energies to a Simple Model SiteJournal of Molecular Biology, 2007
- Structural and thermodynamic studies of simple aldose reductase–inhibitor complexesBioorganic Chemistry, 2006
- Expect the Unexpected or Caveat for Drug Designers: Multiple Structure Determinations Using Aldose Reductase Crystals Treated under Varying Soaking and Co-crystallisation ConditionsJournal of Molecular Biology, 2006
- A novel serine protease inhibition motif involving a multi-centered short hydrogen bonding network at the active site11Edited by D. ReesJournal of Molecular Biology, 2001
- The Protein Data BankNucleic Acids Research, 2000
- Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94Journal of Computational Chemistry, 1996
- Path integral for a class of nonlinear Fokker-Planck processesPhysica, 1967
- Functional integration and the brownian movement in a field of forcePhysica, 1967