Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Skip to main content

Advertisement

Enabling the hypothesis-driven prioritization of ligand candidates in big databases: Screenlamp and its application to GPCR inhibitor discovery for invasive species control

  • Published:
Journal of Computer-Aided Molecular Design Aims and scope Submit manuscript

Abstract

While the advantage of screening vast databases of molecules to cover greater molecular diversity is often mentioned, in reality, only a few studies have been published demonstrating inhibitor discovery by screening more than a million compounds for features that mimic a known three-dimensional (3D) ligand. Two factors contribute: the general difficulty of discovering potent inhibitors, and the lack of free, user-friendly software to incorporate project-specific knowledge and user hypotheses into 3D ligand-based screening. The Screenlamp modular toolkit presented here was developed with these needs in mind. We show Screenlamp’s ability to screen more than 12 million commercially available molecules and identify potent in vivo inhibitors of a G protein-coupled bile acid receptor within the first year of a discovery project. This pheromone receptor governs sea lamprey reproductive behavior, and to our knowledge, this project is the first to establish the efficacy of computational screening in discovering lead compounds for aquatic invasive species control. Significant enhancement in activity came from selecting compounds based on one of the hypotheses: that matching two distal oxygen groups in the 3D structure of the pheromone is crucial for activity. Six of the 15 most active compounds met these criteria. A second hypothesis—that presence of an alkyl sulfate side chain results in high activity—identified another 6 compounds in the top 10, demonstrating the significant benefits of hypothesis-driven screening.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save

Springer+ Basic
$34.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9

Similar content being viewed by others

Abbreviations

2D:

Two-dimensional

3D:

Three-dimensional

3kPZS:

3-keto petromyzonol sulfate

3sPZS:

Trisulfated petromyzonol sulfate

GLL:

GPCR Ligand Library

GPCR:

G protein-coupled receptor

Mgβ1AR:

β1-Adrenergic receptor from Meleagris gallopavo

PAIN:

Pan-assay interference compound

PZS:

Petromyzonol sulfate

SQL:

Structured Query Language

SLOR1:

Sea lamprey olfactory receptor 1

TFM:

Trifluoromethyl nitrophenol

EOG:

Electro-olfactogram

References

  1. Shoichet BK (2004) Virtual screening of chemical libraries. Nature 432:862–865

    Article  CAS  Google Scholar 

  2. Zavodszky MI, Rohatgi A, Van Voorst JR, Yan H, Kuhn LA (2009) Scoring ligand similarity in structure-based virtual screening. J Mol Recognit 22:280–292

    Article  CAS  Google Scholar 

  3. Ferrara P, Gohlke H, Price DJ, Klebe G, Brooks CL (2004) Assessing scoring functions for protein-ligand interactions. J Med Chem 47:3032–3047

    Article  CAS  Google Scholar 

  4. Cozzini P, Kellogg GE, Spyrakis F, Abraham DJ, Costantino G, Emerson A, Fanelli F, Gohlke H, Kuhn LA, Morris GM, Orozco M, Pertinhez TA, Rizzi M, Sotriffer CA (2008) Target flexibility: an emerging consideration in drug discovery and design. J Med Chem 51:6237–6255

    Article  CAS  Google Scholar 

  5. Capuccini M, Ahmed L, Schaal W, Laure E, Spjuth O (2017) Large-scale virtual screening on public cloud resources with Apache Spark. J Cheminform 9:15

    Article  Google Scholar 

  6. Merz KM Jr. (2010) Limits of free energy computation for protein-ligand interactions. J Chem Theory Comput 6:1769–1776

    Article  CAS  Google Scholar 

  7. Hou T, Wang J, Li Y, Wang W (2011) Assessing the performance of the MM/PBSA and MM/GBSA methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations. J Chem Inf Model 51:69–82

    Article  CAS  Google Scholar 

  8. Drwal MN, Griffith R (2013) Combination of ligand- and structure-based methods in virtual screening. Drug Discov Today Technol 10:395–401

    Article  Google Scholar 

  9. McGaughey GB, Sheridan RP, Bayly CI, Culberson JC, Kreatsoulas C, Lindsley S, Maiorov V, Truchon J-F, Cornell WD (2007) Comparison of topological, shape, and docking methods in virtual screening. J Chem Inf Model 47:1504–1519

    Article  CAS  Google Scholar 

  10. Hawkins PCD, Skillman AG, Nicholls A (2007) Comparison of shape-matching and docking as virtual screening tools. J Med Chem 50:74–82

    Article  CAS  Google Scholar 

  11. Hu G, Kuang G, Xiao W, Li W, Liu G, Tang Y (2012) Performance evaluation of 2D fingerprint and 3D shape similarity methods in virtual screening. J Chem Inf Model 52:1103–1113

    Article  CAS  Google Scholar 

  12. Rush TS III, Grant JA, Mosyak L, Nicholls A (2005) A shape-based 3-D scaffold hopping method and its application to a bacterial protein-protein interaction. J Med Chem 1489–1495

  13. Muegge I, Mukherjee P (2016) An overview of molecular fingerprint similarity search in virtual screening. Expert Opin Drug Discov 11:137–148

    Article  CAS  Google Scholar 

  14. Zoete V, Daina A, Bovigny C, Michielin O (2016) SwissSimilarity: a web tool for low to ultra high throughput ligand-based virtual screening. J Chem Inf Model 56:1399–1404

    Article  CAS  Google Scholar 

  15. Dixon SL, Smondyrev AM, Knoll EH, Rao SN, Shaw DE, Friesner RA (2006) PHASE: a new engine for pharmacophore perception, 3D QSAR model development, and 3D database screening: 1. Methodology and preliminary results. J Comput Aided Mol Des 20:647–671

    Article  CAS  Google Scholar 

  16. Nagamine N, Shirakawa T, Minato Y, Torii K, Kobayashi H, Imoto M, Sakakibara Y (2009) Integrating statistical predictions and experimental verifications for enhancing protein-chemical interaction predictions in virtual screening. PLoS Comput Biol 5:e1000397

    Article  Google Scholar 

  17. Koes D, Khoury K, Huang Y, Wang W, Bista M, Popowicz GM, Wolf S, Holak TA, Dömling A, Camacho CJ (2012) Enabling large-scale design, synthesis and validation of small molecule protein-protein antagonists. PLoS ONE 7:e32839

    Article  CAS  Google Scholar 

  18. Miller BR, Roitberg AE (2013) Design of e-pharmacophore models using compound fragments for the trans-sialidase of Trypanosoma cruzi: screening for novel inhibitor scaffolds. J Mol Graph Model 45:84–97

    Article  CAS  Google Scholar 

  19. Murgueitio MS, Henneke P, Glossmann H, Santos-Sierra S, Wolber G (2014) Prospective virtual screening in a sparse data scenario: design of small-molecule TLR2 antagonists. ChemMedChem 9:813–822

    Article  CAS  Google Scholar 

  20. Almela MJ, Lozano S, Lelièvre J, Colmenarejo G, Coterón JM, Rodrigues J, Gonzalez C, Herreros E (2015) A new set of chemical starting points with Plasmodium falciparum transmission-blocking potential for antimalarial drug discovery. PLoS ONE 10:e0135139

    Article  Google Scholar 

  21. Allen BK, Mehta S, Ember SWJ, Schonbrunn E, Ayad N, Schürer SC (2015) Large-scale computational screening identifies first in class multitarget inhibitor of EGFR kinase and BRD4. Sci Rep 5:16924.

    Article  CAS  Google Scholar 

  22. Mirza SB, Salmas RE, Fatmi MQ, Durdagi S (2016) Virtual screening of eighteen million compounds against dengue virus: Combined molecular docking and molecular dynamics simulations study. J Mol Graph Model 66:99–107

    Article  CAS  Google Scholar 

  23. Johnson DK, Karanicolas J (2016) Ultra-high-throughput structure-based virtual screening for small-molecule inhibitors of protein-protein interactions. J Chem Inf Model 56:399–411

    Article  CAS  Google Scholar 

  24. Naylor E, Arredouani A, Vasudevan SR, Lewis AM, Parkesh R, Mizote A, Rosen D, Thomas JM, Izumi M, Ganesan A et al. (2009) Identification of a chemical probe for NAADP by virtual screening. Nat Chem Biol 5:220–226

    Article  CAS  Google Scholar 

  25. LaLonde JM, Elban MA, Courter JR, Sugawara A, Soeta T, Madani N, Princiotto AM, Do Kwon Y, Kwong PD, Schön A et al (2011) Design, synthesis and biological evaluation of small molecule inhibitors of CD4-gp120 binding based on virtual screening. Bioorg Med Chem 19:91–101

    Article  CAS  Google Scholar 

  26. Irwin JJ, Shoichet BK (2005) ZINC—a free database of commercially available compounds for virtual screening. J Chem Inf Model 45:177–182

    Article  CAS  Google Scholar 

  27. Li W, Scott AP, Siefkes MJ, Yan H, Liu Q, Yun S-S, Gage DA (2002) Bile acid secreted by male sea lamprey that acts as a sex pheromone. Science 296:138–141

    Article  CAS  Google Scholar 

  28. Hansen GJA, Jones ML (2008) A rapid assessment approach to prioritizing streams for control of Great Lakes sea lampreys (Petromyzon marinus): a case study in adaptive management. Can J Fish Aquat Sci 65:2471–2484

    Article  Google Scholar 

  29. Lucas MC, Bubb DH, Jang MH, Ha K, Masters JEG (2009) Availability of and access to critical habitats in regulated rivers: effects of low-head barriers on threatened lampreys. Freshw Biol 54:621–634

    Article  Google Scholar 

  30. McDonald DG, Kolar CS (2007) Research to guide the use of lampricides for controlling sea lamprey. J Great Lakes Res 33:20–34

    Article  CAS  Google Scholar 

  31. Scott WB, Crossman EJ (1998) Freshwater fishes of Canada. Galt House Pub., Oakville

    Google Scholar 

  32. Becker GC (1983) Fishes of Wisconsin. University of Wisconsin Press, Madison

    Google Scholar 

  33. Boogaard MA, Bills TD, Johnson DA (2003) Acute toxicity of TFM and a TFM/niclosamide mixture to selected species of fish, including lake sturgeon (Acipenser fulvescens) and mudpuppies (Necturus maculosus), in laboratory and field exposures. J Great Lakes Res 29:529–541

    Article  CAS  Google Scholar 

  34. Johnson NS, Yun S-S, Thompson HT, Brant CO, Li W (2009) A synthesized pheromone induces upstream movement in female sea lamprey and summons them into traps. Proc Natl Acad Sci 106:1021–1026

    Article  CAS  Google Scholar 

  35. Johnson NS, Siefkes MJ, Wagner CM, Bravener G, Steeves T, Twohey M, Li W (2015) Factors influencing capture of invasive sea lamprey in traps baited with a synthesized sex pheromone component. J Chem Ecol 41:913–923

    Article  CAS  Google Scholar 

  36. Libants S, Carr K, Wu H, Teeter JH, Chung-Davidson Y-W, Zhang Z, Wilkerson C, Li W (2009) The sea lamprey Petromyzon marinus genome reveals the early origin of several chemosensory receptor families in the vertebrate lineage. BMC Evol Biol 9:180

    Article  Google Scholar 

  37. Katritch V, Cherezov V, Stevens RC (2012) Diversity and modularity of G protein-coupled receptor structures. Trends Pharmacol Sci 33:17–27

    Article  CAS  Google Scholar 

  38. Lundstrom K (2009) An overview on GPCRs and drug discovery: structure-based drug design and structural biology on GPCRs. Methods Mol Biol 51–66

  39. Chang AJ, Ortega FE, Riegler J, Madison DV, Krasnow MA (2015) Oxygen control of breathing by an olfactory receptor activated by lactate. Nature 527:240–244

    Article  CAS  Google Scholar 

  40. Bessac BF, Jordt S-E (2010) Sensory detection and responses to toxic gases: mechanisms, health effects, and countermeasures. Proc Am Thorac Soc 7:269–277

    Article  CAS  Google Scholar 

  41. Hussain A, Saraiva LR, Ferrero DM, Ahuja G, Krishna VS, Liberles SD, Korsching SI (2013) High-affinity olfactory receptor for the death-associated odor cadaverine. Proc Natl Acad Sci USA 110:19579–19584

    Article  CAS  Google Scholar 

  42. Milligan G, Ulven T, Murdoch H, Hudson BD (2014) G-protein-coupled receptors for free fatty acids: nutritional and therapeutic targets. Br J Nutr 111:3–7

    Article  Google Scholar 

  43. Liberles SD (2014) Mammalian pheromones. Annu Rev Physiol 76:151–175

    Article  CAS  Google Scholar 

  44. Niimura Y, Niimura Y (2009) On the origin and evolution of vertebrate olfactory receptor genes: Comparative genome analysis among 23 chordate species. Genome Biol Evol 1:34–44

    Article  Google Scholar 

  45. Benton R, Sachse S, Michnick SW, Vosshall LB (2006) Atypical membrane topology and heteromeric function of Drosophila odorant receptors in vivo. PLoS Biol 4:240–257

    Article  CAS  Google Scholar 

  46. Kain P, Boyle SM, Tharadra SK, Guda T, Pham C, Dahanukar A, Ray A (2013) Odour receptors and neurons for DEET and new insect repellents. Nature 502:507–512

    Article  CAS  Google Scholar 

  47. Berman HM (2000) The Protein Data Bank. Nucleic Acids Res 28:235–242

    Article  CAS  Google Scholar 

  48. Sander C, Schneider R (1991) Database of homology-derived protein structures and the structural meaning of sequence alignment. Proteins Struct Funct Bioinform 9:56–68

    Article  CAS  Google Scholar 

  49. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215:403–410

    Article  CAS  Google Scholar 

  50. Rosenbaum DM, Rasmussen SGF, Kobilka BK (2009) The structure and function of G-protein-coupled receptors. Nature 459:356–363

    Article  CAS  Google Scholar 

  51. Warne T, Serrano-Vega MJ, Baker JG, Moukhametzianov R, Edwards PC, Henderson R, Leslie AGW, Tate CG, Schertler GFX (2008) Structure of a beta1-adrenergic G-protein-coupled receptor. Nature 454:486–491

    Article  CAS  Google Scholar 

  52. Webb B, Sali A (2014) Protein structure modeling with MODELLER. In: Kihara D (ed) Protein Structure Prediction. Methods Mol Biol. Humana Press, New York, NY

    Google Scholar 

  53. Hawkins PCD, Skillman AG, Warren GL, Ellingson BA, Stahl MT (2010) Conformer generation with OMEGA: algorithm and validation using high quality structures from the Protein Databank and Cambridge Structural Database. J Chem Inf Model 50:572–584

    Article  CAS  Google Scholar 

  54. Hawkins PCD, Nicholls A (2012) Conformer generation with OMEGA: learning from the data set and the analysis of failures. J Chem Inf Model 52:2919–2936

    Article  CAS  Google Scholar 

  55. Zavodszky MI, Sanschagrin PC, Korde RS, Kuhn LA (2002) Distilling the essential features of a protein surface for improving protein-ligand docking, scoring, and virtual screening. J Comput Aided Mol Des 16:883–902

    Article  CAS  Google Scholar 

  56. Van Der Walt S, Colbert SC, Varoquaux G (2011) The NumPy array: a structure for efficient numerical computation. Comput Sci Eng 13:22–30

    Article  Google Scholar 

  57. McKinney W (2010) Data structures for statistical computing in Python. In: Millman J, van der Walt S (eds) Proceedings of the 9th Python in science conference, pp 51–56

  58. Raschka S (2017) BioPandas: working with molecular structures in pandas DataFrames. J Open Source Softw 2:1–3

    Article  Google Scholar 

  59. Halgren TA (1996) Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94. J Comput Chem 17:490–519

    Article  CAS  Google Scholar 

  60. Jakalian A, Jack DB, Bayly CI (2002) Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation. J Comput Chem 23:1623–1641

    Article  CAS  Google Scholar 

  61. Sheridan RP, McGaughey GB, Cornell WD (2008) Multiple protein structures and multiple ligands: effects on the apparent goodness of virtual screening results. J Comput Aided Mol Des 22:257–265

    Article  CAS  Google Scholar 

  62. Chamberlin DD, Boyce RF (1974) SEQUEL: a structured English query language. In: Altshuler G, Rustin R, Plagman B (eds) Proceedings of the 1974 ACM SIGFIDET (now SIGMOD) workshop on data description, access and control, pp 249–264

  63. Folk M, Heber G, Koziol Q, Pourmal E, Robinson D (2011) An overview of the HDF5 technology suite and its applications. In: Stoyanovich J (ed) Proceedings of the 2011 EDBT/ICDT workshop on array databases, pp 36–47

  64. Mitchell JBO (2014) Machine learning methods in chemoinformatics. Wiley Interdiscip Rev Comput Mol Sci 4:468–481

    Article  CAS  Google Scholar 

  65. Raschka S, Scott AM, Huertas M, Weiming L, Kuhn LA (2018) In: Gore (ed) Methods in molecular biology: computational drug discovery and design. Springer, New York (ISBN: 978-1-4939-7755-0)

    Google Scholar 

  66. Weininger D (1988) SMILES, a chemical language and information system. 1. Introduction to methodology and encoding rules. J Chem Inf Comput Sci 28:31–36

    Article  CAS  Google Scholar 

  67. Gatica EA, Cavasotto CN (2012) Ligand and decoy sets for docking to G protein-coupled receptors. J Chem Inf Model 52:1–6

    Article  CAS  Google Scholar 

  68. Lipinski CA, Lombardo F, Dominy BW, Feeney PJ (1997) Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv Drug Deliv Rev 23:3–25

    Article  CAS  Google Scholar 

  69. Mills JEJ, Dean PM (1996) Three-dimensional hydrogen-bond geometry and probability information from a crystal survey. J Comput Aided Mol Des 10:607–622

    Article  Google Scholar 

  70. Scott JW, Scott-Johnson PE (2002) The electroolfactogram: a review of its history and uses. Microsc Res Tech 58:152–160

    Article  Google Scholar 

  71. Brant CO, Huertas M, Li K, Li W (2016) Mixtures of two bile alcohol sulfates function as a proximity pheromone in sea lamprey. PLoS ONE 11:e0149508

    Article  Google Scholar 

  72. Siefkes MJ, Scott AP, Zielinski B, Yun S-S, Li W (2003) Male sea lampreys, Petromyzon marinus L., excrete a sex pheromone from gill epithelia. Biol Reprod 69:125–132

    Article  CAS  Google Scholar 

  73. DeLano WL (2002) Pymol: an open-source molecular graphics tool. CCP4 Newsl Protein Crystallogr 40:82–92

    Google Scholar 

  74. Hunter JD (2007) Matplotlib: a 2D graphics environment. Comput Sci Eng 9:90–95

    Article  Google Scholar 

  75. Benkert P, Biasini M, Schwede T (2011) Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics 27:343–350

    Article  CAS  Google Scholar 

  76. Zhou H, Zhou Y (2002) Distance-scaled, finite ideal-gas reference state improves structure-derived potentials of mean force for structure selection and stability prediction. Protein Sci 11:2714–2726

    Article  CAS  Google Scholar 

  77. Liu W, Chun E, Thompson AA, Chubukov P, Xu F, Katritch V, Han GW, Roth CB, Heitman LH, IJzerman AP et al (2012) Structural basis for allosteric regulation of GPCRs by sodium ions. Science 337:232–236

    Article  CAS  Google Scholar 

  78. Venkatakrishnan AJ, Deupi X, Lebon G, Tate CG, Schertler GF, Babu MM (2013) Molecular signatures of G-protein-coupled receptors. Nature 494:185–194

    Article  CAS  Google Scholar 

  79. Liu N, Van Voorst JR, Johnston JB, Kuhn LA (2015) CholMine: determinants and prediction of cholesterol and cholate binding across nonhomologous protein structures. J Chem Inf Model 55:747–759

    Article  CAS  Google Scholar 

  80. Ballesteros JA, Weinstein H (1995) Integrated methods for the construction of three-dimensional models and computational probing of structure-function relations in G protein-coupled receptors. Methods Neurosci 25:366–428

    Article  CAS  Google Scholar 

  81. Kabsch W, Sander C (1983) Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 22:2577–2637

    Article  CAS  Google Scholar 

  82. Li W, Sorensen PW, Gallaher DD (1995) The olfactory system of migratory adult sea lamprey (Petromyzon marinus) is specifically and acutely sensitive to unique bile acids released by conspecific larvae. J Gen Physiol 105:569–587

    Article  CAS  Google Scholar 

  83. Burns AC, Sorensen PW, Hoye TR (2011) Synthesis and olfactory activity of unnatural, sulfated 5beta-bile acid derivatives in the sea lamprey (Petromyzon marinus). Steroids 76:291–300

    Article  CAS  Google Scholar 

  84. Hagey LR, Møller PR, Hofmann AF, Krasowski MD (2010) Diversity of bile salts in fish and amphibians: Evolution of a complex biochemical pathway. Physiol Biochem Zool 83:308–321

    Article  CAS  Google Scholar 

  85. Brant CO (2015) Characterization of sea lamprey pheromone components. Ph.D. Thesis. Michigan State University

  86. Chen G, Heim A, Riether D, Yee D, Milgrom Y, Gawinowicz MA, Sames D (2003) Reactivity of functional groups on the protein surface: development of epoxide probes for protein labeling. JACS 125:8130–8133

    Article  CAS  Google Scholar 

  87. Kruse AC, Ring AM, Manglik A, Hu J, Hu K, Eitel K, Hübner H, Pardon E, Valant C, Sexton PM et al (2013) Activation and allosteric modulation of a muscarinic acetylcholine receptor. Nature 504:101–106

    Article  CAS  Google Scholar 

  88. Baell JB, Holloway GA (2010) New substructure filters for removal of pan assay interference compounds (PAINS) from screening libraries and for their exclusion in bioassays. J Med Chem 53:2719–2740

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank Qinghua Yuan for her contributions to the homology modeling of SLOR1 and Stacey Kneeshaw for evaluating protein–ligand energy minimization protocols for SLOR1-3kPZS docking and analyzing charge distributions for matching functional groups. This research was supported by funding from the Great Lakes Fishery Commission from 2012-present (Project ID: 2015_KUH_54031). We gratefully acknowledge OpenEye Scientific Software (Santa Fe, NM) for providing academic licenses for the use of their ROCS, OMEGA, QUACPAC (molcharge), and OEChem toolkit software.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Leslie A. Kuhn.

Additional information

Sebastian Raschka and Anne M. Scott contributed equally to this work.

Electronic supplementary material

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Raschka, S., Scott, A.M., Liu, N. et al. Enabling the hypothesis-driven prioritization of ligand candidates in big databases: Screenlamp and its application to GPCR inhibitor discovery for invasive species control. J Comput Aided Mol Des 32, 415–433 (2018). https://doi.org/10.1007/s10822-018-0100-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s10822-018-0100-7

Keywords