Disentangling bias between Gq, GRK2, and arrestin3 recruitment to the M3 muscarinic acetylcholine receptor
Abstract
G protein-coupled receptors (GPCRs) transmit extracellular signals to the inside by activation of intracellular effector proteins. Different agonists can promote differential receptor-induced signaling responses – termed bias – potentially by eliciting different levels of recruitment of effector proteins. As activation and recruitment of effector proteins might influence each other, thorough analysis of bias is difficult. Here, we compared the efficacy of seven agonists to induce G protein, G protein-coupled receptor kinase 2 (GRK2), as well as arrestin3 binding to the muscarinic acetylcholine receptor M3 by utilizing FRET-based assays. In order to avoid interference between these interactions, we studied GRK2 binding in the presence of inhibitors of Gi and Gq proteins and analyzed arrestin3 binding to prestimulated M3 receptors to avoid differences in receptor phosphorylation influencing arrestin recruitment. We measured substantial differences in the agonist efficacies to induce M3R-arrestin3 versus M3R-GRK2 interaction. However, the rank order of the agonists for G protein- and GRK2-M3R interaction was the same, suggesting that G protein and GRK2 binding to M3R requires similar receptor conformations, whereas requirements for arrestin3 binding to M3R are distinct.
Editor's evaluation
This paper investigates the molecular mechanism of ligand bias in G protein-coupled receptors, specifically the M3 muscarinic receptor, which is the property that different receptor agonists favor activation of G protein signaling vs. arrestin-mediated signaling. The interaction with arrestin is promoted by receptor phosphorylation by G protein-coupled receptor kinases (GRK), and GRK recruitment is expected to influence arrestin-recruiting activity of a particular ligand. However, the possibility of a distinct agonist-dependent receptor conformation on GRK association has not been investigated. This study demonstrates that the G protein Gq and GRK2 appear to interact with a similar M3 receptor conformation, whereas arrestin3 interacts with a distinct conformation. This represents a significant advance in understanding the mechanism of ligand bias in G protein-coupled receptors.
https://doi.org/10.7554/eLife.58442.sa0Introduction
G protein-coupled receptors (GPCRs) are membrane-spanning proteins that convert extracellular to intracellular signals. GPCRs can sense a wide variety of different agents, from single photons in the case of rhodopsin to small proteins that activate chemokine receptors (Shichida and Imai, 1998; Griffith et al., 2014). The canonical sequence of events during receptor activation (not taking into account precoupling of G protein) is the following: (1) an agonist, increasing the population of active receptor conformations, binds to the orthosteric pocket; this leads to a further conformational change of the transmembrane helix bundle that results in an opening of the intracellular effector-binding cavity; (2) G proteins bind to the intracellular pocket of the receptors, are thereby activated and enabled to transmit and modulate a multitude of signals; (3) G protein-coupled receptor kinases (GRKs) are recruited to the activated receptors, leading to receptor phosphorylation; and (4) arrestins bind to the active conformation of these receptors after being activated by an initial binding to phosphorylated residues at the intracellular end of the receptors. Subsequently, arrestins induce internalization and desensitization of the receptors (Ritter and Hall, 2009; Martí-Solano et al., 2016; Scheerer and Sommer, 2017; Kenakin, 2019). In the sequence of events described above, the three GPCR-binding proteins investigated in this work (G proteins, GRKs, and arrestins) recognize agonist-induced receptor conformations independently of each other, as shown in Figure 1—figure supplement 1. There is strong evidence that different agonists can stabilize distinct receptor conformations and thereby influence the recruitment of effector proteins leading to ‘functional selectivity’ or, more specifically, ‘ligand bias’ to better reflect the nomenclature in the upcoming BJP guidelines (Costa-Neto et al., 2016; Martí-Solano et al., 2016; Wingler et al., 2019; Xu et al., 2019; Kenakin, 2019). It has been shown before that receptor conformations exist that are preferentially recognized by either arrestins or G proteins. What these different receptor conformations might look like is a hot topic of current research (Wingler et al., 2020; Suomivuori et al., 2020). The recent study of Stoeber et al., 2020 has shown that agonist-dependent recruitment of GRK2 to an opioid receptor exists. This leads to the question whether a functional bias at the level of arrestin recruitment (Carr et al., 2016) is affected by both receptor recruitment rates of GRKs and arrestins. Experimentally it has been difficult to distinguish between both events. Here, we set out to address agonist selectivity independently for all three receptor interaction events by comparing agonist-induced binding of Gq proteins, GRK2, as well as arrestin3 to muscarinic acetylcholine M3 receptors (M3R) as a model. These effector proteins and their recruitment were monitored by means of single-cell Förster resonance energy transfer (FRET) imaging under conditions specifically optimized to minimize interference of upstream events (G protein activation in case of GRK2 and prephosphorylation in case of arrestin). By investigating the effect of seven different agonists on these effectors, we answered the following two questions: firstly, can we detect biased recruitment of effector proteins at the M3R between these seven agonists? Secondly, is a potential bias in arrestin3 recruitment caused by a bias in GRK2 recruitment or can both arrestin3 and GRK2 recruitment be biased differently in comparison to Gq binding and activation?
Our results show that ligand-induced biased recruitment can indeed be detected at the M3R. Assessing bias as a change in sequence when ranking the agonists with respect to their efficacy to induce Gq, GRK2, or arrestin3 recruitment, we found a difference between arrestin3 recruitment and the other two effector proteins. In contrast, no such difference was found between GRK2 recruitment and G protein binding and activation, suggesting that very similar active receptor conformations are required for Gq activation and GRK2 recruitment.
Results
Assays measuring agonist bias at the M3R
We selected seven agonists with pharmacological relevance or structural similarity to a pharmacologically relevant agonist (Figure 1—figure supplement 2): acetylcholine (ACh), arecoline (Are), methacholine (Metha), pilocarpine (Pilo), guvacoline (Guva), methyl-3-dimethylaminopropionate (Mda), and 5-methyl-furmethiodide (Fur). The ability of every agonist to enable Gq, GRK2, or arrestin3 binding to the M3R was investigated with single-cell FRET-based measurements in human embryonic kidney 293T (HEK293T) cells. As expected, receptor stimulation with 30 µM ACh resulted in a reversible increase in YFP fluorescence and a corresponding decrease in mTurq fluorescence, reflecting the FRET development due to Gq, GRK2, or arrestin3 binding to the M3R, respectively (Figure 1—figure supplement 3). Thus, we were able to compare multiple agonists in the same cell under the same conditions.
Agonist efficacy is reflected in the amplitudes of M3R-Gq protein binding
For many receptors including M3R, it is well known that signal amplification at the level of G protein activation leads to a spare receptor phenomenon, which allows only indirect determination of agonist efficacy by means of comparing of concentration-response curves with those of ligand binding. Therefore, we determined agonist evoked binding of Gq to M3R by means of single-cell FRET recording as a direct measurement of agonist efficacy (Figure 1A and B). During laminar superfusion of cells, we sequentially switched between buffer and all seven agonists used in the study. Due to fast kinetics, we were able to wash out the respective agonists before we added the next. In order to get a robust assessment of agonist efficacy, we used high agonist concentrations for all agonists and kept them constant for comparison in the GRK2 and arrrestin3 recruitment assays. At the end of the protocol, we reapplied 30 µM ACh in order to test for the stability of the signal. Our results showed the highest recruitment of Gq for ACh, followed by Metha and Fur (∼60% of the ACh signal), and much weaker recruitment for Guva and Are (∼20%) followed by even lower values for Pilo and Mda (∼8%) (Figure 1B), leading to the following ranking in efficacy: ACh>Metha≈Fur>Are≈Guva>Pilo>Mda. We also performed bioluminescence resonance energy transfer (BRET)-based assays to study Gq activation by these agonists (Figure 1—figure supplement 4), with the following similar ranking in potency: ACh>Metha≈Fur>Are>Guva>Pilo>Mda. In addition, radioligand displacement assays were performed (Figure 1—figure supplement 5) in order to determine the rank order of agonist affinity: Pilo>ACh>Metha≈Fur>Are>Guva»Mda. Mda failed to displace the radioligand at concentrations up to 100 mM, indicating a very low affinity for the orthosteric binding pocket of M3R or – despite its structural similarity to ACh – a noncompetitive binding mode.

Measuring G protein binding to the M3R upon stimulation with muscarinic receptor agonists (A).
(B) G protein-M3R interaction was measured as described in Figure 1—figure supplement 3 by means of single-cell Förster resonance energy transfer (FRET) recording. The cells were stimulated with distinct concentrations of each agonist as indicated. The amplitudes were normalized to the ACh amplitude of every cell. An ordinary one-way ANOVA (****p<0.0001) with Tukey’s multiple comparison test, and additionally for paired measurements a paired Student’s t-test, was conducted. Only the significant differences (ANOVA: black, t-test: purple) in terms of biased recruitment to the M3R are shown. All data are plotted as mean values ± SD for each condition.
-
Figure 1—source data 1
Source data related to Figure 1B.
- https://cdn.elifesciences.org/articles/58442/elife-58442-fig1-data1-v2.xlsx
G protein-independent GRK2 recruitment to the M3R
Our aim to measure and compare the efficacy of the agonists also at the level of GRK2 recruitment to M3R was complicated by the fact that GRK2 is activated by G protein-induced recruitment to the plasma membrane. The rank order of the efficacy of the different agonists to recruit GRK2 to M3R, determined by single-cell FRET assays, was the same as determined for Gq recruitment (Figure 2—figure supplement 1A and B in comparison with Figure 1A and B). However, due to possible interference with Gq in the recruitment of M3R, these data are difficult to interpret. Indeed, when Gq activation was completely abolished by pretreatment with FR900359, no recruitment of GRK2 to M3R was observed in cells transfected with Gq and fluorescently labeled M3R and GRK2 (Figure 2—figure supplement 1C). Likewise, preincubation with PTX completely abolished agonist-evoked FRET between Gβγ and GRK2 (Figure 2—figure supplement 1D). Interestingly, without overexpression of Gαq proteins, preincubation with PTX and FR900359 pretreatment did not lead to a full inhibition of agonist evoked FRET between Gβγ and GRK2 (Figure 2—figure supplement 1E), presumably due to bystander FRET upon G protein-independent translocation of GRK2 to M3R. Kinetic analysis of agonist-evoked FRET between Gβγ and GRK2 with or without pretreatment with PTX and FR900359 showed fast one-phased decays for signals measured with pretreatment and slower two-phased decays for signals measured without pretreatment (Figure 2—figure supplement 1F and G). In order to abolish the possible involvement of Gβγ and Gαq proteins in the recruitment process of GRK2 to M3R (Wolters et al., 2015) further, we aimed to bypass Gq-dependent translocation of GRK2 to the plasma membrane by introducing a CAAX-box at the C-terminus of GRK2-YFP to allow for G protein-independent plasma membrane localization of GRK2 (Inglese et al., 1992). In addition, we mutated D110 to alanine, which attenuates Gαq binding to GRK2 (Wolters et al., 2015). Under conditions of complete Gq and Gi inhibition (see for control Figure 2—figure supplement 1D and E), amplitudes of agonist-evoked FRET increases were reduced and offset kinetics after agonist withdrawal were substantially faster (compare Figure 2A and Figure 2—figure supplement 1A). However, we were still able to measure agonist-specific interactions between (D110A)GRK2-CAAX and M3R (Figure 2A and B, Figure 2—figure supplement 1H). Our measurements of the relative efficacy of the agonists to recruit this GRK2 construct without the influence of G proteins resulted in the following ranking: ACh>Metha≈Fur>Guva≈Are>Pilo (Figure 2B). Mda did not lead to detectable responses. Even though the relative responses induced by the weak partial agonists were substantially reduced compared to the G protein-dependent recruitment of GRK2 (Figure 2—figure supplement 1A and B), the rank order of the agonists was essentially unchanged. We ensured that the color switch of the attached fluorophores had no effect on the relative agonist efficacies (Figure 2—figure supplement 1I). Control experiments without overexpression of Gαq, but in the context of non-mutated fluorescent GRK2, showed a similar efficacy of the partial agonists Pilo and Are for GRK2 binding to M3R compared to those measured under complete G protein inhibition with the G protein-insensitive GRK2 (Figure 2—figure supplement 1J and K). These results demonstrate, as previously reported (Wolters et al., 2015), that the overexpression of Gq proteins leads to an enhanced affinity of GKR2 for the M3R.

Bias in arrestin3 recruitment is distinct from bias in GRK2 recruitment.
(A, B) M3R-(D110A)GRK2-CAAX interaction without Gαq overexpression was measured with a single-cell Förster resonance energy transfer (FRET)-based assay, after cells were preincubated with 50 ng/mL of PTX overnight and additionally for 10 min before the start of the measurement, while continuously superfused with FR900359 at a concentration of 1 µM during the whole measurement. (C–F) M3R-arrestin3 interaction was measured with a single-cell FRET-based assay. Each cell was stimulated for 30 s with saturating concentrations of each indicated agonist. (C) Stimulation with agonist (first stimulation), thereafter with ACh and after that with agonist again (second stimulation) shows that M3R-arrestin3 interaction increases with repeated stimulation by an agonist in the same cell. (D) Mean amplitude values of first and second stimulation were statistically analyzed with a paired t-test. (B, D, F) All amplitudes were normalized to the ACh amplitude of the same cell (norm. to max. ACh). One-way ANOVA (****p<0.0001) with Tukey’s multiple comparison test, and additionally for paired measurements a paired Student’s t-test, was conducted. Only the significant differences (ANOVA: black, t-test: purple) in terms of biased recruitment to the M3R are shown. Data are shown as mean ± SD, *p<0.05, **p<0.001, ***p<0.001, ****p<0.0001. The number of experiments is indicated in parentheses.
-
Figure 2—source data 1
Source data related to Figure 2B, D and F.
- https://cdn.elifesciences.org/articles/58442/elife-58442-fig2-data1-v2.xlsx
Agonist-induced arrestin3 recruitment to prestimulated M3R
Previous studies have shown that GRK2-mediated phosphorylation of the receptor is initially the rate-limiting step of arrestin binding and that β2-adrenergic receptor-arrestin interaction and M3R-arrestin interaction accelerate with repeated stimulation of the same cell (Krasel et al., 2005; Wolters et al., 2015). Since GRK-mediated receptor phosphorylation is in most cases a prerequisite for agonist-evoked arrestin3 binding (Choi et al., 2018), it is difficult to selectively measure the efficacy of agonists to induce binding of arrestin3 to (phosphorylated) receptors. The dependency on agonist-induced phosphorylation for arrestin binding is exemplarily shown in Figure 2C and D as Metha-induced arrestin3 recruitment is significantly enhanced after a prepulse with the full agonist ACh.
To minimize the influence on arrestin3 recruitment by each agonist’s potentially different ability to activate G proteins and subsequently GRK2, we analyzed the time dependency and GRK2 dependency of the arrestin3 recruitment (Figure 2—figure supplement 2A–G) and developed a prestimulation protocol. The amplitudes measured for prestimulated and not prestimulated M3R differed for all seven agonists except for ACh (Figure 2—figure supplement 2D). This indicates that the increase in amplitude after prestimulation is specific for each agonist (Figure 2—figure supplement 2E), presumably due to differences in their capability to induce phosphorylation of M3R (Butcher et al., 2011). The prestimulation protocol contained a first stimulation with saturating concentrations of ACh, followed by a brief washout, immediately followed by test pulses (Figure 2E). We verified the reliability of the protocol by comparing relative FRET amplitudes, evoked by the test compounds when applied in different order, with the amplitudes obtained directly after prestimulation with ACh (Figure 2—figure supplement 3). By applying each of the seven agonists to the M3R subsequent to the prestimulation with ACh, we obtained results that show marked differences in the ability of the individual agonists to recruit arrestin3 to M3R (Figure 2E and F). This led to the following efficacy ranking of the agonists: ACh>Metha≈Fur>Pilo≈Guva>Are>Mda. It is important to note that the rank order for the efficacy to recruit arrestin3 to M3R differs in comparison to the one obtained for Gq recruitment and activation.
Agonist-induced M3R-arrestin3 recruitment can be different from M3R-Gq and M3R-GRK2 recruitment
In order to compare the relative efficacies (normalized to ACh) of all agonists tested for binding of the three different effector proteins to M3R as determined from experiments shown in Figures 1B, 2B and F, we used 2D plots and compared all three modalities pairwise, as depicted in Figure 3A–C. As shown in Figure 3A, the agonist efficacies to recruit Gq to M3R in comparison to D110A-GRK-CAAX were very similar, showing nearly equal efficacy for each agonist. In contrast, both 2D plots comparing agonist efficacies for arrestin3 binding to M3R with those for either Gq or D110A-GRK2-CAAX showed higher relative agonist efficacies for arrestin binding. Most importantly, the rank order of agonist efficacies was different for arrestin3 binding as Pilo was more efficient than Are and similarly efficient as Guva to induce arrestin3 recruitment, whereas for both GRK2 recruitment and Gq binding to M3R, Pilo was less efficient than Are and Guva. We further determined efficacies of the agonists with respect to G protein activation (Figure 1—figure supplement 4), GRK2 recruitment with Gq expression and without G protein inhibitors (Figure 3—figure supplement 1B, Figure 3—figure supplement 2) as well as arrestin3 recruitment (Figure 3—figure supplement 1A, Figure 3—figure supplement 2). Using the determined binding affinities for the agonists (Figure 1—figure supplement 5), the obtained curves were fitted according to the operational model of Black and Leff, 1997, yielding their respective τ values. The values of the selected agonists relative to ACh in all three assays were then plotted in a radar plot (Figure 3—figure supplement 2). The results confirmed that agonist efficacies for GRK2 recruitment and Gq activation followed a similar pattern with the identical rank order of agonists, whereas in the case of arrestin3 binding, at least Pilo exhibited a reverted rank order relative to Are and Guva.

Measuring GRK2 recruitment disentangled from Gq reveals no bias between Gq and GRK2 recruitment to M3R.
Recruitment of Gq (data from Figure 1) was plotted in relation to (A) recruitment of (D110A)GRK2-CAAX and (B) recruitment of arrestin3 (data from Figure 2) for six agonists normalized to the maximum amplitude of ACh at a concentration of 30 µM. (C) The normalized M3R-arrestin3 recruitment was plotted in relation to the M3R-(D110A)GRK2-CAAX recruitment.
Discussion
The notion that different agonists can stabilize different conformational states of the receptor is well accepted within the field of GPCR pharmacology, and numerous studies indicate that selective activation of signaling pathways might at least in part be induced by these different conformations (Wingler et al., 2019). So far, a direct comparison of the bias of agonists at the level of binding of G proteins, GRKs, and arrestins to receptors has been lacking. Instead of analyzing downstream signaling, we directly investigated the recruitment of Gq, GRK2, and arrestin3 to the M3R using FRET-based approaches. The specific aim and, at the same time, difficulty was to analyze these events, which are ultimately interconnected, in a mutually independent way in order to differentiate between their conformational requirements. Specifically, there are two main dependencies that need special attention: (1) recruitment of GRK2 to the plasma membrane is enhanced by Gαq and Gβγ subunits upon G protein activation; and (2) arrestin binding to receptors requires both GRK-mediated phosphorylation of receptors and receptor activation. Our experimental design accounted for both processes.
(1) We measured GRK2 recruitment under conditions of complete Gi/o and Gq inhibition by preincubation with PTX and FR900359 (Figure 2A and B). To circumvent G protein-mediated membrane recruitment of GRK2, we included a CAAX-motif at the C-terminus of the fluorophore to induce isoprenylation (Inglese et al., 1992), and in addition a D110A mutation to attenuate interaction with Gαq (Wolters et al., 2015). However, we had to omit coexpression of Gq in FR900359-treated cells, otherwise the interaction between M3R and GRK2 was completely blocked (Figure 2—figure supplement 1C), presumably due to occupation of receptors by inactivated G proteins, and thereby blockade of the interaction site for GRK2 binding. Strikingly, the rank order of agonist efficacies to recruit GRK2 was the same, independent of whether G proteins were inhibited or not. However, the overall FRET amplitudes and also the agonist responses normalized to the one of ACh were considerably larger and their overall kinetics considerably slower in the absence of G protein inhibition, demonstrating the contribution of G proteins to translocating GRK2 to the plasma membrane. Importantly, the rank order of the efficacy of Metha, Are, and Pilo to recruit GRK2 to the M3R (Figure 2A and B, Figure 2—figure supplement 1A and B) is the same as observed for M3R phosphorylation determined by Butcher et al. (Figure 6 in Butcher et al., 2011), indicating that GRK2 recruitment to the M3R correlates with functional phosphorylation.
(2) To avoid an influence of the different efficacies of agonists when recruiting GRK2 and thus the receptor being phosphorylated, we established a specific single-cell FRET protocol. We confirmed arrestin recruitment to the M3R to be dependent on GRK2-mediated phosphorylation (Figure 2—figure supplement 2), as in the absence of coexpression of GRK2 or upon expression of the dominant negative GRK2(K220R), a decrease of the arrestin recruitment was detected as expected based on previous observations at the β2-adrenergic receptor (Krasel et al., 2005). All cells were preexposed to the full agonist ACh for a defined period of time, in order to allow for phosphorylation of receptors before ACh was withdrawn, followed by the measurement of arrestin recruitment induced by the seven different agonists (Figure 2E and F) to the thus prephosphorylated receptors. As the dephosphorylation of receptors typically occurs on the order of at least several minutes, our experimental design allowed the comparison of agonist efficacies to prephosphorylated receptors, minimizing any influence of agonist efficacy towards G protein activation and GRK2 binding. Our results show that unlike for G protein and GRK2 recruitment, Pilo induces significantly more arrestin3 recruitment compared to Are. These results are in line with a potential bias of Pilo towards arrestin that was reported by Pronin et al., 2017. However, here we were able to observe that the bias towards arrestin recruitment was not introduced by a bias towards GRK2 recruitment. In contrast, no detectable bias was observed between Gq recruitment and GRK2 recruitment for any of these agonists. Our findings disentangle bias between Gq, GRK2, and arrestin3 recruitment to the M3R. Arrestin3 recruitment was tested in an overexpression system for all agonists under exactly the same conditions regarding prephosphorylation and kinase expression. The influence of endogenous GRK and arrestin isoforms is therefore negligible. However, we cannot rule out that the overall agonist-dependent recruitment might as well depend on the expression pattern of GRKs. Since M3R-mediated signaling is regulated by distinct mechanisms (Luo et al., 2008), our results represent the first step on the way to completely disentangle bias signaling. Measurement of agonist-induced G protein recruitment to receptors by means of FRET is a much more direct way to study the efficacy of agonists to activate receptors compared to the measurement of G protein activation or downstream signaling due to its lack of signal amplification, which typically leads to the phenomenon of signal saturation at submaximal levels of receptor activation. However, our determination of the efficacy of the different agonists to activate G proteins by calculating resulted in the same rank order for the agonists, confirming the applicability of receptor G protein interaction assay. We determined also saturating agonist concentration for our GRK2 and arrestin3 recruitment measurements. These concentrations were used for our measurement of agonist efficacies, and, to ensure comparable conditions, to measure G protein binding as well. Having taken very serious care of measuring all three binding events completely independent of each other, a comparison of the agonist efficacies for all three events was possible. Metha and Fur recruited arrestin3, GRK2, and Gq to a similar extent compared to the full agonist ACh, whereas the partial agonists Are, Guva, and Pilo differed substantially in their relative efficacy to recruit arrestin3 versus GRK2 and Gq (Figures 1B, 2B and F). The strongest deviation in the rank order was observed for Pilo, which from the set of partial agonists was the poorest GRK2 and Gq recruiter, but the best arrestin3 recruiter, indicating a relative bias towards arrestin3 binding. Mda did not show a detectable effect in the radioligand replacement assay and in the recruitment of (D110A)GRK2-CAAX. Furthermore, it showed only very low levels of recruitment with respect to Gq, low level of recruitment with respect to arrestin3, and a medium level of Gq activation and GRK2 recruitment in the presence of intact Gq. Therefore, it is difficult to ascertain its precise mechanism of action. Even though we did not directly measure GRK2-mediated phosphorylation, we found that the slow component of arrestin3 recruitment upon first stimulation with agonist was dependent on the expression of catalytically active GRK2, similar to what has been described for β2-adrenergic receptors (Krasel et al., 2005).
However, this slow component, even though it is smaller in amplitude, was also visible during subsequent applications of agonist. This could either indicate incomplete phosphorylation at the time points of the repetitive agonist applications or could reflect phosphorylation-independent effects, such as the interaction of arrestin3 with the membrane and a thus prolonged residence time of arrestin3 at the membrane after the first recruitment to the receptor (Lally et al., 2017). Our finding that agonist bias with respect to the ability to activate G proteins versus the ability to induce arrestin recruitment to the receptor can be detected correlates nicely with recent advancements in the determinnation of conformational differences of the angiotensin-II receptor type 1 when bound to Gq or arrestin (Wingler et al., 2020; Suomivuori et al., 2020). Furthermore, our results demonstrating the existence of functional agonist-dependent patterns of differential recruitment of arrestin3 and GRK2 to the M3R are supported by a recent study demonstrating agonist-dependent patterns of recruitment of G proteins, GRK2, and a conformationally selective nanobody directed against the intracellular cavity of activated opioid receptors (Stoeber et al., 2020). The finding that different receptor ligands can differentially affect binding of G proteins, arrestins, and GRKs opens up the potential for future drug development to specifically direct signaling in one or the other direction. In summary, our findings highlight the existence of ligand-induced bias at the M3R and the importance of understanding GRK2 recruitment and its role for the subsequent arrestin3 recruitment in order to fully differentiate between a bias in GRK2 recruitment and a bias in arrestin3 recruitment. While developed for the M3R, the methodology is generally applicable to all GPCRs, G proteins, GRKs, and arrestins, and therefore offers new possibilities to disentangle biased effector pathway engagement at the level of effector protein recruitment and activation.
Materials and methods
Reagent type (species) or resource | Designation | Source or reference | Identifiers | Additional information |
---|---|---|---|---|
Chemical compound, drug | Acetylcholine iodide | Sigma-Aldrich | CAS number: 2260-50-6 | |
Chemical compound, drug | Arecoline hydrobromide | TCI Chemicals | CAS number: 300-08-3 | |
Chemical compound, drug | Methacholine chloride | TCI Chemicals | CAS number: 62-51-1 | |
Chemical compound, drug | Guvacoline hydrobromide | TRC Canada | CAS number: 17210-51-4 | |
Chemical compound, drug | Methyl-3-(dimethylamino) propionate | Sigma-Aldrich | CAS number: 3853-06-3 | |
Chemical compound, drug | Pilocarpine hydrochloride | TCI Chemicals | CAS number: 54-71-7 | |
Chemical compound, drug | Trimethyl-(5-methyl-furan-2-ylmethyl)-ammonium iodide | Sigma-Aldrich | CAS number: 1197-60-0 |
Reagents
Coelenterazine h was obtained from NanoLight, Pinetop, AZ. Dulbecco’s Modified Eagle’s Medium (DMEM), PBS, penicillin/streptomycin, and Trypsin-EDTA were from Capricorn Scientific Gmbh, Ebsdorfergrund, Germany. Poly-L-lysine hydrobromide, FCS, L-glutamine, PEI, acetlycholine iodide, methyl-3-(dimethylamino)propionate (Mda), and trimethyl-(5-methyl-furan-2-ylmethyl)-ammonium iodide (Fur) were obtained from Sigma-Aldrich, Merck KGaA, Darmstadt, Germany. Arecoline hydrobromide, pilocarpine hydrochloride, and methacholine chloride were from TCI Chemicals, Eschborn, Germany, and guvacoline hydrobromide from TRC Canada, Toronto, Canada. Pertussis toxin was purchased from EMD Millipore Corp., Merck KGaA. FR900359 was isolated from Ardisia crenata leaves as previously described in Schrage et al., 2015.
Plasmids
cDNAs for Gαq, Gαq-yellow fluorescent protein (YFP), where YFP was inserted between residues 124 and 125 of Gαq (Hughes et al., 2001), Gβ1, Gγ2 (Bünemann et al., 2003), Gαo, mTurquoise2-Gγ2, where mTurquoise2 was fused to the N-terminus of Gγ2 (Jelinek et al., 2021), M3R-YFP, where YFP was fused C-terminally to M3R (Hoffmann et al., 2012), GRK2 (Winstel et al., 1996), GRK2-mTurquoise, and GRK2-YFP, where mTurquoise or YFP were fused to the GRK2 C-terminus (Wolters et al., 2015), and arrestin3-mTurquoise, where mTurquoise was fused to the C-terminus of arrestin3 (Krasel et al., 2005; Miess et al., 2018; Roseberry et al., 2001), were described previously. M3R-mTurquoise was cloned by an exchange of mCit to mTurquoise of M3R-mCit (Jelinek et al., 2021), where mTurquoise was fused at the C-terminus. The cDNA for the M3R was obtained from the Missouri S&T cDNA Resource Center and pcDNA3 from Invitrogen. cDNA containing pNluc was kindly provided by Dr. N. Lambert (Augusta Medical College, GA), and pNluc-Gγ2, where Nluc was fused to the N-Terminus of GGγ2, was kindly provided by Dr. C. Krasel (Philipps-Universität Marburg, Germany). GRK2-YFP-CAAX was cloned by amplifying the cDNA for YFP or mTurquoise by PCR with the forward primer AAA AAA TCT AGA GTG AGC AAG GGC GAG G and the reverse primerAAAAAAGCGGCCGCCTAggagagcacacacttgcagctcatgcagcccgggccactctcatcaggagggttCTTGTACAGCTCGTCCATGC. The reverse primer attaches the last 18 amino acids of H-Ras to the C-terminus of the fluorescent protein. The resulting PCR product was digested with XbaI and NotI and cloned into pcGRK2-YFP (Wolters et al., 2015) that had been digested with the same enzymes, replacing the YFP with the modified fluorescent protein. (D110A)-GRK2-YFP-CAAX was generated from GRK2-YFP-CAAX by mutagenesis analogously to (D110)-GRK2-YFP in Wolters et al., 2015 using the primer: CCG GGA GAT CTT CGC CTC ATA CAT CAT G.
Cell culture and transfection
Request a detailed protocolFor experiments HEK tsA 201 cell line was used, which was a kind gift from the Lohse laboratory, University of Würzburg. Therefore, all experiments were carried out in HEK293T cells. They were cultured in DMEM with high glucose and supplemented with 10% FCS, 2 mM L-glutamine, 100 U/mL penicillin, and 0.1 mg/mL streptomycin. Cells were transiently transfected 24 hr after seeding (6 cm dish) with linear polyethylenimine (PEI) 25 kDa as transfecting agent. For M3R-G protein interaction experiments, cells were transfected with cDNAs as follows: 1.5 µg M3R -YFP, 2.4 µg Gαq, 0.75 µg Gβ1, and 0.3 µg mTurquoise2-Gγ2. For M3R GRK2 interaction experiments, cells were transfected with the following cDNAs: 1.5 µg M3R-YFP, 2.4 µg Gαq, 0.75 µg Gβ1, 0.3 µg Gγ2, and 0.75 µg GRK2-mTurquoise. For M3R (D110A)-GRK2-YFP-CAAX interaction experiments, the cells were transfected with 1.5 µg M3R-mTurquoise, 2.4 µg pcDNA3, 0.75 µg Gβ1, 0.3 µg Gγ2, and 0.75 µg (D110A)-GRK2-YFP-CAAX. For M3R GRK2 interaction experiments without Gαq overexpression, the cells were transfected with 1.5 µg M3R-mTurquoise, 2.4 µg pcDNA3, 0.75 µg Gβ1, 0.3 µg Gγ2, and 0.75 µg GRK2-YFP. For M3R-arrestin interaction experiments, cells were transfected with the following amounts of cDNAs: 1.5 µg M3R-YFP, 0.75 µg GRK2, and 1.5 µg arrestin3-mTurquoise. For Gβ1/Gγ2 GRK2 interaction experiments, the cells were transfected with: 1.5 µg M3R, 2.4 µg pcDNA3, 0.75 µg Gβ1, 0.3 µg mTurquoise2-Gγ2, and 0.75 µg GRK2-YFP and at the M2R with 1.5 µg M2R 2.4 µg Gαo, 0.75 µg Gβ1, 0.3 µg mTurquoise2-Gγ2, and 0.75 µg GRK2-YFP. For BRET-based Gαq activation experiments, cells were transfected with cDNAs as follows: 1.5 µg M3R, 2.4 µg Gαq-YFP, 0.75 µg Gβ1, 0.3 µg pNuc-Gγ2, and 0.75 µg GRK2. The mixing ratio of PEI to DNA was 3:1. Per 1 µg DNA, 50 µL DMEM without FCS were added to the DNA and PEI solutions, respectively. Both solutions were mixed, added together and incubated, protected against light, at 20 °C for 30 min. The mix was added to the cells and incubated at 37°C in a humidified atmosphere of 95% air and 5% CO2. For FRET-based experiments, cells were seeded onto 6-well plates with 25 mm coverslips coated with poly-L-lysine after 24 hr. For BRET-based experiments, cells were counted after 24 hr and 16,000 cells/well were seeded in 96-well plates (Greiner 96 Flat White) coated with poly-L-lysine. After 48 hr of transfection, measurements were performed at room temperature.
Single-cell FRET-based measurements
Request a detailed protocolUnless indicated otherwise, FRET-based measurements on single cells were performed as described in Milde et al., 2013. Transiently transfected cells were subjected to single-cell time-resolved FRET imaging with constant superfusion of either buffer (137 mM NaCl, 5.4 mM KCl, 2 mM CaCl2, 1 mM MgCl2, 10 mM HEPES, pH 7.3) or buffer-containing agonist. Cells were measured with an inverted fluorescence microscope (Eclipse Ti; Nikon, Düsseldorf, Germany). Light-emitting diodes (LED) at 425 nm and 500 nm were used for excitation. The intensity of both LEDs of the excitation system (pE-2; CoolLED, Andover, UK) was set to 2%. Fluorescence intensity was measured using the imaging software NIS-Elements advanced research (Nikon Corporation) and recorded at 2 Hz. The FRET ratio was calculated as the fluorescence intensity ratio of YFP and mTurquoise emission after mTurquoise excitation at 425 nm. Stimulation with agonist led to an increase in fluorescence intensity ratio, reflecting the interaction of the M3R with arrestin3 or GRK2. The fluorescence data were corrected for background fluorescence, bleed-through, and false excitation, and then plotted over time. The presented data were baseline-corrected to account for photobleaching. Time intervals of every measurement are indicated in the corresponding figure legend. For averaging the single FRET-based measurements, each measurement was normalized to the induced individual maximum response of saturating concentrations of ACh (30 µM concentration) if not indicated otherwise. Amplitudes were determined by calculating the difference between the mean of FRET ratio values 5 s before stimulation with the agonist and 5 s before withdrawal of the agonist. For FRET-based measurements analyzing the increase of receptor-arrestin3 interaction upon repeated stimulation by the same agonist, each cell was stimulated with agonist for 30 s, followed by a stimulation with ACh to make sure that the maximum of the effect was reached. Afterward, another stimulation of 30 s with agonist followed, termed the second stimulation. Pilo and ACh were measured as additional controls in every cell, and the amplitude of the stimulation with ACh was used for normalization. The normalized M3R -arrestin3 interaction level of each agonist was compared between first and second stimulation. For FRET-based experiments using pertussis toxin for Gαi/o inhibition, the cells were pretreated overnight with a resulting concentration of 50 ng/mL PTX. For inhibition of Gαq using FR900359, the cells were preincubated 10 min before the start of the measurement and superfused during the whole measurement as buffer and agonist solutions contained 1 µM FR900359.
BRET-based measurement
Request a detailed protocolBioluminescence resonance energy transfer (BRET)-based measurements were conducted with the luciferase reporter Nluc (Hall et al., 2012). Transiently transfected adherent cells were measured with a Spark 20M Multimode Microplate Reader from Tecan. Fluorescence and luminescence intensities were obtained using the SparkControl application (Tecan). Cells were carefully washed with buffer (as described in single-cell FRET-based measurements). Afterward, coelenterazine h in buffer was added to the cells. Every well contained a volume of 80 µL and a final concentration of 3.07 µM coelenterazine h. The BRET ratio was calculated as the quotient of the YFP signal (light emitted between 520 nm and 700 nm) and the nanoLuc signal (light emitted between 415 nm and 485 nm). After 10 min incubation with coelenterazine h, measurement of the baseline BRET ratio took place for 10 cycles (baseline phase). One cycle had a duration of about 44 s. After 10 cycles, this led to a total of 6.5 min of baseline measurement. The measurement was paused and 20 µL buffer (negative control) or agonist was added to the cells. The development of BRET was measured for 10 cycles (agonist phase) and then paused again for the addition of 20 µL of a solution containing ACh in buffer. The final concentration of ACh in each well was at least 50 µM, aimed at creating saturating conditions and measuring the maximum response obtainable in each individual well. Afterward, 10 more cycles were measured (saturation phase). Every phase was fitted by a line fit, and the change in BRET signal was calculated as the distance of the last time point of the previous phase and the first point of its following phase. Hence, the agonist-induced change in BRET signal was calculated as the distance between the baseline phase and the agonist phase. The additional change induced by saturating concentrations of ACh was calculated as the distance between the agonist phase and the saturation phase. The maximum change in BRET signal was then calculated as the sum of the agonist-induced change and the additional change upon application of saturating concentrations of ACh. The change in BRET signal was normalized to the maximum change in BRET signal for every well.
Radioligand displacement assay
Request a detailed protocolThe human M3R antagonist radioligand receptor binding assay was performed by Eurofins Cerep, Celle-Levescault, France (item 95). Competition between [3H]4-DAMP (0.2 nM) and increasing concentrations of each compound for M3R in CHO whole cells was measured after incubation for 60 min at room temperature. Nonspecific binding was determined in the presence of 1 µM atropine, and 4-DAMP was used as a reference compound. Radioactivity was quantified by scintillation counting. Results are expressed as percent of control specific binding, that is, (specific binding/control-specific binding) × 100 obtained in the presence of the test compound.
Statistics and data analysis
Request a detailed protocolThe radioligand displacement measurements performed by Eurofins are presented as mean ± SD of n experiments with two independent experiments. All other measurements are presented as mean ± SD of n experiments with at least three independent transfections. For statistical analysis, a paired Student’s t-test as well as a one-way ANOVA with Tukey’s multiple comparison test were conducted using GraphPad Prism 6.01. EC50 values were obtained by using GraphPad Prism’s nonlinear regression curve fit for the concentration-response curves, fitting the Hill equation with four parameters, variable slope, and by setting the bottom parameter to zero. For the radioligand displacement assay, IC50 values were obtained from the inhibition curves and Ki values were calculated using the Cheng and Prusoff method (Cheng, 2001).
For the evaluation of the offset kinetics, each cell was normalized to the maximum response of agonist, and a one-phase exponential decay for the treatment with ACh and a two-phase decay for the treatment together with the inhibitors were fitted using GraphPad Prism 8.3.0. The following equation was used for fitting the one-phase decay:
.
The following equation was used for fitting the two-phase decay:
where Y0 is the y value when x is zero and Plateau is the y value at infinite times. For both conditions, the Y0 was constrained to 1 and the Plateau to 0.
The concentration-response data for Gq activation, arrestin3 recruitment, and GRK2 recruitment to the M3R were analyzed using the operational model of agonism by Black and Leff, 1997. The following equation was used for fitting and calculation of τ:
.
where is the molar concentration of agonist, is the theoretical maximum response of the system, is the equilibrium dissociation constant of the agonist-receptor complex, represents a slope factor, and τ is the operational factor of efficacy. The following equation was implemented in GraphPad Prism:
where is the log of the molar concentration of agonist [A]. Each value was set individually to the value determined by the competition binding experiment for each agonist. Because the response of GRK2 and arrestin3 recruitment for each agonist was normalized to the maximum response of ACh, was set to 1. For Gq activation, where each agonist concentration was normalized to the maximum response of ACh in the same well, was set to the maximum response of ACh, to 98.63. One value was determined for all agonists, and and τ were constrained to be greater than 0. The fitted value of τ was used to calculate for each agonist. ACh was used as a reference compound that all other agonists were scaled to, thus allowing the calculation of .
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
References
-
Operational models of pharmacological agonismProceedings of the Royal Society of London. Series B. Biological Sciences 220:141–162.https://doi.org/10.1098/rspb.1983.0093
-
Differential G-protein-coupled receptor phosphorylation provides evidence for a signaling bar codeThe Journal of Biological Chemistry 286:11506–11518.https://doi.org/10.1074/jbc.M110.154526
-
A Pluridimensional View of Biased AgonismMolecular Pharmacology 90:587–595.https://doi.org/10.1124/mol.116.105940
-
Chemokines and Chemokine Receptors: Positioning Cells for Host Defense and ImmunityAnnual Review of Immunology 32:659–702.https://doi.org/10.1146/annurev-immunol-032713-120145
-
Visualization of a functional Galpha q-green fluorescent protein fusion in living cells Association with the plasma membrane is disrupted by mutational activation and by elimination of palmitoylation sites, but not be activation mediated by receptors or AlF4The Journal of Biological Chemistry 276:4227–4235.https://doi.org/10.1074/jbc.M007608200
-
Biased Receptor Signaling in Drug DiscoveryPharmacological Reviews 71:267–315.https://doi.org/10.1124/pr.118.016790
-
A beta-adrenergic receptor kinase dominant negative mutant attenuates desensitization of the beta 2-adrenergic receptorThe Journal of Biological Chemistry 269:13084–13087.
-
Beta-arrestin binding to the beta2-adrenergic receptor requires both receptor phosphorylation and receptor activationThe Journal of Biological Chemistry 280:9528–9535.https://doi.org/10.1074/jbc.M413078200
-
C-edge loops of arrestin function as a membrane anchorNature Communications 8:14258.https://doi.org/10.1038/ncomms14258
-
M3 muscarinic acetylcholine receptor-mediated signaling is regulated by distinct mechanismsMolecular Pharmacology 74:338–347.https://doi.org/10.1124/mol.107.044750
-
Fine-tuning of GPCR activity by receptor-interacting proteinsNature Reviews. Molecular Cell Biology 10:819–830.https://doi.org/10.1038/nrm2803
-
Structural mechanism of arrestin activationCurrent Opinion in Structural Biology 45:160–169.https://doi.org/10.1016/j.sbi.2017.05.001
-
The experimental power of FR900359 to study Gq-regulated biological processesNature Communications 6:1–17.https://doi.org/10.1038/ncomms10156
-
Visual pigment: G-protein-coupled receptor for light signalsCellular and Molecular Life Sciences 54:1299–1315.https://doi.org/10.1007/s000180050256
Decision letter
-
William I WeisReviewing Editor; Stanford University School of Medicine, United States
-
Richard W AldrichSenior Editor; The University of Texas at Austin, United States
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Decision letter after peer review:
Thank you for submitting your article "Disentangling bias between Gq activation, GRK2 recruitment, and arrestin3 recruitment at the M3R" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Richard Aldrich as the Senior Editor. The reviewers have opted to remain anonymous.
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
As the editors have judged that your manuscript is of interest, but as described below that additional experiments are required before it is published, we would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). First, because many researchers have temporarily lost access to the labs, we will give authors as much time as they need to submit revised manuscripts. We are also offering, if you choose, to post the manuscript to bioRxiv (if it is not already there) along with this decision letter and a formal designation that the manuscript is "in revision at eLife". Please let us know if you would like to pursue this option. (If your work is more suitable for medRxiv, you will need to post the preprint yourself, as the mechanisms for us to do so are still in development.)
Summary:
This paper uses fluorescent probes attached to different the M3 muscarinic receptor (M3R) interacting proteins Gq, arrestin and GRK2 in cellular BRET assays to examine bias among different M3R ligands. Bias is typically discussed in terms of preferential activation of the partner G protein vs. arrestin. GRK recruitment/activity is important for arrestin binding and is expected to influence arrestin-recruiting activity of a particular ligand. However, the possibility of a distinct agonist-dependent receptor conformation on GRK association has not been investigated. This study attempts to answer this question.
Essential revisions:
While the work addresses a very important problem and is technically excellent, all of the reviewers agree that the data do not decisively establish that the receptor adopts a conformation distinct from that which binds arrestin or G protein. As you will see in points 1-3 below, the major concern is that the experiments have not completely disentangled GRK and G protein interactions. If you are able to obtain data that convincingly shows a GRK interaction without the influence of G proteins (which would require mutants or some other way of getting rid of/maximally activating G proteins), then the work would provide new and important information regarding bias between GRK- and arrestin-preferring conformations. The reviewers recognize that this may be very hard to do, but there are some suggestions along these lines as outlined below.
1. The main difficulty encountered by the authors relates to the fact, as the same group showed convincingly in Wolters et al., 2015, that agonist-dependent FRET between active M3 receptors and GRK2 depends almost entirely on the interactions of GRK2 with Gqα and betagamma subunits, which translocate GRK2 to the plasma membrane. The authors make the argument that under the conditions of their experiments Gq is maximally activated by all of the agonists they used, therefore any differences in GRK2 recruitment efficacy must be due to differences in the direct interaction between GRK2 and activated receptors. Setting aside possible differences in G protein activation due to the presence of inserted fluorescent tags and the different preparations this is a reasonable argument. However, there are a number of alternative approaches that could be used to validate the findings. Adding a membrane anchor to GRK2, overexpressing G betagamma subunits alone, blocking Gq activation with small molecules or gene deletion, maximally activating Gq with GTPgammaS, or a combination of these approaches, would greatly increase confidence that the result is fundamentally correct. Since the authors observe essentially no bias between Gq activation and GRK2 recruitment with the ligands they tested (Figure 7C), this could be because (as the authors conclude) the receptor conformations that activate Gq are the same as those that recruit GRK2, or alternatively that the GRK2 signals they measure still reflect primarily the efficiency of Gq activation. This is indeed a very important question because of the possibility that arrestin bias is actually a reflection of GRK bias, therefore the equivalency of G protein activation and GRK recruitment (the main novel conclusion of the paper) deserves very stringent testing.
2. GRK2 also binds to activated Galphaq subunits, which may have an effect on the kinetics of G protein activation and/or deactivation (along with the Gbg binding). This seems to be the case because the presence of overexpressed GRK2 left-shifts the G protein activation dose response curves in Figure 3.
3. Comparison of the data from Figure 3A with Figure 7B, which is used by the authors to argue that they don't need to worry about G protein-dependent effects on GRK2 recruitment seems invalid because GRKs likely have less affinity for phosphorylated receptors (Pulvermuller et al. 1993), as is the case in 7B, and they would compete poorly with arretin-bound GPCRs (endogenous or not). Regardless, also not sure if it is valid to say that there should be no coupling between G protein activation efficacy versus GRK recruitment because their EC50 curves are not similar. The affinity is what it is, and could still be dependent the amount of G protein subunits liberated. Issues like this is what makes interpreting this data very complicated.
4. How do the expression levels of tagged proteins compare to the endogenous levels? Do the M3R-YFP receptors used in the GRK2 experiments express at the same level as the unlabeled M3R receptors used in the G protein activation experiments? Activating the Gq BRET constructs with M3R-YFP might also be useful to address this point. Can the authors comment on how GRK2 overexpression increases the potency of Gq activation?
5. The authors have not tested or understood how the endogenous GRK and arrestin isoforms might be adding to or competing with the tagged proteins, raising the possibility that bias is instead for a specific GRK or arrestin 2 versus arrestin 3. The activated M3R is phosphorylated by several GRKs other than GRK2, such as GRK3 and GRK6 (J Luo et al. M3 muscarinic acetylcholine receptor-mediated signaling is regulated by distinct mechanisms. Mol. Pharmacol. 2008). This latter work was also carried out with HEK293 cells. The Luo et al. paper needs to be cited, and the results of this paper need to be discussed in the context of the new data. How does the Luo et al. study affect the conclusions here?
6. Is there anything known that different muscarinic agonists recruit different GRKs to the receptor, leading to different phosphorylation patterns? This could explain the observed discrepancies in rank order of agonist potencies (GRK2 vs. arrestin3 recruitment). No attempt is made to verify that the M3R is getting phosphorylated at GRK2 specific sites.
7. The potential implications of this study for drug development efforts should be discussed in more detail.
8. The times that people normally cite for full GRK phosphorylation of GPCRs is on the order of minutes unless one is talking about rhodopsin. The full activation of the "GRK" interaction shown after application of agonist suggests that this is really a GRK-G protein interaction, not a GRK-GPCR interaction. If so, how would this affect interpretation? It would have been good to see how a non-Gbg dependent GRK would respond in this assay.
[Editors' note: further revisions were suggested prior to acceptance, as described below.]
Thank you for submitting your article "Disentangling bias between Gq, GRK2 and arrestin3 recruitment at the M3 muscarinic acetylcholine receptor" for consideration by eLife. Your article has been reviewed by 2 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Richard Aldrich as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: John JG Tesmer (Reviewer #2); Nevin A Lambert (Reviewer #3).
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
Essential revisions:
The reviewers agree that the paper has been greatly strengthened by the new experimental data. They do agree, however, that the important message of the paper is weakened by inclusion of the "mixed signal" data in figures 2 and 4A,B (see Reviewer 3 comments below), which by their nature do not cleanly decouple signals from different sources. These data do support the separation of direct receptor vs. G protein-dependent signals and emphasize the importance of the decoupled assays in Fig. 4C,D , so we ask that the paper be revised with those data provided as supplemental rather than main figures. You may also want to consider commenting on the kinetics mentioned by reviewer 3.
Reviewer #2:
I am impressed with the robust experimental response to the two reviews and I think it is an example where things did dramatically improve in terms of rigor. I am satisfied with their response to my comments. I think this is a very interesting comparative study and I cannot think of many better ways to disentangle how bias manifests itself at the M3R. One could quibble onwards because of the complexities of the system and the fact that overexpressed proteins were used, but I think this data should be seen by the community at large now.
Reviewer #3:
The authors have done a very nice job in addressing a most difficult problem, i.e. separating G protein-dependent and -independent GRK2 signals. They did this with a CAAX anchor and PTX/FR900359. I still suggest that Figures 1 and 2 are not necessary and do not add to the story, and in fact confuse the situation by showing "M3R-GRK2 interaction" that is clearly a mixed response that includes bystander FRET- the signal is largely free Gbetagamma-GRK, not receptor-GRK. Figure 4A and 4B are similarly confusing and unnecessary- why show FRET signals that are not interpretable?
The results in 4C and 4D are key, and should take center stage. The authors should also point out that the kinetics of the isolated G protein-independent GRK signals are fast, as would be expected for what is essentially a reporter of receptor conformation. In contrast, the mixed GRK signals in 4A, 4B and Figure 2 have slow decays, presumably because this is rate-limited by G protein deactivation (GTP hydrolysis). This kinetic difference supports their claim that they've successfully isolated the receptor-GRK interaction from G proteins. Along these same lines, it's worth noting that the presumed bystander signal after PTX/FR900358 in Figure 4S1 C is also fast; the slower components due to G proteins get blocked by the PTX/FR. Indeed, the decay of the untreated trace in this figure could be fitted to two exponentials, the slower of the two would be consistent with GTP hydrolysis, the faster with receptor deactivation. A brief kinetic analysis of this sort would strengthen the paper.
https://doi.org/10.7554/eLife.58442.sa1Author response
Essential revisions:
While the work addresses a very important problem and is technically excellent, all of the reviewers agree that the data do not decisively establish that the receptor adopts a conformation distinct from that which binds arrestin or G protein. As you will see in points 1-3 below, the major concern is that the experiments have not completely disentangled GRK and G protein interactions. If you are able to obtain data that convincingly shows a GRK interaction without the influence of G proteins (which would require mutants or some other way of getting rid of/maximally activating G proteins), then the work would provide new and important information regarding bias between GRK- and arrestin-preferring conformations. The reviewers recognize that this may be very hard to do, but there are some suggestions along these lines as outlined below.
1. The main difficulty encountered by the authors relates to the fact, as the same group showed convincingly in Wolters et al., 2015, that agonist-dependent FRET between active M3 receptors and GRK2 depends almost entirely on the interactions of GRK2 with Gq α and betagamma subunits, which translocate GRK2 to the plasma membrane. The authors make the argument that under the conditions of their experiments Gq is maximally activated by all of the agonists they used, therefore any differences in GRK2 recruitment efficacy must be due to differences in the direct interaction between GRK2 and activated receptors. Setting aside possible differences in G protein activation due to the presence of inserted fluorescent tags and the different preparations this is a reasonable argument. However, there are a number of alternative approaches that could be used to validate the findings. Adding a membrane anchor to GRK2, overexpressing G betagamma subunits alone, blocking Gq activation with small molecules or gene deletion, maximally activating Gq with GTPgammaS, or a combination of these approaches, would greatly increase confidence that the result is fundamentally correct. Since the authors observe essentially no bias between Gq activation and GRK2 recruitment with the ligands they tested (Figure 7C), this could be because (as the authors conclude) the receptor conformations that activate Gq are the same as those that recruit GRK2, or alternatively that the GRK2 signals they measure still reflect the primarily the efficiency of Gq activation. This is indeed a very important question because of the possibility that arrestin bias is actually a reflection of GRK bias, therefore the equivalency of G protein activation and GRK recruitment (the main novel conclusion of the paper) deserves very stringent testing.
We thank the reviewer to bring up this important point. We spent a lot of time to clearly address the receptor-GRK2 interaction without the influence of activated G proteins. Since the experiments of M3R-GRK2 interaction by fully activated G proteins in presence of GTPyS, as suggested by the reviewers, did not show reliable alterations in FRET signals, we finally combined numerous approaches. To bypass the G protein dependent translocation (Wolters et al) of GRK2 to the plasma membrane, we introduced a CAAX-box at the C-terminus of the GRK2. Furthermore, we mutated the GRK2 at position D110 to an alanine in order to attenuate Gaq binding (Wolters et al). In addition, we performed the experiments without overexpression of Gaq and preincubated the cells with the Gi/o inhibitor PTX and the Gq inhibitor FR900365 to inhibit endogenous G proteins that could potentially influence GRK2 (Figure 4C,D, supplemental Figure 1A,B,C.). Thus, we removed all possibilities of GRK2-G protein interaction by a combination approach, as indicated by the reviewer. Importantly, even though the agonist-evoked FRET responses under complete Gi and Gq inhibition were markedly reduced, we were still able to resolve agonist-specific interactions between (D110A)GRK2-CAXX and M3R shown in Figure 4 C, D. Therefore, we have specifically measured GRK2 interaction with M3R without any influence of G proteins. Nevertheless, the rank order of agonist efficacy for the GRK2-M3R interaction is the same as for G protein interactions, as shown in Figure 4D.
2. GRK2 also binds to activated Galphaq subunits, which may have an effect on the kinetics of G protein activation and/or deactivation (along with the Gbg binding). This seems to be the case because the presence of overexpressed GRK2 left-shifts the G protein activation dose response curves in Figure 3.
We agree with the reviewer that GRK2 also binds to activated Gaq, which was already shown in Wolters et al., 2015. After our control measurements regarding point 4 (different expression levels of wildtype receptors and YFP-tagged receptors), we decided to measure the direct Gq binding to the M3 receptor, determined by single cell FRET assays, instead of G protein activation (Figure 3). The G protein binding assay allows us to much more directly measure agonist efficacy. Furthermore, as described in response to point 1, we have now measured GRK2 interaction with the M3R in the absence of G protein activation, making it unnecessary to tease apart the exact reason of the left shift of the concentration response curves for of Gq protein activation by GRK2.
3. Comparison of the data from Figure 3A with Figure 7B, which is used by the authors to argue that they don't need to worry about G protein-dependent effects on GRK2 recruitment seems invalid because GRKs likely have less affinity for phosphorylated receptors (Pulvermuller et al. 1993), as is the case in 7B, and they would compete poorly with arretin-bound GPCRs (endogenous or not). Regardless, also not sure if it is valid to say that there should be no coupling between G protein activation efficacy versus GRK recruitment because their EC50 curves are not similar. The affinity is what it is, and could still be dependent the amount of G protein subunits liberated. Issues like this is what makes interpreting this data very complicated.
We totally agree with the reviewer that in the previous version of the manuscript we were not able to completely disentangle G protein activation and GRK2 recruitment, which was also shown in the measurement of different expression levels of wildtype receptor and YFP-tagged receptor (point 4), as suggested by the reviewers. However, in the revised version of our manuscript we have solved this issue by performing GRK2 recruitment assays under conditions of complete G protein inhibition (Figure 4 C, D) (as described above in point 1) to exclude the influence of activated G proteins on the GRK2 recruitment. In addition, we now added the M3R Gq-protein interaction assay as a direct measurement of agonist efficacy, which allows for more consistent comparison of efficacies. These major revisions allow for a much better disentanglement of the influence of G proteins on GRK2 and vice versa. We are not speculating about possible differential dynamics of GRK2 binding to phosphorylated and unphosphorylated receptors. However, as published previously, there is no major difference of amplitude and kinetics of GRK2 binding and dissociation to M3R, whether a dominant negative GRK2 (K220R) is used or the GRK2-wt (Wolters et al., 2015).
4. How do the expression levels of tagged proteins compare to the endogenous levels? Do the M3R-YFP receptors used in the GRK2 experiments express at the same level as the unlabeled M3R receptors used in the G protein activation experiments? Activating the Gq BRET constructs with M3R-YFP might also be useful to address this point. Can the authors comment on how GRK2 overexpression increases the potency of Gq activation?
We are glad to follow the reviewers’ suggestion to address this important point. We performed BRET-based measurements of Gq activation for M3 wt and M3-YFP (Author response image 1). Due to the right shift in the concentration response curves of M3-YFP relative to M3 wt, we can no longer claim that G proteins have been fully activated in our GRK2 recruitment measurements. Therefore, we performed, as described above (point 1 and 2), a G protein-independent GRK2 recruitment assay and the direct G protein binding to the receptor to directly compare agonist efficacies (Figure 3 vs Figure 4 C-D), which is now the primary data set for comparing bias. That being said, the main conclusions about the different behavior of the different ligands stand also with this data set.

Gq activation of M3 wt and M3-YFP was measured with a BRET-based assay in HEK293T cells.
(A) Concentration-response curves of Gq activation with GRK2 overexpression were plotted for M3R wt and M3R-YFP. (B) M3R wt and (C) M3R-YFP Gq activation were measured with distinct concentrations of every agonist as indicated. After stimulating every well with agonist, a saturating concentration of ACh (50 µM) was applied. The BRET ratio after agonist application was normalized to the maximum response of ACh (50 µM) of every well. All data are plotted as mean values ± SD for each condition.
As the rank order of agonist efficacies to activate Gq-proteins is the same for the tagged versus untagged receptor (Author response image 1), and the G protein activation data all moved to the supplemental figures we decided to show this dataset only to the reviewers. If requested, we could show this dataset also in the supplemental figures.
Regarding the reviewers’ question “Can the authors comment on how GRK2 overexpression increases the potency of Gq activation?” we admit that this is a bit complicated: As described by Wolters et al., the FRET decrease is markedly enhanced due to simultaneous binding of Gaq and Gbg to GRK2 and in addition the recovery of the FRET between Gaq and Gbg after agonist withdrawal is a bit slowed. This could lead to an apparent left shift of the concentration-response curve for Gq activation. However, this issue is of much less importance for our conclusions due to the implementation of the direct measurement of Gq interaction with M3R. Therefore, we prefer to not discuss this issue in the manuscript.
5. The authors have not tested or understood how the endogenous GRK and arrestin isoforms might be adding to or competing with the tagged proteins, raising the possibility that bias is instead for a specific GRK or arrestin 2 versus arrestin 3. The activated M3R is phosphorylated by several GRKs other than GRK2, such as GRK3 and GRK6 (J Luo et al. M3 muscarinic acetylcholine receptor-mediated signaling is regulated by distinct mechanisms. Mol. Pharmacol. 2008). This latter work was also carried out with HEK293 cells. The Luo et al. paper needs to be cited, and the results of this paper need to be discussed in the context of the new data. How does the Luo et al. study affect the conclusions here?
We agree with the reviewer that we have not tested how the endogenous GRK and arrestin isoforms add or compete with the tagged proteins, but the influence on our measurements of endogenous isoforms should be negligible through the overexpression of our proteins. Moreover, Hosey et al. (1995) showed that GRK2 and GRK3 were the most effective kinases at the M3 receptor. Additionally, Luo et al. (2008) suggested that it is primarily the GRK2 that regulates the M3 receptor through its ability to interact with Gq. With our data we can assume that the GRK2-M3 receptor interaction is essentially based on the regulation via G proteins. The difference between agonist efficacies of wildtype GRK2-M3 receptor interaction (Figure 4 B) and G protein independent GRK2-M3 receptor interaction (Figure 4 D) show the impact of G proteins on the translocation of the GRK2 to the plasma membrane. It was also shown in Luo et al. that arrestin-3 appears to play a larger role in terminating signaling in response to agonist exposure at the M3 receptor compared to arrestin-2. (We included a sentence indicating this fact: “We used in our study GRK2 and arrestin-3 isoforms, which appear to be the most regulating isoforms at the M3 receptor as shown in Luo et al.”). We discuss this issue and the respective limitations in the Discussion section.
6. Is there anything known that different muscarinic agonists recruit different GRKs to the receptor, leading to different phosphorylation patterns? This could explain the observed discrepancies in rank order of agonist potencies (GRK2 vs. arrestin3 recruitment). No attempt is made to verify that the M3R is getting phosphorylated at GRK2 specific sites.
We thank the reviewer for bringing up this point. It was shown in Butcher et al. (2011) that different agonists can induce phosphorylation of the M3 receptor preferentially at specific sites. However, in our study we avoided complications by agonist-specific phosphorylation patterns or extent by always pre-exposing the cells to the full agonist acetylcholine, just before we applied one of the specific investigated agonists. Thus, our experimental protocol with single cell FRET recording, which allows for a fast solution exchange, allowed to specifically compare agonist-specific recruitment of arrestin 3 to M3R without interference with effects derived from agonist-specific G protein activation or GRK recruitment (Figure 4 A, C). Thus, even without actually measuring bulk receptor phosphorylation (it is in our opinion impossible to mirror our experimental conditions in a phosphorylation assay), we can conclude that the difference in agonist efficacies is not based on different phosphorylation patterns of the M3 receptor.
7. The potential implications of this study for drug development efforts should be discussed in more detail.
We thank the reviewers for their helpful suggestion. Therefore, we stated in the last part of the discussion: "The finding that different receptor ligands can differentially affect binding of G proteins, arrestins and GRKs opens up the potential for future drug development to specifically direct signaling in one or the other direction."
8. The times that people normally cite for full GRK phosphorylation of GPCRs is on the order of minutes unless one is talking about rhodopsin. The full activation of the "GRK" interaction shown after application of agonist suggests that this is really a GRK-G protein interaction, not a GRK-GPCR interaction. If so, how would this affect interpretation? It would have been good to see how a non-Gbg dependent GRK would respond in this assay.
We agree with the reviewer that our GRK2 recruitment measurements were still affected by activated G proteins. As explained above the new measurements of the G protein independent GRK2–M3 interaction showed reduced overall amplitudes and reduced agonist-induced responses normalized to the one of acetylcholine, indicating the impact of G proteins on GRK2 recruitment (Figure 4 A-D). However, having eliminated any influence from the G protein, we still find the same rank order of agonist efficacies to recruit GRK2 to the receptor.
[Editors' note: further revisions were suggested prior to acceptance, as described below.]
Reviewer #3:
The authors have done a very nice job in addressing a most difficult problem, i.e. separating G protein-dependent and -independent GRK2 signals. They did this with a CAAX anchor and PTX/FR900359. I still suggest that Figures 1 and 2 are not necessary and do not add to the story, and in fact confuse the situation by showing "M3R-GRK2 interaction" that is clearly a mixed response that includes bystander FRET- the signal is largely free Gbetagamma-GRK, not receptor-GRK. Figure 4A and 4B are similarly confusing and unnecessary- why show FRET signals that are not interpretable?
We agree with the reviewer that Figures 1 and 2 are not really necessary and do not add to the story. Therefore, we moved both figures to the supplement, resulting in Figure 1—figure supplement 1 and Figure 1—figure supplement 3. Additionally, we fully agree that it is confusing to show the M3R-GRK2 interaction in Figure 2 and in Figure 4A, B. Thus, we exchanged the example measurement of M3R-GRK2 interaction to M3R-(D110A)GRK2-CAAX interaction shown in the new Figure 1—figure supplement 3B. Figure 4A and B were moved to Figure 2—figure supplement 1 A and B.
The results in 4C and 4D are key, and should take center stage. The authors should also point out that the kinetics of the isolated G protein-independent GRK signals are fast, as would be expected for what is essentially a reporter of receptor conformation. In contrast, the mixed GRK signals in 4A, 4B and Figure 2 have slow decays, presumably because this is rate-limited by G protein deactivation (GTP hydrolysis). This kinetic difference supports their claim that they've successfully isolated the receptor-GRK interaction from G proteins. Along these same lines, it's worth noting that the presumed bystander signal after PTX/FR900358 in Figure 4S1 C is also fast; the slower components due to G proteins get blocked by the PTX/FR. Indeed, the decay of the untreated trace in this figure could be fitted to two exponentials, the slower of the two would be consistent with GTP hydrolysis, the faster with receptor deactivation. A brief kinetic analysis of this sort would strengthen the paper.
We are glad to follow the reviewer’s suggestion to point out the differences in offset-kinetics in presence and inhibition of G proteins. We added a half sentence to the results part stating now on page 4 : “… Under conditions of complete Gq and Gi inhibition (see for control Figure 2—figure supplement 1D,E), amplitudes of agonist-evoked FRET increases were reduce and offset kinetics after agonist withdrawal were substantially faster (compare Figure 2A and Figure 2 supplement 1A). …..” Furthermore, as suggested we analysed the offset-kinetics of the Gβγ/GRK2 interaction with and without the inhibitors (Figure 2—figure supplement 1 E, F and G) and stated in the results: “Kinetic analysis of agonist-evoked FRET between Gβγ and GRK2 with or without pre-treatment with PTX and FR900359 showed fast one phased decays for signals measured with pre-treatment and slower two phased decays for signals measured without pre-treatment (Figure 2—figure supplement 1 F, G).” Finally, we completed the sentence in the discussion: “However, the overall FRET amplitudes and also the agonist responses normalized to the one of ACh were considerably larger and their overall kinetics considerably slower in the absence of G protein inhibition, demonstrating the contribution of G proteins to translocating GRK2 to the plasma membrane.”
Again, we would like to thank the reviewers for their many helpful comments and suggestions, which have indeed substantially improved the manuscript and hope that the revised manuscript can be accepted for publication in eLife.
https://doi.org/10.7554/eLife.58442.sa2Article and author information
Author details
Funding
Deutsche Forschungsgemeinschaft (KO-4095/1-1)
- Peter Kolb
Deutsche Forschungsgemeinschaft (KO4095/4-1)
- Peter Kolb
Deutsche Forschungsgemeinschaft (KO4095/5-1)
- Peter Kolb
Deutsche Forschungsgemeinschaft (FOR 2372 290827466)
- Gabriele König
Deutsche Forschungsgemeinschaft (FOR 2372 290847012)
- Evi Kostenis
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Sandra Engel for technical assistance and Daniel Hilger for critical comments on the manuscript. We thank the German Research Foundation DFG for Emmy Noether fellowship KO4095/1-1, Heisenberg professorship grants KO4095/4-1 and KO4095/5-1 (to PK) and grants for the research group FOR 2372, 290827466 (to GK)/290847012 (to EK).
Senior Editor
- Richard W Aldrich, The University of Texas at Austin, United States
Reviewing Editor
- William I Weis, Stanford University School of Medicine, United States
Publication history
- Received: April 30, 2020
- Accepted: November 18, 2021
- Version of Record published: December 1, 2021 (version 1)
- Version of Record updated: December 7, 2021 (version 2)
Copyright
© 2021, Flöser 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
-
- 823
- Page views
-
- 156
- Downloads
-
- 7
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
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
Mitochondrial biogenesis requires the import of >1,000 mitochondrial preproteins from the cytosol. Most studies on mitochondrial protein import are focused on the core import machinery. Whether and how the biophysical properties of substrate preproteins affect overall import efficiency is underexplored. Here, we show that protein traffic into mitochondria can be disrupted by amino acid substitutions in a single substrate preprotein. Pathogenic missense mutations in ADP/ATP translocase 1 (ANT1), and its yeast homolog ADP/ATP carrier 2 (Aac2), cause the protein to accumulate along the protein import pathway, thereby obstructing general protein translocation into mitochondria. This impairs mitochondrial respiration, cytosolic proteostasis, and cell viability independent of ANT1’s nucleotide transport activity. The mutations act synergistically, as double mutant Aac2/ANT1 causes severe clogging primarily at the translocase of the outer membrane (TOM) complex. This confers extreme toxicity in yeast. In mice, expression of a super-clogger ANT1 variant led to neurodegeneration and an age-dependent dominant myopathy that phenocopy ANT1-induced human disease, suggesting clogging as a mechanism of disease. More broadly, this work implies the existence of uncharacterized amino acid requirements for mitochondrial carrier proteins to avoid clogging and subsequent disease.
-
- Biochemistry and Chemical Biology
Metabolic scaling, the inverse correlation of metabolic rates to body mass, has been appreciated for more than 80 years. Studies of metabolic scaling have largely been restricted to mathematical modeling of caloric intake and oxygen consumption, and mostly rely on computational modeling. The possibility that other metabolic processes scale with body size has not been comprehensively studied. To address this gap in knowledge, we employed a systems approach including transcriptomics, proteomics, and measurement of in vitro and in vivo metabolic fluxes. Gene expression in livers of five species spanning a 30,000-fold range in mass revealed differential expression according to body mass of genes related to cytosolic and mitochondrial metabolic processes, and to detoxication of oxidative damage. To determine whether flux through key metabolic pathways is ordered inversely to body size, we applied stable isotope tracer methodology to study multiple cellular compartments, tissues, and species. Comparing C57BL/6 J mice with Sprague-Dawley rats, we demonstrate that while ordering of metabolic fluxes is not observed in in vitro cell-autonomous settings, it is present in liver slices and in vivo. Together, these data reveal that metabolic scaling extends beyond oxygen consumption to other aspects of metabolism, and is regulated at the level of gene and protein expression, enzyme activity, and substrate supply.