Estimates of ligand-binding affinities supported by quantum mechanical methods

In this paper, we review our efforts to use quantum mechanical (QM) methods to improve free-energy estimates of the binding of drug candidates to their receptor proteins. First, we have tested the influence of various implicit solvation models on predictions of the ligand-binding affinity. The accuracy of implicit solvation models strongly depend on the parameterisation, but also on the magnitude of the solvation energy (i.e. their accuracy should be discussed in relative terms). However, if only relative solvation energies within a series of similar drug molecules with the same net charge are considered, nearly all methods tested give a comparable accuracy of 2–5 kJ/mol. Second, we have studied the conformational dependence of QM charges and their influence on ligand-binding affinities. The conformational dependence is significant, but it is to a large extent cancelled by solvation energies. Third, we have estimated the effect and range of electrostatic interactions beyond a point-charge model. The results show that multipoles up to octupoles and anisotropic polarisabilities have a significant influence on energies for residues up to 10–15 °A from the ligand and that different sets of point-charge models may give strongly varying results. However, if only relative energies are considered, the effect is to a large extent cancelled. Fourth, we have tried to develop an accurate QM-based molecular mechanics potential, in which not only the electrostatic terms are improved, but also the dispersion and repulsion. However, even with quite sophisticated expressions, it seems difficult to reduce the average error below 2–3 kJ/mol per interaction (e.g. a hydrogen bond), compared to the full QM treatment. Finally, we have developed a new method, PMISP (polarised multipole interaction with supermolecular pairs), for the calculation of accurate interaction energies. It employs an accurate force field for electrostatics and induction, including multipoles up to octupoles and anisotropic polarisabilities calculated by QM methods on amino-acid fragments of the protein in each conformation observed in snapshots from a molecular dynamics simulation, whereas short-range interactions are estimated by high-level QM calculations for all pairs of the ligand with near-by residues. We show that this approach allows us to go far beyond the current accuracy of molecular mechanics methods, down to an error of 5–10 kJ/mol for a full protein-ligand complex. It can be combined with estimates of solvation, entropy, and dynamic effects to give estimates of binding affinities. However, several problems remain to be solved before any significant improvement in the accuracy can be seen.

[1]  William H. Fink,et al.  Frozen fragment reduced variational space analysis of hydrogen bonding interactions. Application to the water dimer , 1987 .

[2]  Roland Lindh,et al.  Accuracy of distributed multipoles and polarizabilities: Comparison between the LoProp and MpProp models , 2007, J. Comput. Chem..

[3]  Samuel Genheden,et al.  A comparison of different initialization protocols to obtain statistically independent molecular dynamics simulations , 2011, J. Comput. Chem..

[4]  Fumio Hirata,et al.  Theory of Molecular Liquids , 2004 .

[5]  Ulf Ryde,et al.  Accurate metal-site structures in proteins obtained by combining experimental data and quantum chemistry. , 2007, Dalton transactions.

[6]  M. Alderton,et al.  Distributed multipole analysis , 2006 .

[7]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[8]  David L Mobley,et al.  Predicting small-molecule solvation free energies: an informal blind test for computational chemistry. , 2008, Journal of medicinal chemistry.

[9]  Arieh Warshel,et al.  Polarizable Force Fields:  History, Test Cases, and Prospects. , 2007, Journal of chemical theory and computation.

[10]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[11]  On the coupling of intermolecular polarization and repulsion through pseudo-potentials , 2009 .

[12]  Arieh Warshel,et al.  Calculations of chemical processes in solutions , 1979 .

[13]  T. Halgren,et al.  Polarizable force fields. , 2001, Current opinion in structural biology.

[14]  Nohad Gresh,et al.  Binding of 5‐phospho‐D‐arabinonohydroxamate and 5‐phospho‐D‐arabinonate inhibitors to zinc phosphomannose isomerase from Candida albicans studied by polarizable molecular mechanics and quantum mechanics , 2007, J. Comput. Chem..

[15]  Jonathan W. Essex,et al.  Errors in free-energy perturbation calculations due to neglecting the conformational variation of atomic charges , 1992 .

[16]  Frank Neese,et al.  A critical evaluation of DFT, including time-dependent DFT, applied to bioinorganic chemistry , 2006, JBIC Journal of Biological Inorganic Chemistry.

[17]  Pavel Hobza,et al.  Assessment of the MP2 method, along with several basis sets, for the computation of interaction energies of biologically relevant hydrogen bonded and dispersion bound complexes. , 2007, The journal of physical chemistry. A.

[18]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[19]  Kazuo Kitaura,et al.  Molecular recognition mechanism of FK506 binding protein: An all‐electron fragment molecular orbital study , 2007, Proteins.

[20]  S. Durell,et al.  Specificity of acyl transfer from 2-mercaptobenzamide thioesters to the HIV-1 nucleocapsid protein. , 2007, Journal of the American Chemical Society.

[21]  Fahmi Himo,et al.  Quantum chemical modeling of enzyme active sites and reaction mechanisms , 2006 .

[22]  K. Sharp,et al.  Electrostatic interactions in macromolecules: theory and applications. , 1990, Annual review of biophysics and biophysical chemistry.

[23]  Alessandro Laio,et al.  A Hamiltonian electrostatic coupling scheme for hybrid Car-Parrinello molecular dynamics simulations , 2002 .

[24]  John Z H Zhang,et al.  Quantum mechanical map for protein-ligand binding with application to beta-trypsin/benzamidine complex. , 2004, The Journal of chemical physics.

[25]  W. J. Orville-Thomas.  Atoms in Molecules — a Quantum Theory , 1996 .

[26]  Nohad Gresh,et al.  Anisotropic, Polarizable Molecular Mechanics Studies of Inter- and Intramolecular Interactions and Ligand-Macromolecule Complexes. A Bottom-Up Strategy. , 2007, Journal of chemical theory and computation.

[27]  Yun Xiang,et al.  Quantum study of mutational effect in binding of efavirenz to HIV‐1 RT , 2005, Proteins.

[28]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[29]  Jan H. Jensen,et al.  Continuum solvation of large molecules described by QM/MM: a semi-iterative implementation of the PCM/EFP interface , 2003 .

[30]  Bruce L. Bush,et al.  Restrained point‐charge models for disaccharides , 2002, J. Comput. Chem..

[31]  Donald G Truhlar,et al.  Electrostatically Embedded Many-Body Correlation Energy, with Applications to the Calculation of Accurate Second-Order Møller-Plesset Perturbation Theory Energies for Large Water Clusters. , 2007, Journal of chemical theory and computation.

[32]  Samuel Genheden,et al.  How to obtain statistically converged MM/GBSA results , 2009, J. Comput. Chem..

[33]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[34]  Kaori Fukuzawa,et al.  Molecular interactions between estrogen receptor and its ligand studied by the ab initio fragment molecular orbital method. , 2006, The journal of physical chemistry. B.

[35]  Bing Wang,et al.  The role of quantum mechanics in structure-based drug design. , 2007, Drug discovery today.

[36]  Jirí Cerný,et al.  Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs. , 2006, Physical chemistry chemical physics : PCCP.

[37]  Mark S. Gordon,et al.  Evaluation of Charge Penetration Between Distributed Multipolar Expansions , 2000 .

[38]  Giulio Rastelli,et al.  Fast and accurate predictions of binding free energies using MM‐PBSA and MM‐GBSA , 2009, J. Comput. Chem..

[39]  Richard H. Henchman,et al.  Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy. , 2004, Biophysical journal.

[40]  Donald E. Williams,et al.  Alanyl dipeptide potential‐derived net atomic charges and bond dipoles, and their variation with molecular conformation , 1990 .

[41]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[42]  David Chandler,et al.  Optimized Cluster Expansions for Classical Fluids. II. Theory of Molecular Liquids , 1972 .

[43]  Steven M. Bachrach,et al.  Population Analysis and Electron Densities from Quantum Mechanics , 2007 .

[44]  Paul S. Bagus,et al.  A new analysis of charge transfer and polarization for ligand–metal bonding: Model studies of Al4CO and Al4NH3 , 1984 .

[45]  M. Gilson,et al.  Calculation of protein-ligand binding affinities. , 2007, Annual review of biophysics and biomolecular structure.

[46]  Yixiang Cao,et al.  A Polarizable Force Field and Continuum Solvation Methodology for Modeling of Protein-Ligand Interactions. , 2005, Journal of chemical theory and computation.

[47]  U. Ryde,et al.  Comparison of overlap-based models for approximating the exchange-repulsion energy. , 2006, The Journal of chemical physics.

[48]  Pengyu Y. Ren,et al.  Consistent treatment of inter‐ and intramolecular polarization in molecular mechanics calculations , 2002, J. Comput. Chem..

[49]  K. Merz,et al.  Large-scale validation of a quantum mechanics based scoring function: predicting the binding affinity and the binding mode of a diverse set of protein-ligand complexes. , 2005, Journal of medicinal chemistry.

[50]  Pär Söderhjelm,et al.  Accuracy of typical approximations in classical models of intermolecular polarization. , 2008, The Journal of chemical physics.

[51]  J. H. Zhang,et al.  Full quantum mechanical study of binding of HIV‐1 protease drugs , 2005 .

[52]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[53]  Ye Mei,et al.  Quantum and Molecular Dynamics Study for Binding of Macrocyclic Inhibitors to Human α-Thrombin , 2007 .

[54]  Ulf Ryde,et al.  Comparison of methods for deriving atomic charges from the electrostatic potential and moments , 1998 .

[55]  Jacob Kongsted,et al.  An improved method to predict the entropy term with the MM/PBSA approach , 2009, J. Comput. Aided Mol. Des..

[56]  Many-body force field models based solely on pairwise Coulomb screening do not simultaneously reproduce correct gas-phase and condensed-phase polarizability limits. , 2004, The Journal of chemical physics.

[57]  T. Darden,et al.  Intermolecular electrostatic energies using density fitting. , 2005, The Journal of chemical physics.

[58]  Edward I. Solomon,et al.  Computational inorganic and bioinorganic chemistry , 2009 .

[59]  Gregory D. Hawkins,et al.  Parametrized Models of Aqueous Free Energies of Solvation Based on Pairwise Descreening of Solute Atomic Charges from a Dielectric Medium , 1996 .

[60]  B. Kuhn,et al.  Validation and use of the MM-PBSA approach for drug discovery. , 2005, Journal of medicinal chemistry.

[61]  Fumio Hirata,et al.  Potentials of mean force of simple ions in ambient aqueous solution. I. Three-dimensional reference interaction site model approach , 2000 .

[62]  Wei Zhang,et al.  A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations , 2003, J. Comput. Chem..

[63]  J. Tomasi,et al.  Ab initio study of solvated molecules: A new implementation of the polarizable continuum model , 1996 .

[64]  Ryan P. A. Bettens,et al.  On the accurate reproduction of ab initio interaction energies between an enzyme and substrate , 2007 .

[65]  Kazuo Kitaura,et al.  A new energy decomposition scheme for molecular interactions within the Hartree‐Fock approximation , 1976 .

[66]  Xiao He,et al.  Quantum computational analysis for drug resistance of HIV‐1 reverse transcriptase to nevirapine through point mutations , 2005, Proteins.

[67]  U. Ryde.  Quantum Mechanical/Molecular Mechanical (QM/MM) Methods and Applictions in Bioinorganic Chemistry , 2009 .

[68]  Holger Gohlke,et al.  Converging free energy estimates: MM‐PB(GB)SA studies on the protein–protein complex Ras–Raf , 2004, J. Comput. Chem..

[69]  Robert Bell Hermann,et al.  Theory of Hydrophobic Bonding , 1974 .

[70]  Johan Åqvist,et al.  Ligand binding affinity prediction by linear interaction energy methods , 1998, J. Comput. Aided Mol. Des..

[71]  John Z H Zhang,et al.  New Advance in Computational Chemistry:  Full Quantum Mechanical ab Initio Computation of Streptavidin-Biotin Interaction Energy. , 2003, The journal of physical chemistry. B.

[72]  Nohad Gresh,et al.  Improved Formulas for the Calculation of the Electrostatic Contribution to the Intermolecular Interaction Energy from Multipolar Expansion of the Electronic Distribution. , 2003, The journal of physical chemistry. A.

[73]  R. S. Mulliken.  Electronic Population Analysis on LCAO–MO Molecular Wave Functions. I , 1955 .

[74]  Robert B. Hermann,et al.  Theory of hydrophobic bonding. II. Correlation of hydrocarbon solubility in water with solvent cavity surface area , 1972 .

[75]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[76]  A. Warshel,et al.  Calculations of Hydration Entropies of Hydrophobic, Polar, and Ionic Solutes in the Framework of the Langevin Dipoles Solvation Model , 1999 .

[77]  O Engkvist,et al.  Accurate Intermolecular Potentials Obtained from Molecular Wave Functions: Bridging the Gap between Quantum Chemistry and Molecular Simulations. , 2000, Chemical reviews.

[78]  Samuel Genheden,et al.  Transferability of conformational dependent charges from protein simulations , 2012 .

[79]  A. Stone,et al.  Distributed dispersion: A new approach , 2003 .

[80]  M. Probst,et al.  On the performance of molecular polarization methods. II. Water and carbon tetrachloride close to a cation. , 2005, The Journal of chemical physics.

[81]  P. Kollman,et al.  Atomic charges derived from semiempirical methods , 1990 .

[82]  Pär Söderhjelm,et al.  Conformational dependence of charges in protein simulations , 2009, J. Comput. Chem..

[83]  M. Alderton,et al.  Distributed multipole analysis Methods and applications , 1985 .

[84]  Francesco Aquilante,et al.  Calculation of protein-ligand interaction energies by a fragmentation approach combining high-level quantum chemistry with classical many-body effects. , 2009, The journal of physical chemistry. B.

[85]  K. Kitaura,et al.  Fragment molecular orbital method: an approximate computational method for large molecules , 1999 .

[86]  Jacopo Tomasi,et al.  Quantum Mechanical Continuum Solvation Models , 2005 .

[87]  K. Merz,et al.  A quantum mechanics-based scoring function: study of zinc ion-mediated ligand binding. , 2004, Journal of the American Chemical Society.

[88]  Tomasz Borowski,et al.  Modeling enzymatic reactions involving transition metals. , 2006, Accounts of chemical research.

[89]  U. Ryde,et al.  Ligand affinities predicted with the MM/PBSA method: dependence on the simulation method and the force field. , 2006, Journal of Medicinal Chemistry.

[90]  G. Klebe,et al.  Approaches to the Description and Prediction of the Binding Affinity of Small-Molecule Ligands to Macromolecular Receptors , 2002 .

[91]  Mark S. Gordon,et al.  An effective fragment method for modeling solvent effects in quantum mechanical calculations , 1996 .

[92]  J Andrew McCammon,et al.  Generalized Born model with a simple, robust molecular volume correction. , 2007, Journal of chemical theory and computation.

[93]  Roland Lindh,et al.  Local properties of quantum chemical systems: the LoProp approach. , 2004, The Journal of chemical physics.

[94]  Pär Söderhjelm,et al.  Protein Influence on Electronic Spectra Modeled by Multipoles and Polarizabilities. , 2009, Journal of chemical theory and computation.

[95]  Jacopo Tomasi,et al.  A new definition of cavities for the computation of solvation free energies by the polarizable continuum model , 1997 .

[96]  Claude Millot,et al.  Fast and accurate determination of induction energies: reduction of topologically distributed polarizability models , 2001 .

[97]  Pär Söderhjelm,et al.  How accurate can a force field become? A polarizable multipole model combined with fragment-wise quantum-mechanical calculations. , 2009, The journal of physical chemistry. A.

[98]  Nohad Gresh,et al.  Representation of Zn(II) complexes in polarizable molecular mechanics. Further refinements of the electrostatic and short‐range contributions. Comparisons with parallel ab initio computations , 2005, J. Comput. Chem..

[99]  Ray Luo,et al.  Implicit nonpolar solvent models. , 2007, The journal of physical chemistry. B.

[100]  John Z. H. Zhang,et al.  Molecular fractionation with conjugate caps for full quantum mechanical calculation of protein-molecule interaction energy , 2003 .

[101]  Terry R. Stouch,et al.  Conformational dependence of electrostatic potential derived charges of a lipid headgroup: Glycerylphosphorylcholine , 1992 .

[102]  Andriy Kovalenko,et al.  An MM/3D-RISM approach for ligand binding affinities. , 2010, The journal of physical chemistry. B.