High-Performance Drug Discovery: Computational Screening by Combining Docking and Molecular Dynamics Simulations

Virtual compound screening using molecular docking is widely used in the discovery of new lead compounds for drug design. However, this method is not completely reliable and therefore unsatisfactory. In this study, we used massive molecular dynamics simulations of protein-ligand conformations obtained by molecular docking in order to improve the enrichment performance of molecular docking. Our screening approach employed the molecular mechanics/Poisson-Boltzmann and surface area method to estimate the binding free energies. For the top-ranking 1,000 compounds obtained by docking to a target protein, approximately 6,000 molecular dynamics simulations were performed using multiple docking poses in about a week. As a result, the enrichment performance of the top 100 compounds by our approach was improved by 1.6–4.0 times that of the enrichment performance of molecular dockings. This result indicates that the application of molecular dynamics simulations to virtual screening for lead discovery is both effective and practical. However, further optimization of the computational protocols is required for screening various target proteins.

[1]  K. Sharp,et al.  Accurate Calculation of Hydration Free Energies Using Macroscopic Solvent Models , 1994 .

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

[3]  Christopher W Murray,et al.  Fragment-based lead discovery: leads by design. , 2005, Drug discovery today.

[4]  M L Lamb,et al.  Computational approaches to molecular recognition. , 1997, Current opinion in chemical biology.

[5]  Michael K Gilson,et al.  Evaluating the Accuracy of the Quasiharmonic Approximation. , 2005, Journal of chemical theory and computation.

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

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

[8]  M. Karplus,et al.  Evaluation of the configurational entropy for proteins: application to molecular dynamics simulations of an α-helix , 1984 .

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

[10]  F. Lombardo,et al.  Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. , 2001, Advanced drug delivery reviews.

[11]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[12]  M Webb,et al.  Ultra-High Throughput Screen of Two-Million-Member Combinatorial Compound Collection in a Miniaturized, 1536-Well Assay Format , 2000, Journal of biomolecular screening.

[13]  B. Katz,et al.  Structural basis for selectivity of a small molecule, S1-binding, submicromolar inhibitor of urokinase-type plasminogen activator. , 2000, Chemistry & biology.

[14]  Alessandro Curioni,et al.  New Scoring Functions for Virtual Screening from Molecular Dynamics Simulations with a Quantum‐Refined Force‐Field (QRFF‐MD). Application to Cyclin‐Dependent Kinase 2. , 2006 .

[15]  Maurizio Pellecchia,et al.  A fragment-based approach for the discovery of isoform-specific p38alpha inhibitors. , 2007, ACS chemical biology.

[16]  Emil Alexov,et al.  Rapid grid‐based construction of the molecular surface and the use of induced surface charge to calculate reaction field energies: Applications to the molecular systems and geometric objects , 2002, J. Comput. Chem..

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

[18]  Stephen L Hillis,et al.  Monte Carlo validation of the Dorfman-Berbaum-Metz method using normalized pseudovalues and less data-based model simplification. , 2005, Academic radiology.

[19]  Martin Stahl,et al.  Novel dihydrofolate reductase inhibitors. Structure-based versus diversity-based library design and high-throughput synthesis and screening. , 2003, Journal of medicinal chemistry.

[20]  Jane A. Endicott,et al.  Structure-based design of a potent purine-based cyclin-dependent kinase inhibitor , 2002, Nature Structural Biology.

[21]  A. Pope,et al.  A Homogenous 384-Well High Throughput Screen for Novel Tumor Necrosis Factor Receptor: Ligand Interactions Using Time Resolved Energy Transfer , 1999, Journal of biomolecular screening.

[22]  Agrobacterium tumefaciens-mediated transformation of filamentous fungi , 1998, Nature Biotechnology.

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

[24]  J Andrew McCammon,et al.  Studying enzyme binding specificity in acetylcholinesterase using a combined molecular dynamics and multiple docking approach. , 2002, Journal of the American Chemical Society.

[25]  Hiroshi Habe,et al.  Crystal structure of a histidine-tagged serine hydrolase involved in the carbazole degradation (CarC enzyme). , 2003, Biochemical and biophysical research communications.

[26]  P. Fischer,et al.  Protein structures in virtual screening: a case study with CDK2. , 2006, Journal of medicinal chemistry.

[27]  Andrew R. Leach,et al.  Molecular Complexity and Its Impact on the Probability of Finding Leads for Drug Discovery , 2001, J. Chem. Inf. Comput. Sci..

[28]  Jeffry D Madura,et al.  Relative free energy of binding and binding mode calculations of HIV‐1 RT inhibitors based on dock‐MM‐PB/GS , 2004, Proteins.

[29]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[30]  D S Goodsell,et al.  Automated docking of flexible ligands: Applications of autodock , 1996, Journal of molecular recognition : JMR.

[31]  F. J. Luque,et al.  3D structure of Torpedo californica acetylcholinesterase complexed with huprine X at 2.1 A resolution: kinetic and molecular dynamic correlates. , 2002, Biochemistry.

[32]  Ruben Abagyan,et al.  Comparative study of several algorithms for flexible ligand docking , 2003, J. Comput. Aided Mol. Des..

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

[34]  C E Metz,et al.  Variance-component modeling in the analysis of receiver operating characteristic index estimates. , 1997, Academic radiology.

[35]  William A. Goddard,et al.  The MPSim‐Dock hierarchical docking algorithm: Application to the eight trypsin inhibitor cocrystals , 2005, J. Comput. Chem..

[36]  Thomas Lengauer,et al.  A fast flexible docking method using an incremental construction algorithm. , 1996, Journal of molecular biology.

[37]  Ryutaro Himeno,et al.  A 55 TFLOPS simulation of amyloid-forming peptides from yeast prion Sup35 with the special-purpose computer system MDGRAPE-3 , 2006, SC.

[38]  Todd J. A. Ewing,et al.  DOCK 4.0: Search strategies for automated molecular docking of flexible molecule databases , 2001, J. Comput. Aided Mol. Des..

[39]  Ian H. Witten,et al.  Data mining: practical machine learning tools and techniques with Java implementations , 2002, SGMD.

[40]  Matthew P. Repasky,et al.  Glide: a new approach for rapid, accurate docking and scoring. 1. Method and assessment of docking accuracy. , 2004, Journal of medicinal chemistry.

[41]  Ernst-Walter Knapp,et al.  Conformational entropy of biomolecules: beyond the quasi-harmonic approximation. , 2007, Genome informatics. International Conference on Genome Informatics.

[42]  Renxiao Wang,et al.  The PDBbind database: collection of binding affinities for protein-ligand complexes with known three-dimensional structures. , 2004, Journal of medicinal chemistry.

[43]  K Hamasaki,et al.  A high-throughput fluorescence screen to monitor the specific binding of antagonists to RNA targets. , 1998, Analytical biochemistry.

[44]  K. Berbaum,et al.  Receiver operating characteristic rating analysis. Generalization to the population of readers and patients with the jackknife method. , 1992, Investigative radiology.

[45]  J. Åqvist,et al.  Ligand binding affinities from MD simulations. , 2002, Accounts of chemical research.

[46]  B. Shoichet,et al.  Molecular docking and high-throughput screening for novel inhibitors of protein tyrosine phosphatase-1B. , 2002, Journal of medicinal chemistry.

[47]  Nancy A Obuchowski,et al.  A comparison of the Dorfman–Berbaum–Metz and Obuchowski–Rockette methods for receiver operating characteristic (ROC) data , 2005, Statistics in medicine.

[48]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[49]  Chong-Hwan Chang,et al.  Molecular Recognition of Cyclic Urea HIV-1 Protease Inhibitors* , 1998, The Journal of Biological Chemistry.

[50]  Maurizio Pellecchia,et al.  A fragment-based approach for the discovery of isoform-specific p38? inhibitors , 2007 .

[51]  P. Charifson,et al.  Are free energy calculations useful in practice? A comparison with rapid scoring functions for the p38 MAP kinase protein system. , 2001, Journal of medicinal chemistry.

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

[53]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[54]  B G Benson,et al.  Prevention of chemotherapy-induced alopecia in rats by CDK inhibitors. , 2001, Science.

[55]  J. Irwin,et al.  Benchmarking sets for molecular docking. , 2006, Journal of medicinal chemistry.

[56]  J Andrew McCammon,et al.  Optimized Radii for Poisson-Boltzmann Calculations with the AMBER Force Field. , 2005, Journal of chemical theory and computation.

[57]  R. Glen,et al.  Molecular recognition of receptor sites using a genetic algorithm with a description of desolvation. , 1995, Journal of molecular biology.

[58]  Barry Honig,et al.  Extending the Applicability of the Nonlinear Poisson−Boltzmann Equation: Multiple Dielectric Constants and Multivalent Ions† , 2001 .

[59]  Irwin D Kuntz,et al.  Investigation of neuraminidase-substrate recognition using molecular dynamics and free energy calculations. , 2003, Journal of medicinal chemistry.

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

[61]  Piotr Cieplak,et al.  Molecular dynamics and free energy analyses of cathepsin D-inhibitor interactions: insight into structure-based ligand design. , 2002, Journal of medicinal chemistry.

[62]  Hege S. Beard,et al.  Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. , 2004, Journal of medicinal chemistry.

[63]  T Lengauer,et al.  Two-stage method for protein-ligand docking. , 1999, Journal of medicinal chemistry.

[64]  Maria Kontoyianni,et al.  Evaluation of docking performance: comparative data on docking algorithms. , 2004, Journal of medicinal chemistry.

[65]  C A Roe,et al.  Dorfman-Berbaum-Metz method for statistical analysis of multireader, multimodality receiver operating characteristic data: validation with computer simulation. , 1997, Academic radiology.

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

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

[68]  Richard A. Lewis,et al.  Lessons in molecular recognition: the effects of ligand and protein flexibility on molecular docking accuracy. , 2004, Journal of medicinal chemistry.

[69]  S H Kim,et al.  Exploiting chemical libraries, structure, and genomics in the search for kinase inhibitors. , 1998, Science.