Complementarity between in Silico and Biophysical Screening Approaches in Fragment-Based Lead Discovery against the A2A Adenosine Receptor

Fragment-based lead discovery (FBLD) is becoming an increasingly important method in drug development. We have explored the potential to complement NMR-based biophysical screening of chemical libraries with molecular docking in FBLD against the A(2A) adenosine receptor (A(2A)AR), a drug target for inflammation and Parkinson's disease. Prior to an NMR-based screen of a fragment library against the A(2A)AR, molecular docking against a crystal structure was used to rank the same set of molecules by their predicted affinities. Molecular docking was able to predict four out of the five orthosteric ligands discovered by NMR among the top 5% of the ranked library, suggesting that structure-based methods could be used to prioritize among primary hits from biophysical screens. In addition, three fragments that were top-ranked by molecular docking, but had not been picked up by the NMR-based method, were demonstrated to be A(2A)AR ligands. While biophysical approaches for fragment screening are typically limited to a few thousand compounds, the docking screen was extended to include 328,000 commercially available fragments. Twenty-two top-ranked compounds were tested in radioligand binding assays, and 14 of these were A(2A)AR ligands with K(i) values ranging from 2 to 240 μM. Optimization of fragments was guided by molecular dynamics simulations and free energy calculations. The results illuminate strengths and weaknesses of molecular docking and demonstrate that this method can serve as a valuable complementary tool to biophysical screening in FBLD.

[1]  Ajay N. Jain,et al.  Recommendations for evaluation of computational methods , 2008, J. Comput. Aided Mol. Des..

[2]  R. Stevens,et al.  Structural Basis for Allosteric Regulation of GPCRs by Sodium Ions , 2012, Science.

[3]  Olan Dolezal,et al.  Parallel Screening of Low Molecular Weight Fragment Libraries , 2013, Journal of biomolecular screening.

[4]  Anthony Nicholls,et al.  What do we know and when do we know it? , 2008, J. Comput. Aided Mol. Des..

[5]  Leigh A. Stoddart,et al.  Fragment Screening at Adenosine-A3 Receptors in Living Cells Using a Fluorescence-Based Binding Assay , 2012, Chemistry & biology.

[6]  Nathan Robertson,et al.  Article pubs.acs.org/jmc Identification of Novel Adenosine A 2A Receptor Antagonists by Virtual Screening , 2022 .

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

[8]  Avner Schlessinger,et al.  Ligand Discovery from a Dopamine D3 Receptor Homology Model and Crystal Structure , 2011, Nature chemical biology.

[9]  B. Shoichet,et al.  Hierarchical docking of databases of multiple ligand conformations. , 2005, Current topics in medicinal chemistry.

[10]  Jan Schultz,et al.  Integration of fragment screening and library design. , 2007, Drug discovery today.

[11]  Michael M. Mysinger,et al.  Automated Docking Screens: A Feasibility Study , 2009, Journal of medicinal chemistry.

[12]  P. K. Smith,et al.  Measurement of protein using bicinchoninic acid. , 1985, Analytical biochemistry.

[13]  R. Stevens,et al.  Crystal structure-based virtual screening for fragment-like ligands of the human histamine H(1) receptor. , 2011, Journal of medicinal chemistry.

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

[15]  Andrew L. Hopkins,et al.  Screening for GPCR Ligands Using Surface Plasmon Resonance , 2011, ACS medicinal chemistry letters.

[16]  I. Kuntz,et al.  Ligand solvation in molecular docking , 1999, Proteins.

[17]  Amedeo Caflisch,et al.  Library screening by fragment‐based docking , 2009, Journal of molecular recognition : JMR.

[18]  Ana M Sebastião,et al.  Adenosine receptors and the central nervous system. , 2009, Handbook of experimental pharmacology.

[19]  Vadim Cherezov,et al.  Diversity and modularity of G protein-coupled receptor structures. , 2012, Trends in pharmacological sciences.

[20]  Ruben Abagyan,et al.  Structure based prediction of subtype-selectivity for adenosine receptor antagonists , 2011, Neuropharmacology.

[21]  Kenneth Jones,et al.  Use of the X-ray structure of the beta2-adrenergic receptor for drug discovery. Part 2: Identification of active compounds. , 2008, Bioorganic & medicinal chemistry letters.

[22]  Michael J. Keiser,et al.  Complementarity Between a Docking and a High-Throughput Screen in Discovering New Cruzain Inhibitors , 2010, Journal of medicinal chemistry.

[23]  Albert C. Pan,et al.  Pathway and mechanism of drug binding to G-protein-coupled receptors , 2011, Proceedings of the National Academy of Sciences.

[24]  Arieh Warshel,et al.  A surface constrained all‐atom solvent model for effective simulations of polar solutions , 1989 .

[25]  Ryan G. Coleman,et al.  ZINC: A Free Tool to Discover Chemistry for Biology , 2012, J. Chem. Inf. Model..

[26]  Anakarina B. Bedran-Russo,et al.  Identification of active compounds from a PA-based dentin collagen cross-linker , 2010 .

[27]  Gregg Siegal,et al.  Fragment screening of stabilized G-protein-coupled receptors using biophysical methods. , 2011, Methods in enzymology.

[28]  A. Leslie,et al.  Agonist-bound adenosine A2A receptor structures reveal common features of GPCR activation , 2011, Nature.

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

[30]  I. Kuntz,et al.  The maximal affinity of ligands. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[31]  Jonathan S. Mason,et al.  Progress in Structure Based Drug Design for G Protein-Coupled Receptors , 2011, Journal of medicinal chemistry.

[32]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[33]  Saskia Nijmeijer,et al.  Virtual Fragment Screening: Discovery of Histamine H3 Receptor Ligands Using Ligand-Based and Protein-Based Molecular Fingerprints , 2012, J. Chem. Inf. Model..

[34]  B. Fredholm,et al.  International Union of Basic and Clinical Pharmacology. LXXXI. Nomenclature and Classification of Adenosine Receptors—An Update , 2011, Pharmacological Reviews.

[35]  A. Caflisch,et al.  Efficient electrostatic solvation model for protein‐fragment docking , 2001, Proteins.

[36]  Gregg Siegal,et al.  Fragment screening of GPCRs using biophysical methods: identification of ligands of the adenosine A(2A) receptor with novel biological activity. , 2012, ACS chemical biology.

[37]  C. Higgs,et al.  Hydration Site Thermodynamics Explain SARs for Triazolylpurines Analogues Binding to the A2A Receptor. , 2010, ACS medicinal chemistry letters.

[38]  A. Hopkins,et al.  Ligand efficiency: a useful metric for lead selection. , 2004, Drug discovery today.

[39]  K. Jacobson,et al.  Structure-activity relationships of 8-styrylxanthines as A2-selective adenosine antagonists. , 1993, Journal of medicinal chemistry.

[40]  B. Shoichet,et al.  Flexible ligand docking using conformational ensembles , 1998, Protein science : a publication of the Protein Society.

[41]  Roderick E Hubbard,et al.  Informatics and modeling challenges in fragment-based drug discovery. , 2007, Current opinion in drug discovery & development.

[42]  Brian K. Shoichet,et al.  Structure-Based Discovery of A2A Adenosine Receptor Ligands , 2010, Journal of medicinal chemistry.

[43]  Peter Kolb,et al.  Limits of Ligand Selectivity from Docking to Models: In Silico Screening for A1 Adenosine Receptor Antagonists , 2012, PloS one.

[44]  M. Congreve,et al.  Structure of the adenosine A(2A) receptor in complex with ZM241385 and the xanthines XAC and caffeine. , 2011, Structure.

[45]  Yoko Shibata,et al.  Co-evolving stability and conformational homogeneity of the human adenosine A2a receptor , 2008, Proceedings of the National Academy of Sciences.

[46]  J. Ballesteros,et al.  [19] Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors , 1995 .

[47]  Gregg Siegal,et al.  TINS, target immobilized NMR screening: an efficient and sensitive method for ligand discovery. , 2005, Chemistry & biology.

[48]  John P. Overington,et al.  ChEMBL: a large-scale bioactivity database for drug discovery , 2011, Nucleic Acids Res..

[49]  A. P. Ijzerman,et al.  Allosteric modulation of adenosine receptors , 2000, European Neuropsychopharmacology.

[50]  Andrew J. Woodhead,et al.  Recent developments in fragment-based drug discovery. , 2008, Journal of medicinal chemistry.

[51]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[52]  B. Honig,et al.  A rapid finite difference algorithm, utilizing successive over‐relaxation to solve the Poisson–Boltzmann equation , 1991 .

[53]  Arieh Warshel,et al.  A local reaction field method for fast evaluation of long‐range electrostatic interactions in molecular simulations , 1992 .

[54]  J. Åqvist,et al.  Q: a molecular dynamics program for free energy calculations and empirical valence bond simulations in biomolecular systems. , 1998, Journal of molecular graphics & modelling.

[55]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[56]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[57]  Monya Baker,et al.  Fragment-based lead discovery grows up , 2012, Nature Reviews Drug Discovery.

[58]  Jens Carlsson,et al.  Combining docking, molecular dynamics and the linear interaction energy method to predict binding modes and affinities for non-nucleoside inhibitors to HIV-1 reverse transcriptase. , 2008, Journal of medicinal chemistry.

[59]  R. Stevens,et al.  The 2.6 Angstrom Crystal Structure of a Human A2A Adenosine Receptor Bound to an Antagonist , 2008, Science.

[60]  Martin Almlöf,et al.  Free energy calculations and ligand binding. , 2003, Advances in protein chemistry.

[61]  Donald G. Truhlar,et al.  MODEL FOR AQUEOUS SOLVATION BASED ON CLASS IV ATOMIC CHARGES AND FIRST SOLVATION SHELL EFFECTS , 1996 .

[62]  I. Kuntz,et al.  Automated docking with grid‐based energy evaluation , 1992 .

[63]  Brian K. Shoichet,et al.  Rapid Context-Dependent Ligand Desolvation in Molecular Docking , 2010, J. Chem. Inf. Model..

[64]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

[65]  J. Wess,et al.  Site-directed Mutagenesis Identifies Residues Involved in Ligand Recognition in the Human A2a Adenosine Receptor (*) , 1995, The Journal of Biological Chemistry.

[66]  Ruben Abagyan,et al.  Structure-based discovery of novel chemotypes for adenosine A(2A) receptor antagonists. , 2010, Journal of medicinal chemistry.

[67]  P. Hajduk,et al.  A decade of fragment-based drug design: strategic advances and lessons learned , 2007, Nature Reviews Drug Discovery.

[68]  Michael M. Mysinger,et al.  Structure-based ligand discovery for the protein–protein interface of chemokine receptor CXCR4 , 2012, Proceedings of the National Academy of Sciences.

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

[70]  M. Blackburn,et al.  Adenosine receptors and inflammation. , 2009, Handbook of experimental pharmacology.

[71]  Donald G. Truhlar,et al.  New Class IV Charge Model for Extracting Accurate Partial Charges from Wave Functions , 1998 .

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

[73]  Ruben Abagyan,et al.  Optimization of adenosine 5'-carboxamide derivatives as adenosine receptor agonists using structure-based ligand design and fragment screening. , 2012, Journal of medicinal chemistry.

[74]  I. Kuntz,et al.  Matching chemistry and shape in molecular docking. , 1993, Protein engineering.

[75]  O. Berger,et al.  Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. , 1997, Biophysical journal.

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

[77]  Jean-Louis Reymond,et al.  Enumeration of 166 Billion Organic Small Molecules in the Chemical Universe Database GDB-17 , 2012, J. Chem. Inf. Model..