Protein-ligand interaction energies with dispersion corrected density functional theory and high-level wave function based methods.

With dispersion-corrected density functional theory (DFT-D3) intermolecular interaction energies for a diverse set of noncovalently bound protein-ligand complexes from the Protein Data Bank are calculated. The focus is on major contacts occurring between the drug molecule and the binding site. Generalized gradient approximation (GGA), meta-GGA, and hybrid functionals are used. DFT-D3 interaction energies are benchmarked against the best available wave function based results that are provided by the estimated complete basis set (CBS) limit of the local pair natural orbital coupled-electron pair approximation (LPNO-CEPA/1) and compared to MP2 and semiempirical data. The size of the complexes and their interaction energies (ΔE(PL)) varies between 50 and 300 atoms and from -1 to -65 kcal/mol, respectively. Basis set effects are considered by applying extended sets of triple- to quadruple-ζ quality. Computed total ΔE(PL) values show a good correlation with the dispersion contribution despite the fact that the protein-ligand complexes contain many hydrogen bonds. It is concluded that an adequate, for example, asymptotically correct, treatment of dispersion interactions is necessary for the realistic modeling of protein-ligand binding. Inclusion of the dispersion correction drastically reduces the dependence of the computed interaction energies on the density functional compared to uncorrected DFT results. DFT-D3 methods provide results that are consistent with LPNO-CEPA/1 and MP2, the differences of about 1-2 kcal/mol on average (<5% of ΔE(PL)) being on the order of their accuracy, while dispersion-corrected semiempirical AM1 and PM3 approaches show a deviating behavior. The DFT-D3 results are found to depend insignificantly on the choice of the short-range damping model. We propose to use DFT-D3 as an essential ingredient in a QM/MM approach for advanced virtual screening approaches of protein-ligand interactions to be combined with similarly "first-principle" accounts for the estimation of solvation and entropic effects.

[1]  Edward G Hohenstein,et al.  Basis set consistent revision of the S22 test set of noncovalent interaction energies. , 2010, The Journal of chemical physics.

[2]  Kichul Park,et al.  Quantum mechanical scoring for protein docking. , 2009, The Journal of chemical physics.

[3]  S. Grimme.  Density functional theory with London dispersion corrections , 2011 .

[4]  Jacob Kongsted,et al.  Estimates of ligand-binding affinities supported by quantum mechanical methods , 2010, Interdisciplinary Sciences: Computational Life Sciences.

[5]  Parr,et al.  Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. , 1988, Physical review. B, Condensed matter.

[6]  S. Grimme,et al.  A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. , 2010, The Journal of chemical physics.

[7]  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.

[8]  Joachim Sauer,et al.  Application of semiempirical long‐range dispersion corrections to periodic systems in density functional theory , 2008, J. Comput. Chem..

[9]  Claudio N. Cavasotto,et al.  The Challenge of Considering Receptor Flexibility in Ligand Docking and Virtual Screening , 2005 .

[10]  K. Baldridge,et al.  Steric isotope effects gauged by the bowl-inversion barrier in selectively deuterated pentaarylcorannulenes. , 2008, Journal of the American Chemical Society.

[11]  S. Grimme,et al.  IR/UV spectra and quantum chemical calculations of Trp-Ser: stacking interactions between backbone and indole side-chain. , 2008, Physical chemistry chemical physics : PCCP.

[12]  Xiaoqin Zou,et al.  Scoring functions and their evaluation methods for protein-ligand docking: recent advances and future directions. , 2010, Physical chemistry chemical physics : PCCP.

[13]  Stefan Grimme,et al.  Effect of the damping function in dispersion corrected density functional theory , 2011, J. Comput. Chem..

[14]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials , 1995 .

[15]  G. Klebe,et al.  Approaches to the description and prediction of the binding affinity of small-molecule ligands to macromolecular receptors. , 2002, Angewandte Chemie.

[16]  Kenneth M Merz,et al.  Semiempirical Comparative Binding Energy Analysis (SE-COMBINE) of a Series of Trypsin Inhibitors. , 2006, Journal of chemical theory and computation.

[17]  M. Otyepka,et al.  Interaction energies for the purine inhibitor roscovitine with cyclin-dependent kinase 2: correlated ab initio quantum-chemical, DFT and empirical calculations. , 2006, Chemistry.

[18]  S. Grimme,et al.  Density functional theory with dispersion corrections for supramolecular structures, aggregates, and complexes of (bio)organic molecules. , 2007, Organic & biomolecular chemistry.

[19]  Kenneth M Merz,et al.  A Mixed QM/MM Scoring Function to Predict Protein-Ligand Binding Affinity. , 2010, Journal of chemical theory and computation.

[20]  F. Weigend,et al.  Gaussian basis sets of quadruple zeta valence quality for atoms H–Kr , 2003 .

[21]  Kwang Soo Kim,et al.  Molecular Clusters of pi-Systems: Theoretical Studies of Structures, Spectra, and Origin of Interaction Energies. , 2000, Chemical reviews.

[22]  Jindřich Fanfrlík,et al.  Semiempirical Quantum Chemical PM6 Method Augmented by Dispersion and H-Bonding Correction Terms Reliably Describes Various Types of Noncovalent Complexes. , 2009, Journal of chemical theory and computation.

[23]  Sándor Suhai,et al.  Self-consistent-charge density-functional tight-binding method for simulations of complex materials properties , 1998 .

[24]  M. Preuss,et al.  Organic molecule adsorption on solid surfaces: chemical bonding, mutual polarisation and dispersion interaction , 2006 .

[25]  W. Klopper,et al.  Ab Initio Modeling of Methanol Interaction with Single-Walled Carbon Nanotubes , 2007 .

[26]  M. Vincent,et al.  Can the DFT-D method describe the full range of noncovalent interactions found in large biomolecules? , 2007, Physical Chemistry, Chemical Physics - PCCP.

[27]  S. Grimme,et al.  Density functional theory including dispersion corrections for intermolecular interactions in a large benchmark set of biologically relevant molecules. , 2006, Physical chemistry chemical physics : PCCP.

[28]  Stefan Grimme,et al.  Accurate description of van der Waals complexes by density functional theory including empirical corrections , 2004, J. Comput. Chem..

[29]  Kevin E. Riley,et al.  A DFT-D investigation of the mechanisms for activation of the wild-type and S810L mutated mineralocorticoid receptor by steroid hormones. , 2008, Journal of Physical Chemistry B.

[30]  Stefan Grimme,et al.  Semiempirical GGA‐type density functional constructed with a long‐range dispersion correction , 2006, J. Comput. Chem..

[31]  Wang,et al.  Accurate and simple analytic representation of the electron-gas correlation energy. , 1992, Physical review. B, Condensed matter.

[32]  Pietro Cozzini,et al.  Getting it right: modeling of pH, solvent and "nearly" everything else in virtual screening of biological targets. , 2004, Journal of molecular graphics & modelling.

[33]  B. Sumpter,et al.  Density-functional approaches to noncovalent interactions: a comparison of dispersion corrections (DFT-D), exchange-hole dipole moment (XDM) theory, and specialized functionals. , 2011, The Journal of chemical physics.

[34]  J. Perdew,et al.  Density-functional approximation for the correlation energy of the inhomogeneous electron gas. , 1986, Physical review. B, Condensed matter.

[35]  A Srinivas Reddy,et al.  Virtual screening in drug discovery -- a computational perspective. , 2007, Current protein & peptide science.

[36]  C David Sherrill,et al.  Formal Estimation of Errors in Computed Absolute Interaction Energies of Protein-ligand Complexes. , 2011, Journal of chemical theory and computation.

[37]  H. Senn,et al.  QM/MM Methods for Biological Systems , 2006 .

[38]  H. Valdés,et al.  Ab Initio Study on the (OCS)2·CO2 van der Waals Trimers , 2002 .

[39]  Vladislav Vasilyev,et al.  Application of semiempirical quantum chemical methods as a scoring function in docking , 2004 .

[40]  Trygve Helgaker,et al.  Basis set convergence of the interaction energy of hydrogen-bonded complexes , 1999 .

[41]  Jirí Cerný,et al.  Density functional theory augmented with an empirical dispersion term. Interaction energies and geometries of 80 noncovalent complexes compared with ab initio quantum mechanics calculations , 2007, J. Comput. Chem..

[42]  Jeremy C. Smith,et al.  Protein/ligand binding free energies calculated with quantum mechanics/molecular mechanics. , 2005, The journal of physical chemistry. B.

[43]  Lifeng Chi,et al.  Simple and complex lattices of N-alkyl fatty acid amides on a highly oriented pyrolytic graphite surface. , 2005, Langmuir : the ACS journal of surfaces and colloids.

[44]  Kevin E. Riley,et al.  Pairwise decomposition of residue interaction energies using semiempirical quantum mechanical methods in studies of protein-ligand interaction. , 2005, Journal of the American Chemical Society.

[45]  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.

[46]  Frank Neese,et al.  Accurate theoretical chemistry with coupled pair models. , 2009, Accounts of chemical research.

[47]  Pavel Hobza,et al.  A Transferable H-Bonding Correction for Semiempirical Quantum-Chemical Methods. , 2010, Journal of chemical theory and computation.

[48]  Dimitrios G Liakos,et al.  Efficient and accurate approximations to the local coupled cluster singles doubles method using a truncated pair natural orbital basis. , 2009, The Journal of chemical physics.

[49]  J. Brynda,et al.  Interpretation of protein/ligand crystal structure using QM/MM calculations: case of HIV-1 protease/metallacarborane complex. , 2008, The journal of physical chemistry. B.

[50]  F. Neese,et al.  Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method. , 2009, The Journal of chemical physics.

[51]  Vincenzo Barone,et al.  Role and effective treatment of dispersive forces in materials: Polyethylene and graphite crystals as test cases , 2009, J. Comput. Chem..

[52]  G. Scuseria,et al.  Climbing the density functional ladder: nonempirical meta-generalized gradient approximation designed for molecules and solids. , 2003, Physical review letters.

[53]  A. Becke,et al.  Density-functional exchange-energy approximation with correct asymptotic behavior. , 1988, Physical review. A, General physics.

[54]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[55]  Stefan Grimme,et al.  Noncovalent Interactions between Graphene Sheets and in Multishell (Hyper)Fullerenes , 2007 .

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

[57]  Stefan Grimme,et al.  Van der Waals interactions in aromatic systems: structure and energetics of dimers and trimers of pyridine. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.

[58]  L. Cavallo,et al.  Parametrization of an empirical correction term to density functional theory for an accurate description of pi-stacking interactions in nucleic acids. , 2007, The journal of physical chemistry. B.

[59]  Pavel Hobza,et al.  A reliable docking/scoring scheme based on the semiempirical quantum mechanical PM6-DH2 method accurately covering dispersion and H-bonding: HIV-1 protease with 22 ligands. , 2010, The journal of physical chemistry. B.

[60]  Harald Lanig,et al.  Quantum mechanical/molecular mechanical (QM/MM) docking: an evaluation for known test systems , 2003 .

[61]  S. Grimme,et al.  Is spin-component scaled second-order Møller-Plesset perturbation theory an appropriate method for the study of noncovalent interactions in molecules? , 2007, The journal of physical chemistry. A.

[62]  S. Grimme,et al.  Structure and binding energies of the porphine dimer , 2007 .

[63]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[64]  Thomas Frauenheim,et al.  Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment , 2001 .

[65]  Renxiao Wang,et al.  Comparative evaluation of 11 scoring functions for molecular docking. , 2003, Journal of medicinal chemistry.

[66]  Kenneth M Merz,et al.  Limits of Free Energy Computation for Protein-Ligand Interactions. , 2010, Journal of chemical theory and computation.

[67]  Kazumasa Honda,et al.  High-level ab initio computations of structures and interaction energies of naphthalene dimers: origin of attraction and its directionality. , 2004, The Journal of chemical physics.

[68]  Xin Wen,et al.  BindingDB: a web-accessible database of experimentally determined protein–ligand binding affinities , 2006, Nucleic Acids Res..

[69]  Renxiao Wang,et al.  The PDBbind database: methodologies and updates. , 2005, Journal of medicinal chemistry.

[70]  A. Becke.  Density-functional thermochemistry. III. The role of exact exchange , 1993 .

[71]  M. Frisch,et al.  Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields , 1994 .

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

[73]  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.

[74]  Matthew Paul Gleeson,et al.  QM/MM Calculations in Drug Discovery: A Useful Method for Studying Binding Phenomena? , 2009, J. Chem. Inf. Model..

[75]  Cláudio Nahum Alves,et al.  A quantum mechanics/molecular mechanics study of the protein-ligand interaction of two potent inhibitors of human O-GlcNAcase: PUGNAc and NAG-thiazoline. , 2008, The journal of physical chemistry. B.

[76]  Jacob Kongsted,et al.  Ligand Affinities Estimated by Quantum Chemical Calculations. , 2010, Journal of chemical theory and computation.

[77]  Takeshi Ishikawa,et al.  Theoretical study of the prion protein based on the fragment molecular orbital method , 2009, J. Comput. Chem..

[78]  Supot Hannongbua,et al.  Insights into saquinavir resistance in the G48V HIV-1 protease: quantum calculations and molecular dynamic simulations. , 2005, Biophysical journal.

[79]  J. Perdew,et al.  Erratum: Density-functional approximation for the correlation energy of the inhomogeneous electron gas , 1986, Physical review. B, Condensed matter.

[80]  J. H. Zhang,et al.  Quantum study of HIV-1 protease-bridge water interaction. , 2007, The Journal of chemical physics.

[81]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[82]  Stefan Grimme,et al.  A Theoretical Investigation of the Geometries and Binding Energies of Molecular Tweezer and Clip Host-Guest Systems. , 2005, Journal of chemical theory and computation.

[83]  Dimitrios G Liakos,et al.  Weak Molecular Interactions Studied with Parallel Implementations of the Local Pair Natural Orbital Coupled Pair and Coupled Cluster Methods. , 2011, Journal of chemical theory and computation.

[84]  A. Schäfer,et al.  Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr , 1994 .

[85]  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.

[86]  Stefan Grimme,et al.  Van der Waals complexes of polar aromatic molecules: unexpected structures for dimers of azulene. , 2005, Journal of the American Chemical Society.

[87]  Art E. Cho,et al.  Extension of QM/MM docking and its applications to metalloproteins , 2009, J. Comput. Chem..

[88]  Izhar Wallach,et al.  The protein-small-molecule database, a non-redundant structural resource for the analysis of protein-ligand binding , 2009, Bioinform..

[89]  Holger Patzelt,et al.  RI-MP2: optimized auxiliary basis sets and demonstration of efficiency , 1998 .

[90]  Akash Khandelwal,et al.  A combination of docking, QM/MM methods, and MD simulation for binding affinity estimation of metalloprotein ligands. , 2005, Journal of medicinal chemistry.

[91]  Marco Häser,et al.  Auxiliary basis sets to approximate Coulomb potentials (Chem. Phys. Letters 240 (1995) 283-290) , 1995 .

[92]  Gerhard Klebe,et al.  AffinDB: a freely accessible database of affinities for protein–ligand complexes from the PDB , 2005, Nucleic Acids Res..

[93]  N. Burton,et al.  A QM/MM study of fluoroaromatic interactions at the binding site of carbonic anhydrase II, using a DFT method corrected for dispersive interactions. , 2008, Physical chemistry chemical physics : PCCP.

[94]  G. Schaftenaar,et al.  Molden: a pre- and post-processing program for molecular and electronic structures* , 2000, J. Comput. Aided Mol. Des..

[95]  Hans W. Horn,et al.  Fully optimized contracted Gaussian basis sets for atoms Li to Kr , 1992 .

[96]  S. Grimme,et al.  Structures and interaction energies of stacked graphene-nucleobase complexes. , 2008, Physical chemistry chemical physics : PCCP.

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

[98]  Kwang S. Kim,et al.  A theoretical investigation of the nature of the π-H interaction in ethene–H2O, benzene–H2O, and benzene–(H2O)2 , 1999 .

[99]  Richard D. Smith,et al.  Binding MOAD, a high-quality protein–ligand database , 2007, Nucleic Acids Res..