Analyzing kinetic signaling data for G-protein-coupled receptors

Computational modelling of dynamic cAMP responses to GPCR agonists for exploration of GLP-1R ligand effects in pancreatic β-cells and neurons

The glucagon-like peptide-1 receptor (GLP-1R) is a class B G protein-coupled receptor (GPCR) which plays important physiological roles in insulin release and promoting fullness. GLP-1R agonists initiate cellular responses by cyclic AMP (cAMP) pathway signal transduction. Understanding of the potential of GLP-1R agonists in the treatment of type 2 diabetes may be advanced by considering the cAMP dynamics for agonists at GLP-1R in both pancreatic β-cells (important in insulin release) and neurons (important in appetite regulation). Receptor desensitisation in the cAMP pathway is known to be an important regulatory mechanism, with different ligands differentially promoting G protein activation and desensitisation. Here, we use mathematical modelling to quantify and understand experimentally obtained cAMP timecourses for two GLP-1R agonists, exendin-F1 (ExF1) and exendin-D3 (ExD3), which give markedly different signals in β-cells and neurons. We formulate an ordinary differential equation (ODE) model for the dynamics of cAMP signalling in response to G protein-coupled receptor (GPCR) ligands, encompassing ligand binding, receptor activation, G protein activation, desensitisation and second messenger generation. We validate our model initially by fitting to timecourse data for HEK293 cells, then proceed to parameterise the model for β-cells and neurons. Through numerical simulation and sensitivity studies, our analysis adds support to the hypothesis that ExF1 offers more potential glucose regulation benefit than ExD3 over long timescales via signalling in pancreatic β-cells, but that there is little difference between the two ligands in the potential appetite suppression effects offered via long-time signalling in neurons on the same timescales.

Importance of receptor expression in the classification of novel ligands at the M2 muscarinic acetylcholine receptor

BACKGROUND AND PURPOSE

Affinity-based, selective orthosteric ligands for the muscarinic acetylcholine receptors (mAChRs) are difficult to develop due to high sequence homology across the five subtypes. Selectivity can also be achieved via the selective activation of a particular subtype or signalling pathway. Promisingly, a prior study identified compounds 6A and 7A as functionally selective and Gi biased compounds at the M2 mAChR. Here, we have investigated the activation of individual G protein subfamilies and the downstream signalling profiles of 6A and 7A at the M2 mAChR.

EXPERIMENTAL APPROACH

G protein activation was measured with the TRUPATH assay in M2 mAChR FlpIn CHO cells. Activity in downstream signalling pathways was determined using the cAMP CAMYEL BRET sensor and assay of ERK 1/2 phosphorylation.

KEY RESULTS

M2 mAChRs coupled to Gɑi1, GɑoA and Gɑs, but not Gɑq, in response to canonical orthosteric agonists. Compounds 6A and 7A did not elicit any G protein activation, cAMP inhibition or stimulation, or ERK 1/2 phosphorylation. Instead, a Schild analysis indicates a competitive, antagonistic interaction of compounds 6A and 7A with ACh in the Gɑi1 activation assay. Overexpression of the M2 mAChR may suggest an expression-dependent activation profile of compounds 6A and 7A.

CONCLUSIONS AND IMPLICATIONS

These data confirm that the M2 mAChR preferentially couples to Gɑi/o and to a lesser extent to Gɑs in response to canonical orthosteric ligands. However, this study was not able to detect Gɑi bias of compounds 6A and 7A, highlighting the importance of cellular background when classifying new ligands.

LINKED ARTICLES

This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

Time-dependent ligand-receptor binding kinetics and functionality in a heterodimeric receptor model

GPCRs heteromerize both in CNS and non-CNS regions. The cell uses receptor heteromerization to modulate receptor functionality and to provide fine tuning of receptor signaling. In order for pharmacologists to explore these mechanisms for therapeutic purposes, quantitative receptor models are needed. We have developed a time-dependent model of the binding kinetics and functionality of a preformed heterodimeric receptor involving two drugs. Two cases were considered: both or only one of the drugs are in excess with respect to the total concentration of the receptor. The latter case can be applied to those situations in which a drug causes unwanted side effects that need to be reduced by decreasing its concentration. The required efficacy can be maintained by the allosteric effects mutually exerted by the two drugs in the two-drug combination system. We discuss this concept assuming that the drug causing unwanted side effects is an opioid and that analgesia is the therapeutic effect. As additional points, allosteric modulation by endogenous compounds and synthetic bivalent ligands was included in the study. Receptor heteromerization offers a mechanistic understanding and quantification of the pharmacological effects elicited by combinations of two drugs at different doses and with different efficacies and cooperativity effects, thus providing a conceptual framework for drug combination therapy.

Positive allosteric modulation of the cannabinoid CB1 receptor potentiates endocannabinoid signalling and changes ERK1/2 phosphorylation kinetics

BACKGROUND AND PURPOSE

Activation of CB1 by exogenous agonists causes adverse effects in vivo. Positive allosteric modulation may offer improved therapeutic potential and a reduced on-target adverse effect profile compared with orthosteric agonists, due to reduced desensitisation/tolerance, but this has not been directly tested. This study investigated the ability of PAMs/ago-PAMs to induce receptor regulation pathways, including desensitisation and receptor internalisation.

EXPERIMENTAL APPROACH

Bioluminescence resonance energy transfer (BRET) assays in HEK293 cells were performed to investigate G protein dissociation, ERK1/2 phosphorylation and β-arrestin 2 translocation, while immunocytochemistry was performed to measure internalisation of CB1 in response to the PAMs ZCZ011, GAT229 and ABD1236 alone and in combination with the orthosteric agonists AEA, 2-AG, and AMB-FUBINACA.

KEY RESULTS

ZCZ011, GAT229 and ABD1236 were allosteric agonists in all pathways tested. The ago-PAM ZCZ011 induced a biphasic ERK1/2 phosphorylation time course compared to transient activation by orthosteric agonists. In combination with 2-AG but not AEA or AMB-FUBINACA, ZCZ011 and ABD1236 caused the transient peak of ERK1/2 phosphorylation to become sustained. All PAMs increased the potency and efficacy of AEA-induced signalling in all pathways tested; however, no notable potentiation of 2-AG or AMB-FUBINACA was observed.

CONCLUSION AND IMPLICATIONS

Ago-PAMs can potentiate endocannabinoid CB1 agonism by AEA to a larger extent compared with 2-AG. However, all compounds were found to be allosteric agonists and induce activation of CB1 in the absence of endocannabinoid, including β-arrestin 2 recruitment and internalisation. Thus, the spatiotemporal signalling of endogenous cannabinoids will not be retained in vivo.

GPCR kinase subtype requirements for arrestin-2 and -3 translocation to the cannabinoid CB1 receptor and the consequences on G protein signalling

Arrestins are key negative regulators of G Protein-Coupled Receptors (GPCRs) through mediation of G protein desensitisation and receptor internalisation. Arrestins can also contribute to signal transduction by scaffolding downstream signalling effectors for activation. GPCR kinase (GRK) enzymes phosphorylate the intracellular C-terminal domain, or intracellular loop regions of GPCRs to promote arrestin interaction. There are seven different GRK subtypes, which may uniquely phosphorylate the C-terminal tail in a type of 'phosphorylation barcode,' potentially differentially contributing to arrestin translocation and arrestin-dependent signalling. Such contributions may be exploited to develop arrestin-biased ligands. Here, we examine the effect of different GRK subtypes on the ability to promote translocation of arrestin-2 and arrestin-3 to the cannabinoid CB1 receptor (CB1) with a range of ligands. We find that most GRK subtypes (including visual GRK1) can enhance arrestin-2 and -3 translocation to CB1, and that GRK-dependent changes in arrestin-2 and arrestin-3 translocation were broadly shared for most agonists tested. GRK2/3 generally enhanced arrestin translocation more than the other GRK subtypes, with some small differences between ligands. We also explore the interplay between G protein activity and GRK2/3-dependent arrestin translocation, highlighting that high-efficacy G protein agonists will cause GRK2/3 dependent arrestin translocation. This study supports the hypothesis that arrestin-biased ligands for CB1 must engage GRK5/6 rather than GRK2/3, and G protein-biased ligands must have inherently low efficacy.

Bias translation: The final frontier?

Biased signalling is a natural result of GPCR allosteric function and should be expected from any and all synthetic and natural agonists. Therefore, it may be encountered in all agonist discovery projects and must be considered as a beneficial (or possible detrimental) feature of new candidate molecules. While bias is detected easily, the synoptic nature of GPCR signalling makes translation of simple in vitro bias to complex in vivo systems problematic. The practical outcome of this is a difficulty in predicting the therapeutic value of biased signalling due to the failure of translation of identified biased signalling to in vivo agonism. This is discussed in this review as well as some new ways forward to improve this translation process and better exploit this powerful pharmacologic mechanism.

Frizzleds act as dynamic pharmacological entities

The Frizzled family of transmembrane receptors (FZD1-10) belongs to the class F of G protein-coupled receptors (GPCRs). FZDs bind to and are activated by Wingless/Int1 (WNT) proteins. The WNT/FZD signaling system regulates crucial aspects of developmental biology and stem-cell regulation. Dysregulation of WNT/FZD communication can lead to developmental defects and diseases such as cancer and fibrosis. Recent insight into the activation mechanisms of FZDs has underlined that protein dynamics and conserved microswitches are essential for FZD-mediated information flow and build the basis for targeting these receptors pharmacologically. In this review, we summarize recent advances in our understanding of FZD activation, and how novel concepts merge and collide with existing dogmas in the field.

Understanding the effects of oxytocin receptor variants on OXT-OXT receptor binding: A mathematical model

Approximately half of U.S. women giving birth annually receive Pitocin, the synthetic form of oxytocin (OXT), yet its effective dose can vary significantly. This variability presents safety concerns due to unpredictable responses, which may lead to adverse outcomes for both mother and baby. To address the need for improved dosing, we developed a data-driven mathematical model to predict OXT receptor (OXTR) binding. Our study focuses on five prevalent OXTR variants (V45L, P108A, L206V, V281M, and E339K) and their impact on OXT-OXTR binding dynamics in two distinct cell types: human embryonic kidney cells (HEK293T), commonly used in experimental systems, and human myometrial smooth muscle cells, containing endogenous OXTR. We parameterized the model with cell-specific OXTR surface localization measurements. To strengthen the robustness of our study, we conducted a comprehensive meta-analysis of OXT- OXTR binding, enabling parameterization of our model with cell-specific OXT-OXTR binding kinetics (myometrial OXT-OXTR K d = 1.6 nM, kon = 6.8 × 10 5 M -1 min -1 , and koff = 0.0011 min -1 ). Our meta-analysis revealed significant homogeneity in OXT-OXTR affinity across experiments and species with a K d = 0.52 - 9.32 nM and mean K d = 1.48 ± 0.36 nM. Our model achieves several valuable insights into designing dosage strategies. First, we predicted that the OXTR complex reaches maximum occupancy at 10 nM OXT in myometrial cells and at 1 µM in HEK293T cells. This information is pivotal for guiding experimental design and data interpretation when working with these distinct cell types, emphasizing the need to consider effects for specific cell types when choosing OXTR-transfected cell lines. Second, our model recapitulated the significant effects of genetic variants for both experimental and physiologically relevant systems, with V281M and E339K substantially compromising OXT-OXTR binding capacity. These findings suggest the need for personalized oxytocin dosing based on individual genetic profiles to enhance therapeutic efficacy and reduce risks, especially in the context of labor and delivery. Third, we demonstrated the potential for rescuing the attenuated cell response observed in V281M and E339K variants by increasing the OXT dosage at specific, early time points. Cellular responses to OXT, including Ca 2+ release, manifest within minutes. Our model indicates that providing V281M- and E339K-expressing cells with doubled OXT dose during the initial minute of binding can elevate OXT-OXTR complex formation to levels comparable to wild-type OXTR. In summary, our study provides a computational framework for precision oxytocin dosing strategies, paving the way for personalized medicine.

An Algebraic Blueprint for Predicting Turnover Numbers and Endpoints in Photocatalysis

Photocatalysis is a contemporary research field given that the world's fossil energy resources including coal, mineral oil and natural gas are finite. The vast variety of photocatalytic systems demands for standardized protocols facilitating an objective comparison. While there are commonly accepted performance indicators such as the turnover number (TON) that are usually reported, to date there is no unified concept for the determination of TONs and the endpoint of the reaction during continuous measurements. Herein, we propose an algebraic approach using defined parameters and boundary conditions based on partial-least squares regression for generically calculating and predicting the turnover number and the endpoint of a photocatalytic experiment. Furthermore, the impact of the analysis period was evaluated with respect to the fidelity of the obtained TON, and the influence of the data point density along critical segments of the obtained fitting function is demonstrated.

Signaling Specificity and Kinetics of the Human Metabotropic Glutamate Receptors

Metabotropic glutamate receptors (mGluRs) are obligate dimer G protein coupled receptors that can all function as homodimers. Here, each mGluR homodimer was examined for its G protein coupling profile using a bioluminescence resonance energy transfer-based assay that detects the interaction between a split YFP-tagged Gβ 1γ2 and a Nanoluciferase tagged free Gβγ sensor, MAS-GRK3-ct- nanoluciferase with 14 specific Gα proteins heterologously expressed, representing each family. Canonically, the group II and III mGluRs (2 and 3 and 4, 6, 7, and 8, respectively) are thought to couple to Gi/o exclusively. In addition, the group I mGluRs (1 and 5) are known to couple to the Gq/11 family and generally thought to also couple to the pertussis toxin-sensitive Gi/o family some reports have suggested Gs coupling is possible as cAMP elevations have been noted. In this study, coupling was observed with all eight mGluRs through the Gi/o proteins and only mGluR1 and mGluR5 through Gq/11, and, perhaps surprisingly, not G14 None activated any Gs protein. Interestingly, coupling was seen with the group I and II but not the group III mGluRs to G16 Slow but significant coupling to Gz was also seen with the group II receptors. SIGNIFICANCE STATEMENT: Metabotropic glutamate receptor (mGluR)-G protein coupling has not been thoroughly examined, and some controversy remains about whether some mGluRs can activate Gαs family members. Here we examine the ability of each mGluR to activate representative members of every Gα protein family. While all mGluRs can activate Gαi/o proteins, only the group I mGluRs couple to Gαq/11, and no members of the family can activate Gαs family members, including the group I receptors alone or with positive allosteric modulators.

ClickArr: a novel, high-throughput assay for evaluating β-arrestin isoform recruitment

Background: Modern methods for quantifying signaling bias at G protein-coupled receptors (GPCRs) rely on using a single β-arrestin isoform. However, it is increasingly appreciated that the two β-arrestin isoforms have unique roles, requiring the ability to assess β-arrestin isoform preference. Thus, methods are needed to efficiently screen the recruitment of both β-arrestin isoforms as they compete for a target GPCR in cells. Methods: We used molecular cloning to develop fusion proteins of the δ-opioid receptor (δOR), β-arrestin 1, and β-arrestin 2 to fragments of click beetle green and click beetle red luciferases. In this assay architecture, recruitment of either β-arrestin 1 or 2 to the δOR generates a spectrally distinct bioluminescent signal, allowing us to co-transfect all three constructs into cells prior to agonist challenge. Results: We demonstrate that our new assay, named "ClickArr," is a live-cell assay that simultaneously reports the recruitment of both β-arrestin isoforms as they compete for interaction with the δOR. We further find that the partial δOR agonist TAN67 has a significant efficacy bias for β-arrestin 2 over β-arrestin 1 when recruitment is normalized to the reference agonist leu-enkephalin. We confirm that ClickArr reports this bias when run either as a high-throughput endpoint or high-throughput kinetic assay, and cross-validate this result using the PathHunter assay, an orthogonal commercial assay for reporting β-arrestin recruitment to the δOR. Conclusion: Our results suggest that agonist:GPCR complexes can have relative β-arrestin isoform bias, a novel signaling bias that may potentially open up a new dimension for drug development.

REGA-SIGN: Development of a Novel Set of NanoBRET-Based G Protein Biosensors

Despite G protein-coupled receptors (GPCRs) being important theapeutic targets, the signaling properties of many GPCRs remain poorly characterized. GPCR activation primarily initiates heterotrimeric G protein signaling. To detect ligand-induced G protein activation, Bioluminescence Resonance Energy Transfer (BRET)-based biosensors were previously developed. Here, we designed a novel set of Nanoluciferase (NLuc) BRET-based biosensors (REGA-SIGN) that covers all Gα protein families (i.e., Gαi/o, GαSs/L, Gα12/13 and Gαq/15). REGA-SIGN uses NLuc as a bioluminescent donor and LSS-mKATE2, a red-shifted fluorophore, as an acceptor. Due to the enhanced spectral separation between donor and acceptor emission and the availability of a stable substrate for NLuc, this donor-acceptor pair enables sensitive kinetic assessment of G protein activity. After optimization, the NLuc integration sites into the Gα subunit largely corresponded with previously reported integration sites, except for GαSs/L for which we describe an alternative NLuc insertion site. G protein rescue experiments validated the biological activity of these Gα donor proteins. Direct comparison between EGFP and LSS-mKATE2 as acceptor fluorophores revealed improved sensitivity for nearly all G protein subtypes when using the latter one. Hence, REGA-SIGN can be used as a panel of kinetic G protein biosensors with high sensitivity.

Signaling specificity and kinetics of the human metabotropic glutamate receptors

Metabotropic glutamate receptors (mGluRs) are obligate dimer G protein coupled receptors that can all function as homodimers. Here, each mGluR homodimer was examined for its G protein coupling profile using a BRET based assay that detects the interaction between a split YFP-tagged Gβ1γ2 and a Nanoluc tagged free Gβγ sensor, MAS-GRK3-ct-NLuc with 14 specific Ga proteins heterologously expressed, representing each family. Canonically, the group II and III mGluRs (2&3, and 4, 6, 7&8, respectively) are thought to couple to Gi/o exclusively. In addition, the group I mGluRs (1&5) are known to couple to the Gq/11 family, and generally thought to also couple to the PTX-sensitive Gi/o family; some reports have suggested Gs coupling is possible as cAMP elevations have been noted. In this study, coupling was observed with all 8 mGluRs through the Gi/o proteins, and only mGluR1&5 through Gq/11, and perhaps surprisingly, not G14. None activated any Gs protein. Interestingly, coupling was seen with the group I and II, but not the group III mGluRs to G16. Slow but significant coupling to Gz was also seen with the group II receptors.

Investigation of adenosine A1 receptor-mediated β-arrestin 2 recruitment using a split-luciferase assay

Background: Adenosine A1 receptor (A₁AR) plays a prominent role in neurological and cardiac diseases and inflammatory processes. Its endogenous ligand adenosine is known to be one of the key players in the sleep-wake cycle. Like other G protein-coupled receptors (GPCRs), stimulation of A₁AR leads to the recruitment of arrestins in addition to the activation of G proteins. So far, little is known about the role of these proteins in signal transduction and regulation of A₁AR compared to the activation of G proteins. In this work, we characterized a live cell assay for A₁AR-mediated β-arrestin 2 recruitment. We have applied this assay to a set of different compounds that interact with this receptor. Methods: Based on NanoBit^® technology, a protein complementation assay was developed in which the A₁AR is coupled to the large part of the nanoluciferase (LgBiT), whereas its small part (SmBiT) is fused to the N-terminus of β-arrestin 2. Stimulation of A₁AR results in the recruitment of β-arrestin 2 and subsequent complementation of a functional nanoluciferase. For comparison, corresponding data on the effect of receptor stimulation on intracellular cAMP levels were collected for some data sets using the GloSensor™ assay. Results: The assay gives highly reproducible results with a very good signal-to-noise ratio. Capadenoson, in contrast to adenosine, CPA, or NECA, shows only partial agonism in this assay with respect to the recruitment of β-arrestin 2, whereas it shows full agonism in the case of the inhibitory effect of A₁AR on cAMP production. By using a GRK2 inhibitor, it becomes clear that the recruitment is at least partially dependent on the phosphorylation of the receptor by this kinase. Interestingly, this was also the first time that we demonstrate the A₁AR-mediated recruitment of β-arrestin 2 by stimulation with a valerian extract. Conclusion: The presented assay is a useful tool for the quantitative study of A₁AR-mediated β-arrestin 2 recruitment. It allows data collection for stimulatory, inhibitory, and modulatory substances and is also suitable for more complex substance mixtures such as valerian extract.

Kinetic analysis of endogenous β2 -adrenoceptor-mediated cAMP GloSensor™ responses in HEK293 cells

BACKGROUND AND AIM

Standard pharmacological analysis of agonist activity utilises measurements of receptor-mediated responses at a set time-point, or at the peak response level, to characterise ligands. However, the occurrence of non-equilibrium conditions may dramatically impact the properties of the response being measured. Here we have analysed the initial kinetic phases of cAMP responses to β2 -adrenoceptor agonists in HEK293 cells expressing the endogenous β2 -adrenoceptor at extremely low levels.

EXPERIMENTAL APPROACH

The kinetics of β2 -adrenoceptor agonist-stimulated cAMP responses were monitored in real-time, in the presence and absence of antagonists, in HEK293 cells expressing the cAMP GloSensor™ biosensor. Potency (EC50 ) and efficacy (Emax ) values were determined at the peak of the agonist GloSensor™ response and compared to kinetic parameters L50 and IRmax values derived from initial response rates.

KEY RESULTS

The partial agonists salbutamol and salmeterol displayed reduced relative IRmax values (with respect to isoprenaline) when compared with their Emax values. Except for the fast dissociating bisoprolol, preincubation with β2 -adrenoceptor antagonists produced a large reduction in the isoprenaline peak response due to a state of hemi-equilibrium in this low receptor reserve system. This effect was exacerbated when IRmax parameters were measured. Furthermore, bisoprolol produced a large reduction in isoprenaline IRmax consistent with its short residence time.

CONCLUSIONS AND IMPLICATIONS

Kinetic analysis of real-time signalling data can provide valuable insights into the hemi-equilibria that can occur in low receptor reserve systems with agonist-antagonist interactions, due to incomplete dissociation of antagonist whilst the peak agonist response is developing.

Cannabinoid 1 (CB1 ) receptor arrestin subtype-selectivity and phosphorylation dependence

BACKGROUND AND PURPOSE

Arrestin or G protein bias may be desirable for novel cannabinoid therapeutics. Arrestin-2 and arrestin-3 translocation to CB₁ receptor have been suggested to mediate different functions that may be exploited with biased ligands. Here, the requirement of a recently described phosphorylation motif 'pxxp' (where 'p' denotes phosphorylatable serine or threonine and 'x' denotes any other amino acid) within the CB₁ receptor C-terminus for interaction with different arrestin subtypes was examined.

EXPERIMENTAL APPROACH

Site-directed mutagenesis was conducted to generate nine different phosphorylation-impaired CB₁ receptor C-terminal mutants. Bioluminescence resonance energy transfer (BRET) was employed to measure arrestin-2/3 translocation and G protein dissociation of a high efficacy agonist for each mutant. Immunocytochemistry was used to quantify receptor expression.

KEY RESULTS

The effects of each mutation were shared for arrestin-2 and arrestin-3 translocation to CB₁ receptor pxxp motifs are partially required for arrestin-2/3 translocation, but translocation was not completely inhibited until all phosphorylation sites were mutated. The rate of arrestin translocation was reduced with simultaneous mutation of S425 and S429. Desensitisation of G protein dissociation was inhibited in different mutants proportional to the extent of their respective loss of arrestin translocation.

CONCLUSIONS AND IMPLICATIONS

These data do not support the existence of an 'essential' pxxp motif for arrestin translocation to CB₁ receptor. These data also identify that arrestin-2 and arrestin-3 have equivalent phosphorylation requirements within the CB₁ receptor C-terminus, suggesting arrestin subtype-selective biased ligands may not be viable and that different regions of the C-terminus contribute differently to arrestin translocation.

Thermodynamic architecture and conformational plasticity of GPCRs

G-protein-coupled receptors (GPCRs) are ubiquitous integral membrane proteins involved in diverse cellular signaling processes. Here, we carry out a large-scale ensemble thermodynamic study of 45 ligand-free GPCRs employing a structure-based statistical mechanical framework. We find that multiple partially structured states co-exist in the GPCR native ensemble, with the TM helices 1, 6 and 7 displaying varied folding status, and shaping the conformational landscape. Strongly coupled residues are anisotropically distributed, accounting for only 13% of the residues, illustrating that a large number of residues are inherently dynamic. Active-state GPCRs are characterized by reduced conformational heterogeneity with altered coupling-patterns distributed throughout the structural scaffold. In silico alanine-scanning mutagenesis reveals that extra- and intra-cellular faces of GPCRs are coupled thermodynamically, highlighting an exquisite structural specialization and the fluid nature of the intramolecular interaction network. The ensemble-based perturbation methodology presented here lays the foundation for understanding allosteric mechanisms and the effects of disease-causing mutations in GCPRs.

Development of fluorescent peptide G protein-coupled receptor activation biosensors for NanoBRET characterization of intracellular allosteric modulators

G protein-coupled receptors (GPCRs) are widely therapeutically targeted, and recent advances in allosteric modulator development at these receptors offer further potential for exploitation. Intracellular allosteric modulators (IAM) represent a class of ligands that bind to the receptor-effector interface (e.g., G protein) and inhibit agonist responses noncompetitively. This potentially offers greater selectivity between receptor subtypes compared to classical orthosteric ligands. However, while examples of IAM ligands are well described, a more general methodology for assessing compound interactions at the IAM site is lacking. Here, fluorescent labeled peptides based on the Gα peptide C terminus are developed as novel binding and activation biosensors for the GPCR-IAM site. In TR-FRET binding studies, unlabeled peptides derived from the Gαs subunit were first characterized for their ability to positively modulate agonist affinity at the β₂ -adrenoceptor. On this basis, a tetramethylrhodamine (TMR) labeled tracer was synthesized based on the 19 amino acid Gαs peptide (TMR-Gαs19cha18, where cha = cyclohexylalanine). Using NanoBRET technology to detect binding, TMR-Gαs19cha18 was recruited to Gs coupled β₂ -adrenoceptor and EP₂ receptors in an agonist-dependent manner, but not the Gi-coupled CXCR2 receptor. Moreover, NanoBRET competition binding assays using TMR-Gαs19cha18 enabled direct assessment of the affinity of unlabeled ligands for β₂ -adrenoceptor IAM site. Thus, the NanoBRET platform using fluorescent-labeled G protein peptide mimetics offers novel potential for medium-throughput screens to identify IAMs, applicable across GPCRs coupled to a G protein class. Using the same platform, Gs peptide biosensors also represent useful tools to probe orthosteric agonist efficacy and the dynamics of receptor activation.

The Impact of the Secondary Binding Pocket on the Pharmacology of Class A GPCRs

G-protein coupled receptors (GPCRs) are considered important therapeutic targets due to their pathophysiological significance and pharmacological relevance. Class A receptors represent the largest group of GPCRs that gives the highest number of validated drug targets. Endogenous ligands bind to the orthosteric binding pocket (OBP) embedded in the intrahelical space of the receptor. During the last 10 years, however, it has been turned out that in many receptors there is secondary binding pocket (SBP) located in the extracellular vestibule that is much less conserved. In some cases, it serves as a stable allosteric site harbouring allosteric ligands that modulate the pharmacology of orthosteric binders. In other cases it is used by bitopic compounds occupying both the OBP and SBP. In these terms, SBP binding moieties might influence the pharmacology of the bitopic ligands. Together with others, our research group showed that SBP binders contribute significantly to the affinity, selectivity, functional activity, functional selectivity and binding kinetics of bitopic ligands. Based on these observations we developed a structure-based protocol for designing bitopic compounds with desired pharmacological profile.

Quantifying the Kinetics of Signaling and Arrestin Recruitment by Nervous System G-Protein Coupled Receptors

Neurons integrate inputs over different time and space scales. Fast excitatory synapses at boutons (ms and μm), and slow modulation over entire dendritic arbors (seconds and mm) are all ultimately combined to produce behavior. Understanding the timing of signaling events mediated by G-protein-coupled receptors is necessary to elucidate the mechanism of action of therapeutics targeting the nervous system. Measuring signaling kinetics in live cells has been transformed by the adoption of fluorescent biosensors and dyes that convert biological signals into optical signals that are conveniently recorded by microscopic imaging or by fluorescence plate readers. Quantifying the timing of signaling has now become routine with the application of equations in familiar curve fitting software to estimate the rates of signaling from the waveform. Here we describe examples of the application of these methods, including (1) Kinetic analysis of opioid signaling dynamics and partial agonism measured using cAMP and arrestin biosensors; (2) Quantifying the signaling activity of illicit synthetic cannabinoid receptor agonists measured using a fluorescent membrane potential dye; (3) Demonstration of multiplicity of arrestin functions from analysis of biosensor waveforms and quantification of the rates of these processes. These examples show how temporal analysis provides additional dimensions to enhance the understanding of GPCR signaling and therapeutic mechanisms in the nervous system.

Metadynamics simulations leveraged by statistical analyses and artificial intelligence-based tools to inform the discovery of G protein-coupled receptor ligands

G Protein-Coupled Receptors (GPCRs) are a large family of membrane proteins with pluridimensional signaling profiles. They undergo ligand-specific conformational changes, which in turn lead to the differential activation of intracellular signaling proteins and the consequent triggering of a variety of biological responses. This conformational plasticity directly impacts our understanding of GPCR signaling and therapeutic implications, as do ligand-specific kinetic differences in GPCR-induced transducer activation/coupling or GPCR-transducer complex stability. High-resolution experimental structures of ligand-bound GPCRs in the presence or absence of interacting transducers provide important, yet limited, insights into the highly dynamic process of ligand-induced activation or inhibition of these receptors. We and others have complemented these studies with computational strategies aimed at characterizing increasingly accurate metastable conformations of GPCRs using a combination of metadynamics simulations, state-of-the-art algorithms for statistical analyses of simulation data, and artificial intelligence-based tools. This minireview provides an overview of these approaches as well as lessons learned from them towards the identification of conformational states that may be difficult or even impossible to characterize experimentally and yet important to discover new GPCR ligands.

Prior Acoustic Trauma Alters Type II Afferent Activity in the Mouse Cochlea

Auditory stimuli travel from the cochlea to the brainstem through type I and type II cochlear afferents. While type I afferents convey information about the frequency, intensity, and timing of sounds, the role of type II afferents remains unresolved. Limited recordings of type II afferents from cochlear apex of prehearing rats reveal they are activated by widespread outer hair cell stimulation, ATP, and by the rupture of nearby outer hair cells. Altogether, these lines of evidence suggest that type II afferents sense loud, potentially damaging levels of sound. To explore this hypothesis further, calcium imaging was used to determine the impact of acoustic trauma on the activity of type II cochlear afferents of young adult mice of both sexes. Two known marker genes (Th, Drd2) and one new marker gene (Tac1), expressed in type II afferents and some other cochlear cell types, drove GCaMP6f expression to reveal calcium transients in response to focal damage in the organ of Corti in all turns of the cochlea. Mature type II afferents responded to acute photoablation damage less often but at greater length compared with prehearing neurons. In addition, days after acoustic trauma, acute photoablation triggered a novel response pattern in type II afferents and surrounding epithelial cells, delayed bursts of activity occurring minutes after the initial response subsided. Overall, calcium imaging can report type II afferent responses to damage even in mature and noise-exposed animals and reveals previously unknown tissue hyperactivity subsequent to acoustic trauma.

Circuit-specific enteric glia regulate intestinal motor neurocircuits

Glia in the central nervous system exert precise spatial and temporal regulation over neural circuitry on a synapse-specific basis, but it is unclear if peripheral glia share this exquisite capacity to sense and modulate circuit activity. In the enteric nervous system (ENS), glia control gastrointestinal motility through bidirectional communication with surrounding neurons. We combined glial chemogenetics with genetically encoded calcium indicators expressed in enteric neurons and glia to study network-level activity in the intact myenteric plexus of the proximal colon. Stimulation of neural fiber tracts projecting in aboral, oral, and circumferential directions activated distinct populations of enteric glia. The majority of glia responded to both oral and aboral stimulation and circumferential pathways, while smaller subpopulations were activated only by ascending and descending pathways. Cholinergic signaling functionally specifies glia to the descending circuitry, and this network plays an important role in repressing the activity of descending neural pathways, with some degree of cross-inhibition imposed upon the ascending pathway. Glial recruitment by purinergic signaling functions to enhance activity within ascending circuit pathways and constrain activity within descending networks. Pharmacological manipulation of glial purinergic and cholinergic signaling differentially altered neuronal responses in these circuits in a sex-dependent manner. Collectively, our findings establish that the balance between purinergic and cholinergic signaling may differentially control specific circuit activity through selective signaling between networks of enteric neurons and glia. Thus, enteric glia regulate the ENS circuitry in a network-specific manner, providing profound insights into the functional breadth and versatility of peripheral glia.

High-Content Single-Cell Förster Resonance Energy Transfer Imaging of Cultured Striatal Neurons Reveals Novel Cross-Talk in the Regulation of Nuclear Signaling by Protein Kinase A and Extracellular Signal-Regulated Kinase 1/2

Genetically encoded biosensors can be used to track signaling events in living cells by measuring changes in fluorescence emitted by one or more fluorescent proteins. Here, we describe the use of genetically encoded biosensors based on Förster resonance energy transfer (FRET), combined with high-content microscopy, to image dynamic signaling events simultaneously in thousands of neurons in response to drug treatments. We first applied this approach to examine intercellular variation in signaling responses among cultured striatal neurons stimulated with multiple drugs. Using high-content FRET imaging and immunofluorescence, we identified neuronal subpopulations with unique responses to pharmacological manipulation and used nuclear morphology to identify medium spiny neurons within these heterogeneous striatal cultures. Focusing on protein kinase A (PKA) and extracellular signal-regulated kinase 1/2 (ERK1/2) signaling in the cytoplasm and nucleus, we noted pronounced intercellular differences among putative medium spiny neurons, in both the magnitude and kinetics of signaling responses to drug application. Importantly, a conventional "bulk" analysis that pooled all cells in culture yielded a different rank order of drug potency than that revealed by single-cell analysis. Using a single-cell analytical approach, we dissected the relative contributions of PKA and ERK1/2 signaling in striatal neurons and unexpectedly identified a novel role for ERK1/2 in promoting nuclear activation of PKA in striatal neurons. This finding adds a new dimension of signaling crosstalk between PKA and ERK1/2 with relevance to dopamine D1 receptor signaling in striatal neurons. In conclusion, high-content single-cell imaging can complement and extend traditional population-level analyses and provides a novel vantage point from which to study cellular signaling. SIGNIFICANCE STATEMENT: High-content imaging revealed substantial intercellular variation in the magnitude and pattern of intracellular signaling events driven by receptor stimulation. Since individual neurons within the same population can respond differently to a given agonist, interpreting measures of intracellular signaling derived from the averaged response of entire neuronal populations may not always reflect what happened at the single-cell level. This study uses this approach to identify a new form of cross-talk between PKA and ERK1/2 signaling in the nucleus of striatal neurons.

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.

Conserved residues in the extracellular loop 2 regulate Stachel-mediated activation of ADGRG2

Cleavage and dissociation of a large N-terminal fragment and the consequent unmasking of a short sequence (Stachel) remaining on the N-terminus have been proposed as mechanisms of activation of some members of the adhesion G protein-coupled receptor (aGPCR) family. However, the identity of residues that play a role in the activation of aGPCRs by the cognate Stachel remains largely unknown. Protein sequence alignments revealed a conserved stretch of residues in the extracellular loop 2 (ECL2) of all 33 members of the aGPCR family. ADGRG2, an orphan aGPCR, plays a major role in male fertility, Ewing sarcoma cell proliferation, and parathyroid cell function. We used ADGRG2 as a model aGPCR and generated mutants of the conserved residues in the ECL2 via site-directed mutagenesis. We show that tryptophan and isoleucine in the ECL2 are essential for receptor stability and surface expression in the HEK293 cells. By adjusting the receptor surface expression levels, we show that mutation of these residues of ECL2 ablates the Stachel-mediated activation of multiple signaling pathways of ADGRG2. This study provides a novel understanding of the role of the ECL2 in Stachel-mediated signaling and degradation of ADGRG2, which may lay the foundation for the rational design of therapeutics to target aGPCRs.