An automated system for the analysis of G protein-coupled receptor transmembrane binding pockets: alignment, receptor-based pharmacophores, and their application

G protein-coupled receptors (GPCRs) share a common architecture consisting of seven transmembrane (TM) domains. Various lines of evidence suggest that this fold provides a generic binding pocket within the TM region for hosting agonists, antagonists, and allosteric modulators. Here, a comprehensive and automated method allowing fast analysis and comparison of these putative binding pockets across the entire GPCR family is presented. The method relies on a robust alignment algorithm based on conservation indices, focusing on pharmacophore-like relationships between amino acids. Analysis of conservation patterns across the GPCR family and alignment to the rhodopsin X-ray structure allows the extraction of the amino acids lining the TM binding pocket in a so-called ligand binding pocket vector (LPV). In a second step, LPVs are translated to simple 3D receptor pharmacophore models, where each amino acid is represented by a single spherical pharmacophore feature and all atomic detail is omitted. Applications of the method include the assessment of selectivity issues, support of mutagenesis studies, and the derivation of rules for focused screening to identify chemical starting points in early drug discovery projects. Because of the coarseness of this 3D receptor pharmacophore model, however, meaningful scoring and ranking procedures of large sets of molecules are not justified. The LPV analysis of the trace amine-associated receptor family and its experimental validation is discussed as an example. The value of the 3D receptor model is demonstrated for a class C GPCR family, the metabotropic glutamate receptors.     

Development of selective high affinity antagonists, agonists, and radioligands for the P2Y1 receptor

The P2Y(1) receptor is a member of the P2Y family of nucleotide-activated G protein-coupled receptors, and it is an important therapeutic target based on its broad tissue distribution and essential role in platelet aggregation. We have designed a set of highly selective and diverse pharmacological tools for studying the P2Y(1) receptor using a rational approach to ligand design. Based on the discovery that bisphosphate analogues of the P2Y(1) receptor agonist, ADP, are partial agonists/competitive antagonists of this receptor, an iterative approach was used to develop competitive antagonists with enhanced affinity and selectivity. Halogen substitutions of the 2-position of the adenine ring provided increased affinity while an N(6) methyl substitution eliminated partial agonist activity. Furthermore, various replacements of the ribose ring with symmetrically branched, phosphorylated acyclic structures revealed that the ribose is not necessary for recognition at the P2Y(1) receptor. Finally, replacement of the ribose ring with a five member methanocarba ring constrained in the Northern conformation conferred dramatic increases in affinity to both P2Y(1) receptor antagonists as well as agonists. These combined structural modifications have resulted in a series of selective high affinity antagonists of the P2Y(1) receptor, two broadly applicable radioligands, and a high affinity agonist capable of selectively activating the P2Y(1) receptor in human platelets. Complementary receptor modeling and computational ligand docking have provided a putative structural framework for the drug-receptor interactions. A similar rational approach is being applied to develop selective ligands for other subtypes of P2Y receptors.     

High content screening for G protein-coupled receptors using cell-based protein translocation assays

G protein-coupled receptors (GPCRs) have been one of the most productive classes of drug targets for several decades, and new technologies for GPCR-based discovery promise to keep this field active for years to come. While molecular screens for GPCR receptor agonist- and antagonist-based drugs will continue to be valuable discovery tools, the most exciting developments in the field involve cell-based assays for GPCR function. Some cell-based discovery strategies, such as the use of beta-arrestin as a surrogate marker for GPCR function, have already been reduced to practice, and have been used as valuable discovery tools for several years. The application of high content cell-based screening to GPCR discovery has opened up additional possibilities, such as direct tracking of GPCRs, G proteins and other signaling pathway components using intracellular translocation assays. These assays provide the capability to probe GPCR function at the cellular level with better resolution than has previously been possible, and offer practical strategies for more definitive selectivity evaluation and counter-screening in the early stages of drug discovery. The potential of cell-based translocation assays for GPCR discovery is described, and proof-of-concept data from a pilot screen with a CXCR4 assay are presented. This chemokine receptor is a highly relevant drug target which plays an important role in the pathogenesis of inflammatory disease and also has been shown to be a co-receptor for entry of HIV into cells as well as to play a role in metastasis of certain cancer cells.     

Toward deciphering the code to aminergic G protein-coupled receptor drug design

The trace amine-associated receptor 1 (TAAR(1)) is a biogenic amine G protein-coupled receptor (GPCR) that is potently activated by 3-iodothyronamine (1, T(1)AM) in vitro. Compound 1 is an endogenous derivative of the thyroid hormone thyroxine which rapidly induces hypothermia, anergia, and bradycardia when administered to mice. To explore the role of TAAR(1) in mediating the effects of 1, we rationally designed and synthesized rat TAAR(1) superagonists and lead antagonists using the rotamer toggle switch model of aminergic GPCR activation. The functional activity of a ligand is proposed to be correlated to its probable interactions with the rotamer switch residues; agonists allow the rotamer switch residues to toggle to their active conformation, whereas antagonists interfere with this conformational transition. These agonist and antagonist design principles provide a conceptual model for understanding the relationship between the molecular structure of a drug and its pharmacological properties.     

A CE assay for the detection of agonist-stimulated adenylyl cyclase activity

A CE assay was developed for the detection of adenylyl cyclase (AC) activity stimulated at the AC and G protein-coupled receptor (GPCR) level. In the assay, cell membranes overexpressing GPCR and/or AC were incubated with modulators and substrate ATP to produce cAMP in a dose-dependent manner. In both the CE-UV and a radiochemical assay, the addition of forskolin (FSK) resulted in a two- to three-fold maximum increase in AC activity with EC50s of 4.2 +/- 0.7 and 2.4 +/- 0.7 microM, respectively, demonstrating that similar results were obtained by both assays. GPCR activation was also detected using cell membranes overexpressing AC and the beta2-adrenergic receptor (beta2AR) fused to the stimulatory G protein. Terbutaline (beta2AR agonist) increased the basal rate of cAMP formation 1.7 +/- 0.1-fold resulting in an EC50 of 62 +/- 10 nM. The assay's ability to detect antagonists is demonstrated by the expected right-shifted EC50 of terbutaline by the beta2AR antagonist propranolol. The CE-UV assay offers advantages over the traditional radioactivity assay in terms of safety and labor.     

Mechanistic Studies on the Stereoselectivity of the Serotonin 5‐HT1A Receptor

G‐protein‐coupled receptors (GPCRs) are involved in a wide range of physiological processes, and they have attracted considerable attention as important targets for developing new medicines. A central and largely unresolved question in drug discovery, which is especially relevant to GPCRs, concerns ligand selectivity: Why do certain molecules act as activators (agonists) whereas others, with nearly identical structures, act as blockers (antagonists) of GPCRs? To address this question, we employed all‐atom, long‐timescale molecular dynamics simulations to investigate how two diastereomers (epimers) of dihydrofuroaporphine bind to the serotonin 5‐HT_1A receptor and exert opposite effects. By using molecular interaction fingerprints, we discovered that the agonist could mobilize nearby amino acid residues to act as molecular switches for the formation of a continuous water channel. In contrast, the antagonist epimer remained firmly stabilized in the binding pocket.     

Biophysical approaches to G protein-coupled receptors: Structure, function and dynamics

G protein-coupled receptors (GPCR) represent a large family of drug targets for which there is no high-resolution structural information. In order to understand the mechanisms of ligand recognition and receptor activation, there is a strong need for novel biophysical methods. In this Perspective we provide an overview of recent experimental approaches used to explore the molecular architecture and dynamics of GPCR and their interactions with ligands and G proteins using biophysical, non-crystallographic, methods.     

Ligand and decoy sets for docking to G protein-coupled receptors

We compiled a G protein-coupled receptor (GPCR) ligand library (GLL) for 147 targets, selecting for each ligand 39 decoy molecules, collected in the GPCR Decoy Database (GDD). Decoys were chosen ensuring a ligand-decoy similarity of six physical properties, while enforcing ligand-decoy chemical dissimilarity. The performance in docking of the GDD was evaluated on 19 GPCRs, showing a marked decrease in enrichment compared to bias-uncorrected decoy sets. Both the GLL and GDD are freely available for the scientific community.     

Matching cavities in g protein-coupled receptors to infer ligand-binding sites

To understand the activity and cross reactivity of ligands and G protein-coupled receptors, we take stock of relevant existing receptor mutation, sequence, and structural data to develop a statistically robust and transparent scoring system. Our method evaluates the viability of binding of any ligand for any GPCR sequence of amino acids. This enabled us to explore the binding repertoire of both receptors and ligands, relying solely on correlations between carefully identified receptor features and without requiring any chemical information about ligands. This study suggests that sequence similarity at specific binding pockets can predict relative affinity of ligands; enabling recovery of over 80% of known ligands for a withheld receptor and almost 80% of known receptors for a ligand. The method enables qualitative prediction of ligand binding for all nonredundant human G protein-coupled receptors.     

Small Molecule Tools to Study Cellular Target Engagement for the Intracellular Allosteric Binding Site of GPCRs

A conserved intracellular allosteric binding site (IABS) has recently been identified at several G protein-coupled receptors (GPCRs). Ligands targeting the IABS, so-called intracellular allosteric antagonists, are highly promising compounds for pharmaceutical intervention and currently evaluated in several clinical trials. Beside co-crystal structures that laid the foundation for the structure-based development of intracellular allosteric GPCR antagonists, small molecule tools that enable an unambiguous identification and characterization of intracellular allosteric GPCR ligands are of utmost importance for drug discovery campaigns in this field. Herein, we discuss recent approaches that leverage cellular target engagement studies for the IABS and thus play a critical role in the evaluation of IABS-targeted ligands as potential therapeutic agents.     

A virtual screen for diverse ligands: discovery of selective G protein-coupled receptor antagonists

Virtual screening has become a major focus of bioactive small molecule lead identification, and reports of agonists and antagonists discovered via virtual methods are becoming more frequent. G protein-coupled receptors (GPCRs) are the one class of protein targets for which success with this approach has been limited. This is likely due to the paucity of detailed experimental information describing GPCR structure and the intrinsic function-associated structural flexibility of GPCRs which present major challenges in the application of receptor-based virtual screening. Here we describe an in silico methodology that diminishes the effects of structural uncertainty, allowing for more inclusive representation of a potential docking interaction with exogenous ligands. Using this approach, we screened one million compounds from a virtual database, and a diverse subgroup of 100 compounds was selected, leading to experimental identification of five structurally diverse antagonists of the thyrotropin-releasing hormone receptors (TRH-R1 and TRH-R2). The chirality of the most potent chemotype was demonstrated to be important in its binding affinity to TRH receptors; the most potent stereoisomer was noted to have a 13-fold selectivity for TRH-R1 over TRH-R2. A comprehensive mutational analysis of key amino acid residues that form the putative binding pocket of TRH receptors further verified the binding modality of these small molecule antagonists. The described virtual screening approach may prove applicable in the search for novel small molecule agonists and antagonists of other GPCRs.     

Atomistic insights into rhodopsin activation from a dynamic model

Rhodopsin, the light sensitive receptor responsible for blue-green vision, serves as a prototypical G protein-coupled receptor (GPCR). Upon light absorption, it undergoes a series of conformational changes that lead to the active form, metarhodopsin II (META II), initiating a signaling cascade through binding to the G protein transducin (G(t)). Here, we first develop a structural model of META II by applying experimental distance restraints to the structure of lumi-rhodopsin (LUMI), an earlier intermediate. The restraints are imposed by using a combination of biased molecular dynamics simulations and perturbations to an elastic network model. We characterize the motions of the transmembrane helices in the LUMI-to-META II transition and the rearrangement of interhelical hydrogen bonds. We then simulate rhodopsin activation in a dynamic model to study the path leading from LUMI to our META II model for wild-type rhodopsin and a series of mutants. The simulations show a strong correlation between the transition dynamics and the pharmacological phenotypes of the mutants. These results help identify the molecular mechanisms of activation in both wild type and mutant rhodopsin. While static models can provide insights into the mechanisms of ligand recognition and predict ligand affinity, a dynamic model of activation could be applicable to study the pharmacology of other GPCRs and their ligands, offering a key to predictions of basal activity and ligand efficacy.     

Design of Drug Efficacy Guided by Free Energy Simulations of the β2-Adrenoceptor

G-protein-coupled receptors (GPCRs) play important roles in physiological processes and are modulated by drugs that either activate or block signaling. Rational design of the pharmacological efficacy profiles of GPCR ligands could enable the development of more efficient drugs, but is challenging even if high-resolution receptor structures are available. We performed molecular dynamics simulations of the β₂ adrenergic receptor in active and inactive conformations to assess if binding free energy calculations can predict differences in ligand efficacy for closely related compounds. Previously identified ligands were successfully classified into groups with comparable efficacy profiles based on the calculated shift in ligand affinity upon activation. A series of ligands were then predicted and synthesized, leading to the discovery of partial agonists with nanomolar potencies and novel scaffolds. Our results demonstrate that free energy simulations enable design of ligand efficacy and the same approach can be applied to other GPCR drug targets.     

INPHARMA-Based Determination of Ligand Binding Modes at α1-Adrenergic Receptors Explains the Molecular Basis of Subtype Selectivity

The structural poses of ligands that bind weakly to protein receptors are challenging to define. In this work we have studied ligand interactions with the adrenoreceptor (AR) subtypes, α_1A-AR and α_1B-AR, which belong to the G protein-coupled receptor (GPCR) superfamily, by employing the solution-based ligand-observed NMR method interligand NOEs for pharmacophore mapping (INPHARMA). A lack of receptor crystal structures and of subtype-selective drugs has hindered the definition of the physiological roles of each subtype and limited drug development. We determined the binding pose of the weakly binding α_1A-AR-selective agonist A-61603 relative to an endogenous agonist, epinephrine, at both α_1A-AR and α_1B-AR. The NMR experimental data were quantitatively compared, by using SpINPHARMA, to the back-calculated spectra based on ligand poses obtained from all-atom molecular dynamics simulations. The results helped mechanistically explain the selectivity of (R)-A-61603 towards α_1A-AR, thus demonstrating an approach for targeting subtype selectivity in ARs.     

Ligand and structure-based methodologies for the prediction of the activity of G protein-coupled receptor ligands

Accurate in silico models for the quantitative prediction of the activity of G protein-coupled receptor (GPCR) ligands would greatly facilitate the process of drug discovery and development. Several methodologies have been developed based on the properties of the ligands, the direct study of the receptor-ligand interactions, or a combination of both approaches. Ligand-based three-dimensional quantitative structure-activity relationships (3D-QSAR) techniques, not requiring knowledge of the receptor structure, have been historically the first to be applied to the prediction of the activity of GPCR ligands. They are generally endowed with robustness and good ranking ability; however they are highly dependent on training sets. Structure-based techniques generally do not provide the level of accuracy necessary to yield meaningful rankings when applied to GPCR homology models. However, they are essentially independent from training sets and have a sufficient level of accuracy to allow an effective discrimination between binders and nonbinders, thus qualifying as viable lead discovery tools. The combination of ligand and structure-based methodologies in the form of receptor-based 3D-QSAR and ligand and structure-based consensus models results in robust and accurate quantitative predictions. The contribution of the structure-based component to these combined approaches is expected to become more substantial and effective in the future, as more sophisticated scoring functions are developed and more detailed structural information on GPCRs is gathered.     

A ligand-based approach to mining the chemogenomic space of drugs

The practical implementation and validation of a ligand-based approach to mining the chemogenomic space of drugs is presented and applied to the in silico target profiling of 767 drugs against 684 targets of therapeutic relevance. The results reveal that drugs targeting aminergic G protein-coupled receptors (GPCRs) show the most promiscuous pharmacological profiles. The detection of cross-pharmacologies between aminergic GPCRs and the opioid, sigma, NMDA, and 5-HT3 receptors aggravate the potential promiscuity of those drugs, predominantly including analgesics, antidepressants, and antipsychotics.     

G protein- and agonist-bound serotonin 5-HT2A receptor model activated by steered molecular dynamics simulations

A 5-HT(2A) receptor model was constructed by homology modeling based on the β(2)-adrenergic receptor and the G protein-bound opsin crystal structures. The 5-HT(2A) receptor model was transferred into an active conformation by an agonist ligand and a G(αq) peptide in four subsequent steered molecular dynamics (MD) simulations. The driving force for the transformation was the addition of several known intermolecular and receptor interhelical hydrogen bonds enforcing the necessary helical and rotameric movements. Subsquent MD simulations without constraints confirmed the stability of the activated receptor model as well as revealed new information about stabilizing residues and bonds. The active 5-HT(2A) receptor model was further validated by retrospective ligand screening of more than 9400 compounds, whereof 182 were known ligands. The results show that the model can be used in drug discovery for virtual screening and structure-based ligand design as well as in GPCR activation studies.     

Label-free cell-based assays for GPCR screening

G protein-coupled receptors (GPCRs) have been proven to be the largest family of druggable targets in the human genome. Given the importance of GPCRs as drug targets and the de-orphanization of novel targets, GPCRs are likely to remain the frequent targets of many drug discovery programs. With recent advances in instrumentation and understanding of cellular mechanisms for the signals measured, biosensor-centered label-free cell assay technologies become a very active area for GPCR screening. This article reviews the principles and potential of current label-free cell assay technologies in GPCR drug discovery.     

Monitoring the Reversibility of GPCR Signaling by Combining Photochromic Ligands with Label-free Impedance Analysis

G protein-coupled cell surface receptors (GPCR) trigger complex intracellular signaling cascades upon agonist binding. Classic pharmacological assays provide information about binding affinities, activation or blockade at different stages of the signaling cascade, but real time dynamics and reversibility of these processes remain often disguised. We show that combining photochromic NPY receptor ligands, which can be toggled in their receptor activation ability by irradiation with light of different wavelengths, with whole cell label-free impedance assays allows observing the cell response to receptor activation and its reversibility over time. The concept demonstrated on NPY receptors may be well applicable to many other GPCRs providing a deeper insight into the time course of intracellular signaling processes.