Effector membrane translocation biosensors reveal G protein and βarrestin coupling profiles of 100 therapeutically relevant GPCRs
Abstract
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.
Editor's evaluation
The authors describe a novel set of biosensors to assess the coupling specificity of 100 therapeutically relevant G proteins-coupled receptors (GPCRs) to various G proteins. The utility of the assay system is well-supported by the data. These tools are likely to be useful for many specific studies of individual receptors, including efforts to discover ligands that display functional selectivity (bias) between G protein pathways or between G proteins and arrestins. The work provides a rich repository of data informing on the possible effector coupling of 100 GPCRs and a set of analytical tools that could guide the development of new drugs, including efforts to discover ligands that display functional selectivity (bias) between G protein pathways or between G proteins and arrestins.
https://doi.org/10.7554/eLife.74101.sa0Introduction
G protein-coupled receptors (GPCRs) play crucial roles in the regulation of a wide variety of physiological processes and represent one-third of clinically prescribed drugs (Hauser et al., 2017). Classically, GPCR-mediated signal transduction was believed to rely on linear signaling pathways whereby a given GPCR selectively activates a single G protein family, defined by the nature of its Gα subunit (Oldham and Hamm, 2008). Gα proteins are divided into four major families (Gs, Gi/o, Gq/11, and G12/13) encoded by 16 human genes. Once activated, these proteins each trigger different downstream effectors yielding different biological outcomes. It has now become evident that many GPCRs can couple to more than one G protein family and that ligands can selectively promote the activation of different subsets of these pathways (Namkung et al., 2018; Quoyer et al., 2013). These observations extended the concept of ligand-biased signaling, which was first established for ligand-directed selectivity between βarrestin and G protein (Azzi et al., 2003; Wei et al., 2003), to functional selectivity between G proteins. Ligand-directed functional selectivity represents a promising avenue for GPCRs drug discovery since it offers the opportunity of activating pathways important for therapeutic efficacy while minimizing activation of pathways responsible for undesirable side effects (Galandrin et al., 2007; Kenakin, 2019).
To fully explore the potential of functional selectivity, it is essential to have an exhaustive description of the signaling partners that can be activated by a given receptor, providing receptor- and ligand-specific signaling signatures. Currently, few assays allow for an exhaustive pathway-specific analysis of GPCR signaling; these include BRET-based G protein activation sensors platforms (Galés et al., 2005; Masuho et al., 2015; Maziarz et al., 2020; Mende et al., 2018; Olsen et al., 2020) and the TGF-α shedding assay (Inoue et al., 2019). However, several of these platforms require modification of G protein subunits that may create functional distortions. Moreover, these assays may detect non-productive conformational rearrangements of the G protein heterotrimer as was recently reported for G12 (Okashah et al., 2020).
Here, we describe unique sensors that do not require modification of receptors or G proteins (except for Gs) for interrogating the signaling profiles of GPCRs. The platform includes 15 pathway-selective enhanced bystander bioluminescence resonance energy transfer (ebBRET) biosensors monitoring the translocation of downstream effectors to the plasma membrane for Gi/o, Gq/11, and G12/13, the dissociation of the Gα subunit from the plasma membrane for Gs and the recruitment of βarrestin to the plasma membrane. Overall, the new ebBRET-based Effector Membrane Translocation Assays, named EMTA, provide a readily accessible large scale and comprehensive platform to study constitutive and ligand-directed GPCR signaling. The signaling signatures of 100 GPCRs using the EMTA platform also provides a rich source of information to explore the principles underlying receptor/G protein/βarrestin coupling selectivity relationships. It thus provides a unique set of tools that is complementary to previously described platforms and existing datasets, and offers a map of the coupling potentials for individual GPCR that will stimulate future studies investigating the relevance of these couplings in different physiological systems.
Results
ebBRET-based G protein effector membrane translocation assay (EMTA) allows detection of each Gα protein subunit activation
To detect the activation of Gα subtypes, we created an EMTA biosensor platform based on ebBRET (Namkung et al., 2016; Figure 1A). The biosensors at the heart of EMTA consist of sub-domains of the G protein-effector proteins p63-RhoGEF, Rap1GAP and PDZ-RhoGEF that selectively interact with activated Gq/11, Gi/o, or G12/13, respectively. These domains were fused at their C-terminus to Renilla luciferase (RlucII) and co-expressed with different unmodified receptor and Gα protein subtypes. Upon GPCR activation, the energy donor-fused effectors translocate to the plasma membrane to bind activated Gα proteins, bringing RlucII in close proximity to the energy acceptor, Renilla green fluorescent protein, targeted to the plasma membrane through a CAAX motif (rGFP-CAAX), thus leading to an increase in ebBRET. The same plasma membrane translocation principle is used to measure βarrestin recruitment (Namkung et al., 2016; Figure 1B, top). Because no selective soluble downstream effector of Gs exists, the assay was modified taking advantage of Gαs dissociation from the plasma membrane following its activation (Wedegaertner et al., 1996). In this configuration, the RlucII is directly fused to Gαs (Carr et al., 2014). Its activation upon GPCR stimulation leads to its dissociation from the plasma membrane (Martin and Lambert, 2016), resulting in a reduction in ebBRET (Figure 1B, bottom).
The sensitivity and selectivity of the newly created G protein EMTA biosensors, were validated using prototypical GPCRs known to activate specific Gα subtypes. The responses were monitored upon heterologous expression of specific Gα subunits belonging to Gi/o, Gq/11, or G12/13 families in the absence or presence of pharmacological inhibitors and using engineered cells lacking selected Gα subtypes. The dopamine D2 receptor was used to validate the ability of the Gi/o binding domain of Rap1GAP (Jordan et al., 1999; Meng et al., 1999) to selectively detect Gi/o activation. The dopamine-promoted increase in ebBRET between Rap1GAP-RlucII and rGFP-CAAX in the presence of Gαi/o subunits was not affected by the Gq/11-selective inhibitor UBO-QIC (a.k.a., FR900359; Schrage et al., 2015; Figure 2A, left), whereas the Gαi/o family inhibitor, pertussis toxin (PTX), completely blocked the response for all members of Gαi/o family except for Gαz, known to be insensitive to PTX (Casey et al., 1990; Figure 2A, right). Gonadotropin-releasing hormone (GnRH) stimulation of the GnRH receptor (GnRHR), used as a prototypical Gq/11-coupled receptor, promoted ebBRET between the RlucII-fused Gq/11 binding domain of p63-RhoGEF (p63-RhoGEF-RlucII; Lutz et al., 2007; Rojas et al., 2007) and rGFP-CAAX. The ebBRET increase observed in the presence of different Gαq/11 subunits was not significantly (p = 0.077, 0.0636 and 0.073 for Gq, G11, and G14, respectively) affected by PTX (Figure 2B, right), whereas UBO-QIC completely blocked the response for all members of Gαq/11 family except for Gα15, known to be insensitive to UBO-QIC (Schrage et al., 2015; Figure 2B, left). These two G protein-specific EMTA were sensitive enough to detect responses elicited by endogenous G proteins since deletion of Gi/o (ΔGi/o) or Gq/11 (ΔGq/11) subtypes completely abolished the responses induced by D2 or GnRHR activation in the absence of heterologously expressed G proteins (Figure 2—figure supplement 1I). It should however be noted that relying on endogenous proteins does not allow the identification of specific members of Gi/o (i.e.: Gi1, Gi2, Gi3, GoA, GoB, or Gz) or Gq/11 (i.e.: Gq, G11, G14, or G15) families.
The selectivity of the G12/13 binding domain of PDZ-RhoGEF (Fukuhara et al., 2001) was confirmed using the cannabinoid receptor type 1 (CB1). The ebBRET between PDZ-RhoGEF-RlucII and rGFP-CAAX in the presence of Gα12 or Gα13 promoted by the cannabinoid agonist WIN-55,212–2 was not affected by UBO-QIC (Figure 2C, top left), nor PTX (Figure 2C, top right). Given the lack of selective G12/13 pharmacological inhibitor, we used HEK293 cells genetically deleted for Gα12 and Gα13 proteins (ΔG12/13) to further confirm the response selectivity. As expected, PDZ-RhoGEF-RlucII/rGFP-CAAX ebBRET was observed only following reintroduction of either Gα12 (ΔG12/13_+G12) or Gα13 (ΔG12/13_+G13) (Figure 2C, bottom left). The G12/13 coupling of CB1 was further confirmed by monitoring the recruitment of PKN to the plasma membrane (Figure 2C, bottom right) in agreement with previous reports (Inoue et al., 2019).
To further assess the selectivity of each EMTA biosensor, we took advantage of the fact that the endothelin-1 receptor (ETA) can activate Gq/11, Gi/o, and G12/13 family members. As shown in Figure 2—figure supplement 2, only over-expression of the Gα family members corresponding to their selective effectors (Rap1GAP for Gi/o, p63-RhoGEF for Gq/11, and PDZ-RhoGEF for G12/13) significantly increased the recruitment of the effector-RlucII to the plasma membrane. A recent study (Chandan et al., 2021) showed that Gi/o can also activate full length PDZ-RhoGEF. Although the domain of PDZ-RhoGEF required for this activation has not been identified yet, the selectivity of our PDZ-RhoGEF sensor for G12/13 vs. all other G protein families most likely results from the fact that we used a truncated version of PDZ-RhoGEF that only contains the G12/13 binding domain and lacks the PDZ domain involved in protein-protein interaction, the actin-binding domain and the DH/PH domains involved in GEF activity and RhoA activation (Aittaleb et al., 2010).
It should be noted that in the heterologous expression configuration, competition with endogenous G proteins did not occur to a significant extent since the potencies of the responses to a given G protein subtype were not affected by genetic deletion of the different G protein family members (Figure 2—figure supplement 1 and Supplementary file 1A). Similarly, overexpression of G proteins, GPCRs or effectors-RlucII did not affect the potencies of the responses observed (Figure 2—figure supplement 3 and Supplementary file 1B-D), indicating that, in our experimental conditions, overexpression of the different components of EMTA sensors must likely not bias the coupling response. In addition to spectrometric assessment of coupling selectivity (above) and activation kinetics (Figure 2—figure supplement 4), EMTA allows to image the real-time recruitment of the G protein effectors to the plasma membrane (Videos 1–3) thus providing spatiotemporal resolution for the imaging detection of Gαi/o, Gαq/11, and Gα12/13 activation.
The sensitivity of the EMTA platform is illustrated by a direct side-by-side comparison of the signals detected with EMTA vs. BRET assays based on Gαβγ dissociation (Gαβγ) (Galés et al., 2005; Galés et al., 2006; Olsen et al., 2020), that reveals a significantly larger assay windows for EMTA for the 6 Gα subunits tested for eight selected receptors, (Figure 2—figure supplement 5).
For the Gαs translocation biosensor, the bile acid receptor (GPBA) was chosen for validation (Kawamata et al., 2003). As expected, lithocholic acid stimulation resulted in a concentration-dependent decrease in ebBRET between Gαs-RlucII and rGFP-CAAX (Figure 2D, left). Cholera toxin (CTX), which directly activates Gαs (De Haan and Hirst, 2004), led to a decrease in ebBRET (Figure 2D, center), confirming that loss of Gαs plasma membrane localization results from its activation. The potency of lithocholic acid to promote Gs dissociation from the plasma membrane was well in line with its potency to increase cAMP production as assessed using a BRET²-based EPAC biosensor (Leduc et al., 2009; Figure 2D, right). The Gs-plasma membrane dissociation ebBRET signal was not affected by UBO-QIC or PTX (Figure 2D, left), confirming the selectivity of the biosensor.
Signaling signatures of one hundred therapeutically relevant receptors reveals distinct G protein and βarrestin selectivity profiles
We used EMTA to assess the signaling signature of a panel of 100 human GPCRs that are either already the target of clinically used drugs (74 receptors), considered for pre- or clinical drug development (6 receptors), or pathophysiologically relevant (Supplementary file 2A). To establish the coupling potentials for each receptor, we quantified its ability to activate 15 pathways: Gαs, Gαi1, Gαi2, GαoA, GαoB, Gαz, Gα12, Gα13, Gαq, Gα11, Gα14, Gα15 and βarrestin 2 as well as βarrestin 1 and 2 in the presence of GRK2 (Supplementary file 3). Emax and pEC50 values were determined (Supplementary file 2) and, based on the pre-determined threshold criteria (Emax ≥mean of vehicle-stimulated +2*SD; see Materials and methods), a ‘yes or no’ agonist-dependent activation was assigned to each signaling pathway and summarized using radial graph representations (Figure 3—figure supplement 1). To assess whether endogenous receptors could contribute to the observed responses, assays were also carried out in cells not transfected with the studied receptor (Figure 3—figure supplement 2). When an agonist-promoted response was observed in non-transfected parental HEK293 cells, this response was not considered as a receptor-specific response (see Materials and methods).
To compare the signaling profiles across all receptors and pathways and to overcome differences in receptor expression levels and individual biosensor dynamic windows, we first min-max normalized Emax and pEC50 values (between 0 and 1) across receptors as a function of a reference receptor yielding the largest response for a given pathway (Figure 3A, left). Then, these values were again min-max normalized (between 0 and 1) for the same receptor across pathways, using the pathway with the largest response for this receptor as the reference (Figure 3A, right; see description in Materials and methods). Such double normalization allows direct comparison of the coupling efficiency to different G proteins for a given receptor and across receptors for a given G protein. This coupling efficiency is summarized as heatmaps (Figure 3B) that reveals a high diversity of signaling profiles. The selectivity toward the different G protein families varies considerably among GPCRs (Figure 4). In our dataset, which is the first using unmodified GPCRs and Gα proteins (except for Gs), 29% of the receptors coupled to only one family, whereas others displayed more promiscuity by coupling to 2, 3, or 4 families (36%, 25%, and 10%, respectively). Receptors coupling to a single G protein family favored the members of the Gi/o family. Indeed, 27% of the receptors coupling to Gi/o only activated this subtype family in comparison to 0, 2.4 and 9.1% for receptors activating G12/13, Gq/11, and Gs, respectively, thus displaying more promiscuous coupling. A detailed comparative analysis of the selectivity profiles that we observed using the EMTA sensors with that of the chimeric G protein-based assay developed by Inoue et al., 2019 and the IUPHAR/BPS Guide to Pharmacology database (GtP; https://www.guidetopharmacology.org/) is presented in the accompanying paper (Hauser et al., 2022). Supplementary file 2C allows a direct comparison of the relative potency determined using EMTA for both the new and the already known (i.e.: identified in GtP database) couplings. As can be seen in the table, although in many cases the potency for the novel couplings is lower, this is not a universal finding since for some receptors, the pEC50s for the new couplings are similar (ex: G12 for CB1; G13 for serotonin 5-HT2C; G12/13 for adenosine 2A (A2A) and prostaglandin E1 (EP1) receptors; Gi/o for corticotropin-releasing hormone receptor 1 (CRFR1), ETA and G protein-coupled receptor 39 (GPR39)) or higher (ex: Gz for serotonin 5-HT2B; G15 for adenosine 3 (A3) and melanocortin 3 (MC3R) receptors; G12 for bradykinin 2 (B2), cholecystokinin A (CCK1), chemokine receptor 6 (CCR6) and ETA receptors; G12/13 for CRFR1 and GPR68) than those for the canonical ones. Interestingly, in many instances the potency for the newly uncovered couplings are similar to those for βarrestins, which is generally lower than for their canonical G proteins, a finding consistent with the role of βarrestins in signaling arrest at the plasma membrane. The potency differences observed for the activation of different G protein subtypes by a given receptor may lead to preferential activation of some pathways over others. This relative selectivity is likely to be influenced by tissue-dependent G protein subtype expression levels. The physiological consequences of such selectivity remain to be investigated.
When examining the frequency of coupling for each Gα subunit family (Figure 4C), the Gi/o family members were the most commonly activated, with 89% of the tested receptors activating a Gi/o family member. In contrast, only 33%, 49%, and 45% of the receptors activate Gs, G12/13, or Gq/11 (excluding Gα15) family members, respectively. Not surprisingly, and consistent with its reported promiscuous coupling, Gα15 was found to be activated by 81% of the receptors. For some receptors, we also observed preferential coupling of distinct members within a subtype family (Figure 3—figure supplement 1). For instance, 33% of Gi/o-coupled receptors can couple to only a subpopulation of the family (Figure 4—figure supplement 1A). For the Gq/11 family, only 44% activate all family members with 45% activating only Gα15 and 11% engaging only two or three members of the family. A matrix expressing the % of receptors engaging a specific Gα subtype that also activated another subtype, is illustrated in Figure 4—figure supplement 1B. When considering individual families, considerable variation within the Gi/o family was observed. The greatest similarities were observed between GαoB and either GαoA or Gαz, and the lowest between Gαi1 and Gαz. A striking example of intra-family coupling selectivity is the serotonin 5-HT2B that activates only GαoB and Gαz and GPR65 that selectively activates GαoB. Similarly, when considering the ligand-promoted responses above our threshold criteria (see Materials and methods), histamine H2 and MC3R receptors show preferred activation of GαoB and Gαz, whereas the prostaglandin F (FP) and neuropeptide Y5 (Y5) receptors preferentially activate GαoB, GαoA, and Gαz. Even when all members of a given family are found to be activated, some receptors activate specific family members with greater potencies (Supplementary file 2C).
When considering βarrestin recruitment, our analysis shows that 22% of receptors did not recruit βarrestin 1 or 2, even in the presence of overexpressed GRK2 (Figure 4D). Among the receptors able to recruit βarrestins, only a very small number selectively recruited βarrestin1 (1.3%) or βarrestin2 (6.4%), most of them recruiting both βarrestins in the presence of GRK2 (92.3%) (Figure 4D). Overexpression of GRK2 potentiated the recruitment of βarrestin2 for 68% of receptors highlighting the importance of GRK2 expression level in determining βarrestin activation (Supplementary files 3 and 2).
Comparison with previous datasets reveals commonalities and crucial differences
We compared the signaling profiles obtained here with those presented by Inoue et al., 2019 and the GtP dataset. Of note, this comparison only considers the final reported couplings that in the Inoue’s study were based on the criteria of positive coupling if LogRAi ≥ –1 and negative coupling if LogRAi ≤ –1, and is influenced by the different cut-offs and normalization used in the two studies. A comparison of couplings using common Emax standard deviation cut-off, quantitative normalization and aggregation of G proteins into families is provided in the accompanying paper (Hauser et al., 2022). As can be seen in Supplementary file 4A, among the 70 receptors common to both studies, less couplings were detected in our study than reported in Inoue et al. for Gαs (21 vs. 28), Gαi1 (54 vs. 56), Gαq (31 vs. 34), and Gα14 (36 vs. 40). In contrast, more receptors activating Gα12 (29 vs. 23), Gαo (59 vs. 41), Gα13 (30 vs. 15), Gαz (52 vs. 37), and Gα15 (62 vs. 15) were detected in our study. When comparing with data collected in GtP, that reports couplings grouped for G protein families (i.e.: Gs, Gi/o, Gq/11, or G12/13) and not at the single G protein subtype level, we detected less couplings than what was reported in GtP for Gαs (32 vs. 37), but more for Gαi/o (89 vs. 69), Gαq/11 (81 vs. 48), and Gα12/13 (47 vs. 10), among the 99 receptors common to both datasets (Supplementary file 4B).
Altogether, the comparative analysis reveals 64% and 69% identity of couplings between the EMTA and Inoue’s or GtP datasets, respectively. Each dataset reporting unique couplings and missing couplings found in the other two datasets. The reasons for these differences are plausibly due to intrinsic differences in the assays used. For instance, for G12/13 and G15 specifically, the difference with the GtP dataset most likely results from the fact that in most cases G12/13 or G15 activation were determined indirectly since, until their recent description (G12/13: Quoyer et al., 2013; Schrage et al., 2015; G15:Inoue et al., 2019; Olsen et al., 2020), no robust readily available assay existed to monitor the activation of these G proteins.
Validation of newly identified G12/13 and G15 couplings
Given the overrepresentation of both G12/13 and G15 couplings, obtained with the EMTA assays vs. those reported by Inoue et al. and the GtP datasets, the validity of the EMTA assay to detect real productive couplings, was confirmed using orthogonal assays for selective examples not reported in the two other datasets. For G12/13, we used the PKN-based BRET biosensor detecting Rho activation downstream of either G12/13 or Gq/11 (Namkung et al., 2018) and the MyrPB-Ezrin-based BRET biosensor detecting the activation of Ezrin downstream of G12/13 (Leguay et al., 2021), both in the absence of heterologously expressed G proteins. Ligand stimulation of FP and CysLT2 receptors led to Rho and ezrin activation (Figure 3—figure supplement 3A), that were insensitive to the Gq/11 inhibitor YM-254890, confirming that these receptors activate Gα12/13.
For newly identified G15 couplings, we took advantage of the lack of Gα15 in HEK293 cells and assessed the impact of Gα15 heterologous expression on receptor-mediated calcium responses (Figure 3—figure supplement 3B). For prostaglandin E2 (EP2) and κ-opioid (κOR) receptors, which couple to G15 but no other Gq/11 members, expression of Gα15 significantly increased the PGE2- and Dynorphin A- promoted calcium responses. For α2A adrenergic (α2AAR) and vasopressin 2 (V2) receptors that couple other Gq/11 family members, treatment with YM-254890 completely abolished the agonist-promoted calcium response in the absence of Gα15. In contrast, the calcium response evoked by α2AAR and V2 agonists following Gα15 expression was completely insensitive to YM-254890 (Figure 3—figure supplement 3B), confirming that these receptors can activate this YM-254890-insensitive G protein subtype (Takasaki et al., 2004).
EMTA platform detects constitutive receptor activity and biased signaling
We went on to assess the ability of the EMTA platform to detect receptor constitutive activity. Transfection of increasing amounts of adenosine A1 receptor (A1) led to a receptor-dependent increase in basal ebBRET of the Gαi2-activation sensor (Figure 5A, left), reflecting A1 constitutive activity. The A1 inverse agonist DPCPX (Lu et al., 2014) dose-dependently decreased the constitutive A1-mediated activation of Gαi2 (Figure 5A, left), indicating that EMTA can detect inverse agonism. Although we can not exclude that the high basal activity resulted from activation by adenosine in the cell culture medium, the fact that high basal activity was observed for A1 but not A3, despite a similar potency of adenosine to activate these two receptors subtypes (see Figure 5—figure supplement 1A), supports the notion that the increased basal activity reflects A1 constitutive activity.
To further confirm that the platform can adequately detect inverse agonism, a second receptor for which no endogenous ligand should be present in the media, the CB1 receptor, was used. As illustrated in Figure 5A (right), increase CB1 expression led to a ligand-independent constitutive activation of Gz, that could be completely blocked by the CB1 inverse agonist rimonabant.
EMTA also faithfully detected biased signaling. Indeed, as previously reported (Namkung et al., 2018; Wei et al., 2003), angiotensin analogs such as SII, saralasin or TRV027 displayed biased signaling by promoting efficient βarrestin2 recruitment but marginal or no Gαq, Gαi2, or Gα13 activation as compared to angiotensin II that activated all G proteins and βarrestin2 (Figure 5B). The platform was also used to identify biased signaling resulting from single nucleotide polymorphisms. As shown in Figure 5C, two naturally occurring variants of human GPR17 (isoform 2) localised in the TM3 E/DRY motif resulted in altered functional selectivity profiles. Whereas the Asp128Asn variant displayed WT-like activity on Gαi2, it lost the ability to activate Gαq and βarrestin2. In contrast, variant Arg129His at the neighboring position resulted in an increased constitutive βarrestin2 recruitment and a loss of Gαi2 and Gαq protein signaling.
Combining Gz and G15 biosensors for safety panels and systems pharmacology
The G protein coupling profiles obtained for the 100 GPCRs revealed that 95% of receptors activate either Gαz (73%) or Gα15 (81%). Measuring activation of both pathways simultaneously provides an almost universal sensor applicable to screening. Combining the two sensors (Rap1GAP-RlucII/p63-RhoGEF-RlucII/rGFP-CAAX) in the same cells allowed to detect ligand concentration-dependent activation of a safety panel of 24 GPCRs, that are well established as contributors to clinical adverse drug reactions (Bowes et al., 2012; Figure 6—figure supplement 1). Indeed, the Gz/G15 sensor captured the activation of receptors largely or uniquely coupled to either Gαz (e.g. CB2) or Gα15 (e.g. A2A and A2B), as well as receptors coupled (to varying degrees) to both pathways. The usefulness of the Gz/G15 combined sensor to detect off-target ligand activity is illustrated in Figure 6A. Most ligands tested were specific for their primary target(s). However, certain ligands displayed functional cross-reactivity with GPCRs other than their cognate targets. These included the activation of the α2AAR by dopamine and serotonin, the D2 by noradrenaline and serotonin, and of the CB1 and CB2 receptors by acetylcholine (Figure 6B–C). The activation of D2 by noradrenaline and serotonin was confirmed by the ability of the D2-family selective antagonist eticlopride to block the dopamine-, serotonin-, and noradrenaline-promoted responses detected using the combined Gz/G15 or the Gi2- and GoB-selective sensors and βarrestin2 sensor (Figure 6B, top). Similarly, use of the α2AR selective antagonist, WB4101, allowed to confirm that dopamine can activate Gαi2, GαoB and βarrestin2 through the α2AAR (Figure 6B, bottom). Such pleiotropic activation of different monoaminergic receptors by catecholamines and serotonin has been previously observed (Roth et al., 2004; Sánchez-Soto et al., 2016; Sunahara et al., 1991). Direct activation of the α2AAR by dopamine was confirmed by showing that treatment with the D2-family receptor selective antagonist eticlopride had negligible effect on dopamine-mediated activation of Gαi2 and GαoB in cells heterologously expressing α2AAR, confirming that the response did not result from the activation of endogenously expressed dopamine receptor. In contrast, eticlopride blocked the activation of Gαi2 and GαoB in cells heterologously expressing D2 (Figure 6—figure supplement 2).
These cross-reactivity may be direct (i.e. via direct binding of a ligand to its non-cognate receptor) as suggested above, or indirect (e.g. ‘trans’, via ligand activation of its canonical receptor, leading to subsequent secretion of factors that activate the non-canonical target). One such example of trans-activation is provided by the activation of cannabinoid CB1 and CB2 receptors by acetylcholine (detected by the Gz/15 and confirmed with the GoB sensors; Figure 6A and C). Indeed, the activation was completely inhibited by both the CB inverse agonist AM-630 and by the cholinergic antagonist atropine (Figure 6C, left). Yet the response evoked by the CB selective agonist WIN55,212 2 was not blocked by atropine (Figure 6C, center). GαoB activation by acetylcholine did not result from direct activation of endogenous muscarinic receptors since no GαoB response was observed in parental cells (Figure 3—figure supplement 2). Given that the M3 muscarinic receptor, which is endogenously expressed at relatively high levels in HEK293 cells (Atwood et al., 2011), is strongly coupled to the Gq/11, CB1-expressing cells were pretreated with Gq/11/14 inhibitor UBO-QIC prior to stimulation with acetylcholine. UBO-QIC pre-treatment blocked acetylcholine- but not WIN55,212–2-mediated GαoB activation (Figure 6C, right). These results demonstrate that CB1 activation by acetylcholine is indirect and potentially involves the secretion of an endogenous CBR ligand following activation of Gq/11 by endogenous muscarinic acetylcholine receptors. The combined Gz/G15 sensor is therefore a useful tool to identify interplay between receptors and to explore systems pharmacology resulting from such cross-talks.
Discussion
This study describes the development and validation of a genetically encoded ebBRET-based biosensor platform allowing live-cell mapping of GPCR-G protein coupling preferences covering 12 heterotrimeric G proteins. The novel EMTA biosensors were combined with previously described ebBRET-based βarrestin trafficking sensors (Namkung et al., 2016), providing an unprecedented description of GPCR signaling partner couplings. In addition to providing a resource to study GPCR functional selectivity (Pándy-Szekeres et al., 2022) , the sensors provide versatile and readily usable tools to study, on a large-scale, pharmacological processes such as constitutive activity, inverse agonism, ligand-biased signaling, and signaling cross-talk.
Our EMTA-based biosensor platform offers several advantages relative to other available approaches. First, EMTA provides direct real-time measurement of proximal signaling events following GPCR activation (i.e. Gα protein activation and βarrestin recruitment) and resulting in lower level of amplification than those of assays relying on enzymatic activity of downstream effectors (i.e.: adenylyl cyclase or phospholipase C) or artificial detection systems (i.e.: gene-reporter or TGF-α shedding assays) that measure signal accumulation sometimes following extended incubation times. In addition, measuring proximal activity reduces the risk of cross-talks between pathways that may complicate data interpretation when considering downstream signaling as the readout (Mancini et al., 2015).
Second, EMTA uses native untagged GPCRs and G protein subunits (except for Gs), contrary to protein complementation (Laschet et al., 2019), FRET/BRET-based Gαβγ dissociation/receptor-G protein interaction (Bünemann et al., 2003; Galés et al., 2005; Galés et al., 2006; Hoffmann et al., 2005; Namkung et al., 2018; Olsen et al., 2020) or TGF-α shedding (Inoue et al., 2019) assays. Modifying these core-signaling components could alter responses, complicate interpretation and explain some of the discrepancies observed between the EMTA platform and other approaches used to study G protein activation. Moreover, the ability to work with unmodified receptors and G proteins (except for Gs) offers numerous advantages. First, it allows for the detection of endogenous GPCR signaling in either generic HEK293 cells (Figure 3—figure supplement 2) or more physiologically relevant cell lines such as induced pluripotent stem cell (iPSC)-derived cardiomyocytes (Figure 7A) and promyelocytic HL-60 cells (Figure 7B). Further it allows, in cells expressing sufficient endogenous level of the G proteins of interest, to detect activation of both native receptor and G proteins with no need of overexpression (Figure 7C–D). This is illustrated by the ability to detect the recruitment of Rap1GAP upon activation of the endogenous Gi/o family members by the formyl peptide receptor 2 (FPR2) in HL-60 cells (Figure 7C) or protease-activated receptor-2 (PAR2) in HEK293 cells (Figure 7D). The ability to detect the activation of endogenous G protein was also illustrated in Figure 2—figure supplement 1I, where the responses elicited by agonist stimulation were lost in cells genetically deleted of the G protein engaged by the studied receptor (i.e.: Gq/11 or Gi/o families). Recently, another BRET-based approach (Maziarz et al., 2020), taking advantage of a synthetic peptide recognizing the GTP-bound form of Gα subunits, also allows the detection of native G protein activation, offering alternative means to probe coupling selectivity profiles for both endogenously and heterologously expressed GPCRs.
Finally, similarly to BERKY, the EMTA assay platform detects the active form of the Gα subunits rather than the surrogate measurement of Gαβγ dissociation (Galés et al., 2005; Masuho et al., 2015; Maziarz et al., 2020; Mende et al., 2018), which can also detect non-productive binding as recently described for the V2 engagement of G12 (Okashah et al., 2020).
A potential caveat of EMTA is the use of common downstream effectors for all members of a given G protein family. Indeed, one cannot exclude that distinct members of a given family may display different relative affinities for their common effector. However, such differences are compensated by our data normalization that establishes the maximal response observed for a given subtype as the reference for this pathway (Figure 3A), as long as the number of the diversity of receptors included in the analysis is sufficient.
A second potential caveat of EMTA is that, when using heterologously expressed GPCRs and G proteins, some of the responses could result from favorable stoichiometries that may not exist under physiological conditions. It follows that such profiling represents the coupling possibilities of a given GPCR and not necessarily the coupling that will be observed in all cell types. Any couplings observed in such high-throughput studies requires further validation to conclude on their physiological relevance in cells or tissues of interest, and to form hypothesis for futures studies. Because we elected to use unmodified receptors (i.e.: not bearing any tags), the expression level of receptors could not be directly monitored. However, the double normalization method developed (see Materials and methods) allows quantitative comparison of coupling preferences across different receptors curtailing the influence of the assay response windows as well as receptor expression levels. Indeed, the double normalization allows ranking the coupling propensity of the receptors first as a function of the receptor which shows the strongest coupling to a specific G protein subtype, and then establishing the maximal response observed for a given G protein subtype as the reference for all G protein activated by a given receptor. In addition, as illustrated using the ETA receptor as example, titrating receptor levels did not influence the pEC50 for the activation of the different G protein coupled to this receptor (Figure 2—figure supplement 3B and Supplementary file 1C). Similarly, the pEC50 was not affected when titrating the amount of G protein subtype expressed (Figure 2—figure supplement 3A and Supplementary file 1B). As expected, only the amplitude of the response was affected.
It could be argued that overexpressing the G protein effectors (i.e.: p63-RhoGEF, Rap1GAP or PDZ-RhoGEF) used as sensors could influence the couplings observed. This potential caveat is mitigated by the fact that we used truncated part and/or modified versions of these effectors that limit the possibilities of interference with other components of the signaling machinery, and served essentially as a binding detector of the active forms of the G proteins (see Materials and methods). Supporting this notion, titrating the amount of the biosensor effector component did not affect the pEC50 of G protein activation (Figure 2—figure supplement 3C and Supplementary file 1D).
Another limitation of the EMTA platform is the lack of a soluble effector protein selective for activated Gαs thus requiring tagging of the Gαs subunit (Figure 1B, bottom) and monitoring its dissociation from the plasma membrane. Yet, our data show that this translocation reflects Gs activation state, justifying its use in a G protein activation detection platform.
Finally, because EMTA is able to detect constitutive activity, high receptor expression levels may lead to an elevated basal signal level that may obscure an agonist-promoted response. Such an example can be appreciated for the A1 receptor for which the agonist-promoted Gαi2 response did not reach the activation threshold criteria because of a very high constitutive activity level (Figure 5A). The impact of receptor expression on the constitutive activity and the narrowing on the agonist-promoted response is illustrated for Gαq activation by the 5-HT2C (Figure 5—figure supplement 1B).
A limitation of any large-scale signaling study and drug discovery program is that ligands may elicit responses downstream of receptors other than the one under study. The development of a Gz/G15 quasi-universal biosensor enables efficient screening and detection of such polypharmacology and cross-talk. Using a combination of EMTA and appropriate pharmacological tools, we also proposed a systematic approach to distinguish off-target action of ligands from cross-talk. Interestingly, the cross-talk between the M3 and CB receptors detected (Figure 6) may have physiological relevance since activation of muscarinic acetylcholine receptors has been shown to enhance the release of endocannabinoids in the hippocampus (Kim et al., 2002). The combined Gz/G15 biosensor should be particularly useful for early profiling of compound activity on safety panels and for the design of drugs displaying polypharmacology, an approach that is increasingly considered for the development of neuropsychiatric drugs (Roth et al., 2004).
The EMTA platform undoubtedly represents a novel tool-set that could be amenable for high throughput screening of small molecules and biologics across an array of signaling pathways, allowing for the discovery of functionally selective molecules or for GPCR deorphanization campaigns. The ability of the EMTA platform to quantitatively assess G protein coupling selectivity firmly expands the concept of functional selectivity and potential ligand bias beyond the dichotomic G protein vs. βarrestin view and provides plausible functional selectivity profiles that could be tested for their biological and pharmacological outcomes.
Materials and methods
Cells
HEK293 clonal cell line (HEK293SL cells), hereafter referred as HEK293 cells, were a gift from S. Laporte (McGill University, Montreal, Quebec, Canada) and previously described (Namkung et al., 2016). HEK293 cells devoid of functional Gαs (ΔGs), Gα12 and Gα13 (ΔG12/13), Gαq, Gα11, Gα14 and Gα15 (ΔGq/11) and, Gαi, and Gαo (ΔGi/o) proteins were a gift from Dr. A. Inoue (Tohoku University, Sendai, Miyagi, Japan) and previously described (Devost et al., 2017; Namkung et al., 2018; Schrage et al., 2015; Stallaert et al., 2017). Cells were maintained in Dulbecco’s Modified Eagle Medium (DMEM, Wisent, Saint-Jean-Baptiste, QC, Canada) supplemented with 10% fetal bovine serum (FBS, Wisent) and 1% antibiotics (100 U/mL penicillin and 100 μg/mL streptomycin (PS); Wisent). HL-60 cells were obtained from ATCC and maintained in RPMI 1640 medium containing L-Glutamine and 25 mM HEPES (Gibco) supplemented with 20% FBS (Wisent) and 1/100 volume PS (Wisent). Differentiation of HL-60 cells into neutrophil-like cells was induced by maintaining the cells in growth medium containing 1.3% DMSO (Bioshop) during 5 days. Cardiomyocytes derived from induced pluripotent stem cells (iPSCs; iCell Cardiomyocytes) were obtained from FUJIFILM Cellular Dynamics (Madison, WI, USA) and maintained in maintenance medium provided with the cells (special formulation by FujiFilm). Cells were grown at 37 °C in 5% CO2 and 90% humidity and checked for mycoplasma contamination.
Plasmids and ebBRET biosensor constructs
Request a detailed protocolOnly human GPCRs and human Gα subunits were used in this study. An open reading frame of each full-length GPCR was cloned into pcDNA3.1(+) expression plasmid. Except when otherwise specified, GPCRs sequences were devoid of epitope tags.
Gαs-67-RlucII (Carr et al., 2014), Gαi1-loop-RlucII and GFP10-Gγ1 (Armando et al., 2014), Gαi2-loop-RlucII and βarrestin2-RlucII (Quoyer et al., 2013), GαoB-99-RlucII (Mende et al., 2018), Gαq-118-RlucII (Breton et al., 2010), Gα12-136-RlucII and PKN-RBD-RlucII (Namkung et al., 2018), Gα13-130-RlucII (Avet et al., 2020), GFP10-Gγ2 (Galés et al., 2006), βarrestin1-RlucII (Zimmerman et al., 2012), rGFP-CAAX (Namkung et al., 2016), EPAC (Leduc et al., 2009), MyrPB-Ezrin-RlucII (Leguay et al., 2021), HA-β2AR (Lavoie et al., 2002), signal peptide-Flag-AT1 (Goupil et al., 2015), and EAAC-1 (Brabet et al., 1998) were previously described. Full-length, untagged Gα subunits, Gβ1 and Gγ9 were purchased from cDNA Resource Center. GRK2 was generously provided by Dr. Antonio De Blasi (Istituto Neurologico Mediterraneo Neuromed, Pozzilli, Italy).
To selectively detect Gi/o activation, a construct coding for aa 1–442 of Rap1 GTPase-activating protein (comprising a Gi/o binding domain) fused to Rluc8, was sequence-optimized, synthetized and subcloned at TopGenetech (St-Laurent, QC, Canada). From this construct, a RlucII-tagged version of Rap1GAP (1-442) with a linker sequence (GSAGTGGRAIDIKLPAT) between Rap1GAP and RlucII was created by Gibson assembly in pCDNA3.1_Hygro (+) GFP10-RlucII, replacing GFP10. Three substitutions (i.e. S437A/S439A/S441A) were introduced into the Rap1GAP sequence by PCR-mediated mutagenesis. These putative (S437 and S439) and documented (S441) (McAvoy et al., 2009) protein kinase A phosphorylation sites were removed in order to eliminate any Gs-mediated Rap1GAP recruitment to the plasma-membrane.
To selectively detect Gq/11 activation, a construct encoding the Gq binding domain of the human p63 Rho guanine nucleotide exchange factor (p63RhoGEF; residues: 295–502) tagged with RlucII was done from IMAGE clones (OpenBiosystems; Burlington, ON, Canada) and subcloned by Gibson assembly in pCDNA3.1_Hygro (+) GFP10-RlucII, replacing GFP10. The Gq binding domain of p63RhoGEF and RlucII were separated by the peptidic linker ASGSAGTGGRAIDIKLPAT. N-term part containing palmitoylation sites maintaining p63 to plasma membrane and part of its DH domain involved in RhoA binding/activation (Aittaleb et al., 2010; Aittaleb et al., 2011) are absent of the sensor.
To selectively detect G12/13 activation, a construct encoding the G12/13 binding domain of the human PDZ-RhoGEF (residues: 281–483) tagged with RlucII was done by PCR amplification from IMAGE clones (OpenBiosystems) and subcloned by Gibson assembly in pCDNA3.1_Hygro (+) GFP10-RlucII, replacing GFP10. The peptidic linker GIRLREALKLPAT is present between RlucII and the G12/13 binding domain of PDZ-RhoGEF. The sensor is lacking the PDZ domain of PDZ-RhoGEF involved in protein-protein interaction, as well as actin-binding domain and DH/PH domains involved in GEF activity and RhoA activation (Aittaleb et al., 2010).
The sequence of each EMTA biosensors is provided in the Supplementary file 5.
Transfection
Request a detailed protocolFor BRET experiments, HEK293 cells (1.2 mL at 3.5 × 105 cells per mL) were transfected with a fixed final amount of pre-mixed biosensor-encoding DNA (0.57 μg, adjusted with salmon sperm DNA; Invitrogen) and human receptor DNA. Transfections were performed using a polyethylenimine solution (PEI, 1 mg/mL; Polysciences, Warrington, PA, USA) diluted in NaCl (150 mM, pH 7.0; 3:1 PEI/DNA ratio). Gelatin solution (1%; Sigma-Aldrich, Saint-Louis, Missouri) was used to stabilize DNA/PEI transfection mixes. Following addition of cells to the stabilized DNA/PEI transfection mix, cells were immediately seeded (3.5 × 104 cells/well) into 96-well white microplates (Greiner Bio-one; Monroe, NC, USA) and maintained in culture for the next 48 hr in DMEM containing 2% FBS and 1% PS. DMEM medium without L-glutamine (Wisent) was used for transfection of cells with mGluR to avoid receptor activation and desensitization. For Neutrophil-like differentiated HL-60 cells, cells were resuspended in electroporation medium (growth medium containing an extra 15 mM of HEPES pH 7.0) at 25 × 106 cells/mL. Electroporation reactions were prepared by adding 50 µL of DNA mastermix (20 µg total of DNA adjusted with salmon sperm DNA, supplemented with 210 mM NaCl) to 200 µL of cell suspension and transferring into 0.4 cm gap electroporation cuvettes (Bio-Rad). The cells were electroporated at 350 µF/400 V using a Bio-Rad Gene Pulser II electroporation system, washed in electroporation medium, and seeded in 96-well plates at 0.8 × 106 cells/well in 200 µL of growth medium. BRET assays were performed 6 hr post-electroporation. For iPSC Cardiomyocytes, cells were seeded in 96-well plates pretreated with fibronectin (10 µg/ml 60 min; Sigma-Aldrich) at 3.5 × 104 cells /well. After 48 hr, attached iPSCs cells were transfected with the indicated biosensor components, using TransIT-LT1 reagent (Mirus; Madison, WI, USA), according to manufacturer recommendation. BRET assays were performed 48 hr after transfection.
For Ca2+ experiments, cells (3.5 × 104 cells/well) were co-transfected with the indicated receptor, with or without Gα15 protein, using PEI and seeded in poly-ornithine-coated 96-well clear-bottom black microplates (Greiner Bio-one) and maintained in culture for the next 48 hr.
For BRET-based imagery, cells (4 × 105 cells/dish) were seeded into 35 mm poly-d-lysine-coated glass-bottom culture dishes (Mattek Corporation; Ashland, MA, USA) in 2 ml of fresh medium and incubated at 37 °C in 5% CO2, 3 day before imaging experiments. Twenty-four hours later, cells were transfected with EMTA ebBRET biosensors and the indicated receptor (i.e. p63-RhoGEF-RlucII/rGFP-CAAX + Gαq and AT1, Rap1GAP-RlucII/rGFP-CAAX + Gαi2 and D2 or PDZ-RhoGEF-RlucII/rGFP-CAAX + Gα13 and TPαR) using X-tremeGENE 9 DNA transfection reagent (3:1 reagent/DNA ratio; Roche) diluted in OptiMEM (Gibco) and maintained in culture for the next 48 hr in DMEM containing 10% FBS and 1% PS.
Bioluminescence resonance energy transfer measurement
Request a detailed protocolEnhanced bystander BRET (ebBRET) was used to monitor the activation of each Gα protein, as well as βarrestin 1 and 2 recruitment to the plasma membrane. Gαs protein activation was measured between the plasma membrane marker rGFP-CAAX and human Gαs-RlucII in the presence of human Gβ1, Gγ9 and the tested receptor. Gαs downstream cAMP production was determined using the EPAC biosensor and GPBA receptor. Gαi/o protein family activation was followed using the selective-Gi/o effector Rap1GAP-RlucII and rGFP-CAAX along with the human Gαi1, Gαi2, GαoA, GαoB, or Gαz subunits and the tested receptor. Gαq/11 protein family activation was determined using the selective-Gq/11 effector p63-RhoGEF-RlucII and rGFP-CAAX along with the human Gαq, Gα11, Gα14, or Gα15/16 subunits and the tested receptor. Gα12/13 protein family activation was monitored using the selective-G12/13 effector PDZ-RhoGEF-RlucII and rGFP-CAAX in the presence of either Gα12 or Gα13 and the tested receptor. The expression level of the Gα subunits was monitored by western blot in HEK293 cells that endogenously expressed Gαi1, Gαi2, Gα12, Gα13, Gαq, Gα11, Gα14, and Gαs but not GαoA, GαoB, Gαz, and Gα15 (Figure 2—figure supplement 6). Gα12/13-downstream activation of the Rho pathway was measured using PKN-RBD-RlucII or Ezrin-RlucII and rGFP-CAAX with the indicated receptor. βarrestin recruitment to the plasma membrane was determined using DNA mix containing rGFP-CAAX and βarrestin1-RlucII with GRK2 or βarrestin2-RlucII alone or with GRK2 and the tested receptor. Glutamate transporters EAAC-1 and EAAT-1 were systematically co-transfected with the mGluR to prevent receptor activation and desensitization by glutamate secreted in the medium by the cells (Brabet et al., 1998). All ligands were also tested for potential activation of endogenous receptors by transfecting the biosensors without receptor DNA. The Gz/G15 biosensor consists of a combination of the following plasmids: rGFP-CAAX, Rap1GAP-RlucII, Gαz, p63-RhoGEF-RlucII and Gα15. For G protein activation detection using the BRET-based Gαβγ dissociation sensors, cells were co-transfected with untagged Gβ1 and Gαq-118-RlucII, Gα12-136-RlucII or Gα13-130-RlucII with GFP10-Gγ1, or Gαi1-loop-RlucII, Gαi2-loop-RlucII or GαoB-99-RlucII with GFP10-Gγ2, along with the indicated receptor.
The day of the BRET experiment, cells were incubated in HBSS for 1 hr at room temperature (RT). Cells were then co-treated with increasing concentrations of ligand (see Appendix 1—key resources table and Supplementary file 2 for details) and the luciferase substrate coelenterazine prolume purple (1 µM, NanoLight Technologies; Pinetop, AZ, USA) for 10 min at RT. Plates were read on a Synergy Neo microplate reader (BioTek Instruments, Inc; Winooski, VT, USA) equipped with 410 ± 80 nm donor and 515 ± 30 nm acceptor filters or with a Spark microplate reader (Tecan; Männedorf, Switzerland) using the BRET2 manufacturer settings. The BRET signal (BRET²) was determined by calculating the ratio of the light intensity emitted by the acceptor over the light intensity emitted by the donor. To validate the specificity of the biosensor responses, cells were pretreated in the absence or presence of either the Gαq inhibitor UBO-QIC (100 nM, 30 min; Institute for Pharmaceutical Biology of the University of Bonn, Germany), the Gαi/o inhibitor PTX (100 ng/mL, 18 hr; List Biological Laboratories, Campbell, California, USA) or the Gαs activator CTX (0–200 ng/mL, 4 hr; Sigma-Aldrich) before stimulation with agonist. For inverse agonist activity detection of A1 or CB1 receptors, cells were stimulated during 10 min with increasing concentrations of DPCPX or rimonabant, respectively. For ligand-cross receptor activation experiments, cells were pretreated for 10 min with increasing concentrations of antagonists or inverse agonist (eticlopride for D2, WB4101 for α2AAR, atropine for muscarinic receptors and AM-630 for CB1) before a 10 min stimulation with an EC80 concentration of the indicated agonist. BRET was measured as described above. For the safety target panel ligand screen using the combined Gz/G15 sensor, basal ebBRET level was first measured 10 min following the addition of coelenterazine prolume purple (1 µM) and ebBRET level was measured again following a 10 min stimulation with a single dose of the indicated ligand (1 μM for endothelin-1 and 10 μM for all other ligands). Technical replicates for each receptor were included on the same 96-well plate. For kinetics experiment of ETA activation, basal BRET was measured during 150 s before cells stimulation with either vehicle (DMSO) or 1 µM of endothelin-1 (at time 0 sec) and BRET signal was recorded each 30 s during 3570 s. For the validation of G12/13-mediated signal by new identified G12/13-coupled receptor using PKN- or Ezrin-based BRET biosensors, cells were pretreated or not with the Gαq inhibitor YM-254890 (1 µM, 30 min; Wako Pure Chemical Industries (Fujifilm), Osaka, Japan) before agonist stimulation for 10 min. For G protein activation detection using the BRET-based Gαβγ dissociation sensors, and for titration experiments of either Gα proteins subunit with GEMTA sensors, GPCRs with GEMTA sensors or Effector-RlucII (p63-RhoGEF-RlucII for Gαq/11, Rap1GAP-RlucII for Gαi/o or PDZ-RhoGEF-RlucII for Gα12/13) from GEMTA sensors, cells were stimulated with increasing concentrations of the indicated agonist in the presence of prolume purple for 10 min before BRET measurement. For BRET in iPSC cardiomyocytes and HL-60 cells, cells were incubated in Tyrode Hepes buffer (137 mM NaCl, 0.9 mM KCl, 1 mM MgCl2, 11.9 mM NaHCO3, 3.6 mM NaH2PO4, 25 mM HEPES, 5.5 mM D-Glucose and 1 mM CaCl2, pH 7.4) 30 min at RT before being treated with increasing concentrations of agonist for 15 min, using prolume purple (2 µM) as luciferase substrate, and BRET measured.
BRET data analyses and coupling efficiency evaluation
Request a detailed protocolAll BRET ratios were standardized using the equation below and represented as universal BRET (uBRET) values: uBRET = ((BRET ratio – A)/(B-A)) * 10,000. Constants A and B correspond to the following values:
A = pre-established BRET ratio obtained from transfection of negative control (vector coding for RlucII alone).
B = pre-established BRET ratio obtained from transfection of positive control (vector coding for a GFP10-RlucII fusion protein).
For a given signaling pathway, uBRET values at each agonist concentration were normalized as the % of the response obtained in the absence of agonist (vehicle) and concentration-response curves were fitted in GraphPad Prism 8.3 software using a four-parameter logistic nonlinear regression model. Results are expressed as mean ± SEM of at least three independent experiments.
A ligand-promoted response was considered real when the Emax value was ≥to the mean + 2*SD of the response obtained in vehicle condition and that a pEC50 value could be determined in the agonist concentration range used to stimulate the receptor. Consequently, a score of 0 or 1 was assigned to each signaling pathway depending on an agonist’s ability to activate the tested pathway (0 = no activation; 1 = activation). In the case were responses associated to endogenous receptor were detectable, we considered as ‘distorted’ and exclude all the responses observed in the presence of transfected receptor for which Emax was ≤to 2*mean of the Emax value obtained with endogenous receptors or pEC50 was ≥to 2*mean of the pEC50 value obtained with endogenous receptors. Consequently, a score of 0 was assigned for these distorted responses in radial graph representation (Figure 3—figure supplement 1) and concentration-response curves were placed on a gray background in signaling signature profile panels (Supplementary file 3). Whenever transfected receptors produced an increase in Emax or a left-shift in pEC50 values compared to endogenous receptors, responses were considered ‘true’ and were assigned with a score of 1 for radial graph representation (Figure 3—figure supplement 1) and concentration-response curves were placed on a yellow background in signaling signature profile panels to indicate a partial contribution of endogenous receptors (Supplementary file 3).
We used a double normalization of Emax and pEC50 values to compare the signaling efficiency obtained for the 100 GPCRs across all receptors and pathways. Emax and pEC50 values deduced from concentration-response curves were first normalized between 0 and 1 across receptors by ranking the receptors as a function of the receptor that most efficiently activate a given pathway and then using the activation value for the pathway (including G protein and βarrestin subtypes) that a given receptor most efficiently activate as a reference for the other pathways that can be activated by this receptor. This double normalization can be translated in the following formalized equation:
STEP1: For each receptor and for each pathway:
Pathway A = Pathway specific normalized score for GPCRx on pathway A ([PSNS GPCRx]Pathway A)
where: GPCRx is receptor being analyzed, GPCRRef is the receptor giving greatest Emax on pathway A of all receptors studied (i.e. reference receptor for pathway A). A PSNS was determined for every receptor and every pathway coupled to that receptor.
STEP2: For any given receptor:
= Normalized pathway A coupling score for GPCRx
where: [PSNS GPCRx] Pathway A is the pathway specific normalized score for GPCRx on pathway A, and [PSNS GPCRx] Ref pathway is the pathway specific normalized score for the pathway giving the highest PSNS for GPCRx (i.e., reference pathway for GPCRx).
For the safety target panel ligand screen using the combined Gz/G15 sensor, the fold ligand-induced stimulation was calculated for each receptor by dividing the BRET ratio after ligand addition (measured at 10 min post stimulation) by the basal BRET ratio prior to receptor stimulation. Activation thresholds were defined as the mean + 2*SD of the ligand-stimulated response obtained with receptor-null cells expressing only the combined Gz/G15 sensor.
Ca2+ mobilization assay
Request a detailed protocolThe day of the experiment, cells were incubated with 100 μL of a Ca2+-sensitive dye-loading buffer (FLIPR calcium five assay kit, Molecular Devices; Sunnyvale, CA, USA) containing 2.5 mM probenecid (Sigma-Aldrich) for 1 hr at 37 °C in a 5% CO2 incubator. During a data run, cells in individual wells were exposed to an EC80 concentration of agonist, and fluorescent signals were recorded every 1.5 s for 3 min using the FlexStation II microplate reader (Molecular Devices). For receptors that also activate other Gq/11 family members, cells were pretreated with the Gq/11 inhibitor YM-254890 (1 µM, 30 min) before agonist stimulation. Gα15 is resistant to inhibition by YM-254890, thus allowing to measure Ca2+ responses generated specifically by Gα15.
BRET-based imaging
Request a detailed protocolBRET images were obtained as previously described (Kobayashi et al., 2019). Briefly, the day of imaging experiment, cells were carefully rinsed with HBSS, and images were acquired before and after agonists addition (100 nM for Angiotensin II and U46619, and 1 µM for dopamine) diluted in HBSS in the presence of the luciferase substrate coelenterazine prolume purple (20 µM).
Images were recorded using an inverted microscope (Nikon Eclipse Ti-U) equipped with x60 objective lens (Nikon CFI Apochromat TIRF) and EM-CCD camera (Nuvu HNu 512). Measurements were carried out in photon counting mode with EM gain 3000. Exposure time of each photon counting was 100ms. Successive 100 frames were acquired alternatively with 480 nm longpass filter (acceptor frames) or without filter (total luminescence frames), and integrated. Image integrations were repeated 10 times and 1000 frames (Video 1) or 5 times and 500 frames (Videos 2 and 3) of acceptor and total luminescence were used to generate each image.
BRET values were obtained by dividing acceptor counts by total luminescence counts pixelwise. BRET values from 0.0 to 0.8 (Video 1) or 0.0–0.5 (Videos 2 and 3) were allocated to ‘jet’ heatmap array using MATLAB 2019b. Brightness of each pixel was mapped from the signal level of total luminescence image. 0% and 99.9% signal strength were allocated to the lowest and highest brightness to exclude the influence of defective pixels with gamma correction factor of 2.0.
The movies were generated using ImageJ 1.52 a. Frame rate is 10 (Video 1) or 3 (Videos 2 and 3) frames/s, and frame interval is 20 or 100 s for Videos 1 and 2–3, respectively. The field of view of the movie is 137 µm x 137 µm.
Western blot analysis
Request a detailed protocolCells were transfected or not with the indicated biosensors mix as previously described and whole-cell extracts were prepared 48 hr later. Briefly, cells were washed with ice-cold PBS and lysed in a buffer containing 10 mM Tris buffer (pH 7.4), 100 mM NaCl, 1 mM EDTA, 1 mM EGTA, 0.1% SDS, 1% Triton X-100, 10% Glycerol supplemented with protease inhibitors cocktails (Thermo Fisher Scientific). Cell lysates were centrifuged at 13,000 × g for 30 min at 4 °C. Equal amounts of proteins were separated by SDS-PAGE and transferred onto polyvinylidene fluoride membrane. The membranes were blocked (1 hr incubation at RT in TBS, 0.1% Tween-20, 5% BSA) and successively probed with primary antibody and appropriate goat secondary antibodies coupled to horseradish peroxidase (see Appendix 1-Key Resources Table). Western blots were visualized using enhanced chemiluminescence and detection was performed using a ChemiDoc MP Imaging System (BioRad). Relative densitometry analysis on protein bands was performed using MultiGauge software (Fujifilm). Results were normalized against control bands.
Statistical analyses
Request a detailed protocolCurve fitting and statistical analyses were performed using GraphPad Prism 8.3 software and methods are described in the legends of the figures. Significance was determined as p < 0.05.
Appendix 1
Key Resources Table
Data availability
All data generated or analysed during this study are included in the manuscript and supporting file; Source Data files have been provided for Figures 2, 4, 5, 6 and 7 and associated figure supplements. Supplementary file 1 contains the numerical data used to generate Figure 2-figure supplement 1 and Figure 2-figure supplement 3. Supplementary File 2 contains the numerical data used to generate figure 3. Further data are available in the companion paper co-published in eLife: https://doi.org/10.7554/eLife.74107.
References
-
Plasma membrane association of p63 Rho guanine nucleotide exchange factor (p63RhoGEF) is mediated by palmitoylation and is required for basal activity in cellsThe Journal of Biological Chemistry 286:34448–34456.https://doi.org/10.1074/jbc.M111.273342
-
Reducing safety-related drug attrition: the use of in vitro pharmacological profilingNature Reviews. Drug Discovery 11:909–922.https://doi.org/10.1038/nrd3845
-
Multiplexing of multicolor bioluminescence resonance energy transferBiophysical Journal 99:4037–4046.https://doi.org/10.1016/j.bpj.2010.10.025
-
Development and characterization of pepducins as Gs-biased allosteric agonistsThe Journal of Biological Chemistry 289:35668–35684.https://doi.org/10.1074/jbc.M114.618819
-
Gz, a guanine nucleotide-binding protein with unique biochemical propertiesThe Journal of Biological Chemistry 265:2383–2390.
-
Cholera toxin: a paradigm for multi-functional engagement of cellular mechanisms (Review)Molecular Membrane Biology 21:77–92.https://doi.org/10.1080/09687680410001663267
-
Conformational Profiling of the AT1 Angiotensin II Receptor Reflects Biased AgonismG Protein Coupling, and Cellular Context. J Biol Chem 292:5443–5456.https://doi.org/10.1074/jbc.M116.763854
-
The evasive nature of drug efficacy: implications for drug discoveryTrends in Pharmacological Sciences 28:423–430.https://doi.org/10.1016/j.tips.2007.06.005
-
Probing the activation-promoted structural rearrangements in preassembled receptor-G protein complexesNature Structural & Molecular Biology 13:778–786.https://doi.org/10.1038/nsmb1134
-
Angiotensin II type I and prostaglandin F2alpha receptors cooperatively modulate signaling in vascular smooth muscle cellsThe Journal of Biological Chemistry 290:3137–3148.https://doi.org/10.1074/jbc.M114.631119
-
Trends in GPCR drug discovery: new agents, targets and indicationsNature Reviews. Drug Discovery 16:829–842.https://doi.org/10.1038/nrd.2017.178
-
Modulation of rap activity by direct interaction of Galpha(o) with Rap1 GTPase-activating proteinThe Journal of Biological Chemistry 274:21507–21510.https://doi.org/10.1074/jbc.274.31.21507
-
A G protein-coupled receptor responsive to bile acidsThe Journal of Biological Chemistry 278:9435–9440.https://doi.org/10.1074/jbc.M209706200
-
Biased Receptor Signaling in Drug DiscoveryPharmacological Reviews 71:267–315.https://doi.org/10.1124/pr.118.016790
-
Activation of muscarinic acetylcholine receptors enhances the release of endogenous cannabinoids in the hippocampusThe Journal of Neuroscience 22:10182–10191.https://doi.org/10.1523/JNEUROSCI.22-23-10182.2002
-
A dynamic and screening-compatible nanoluciferase-based complementation assay enables profiling of individual GPCR-G protein interactionsThe Journal of Biological Chemistry 294:4079–4090.https://doi.org/10.1074/jbc.RA118.006231
-
Beta 1/beta 2-adrenergic receptor heterodimerization regulates beta 2-adrenergic receptor internalization and ERK signaling efficacyThe Journal of Biological Chemistry 277:35402–35410.https://doi.org/10.1074/jbc.M204163200
-
Functional selectivity of natural and synthetic prostaglandin EP4 receptor ligandsThe Journal of Pharmacology and Experimental Therapeutics 331:297–307.https://doi.org/10.1124/jpet.109.156398
-
Development of conformational BRET biosensors that monitor ezrin, radixin and moesin activation in real timeJournal of Cell Science 134:jcs255307.https://doi.org/10.1242/jcs.255307
-
Structure of Galphaq-p63RhoGEF-RhoA complex reveals a pathway for the activation of RhoA by GPCRsScience (New York, N.Y.) 318:1923–1927.https://doi.org/10.1126/science.1147554
-
Exploring the Technology Landscape of 7TMR Drug Signaling ProfilingCurrent Topics in Medicinal Chemistry 15:2528–2542.https://doi.org/10.2174/1568026615666150701113344
-
Activated G Protein Galphas Samples Multiple Endomembrane CompartmentsThe Journal of Biological Chemistry 291:20295–20302.https://doi.org/10.1074/jbc.M116.729731
-
Functional interaction between Galpha(z) and Rap1GAP suggests a novel form of cellular cross-talkThe Journal of Biological Chemistry 274:36663–36669.https://doi.org/10.1074/jbc.274.51.36663
-
Heterotrimeric G protein activation by G-protein-coupled receptorsNature Reviews. Molecular Cell Biology 9:60–71.https://doi.org/10.1038/nrm2299
-
TRUPATH, an open-source biosensor platform for interrogating the GPCR transduceromeNature Chemical Biology 16:841–849.https://doi.org/10.1038/s41589-020-0535-8
-
The G protein database, GproteinDbNucleic Acids Research 50:D518–D525.https://doi.org/10.1093/nar/gkab852
-
Galphaq directly activates p63RhoGEF and Trio via a conserved extension of the Dbl homology-associated pleckstrin homology domainThe Journal of Biological Chemistry 282:29201–29210.https://doi.org/10.1074/jbc.M703458200
-
Magic shotguns versus magic bullets: selectively non-selective drugs for mood disorders and schizophreniaNature Reviews. Drug Discovery 3:353–359.https://doi.org/10.1038/nrd1346
-
The experimental power of FR900359 to study Gq-regulated biological processesNature Communications 6:10156.https://doi.org/10.1038/ncomms10156
-
A novel Galphaq/11-selective inhibitorThe Journal of Biological Chemistry 279:47438–47445.https://doi.org/10.1074/jbc.M408846200
-
Activation-induced subcellular redistribution of Gs alphaMolecular Biology of the Cell 7:1225–1233.https://doi.org/10.1091/mbc.7.8.1225
Article and author information
Author details
Funding
Canada Research Chairs
- Michel Bouvier
Canadian Institutes of Health Research (FDN-148431)
- Michel Bouvier
Lundbeckfonden (R218-2016-1266)
- David E Gloriam
Lundbeckfonden (R313-2019-526)
- David E Gloriam
Novo Nordisk Fonden (NNF18OC0031226)
- David E Gloriam
Lundbeckfonden (R278-2018-180)
- Alexander S Hauser
Bristol-Myers Squibb
- Michel Bouvier
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Shane C Wright for scientific discussion and Monique Lagacé for critical reading of the manuscript. The authors are grateful to the funding from Bristol-Myers Squibb that supported the detection of Gα proteins by endogenous receptors in iPSC cardiomyocytes and promyelocytic HL-60 cells.
Copyright
© 2022, Avet et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 8,085
- views
-
- 1,375
- downloads
-
- 142
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
The conformational ensemble and function of intrinsically disordered proteins (IDPs) are sensitive to their solution environment. The inherent malleability of disordered proteins, combined with the exposure of their residues, accounts for this sensitivity. One context in which IDPs play important roles that are concomitant with massive changes to the intracellular environment is during desiccation (extreme drying). The ability of organisms to survive desiccation has long been linked to the accumulation of high levels of cosolutes such as trehalose or sucrose as well as the enrichment of IDPs, such as late embryogenesis abundant (LEA) proteins or cytoplasmic abundant heat-soluble (CAHS) proteins. Despite knowing that IDPs play important roles and are co-enriched alongside endogenous, species-specific cosolutes during desiccation, little is known mechanistically about how IDP-cosolute interactions influence desiccation tolerance. Here, we test the notion that the protective function of desiccation-related IDPs is enhanced through conformational changes induced by endogenous cosolutes. We find that desiccation-related IDPs derived from four different organisms spanning two LEA protein families and the CAHS protein family synergize best with endogenous cosolutes during drying to promote desiccation protection. Yet the structural parameters of protective IDPs do not correlate with synergy for either CAHS or LEA proteins. We further demonstrate that for CAHS, but not LEA proteins, synergy is related to self-assembly and the formation of a gel. Our results suggest that functional synergy between IDPs and endogenous cosolutes is a convergent desiccation protection strategy seen among different IDP families and organisms, yet the mechanisms underlying this synergy differ between IDP families.
-
- Biochemistry and Chemical Biology
- Stem Cells and Regenerative Medicine
Human induced pluripotent stem cells (hiPSCs) have great potential to be used as alternatives to embryonic stem cells (hESCs) in regenerative medicine and disease modelling. In this study, we characterise the proteomes of multiple hiPSC and hESC lines derived from independent donors and find that while they express a near-identical set of proteins, they show consistent quantitative differences in the abundance of a subset of proteins. hiPSCs have increased total protein content, while maintaining a comparable cell cycle profile to hESCs, with increased abundance of cytoplasmic and mitochondrial proteins required to sustain high growth rates, including nutrient transporters and metabolic proteins. Prominent changes detected in proteins involved in mitochondrial metabolism correlated with enhanced mitochondrial potential, shown using high-resolution respirometry. hiPSCs also produced higher levels of secreted proteins, including growth factors and proteins involved in the inhibition of the immune system. The data indicate that reprogramming of fibroblasts to hiPSCs produces important differences in cytoplasmic and mitochondrial proteins compared to hESCs, with consequences affecting growth and metabolism. This study improves our understanding of the molecular differences between hiPSCs and hESCs, with implications for potential risks and benefits for their use in future disease modelling and therapeutic applications.