Predictive Power of Molecular Dynamics Receptor Structures in Virtual Screening

Molecular dynamics (MD) simulation is a well-established method for understanding protein dynamics. Conformations from unrestrained MD simulations have yet to be assessed for blind virtual screening (VS) by docking. This study presents a critical analysis of the predictive power of MD snapshots to this regard, evaluating two well-characterized systems of varying flexibility in ligand-bound and unbound configurations. Results from such VS predictions are discussed with respect to experimentally determined structures. In all cases, MD simulations provide snapshots that improve VS predictive power over known crystal structures, possibly due to sampling more relevant receptor conformations. Additionally, MD can move conformations previously not amenable to docking into the predictive range.

[1]  Erin S. Bolstad,et al.  In pursuit of virtual lead optimization: The role of the receptor structure and ensembles in accurate docking , 2008, Proteins.

[2]  D. Goodin,et al.  Artificial protein cavities as specific ligand-binding templates: characterization of an engineered heterocyclic cation-binding site that preserves the evolved specificity of the parent protein. , 2002, Journal of molecular biology.

[3]  J. Bajorath,et al.  Docking and scoring in virtual screening for drug discovery: methods and applications , 2004, Nature Reviews Drug Discovery.

[4]  Woody Sherman,et al.  Improving database enrichment through ensemble docking , 2008, J. Comput. Aided Mol. Des..

[5]  J. Mongan,et al.  Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules. , 2004, The Journal of chemical physics.

[6]  William L. Jorgensen,et al.  Discovery of Wild-Type and Y181C Mutant Non-nucleoside HIV-1 Reverse Transcriptase Inhibitors Using Virtual Screening with Multiple Protein Structures , 2009, J. Chem. Inf. Model..

[7]  Markus Christen,et al.  Architecture, implementation and parallelisation of the GROMOS software for biomolecular simulation , 2012, Comput. Phys. Commun..

[8]  Kenneth M. Merz,et al.  Drug Design : Structure-and Ligand-Based Approaches , 2017 .

[9]  Gisbert Schneider,et al.  Virtual screening: an endless staircase? , 2010, Nature Reviews Drug Discovery.

[10]  Tom Fawcett,et al.  An introduction to ROC analysis , 2006, Pattern Recognit. Lett..

[11]  E. De Clercq,et al.  Non-nucleoside reverse transcriptase inhibitors (NNRTIs): past, present, and future. , 2004, Chemistry & biodiversity.

[12]  Chris Oostenbrink,et al.  A biomolecular force field based on the free enthalpy of hydration and solvation: The GROMOS force‐field parameter sets 53A5 and 53A6 , 2004, J. Comput. Chem..

[13]  J. Andrew McCammon,et al.  Mapping the Druggable Allosteric Space of G-Protein Coupled Receptors: a Fragment-Based Molecular Dynamics Approach , 2010, Chemical biology & drug design.

[14]  J. Mccammon,et al.  Computational drug design accommodating receptor flexibility: the relaxed complex scheme. , 2002, Journal of the American Chemical Society.

[15]  Ruben Abagyan,et al.  Recipes for the Selection of Experimental Protein Conformations for Virtual Screening , 2010, J. Chem. Inf. Model..

[16]  D E McRee,et al.  The role of aspartate‐235 in the binding of cations to an artificial cavity at the radical site of cytochrome c peroxidase , 1995, Protein science : a publication of the Protein Society.

[17]  Roger S Armen,et al.  An Evaluation of Explicit Receptor Flexibility in Molecular Docking Using Molecular Dynamics and Torsion Angle Molecular Dynamics. , 2009, Journal of chemical theory and computation.

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

[19]  Chung F Wong,et al.  Flexible ligand-flexible protein docking in protein kinase systems. , 2008, Biochimica et biophysica acta.

[20]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[21]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[22]  Alexander D. MacKerell,et al.  Computational evaluation of protein-small molecule binding. , 2009, Current opinion in structural biology.

[23]  J. Pin,et al.  Virtual screening workflow development guided by the "receiver operating characteristic" curve approach. Application to high-throughput docking on metabotropic glutamate receptor subtype 4. , 2005, Journal of medicinal chemistry.

[24]  Alan E. Mark,et al.  The GROMOS96 Manual and User Guide , 1996 .

[25]  J Andrew McCammon,et al.  (Thermo)dynamic role of receptor flexibility, entropy, and motional correlation in protein-ligand binding. , 2008, ChemPhysChem.

[26]  J Andrew McCammon,et al.  Dynamics, hydration, and motional averaging of a loop-gated artificial protein cavity: the W191G mutant of cytochrome c peroxidase in water as revealed by molecular dynamics simulations. , 2007, Biochemistry.

[27]  Adrian H Elcock,et al.  Structure selection for protein kinase docking and virtual screening: homology models or crystal structures? , 2006, Current protein & peptide science.

[28]  Martha S. Head,et al.  Drug Design: Docking: a domesday report , 2010 .

[29]  Yvonne Jones,et al.  High resolution structures of HIV-1 RT from four RT–inhibitor complexes , 1995, Nature Structural Biology.

[30]  M. M. Fitzgerald,et al.  A ligand-gated, hinged loop rearrangement opens a channel to a buried artificial protein cavity , 1996, Nature Structural Biology.

[31]  B. Shoichet,et al.  Information decay in molecular docking screens against holo, apo, and modeled conformations of enzymes. , 2003, Journal of medicinal chemistry.

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

[33]  Jianping Ding,et al.  Targeting HIV reverse transcriptase for anti-AIDS drug design: structural and biological considerations for chemotherapeutic strategies. , 1996, Drug design and discovery.

[34]  H. Carlson.  Protein flexibility and drug design: how to hit a moving target. , 2002, Current opinion in chemical biology.

[35]  Heather A Carlson,et al.  Protein flexibility and species specificity in structure-based drug discovery: dihydrofolate reductase as a test system. , 2007, Journal of the American Chemical Society.

[36]  H. Berendsen,et al.  A consistent empirical potential for water–protein interactions , 1984 .

[37]  Rommie E. Amaro,et al.  An improved relaxed complex scheme for receptor flexibility in computer-aided drug design , 2008, J. Comput. Aided Mol. Des..

[38]  Jonathan W. Essex,et al.  Ensemble Docking into Multiple Crystallographically Derived Protein Structures: An Evaluation Based on the Statistical Analysis of Enrichments , 2010, J. Chem. Inf. Model..

[39]  B. Shoichet,et al.  Probing molecular docking in a charged model binding site. , 2006, Journal of molecular biology.

[40]  Michael K. Gilson,et al.  Screening Drug-Like Compounds by Docking to Homology Models: A Systematic Study , 2006, J. Chem. Inf. Model..

[41]  Erik De Clercq,et al.  Non‐Nucleoside Reverse Transcriptase Inhibitors (NNRTIs): Past, Present, and Future , 2004, Chemistry & biodiversity.

[42]  Jung-Hsin Lin,et al.  The relaxed complex method: Accommodating receptor flexibility for drug design with an improved scoring scheme. , 2003, Biopolymers.

[43]  X. Barril,et al.  Unveiling the full potential of flexible receptor docking using multiple crystallographic structures. , 2005, Journal of medicinal chemistry.

[44]  S. Ozkan,et al.  A flexible docking scheme to explore the binding selectivity of PDZ domains , 2010, Protein science : a publication of the Protein Society.

[45]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[46]  D. Case,et al.  Rescoring docking hit lists for model cavity sites: predictions and experimental testing. , 2008, Journal of molecular biology.

[47]  Andrew E. Torda,et al.  The GROMOS biomolecular simulation program package , 1999 .

[48]  J Andrew McCammon,et al.  Elucidating the inhibition mechanism of HIV-1 non-nucleoside reverse transcriptase inhibitors through multicopy molecular dynamics simulations. , 2009, Journal of molecular biology.

[49]  X. Zou,et al.  Ensemble docking of multiple protein structures: Considering protein structural variations in molecular docking , 2006, Proteins.

[50]  Markus Christen,et al.  The GROMOS software for biomolecular simulation: GROMOS05 , 2005, J. Comput. Chem..

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