Cryo-electron microscopy structure of the glucagon receptor with a dual-agonist peptide

Molecular features of the ligand-free GLP-1R, GCGR and GIPR in complex with Gs proteins

Class B1 G protein-coupled receptors (GPCRs) are important regulators of many physiological functions such as glucose homeostasis, which is mainly mediated by three peptide hormones, i.e., glucagon-like peptide-1 (GLP-1), glucagon (GCG), and glucose-dependent insulinotropic polypeptide (GIP). They trigger a cascade of signaling events leading to the formation of an active agonist-receptor-G protein complex. However, intracellular signal transducers can also activate the receptor independent of extracellular stimuli, suggesting an intrinsic role of G proteins in this process. Here, we report cryo-electron microscopy structures of the human GLP-1 receptor (GLP-1R), GCG receptor (GCGR), and GIP receptor (GIPR) in complex with Gs proteins without the presence of cognate ligands. These ligand-free complexes share a similar intracellular architecture to those bound by endogenous peptides, in which, the Gs protein alone directly opens the intracellular binding cavity and rewires the extracellular orthosteric pocket to stabilize the receptor in a state unseen before. While the peptide-binding site is partially occupied by the inward folded transmembrane helix 6 (TM6)-extracellular loop 3 (ECL3) juncture of GIPR or a segment of GCGR ECL2, the extracellular portion of GLP-1R adopts a conformation close to the active state. Our findings offer valuable insights into the distinct activation mechanisms of these three important receptors. It is possible that in the absence of a ligand, the intracellular half of transmembrane domain is mobilized with the help of Gs protein, which in turn rearranges the extracellular half to form a transitional conformation, facilitating the entry of the peptide N-terminus.

Hormone Analogues with Unique Signaling Profiles from Replacement of α-Residue Triads with β/γ Diads

We have applied an underexplored backbone modification strategy to generate new analogues of peptides that activate two clinically important class B1 G protein-coupled receptors (GPCRs). Most peptide modification strategies involve changing side chains or, less commonly, changing the configuration at side chain-bearing carbons (i.e., l residues replaced by d residues). In contrast, backbone modifications alter the number of backbone atoms and the identities of backbone atoms relative to a poly-α-amino acid backbone. Starting from the peptide agonists PTH(1-34) (the first 34 residues of the parathyroid hormone, used clinically as the drug teriparatide) and glucagon-like peptide-1 (7-36) (GLP-1(7-36)), we replaced native α-residue triads with a diad composed of a β-amino acid residue and a γ-amino acid residue. The β/γ diad retains the number of backbone atoms in the ααα triad. Because the β and γ residue each bear a single side chain, we implemented ααα→βγ replacements at sites that contained a Gly residue (i.e., at α-residue triads that presented only two side chains). All seven of the α/β/γ-peptides derived from PTH(1-34) or GLP-1(7-36) bind to the cognate receptor (the PTHR1 or the GLP-1R), but they vary considerably in their activity profiles. Outcomes include functional mimicry of the all-α agonist, receptor-selective agonist activity, biased agonism, or strong binding with weak activation, which could lead to antagonist development. Collectively, these findings demonstrate that ααα→βγ replacements, which are easily implemented via solid-phase synthesis, can generate peptide hormone analogues that display unique and potentially useful signaling behavior.

Structural analysis of the dual agonism at GLP-1R and GCGR

Glucagon-like peptide-1 receptor (GLP-1R) and glucagon receptor (GCGR), two members of class B1 G protein-coupled receptors, play important roles in glucose homeostasis and energy metabolism. They share a high degree of sequence homology but have different functionalities. Unimolecular dual agonists of both receptors developed recently displayed better clinical efficacies than that of monotherapy. To study the underlying molecular mechanisms, we determined high-resolution cryo-electron microscopy structures of GLP-1R or GCGR in complex with heterotrimeric Gs protein and three GLP-1R/GCGR dual agonists including peptide 15, MEDI0382 (cotadutide) and SAR425899 with variable activating profiles at GLP-1R versus GCGR. Compared with related structures reported previously and supported by our published pharmacological data, key residues responsible for ligand recognition and dual agonism were identified. Analyses of peptide conformational features revealed a difference in side chain orientations within the first three residues, indicating that distinct engagements in the deep binding pocket are required to achieve receptor selectivity. The middle region recognizes extracellular loop 1 (ECL1), ECL2, and the top of transmembrane helix 1 (TM1) resulting in specific conformational changes of both ligand and receptor, especially the dual agonists reshaped ECL1 conformation of GLP-1R relative to that of GCGR, suggesting an important role of ECL1 interaction in executing dual agonism. Structural investigation of lipid modification showed a better interaction between lipid moiety of MEDI0382 and TM1-TM2 cleft, in line with its increased potency at GCGR than SAR425899. Together, the results provide insightful information for the design and development of improved therapeutics targeting these two receptors simultaneously.

Differential Responses of the GLP-1 and GLP-2 Receptors to N-Terminal Modification of a Dual Agonist

Class B1 G protein-coupled receptors (GPCRs), collectively, respond to a diverse repertoire of extracellular polypeptide agonists and transmit the encoded messages to cytosolic partners. To fulfill these tasks, these highly mobile receptors must interconvert among conformational states in response to agonists. We recently showed that conformational mobility in polypeptide agonists themselves plays a role in activation of one class B1 GPCR, the receptor for glucagon-like peptide-1 (GLP-1). Exchange between helical and nonhelical conformations near the N-termini of agonists bound to the GLP-1R was revealed to be critical for receptor activation. Here, we ask whether agonist conformational mobility plays a role in the activation of a related receptor, the GLP-2R. Using variants of the hormone GLP-2 and the designed clinical agonist glepaglutide (GLE), we find that the GLP-2R is quite tolerant of variations in α-helical propensity near the agonist N-terminus, which contrasts with signaling at the GLP-1R. A fully α-helical conformation of the bound agonist may be sufficient for GLP-2R signal transduction. GLE is a GLP-2R/GLP-1R dual agonist, and the GLE system therefore enables direct comparison of the responses of these two GPCRs to a single set of agonist variants. This comparison supports the conclusion that the GLP-1R and GLP-2R differ in their response to variations in helical propensity near the agonist N-terminus. The data offer a basis for development of new hormone analogues with distinctive and potentially useful activity profiles; for example, one of the GLE analogues is a potent agonist of the GLP-2R but also a potent antagonist of the GLP-1R, a novel form of polypharmacology.

Effects of site-directed mutagenesis of GLP-1 and glucagon receptors on signal transduction activated by dual and triple agonists

The paradigm of one drug against multiple targets, known as unimolecular polypharmacology, offers the potential to improve efficacy while overcoming some adverse events associated with the treatment. This approach is best exemplified by targeting two or three class B1 G protein-coupled receptors, namely, glucagon-like peptide-1 receptor (GLP-1R), glucagon receptor (GCGR) and glucose-dependent insulinotropic polypeptide receptor for treatment of type 2 diabetes and obesity. Some of the dual and triple agonists have already shown initial successes in clinical trials, although the molecular mechanisms underlying their multiplexed pharmacology remain elusive. In this study we employed structure-based site-directed mutagenesis together with pharmacological assays to compare agonist efficacy across two key signaling pathways, cAMP accumulation and ERK1/2 phosphorylation (pERK1/2). Three dual agonists (peptide 15, MEDI0382 and SAR425899) and one triple agonist (peptide 20) were evaluated at GLP-1R and GCGR, relative to the native peptidic ligands (GLP-1 and glucagon). Our results reveal the existence of residue networks crucial for unimolecular agonist-mediated receptor activation and their distinct signaling patterns, which might be useful to the rational design of biased drug leads.

Structural Understanding of Peptide-Bound G Protein-Coupled Receptors: Peptide-Target Interactions

The activation of G protein-coupled receptors (GPCRs) is triggered by ligand binding to their orthosteric sites, which induces ligand-specific conformational changes. Agonists and antagonists bound to GPCR orthosteric sites provide detailed information on ligand-binding modes. Among these, peptide ligands play an instrumental role in GPCR pharmacology and have attracted increased attention as therapeutic drugs. The recent breakthrough in GPCR structural biology has resulted in the remarkable availability of peptide-bound GPCR complexes. Despite the several structural similarities shared by these receptors, they exhibit distinct features in terms of peptide recognition and receptor activation. From this perspective, we have summarized the current status of peptide-bound GPCR structural complexes, largely focusing on the interactions between the receptor and its peptide ligand at the orthosteric site. In-depth structural investigations have yielded valuable insights into the molecular mechanisms underlying peptide recognition. This study would contribute to the discovery of GPCR peptide drugs with improved therapeutic effects.

Lipid Modulation of a Class B GPCR: Elucidating the Modulatory Role of PI(4,5)P2 Lipids

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) lipids have been shown to stabilize an active conformation of class A G-protein coupled receptors (GPCRs) through a conserved binding site, not present in class B GPCRs. For class B GPCRs, previous molecular dynamics (MD) simulation studies have shown PI(4,5)P₂ interacting with the Glucagon receptor (GCGR), which constitutes an important target for diabetes and obesity therapeutics. In this work, we applied MD simulations supported by native mass spectrometry (nMS) to study lipid interactions with GCGR. We demonstrate how tail composition plays a role in modulating the binding of PI(4,5)P₂ lipids to GCGR. Specifically, we find the PI(4,5)P₂ lipids to have a higher affinity toward the inactive conformation of GCGR. Interestingly, we find that in contrast to class A GPCRs, PI(4,5)P₂ appear to stabilize the inactive conformation of GCGR through a binding site conserved across class B GPCRs but absent in class A GPCRs. This suggests differences in the regulatory function of PI(4,5)P₂ between class A and class B GPCRs.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

Understanding the Allosteric Modulation of PTH1R by a Negative Allosteric Modulator

The parathyroid hormone type 1 receptor (PTH1R) acts as a canonical class B G protein-coupled receptor, regulating crucial functions including calcium homeostasis and bone formation. The identification and development of PTH1R non-peptide allosteric modulators have obtained widespread attention. It has been found that a negative allosteric modulator (NAM) could inhibit the activation of PTH1R, but the implied mechanism remains unclear. Herein, extensive molecular dynamics simulations together with multiple analytical approaches are utilized to unravel the mechanism of PTH1R allosteric inhibition. The results suggest that the binding of NAM destabilizes the structure of the PTH1R-PTH-spep/qpep (the C terminus of Gs/Gq proteins) complexes. Moreover, the presence of NAM weakens the binding of PTH/peps (spep and qpep) and PTH1R. The intra- and inter-molecular couplings are also weakened in PTH1R upon NAM binding. Interestingly, compared with our previous study of the positive allosteric effects induced by extracellular Ca²⁺, the enhanced correlation between the PTH and G-protein binding sites is significantly reduced by the replacement of this negative allosteric regulator. Our findings might contribute to the development of new therapeutic agents for diseases caused by the abnormal activation of PTH1R.

Targeting VIP and PACAP Receptor Signaling: New Insights into Designing Drugs for the PACAP Subfamily of Receptors

Pituitary Adenylate Cyclase-Activating Peptide (PACAP) and Vasoactive Intestinal Peptide (VIP) are neuropeptides involved in a diverse array of physiological and pathological processes through activating the PACAP subfamily of class B1 G protein-coupled receptors (GPCRs): VIP receptor 1 (VPAC1R), VIP receptor 2 (VPAC2R), and PACAP type I receptor (PAC1R). VIP and PACAP share nearly 70% amino acid sequence identity, while their receptors PAC1R, VPAC1R, and VPAC2R share 60% homology in the transmembrane regions of the receptor. PACAP binds with high affinity to all three receptors, while VIP binds with high affinity to VPAC1R and VPAC2R, and has a thousand-fold lower affinity for PAC1R compared to PACAP. Due to the wide distribution of VIP and PACAP receptors in the body, potential therapeutic applications of drugs targeting these receptors, as well as expected undesired side effects, are numerous. Designing selective therapeutics targeting these receptors remains challenging due to their structural similarities. This review discusses recent discoveries on the molecular mechanisms involved in the selectivity and signaling of the PACAP subfamily of receptors, and future considerations for therapeutic targeting.

Critical Extracellular Ca2+ Dependence of the Binding between PTH1R and a G-Protein Peptide Revealed by MD Simulations

The parathyroid hormone type 1 receptor (PTH1R), a canonical class B GPCR, is regulated by a positive allosteric modulator, extracellular Ca²⁺. Calcium ions prolong the residence time of PTH on the PTH1R, leading to increased receptor activation and duration of cAMP signaling. But the essential mechanism of the allosteric behavior of PTH1R is not fully understood. Here, extensive molecular dynamics (MD) simulations are performed for the PTH1R-G-protein combinations with and without Ca²⁺ to describe how calcium ions allosterically engage receptor-G-protein coupling. We find that the binding of Ca²⁺ stabilizes the conformation of the PTH1R-PTH-spep (the α5 helix of Gs protein) complex, especially the extracellular loop 1 (ECL1). Moreover, the MM-GBSA result indicates that Ca²⁺ allosterically promotes the interaction between PTH1R and spep, consistent with the observation of steered molecular dynamics (SMD) simulations. We further illuminate the possible allosteric signaling pathway from the stable Ca²⁺-coupling site to the intracellular G-protein binding site. These results unveil structural determinants for Ca²⁺ allosterism in the PTH1R-PTH-spep complex and give insights into pluridimensional GPCR signaling regulated by calcium ions.

Structural insights into multiplexed pharmacological actions of tirzepatide and peptide 20 at the GIP, GLP-1 or glucagon receptors

Glucose homeostasis, regulated by glucose-dependent insulinotropic polypeptide (GIP), glucagon-like peptide-1 (GLP-1) and glucagon (GCG) is critical to human health. Several multi-targeting agonists at GIPR, GLP-1R or GCGR, developed to maximize metabolic benefits with reduced side-effects, are in clinical trials to treat type 2 diabetes and obesity. To elucidate the molecular mechanisms by which tirzepatide, a GIPR/GLP-1R dual agonist, and peptide 20, a GIPR/GLP-1R/GCGR triagonist, manifest their multiplexed pharmacological actions over monoagonists such as semaglutide, we determine cryo-electron microscopy structures of tirzepatide-bound GIPR and GLP-1R as well as peptide 20-bound GIPR, GLP-1R and GCGR. The structures reveal both common and unique features for the dual and triple agonism by illustrating key interactions of clinical relevance at the near-atomic level. Retention of glucagon function is required to achieve such an advantage over GLP-1 monotherapy. Our findings provide valuable insights into the structural basis of functional versatility of tirzepatide and peptide 20.

Multi-targeting agonists at GIPR, GLP-1R or GCGR are pursued vigorously. Here, the authors report cryo-EM structures of tirzepatide-bound GIPR and GLP-1R, peptide 20-bound GIPR, GLP-1R and GCGR, revealing the molecular basis of their multiplexed pharmacological actions.

Dynamics of GLP-1R peptide agonist engagement are correlated with kinetics of G protein activation

The glucagon-like peptide-1 receptor (GLP-1R) has broad physiological roles and is a validated target for treatment of metabolic disorders. Despite recent advances in GLP-1R structure elucidation, detailed mechanistic understanding of how different peptides generate profound differences in G protein-mediated signalling is still lacking. Here we combine cryo-electron microscopy, molecular dynamics simulations, receptor mutagenesis and pharmacological assays, to interrogate the mechanism and consequences of GLP-1R binding to four peptide agonists; glucagon-like peptide-1, oxyntomodulin, exendin-4 and exendin-P5. These data reveal that distinctions in peptide N-terminal interactions and dynamics with the GLP-1R transmembrane domain are reciprocally associated with differences in the allosteric coupling to G proteins. In particular, transient interactions with residues at the base of the binding cavity correlate with enhanced kinetics for G protein activation, providing a rationale for differences in G protein-mediated signalling efficacy from distinct agonists.

The glucagon-like peptide-1 receptor (GLP-1R) can be targeted in the treatment of diabetes, obesity and other metabolic disorders. Here, the authors assess the molecular mechanisms of peptide agonists binding to GLP-1R and the responses elucidated by these ligands, including distinct kinetics of G protein activation.

The current significance and prospects for the use of dual receptor agonism GLP-1/Glucagon

The therapeutic arsenal for treating type 2 diabetes mellitus (T2DM) has been enriched recently with the inclusion of type 1 glucagon-like peptide (GLP-1). GLP-1 receptor agonists (RA) secondarily reduce appetite, decrease gastric emptying, and reduce body weight. This effect has been used to treat overweight/obesity, especially with comorbidities associated with T2DM. However, the first formulations and adverse effects gradually gave way to new formulations with fewer unpleasant effects and a more extended period of action (weekly subcutaneous administration and even oral administration), which improved the acceptance and adherence to the treatment. Therefore, titration of GLP-1RA should be done gradually. Furthermore, when side effects are consistent and intolerable after weeks/months of titration, a lower dose or a combination of antidiabetic therapies should be implemented, avoiding treatment interruption. The effort to produce increasingly powerful molecules with fewer side effects is the driving force behind the pharmaceutical industry. The unimolecular dual agonism GLP-1RA plus glucagon receptor agonism (GRA) represents an updated pharmacological indication for controlling blood glucose levels in treating T2DM and its comorbidities, showing better effects with less adverse impact than mono GLP-1RA. There are currently different proposals in this way by different laboratories. Nevertheless, the experimental results are promising and show that soon, we will have the contribution of new drugs for the treatment of T2DM.

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 (D63^ECDN, P86^ECDS, V96^ECDE, G125^ECDC, R225^3.30H, R308^5.40W, V368^6.59M, and R378^7.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.

Partial agonism improves the anti-hyperglycaemic efficacy of an oxyntomodulin-derived GLP-1R/GCGR co-agonist

OBJECTIVE

Glucagon-like peptide-1 and glucagon receptor (GLP-1R/GCGR) co-agonism can maximise weight loss and improve glycaemic control in type 2 diabetes and obesity. In this study, we investigated the cellular and metabolic effects of modulating the balance between G protein and β-arrestin-2 recruitment at GLP-1R and GCGR using oxyntomodulin (OXM)-derived co-agonists. This strategy has been previously shown to improve the duration of action of GLP-1R mono-agonists by reducing target desensitisation and downregulation.

METHODS

Dipeptidyl dipeptidase-4 (DPP-4)-resistant OXM analogues were generated and assessed for a variety of cellular readouts. Molecular dynamic simulations were used to gain insights into the molecular interactions involved. In vivo studies were performed in mice to identify the effects on glucose homeostasis and weight loss.

RESULTS

Ligand-specific reductions in β-arrestin-2 recruitment were associated with slower GLP-1R internalisation and prolonged glucose-lowering action in vivo. The putative benefits of GCGR agonism were retained, with equivalent weight loss compared to the GLP-1R mono-agonist liraglutide despite a lesser degree of food intake suppression. The compounds tested showed only a minor degree of biased agonism between G protein and β-arrestin-2 recruitment at both receptors and were best classified as partial agonists for the two pathways measured.

CONCLUSIONS

Diminishing β-arrestin-2 recruitment may be an effective way to increase the therapeutic efficacy of GLP-1R/GCGR co-agonists. These benefits can be achieved by partial rather than biased agonism.

Insights into agonist-elicited activation of the human glucose-dependent insulinotropic polypeptide receptor

Glucose-dependent insulinotropic polypeptide (GIP) and its receptor (GIPR) are part of the incretin system that regulates glucose homeostasis. A series of GIPR residues putatively important for ligand binding and receptor activation were mutated and pharmacologically evaluated using GIPR selective agonists in cAMP accumulation, ERK1/2 phosphorylation (pERK1/2) and β-arrestin 2 recruitment assays. The impact of mutation on ligand efficacy was determined by operational modelling of experimental data for each mutant, with results mapped onto the full-length, active-state GIPR structure. Two interaction networks, comprising transmembrane helix (TM) 7, TM1 and TM2, and extracellular loop (ECL) 2, TM5 and ECL3 were revealed, respectively. Both networks were critical for Gα_s-mediated cAMP accumulation and the recruitment of β-arrestin 2, however, cAMP response was more sensitive to alanine substitution, with most mutated residues displaying reduced signaling. Unlike the other two assays, activation of ERK1/2 was largely independent of the network involving ECL2, TM5 and ECL3, indicating that pERK1/2 is at least partially distinct from Gα_s or β-arrestin pathways and this network is also crucial for potential biased agonism at GIPR. Collectively, our work advances understanding of the structure-function relationship of GIPR and provides a framework for the design and/or interpretation of GIP analogues with unique signaling profiles.

Structural insights into hormone recognition by the human glucose-dependent insulinotropic polypeptide receptor

Glucose-dependent insulinotropic polypeptide (GIP) is a peptide hormone that exerts crucial metabolic functions by binding and activating its cognate receptor, GIPR. As an important therapeutic target, GIPR has been subjected to intensive structural studies without success. Here, we report the cryo-EM structure of the human GIPR in complex with GIP and a Gs heterotrimer at a global resolution of 2.9 Å. GIP adopts a single straight helix with its N terminus dipped into the receptor transmembrane domain (TMD), while the C terminus is closely associated with the extracellular domain and extracellular loop 1. GIPR employs conserved residues in the lower half of the TMD pocket to recognize the common segments shared by GIP homologous peptides, while uses non-conserved residues in the upper half of the TMD pocket to interact with residues specific for GIP. These results provide a structural framework of hormone recognition and GIPR activation.

Evolving cryo-EM structural approaches for GPCR drug discovery

G protein-coupled receptors (GPCRs) are the largest class of cell surface drug targets. Advances in stabilization of GPCR:transducer complexes, together with improvements in cryoelectron microscopy (cryo-EM) have recently been applied to structure-assisted drug design for GPCR agonists. Nonetheless, limitations in the commercial application of these approaches, including the use of nanobody 35 (Nb35) to aid complex stabilization and the high cost of 300 kV imaging, have restricted broad application of cryo-EM in drug discovery. Here, using the PF 06882961-bound GLP-1R as exemplar, we validated the formation of stable complexes with a modified Gs protein in the absence of Nb35. In parallel, we compare 200 versus 300 kV image acquisition using a Falcon 4 or K3 direct electron detector. Moreover, the 200 kV Glacios-Falcon 4 yielded a 3.2 Å map with clear density for bound drug and multiple structurally ordered waters. Our work paves the way for broader commercial application of cryo-EM for GPCR drug discovery.

Ligand-Receptor Interactions and Machine Learning in GCGR and GLP-1R Drug Discovery

The large amount of data that has been collected so far for G protein-coupled receptors requires machine learning (ML) approaches to fully exploit its potential. Our previous ML model based on gradient boosting used for prediction of drug affinity and selectivity for a receptor subtype was compared with explicit information on ligand-receptor interactions from induced-fit docking. Both methods have proved their usefulness in drug response predictions. Yet, their successful combination still requires allosteric/orthosteric assignment of ligands from datasets. Our ligand datasets included activities of two members of the secretin receptor family: GCGR and GLP-1R. Simultaneous activation of two or three receptors of this family by dual or triple agonists is not a typical kind of information included in compound databases. A precise allosteric/orthosteric ligand assignment requires a continuous update based on new structural and biological data. This data incompleteness remains the main obstacle for current ML methods applied to class B GPCR drug discovery. Even so, for these two class B receptors, our ligand-based ML model demonstrated high accuracy (5-fold cross-validation Q² > 0.63 and Q² > 0.67 for GLP-1R and GCGR, respectively). In addition, we performed a ligand annotation using recent cryogenic-electron microscopy (cryo-EM) and X-ray crystallographic data on small-molecule complexes of GCGR and GLP-1R. As a result, we assigned GLP-1R and GCGR actives deposited in ChEMBL to four small-molecule binding sites occupied by positive and negative allosteric modulators and a full agonist. Annotated compounds were added to our recently released repository of GPCR data.

Engineering of Challenging G Protein-Coupled Receptors for Structure Determination and Biophysical Studies

Membrane proteins such as G protein-coupled receptors (GPCRs) exert fundamental biological functions and are involved in a multitude of physiological responses, making these receptors ideal drug targets. Drug discovery programs targeting GPCRs have been greatly facilitated by the emergence of high-resolution structures and the resulting opportunities to identify new chemical entities through structure-based drug design. To enable the determination of high-resolution structures of GPCRs, most receptors have to be engineered to overcome intrinsic hurdles such as their poor stability and low expression levels. In recent years, multiple engineering approaches have been developed to specifically address the technical difficulties of working with GPCRs, which are now beginning to make more challenging receptors accessible to detailed studies. Importantly, successfully engineered GPCRs are not only valuable in X-ray crystallography, but further enable biophysical studies with nuclear magnetic resonance spectroscopy, surface plasmon resonance, native mass spectrometry, and fluorescence anisotropy measurements, all of which are important for the detailed mechanistic understanding, which is the prerequisite for successful drug design. Here, we summarize engineering strategies based on directed evolution to reduce workload and enable biophysical experiments of particularly challenging GPCRs.

Design and preparation of the class B G protein-coupled receptors GLP-1R and GCGR for 19 F-NMR studies in solution

The human glucagon-like peptide-1 receptor (GLP-1R) and the glucagon receptor (GCGR) are class B G protein-coupled receptors (GPCRs) that are activated by interactions with, respectively, the glucagon-like peptide-1 (GLP-1) and glucagon (GCG). These polypeptide hormones are involved in the regulation of lipid and cholic acid metabolism, and thus play an important role in the pathogenesis of glucose metabolism and diabetes mellitus, which attracts keen interest of these GPCRs as drug targets. GLP-1R and GCGR have therefore been extensively investigated by X-ray crystallography and cryo-electron microscopy (cryo-EM), so that their structures are well known. Here, we present the groundwork for using nuclear magnetic resonance (NMR) spectroscopy in solution to complement the molecular architectures with information on intramolecular dynamics and on the thermodynamics and kinetics of interactions with physiological ligands and extrinsic drug candidates. This includes the generation of novel, near-wild-type constructs of GLP-1R and GCGR, optimization of the solution conditions for NMR studies in detergent micelles and in nanodiscs, post-translational chemical introduction of fluorine-19 NMR probes, and sequence-specific assignments of the ¹⁹ F-labels attached to indigenous cysteines. Addition of the negative allosteric modulator (NAM) NNC0640 was critically important for obtaining the long-time stability needed for our NMR experiments, and we report on novel insights into the allosteric effects arising from binding of NNC0640 to the transmembrane domain of GLP-1R (GLP-1R[TMD]).

GIP receptor suppresses PAC1receptor-mediated neuronal differentiation via formation of a receptor heterocomplex

Pituitary adenylate cyclase-activating peptide (PACAP) receptor (PAC1R) is a class B Gprotein-coupled receptor (GPCR) that is widely expressed in the human body and is involved in neuronal differentiation. As class B GPCRs are known to form heterocomplexes with family members, we hypothesized that PAC1R mediates neuronal differentiation through interaction with a class B GPCR. We used the BRET assay to identify potential interactions between PAC1R and 11 class B GPCRs. Gastric inhibitory polypeptide receptor (GIPR) and secretin receptor were identified as putative binding partners of PAC1R. The effect of heterocomplex formation by PAC1R on receptor activation was evaluated with the cyclic (c)AMP, luciferase reporter, and calcium signaling assays; and the effects on receptor internalization and subcellular localization were examined by confocal microscopy. The results suggested he PAC1R/GIPR heterocomplex suppressed signaling events downstream of PAC1R, including cAMP production, serum response element and calcium signaling, and β-arrestin recruitment. Protein-protein interaction was analyzed in silico, and induction of neuronal differentiation by the PAC1R heterocomplex was assessed in SH-SY5Y neuronal cells by measure the morphological changes and marker genes expression by real-time quantitative PCR and western blot. Over-expression of GIPR suppressed PACAP/PAC1R-mediated neuronal differentiation and the differentiation markers expression in SH-SY5Y cells. GIPR regulates neuronal differentiation through heterocomplex formation with PAC1R.