Discovery of Human Signaling Systems: Pairing Peptides to G Protein-Coupled Receptors.

The peptidergic system is the most abundant network of ligand-receptor-mediated signaling in humans. However, the physiological roles remain elusive for numerous peptides and more than 100 G protein-coupled receptors (GPCRs). Here we report the pairing of cognate peptides and receptors. Integrating comparative genomics across 313 species and bioinformatics on all protein sequences and structures of human class A GPCRs, we identify universal characteristics that uncover additional potential peptidergic signaling systems. Using three orthogonal biochemical assays, we pair 17 proposed endogenous ligands with five orphan GPCRs that are associated with diseases, including genetic, neoplastic, nervous and reproductive system disorders. We also identify additional peptides for nine receptors with recognized ligands and pathophysiological roles. This integrated computational and multifaceted experimental approach expands the peptide-GPCR network and opens the way for studies to elucidate the roles of these signaling systems in human physiology and disease. VIDEO ABSTRACT.

Pharmacogenomics of GPCR Drug Targets.

Natural genetic variation in the human genome is a cause of individual differences in responses to medications and is an underappreciated burden on public health. Although 108 G-protein-coupled receptors (GPCRs) are the targets of 475 (∼34%) Food and Drug Administration (FDA)-approved drugs and account for a global sales volume of over 180 billion US dollars annually, the prevalence of genetic variation among GPCRs targeted by drugs is unknown. By analyzing data from 68,496 individuals, we find that GPCRs targeted by drugs show genetic variation within functional regions such as drug- and effector-binding sites in the human population. We experimentally show that certain variants of µ-opioid and Cholecystokinin-A receptors could lead to altered or adverse drug response. By analyzing UK National Health Service drug prescription and sales data, we suggest that characterizing GPCR variants could increase prescription precision, improving patients' quality of life, and relieve the economic and societal burden due to variable drug responsiveness. VIDEO ABSTRACT.

G protein-coupled receptor (GPCR) gene variants and human genetic disease.

Genetic variations in the genes encoding G protein-coupled receptors (GPCRs) can disrupt receptor structure and function, which can result in human genetic diseases. Disease-causing mutations have been reported in at least 55 GPCRs for more than 66 monogenic diseases in humans. The spectrum of pathogenic and likely pathogenic variants includes loss of function variants that decrease receptor signaling on one extreme and gain of function that may result in biased signaling or constitutive activity, originally modeled on prototypical rhodopsin GPCR variants identified in retinitis pigmentosa, on the other. GPCR variants disrupt ligand binding, G protein coupling, accessory protein function, receptor desensitization and receptor recycling. Next generation sequencing has made it possible to identify variants of uncertain significance (VUS). We discuss variants in receptors known to result in disease and in silico strategies for disambiguation of VUS such as sorting intolerant from tolerant and polymorphism phenotyping. Modeling of variants has contributed to drug development and precision medicine, including drugs that target the melanocortin receptor in obesity and interventions that reverse loss of gonadotropin-releasing hormone receptor from the cell surface in idiopathic hypogonadotropic hypogonadism. Activating and inactivating variants of the calcium sensing receptor (CaSR) gene that are pathogenic in familial hypocalciuric hypercalcemia and autosomal dominant hypocalcemia have enabled the development of calcimimetics and calcilytics. Next generation sequencing has continued to identify variants in GPCR genes, including orphan receptors, that contribute to human phenotypes and may have therapeutic potential. Variants of the CaSR gene, some encoding an arginine-rich region that promotes receptor phosphorylation and intracellular retention, have been linked to an idiopathic epilepsy syndrome. Agnostic strategies have identified variants of the pyroglutamylated RF amide peptide receptor gene in intellectual disability and G protein-coupled receptor 39 identified in psoriatic arthropathy. Coding variants of the G protein-coupled receptor L1 (GPR37L1) orphan receptor gene have been identified in a rare familial progressive myoclonus epilepsy. The study of the role of GPCR variants in monogenic, Mendelian phenotypes has provided the basis of modeling the significance of more common variants of pharmacogenetic significance.

The G protein database, GproteinDb.

Two-thirds of signaling substances, several sensory stimuli and over one-third of drugs act via receptors coupling to G proteins. Here, we present an online platform for G protein research with reference data and tools for analysis, visualization and design of scientific studies across disciplines and areas. This platform may help translate new pharmacological, structural and genomic data into insights on G protein signaling vital for human physiology and medicine. The G protein database is accessible at https://gproteindb.org.

Molecular and in vivo phenotyping of missense variants of the human glucagon receptor.

Naturally occurring missense variants of G protein-coupled receptors with loss of function have been linked to metabolic disease in case studies and in animal experiments. The glucagon receptor, one such G protein-coupled receptor, is involved in maintaining blood glucose and amino acid homeostasis; however, loss-of-function mutations of this receptor have not been systematically characterized. Here, we observed fewer glucagon receptor missense variants than expected, as well as lower allele diversity and fewer variants with trait associations as compared with other class B1 receptors. We performed molecular pharmacological phenotyping of 38 missense variants located in the receptor extracellular domain, at the glucagon interface, or with previously suggested clinical implications. These variants were characterized in terms of cAMP accumulation to assess glucagon-induced Gαs coupling, and of recruitment of ß-arrestin-1/2. Fifteen variants were impaired in at least one of these downstream functions, with six variants affected in both cAMP accumulation and ß-arrestin-1/2 recruitment. For the eight variants with decreased Gαs signaling (D63ECDN, P86ECDS, V96ECDE, G125ECDC, R2253.30H, R3085.40W, V3686.59M, and R3787.35C) binding experiments revealed preserved glucagon affinity, although with significantly reduced binding capacity. Finally, using the UK Biobank, we found that variants with wildtype-like Gαs signaling did not associate with metabolic phenotypes, whereas carriers of cAMP accumulation-impairing variants displayed a tendency toward increased risk of obesity and increased body mass and blood pressure. These observations are in line with the essential role of the glucagon system in metabolism and support that Gαs is the main signaling pathway effecting the physiological roles of the glucagon receptor.

Selectivity determinants of GPCR-G-protein binding.

The selective coupling of G-protein-coupled receptors (GPCRs) to specific G proteins is critical to trigger the appropriate physiological response. However, the determinants of selective binding have remained elusive. Here we reveal the existence of a selectivity barcode (that is, patterns of amino acids) on each of the 16 human G proteins that is recognized by distinct regions on the approximately 800 human receptors. Although universally conserved positions in the barcode allow the receptors to bind and activate G proteins in a similar manner, different receptors recognize the unique positions of the G-protein barcode through distinct residues, like multiple keys (receptors) opening the same lock (G protein) using non-identical cuts. Considering the evolutionary history of GPCRs allows the identification of these selectivity-determining residues. These findings lay the foundation for understanding the molecular basis of coupling selectivity within individual receptors and G proteins.

Receptor selectivity between the G proteins Gα12 and Gα13 is defined by a single leucine-to-isoleucine variation.

Despite recent advances in structural definition of GPCR-G protein complexes, the basis of receptor selectivity between G proteins remains unclear. The Gα12 and Gα13 subtypes together form the least studied group of heterotrimeric G proteins. G protein-coupled receptor 35 (GPR35) has been suggested to couple efficiently to Gα13 but weakly to Gα12. Using combinations of cells genome-edited to not express G proteins and bioluminescence resonance energy transfer-based sensors, we confirmed marked selectivity of GPR35 for Gα13. Incorporating Gα12/Gα13 chimeras and individual residue swap mutations into these sensors defined that selectivity between Gα13 and Gα12 was imbued largely by a single leucine-to-isoleucine variation at position G.H5.23. Indeed, leucine could not be substituted by other amino acids in Gα13 without almost complete loss of GPR35 coupling. The critical importance of leucine at G.H5.23 for GPR35-G protein interaction was further demonstrated by introduction of this leucine into Gαq, resulting in the gain of coupling to GPR35. These studies demonstrate that Gα13 is markedly the most effective G protein for interaction with GPR35 and that selection between Gα13 and Gα12 is dictated largely by a single conservative amino acid variation.-Mackenzie, A. E., Quon, T., Lin, L.-C., Hauser, A. S., Jenkins, L., Inoue, A., Tobin, A. B., Gloriam, D. E., Hudson, B. D., Milligan, G. Receptor selectivity between the G proteins Gα12 and Gα13 is defined by a single leucine-to-isoleucine variation.

GPCRdb in 2018: adding GPCR structure models and ligands.

G protein-coupled receptors are the most abundant mediators of both human signalling processes and therapeutic effects. Herein, we report GPCRome-wide homology models of unprecedented quality, and roughly 150 000 GPCR ligands with data on biological activities and commercial availability. Based on the strategy of 'Less model - more Xtal', each model exploits both a main template and alternative local templates. This achieved higher similarity to new structures than any of the existing resources, and refined crystal structures with missing or distorted regions. Models are provided for inactive, intermediate and active states-except for classes C and F that so far only have inactive templates. The ligand database has separate browsers for: (i) target selection by receptor, family or class, (ii) ligand filtering based on cross-experiment activities (min, max and mean) or chemical properties, (iii) ligand source data and (iv) commercial availability. SMILES structures and activity spreadsheets can be downloaded for further processing. Furthermore, three recent landmark publications on GPCR drugs, G protein selectivity and genetic variants have been accompanied with resources that now let readers view and analyse the findings themselves in GPCRdb. Altogether, this update will enable scientific investigation for the wider GPCR community. GPCRdb is available at http://www.gpcrdb.org.

GPCRdb: an information system for G protein-coupled receptors.

Recent developments in G protein-coupled receptor (GPCR) structural biology and pharmacology have greatly enhanced our knowledge of receptor structure-function relations, and have helped improve the scientific foundation for drug design studies. The GPCR database, GPCRdb, serves a dual role in disseminating and enabling new scientific developments by providing reference data, analysis tools and interactive diagrams. This paper highlights new features in the fifth major GPCRdb release: (i) GPCR crystal structure browsing, superposition and display of ligand interactions; (ii) direct deposition by users of point mutations and their effects on ligand binding; (iii) refined snake and helix box residue diagram looks; and (iii) phylogenetic trees with receptor classification colour schemes. Under the hood, the entire GPCRdb front- and back-ends have been re-coded within one infrastructure, ensuring a smooth browsing experience and development. GPCRdb is available at http://www.gpcrdb.org/ and it's open source code at https://bitbucket.org/gpcr/protwis.

A new crystal structure fragment-based pharmacophore method for G protein-coupled receptors.

We have developed a new method for the building of pharmacophores for G protein-coupled receptors, a major drug target family. The method is a combination of the ligand- and target-based pharmacophore methods and founded on the extraction of structural fragments, interacting ligand moiety and receptor residue pairs, from crystal structure complexes. We describe the procedure to collect a library with more than 250 fragments covering 29 residue positions within the generic transmembrane binding pocket. We describe how the library fragments are recombined and inferred to build pharmacophores for new targets. A validating retrospective virtual screening of histamine H1 and H3 receptor pharmacophores yielded area-under-the-curves of 0.88 and 0.82, respectively. The fragment-based method has the unique advantage that it can be applied to targets for which no (homologous) crystal structures or ligands are known. 47% of the class A G protein-coupled receptors can be targeted with at least four-element pharmacophores. The fragment libraries can also be used to grow known ligands or for rotamer refinement of homology models. Researchers can download the complete fragment library or a subset matching their receptor of interest using our new tool in GPCRDB.

Differential G protein activation by the long and short isoforms of the dopamine D2 receptor.

BACKGROUND AND PURPOSE: The dopamine D2 receptor is expressed as a short (D2S) and a long (D2L) isoform with 29 additional amino acids in the third intracellular loop. The D2S isoform shows higher presynaptic expression than the D2L isoform, and decreased D2S expression has recently been linked to an increased risk for schizophrenia. Here, we present the first investigation, at receptor isoform level, of kinetic differences in the G protein activation profiles of the D2S, compared with the D2L isoform. EXPERIMENTAL APPROACH: We employed a NanoBRET-based approach to G protein dissociation to interrogate the time-resolved coupling profile of 3×HA-tagged D2L and D2S to Gαi/o/z proteins in vitro. KEY RESULTS: Using dopamine as a D2 receptor agonist, we observed a more pronounced activation of Gαo and Gαz than Gαi proteins by D2L compared with D2S. This differentiation was not observed for D2S, which activated Gαo and Gαz with lower efficacy than D2L. These signalling differences were preserved on second messenger level and were not due to differences in receptor expression. Expanding to a set of seven full and partial D2 receptor agonists showed these effects were not restricted to dopamine but rather a mutual, receptor-associated property. Contrasting this trend, we found that D2S activated G proteins faster than D2L upon full receptor activation. CONCLUSION AND IMPLICATIONS: The findings highlight that both D2L and D2S are mechanistically able to activate all non-visual Gαi/o proteins. Thereby, they add to previous reports about isoform-specificity to certain Gαi/o proteins observed in specific cell types.

Characterization of genetic variants of GIPR reveals a contribution of ß-arrestin to metabolic phenotypes.

Incretin-based therapies are highly successful in combatting obesity and type 2 diabetes1. Yet both activation and inhibition of the glucose-dependent insulinotropic polypeptide (GIP) receptor (GIPR) in combination with glucagon-like peptide-1 (GLP-1) receptor (GLP-1R) activation have resulted in similar clinical outcomes, as demonstrated by the GIPR-GLP-1R co-agonist tirzepatide2 and AMG-133 (ref. 3) combining GIPR antagonism with GLP-1R agonism. This underlines the importance of a better understanding of the GIP system. Here we show the necessity of ß-arrestin recruitment for GIPR function, by combining in vitro pharmacological characterization of 47 GIPR variants with burden testing of clinical phenotypes and in vivo studies. Burden testing of variants with distinct ligand-binding capacity, Gs activation (cyclic adenosine monophosphate production) and ß-arrestin 2 recruitment and internalization shows that unlike variants solely impaired in Gs signalling, variants impaired in both Gs and ß-arrestin 2 recruitment contribute to lower adiposity-related traits. Endosomal Gs-mediated signalling of the variants shows a ß-arrestin dependency and genetic ablation of ß-arrestin 2 impairs cyclic adenosine monophosphate production and decreases GIP efficacy on glucose control in male mice. This study highlights a crucial impact of ß-arrestins in regulating GIPR signalling and overall preservation of biological activity that may facilitate new developments in therapeutic targeting of the GIPR system.

GLP-1R signaling neighborhoods associate with the susceptibility to adverse drug reactions of incretin mimetics.

G protein-coupled receptors are important drug targets that engage and activate signaling transducers in multiple cellular compartments. Delineating therapeutic signaling from signaling associated with adverse events is an important step towards rational drug design. The glucagon-like peptide-1 receptor (GLP-1R) is a validated target for the treatment of diabetes and obesity, but drugs that target this receptor are a frequent cause of adverse events. Using recently developed biosensors, we explored the ability of GLP-1R to activate 15 pathways in 4 cellular compartments and demonstrate that modifications aimed at improving the therapeutic potential of GLP-1R agonists greatly influence compound efficacy, potency, and safety in a pathway- and compartment-selective manner. These findings, together with comparative structure analysis, time-lapse microscopy, and phosphoproteomics, reveal unique signaling signatures for GLP-1R agonists at the level of receptor conformation, functional selectivity, and location bias, thus associating signaling neighborhoods with functionally distinct cellular outcomes and clinical consequences.

Rare Heterozygous Loss-of-Function Variants in the Human GLP-1 Receptor Are Not Associated With Cardiometabolic Phenotypes.

CONTEXT: Lost glucagon-like peptide 1 receptor (GLP-1R) function affects human physiology. OBJECTIVE: This work aimed to identify coding nonsynonymous GLP1R variants in Danish individuals to link their in vitro phenotypes and clinical phenotypic associations. METHODS: We sequenced GLP1R in 8642 Danish individuals with type 2 diabetes or normal glucose tolerance and examined the ability of nonsynonymous variants to bind GLP-1 and to signal in transfected cells via cyclic adenosine monophosphate (cAMP) formation and ß-arrestin recruitment. We performed a cross-sectional study between the burden of loss-of-signaling (LoS) variants and cardiometabolic phenotypes in 2930 patients with type 2 diabetes and 5712 participants in a population-based cohort. Furthermore, we studied the association between cardiometabolic phenotypes and the burden of the LoS variants and 60 partly overlapping predicted loss-of-function (pLoF) GLP1R variants found in 330 566 unrelated White exome-sequenced participants in the UK Biobank cohort. RESULTS: We identified 36 nonsynonymous variants in GLP1R, of which 10 had a statistically significant loss in GLP-1-induced cAMP signaling compared to wild-type. However, no association was observed between the LoS variants and type 2 diabetes, although LoS variant carriers had a minor increased fasting plasma glucose level. Moreover, pLoF variants from the UK Biobank also did not reveal substantial cardiometabolic associations, despite a small effect on glycated hemoglobin A1c. CONCLUSION: Since no homozygous LoS nor pLoF variants were identified and heterozygous carriers had similar cardiometabolic phenotype as noncarriers, we conclude that GLP-1R may be of particular importance in human physiology, due to a potential evolutionary intolerance of harmful homozygous GLP1R variants.

Common coupling map advances GPCR-G protein selectivity.

Two-thirds of human hormones and one-third of clinical drugs act on membrane receptors that couple to G proteins to achieve appropriate functional responses. While G protein transducers from literature are annotated in the Guide to Pharmacology database, two recent large-scale datasets now expand the receptor-G protein 'couplome'. However, these three datasets differ in scope and reported G protein couplings giving different coverage and conclusions on G protein-coupled receptor (GPCR)-G protein signaling. Here, we report a common coupling map uncovering novel couplings supported by both large-scale studies, the selectivity/promiscuity of GPCRs and G proteins, and how the co-coupling and co-expression of G proteins compare to the families from phylogenetic relationships. The coupling map and insights on GPCR-G protein selectivity will catalyze advances in receptor research and cellular signaling toward the exploitation of G protein signaling bias in design of safer drugs.

Ligand-directed bias of G protein signaling at the dopamine D2 receptor.

G-protein-coupled receptors (GPCRs) represent the largest family of drug targets. Upon activation, GPCRs signal primarily via a diverse set of heterotrimeric G proteins. Most GPCRs can couple to several different G protein subtypes. However, how drugs act at GPCRs contributing to the selectivity of G protein recognition is poorly understood. Here, we examined the G protein selectivity profile of the dopamine D2 receptor (D2), a GPCR targeted by antipsychotic drugs. We show that D2 discriminates between six individual members of the Gi/o family, and its profile of functional selectivity is remarkably different across its ligands, which all engaged D2 with a distinct G protein coupling pattern. Using structural modeling, receptor mutagenesis, and pharmacological evaluation, we identified residues in the D2 binding pocket that shape these ligand-directed biases. We further provide pharmacogenomic evidence that natural variants in D2 differentially affect its G protein biases in response to different ligands.

Psychopharmacological treatment patterns prior to a schizophrenia diagnosis: A Danish nationwide study.

BACKGROUND: Many patients with schizophrenia experience psychiatric symptoms long before being diagnosed. We investigated patterns of pre-diagnostic psychopharmacological treatment in individuals diagnosed with first-episode schizophrenia during the last two decades. DESIGN: Using Danish nationwide healthcare registers, we identified all individuals aged ≥10 years registered with their first ICD-10 schizophrenia diagnosis between January 1999 and March 2019. For each calendar year from 1999 to 2019, we calculated the proportion of patients - among those having received their first schizophrenia diagnosis in the respective calendar year - who redeemed prescriptions for various psychotropics in the two years preceding the schizophrenia diagnosis. We calculated proportions of all pre-diagnostic prescriptions since 1995 for a sub-population diagnosed 2011-2019 and for an age- and sex-matched reference group without schizophrenia. RESULTS: Among 33,361 individuals with schizophrenia (58 % males), the schizophrenia incidence rate was stable during the study period but the mean age at diagnosis decreased by >10 years. In the two pre-diagnostic years, 69 % received psychopharmacological treatment (52 % antipsychotics, 40 % antidepressants). This was stable between 1999 and 2019. Among 14,425 individuals diagnosed 2011-2019, psychotropic drug use was observed among 14-20 % between 24 and 10 years before the diagnosis, being four times higher than the reference group. Particularly antipsychotic and antidepressant drug use increased steadily during the ten pre-diagnostic years. CONCLUSIONS: Pre-diagnostic psychotropic drug use in schizophrenia was frequent but stable between 1999 and 2019 despite an earlier identification of schizophrenia patients. Our findings emphasize the continued importance of thorough diagnostic interviews, particularly among patients in need of antipsychotic treatment.

Effector membrane translocation biosensors reveal G protein and ßarrestin coupling profiles of 100 therapeutically relevant GPCRs.

The recognition that individual GPCRs can activate multiple signaling pathways has raised the possibility of developing drugs selectively targeting therapeutically relevant ones. This requires tools to determine which G proteins and ßarrestins are activated by a given receptor. Here, we present a set of BRET sensors monitoring the activation of the 12 G protein subtypes based on the translocation of their effectors to the plasma membrane (EMTA). Unlike most of the existing detection systems, EMTA does not require modification of receptors or G proteins (except for Gs). EMTA was found to be suitable for the detection of constitutive activity, inverse agonism, biased signaling and polypharmacology. Profiling of 100 therapeutically relevant human GPCRs resulted in 1500 pathway-specific concentration-response curves and revealed a great diversity of coupling profiles ranging from exquisite selectivity to broad promiscuity. Overall, this work describes unique resources for studying the complexities underlying GPCR signaling and pharmacology.