The recent determination of X-ray structures of pharmacologically relevant GPCRs has made these targets accessible to structure-based ligand discovery. Here we explore whether novel chemotypes may be discovered for the A(2A) adenosine receptor, based on complementarity to its recently determined structure. The A(2A) adenosine receptor signals in the periphery and the CNS, with agonists explored as anti-inflammatory drugs and antagonists explored for neurodegenerative diseases. We used molecular docking to screen a 1.4 million compound database against the X-ray structure computationally and tested 20 high-ranking, previously unknown molecules experimentally. Of these 35% showed substantial activity with affinities between 200 nM and 9 microM. For the most potent of these new inhibitors, over 50-fold specificity was observed for the A(2A) versus the related A(1) and A(3) subtypes. These high hit rates and affinities at least partly reflect the bias of commercial libraries toward GPCR-like chemotypes, an issue that we attempt to investigate quantitatively. Despite this bias, many of the most potent new ligands were novel, dissimilar from known ligands, providing new lead structures for modulation of this medically important target.

Structure-based docking screens of large compound libraries have become common in early drug and probe discovery. As computer efficiency has improved and compound libraries have grown, the ability to screen hundreds of millions, and even billions, of compounds has become feasible for modest-sized computer clusters. This allows the rapid and cost-effective exploration and categorization of vast chemical space into a subset enriched with potential hits for a given target. To accomplish this goal at speed, approximations are used that result in undersampling of possible configurations and inaccurate predictions of absolute binding energies. Accordingly, it is important to establish controls, as are common in other fields, to enhance the likelihood of success in spite of these challenges. Here we outline best practices and control docking calculations that help evaluate docking parameters for a given target prior to undertaking a large-scale prospective screen, with exemplification in one particular target, the melatonin receptor, where following this procedure led to direct docking hits with activities in the subnanomolar range. Additional controls are suggested to ensure specific activity for experimentally validated hit compounds. These guidelines should be useful regardless of the docking software used. Docking software described in the outlined protocol (DOCK3.7) is made freely available for academic research to explore new hits for a range of targets.

G protein-coupled receptors (GPCRs) are intensely studied as drug targets and for their role in signaling. With the determination of the first crystal structures, interest in structure-based ligand discovery increased. Unfortunately, for most GPCRs no experimental structures are available. The determination of the D(3) receptor structure and the challenge to the community to predict it enabled a fully prospective comparison of ligand discovery from a modeled structure versus that of the subsequently released crystal structure. Over 3.3 million molecules were docked against a homology model, and 26 of the highest ranking were tested for binding. Six had affinities ranging from 0.2 to 3.1 μM. Subsequently, the crystal structure was released and the docking screen repeated. Of the 25 compounds selected, five had affinities ranging from 0.3 to 3.0 μM. One of the new ligands from the homology model screen was optimized for affinity to 81 nM. The feasibility of docking screens against modeled GPCRs more generally is considered.

Liver receptor homolog 1 (nuclear receptor LRH-1, NR5A2) is an essential regulator of gene transcription, critical for maintenance of cell pluripotency in early development and imperative for the proper functions of the liver, pancreas, and intestines during the adult life. Although physiological hormones of LRH-1 have not yet been identified, crystallographic and biochemical studies demonstrated that LRH-1 could bind regulatory ligands and suggested phosphatidylinositols as potential hormone candidates for this receptor. No synthetic antagonists of LRH-1 are known to date. Here, we identify the first small molecule antagonists of LRH-1 activity. Our search for LRH-1 modulators was empowered by screening of 5.2 million commercially available compounds via molecular docking followed by verification of the top-ranked molecules using in vitro direct binding and transcriptional assays. Experimental evaluation of the predicted ligands identified two compounds that inhibit the transcriptional activity of LRH-1 and diminish the expression of the receptor's target genes. Among the affected transcriptional targets are co-repressor SHP (small heterodimer partner) as well as cyclin E1 (CCNE1) and G0S2 genes that are known to regulate cell growth and proliferation. Treatments of human pancreatic (AsPC-1), colon (HT29), and breast adenocarcinoma cells T47D and MDA-MB-468 with the LRH-1 antagonists resulted in the receptor-mediated inhibition of cancer cell proliferation. Our data suggest that specific antagonists of LRH-1 could be used as specific molecular probes for elucidating the roles of the receptor in different types of malignancies.