A prospective cross-screening study on G-protein-coupled receptors: lessons learned in virtual compound library design

Structure-Guided Design of G-Protein-Coupled Receptor Polypharmacology

Many diseases are polygenic and can only be treated efficiently with drugs that modulate multiple targets. However, rational design of compounds with multi-target profiles is rarely pursued because it is considered too difficult, in particular if the drug must enter the central nervous system. Here, a structure-based strategy to identify dual-target ligands of G-protein-coupled receptors is presented. We use this approach to design compounds that both antagonize the A_2A adenosine receptor and activate the D₂ dopamine receptor, which have excellent potential as antiparkinson drugs. Atomic resolution models of the receptors guided generation of a chemical library with compounds designed to occupy orthosteric and secondary binding pockets in both targets. Structure-based virtual screens identified ten compounds, of which three had affinity for both targets. One of these scaffolds was optimized to nanomolar dual-target activity and showed the predicted pharmacodynamic effect in a rat model of Parkinsonism.

Accelerating the Throughput of Affinity Mass Spectrometry-Based Ligand Screening toward a G Protein-Coupled Receptor

Affinity mass spectrometry (MS) enables rapid screening of compound mixtures for ligands bound to a specific protein target, yet its current throughput is limited to individually assay pools of 400-2000 compounds. Typical affinity MS screens implemented in pharmaceutical industry laboratories identify putative ligands based on qualitative analysis of compound binding to the target whereas no quantitative information is acquired to discriminate high- and low-affinity ligands in the screening phase. Furthermore, these screens require purification of a stabilized form of the protein target, which poses a great challenge for membrane receptor targets. Here, we describe a new, potentially general affinity MS strategy that allows screening of 20,000 compounds in one pool for highly efficient ligand discovery toward a G protein-coupled receptor (GPCR) target. Quantitative measurement of compound binding to the receptor enables high-affinity ligand selection using both the purified receptor and receptor-embedded cell membranes. This high-throughput, label-free and quantitative affinity MS screen resulted in discovery of three new antagonists of the A_2A adenosine receptor.

Drug Repositioning Inferred from E2F1-Coregulator Interactions Studies for the Prevention and Treatment of Metastatic Cancers

Metastasis management remains a long-standing challenge. High abundance of E2F1 triggers tumor progression by developing protein-protein interactions (PPI) with coregulators that enhance its potential to activate a network of prometastatic transcriptional targets.

Methods: To identify E2F1-coregulators, we integrated high-throughput Co-immunoprecipitation (IP)/mass spectometry, GST-pull-down assays, and structure modeling. Potential inhibitors of PPI discovered were found by bioinformatics-based pharmacophore modeling, and transcriptome profiling was conducted to screen for coregulated downstream targets. Expression and target gene regulation was validated using qRT-PCR, immunoblotting, chromatin IP, and luciferase assays. Finally, the impact of the E2F1-coregulator complex and its inhibiting drug on metastasis was investigated in vitro in different cancer entities and two mouse metastasis models.

Results: We unveiled that E2F1 forms coactivator complexes with metastasis-associated protein 1 (MTA1) which, in turn, is directly upregulated by E2F1. The E2F1:MTA1 complex potentiates hyaluronan synthase 2 (HAS2) expression, increases hyaluronan production and promotes cell motility. Disruption of this prometastatic E2F1:MTA1 interaction reduces hyaluronan synthesis and infiltration of tumor-associated macrophages in the tumor microenvironment, thereby suppressing metastasis. We further demonstrate that E2F1:MTA1 assembly is abrogated by small-molecule, FDA-approved drugs. Treatment of E2F1/MTA1-positive, highly aggressive, circulating melanoma cells and orthotopic pancreatic tumors with argatroban prevents metastasis and cancer relapses in vivo through perturbation of the E2F1:MTA1/HAS2 axis.

Conclusion: Our results propose argatroban as an innovative, E2F-coregulator-based, antimetastatic drug. Cancer patients with the infaust E2F1/MTA1/HAS2 signature will likely benefit from drug repositioning.

A Structural Framework for GPCR Chemogenomics: What's In a Residue Number?

The recent surge of crystal structures of G protein-coupled receptors (GPCRs), as well as comprehensive collections of sequence, structural, ligand bioactivity, and mutation data, has enabled the development of integrated chemogenomics workflows for this important target family. This chapter will focus on cross-family and cross-class studies of GPCRs that have pinpointed the need for, and the implementation of, a generic numbering scheme for referring to specific structural elements of GPCRs. Sequence- and structure-based numbering schemes for different receptor classes will be introduced and the remaining caveats will be discussed. The use of these numbering schemes has facilitated many chemogenomics studies such as consensus binding site definition, binding site comparison, ligand repurposing (e.g. for orphan receptors), sequence-based pharmacophore generation for homology modeling or virtual screening, and class-wide chemogenomics studies of GPCRs.

Recent Advances and Applications of Molecular Docking to G Protein-Coupled Receptors

The growing number of studies on G protein-coupled receptors (GPCRs) family are a source of noticeable improvement in our understanding of the functioning of these proteins. GPCRs are responsible for a vast part of signaling in vertebrates and, as such, invariably remain in the spotlight of medicinal chemistry. A deeper insight into the underlying mechanisms of interesting phenomena observed in GPCRs, such as biased signaling or allosteric modulation, can be gained with experimental and computational studies. The latter play an important role in this process, since they allow for observations on scales inaccessible for most other methods. One of the key steps in such studies is proper computational reconstruction of actual ligand-receptor or protein-protein interactions, a process called molecular docking. A number of improvements and innovative applications of this method were documented recently. In this review, we focus particularly on innovations in docking to GPCRs.

In search of novel ligands using a structure-based approach: a case study on the adenosine A2A receptor

In this work, we present a case study to explore the challenges associated with finding novel molecules for a receptor that has been studied in depth and has a wealth of chemical information available. Specifically, we apply a previously described protocol that incorporates explicit water molecules in the ligand binding site to prospectively screen over 2.5 million drug-like and lead-like compounds from the commercially available eMolecules database in search of novel binders to the adenosine A_2A receptor (A_2AAR). A total of seventy-one compounds were selected for purchase and biochemical assaying based on high ligand efficiency and high novelty (Tanimoto coefficient ≤0.25 to any A_2AAR tested compound). These molecules were then tested for their affinity to the adenosine A_2A receptor in a radioligand binding assay. We identified two hits that fulfilled the criterion of ~50 % radioligand displacement at a concentration of 10 μM. Next we selected an additional eight novel molecules that were predicted to make a bidentate interaction with Asn253^6.55, a key interacting residue in the binding pocket of the A_2AAR. None of these eight molecules were found to be active. Based on these results we discuss the advantages of structure-based methods and the challenges associated with finding chemically novel molecules for well-explored targets.

Function-specific virtual screening for GPCR ligands using a combined scoring method

The ability of scoring functions to correctly select and rank docking poses of small molecules in protein binding sites is highly target dependent, which presents a challenge for structure-based drug discovery. Here we describe a virtual screening method that combines an energy-based docking scoring function with a molecular interaction fingerprint (IFP) to identify new ligands based on G protein-coupled receptor (GPCR) crystal structures. The consensus scoring method is prospectively evaluated by: 1) the discovery of chemically novel, fragment-like, high affinity histamine H₁ receptor (H₁R) antagonists/inverse agonists, 2) the selective structure-based identification of ß₂-adrenoceptor (ß₂R) agonists, and 3) the experimental validation and comparison of the combined and individual scoring approaches. Systematic retrospective virtual screening simulations allowed the definition of scoring cut-offs for the identification of H₁R and ß₂R ligands and the selection of an optimal ß-adrenoceptor crystal structure for the discrimination between ß₂R agonists and antagonists. The consensus approach resulted in the experimental validation of 53% of the ß₂R and 73% of the H₁R virtual screening hits with up to nanomolar affinities and potencies. The selective identification of ß₂R agonists shows the possibilities of structure-based prediction of GPCR ligand function by integrating protein-ligand binding mode information.

Structurally Enabled Discovery of Adenosine A2A Receptor Antagonists

Over the past decade there has been a revolution in the field of G protein-coupled receptor (GPCR) structural biology. Many years of innovative research from different areas have come together to fuel this significant change in the fortunes of this field, which for many years was characterized by the paucity of high-resolution structures. The determination to succeed has been in part due to the recognized importance of these proteins as drug targets, and although the pharmaceutical industry has been focusing on these receptors, it can be justifiably argued and demonstrated that many of the approved and commercially successful GPCR drugs can be significantly improved to increase efficacy and/or reduce undesired side effects. In addition, many validated targets in this class remain to be drugged. It is widely recognized that application of structure-based drug design approaches can help medicinal chemists a long way toward discovering better drugs. The achievement of structural biologists in providing high-resolution insight is beginning to transform drug discovery efforts, and there are a number of GPCR drugs that have been discovered by use of structural information that are in clinical development. This review aims to highlight the key developments that have brought success to GPCR structure resolution efforts and exemplify the practical application of structural information for the discovery of adenosine A_2A receptor antagonists that have potential to treat multiple conditions.

Hierarchical virtual screening approaches in small molecule drug discovery

Virtual screening has played a significant role in the discovery of small molecule inhibitors of therapeutic targets in last two decades. Various ligand and structure-based virtual screening approaches are employed to identify small molecule ligands for proteins of interest. These approaches are often combined in either hierarchical or parallel manner to take advantage of the strength and avoid the limitations associated with individual methods. Hierarchical combination of ligand and structure-based virtual screening approaches has received noteworthy success in numerous drug discovery campaigns. In hierarchical virtual screening, several filters using ligand and structure-based approaches are sequentially applied to reduce a large screening library to a number small enough for experimental testing. In this review, we focus on different hierarchical virtual screening strategies and their application in the discovery of small molecule modulators of important drug targets. Several virtual screening studies are discussed to demonstrate the successful application of hierarchical virtual screening in small molecule drug discovery.

A structural chemogenomics analysis of aminergic GPCRs: lessons for histamine receptor ligand design

BACKGROUND AND PURPOSE

Chemogenomics focuses on the discovery of new connections between chemical and biological space leading to the discovery of new protein targets and biologically active molecules. G-protein coupled receptors (GPCRs) are a particularly interesting protein family for chemogenomics studies because there is an overwhelming amount of ligand binding affinity data available. The increasing number of aminergic GPCR crystal structures now for the first time allows the integration of chemogenomics studies with high-resolution structural analyses of GPCR-ligand complexes.

EXPERIMENTAL APPROACH

In this study, we have combined ligand affinity data, receptor mutagenesis studies, and amino acid sequence analyses to high-resolution structural analyses of (hist)aminergic GPCR-ligand interactions. This integrated structural chemogenomics analysis is used to more accurately describe the molecular and structural determinants of ligand affinity and selectivity in different key binding regions of the crystallized aminergic GPCRs, and histamine receptors in particular.

KEY RESULTS

Our investigations highlight interesting correlations and differences between ligand similarity and ligand binding site similarity of different aminergic receptors. Apparent discrepancies can be explained by combining detailed analysis of crystallized or predicted protein-ligand binding modes, receptor mutation studies, and ligand structure-selectivity relationships that identify local differences in essential pharmacophore features in the ligand binding sites of different receptors.

CONCLUSIONS AND IMPLICATIONS

We have performed structural chemogenomics studies that identify links between (hist)aminergic receptor ligands and their binding sites and binding modes. This knowledge can be used to identify structure-selectivity relationships that increase our understanding of ligand binding to (hist)aminergic receptors and hence can be used in future GPCR ligand discovery and design.

Structure-based drug design for G protein-coupled receptors

Our understanding of the structural biology of G protein-coupled receptors has undergone a transformation over the past 5 years. New protein-ligand complexes are described almost monthly in high profile journals. Appreciation of how small molecules and natural ligands bind to their receptors has the potential to impact enormously how medicinal chemists approach this major class of receptor targets. An outline of the key topics in this field and some recent examples of structure- and fragment-based drug design are described. A table is presented with example views of each G protein-coupled receptor for which there is a published X-ray structure, including interactions with small molecule antagonists, partial and full agonists. The possible implications of these new data for drug design are discussed.

Homology modeling of human muscarinic acetylcholine receptors

We have developed homology models of the acetylcholine muscarinic receptors M₁R-M₅R, based on the β₂-adrenergic receptor crystal as the template. This is the first report of homology modeling of all five subtypes of acetylcholine muscarinic receptors with binding sites optimized for ligand binding. The models were evaluated for their ability to discriminate between muscarinic antagonists and decoy compounds using virtual screening using enrichment factors, area under the ROC curve (AUC), and an early enrichment measure, LogAUC. The models produce rational binding modes of docked ligands as well as good enrichment capacity when tested against property-matched decoy libraries, which demonstrates their unbiased predictive ability. To test the relative effects of homology model template selection and the binding site optimization procedure, we generated and evaluated a naïve M₂R model, using the M₃R crystal structure as a template. Our results confirm previous findings that binding site optimization using ligand(s) active at a particular receptor, i.e. including functional knowledge into the model building process, has a more pronounced effect on model quality than target-template sequence similarity. The optimized M₁R-M₅R homology models are made available as part of the Supporting Information to allow researchers to use these structures, compare them to their own results, and thus advance the development of better modeling approaches.

Current progress in Structure-Based Rational Drug Design marks a new mindset in drug discovery

The past decade has witnessed a paradigm shift in preclinical drug discovery with structure-based drug design (SBDD) making a comeback while high-throughput screening (HTS) methods have continued to generate disappointing results. There is a deficit of information between identified hits and the many criteria that must be fulfilled in parallel to convert them into preclinical candidates that have a real chance to become a drug. This gap can be bridged by investigating the interactions between the ligands and their receptors. Accurate calculations of the free energy of binding are still elusive; however progresses were made with respect to how one may deal with the versatile role of water. A corpus of knowledge combining X-ray structures, bioinformatics and molecular modeling techniques now allows drug designers to routinely produce receptor homology models of increasing quality. These models serve as a basis to establish and validate efficient rationales used to tailor and/or screen virtual libraries with enhanced chances of obtaining hits. Many case reports of successful SBDD show how synergy can be gained from the combined use of several techniques. The role of SBDD with respect to two different classes of widely investigated pharmaceutical targets: (a) protein kinases (PK) and (b) G-protein coupled receptors (GPCR) is discussed. Throughout these examples prototypical situations covering the current possibilities and limitations of SBDD are presented.

X-ray structural information of GPCRs in drug design: what are the limitations and where do we go?

INTRODUCTION

In 2007, the X-ray structural determination of non-rhodopsin G-Protein coupled receptors (GPCRs), considered the most extensively targeted protein class for marketed drugs, commenced. With the relatively rapid availability of additional structures, an assessment of the progression made is needed in addition to the assessment of the understandings gleaned, deployment successes and forthcoming prospects.

AREAS COVERED

The author reviews the approaches and tools that have made it possible to determine the three dimensional structures of GPCRs using X-ray crystallography. Furthermore, the author describes the methods suited for crystallization of membrane bound GPCR proteins including the lipidic cubic phase and various protein modification approaches. The author also provides highlights, from the literature, of the structures determined to date including targets solved, the nature of the content provided (such as selectivity, activating vs. inactivating determinants) and how these structural features relate to drug design strategies.

EXPERT OPINION

The GPCR X-ray structures that have been so far determined have yielded significant information. This has presented dramatic evidence concerning their ability to impact the discovery of compounds through their action as traditional, orthosteric modulators. It is, however, noted that more challenging design strategies, such as identifying biased agonists and the use of sites remote from the orthosteric site for allosteric modulation, are still in their infancy.

Chemogenomics of allosteric binding sites in GPCRs

Chemogenomic techniques connect the chemical and biological domains to establish ligand and target relationships not evident from the individual disciplines. Chemogenomics has been applied in lead generation, target classification, focused library design as well as selectivity and polypharmacology profiling. This review describes recent developments structured into ligand-, target- and combined chemogenomic techniques and applications to allosteric GPCR ligands. It also outlines relative strengths and limitations of these techniques and the impact of the increasing crystallographic data.

Small and colorful stones make beautiful mosaics: fragment-based chemogenomics

Smaller stones with a wide variety of colors make a higher resolution mosaic. In much the same way, smaller chemical entities that are structurally diverse are better able to interrogate protein binding sites. This feature article describes the construction of a diverse fragment library and an analysis of the screening of six representative protein targets belonging to three diverse target classes (G protein-coupled receptors ADRB2, H1R, H3R, and H4R, the ligand-gated ion channel 5-HT3R, and the kinase PKA) using chemogenomics approaches. The integration of experimentally determined bioaffinity profiles across related and unrelated protein targets and chemogenomics analysis of fragment binding and protein structure allow the identification of: (i) unexpected similarities and differences in ligand binding properties, and (ii) subtle ligand affinity and selectivity cliffs. With a wealth of fragment screening data being generated in industry and academia, such approaches will contribute to a more detailed structural understanding of ligand-protein interactions.