Virtual screening using molecular simulations

Effective virtual screening relies on our ability to make accurate prediction of protein‐ligand binding, which remains a great challenge. In this work, utilizing the molecular‐mechanics Poisson‐Boltzmann (or Generalized Born) surface area approach, we have evaluated the binding affinity of a set of 156 ligands to seven families of proteins, trypsin β, thrombin α, cyclin‐dependent kinase (CDK), cAMP‐dependent kinase (PKA), urokinase‐type plasminogen activator, β‐glucosidase A, and coagulation factor Xa. The effect of protein dielectric constant in the implicit‐solvent model on the binding free energy calculation is shown to be important. The statistical correlations between the binding energy calculated from the implicit‐solvent approach and experimental free energy are in the range of 0.56–0.79 across all the families. This performance is better than that of typical docking programs especially given that the latter is directly trained using known binding data whereas the molecular mechanics is based on general physical parameters. Estimation of entropic contribution remains the barrier to accurate free energy calculation. We show that the traditional rigid rotor harmonic oscillator approximation is unable to improve the binding free energy prediction. Inclusion of conformational restriction seems to be promising but requires further investigation. On the other hand, our preliminary study suggests that implicit‐solvent based alchemical perturbation, which offers explicit sampling of configuration entropy, can be a viable approach to significantly improve the prediction of binding free energy. Overall, the molecular mechanics approach has the potential for medium to high‐throughput computational drug discovery. Proteins 2011; © 2011 Wiley‐Liss, Inc.

[1]  A. Li,et al.  Screening for human ADME/Tox drug properties in drug discovery. , 2001, Drug discovery today.

[2]  M Rarey,et al.  Detailed analysis of scoring functions for virtual screening. , 2001, Journal of medicinal chemistry.

[3]  Margaret E. Johnson,et al.  Current status of the AMOEBA polarizable force field. , 2010, The journal of physical chemistry. B.

[4]  Huan-Xiang Zhou,et al.  Comparison of calculation and experiment implicates significant electrostatic contributions to the binding stability of barnase and barstar. , 2003, Biophysical journal.

[5]  Jacek Otlewski,et al.  Trypsin specificity as elucidated by LIE calculations, X‐ray structures, and association constant measurements , 2004, Protein science : a publication of the Protein Society.

[6]  A. Cavalli,et al.  Single-molecule pulling simulations can discern active from inactive enzyme inhibitors. , 2010, Journal of the American Chemical Society.

[7]  J J Baldwin,et al.  Application of the three-dimensional structures of protein target molecules in structure-based drug design. , 1994, Journal of medicinal chemistry.

[8]  Michael R. Shirts,et al.  Equilibrium free energies from nonequilibrium measurements using maximum-likelihood methods. , 2003, Physical review letters.

[9]  Michael K. Gilson,et al.  Theory of Free Energy and Entropy in Noncovalent Binding , 2009 .

[10]  A. Warshel,et al.  What are the dielectric “constants” of proteins and how to validate electrostatic models? , 2001, Proteins.

[11]  J. Janc,et al.  A novel serine protease inhibition motif involving a multi-centered short hydrogen bonding network at the active site. , 2002, Journal of molecular biology.

[12]  Tingjun Hou,et al.  Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking , 2011, J. Comput. Chem..

[13]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[14]  A. Mark,et al.  Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , 1994 .

[15]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

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

[17]  M Karplus,et al.  The contribution of vibrational entropy to molecular association. The dimerization of insulin. , 1994, Journal of molecular biology.

[18]  Thomas A. Halgren,et al.  The representation of van der Waals (vdW) interactions in molecular mechanics force fields: potential form, combination rules, and vdW parameters , 1992 .

[19]  Ronald M. Levy,et al.  The SGB/NP hydration free energy model based on the surface generalized born solvent reaction field and novel nonpolar hydration free energy estimators , 2002, J. Comput. Chem..

[20]  Jan H. Jensen,et al.  Very fast prediction and rationalization of pKa values for protein–ligand complexes , 2008, Proteins.

[21]  M. L. Connolly.  Analytical molecular surface calculation , 1983 .

[22]  W. L. Jorgensen.  The Many Roles of Computation in Drug Discovery , 2004, Science.

[23]  A. Sali,et al.  Modeller: generation and refinement of homology-based protein structure models. , 2003, Methods in enzymology.

[24]  D. M. Ryan,et al.  Rational design of potent sialidase-based inhibitors of influenza virus replication , 1993, Nature.

[25]  David A. Case,et al.  Effective Born radii in the generalized Born approximation: The importance of being perfect , 2002, J. Comput. Chem..

[26]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[27]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[28]  Gisbert Schneider,et al.  Computer-based de novo design of drug-like molecules , 2005, Nature Reviews Drug Discovery.

[29]  R. Hilgenfeld,et al.  Modern Methods of Drug Discovery , 2003, EXS.

[30]  N M Dean,et al.  Functional genomics and target validation approaches using antisense oligonucleotide technology. , 2001, Current opinion in biotechnology.

[31]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[32]  D. Case.  Normal mode analysis of protein dynamics , 1994 .

[33]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

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

[35]  Charles C. Persinger,et al.  How to improve R&D productivity: the pharmaceutical industry's grand challenge , 2010, Nature Reviews Drug Discovery.

[36]  C. Tanford,et al.  Empirical correlation between hydrophobic free energy and aqueous cavity surface area. , 1974, Proceedings of the National Academy of Sciences of the United States of America.

[37]  Jaime Rubio-Martinez,et al.  Comparative Evaluation of MMPBSA and XSCORE To Compute Binding Free Energy in XIAP-Peptide Complexes , 2007, J. Chem. Inf. Model..

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

[39]  Jeremy C. Smith,et al.  Can the calculation of ligand binding free energies be improved with continuum solvent electrostatics and an ideal‐gas entropy correction? , 2002, J. Comput. Chem..

[40]  Sapan Shah,et al.  Non‐Boltzmann thermodynamic integration (NBTI) for macromolecular systems: Relative free energy of binding of trypsin to benzamidine and benzylamine , 1999, Proteins.

[41]  Sandeep Grover,et al.  Topical dorzolamide for the treatment of cystoid macular edema in patients with retinitis pigmentosa. , 2006, American journal of ophthalmology.

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

[43]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[44]  Scott P. Brown,et al.  Large-scale application of high-throughput molecular mechanics with Poisson-Boltzmann surface area for routine physics-based scoring of protein-ligand complexes. , 2009, Journal of medicinal chemistry.

[45]  K. Dill,et al.  Binding of small-molecule ligands to proteins: "what you see" is not always "what you get". , 2009, Structure.

[46]  B. Shoichet,et al.  Molecular docking and ligand specificity in fragment-based inhibitor discovery. , 2009, Nature chemical biology.

[47]  B. Roux,et al.  Implicit solvent models. , 1999, Biophysical chemistry.

[48]  B Honig,et al.  Correlating solvation free energies and surface tensions of hydrocarbon solutes. , 1994, Biophysical chemistry.

[49]  William L Jorgensen,et al.  Efficient drug lead discovery and optimization. , 2009, Accounts of chemical research.

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

[51]  Cen Gao,et al.  Accounting for ligand conformational restriction in calculations of protein-ligand binding affinities. , 2010, Biophysical journal.

[52]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[53]  Ray Luo,et al.  Accelerated Poisson–Boltzmann calculations for static and dynamic systems , 2002, J. Comput. Chem..

[54]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[55]  C. E. Peishoff,et al.  A critical assessment of docking programs and scoring functions. , 2006, Journal of medicinal chemistry.

[56]  E. M.,et al.  Statistical Mechanics , 2021, Manual for Theoretical Chemistry.

[57]  Michael G. Lerner,et al.  Binding MOAD (Mother Of All Databases) , 2005, Proteins.

[58]  Nathan A. Baker,et al.  Assessing implicit models for nonpolar mean solvation forces: the importance of dispersion and volume terms. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[59]  A. Szabó,et al.  Langevin modes of macromolecules , 1986 .

[60]  M. Congreve,et al.  Recent developments in fragment-based drug discovery. , 2008, Journal of medicinal chemistry.

[61]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[62]  B. Druker,et al.  Specific Targeted Therapy of Chronic Myelogenous Leukemia with Imatinib , 2003, Pharmacological Reviews.

[63]  Jay W Ponder,et al.  Polarizable Atomic Multipole Solutes in a Generalized Kirkwood Continuum. , 2007, Journal of chemical theory and computation.

[64]  Thomas Simonson,et al.  A POISSON-BOLTZMANN STUDY OF CHARGE INSERTION IN AN ENZYME ACTIVE SITE : THE EFFECT OF DIELECTRIC RELAXATION , 1999 .

[65]  Tingjun Hou,et al.  Assessing the Performance of the MM/PBSA and MM/GBSA Methods. 1. The Accuracy of Binding Free Energy Calculations Based on Molecular Dynamics Simulations , 2011, J. Chem. Inf. Model..

[66]  Laurent Vial,et al.  Dynamic combinatorial chemistry. , 2006, Chemical reviews.

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

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

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

[70]  Shaomeng Wang,et al.  An Extensive Test of 14 Scoring Functions Using the PDBbind Refined Set of 800 Protein-Ligand Complexes , 2004, J. Chem. Inf. Model..

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

[72]  R Talhout,et al.  Thermodynamic analysis of binding of p-substituted benzamidines to trypsin. , 2001, European journal of biochemistry.

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

[74]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[75]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[76]  William L. Jorgensen,et al.  Drug discovery: Pulled from a protein's embrace , 2010, Nature.

[77]  W. C. Still,et al.  The GB/SA Continuum Model for Solvation. A Fast Analytical Method for the Calculation of Approximate Born Radii , 1997 .

[78]  Pengyu Y. Ren,et al.  Calculation of protein–ligand binding free energy by using a polarizable potential , 2008, Proceedings of the National Academy of Sciences.

[79]  D. Case,et al.  Langevin modes of macromolecules: Applications to crambin and DNA hexamers , 1990, Biopolymers.

[80]  Guohui Li,et al.  Trypsin‐ligand binding free energies from explicit and implicit solvent simulations with polarizable potential , 2009, J. Comput. Chem..

[81]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .