- Research article
- Open Access
Cross-communication between Gi and Gs in a G-protein-coupled receptor heterotetramer guided by a receptor C-terminal domain
- Gemma Navarro1, 2, 3,
- Arnau Cordomí4,
- Marc Brugarolas1, 2, 3,
- Estefanía Moreno1, 2, 3,
- David Aguinaga1, 2, 3,
- Laura Pérez-Benito4,
- Sergi Ferre5,
- Antoni Cortés1, 2, 3,
- Vicent Casadó1, 2, 3,
- Josefa Mallol1, 2, 3,
- Enric I. Canela1, 2, 3,
- Carme Lluís1, 2, 3,
- Leonardo Pardo4Email author,
- Peter J. McCormick1, 2, 3, 6Email author and
- Rafael Franco1, 2, 3Email authorView ORCID ID profile
© Franco et al. 2018
- Received: 19 October 2017
- Accepted: 22 January 2018
- Published: 28 February 2018
G-protein-coupled receptor (GPCR) heteromeric complexes have distinct properties from homomeric GPCRs, giving rise to new receptor functionalities. Adenosine receptors (A1R or A2AR) can form A1R-A2AR heteromers (A1-A2AHet), and their activation leads to canonical G-protein-dependent (adenylate cyclase mediated) and -independent (β-arrestin mediated) signaling. Adenosine has different affinities for A1R and A2AR, allowing the heteromeric receptor to detect its concentration by integrating the downstream Gi- and Gs-dependent signals. cAMP accumulation and β-arrestin recruitment assays have shown that, within the complex, activation of A2AR impedes signaling via A1R.
We examined the mechanism by which A1-A2AHet integrates Gi- and Gs-dependent signals. A1R blockade by A2AR in the A1-A2AHet is not observed in the absence of A2AR activation by agonists, in the absence of the C-terminal domain of A2AR, or in the presence of synthetic peptides that disrupt the heteromer interface of A1-A2AHet, indicating that signaling mediated by A1R and A2AR is controlled by both Gi and Gs proteins.
We identified a new mechanism of signal transduction that implies a cross-communication between Gi and Gs proteins guided by the C-terminal tail of the A2AR. This mechanism provides the molecular basis for the operation of the A1-A2AHet as an adenosine concentration-sensing device that modulates the signals originating at both A1R and A2AR.
- C-terminal domain
- Molecular modeling
Adenosine is a purine nucleoside whose relevance in the central nervous system is mainly due to its role in regulating neurotransmitter release . The effects of adenosine are mediated by specific G-protein-coupled receptors (GPCRs) that are coupled to either Gs or Gi heterotrimeric Gαβγ proteins. The endogenous adenosine acts on four receptor subtypes – A1R, A2AR, A2BR, and A3R. Convergent and compelling evidence shows that GPCRs may form complexes constituted by a number of equal (homo) or different (hetero) receptor protomers . As agreed in the field, a GPCR heteromer displays characteristics that are different from those of the constituting protomers, thus giving rise to novel functional entities . Adenosine receptors have been used as a paradigm in the study of receptor homo- and heteromerization. For instance, A1R, which is Gi coupled, and A2AR, which is Gs coupled, form a functional heteromer .
The A1R-A2AR heteromer (A1-A2AHet) is found presynaptically in, inter alia, cortical glutamatergic terminals innervating the striatum and functions as a switch that differentially senses high and low concentrations of adenosine in the inter-synaptic space. Since adenosine has higher affinity for A1R than for A2AR, low concentrations predominantly activate A1R, engaging a Gi-mediated signaling, whereas higher adenosine concentrations also activate A2AR, engaging a Gs-mediated signaling . The physiological role of such a concentration-sensing device is remarkable as it allows adenosine to fine-tune modulate the release of neurotransmitters from presynaptic terminals. However, the mechanism by which A1-A2AHet integrates both Gi- and Gs-dependent signals is not yet understood. We have recently shown, using a combination of single-particle tracking experiments, bioluminescence resonance energy transfer (BRET) assays, and computer modeling, that the (minimal) functional A1-A2AHet/G protein unit is composed by a compact rhombus-shaped heterotetramer (with A1R and A2AR homodimers) bound to two different interacting heterotrimeric G proteins (Gs and Gi) . In the present study, we aim to understand the molecular intricacies underlying the signaling mediated by A1-A2AHet, in which (1) both receptors constituting the heteromer are activated by the same endogenous agonist and (2) is coupled to two different G proteins with opposite effects, i.e., one mediating the inhibition of the adenylate cyclase (Gi) and another mediating the activation of the enzyme (Gs). Our data identifies a new mechanism of signal transduction and provides the molecular basis to understand the unique properties of this heteromer, in which the C-terminal tail of the A2AR influences the Gi-mediated signaling of the partner A1R receptor.
Homodimerization of A1R and A2AR occurs through the transmembrane (TM) 4/5 interface and heterodimerization via the TM5/6 interface in the A1-A2AHet
The complex formed by Gs, Gi, and the A1-A2AHet as a signal transduction unit
In order to test the ability of Gs and Gi proteins to interact with the A1-A2AHet, we used BRET assays . Cells were transfected with cDNAs of A1R-nYFP and A2AR-cYFP, which only upon complementation can act as a BRET acceptor (YFP), and Renilla luciferase (Rluc) as a BRET donor fused to either Gi (Gi-Rluc) or Gs (Gs-Rluc). We observed significant energy transfer (Additional file 1: Figure S1C), indicating that Gi and Gs are bound to their respective receptors in the A1-A2AHet.
To further test for a cross-communication between G proteins in the Gs-Gi-heterotetramer signaling unit, we resolved the real-time signaling signature by using a label-free method, based on optical detection of dynamic changes in cellular density following receptor activation . The magnitude of the signaling by CPA or by CGS 21680 significantly decreased when cells co-expressing both receptors were pre-treated with either PTX or CTX (Fig. 3d). This phenomenon was not observed in cells expressing only A1R (Additional file 1: Figure S2G) or A2AR (Additional file 1: Figure S2H). Again, these results indicate the simultaneous coupling of interacting Gs and Gi proteins within the A1-A2AHet.
Simultaneous activation of both A1R and A2AR with CPA and CGS21680 increased cAMP to similar levels to those obtained with CGS21680 alone and the signal of co-activated receptors was inhibited by both PTX and CTX (Fig. 3c). Therefore, A1R agonist was able to decrease forskolin-induced cAMP (Fig. 3a, b) and yet was unable to decrease A2AR-mediated increases of cAMP (Fig. 3c). Consequently, when both receptors are co-activated in the heterotetramer, only the A2AR-mediated, but not the A1R-mediated signaling occurs. This finding was confirmed in label-free experiments, showing that receptor co-activation with CPA and CGS 21680 did not increase the time-response curve with respect to the activation with CGS 21680 alone (Fig. 3d green and yellow lines, respectively).
It has been shown that the mechanism for receptor-catalyzed nucleotide exchange in G proteins involves a large-scale opening of the α-helical domain (αAH) of the α-subunit, from the Ras domain, allowing GDP to freely dissociate [10–13]. Notably, our proposed model of the A1-A2AHet positions the αiAH and αsAH domains facing each other (Fig. 3e). The fact that both Gs- and Gi-specific toxins and Gs- and Gi-specific minigenes affect both Gs- and Gi-mediated coupling in the A1-A2AHet suggests that the proposed large-scale conformational changes of αAH domains is mutually dependent. We used molecular dynamics (MD) simulations of the A1-A2AHet in complex with Gs and Gi to evaluate intermolecular distances between the αsAH and αiAH domains when αiAH is in the closed conformation and αsAH is either in the open (Fig. 3e) or in the closed conformation (Fig. 3ef). In a previous report, double electron–electron resonance (DEER) distance distributions between spin labels attached to Arg90 (αiAH domain) and Glu238 (Ras domain) of Gi (the distance between Cα atoms is termed d[Arg90αi-Glu238αi] in the manuscript) or Asn112 (αsAH) and Asn261 (Ras) of Gs (d[Asn112αs-Asn261αs]) permitted to faithfully monitor the equilibrium within the open (distance of ~40 Å) and closed (~20 Å) conformation of the αAH domain . Here, we measured the intermolecular distance between the αsAH and αiAH domains using Cα atoms of Arg90 of αi and Asn112 of αs (d[Arg90αi-Asn112αs]). This d[Arg90αi-Asn112αs] intermolecular distance between αiAH in the closed conformation (d[Arg90αi-Glu238αi]: 11 Å, yellow line in Fig. 3e) and αsAH in the closed conformation (d[Asn112αs-Asn261αs]: 14 Å, green line in Fig. 3e) has an average value of 108 Å for inactive A1-A2AHet (Fig. 3e, dark red line). Activation of A2AR would trigger the opening of αsAH (d[Asn112αs-Asn261αs]: 52 Å; Fig. 3f, green line), necessary for GDP/GTP exchange, decreasing the d[Arg90αi-Asn112αs] distance between αiAH and αsAH to 60 Å (Fig. 3f, dark red line). Although the results are based on a single trajectory, it is unlikely that additional replicates would change, in a significant manner, the distances reported from the simulations. Moreover, the differences between the distances are so substantial that results from more simulations would not have a significant impact. We hypothesize that a similar change occurs with activation of A1R. This indicates that both receptors can signal via their cognate G protein by opening their αAH domain. However, in the compact rhombus-shaped A1-A2AHet model, simultaneous opening of both αAH domains (co-activation with CPA and CGS 21680) would not be possible due to a steric clash in such open conformations (Fig. 3g). Due to this steric clash, MD simulations of this open αiAH-open αsAH conformation in the absence of interference peptides (see below) cannot be performed.
Altering the heteromer interface of A1-A2AHet enables simultaneous Gi and Gs signaling
A1-A2AHet as an adenosine concentration-sensing device
Recruitment of β-arrestin-2 by the A1-A2AHet
Using the TAT-fused synthetic peptides we investigated whether the quaternary structure of the A1-A2AHet determines its putative selective A2AR-dependent β-arrestin-2 recruitment. As a negative control, we first corroborated that TM4, TM5, and TM6 peptides of A2AR do not interfere with A1R-mediated signaling (Additional file 1: Figure S3C). Pretreatment of cells expressing Arr-Rluc, A2AR-YFP and non-fused A1R (Fig. 6a), or Arr-Rluc, A1R-YFP and non-fused A2AR (Fig. 6b) with TM4, TM5, and TM6 peptides, but not in the absence of peptides (control) or with the TM7 peptide (negative control), allowed the detection of positive BRET (recruitment of β-arrestin-2) not only when cells were treated with the A2AR-selective agonist CGS-21680 (white bars), but also when treated with the A1R-selective agonist CPA (black bars) (Figs. 6a, b). Importantly, when cells expressing Arr-Rluc, A2AR-YFP, and non-fused A1R were co-activated by CPA and CGS-21680 (striped bars), BRET measurement in the presence of TM4, TM5, or TM6 peptides, but neither in the absence of peptides nor in the presence of TM7 peptide, significantly increased relative to the values obtained by the action of a single agonist (Fig. 6a). The trend is similar in cells expressing Arr-Rluc, A1R-YFP, and non-fused A2AR, but not statistically significant (Fig. 6b). These results indicate that alteration of the A1R-A2AR heteromer interface within the A1-A2AHet allows simultaneous recruitment of β-arrestin-2 to A1R and A2AR when both receptors are activated. Interference peptides abolish cross-communication of G proteins, permitting CPA to activate Gi (Gβγ moving away from Gαi) and recruitment of β-arrestin-2 to A1R, as well as Gs activation by CGS-21680 (Gβγ moving away from Gαs) and simultaneous recruitment of β-arrestin-2 to A2AR.
The C-terminal domain of A2AR is responsible for the dominant A2AR-mediated signaling
We measured cAMP production in cells expressing A1R and wild-type or truncated A2AR receptors (Fig. 7c). Truncated A2AR were able to signal as wild-type receptors. Interestingly, the dominant Gs-mediated signaling when A1R and A2AR were co-activated decreased progressively with the shortening of the A2AR C-tail (Fig. 7c, striped bars). In fact, CPA inhibited CGS-21680-induced cAMP accumulation when truncated receptors were expressed, showing that, in these heteromers, A1R were functional (Additional file 1: Figure S5). Figure 7e shows a detailed view of the orientation of the C-tail (102 amino acids, Gln311-Ser412) of both A2AR protomers in the A1-A2AHet, which was modeled as suggested for the OXER , together with the structure of β-arrestin-2 in complex with V2 vasopressin receptor . It is important to note that the exact conformation of the A2AR C-tail cannot unambiguously be determined, thus, we only predict its orientation as explained in detail in Additional file 1: Figure S4. The fact that the C-tail of the αs-unbound A2AR protomer points toward the αsAH domain suggests that this C-tail is influencing the conformational changes required to open the αsAH, and thus controlling the balance between Gs and Gi activation. Next, we measured β-arrestin-2 recruitment by BRET assays in cells expressing A1R and wild-type or truncated A2AR receptors. In cells expressing non-fused A1R, Arr-Rluc and A2AR-YFP, A2AΔ40R-YFP, or A2AΔCTR-YFP, the A1R agonist CPA could increase BRET values only when the heteromer is formed with A2AR-truncated receptors. In these conditions, co-activation with CPA and CGS-21860 induced a BRET increase higher than the one obtained with CGS-21680 alone (Fig. 7d). These results indicate that the selective A2AR-dependent β-arrestin-2 recruitment in the A1-A2AHet decreases progressively with the shortening of the A2AR C-tail (Fig. 7d).
As previously reviewed [2, 3, 22], the intercommunication between protomers of a GPCR heteromer can be observed at the level of agonist binding, ligand-induced cross-conformational changes between receptor protomers, and the binding of GPCR-associated proteins, including heterotrimeric G proteins and β-arrestins. The intercommunication between protomers is a consequence of a defined quaternary structure that is responsible for the specific functional characteristics of the heteromer. For GPCR heteromers, such as A1-A2AHet, constituted by receptors sensing the same hormone but producing opposite signaling effects, it is not obvious how a defined quaternary structure achieves this dual behavior. A1-A2AHet acts as a concentration-sensing device that allows adenosine to signal by one or the other coupled G protein (Gs or Gi) to fine-tune modulate the release of neurotransmitters from presynaptic terminals. In the present study, we solved this question by discovering a new mechanism of signal transduction, a cross-communication between Gi and Gs in the A1-A2AHet guided by the A2AR C-terminal domain.
We have shown that cross-communication between Gi and Gs proteins involves the formation of a GPCR heterotetramer (i.e., one homodimer of A1R and one of A2AR) that has a 2:2:1:1 (A2AR:A1R:Gs:Gi) stoichiometry. From our data, it is deduced that the cross-talk between Gi and Gs resides on the structural constraints surrounding the mechanism for GDP/GTP exchange, which involves the opening of the αAH domain of the α-subunit of any given G protein. We propose that cross-communication in the Gs-Gi-heterotetramer signaling unit is a property associated with a specific quaternary structure, the compact rhombus-shaped A1-A2AHet (the TM4/5 interface for homodimerization and the TM5/6 interface for heterodimerization), which positions the αiAH and αsAH domains in close proximity, making their conformational changes mutually dependent in a way that simultaneous opening of both αAH domains would not be possible due to a steric clash in such open conformations. Alterations of this quaternary structure of the A1-A2AHet by insertion of synthetic peptides between A1R and A2AR blocks this cross-communication without disrupting the heteromer and permits simultaneous activation of Gi and Gs in the heteromer. Since the cross-talk between Gi and Gs resides on the structural constraints imposed by defined TM interfaces in the heteromer, it is important to note that other heterotetramers, mainly those sensing different hormones and with a different quaternary structure, might not display this cross-communication among G proteins. Moreover, although, from a structural point of view, the A1-A2AHet is capable to recruit not only two G proteins but also two β-arrestins, the cross-talk between Gi and Gs, in which Gs activation inhibits the simultaneous activation of Gi, blocks A1R agonist-promoted arrestin recruitment. Alteration of the A1-A2AHet by insertion of synthetic peptides between A1R and A2AR facilitates simultaneous activation of Gi and Gs and the corresponding binding of two β-arrestins to A1R and A2AR. Our finding that Gi is dependent on Gs-mediated signaling strengthens the conclusion that cross-talk across G proteins is a potentially important functional property of GPCR heteromers. Remarkably, when both receptors are co-activated in this heterotetramer, only A2AR-mediated, but not A1R-mediated signaling occurs. We show that the ability of blunting A1R-mediated signaling when Gs is engaged is dependent of the long C-terminus of the A2AR. In the absence of A2AR activation by agonists, or in the absence of the C-terminal domain of A2AR, the A1R-mediated signaling via Gi is totally functional. The most straightforward hypothesis is that the opening of αsAH parallels a movement of the C-tail to block the opening of αiAH.
Adenosinergic signaling in mammalians is important for energy and temperature homeostasis and for neuroregulation. Multiplicity of adenosine actions is due to a balance between the expression of specific receptors and producing/degrading enzymes and to the biological diversity due to a membrane network established by the interaction among purinergic receptors . Ciruela et al.  first identified the occurrence of heteromers formed by A1R-Gi- and A2A-Gs-coupled adenosine receptors that participate in the regulation of glutamate release by neurons projecting from the cortex to the striatum. The same A1-A2AHet can be found in astrocytes modulating the transport of γ-amino butyric acid (GABA) . Differently from the modulation of neuronal glutamate release, the A1R-Gi-coupled receptor activates and the A2AR-Gs-coupled receptor inhibits the modulation of GABA transport. Under conditions of high extracellular adenosine concentrations, such as hypoxic conditions , the nucleoside will bind to both the high (A1R) and the low (A2AR) affinity receptors in the heteromer, and the predominant A2AR-mediated signaling via Gs will result in counteraction of astrocytic GABA transport. Our results show that the asymmetric signaling is possible because the long C-terminus of A2AR blunts Gi-mediated signaling. We have therefore elucidated the mechanism by which the A1-A2AHet functions as an adenosine concentration-sensing device that can promote even opposite signaling responses depending on the extracellular concentration of adenosine. The molecular mechanism involves the C-terminal domain of the activated Gs-coupled A2AR, which hinders the activation of A1R coupled to Gi.
Using a convergent approach including biochemical, biophysical, cell biology, and molecular biology techniques, together with in silico molecular models, we here provide the mode of action of a membrane receptor complex that responds depending on the concentration of adenosine, a hormone and a neuroregulatory molecule. The concentration sensor is a heteromer composed of four adenosine receptors (two A1 and two A2A) and two G proteins (Gi and Gs). Despite Gi sits underneath the A1 receptor dimer and Gs sits underneath the A2A receptor dimer, both G proteins do interact and are able to convey allosteric regulation depending on how the functional unit is activated. At low adenosine concentrations Gi is engaged via A1 activation without affecting/engaging Gs signaling. At higher concentrations Gs is engaged via A2A activation, and this engagement blocks Gi-mediated signaling. The reason why a rhombus-shaped apparently symmetric structure results in asymmetric signaling is due to the long C-terminal tail of the A2A receptor. In fact, both deletion of the C-terminal end or treatment with interfering peptides derived from the sequence of TM segments of the receptors impair allosteric cross-interaction between receptors and G proteins within the macromolecule, and the device loses its concentration sensing properties.
Cell culture and transient transfection
HEK-293 T cells were grown at 37 °C in in Dulbecco’s modified Eagle’s medium (DMEM) (Gibco) supplemented with 2 mM L-glutamine, 100 U/mL penicillin/streptomycin, and 5% (v/v) heat inactivated fetal bovine serum (all supplements were from Invitrogen, Paisley, Scotland, UK). Cells were transiently transfected with cDNA corresponding to receptors, fusion proteins, A2AR mutant constructs, or minigene vectors using polyethylenimine (Sigma-Aldrich, Cerdanyola del Vallés, Spain) as described elsewhere .
Expression vectors, A2AR mutants and minigenes
Sequences encoding amino acid residues 1–155 or 155–238 of YFP-Venus protein, were subcloned in pcDNA3.1 to obtain the YFP Venus hemi-truncated proteins (nYFP and cYFP). The human cDNAs for A2AR, mutant A2AR, A1R, and Gi or Gs proteins cloned into pcDNA3.1, were amplified without their stop codons using sense and antisense primers harboring unique EcoRI and BamHI sites to subclone receptors in pcDNA3.1RLuc vector (pRLuc-N1 PerkinElmer, Wellesley, MA, USA) and EcoRI and KpnI to subclone receptors in pEYFP-N1 (enhanced yellow variant of GFP; Clontech, Heidelberg, Germany), pcDNA3.1-nVenus, or pcDNA3.1-cVenus vectors. The amplified fragments were subcloned to be in-frame with restriction sites of the corresponding vectors to give the plasmids that express receptors fused to RLuc, YFP, nYFP or cYFP on the C-terminal end (A1R-Rluc, A2AR-Rluc, Gi-RLuc, Gs-RLuc, A1R-YFP, A2AR-YFP, A2AΔ40R-YFP, A2AΔCTR-YFP, A1R-nYFP, A2A-nYFP, and A2A-cYFP). Expression of constructs was tested by confocal microscopy and the receptor-fusion protein functionality by second messengers, ERK1/2 phosphorylation and cAMP production as described previously [4, 26–28]. Mutants with a deletion of aa 372 to aa 412 (A2AΔ40R) or aa 321 to aa 412 (A2AΔCTR) on the C-terminal domain of A2AR were generated as previously described . “Minigene” plasmid vectors are constructs designed to express relatively short polypeptide sequences following their transfection into mammalian cells. Here, we used minigene constructs encoding 11 amino acid residues from the C-terminus sequence of α subunit of Gi1/2 or Gs. The peptide coded by every minigene inhibits the coupling of the G (Gi1/2 or Gs) protein to the receptor and, consequently, it inhibits the G-protein-mediated cellular response, as previously described . The cDNA encoding the last 11 amino acids of human Gα subunit corresponding to Gi1/2 (IKNNLKDCGLF) or Gs (QRMHLRQYELL), inserted in a pcDNA 3.1 plasmid vector, was generously provided by Dr. Heidi Hamm.
MEYMVYFNFFVWVLPPLLLMVLIYLYGRKKRRQRRR for TM5 of A1R, RRRQRRKKRGYLALILFLFALSWLPLHILNCITLF for TM6 of A1R, ILTYIAIFLTHGNSAMNPIVYAFRIYGRKKRRQRRR for TM7 of A1R, VYITVELAIAVLAILGNVLVCWAVWYGRKKRRQRRR for TM1 of A2AR, YGRKKRRQRRRYFVVSLAAADIAVGVLAIPFAITI for TM2 of A2AR, LFIACFVLVLTQSSIFSLLAIAIYGRKKRRQRRR for TM3 of A2AR, YGRKKRRQRRRAKGIIAICWVLSFAIGLTPMLGW for TM4 of A2AR, MNYMVYFNFFACVLVPLLLMLGVYLYGRKKRRQRRR for TM5 of A2AR, YGRKKRRQRRRLAIIVGLFALCWLPLHIINCFTFF for TM6 of A2AR, LWLMYLAIVLSHTNSVVNPFIYAYYGRKKRRQRRR for TM7 of A2AR.
YGRKKRRQRRRILGIWAVSLAIMVPQAAVME for TM4 of OX1R, SSFFIVTYLAPLGLMAMAYFQIFYGRKKRRQRRR for TM5 of OX1R, YASFTFSHWLVYANSAANPIIYNFYGRKKRRQRRR for TM7 of OX1R
Bimolecular fluorescence complementation assay (BiFC)
HEK-293 T cells were transiently transfected with equal amounts of the cDNA for fusion proteins of the hemi-truncated Venus (1 μg of each cDNA). At 48 h after transfection, cells were treated for 4 h at 37° with medium or TAT peptides (4 μM) before plating 20 μg of protein in 96-well black microplates (Porvair, King’s lynn, UK). To quantify reconstituted YFP Venus expression, fluorescence at 530 nm was read in a Fluoro Star Optima Fluorimeter (BMG Labtechnologies, Offenburg, Germany) equipped with a high-energy xenon flash lamp, using a 10 nm bandwidth excitation filter at 400 nm reading. Protein fluorescence expression was determined as fluorescence of the sample minus the fluorescence of cells not expressing the fusion proteins (basal). Cells expressing receptor-cVenus and nVenus or receptor-nVenus and cVenus showed similar fluorescence levels than untransfected cells.
Bioluminescence resonance energy transfer (BRET)
HEK-293 T cells were transiently transfected with a constant amount of cDNA for Rluc fusion proteins and increasing amounts of cDNA for YFP fusion proteins. At 48 h after transfection, 20 μg of cell suspension were plated in 96-well black microplates for fluorescence detection or in 96-well white microplates for BRET readings and Rluc quantification. YFP fluorescence at 530 nm was quantified in a Fluoro Star Optima Fluorimeter as described above. BRET signal was collected 1 min after addition of 5 μM coelenterazine H (Molecular Probes, Eugene, OR, USA) using a Mithras LB 940. The integration of the signals detected in the short-wavelength filter at 485 nm and the long-wavelength filter at 530 nm was recorded. To quantify protein-RLuc expression, luminescence readings were also performed after 10 minutes of adding 5 μM coelenterazine H. The net BRET is defined as (long-wavelength emission/short-wavelength emission)–Cf, where Cf corresponds to long-wavelength emission/short-wavelength emission for the donor construct expressed alone in the same experiment. BRET is expressed as milli-BRET units (net BRET × 1000). To calculate maximum BRET (BRETmax) from saturation curves, data were fitted to a nonlinear regression equation, assuming a single-phase saturation curve with GraphPad Prism software (San Diego, CA, USA).
Proximity ligation assay (PLA)
HEK293T cells were grown on glass coverslips and fixed in 4% paraformaldehyde for 15 min, washed with phosphate-buffered saline containing 20 mM glycine, permeabilized with the same buffer containing 0.05% Triton X-100, and successively washed with tris-buffered saline. Heteromers were detected using the Duolink II in situ PLA detection Kit (OLink; Bioscience, Uppsala, Sweden) following supplier’s instructions. A mixture of the primary antibodies (mouse anti-A2AR antibody (1:100; 05-717, Millipore, Darmstadt, Germany; RRID:AB_309931) and rabbit anti-A1R antibody (1:100; ab82477, Abcam, Bristol, UK; RRID: AB_2049141)) was used to detect A1-A2AHet together with PLA probes detecting mouse or rabbit antibodies. Then, samples were processed for ligation and amplification with a Detection Reagent Red and were mounted using a DAPI-containing mounting medium. Samples were analyzed in a Leica SP2 confocal microscope (Leica Microsystems, Mannheim, Germany) equipped with 405 nm and 561 nm laser lines. For each field of view, a stack of two channels (one per staining) and 4–6 Z-stacks with a step size of 1 μm were acquired. Images were opened and processed with Image J software (National Institutes of Health, Bethesda, MD, USA).
cAMP determination assays
HEK-293 T cells expressing adenosine receptors were incubated for 4 h in serum-free medium containing 50 μM zardeverine. Cells were plated in 384-well white microplates (1500 cells/well), pre-treated with toxins or the corresponding vehicle for the indicated time, stimulated with agonists for 15 min before adding medium or 0.5 μM forskolin, and incubated for an additional 15 min. cAMP production was quantified by a TR-FRET (Time-Resolved Fluorescence Resonance Energy Transfer) methodology using the LANCE Ultra cAMP kit (PerkinElmer) and fluorescence at 665 nm was analyzed on a Pherastar Flagship Microplate Reader (BMG Labtech, Ortenberg, Germany).
Dynamic mass redistribution (DMR) assays
The heteromer-induced cell signaling signature was determined using an EnSpire® Multimode Plate Reader (PerkinElmer, Waltham, MA, USA) by a label-free technology. Refractive waveguide grating optical biosensors, integrated in 384-well microplates, allow extremely sensitive measurements of changes in local optical density in a detecting zone up to 150 nm above the surface of the sensor. Cellular mass movements induced upon receptor activation were detected by illuminating the underside of the biosensor with polychromatic light and measured as changes in wavelength of the reflected monochromatic light, which is a sensitive function of the index of refraction. The magnitude of this wavelength shift (in picometers) is directly proportional to the amount of DMR. Briefly, 24 h before the assay, cells were seeded at a density of 7500 cells per well in 384-well sensor microplates with 40 μL growth medium and cultured for 24 h (37 °C, 5% CO2) to obtain 70–80% confluent monolayers. Previous to the assay, cells were pre-treated with medium or toxins as indicated and incubated for 2 h in 40 μL per well of assay-buffer (HBSS with 20 mM HEPES, pH 7.15) in the reader at 24 °C. Thereafter, the sensor plate was scanned and a baseline optical signature was recorded prior to addition of 10 μL of receptor agonist dissolved in assay buffer containing 0.1% DMSO. DMR responses were monitored for at least 8000 s and data were analyzed using EnSpire Workstation Software v. 4.10.
The structural model of the A1-A2AHet bound to Gs (closed αsAH domain) and Gi (closed αiAH domain) was taken from our previous work . This previous structural model contains a A2AR-based homology model of A1R. The structure of the adenosine A1R has recently been revealed , showing a remarkably similar structure (Additional file 1: Figure S6A). This structure of A1R contains a TM4/5 dimer interface that is in close agreement with our model (Additional file 1: Figure S6B). An intermediate conformation (obtained using the g_morph tool of the GROMACS package ) between the closed αAH domain (PDB id 1AZT) and the conformation observed in the crystal structure of the β2-AR in complex with Gs (PDB id 3SN6) was used to model the open αAH domain (Additional file 1: Figure S6C). This conformation is supported by DEER spectroscopy, deuterium-exchange and electron microscopy data [11–13]. The active state of β-arrestin-2 was built using a multi-template alignment combining the structure of the active β-arrestin-1 (PDB id 4JQI)  and the structure of rhodopsin in complex with visual β-arrestin (PDB id 4ZWJ) . Structural models of the A1-A2AHet bound to β-arrestin-2 were modeled using the crystal structure of rhodopsin bound to β-arrestin (PDB id 4ZWJ) . The structure of TM6 of A2AR fused to the cell-penetrating TAT peptide was modeled from the structure of A2AR. Molecular models of the A1-A2AHet with the TAT-fused TM6 peptide, disrupting the heteromer interface between A1R and A2AR, in complex with Gs (open αsAH domain) and Gi (open αiAH domain), was built from the structure of A1-A2AHet. The conformation of the proximal C-tail of A2AR (Ser305-Ala317) was modeled based on squid rhodopsin . The remaining part of the C-tail (1Gly319–Ser412), cannot be unambiguously determined and it was modeled as suggested for the oxoeicosanoid receptor (OXER) , together with the structure derived from the human V2 vasopressin receptor in complex with β-arrestin-2  (see Additional file 1: Figure S4 for details). Additional file 2: Table S1 shows the template structures used in the protein models. Modeller 9.12 was used to build these models . The molecular models of A1-A2AHet in complex with Gs and Gi or β-arrestin, in the absence or presence of the TAT-fused TM6 peptide, were embedded in a pre-equilibrated box containing a lipid bilayer (~800 POPC molecules) with explicit solvent (~110,000 waters) and 0.15 M concentration of Na+ and Cl– (~1800 ions). These initial complexes were energy-minimized and subsequently subjected to a 21 ns MD equilibration, with positional restraints on protein coordinates. These restraints were released and 500 ns of MD trajectory were produced at constant pressure and temperature. Computer simulations were performed with the GROMACS 4.6.3 simulation package , using the AMBER99SB force field as implemented in GROMACS and Berger parameters for POPC lipids. This procedure has been previously validated .
We would like to thank Jasmina Jiménez for technical help (University of Barcelona). RF, PJM and LP participate in the European COST Action CM1207 (GLISTEN). Authors gratefully acknowledge the computer resources provided by the Barcelona Supercomputing Center - Centro Nacional de Supercomputación.
This study was supported by grants from the Spanish Ministerio de Economía y Competitividad (SAF2015-74627-JIN, BFU2015-64405-R and SAF2016-77830-R; they may contain FEDER funds) and by the intramural funds of the National Institute on Drug Abuse to SF. RF, PJM, and LP participate in the European COST Action CM1207 (GLISTEN).
GN performed the molecular biology experiments. GN, MB, EM, and DA performed BRET experiments. AC and LP-B performed molecular modeling studies. SF, AC, VC, JM, and EIC analyzed the data. CL, LP, PJM, and RF designed the experiments, supervised the work in the respective laboratories and wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Thompson SM, Haas HL, Gahwiler BH. Comparison of the actions of adenosine at pre- and postsynaptic receptors in the rat hippocampus in vitro. J Physiol. 1992;451:347–63.View ArticlePubMedPubMed CentralGoogle Scholar
- Franco R, Martinez-Pinilla E, Lanciego JL, Navarro G. Basic pharmacological and structural evidence for class A G-protein-coupled receptor heteromerization. Front Pharmacol. 2016;7:76.View ArticlePubMedPubMed CentralGoogle Scholar
- Ferre S, Baler R, Bouvier M, Caron MG, Devi LA, Durroux T, Fuxe K, George SR, Javitch JA, Lohse MJ, et al. Building a new conceptual framework for receptor heteromers. Nat Chem Biol. 2009;5:131–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Ciruela F, Casado V, Rodrigues RJ, Lujan R, Burgueno J, Canals M, Borycz J, Rebola N, Goldberg SR, Mallol J, et al. Presynaptic control of striatal glutamatergic neurotransmission by adenosine A1-A2A receptor heteromers. J Neurosci. 2006;26:2080–7.View ArticlePubMedGoogle Scholar
- Navarro G, Cordomi A, Zelman-Femiak M, Brugarolas M, Moreno E, Aguinaga D, Perez-Benito L, Cortes A, Casado V, Mallol J, et al. Quaternary structure of a G-protein-coupled receptor heterotetramer in complex with Gi and Gs. BMC Biol. 2016;14:26.View ArticlePubMedPubMed CentralGoogle Scholar
- Schwarze SR, Ho A, Vocero-Akbani A, Dowdy SF. In vivo protein transduction: delivery of a biologically active protein into the mouse. Science. 1999;285:1569–72.View ArticlePubMedGoogle Scholar
- Carriba P, Navarro G, Ciruela F, Ferre S, Casado V, Agnati L, Cortes A, Mallol J, Fuxe K, Canela EI, et al. Detection of heteromerization of more than two proteins by sequential BRET-FRET. Nat Methods. 2008;5:727–33.View ArticlePubMedGoogle Scholar
- Gilchrist A, Li A, Hamm HE. G alpha COOH-terminal minigene vectors dissect heterotrimeric G protein signaling. Sci STKE. 2002;2002:pl1.PubMedGoogle Scholar
- Schroder R, Schmidt J, Blattermann S, Peters L, Janssen N, Grundmann M, Seemann W, Kaufel D, Merten N, Drewke C, et al. Applying label-free dynamic mass redistribution technology to frame signaling of G protein-coupled receptors noninvasively in living cells. Nat Protoc. 2011;6:1748–60.View ArticlePubMedGoogle Scholar
- Van Eps N, Preininger AM, Alexander N, Kaya AI, Meier S, Meiler J, Hamm HE, Hubbell WL. Interaction of a G protein with an activated receptor opens the interdomain interface in the alpha subunit. Proc Natl Acad Sci U S A. 2011;108:9420–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Chung KY, Rasmussen SG, Liu T, Li S, DeVree BT, Chae PS, Calinski D, Kobilka BK, Woods VL Jr, Sunahara RK. Conformational changes in the G protein Gs induced by the beta2 adrenergic receptor. Nature. 2011;477:611–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Westfield GH, Rasmussen SG, Su M, Dutta S, DeVree BT, Chung KY, Calinski D, Velez-Ruiz G, Oleskie AN, Pardon E, et al. Structural flexibility of the G alpha s alpha-helical domain in the beta2-adrenoceptor Gs complex. Proc Natl Acad Sci U S A. 2011;108:16086–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Dror RO, Mildorf TJ, Hilger D, Manglik A, Borhani DW, Arlow DH, Philippsen A, Villanueva N, Yang Z, Lerch MT, et al. SIGNAL TRANSDUCTION. Structural basis for nucleotide exchange in heterotrimeric G proteins. Science. 2015;348:1361–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Kang Y, Zhou XE, Gao X, He Y, Liu W, Ishchenko A, Barty A, White TA, Yefanov O, Han GW, et al. Crystal structure of rhodopsin bound to arrestin by femtosecond X-ray laser. Nature. 2015;523:561–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Piersen CE, True CD, Wells JN. A carboxyl-terminally truncated mutant and nonglycosylated A2a adenosine receptors retain ligand binding. Mol Pharmacol. 1994;45:861–70.PubMedGoogle Scholar
- Canals M, Burgueno J, Marcellino D, Cabello N, Canela EI, Mallol J, Agnati L, Ferre S, Bouvier M, Fuxe K, et al. Homodimerization of adenosine A2A receptors: qualitative and quantitative assessment by fluorescence and bioluminescence energy transfer. J Neurochem. 2004;88:726–34.View ArticlePubMedGoogle Scholar
- Palmer TM, Stiles GL. Identification of an A2a adenosine receptor domain specifically responsible for mediating short-term desensitization. Biochemistry. 1997;36:832–8.View ArticlePubMedGoogle Scholar
- Klinger M, Kuhn M, Just H, Stefan E, Palmer T, Freissmuth M, Nanoff C. Removal of the carboxy terminus of the A2A-adenosine receptor blunts constitutive activity: differential effect on cAMP accumulation and MAP kinase stimulation. Naunyn Schmiedebergs Arch Pharmacol. 2002;366:287–98.View ArticlePubMedGoogle Scholar
- Schroder R, Merten N, Mathiesen JM, Martini L, Kruljac-Letunic A, Krop F, Blaukat A, Fang Y, Tran E, Ulven T, et al. The C-terminal tail of CRTH2 is a key molecular determinant that constrains Galphai and downstream signaling cascade activation. J Biol Chem. 2009;284:1324–36.View ArticlePubMedGoogle Scholar
- Blattermann S, Peters L, Ottersbach PA, Bock A, Konya V, Weaver CD, Gonzalez A, Schroder R, Tyagi R, Luschnig P, et al. A biased ligand for OXE-R uncouples Galpha and Gbetagamma signaling within a heterotrimer. Nat Chem Biol. 2012;8:631–8.View ArticlePubMedGoogle Scholar
- Shukla AK, Manglik A, Kruse AC, Xiao K, Reis RI, Tseng WC, Staus DP, Hilger D, Uysal S, Huang LY, et al. Structure of active beta-arrestin-1 bound to a G-protein-coupled receptor phosphopeptide. Nature. 2013;497:137–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Maurice P, Kamal M, Jockers R. Asymmetry of GPCR oligomers supports their functional relevance. Trends Pharmacol Sci. 2011;32:514–20.View ArticlePubMedGoogle Scholar
- Schicker K, Hussl S, Chandaka GK, Kosenburger K, Yang JW, Waldhoer M, Sitte HH, Boehm S. A membrane network of receptors and enzymes for adenine nucleotides and nucleosides. Biochim Biophys Acta. 2009;1793:325–34.View ArticlePubMedGoogle Scholar
- Cristovao-Ferreira S, Navarro G, Brugarolas M, Perez-Capote K, Vaz SH, Fattorini G, Conti F, Lluis C, Ribeiro JA, McCormick PJ, et al. A1R-A2AR heteromers coupled to Gs and G i/0 proteins modulate GABA transport into astrocytes. Purinergic Signal. 2013;9:433–49.View ArticlePubMedPubMed CentralGoogle Scholar
- Lopes LV, Sebastiao AM, Ribeiro JA. Adenosine and related drugs in brain diseases: present and future in clinical trials. Curr Top Med Chem. 2011;11:1087–101.View ArticlePubMedGoogle Scholar
- Canals M, Marcellino D, Fanelli F, Ciruela F, de Benedetti P, Goldberg SR, Neve K, Fuxe K, Agnati LF, Woods AS, et al. Adenosine A2A-dopamine D2 receptor-receptor heteromerization: qualitative and quantitative assessment by fluorescence and bioluminescence energy transfer. J Biol Chem. 2003;278:46741–9.View ArticlePubMedGoogle Scholar
- Gonzalez S, Moreno-Delgado D, Moreno E, Perez-Capote K, Franco R, Mallol J, Cortes A, Casado V, Lluis C, Ortiz J, et al. Circadian-related heteromerization of adrenergic and dopamine D(4) receptors modulates melatonin synthesis and release in the pineal gland. PLoS Biol. 2012;10:e1001347.View ArticlePubMedPubMed CentralGoogle Scholar
- Navarro G, Ferre S, Cordomi A, Moreno E, Mallol J, Casado V, Cortes A, Hoffmann H, Ortiz J, Canela EI, et al. Interactions between intracellular domains as key determinants of the quaternary structure and function of receptor heteromers. J Biol Chem. 2010;285:27346–59.View ArticlePubMedPubMed CentralGoogle Scholar
- Burgueno J, Blake DJ, Benson MA, Tinsley CL, Esapa CT, Canela EI, Penela P, Mallol J, Mayor F Jr, Lluis C, et al. The adenosine A2A receptor interacts with the actin-binding protein alpha-actinin. J Biol Chem. 2003;278:37545–52.View ArticlePubMedGoogle Scholar
- He SQ, Zhang ZN, Guan JS, Liu HR, Zhao B, Wang HB, Li Q, Yang H, Luo J, Li ZY, et al. Facilitation of mu-opioid receptor activity by preventing delta-opioid receptor-mediated codegradation. Neuron. 2011;69:120–31.View ArticlePubMedGoogle Scholar
- Glukhova A, Thal DM, Nguyen AT, Vecchio EA, Jorg M, Scammells PJ, May LT, Sexton PM, Christopoulos A. Structure of the adenosine A1 receptor reveals the basis for subtype selectivity. Cell. 2017;168:867–77. e813View ArticlePubMedGoogle Scholar
- Pronk S, Pall S, Schulz R, Larsson P, Bjelkmar P, Apostolov R, Shirts MR, Smith JC, Kasson PM, van der Spoel D, et al. GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics. 2013;29:845–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Murakami M, Kouyama T. Crystal structure of squid rhodopsin. Nature. 2008;453:363–7.View ArticlePubMedGoogle Scholar
- Marti-Renom MA, Stuart AC, Fiser A, Sanchez R, Melo F, Sali A. Comparative protein structure modeling of genes and genomes. Annu Rev Biophys Biomol Struct. 2000;29:291–325.View ArticlePubMedGoogle Scholar
- Cordomi A, Caltabiano G, Pardo L. Membrane protein simulations using AMBER force field and berger lipid parameters. J Chem Theory Comput. 2012;8:948–58.View ArticlePubMedGoogle Scholar