Pharmacological evidence for a metabotropic glutamate receptor heterodimer in neuronal cells

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
References
Article and author information
Metrics

Abstract

Metabotropic glutamate receptors (mGluRs) are mandatory dimers playing important roles in regulating CNS function. Although assumed to form exclusive homodimers, 16 possible heterodimeric mGluRs have been proposed but their existence in native cells remains elusive. Here, we set up two assays to specifically identify the pharmacological properties of rat mGlu heterodimers composed of mGlu2 and 4 subunits. We used either a heterodimer-specific conformational LRET-based biosensor or a system that guarantees the cell surface targeting of the heterodimer only. We identified mGlu2-4 specific pharmacological fingerprints that were also observed in a neuronal cell line and in lateral perforant path terminals naturally expressing mGlu2 and mGlu4. These results bring strong evidence for the existence of mGlu2-4 heterodimers in native cells. In addition to reporting a general approach to characterize heterodimeric mGluRs, our study opens new avenues to understanding the pathophysiological roles of mGlu heterodimers.

https://doi.org/10.7554/eLife.25233.001

Introduction

G-protein-coupled receptors (GPCRs) are essential in cell-cell communication and are considered as major drug targets. Although recognized as activating G proteins in a monomeric form (Whorton et al., 2007), numerous studies revealed their possible association into hetero-oligomers enabling allosteric controls between receptors (Pin et al., 2007; Ferré et al., 2014). The validation of this concept in vivo remains difficult and is a matter of intense debates (Pin et al., 2007; Bouvier and Hébert, 2014; Lambert and Javitch, 2014). The metabotropic glutamate (mGlu) receptors are members of the class C GPCRs activated by the main excitatory neurotransmitter, glutamate. These receptors are strict dimers and have until recently only been considered as homodimers (Romano et al., 1996; Kunishima et al., 2000). However, recent studies revealed the possible existence of heterodimeric mGluRs (Doumazane et al., 2011; Kammermeier, 2012; Yin et al., 2014; Niswender et al., 2016), as observed with other class C GPCRs (Marshall et al., 1999; Zhao et al., 2003; Pin and Bettler, 2016). The mGluRs constitute therefore an interesting model to tackle the issue of heterodimeric GPCRs in vivo.

Among the eight mGluRs, mGlu1 and 5 (group I) are mainly postsynaptic, while mGlu2 and 3 (group II) and mGlu4, 7 and 8 (group III) are predominantly found in presynaptic terminals (Conn and Pin, 1997; Niswender and Conn, 2010). Heterologous expression studies revealed that group-I mGluRs on one hand, and group II and III mGluRs on the other hand could form heterodimers (Doumazane et al., 2011), leading to the possible existence of 16 additional mGluRs with likely specific pharmacological and functional properties. Identifying such properties is a difficult issue to address, although one can expect that they will be essential in identifying the roles of mGlu heterodimers in vivo. What limits such studies is the presence of both homodimers and heterodimers in cells co-expressing both types of mGlu subunits (Kammermeier, 2012; Yin et al., 2014; Niswender et al., 2016).

Among heterodimeric mGluRs, mGlu2-4 was the most studied pair due to its important physiological interest and because different pharmacological tools are available (Kammermeier, 2012; Yin et al., 2014; Niswender et al., 2016). First, in the basal ganglia and the corticostriatal pathway, these two subunits are playing an important role in movement disorders such as Parkinson’s disease (Johnson et al., 2005). Second, previous immunohistochemistry and in situ hybridization studies suggest that mGlu2 and 4 receptors are co-localized in several brain regions (Neki et al., 1996; Bradley et al., 1999). Accordingly, mGlu2 and mGlu4 receptors could co-immunoprecipitate in native rodent tissue (Yin et al., 2014). However, it is difficult to detect pharmacological activation of any heterodimer using single activation protocols due to the co-existence of both homo- and heterodimers. Interestingly, the co-expression of the mGlu2 and mGlu4 subunits was reported to modify the pharmacology of mGlu2 and mGlu4 agonists. In addition, amongst the positive allosteric modulators (PAMs) of mGlu4 receptor, only one is active at the mGlu2-4 heterodimer (Yin et al., 2014; Niswender et al., 2016). The lack of effect of some mGlu4 PAMs in modulating mGlu4-mediated inhibition of cortico-striatal terminals was then used as a first evidence for the existence of mGlu2-4 heterodimers in the brain (Yin et al., 2014). However, in order to discern between homodimers and heterodimers, it is essential to find their specific pharmacological signature.

In the present study, we developed two different innovative approaches to characterize the pharmacological and functional properties of mGlu2-4 heterodimers specifically, without any influence of the co-existing mGlu2 and mGlu4 homodimers. We also used an innovative lanthanide-based time-resolved FRET microscopy approach (Faklaris et al., 2015), to demonstrate mGlu2 and mGlu4 can form heterodimers in transfected neurons. Using our heterodimer selective assays, we identified three pharmacological fingerprints that can be used to identify mGlu2-4 heterodimers in native cells. Such fingerprints could be identified in a neuronal cell line that naturally expresses both mGlu2 and mGlu4 subunits, as well as in the lateral perforant path (LPP) terminals in the hippocampus. These data bring strong evidence for the natural formation of such heterodimeric mGluRs in brain cells. Our observation will then be useful to study the function of mGlu2-4 heterodimers in the brain, and most importantly, to set up the condition to characterize other GPCR heterodimers.

Results

A FRET-based sensor to identify mGlu2-4-specific fingerprints

The co-expression of both rat mGlu2 and mGlu4 subunits led to the surface expression of three types of dimers: mGlu2 and mGlu4 homodimers and mGlu2-4 heterodimers. We then set up the transfection conditions to obtain an optimal expression at the cell surface of the mGlu2-4 heterodimer (Figure 1A; Figure 1—figure supplement 1A–C). To that aim, we co-transfected various amounts of plasmids encoding CLIP-mGlu2 and SNAP-mGlu4, and quantified the surface expression of each dimer population measuring lanthanide resonance energy transfer (LRET) in a time-resolved manner (TR-FRET) between a long life-time donor (Lumi4-Tb) and a fluorescent acceptor covalently attached to N-terminal CLIP (Gautier et al., 2008) or SNAP (Juillerat et al., 2003) tags (Maurel et al., 2008; Doumazane et al., 2011). This approach allowed the orthogonal labeling of the subunits in any dimer combinations. The use of donor and acceptor SNAP substrates allows the specific labeling of mGlu4 homodimers with a TR-FRET pair. Similarly, the use of CLIP substrates allows the measurement of TR-FRET signal originating from mGlu2 homodimers exclusively. Eventually, the use of a combination of SNAP-donor and CLIP-acceptor substrates leads to TR-FRET originating from the heterodimer only, since any homodimers carry either the donor or the acceptor but not both (Figure 1A; Figure 1—figure supplement 1A–C). Under these optimized conditions (40 ng CLIP-mGlu2 and 20 ng SNAP-mGlu4), the measure of the inhibition of cAMP production revealed a partial activity of the mGlu4 agonist L-AP4 and a slight loss in potency of the mGlu2 agonist LY379268 (Figure 1—figure supplement 1D–F), as reported by others (Yin et al., 2014). Of note, the L-AP4 dose-response curve in cells expressing both mGlu2 and 4 subunits can be fitted with a biphasic curve (Figure 1—figure supplement 1F), an effect consistent with the action of L-AP4 on both mGlu2-4 heterodimers and mGlu4 homodimers. This illustrates the difficulty of analyzing the specific properties of the heterodimer under such conditions.

Figure 1 with 3 supplements see all

Download asset Open asset

Phamacological profile of mGlu2, mGlu4 and mGlu2-4 receptors.

(A) Schematic representation of TR-FRET mGlu sensors generated. (**B–D**) Specific effect on VFT rearrangement of the CLIP-CLIP mGlu2 (red), SNAP-CLIP mGlu2-4 (green) or SNAP-SNAP mGlu4 (blue) with increasing concentrations of the indicated compound. (E) Schematic representation of the C1-C2 expression control system used for a specific expression of mGlu2-4 heterodimers at the cell surface. (**F–H**) Specific detection of the inhibition of cAMP pathway using C1-C2 expression control system for mGlu2-4 (green), as well as wild-type mGlu2 (red) and wild-type mGlu4 (blue). Glutamate and the specific agonists of mGlu4 (L–AP4) and mGlu2 (LY354740) present similar pharmacological profile using both techniques. Results are mean ± SEM from three independent experiments performed in triplicates.

https://doi.org/10.7554/eLife.25233.002

To examine the effect of various agonists specifically on the mGlu2-4 heterodimer at the surface of live cells, we took advantage of the large conformational change observed at the level of the extracellular domain of mGlu dimers upon activation. This conformational change led to a drastic decrease in TR-FRET signal (Doumazane et al., 2013) that can be followed specifically in any of the three types of dimers at the surface of cells co-transfected with CLIP-mGlu2 and SNAP-mGlu4 (Figure 1A). Of note, the properties of one dimer combination were then analyzed in the presence of the others.

In this assay, the glutamate potency was similar in the mGlu2-4 heterodimer and mGlu2 homodimer, higher than that on mGlu4 (Figure 1B), as previously reported by others (Yin et al., 2014). However, the potencies of the specific mGlu4 agonists were not increased in the heterodimer (Figure 1C; Figure 1—figure supplement 2A,B) suggesting no change in the affinity of mGlu4 ligands in the heterodimer. However, all of them acted as partial agonists within the mGlu2-4 heterodimer indicating that binding on the mGlu4 protomer only is not sufficient for a full activation of the heterodimer (Figure 1C; Figure 1—figure supplement 2A,B). This partial effect was more pronounced when activating mGlu4 with the partial agonist ACPT-I (Figure 1—figure supplement 2B). On the other hand, when activating the heterodimer with mGlu2 selective agonists, a loss in potency was observed in addition to the partial activity in the heterodimer (Figure 1D; Figure 1—figure supplement 2C,D). Interestingly, the mGlu2 agonist LY354740 displayed a strong loss in potency with a highly reduced slope (n_H = 0.29) on the heterodimer (Figure 1D). The agonist APDC showed a right-shifted curve in the heterodimer in comparison with mGlu2 homodimer and half of the maximal response (Figure 1—figure supplement 2C). Of interest, DCG-IV, a high-affinity mGlu2 agonist and low-affinity mGlu4 antagonist, displayed a biphasic dose-response curve with a reduction of the response at higher concentrations (Figure 1—figure supplement 2D).

Functional characterization of mGlu2-4 heterodimer confirmed the properties of LY354740

We next aimed at verifying that the pharmacological mGlu2-4 properties observed using the TR-FRET conformational sensor correlate with those measured using a functional read out. To that aim, we used a quality control system allowing the cell surface targeting of the mGlu2-4 heterodimer only (Figure 1E). We replaced the C-terminal tails of the SNAP-mGlu4 and CLIP-mGlu2 with a quality control system based on the modified intracellular tails of the GABA_B receptor subunits (called C1 and C2) (Brock et al., 2007; Huang et al., 2011). In that situation both homodimers are retained in the endoplasmic reticulum and do not reach the cell surface, and then are not capable of generating a signal as already reported for mGlu2 (Huang et al., 2011), and mGlu5 (Brock et al., 2007) receptors. In contrast, the coiled coil interaction between the C1 and C2-tails prevents the retention of each subunit, allowing the C1-C2 heterodimer to escape from the endoplasmic reticulum and reach the cell surface (Figure 1E) (see (Huang et al., 2011) for the characterization of the mGlu2-C1 and mGlu2-C2 constructs, and Figure 1—figure supplement 3 for the mGlu4-C1 construct). We set up the transfection conditions to avoid even the minimum leaking of the respective homodimers that might occur during the expression of these constructs and the absence of homodimer formation was checked by TR-FRET (Figure 1—figure supplement 3). The inhibition of forskolin-induced cAMP by mGlu2-4 receptors revealed data that perfectly match those measured with the TR-FRET sensor assay in terms of potency, Hill coefficient and efficacy (Figure 1F–H, Figure 1—figure supplement 2E–H). The potencies of the compounds tested in these different assays are indicated in Table 1.

Table 1

Potencies (pEC50) of the indicated compound on mGlu2, mGlu2-4 and mGlu4 as determined using the TR-FRET-based conformational assay depicted in Figure 1, or the cAMP assay as depicted in Figure 2.

Data are means ± SEM of at least three experiments performed in triplicates.

https://doi.org/10.7554/eLife.25233.006

	Receptor
	mGlu2		mGlu2-4		mGlu4
Compound	pEC₅₀	Emax (%)	pEC₅₀	Emax (%)	pEC₅₀	Emax (%)
TR-FRET conformational sensor assay
Glutamate	5.5 ± 0.04	100 ± 1	5.5 ± 0.04	100 ± 2	5.0 ± 0.05	100 ± 2
LY354740	7.5 ± 0.04	99 ± 1	5.9 ± 0.4	59 ± 9	–	3 ± 1
APDC	5.4 ± 0.07	97 ± 3	4.8 ± 0.1	56 ± 3	–	10 ± 2
DCGIV	6.6 ± 0.06	79 ± 1	6.8 ± 0.2	27 ± 2	4.7 ± 0.2	16 ± 2
L-AP4	–	14 ± 1	6.7 ± 0.1	48 ± 2	6.2 ± 0.1	95 ± 1
LSP4-2022	–	3 ± 1	6.8 ± 0.1	53 ± 2	6.2 ± 0.05	91 ± 1
ACPT-I	–	7 ± 1	5.3 ± 0.2	26 ± 2	5.1 ± 0.1	74 ± 2
cAMP assay
Glutamate	5.6 ± 0.06	100 ± 1	5.6 ± 0.07	100 ± 2	5.1 ± 0.06	100 ± 2
LY354740	7.8 ± 0.04	97 ± 1	5.3 ± 0.2	59 ± 6	–	5 ± 1
APDC	6.3 ± 0.06	96 ± 2	6.0 ± 0.1	50 ± 2	–	9 ± 3
DCGIV	6.9 ± 0.09	69 ± 2	6.6 ± 0.2	42 ± 3	5.5 ± 0.4	−17 ± 2
L-AP4	–	3 ± 2	6.6 ± 0.1	51 ± 1	6.9 ± 0.07	99 ± 2
LSP4-2022	–	3 ± 1	6.9 ± 0.1	48 ± 2	6.9 ± 0.04	99 ± 1
ACPT-I	–	13 ± 7	5.9 ± 0.1	32 ± 1	5.7 ± 0.1	68 ± 3

Taken together, these data revealed a low potency and low Hill coefficient for LY354740, which can be used as a first fingerprint of the mGlu2-4 heterodimer.

Activation of both subunits in mGlu2-4 receptor is required for full activity

To examine the role of each binding site in the activation of an mGlu heterodimer, we examined the effect of mutating either site. The substitution of the conserved Tyr and Asp by Ala in the glutamate-binding site of mGlu receptors (position 216 and 295 in mGlu2), which are called YADA mutants (Kniazeff et al., 2004; Brock et al., 2007; Doumazane et al., 2013), strongly impairs the binding of agonists. Accordingly, mGlu2^YADA and mGlu4^YADA homodimers, as well as the mGlu2^YADA-4^YADA heterodimer could not be activated by glutamate (Figure 2), despite their normal expression at the cell surface (Figure 2—figure supplement 1B,C). However, when a single subunit per heterodimer is mutated, glutamate maximal FRET change was about half the maximal response of the wild-type receptor (Figure 2B), consistent with a full activation requesting both binding sites occupied. As expected, no effect of the selective mGlu2 agonists could be observed in the mGlu2^YADA-4 heterodimer (Figure 2D–F). Similarly, selective mGlu4 agonists had no effect on the heterodimer mutated in the mGlu4 subunit (mGlu2-4^YADA) (Figure 2C, Figure 2—figure supplement 1A), but mGlu2 selective agonists retained their activity (Figure 2D–F) on this mutant heterodimer. Of note, the Hill coefficient of LY354740 was increased to n_H = 0.89 (Figure 2E), and the DCG-IV dose response curve was no longer biphasic, the decreased response obtained at higher dose not being observed in this mutated heterodimer (Figure 2F). These findings demonstrate the importance of an intact mGlu4-binding site in the complex pharmacological effect of the mGlu2 agonists LY354740 and DCG-IV on the mGlu2-4 heterodimer.

Figure 2 with 1 supplement see all

Download asset Open asset

Role of each binding site in agonist-induced activity of mGlu2-4.

(A) Schematic representation of the mGlu2-4 mutants; wild type (green), YADA mutation in mGlu2 (blue), mGlu4 (red) or both (black). (**B–F**) Effect of increasing concentrations of the indicated ligands on the mGlu2-4 TR-FRET sensor. Results are mean± SEM of three independent experiments performed in triplicates.

https://doi.org/10.7554/eLife.25233.007

Cooperativity between the agonist binding sites in the mGlu2-4 heterodimer

The above data prompted us to examine the influence of agonist binding in one subunit on the effect mediated through agonist binding in the second subunit. We observed that agonist binding in mGlu4 receptor increased the potency of mGlu2-specific ligands on the mGlu2-4 heterodimer (Figure 3A–C, Figure 3—figure supplement 1). In the case of LY354740, not only mGlu4 agonists increased its potency, but also restored an nH close to unity, as observed with the TR-FRET conformational sensor (Figure 3A) and cAMP assays (Figure 3B,D). These results revealed a crosstalk between mGlu4 and mGlu2 protomers within the heterodimer, an effect that can be observed both at the level of the ECDs as revealed by the TR-FRET sensor, and at the G-protein coupling site.

Figure 3 with 2 supplements see all

Download asset Open asset

Synergistic action of mGlu2 and mGlu4 agonists in mGlu2-4 heterodimer.

(**A–B**) Dose response curves of the mGlu2 ligand LY354740 in absence or presence of mGlu4 ligands (ACPT-I 10 µM, APcPr 3 µM, L-AP4 3 µM or LSP4-2022 3 µM) on the TR-FRET assay (A) and the inhibition of forskolin-induced cAMP production (B). (**C–D**) LY35740 EC₅₀ (C) or Hill slope (D) in the presence of the indicated concentration of L-AP4. Results are mean ± SEM of three independent experiments performed in triplicates. Curve fitting was performed by using nonlinear regression. p-values<0.05 were considered statistically significant (*).

https://doi.org/10.7554/eLife.25233.009

To quantitatively analyze the effect of mGlu4 agonists on the nH of LY354740 between ligands in the mGlu2-4 heterodimer, mathematical models were developed (Figure 3—figure supplement 2A,B). Assuming that mGlu2 and mGlu4 ligands bind to their respective protomers exclusively (Model 1), the mechanistic model collapses to an empirical model that can be expressed as a Hill equation with a Hill coefficient of 1 (see Appendix 1). This is not consistent with the flat slope curves displayed by LY354740 and the cooperativity exerted by mGlu4 ligands. In a second model, we then assumed that LY354740 could bind to mGlu4 VFT, with a very low affinity in mGlu4 homodimers, but with a higher affinity in the mGlu4 subunit of the mGlu2-4 heterodimer due to its binding on mGlu2 VFT. For a best fitting of the LY354740 curve (Figure 3—figure supplement 2C), we had to assume the closed-closed state is not fully achieved because the ligand behaves as a partial agonist. In addition, two components of the functional activity had to be set up, one related with the binding to a first VFT of the heterodimer and another one related to the binding to the second protomer. Eventually, we had to assume that LY354740 binds the mGlu4 protomer after occupying first the mGlu2 binding site. This model is then consistent with L-AP4 binding in the mGlu4 VFT increasing LY354740 potency and restoring a Hill coefficient close to unity (see Appendix 1).

The synergistic effect between mGlu2 and mGlu4 agonists constitutes a second fingerprint of the mGlu2-4 heterodimer that may be useful for the identification of mGlu2-4 in neurons.

Allosteric modulation of mGlu2-4 heterodimer

Positive allosteric modulators (PAMs) can enhance both agonist affinity and efficacy. They can also have an intrinsic agonist activity on mGluRs (Conn et al., 2014). Using our TR-FRET mGluR conformational sensors, we found that the mGlu2 PAMs, BINA and LY487379 potentiate the effect of glutamate (at its EC₂₀) in the mGlu2 homodimer, but very weakly in the mGlu2-4 heterodimer and not at all in mGlu4 (Figure 4A). Regarding the mGlu4 PAMs, VU0155041 activated both mGlu4 and mGlu2-4 while VU0415374 potentiated mGlu4 homodimer mainly (Figure 4A), as previously reported (Yin et al., 2014; Niswender et al., 2016). Co-application of mGlu2 and mGlu4 PAMs led to specific effects on the mGlu2-4 heterodimers, depending on the PAM used. Neither BINA nor VU0415374 had any effect when applied alone on the glutamate EC₂₀ mediated response (Figure 4B–D). However, a clear and strong potentiation of the glutamate EC₂₀ response was observed when both PAMs were applied together (Figure 4B–D). This synergistic effect was not observed with another combination of PAMs (Figure 4B). The synergistic action of BINA and VU0415374 was also observed in a functional cAMP assay (Figure 4D) and constitutes therefore a third fingerprint for the mGlu2-4 heterodimers.

Figure 4

Download asset Open asset

Synergistic action of mGlu2 and mGlu4 PAMs in mGlu2-4 heterodimer.

(A) Effect of mGlu2 (LY487379 10 µM, BINA 10 µM) and mGlu4 (VU0155041 10 µM, VU0415374 10 µM) PAMs on mGlu2 (red), mGlu4 (blue) and mGlu2-4 (green) TR-FRET sensors in the presence of an EC₂₀ of glutamate. (B) Effect of mGlu2 and/or mGlu4 PAMs on the response mediated by an EC₂₀ concentration of glutamate in mGlu2-4 heterodimer by TR-FRET. The strong synergy between BINA and VU0415374 is highlighted with red bars. (**C–D**) Dose response of BINA and/or VU0415374 in potentiating the effect of EC₂₀ glutamate on TR-FRET sensors (C) and cAMP inhibition (D). Results are mean ± SEM of three independent experiments performed in triplicates.

https://doi.org/10.7554/eLife.25233.012

mGlu2 and mGlu4 can form heterodimers in neurons

In primary cultures of hippocampal neurons, co-expression of CLIP-mGlu2 and SNAP-mGlu4 subunits could be detected at the cell surface through labeling with CLIP and SNAP substrates carrying either Lumi4-Tb or Red. In these neurons, using a lanthanide-based time-resolved FRET microscope that we recently developed (Faklaris et al., 2015), we detected a TR-FRET signal between the CLIP and SNAP subunits equivalent to that measured for homodimeric mGlu2 receptors. This observation is consistent with the formation of mGlu2-4 heterodimers in transfected neurons (Figure 5). In contrast, no TR-FRET could be detected between CLIP-mGlu2 and SNAP-mGlu1 (Figure 5), two subunits reported not to associate into heterodimeric entities (Doumazane et al., 2011; Levitz et al., 2016). Such mGlu2-4 heterodimers are not the consequence of a large over-expression of these tagged subunits, since their quantification using a mGlu2 specific antibody, relative to the endogenous mGlu2 revealed a five times higher expression of the tagged receptor only (Möller et al., manuscript submitted for publication).

Figure 5

Download asset Open asset

TR-FRET detection of mGlu2-4 heterodimers in transfected hippocampal neurons.

Neurons transfected with CLIP-mGlu2 and SNAP-mGlu4 are compared with either CLIP-mGlu2 and SNAP-mGlu1 (negative control) or SNAP-mGlu2 (positive control). The receptors are labeled with Lumi4-Tb as donor and Red as acceptor. (A) Image examples of neurons expressing the three receptor combinations at comparable expression levels showing similar TR-FRET and NFRET (TR-FRET normalized to the expression of donor and acceptor) signals for the mGlu2-4 heterodimer and the mGlu2 homodimer and a low signal for the mGlu1-2 heterodimer. TR-FRET images are corrected for bleedthrough and thresholded to remove background and noise. Images are 138 µm wide. (B) Quantification of NFRET as a function of the expression level of donor and acceptor. Each point is the quantification of one neuron. (C) Scatter plot of NFRET.

https://doi.org/10.7554/eLife.25233.013

Functional mGlu2-4 heterodimers in a neuronal cell line

As a first attempt to identify native mGlu2-4 heterodimers in non-transfected cells, we used the striatal cell line STHdh^Q7 (Trettel et al., 2000), since both mGlu2 and mGlu4 mRNA have been reported in the striatum, although to a low level (Ohishi et al., 1993; Conn et al., 2005; Ferraguti and Shigemoto, 2006; Gu et al., 2008; Beurrier et al., 2009). We first examined whether mGlu2-4 heterodimers can form in these cells upon transfection of CLIP-mGlu2 and SNAP-mGlu4. The formation of mGlu2-4 heterodimers was well illustrated by the cooperativity between the binding sites characteristic of the heterodimer, as revealed with the TR-FRET sensor assay (Figure 6—figure supplement 1). We then examined whether a native mGluR-mediated responses with the characteristics of the mGlu2-4 heterodimer could be detected in these cells. Unfortunately, the endogenous receptor levels are very low and the detection of their activation could not be detected using cAMP assays. However, we achieved measuring responses with mGlu2 and mGlu4 agonists using the xCELLigence technique, a label-free method reporting on small variations in cell shape (cell index). STHdh cells are shrunk when adenylate cyclase is activated (Figure 6A). Taking advantage of this characteristic we observed that mGlu2- and mGlu4-specific ligands impaired the forskolin effect in a pertussis toxin-sensitive way, consistent with the presence of endogenous Gi-coupled mGlu2 and mGlu4 receptors in these cells (Figure 6—figure supplement 2). Both LY354740 and the mGlu4 preferential agonist LSP4-2022 (Goudet et al., 2012) produced a dose-dependent effect with the expected EC₅₀s (Figure 6B).

Figure 6 with 2 supplements see all

Download asset Open asset

Functional evidence for mGlu2-4 heterodimers in a neuronal cell line, STHdh.

(A) Representative image of neuronal cells treated or untreated with forskolin. Cells were transfected with GFP (green) and stained with MAP2 (red) and DAPI (blue). (B) Dose response on LY354740 or LSP4-2022 represented in % variation of cell index versus LY354740. (C) pEC50 values of LY354740 as determined using the change in cell index, in control cells transfected with GFP, and in cells transfected with the SH RNA against mGlu4, under control condition (gray bars) or in the presence of the mGlu4 agonists LSP4-2022 (10 µM). (D) Maximum effect of mGlu2 (BINA 1 µM) and/or mGlu4 PAM (VU0415374 1 µM) in potentiating the effect of low concentration of LY354740 (10 nM) and LSP4-2022 (100 nM) in control STHdh cells infected with GFP vector or silencing shRNA for mGlu4. Data in **B–D** are means ± SEM of three independent experiments performed in triplicates.

https://doi.org/10.7554/eLife.25233.014

As described above with the mGlu2-4 heterodimer (Figure 3), LSP4-2022 increased the potency of the mGlu2 agonist LY354740 (Figure 6B,C). This mGlu4 effect is clearly due to the presence of mGlu4 in these STHdh cells, since depletion of mGlu4 with a lentivirus expressing mGlu4 ShRNA, resulted in a higher potency of LY354740 consistent with its potency on mGlu2 homodimers, with no further effect of the mGlu4 agonist (Figure 6C). Furthermore, the synergistic effect of BINA and VU0415374 observed on the mGlu2-4 heterodimer (Figure 4) could also be observed in the STHdh cells (Figure 6D). Indeed, either BINA (1 µM) or VU0415374 (1 µM) modestly potentiated the effect of agonists (a combination of 10 nM LY354740 and 100 nM LSP4-2022), while a strong potentiation was observed with these two PAMs added simultaneously (Figure 6D). No such effect was observed in the mGlu4-silenced cells (Figure 6D). These data are consistent with the existence of mGlu2-4 heterodimers endogenously expressed in the STHdh cells.

Pharmacological evidence for mGlu2-4 receptors in lateral perforant path terminals

In the hippocampus, mGlu4 receptor expression is prominent in the inner third of the molecular layer of the dentate gyrus (Shigemoto et al., 1997; Corti et al., 2002). mGlu2 receptor is also expressed in the molecular layer of the dentate gyrus (Ohishi et al., 1993; Gu et al., 2008; Wright et al., 2013). As expected, activation of mGlu2 (concentration >300 nM LY354740) or mGlu4 receptor (concentration >5 µM LSP4-2022) inhibited synaptic transmission at the LPP (Figure 7). The effect of LSP4-2022 is absent in slices prepared form mGlu4 KO mice, demonstrating its effect is mediated by mGlu4 in control animals (Figure 7)

Figure 7 with 1 supplement see all

Download asset Open asset

Effect of LY354740 and LSP4-2022 in the LPP of wild-type (WT) mice.

(A) Bar graph illustrating the % inhibition of fEPSPs induced by low (100 nM and 300 nM) and high (5 µM and 10 µM) concentrations of LSP4-2022 in the LPP. Only high concentrations of LSP4-2022 induced a significant decrease of fEPSP amplitude. Green bars indicate data obtained using slices from mGlu4 KO mice. (B) Inhibitory effect of LY354740 on fEPSP amplitude in the LPP. Note that 300 nM or 1 µM LY354740 caused a significant decrease of fEPSP amplitude. Green bars indicate data obtained using slices from mGlu4 KO mice. (C) Representative averaged traces of evoked synaptic activity induced by LPP stimulation in field recording of granular cells from WT mice (*Left*). Bar graph illustrating the % inhibition of fEPSP amplitude by LY354740 (10 nM), LSP4-2022 (100 nM) and LY354740 (10 nM) / LSP4-2022 (100 nM) in the LPP (*Right*). Note that application of LY354740 (10 nM) + LSP4-2022 (100 nM) significantly decreased the fEPSP amplitude. Data in A-C are means ± SEM of (n) independent experiments from at least three different animals. *p<0.05, **p<0.001.

https://doi.org/10.7554/eLife.25233.017

In order to investigate whether mGlu2-4 heterodimers are expressed in these synapses, we applied a low concentration of each agonist that produced no detectable inhibitory effect in the LPP. Only when co-applied, LSP4-2022 (100 nM, 3.73 ± 2.19%, n = 4) and LY354740 (10 nM, 2.62 ± 3.83%, n = 5) induced a significantly large reduction of the fEPSPs in LPP (20.05 ± 5.82%; n = 6; Figure 7C). No such effect was observed in slices prepared from mGlu4 KO mice (Pekhletski et al., 1996) (Figure 7). The synergy of these agonists is consistent with the presence of mGlu2-4 heterodimers in these terminals. However, one cannot exclude a possible synergy between mGlu2 and mGlu4 at the signaling level, rather than within a heterodimer. We think this is unlikely since such strong synergistic effect have not been observed between Gi-coupled receptors. Indeed, upon co-expression of mGlu2 and the Gi-coupled delta opioid receptor (DOR), activation of DOR with SNC162 had no effect on the potency of LY354740 in inhibiting cAMP formation via mGlu2 receptors (Figure 7—figure supplement 1).

Discussion

Despite their description in heterologous cells 5 years ago (Doumazane et al., 2011), evidence for the existence of mGlu heterodimers in vivo remains elusive. Using two different approaches to characterize mGlu2-4 heterodimers specifically, we identified pharmacological fingerprints of such receptors. First, the mGlu2 selective agonist LY354740 behaves differently on the mGlu2-4 heterodimer than on the mGlu2 homodimer, including a lower potency, and a lower Hill coefficient. Such complex properties of LY354740 disappeared in the presence of mGlu4 agonists. Second, mGlu4 agonists largely increase the potency of LY354740. Third, among four mGlu2 or mGlu4 PAMs tested, only VU0155041 potentiated the effect of agonists on the mGlu2-4 heterodimer. Fourth, a combination of two PAMs (BINA and VU0155041) inactive when applied alone enhanced agonist action on the heterodimer when applied together. Such pharmacological fingerprints provide ways for demonstrating the existence of such heterodimers in native cells, as illustrated here with a neuronal cell line and the medial perforant path terminals in the dentate gyrus.

One major difficulty in studying the functional and pharmacological properties of GPCR heterodimers is the ability of each subunit to form functional receptors, making difficult the measurement of signals originating from the heterodimers exclusively. In previous studies, properties of the mGlu2-4 heterodimers were studied in cells co-expressing both mGlu2 and mGlu4 (Yin et al., 2014). However, as illustrated in Figure 1—figure supplement 1, even when conditions were used for an optimal expression of the heterodimer at the cell surface, data obtained were always contaminated by responses mediated by the homodimers. In this study, we developed two different approaches that allowed the analysis of the pharmacological and functional properties of mGlu2-4 heterodimers, both strategies being likely useful for other class C heterodimers. Of note, the TR-FRET sensor assay that relies on the inter-subunit movement in mGlu dimers allows the specific analysis of compounds in any of the three combinations specifically. This allowed us to identify specific properties of the heterodimers that can be useful for the characterization of such receptors in native cells.

Our data show that agonist occupancy of both subunits is required for a full activity of the heterodimer, as well illustrated using specific agonists of one subunit (either mGlu2 or mGlu4), or by mutating the binding site of either subunit. Accordingly, activating either the mGlu2 or mGlu4 VFT in the mGlu2-4 heterodimer leads to a similar partial effect, both at the conformational level of the VFT dimer, as revealed with the TR-FRET sensor assay, and at the signaling level. This finding is consistent with previous studies demonstrating that two agonists are required, with both VFT closed to fully activate mGlu5 homodimers (Kniazeff et al., 2004). However, it was reported that the low activity observed when a single protomer is occupied by an agonist in mGlu2 homodimers possibly results from the spontaneous closure of the second, unliganded VFT (Levitz et al., 2016). It is therefore possible that part of the activity observed with mGlu2 selective agonists on mGlu2-4 is due to a spontaneous closure of the mGlu4 VFT. This is consistent with the biphasic curve of DCG-IV, an mGlu2 agonist that has mGlu4 antagonist activity at high concentration.

The effect of LY354740, a well-known group-II-specific mGluR agonist, appears quite complex on the mGlu2-4 heterodimer. This compound has a partial efficacy, but a low Hill coefficient on the heterodimer. Its mGlu2 potency and normal Hill coefficient are restored when the mGlu4 subunit is either activated or mutated to prevent ligand binding. These findings cannot be explained by a simple model in which both mGlu2 and mGlu4 ligands bind selectively to their respective subunit. Instead, our data suggest that LY354740 can bind to both subunits, its interaction with mGlu2 increasing its affinity to mGlu4 in the heterodimer. Such ligand cooperativity between bindings sites in an mGlu dimer has already been reported in mGlu5 receptors in which one binding site is mutated (Kniazeff et al., 2004; Rovira et al., 2008). The synergistic effect observed between mGlu2 and mGlu4 agonists appears as an interesting property to identify mGlu2-4 heterodimers in neurons.

When examining the effect of two mGlu2 (BINA and LY487379) and two mGlu4 PAMs (VU0155041 and VU0415374), we confirmed that VU0155041 is the only PAM able to potentiate the effect of glutamate on the mGlu2-4 heterodimer (Yin et al., 2014). Although BINA and VU0415374 had very modest effects alone, their co-application largely potentiated the heterodimer. The structural basis for this synergistic effect, especially when considering that a single 7TM is active at a time in mGlu dimers (Hlavackova et al., 2005, 2012), remains unclear and will be the subject of further studies. Whatever the reason, this synergistic effect between these two PAMs offers another way to identify mGlu2-4 heterodimers in neurons.

A recent study revealed that mGlu2-4 heterodimers are likely present in cortico striatal terminals (Yin et al., 2014). This conclusion is based on co-immunoprecipitation data, and on the lack of effect of PHCCC, an mGlu4 PAM devoid of effects in cells co-expressing mGlu2 and mGlu4, while VU0155041 that is active on both mGlu4 and mGlu2-4 heterodimers, potentiated the response (Yin et al., 2014). Our data also revealed that mGlu2-4 can form in transfected neurons, indicating there are no specific mechanisms in neurons that would prevent the formation of such heterodimers. Most importantly, we found that in a neuronal cell line, responses with the pharmacological characteristics of the mGlu2-4 heterodimers can be recorded. Indeed, the synergistic effects of mGlu2 and mGlu4 ligands (both agonists and PAMs) typical of the mGlu2-4 heterodimer were observed in these cells, bringing strong evidence that endogenous mGlu2-4 heterodimers exist in these neuronal cells despite the low expression of both mGlu2 and mGlu4. The synergistic activity of the agonists LY354740 and LSP4-2022 was also observed in the terminals of the medial perforant path in the dentate gyrus where both mGlu2 and mGlu4 subunits are present (Shigemoto et al., 1997). Such a synergistic activity is no longer observed in slices prepared from mGlu4 KO mice, demonstrating the involvement of mGlu4. However, we cannot rule out that such a synergy may come from the signaling cascades activated by both mGlu2 and mGlu4 homodimers. We still think this is unlikely because such a strong synergistic effect has not been observed between Gi-coupled receptors and indeed could not be observed between mGlu2 and the delta opioid receptor co-expressed in the same cells.

It is sometimes argued that GPCR dimers and heterodimers result from the overexpression of the partners. For several reasons, this is unlikely the case for the mGlu2-4 heterodimer. Over-expression is expected to result in larger mGluR complexes since mGluRs are constitutive and covalent dimers (Calebiro et al., 2013), and no proximity could be detected by FRET between mGlu2 and mGlu1 receptors known not to form heterodimers (Doumazane et al., 2011), even in transfected neurons. The FRET efficacy between mGlu2 and mGlu4 largely decrease upon receptor activation, as expected for a correctly assembled dimer. The relative quantification of transfected over endogenous receptors revealed a fivefold only over-expression (Møller et al., manuscript in preparation). Various approaches have been used to estimate the size of the mGlu2-4 complex and all revealed a strict dimer (Doumazane et al., 2011). Eventually, a receptor with the pharmacological fingerprints of the mGlu2-4 heretodimer was observed in STHdh cells where both mGlu2- and mGlu4-mediated responses were difficult to detect, suggesting a low expression level.

Taken together, our data add to previous studies suggesting the existence of mGlu heterodimers in the brain. We show that mGlu2-4 receptors likely exist in the brain and we report innovative approaches that will be useful to confirm the existence of other mGlu heterodimers. For example, one may propose the existence of heterodimers containing an mGlu7 subunit, for which the very low glutamate potency raised a number of questions regarding its roles in vivo. Indeed, mGlu7 can be found with other high-affinity mGluRs, such as mGlu8 (Ferraguti et al., 2005), offering a way to involve mGlu7 in a receptor heterodimer with specific properties. Proteomic experiments also identified mGlu5 as a partner of mGlu1 (Pandya et al., 2016), a finding that likely explains surprising functional studies regarding the effect of specific mGlu1 and mGlu5 inhibitors in the hippocampus (Huber et al., 2001; Volk et al., 2006). Functional studies also suggested mGlu3 receptors could be involved in heterodimeric entities with mGlu2 receptors (Iacovelli et al., 2009). A clear view of such mGlu heterodimers is definitively needed since all possible combinations observed in transfected cells suggest the existence of 16 additional receptor entities in the brain. In addition, first results are already highlighting specific roles played by homo and heterodimers in the actions of drugs with therapeutic potentials. For example, PAMs selective for the homodimeric mGlu4 may be preferred for the treatment of Parkinson disease (Niswender et al., 2016). Our study highlighting techniques to decipher the specific properties of mGlu heterodimers will definitively help solving these important issues.

Materials and methods

Materials

SNAP-Lumi4-Tb, SNAP-Green and CLIP-Green were from Cisbio Bioassays (Codolet, France). SNAP-block and CLIP-block were from New England Biolabs (Ipswich, MA, USA). (2R,4R)-APDC, ACPT-I, L-AP4, BINA, DCG-IV, LY341495, LY354740, LY487379, VU0155041 and SNC162 were purchased from Tocris Bioscience (Bristol, UK). LSP4-2022 was kindly provided by Dr. Francine Acher (Université Paris Descartes, France). VU0415374 was synthesized by Dr. Xavier Gómez and provided by Dr. Amadeu Llebaria (University of Barcelona, Spain). Control GFP and ShmRNA mGlu2 and mGlu4 lentiviral particles were purchased from Santa Cruz Biotechnology Inc. (Dallas, TX, USA). The mGlu4 KO mice were obtained from Dr David Hampson (Toronto, Canada [Pekhletski et al., 1996]), and their genotype determined as reported (Pitsch et al., 2007).

Plasmids

The pRK5 plasmids encoding the wild-type rat mGlu subunits in which the SNAP or CLIP has been inserted at their N-term after the signal peptide, and constructs with YADA mutations in mGlu2 were previously described (Doumazane et al., 2011, 2013). The pRK5 plasmid encoding for the ligand binding deficient SNAP-mGlu4-YADA mutant in which the two residues Y230 and D312 important for agonist binding in the VFT were mutated, was generated by site-directed mutagenesis using QuikChange mutagenesis protocol (Agilent Technologies, Les Ulis, France) using the SNAP-mGlu4 plasmid as a template (Doumazane et al., 2011). The sequence coding C1 (the 47-residue coiled-coil sequence of the C-terminal of GABA_B1), or C2 (the 49-residue coiled-coil region of GABA_B2), followed by the endoplasmic reticulum retention signal KKTN, as previously described (Huang et al., 2011), was used to generate the constructs SNAP-mGlu4-C1 and CLIP-mGlu2-C2. SNAP-mGlu4-C1 was obtained by replacing the last 38 residues in mGlu4 C-term (SNAP-tagged version of mGlu4 was used) by C1KKTN. In this construct, the C-term of SNAP-mGlu4-C1 is …NKFTTGSSTNNNEEEKSRLLEKENRELEKIIAEKEERVSELRHLQSRQQLKKTN (the last residues (up to Thr874) of mGlu4 are underlined, those of C1 are in italic). The C-term sequence of CLIP-mGlu2-C2 was previously described (Huang et al., 2011). The plasmid encoding SNAP-delta opioid receptor was from Cisbio Bioassays.

Cell culture and transfection

Request a detailed protocol

HEK293 cells (ATCC, CRL-1573, lot: 3449904) were cultured in Dulbecco’s modified Eagle’s medium (Thermo Fischer Scientific, Courtaboeuf, France) supplemented with 10% (vol/vol) fetal bovine serum (Sigma Aldrich) in a P2 cell culture room. Absence of mycoplasma was routinely checkedusing the MycoAlert Mycoplasma detection kit (LT07-318 (Lonza, Amboise, France), according to the manufacturer protocol. HEK 293 cells were used after 35 to 40 passages and transfected with a reverse transfection protocol using Lipofectamine 2000 (Thermo Fischer Scientific, Courtaboeuf, France), and finally plated in polyornithine-coated, black-walled, dark-bottom, 96-well plates at 10⁵ cells/well. To avoid too high concentrations of glutamate in the assay medium that could interfere with mGluR activity, cells were cotransfected with the plasmid encoding the glutamate transporter EAAC1 and incubated in DMEM Glutamax medium (Thermo Fischer Scientific) at least 2 hr before the different assays were performed. Frozen-labeled HEK-293 cells were transfected as described above, labeled as described below, then frozen at −80°C with 10% DMSO and fetal bovine serum, and later washed three times in Krebs buffer (10 mM Hepes pH 7.4, 146 mM NaCl, 4.2 mM KCl, 1 mM CaCl₂, 0.5 mM MgCl₂, 5.6 mM glucose, bovine serum albumin (BSA) 0.1%) before use.

In order to optimize the best expression of mGlu2-mGlu4 heteromers, several ratios of mGlu2:mGlu4 were assayed. It was determined by TR-FRET analysis that 2:1 ratio (40 ng CLIP-mGlu2: 20 ng SNAP-mGlu4) was optimal for the detection of all populations (Figure 1—figure supplement 1A–C). Using these conditions, a large batch of cells were transfected, labeled and frozen to perform a complete screening of the different compounds in 384-well plates.

Conditionally immortalized wild-type STHdh^Q7 striatal neuronal progenitor cell line (Trettel et al., 2000) were kindly provided by Dr Sílvia Ginés (University of Barcelona, Spain). These cells nicely differentiated and became MAP2 positive when cultured in a differentiated medium as described (Trettel et al., 2000). We also verified that they were still responsive to dopamine D1 and histamine H3 receptor agonists using the Xcellingence technology. Neuronal cells were grown at 33°C in DMEM (Sigma-Aldrich), supplemented with 10% fetal bovine serum (FBS), 1% streptomycin-penicillin, 2 mM L-glutamine, 1 mM sodium pyruvate, and 400 μg/ml G418 (Geneticin; Invitrogen). Neuronal cells were transfected with Lipofectamine LTX (Thermo Fischer Scientific, Courtaboeuf, France) following the protocol from the provider. To perform silencing of mGlu2 and mGlu4, STHdh cells were infected with control GFP vector, ShmRNA mGlu2 or mGlu4 vector and after 48 hr, infected cells were selected by adding hygromycin-containing medium.

Fluorescence labeling and TR-FRET measurements

Request a detailed protocol

SNAP-tag labeling alone and orthogonal labeling of SNAP- and CLIP-tag were performed as described previously (Scholler et al., 2017). Briefly, for SNAP-tag labeling, 24 hr after transfection, HEK293 cells were incubated at 37°C for 1 hr with a solution of 100 nM of SNAP-Lumi4-Tb, 60 nM of SNAP-Green and 1 μM CLIP-block, in case of FRET detection between SNAP-tag subunits. For CLIP labeling, cells were incubated with 1 μM CLIP-Lumi4-Tb, 800 nM CLIP-Green and 1 μM SNAP-block. For co-labeling of the SNAP- and CLIP tags, cells were incubated at 37°C for 2 hr with a solution of 300 nM SNAP-Lumi4-Tb and 1 μM CLIP-Green. After labeling, cells were washed three times with Krebs buffer, and drugs were added. Then, the TR-FRET measurements were performed on a PHERAstar FS microplate reader (BMG Labtech, Ortenberg, Germany) which is standardly equipped with ‘TR-FRET’ optical modules and two photomultiplier tubes to detect two emission wavelengths representing donor and acceptor emission simultaneously, as previously described (Scholler et al., 2017). To monitor the emissive decay curves, the Lumi4-Tb present in each well was excited using N₂ laser emission line at 337 nm (40 flashes per well for the 96-well plate format, 20 flashes per well for the 384-well plate format). The emission decay was collected during 2500 or 5000 μs with 5 μs or 10 μs steps, respectively, at 620 nm for the donor (Lumi4-Tb) and at 520 nm for Green, as can be indicated in the ‘advanced mode’ option of the plate-reader’s software. For acceptor ratio determination, optimal integration windows were determined as previously reported (Scholler et al., 2017). The acceptor ratio was calculated using the sensitized acceptor signal integrated over the time window [50 μs-100 μs], divided by the sensitized acceptor signal integrated over the time window [800 μs-1200 μs].

cAMP functional assay

Request a detailed protocol

The amount of cAMP was determined using the Glosensor^TM cAMP assay (Promega Corporation, Madison, USA). HEK293 cells were co-transfected with the indicated mGluR plasmids and the pGloSensor-22F plasmid. The day after, cells were starved during 2 hr in serum-free medium and afterwards incubated in Krebs buffer with 450 μg/ml luciferin (Sigma-Aldrich) during 1 hr. The luminescence peak signal was measured on a Mithras microplate reader at 28°C during 8 min since luminescence signal was stable. Then, forskolin plus mGluR ligands were added and luminescence was measured for 30 min.

Label-free impedance assay

Request a detailed protocol

xCELLigence plates were coated with poly-ornithine and laminin during 1 hr, and neuronal cells were seeded at 3 × 10⁴ cells/well and introduced into the incubator at 33°C overnight. Medium was replaced by serum-free medium during 2 hr to reach a stable cell index, and then forskolin and mGluR ligands were added and the signal was followed during at least 2 hr using xCELLigence RTCA DP apparatus (ACEA Bioscience Inc, San Diego, USA). When antagonists or PAMs were used, they were added 20 min before forskolin. Pertussin toxin was added 4 hr after cell plating and incubated overnight.

Neuronal culture and TR-FRET microscopy

Request a detailed protocol

Hippocampi from Sprague-Dawley rat embryos on embryonic day 18 (E18) were dissected, dissociated by treatment with liberase TL (Roche, Boulogne-Billancourt, France), then mechanical triturated and plated on Lab-Tek II chambered cover slides (Thermo Fisher Scientific) coated with poly-L-ornithine and laminin (Sigma-Aldrich) at a density of ~300 neurons/mm². Neurons were cultured in Neurobasal medium (Thermo Fisher Scientific) supplemented with 2% B-27 (Thermo Fisher Scientific), 100 U/ml Penicillin-Streptomycin (Thermo Fisher Scientific), 10 mM HEPES, and 0.5 mM Glutamax medium (Thermo Fisher Scientific). 0.5 mM L-glutamine was added when plating the cells. Half of the medium was exchanged weekly. Neurons were transfected with Lipofectamine 2000 at 10 days in vitro (DIV). The medium was exchanged after 4 hr of incubation with the transfection reagent with half fresh medium and half medium conditioned by incubation with primary neurons. pRK5 plasmids for expression of SNAP- or CLIP-tagged rat mGluRs under control of the CMV promoter were previously described (Doumazane et al., 2011). For increased expression, the CMV promoter was exchanged with the synapsin-1 promoter (gift from B. Bettler) for CLIP-mGlu2 and SNAP-mGlu4. Homo- and heterodimers were expressed by co-transfection with CLIP-mGlu2 (100 ng/well) + SNAP-mGlu4 (200 ng/well), CLIP-mGlu2 (200 ng/well) + SNAP-mGlu1a (100 ng/well), or SNAP-mGlu2 (300 ng/well). For TR-FRET microscopy, 16–17 DIV neurons were labeled with 100 nM SNAP-Lumi4-Tb + 1000 nM CLIP-Red (heterodimers) or 100 nM SNAP-Lumi4-Tb + 500 nM SNAP-Red for 1 hr at 37°C in imaging buffer (10 mM Hepes pH 7.4, 127 mM NaCl, 2.8 mM KCl, 1.1 mM MgCl₂, 1.15 mM CaCl₂, 10 mM glucose) supplemented with 1% BSA followed by a wash in imaging buffer with 1% BSA and three washes in imaging buffer. Cells were imaged in imaging buffer. Images were acquired with a homebuilt TR-FRET microscope (Faklaris et al., 2015). Briefly, the donor was excited with a 349 nm Nd:YLF pulsed laser at 300 Hz with ~68 µJ/pulse followed by collection of either the donor signal using a 550/32 nm bandpass filter or the TR-FRET signal using a 700/75 nm bandpass filter. In both cases, images were acquired with 10 µs delay between excitation and collection of emission, 3 ms acquisition time and 4000 acquisitions. The acceptor was excited with a mercury lamp using a 620/60 nm bandpass filter, and the emission was collected for 300 ms with a 700/75 nm bandpass filter. Time-gated images were shading corrected by dividing the raw image with a background image using ImageJ version 1.51 f (Schneider et al., 2012). Correction for donor bleedthrough (6%) and generation of NFRET images (TR-FRET/(donor × acceptor)^0.5) was done with the PixFRET plugin to ImageJ (Feige et al., 2005). Acceptor bleedthrough and direct acceptor excitation was not detected. For quantification, all non-zero pixels in the NFRET image were selected, pixels not belonging to the cell removed and the modal NFRET value and mean donor signal were measured for this selection.

Electrophysiological recordings

Request a detailed protocol

Acute slices were prepared from adult (P21-P30) control or mGlu4-KO mice following a protocol approved by the European Communities Council Directive and the French low for care and use of experimental animals. Mice were decapitated, and brains quickly removed and chilled in cold artificial cerebro-spinal fluid (ACSF) containing 125 mM NaCl, 2.5 mM KCl, 1.25 mM NaH₂PO₄, 25 mM NaHCO₃, 1 mM MgCl₂, 2 mM CaCl₂, and 10 mM glucose, pH 7.4, equilibrated with 95% O₂ and 5% CO₂. Parasagittal acute 400-µm-thick slices were prepared with a vibratome (7000 smz, Campden Instruments LTD, England) in ice-cold ACSF. Sections were kept at room temperature for at least 1 hr before recording. Slices were transferred to a submersion recording chamber, maintained at 30°C and perfused with oxygenated ACSF at a rate of one chamber volume (1.5 ml) per minute.

fEPSPs were evoked at 0.033 Hz using bipolar stimuling electrode and recorded using glass micropipettes (3–5 MΩ) and filled with 3 M NaCl. Stimuling electrode was placed in the outer thirds of the molecular layer of the dentate gyrus for stimulation of the lateral perforant path (LPP). Correct positioning of electrodes was verified by application of paired-pulse at an interval of 100 ms induces paired-pulse facilitation in the LPP. The effect of paired-pulse stimulation was assessed and only those slices that displayed the correct facilitation in the LPP were used for this study. Input-output curves were generated for each slice, and the stimulation intensity was adjusted to 70% of the maximum response. Baseline fEPSPs were recorded for a minimum of 20 min before bath-application of different agonists or PAMs. Evoked responses were analyzed by measuring the slope of individual fEPSPs. The slopes from two sequential sweeps were averaged. All slopes were normalized to the average slope calculated during the pre-drug period (percentage of baseline). All data were analyzed offline using pClamp 9 (Molecular Devices) and are reported as the mean ± SEM. Statistical comparisons were made using two-tailed unpaired or paired Student’s t-tests. Differences were considered significant at p<0.05. The % inhibition was calculated by the difference of the slope between the baseline and the last 4 min of the drug application.

Curve fitting and data analysis

Request a detailed protocol

SAS/STAT 9.4 (SAS Institute, Cary, NC, USA) statistical package was used for parameter optimization and statistical analyses in mathematical modeling. Curve fitting was performed using nonlinear regression using GraphPad Prism software. p-values<0.05 were considered statistically significant using one-way ANOVA with Bonferroni post-hoc test.

Appendix 1

Analysis of the functional response of the mGlu2-4 system to LY354740 under fixed concentrations of L-AP4

LY354740 is an mGlu2 agonist

L-AP4 is an mGlu4 agonist

Note: LY354740 will be named LY35 for simplification

Concentration-effect results

Concentration-effect results from our experiments were selected for further analysis with a mathematical model.

The model

Model 1

We consider strict mGlu2-4 heterodimers. As we are considering agonists we will focus on the ECDs of the dimeric receptor. Two states either open (O) or closed (C) can be reached by each of the ECDs which can lead to heterodimers arranged as OO, OC or CC dimeric states. We will consider that if all heterodimers were in OO then a functional effect (F) value of 100 would be obtained whereas if they were all in CC an F value of 0 would be observed. Consistently, we assume that OC yields an intermediate value between 0 and 100.

Assumption

LY35 binds exclusively the mGlu2 protomer whereas L-AP4 binds exclusively the mGlu4 protomer.

With

\begin{matrix} K_{1} = \frac{[O^{2} O^{4}] [A]}{[O_{A}^{2} O^{4}]}; X_{1} = \frac{[C^{2} O^{4}]}{[O^{2} O^{4}]}; X_{2} = \frac{[O^{2} C^{4}]}{[O^{2} O^{4}]}; X_{3} = \frac{[C^{2} C^{4}]}{[C^{2} O^{4}]}; X_{4} = \frac{[C^{2} C^{4}]}{[O^{2} C^{4}]}; \\ Y_{1} = \frac{[C_{A}^{2} O^{4}]}{[O_{A}^{2} O^{4}]}; Y_{2} = \frac{[O_{A}^{2} C^{4}]}{[O_{A}^{2} O^{4}]}; Y_{3} = \frac{[C_{A}^{2} C^{4}]}{[C_{A}^{2} O^{4}]}; Y_{4} = \frac{[C_{A}^{2} C^{4}]}{[O_{A}^{2} C^{4}]} \end{matrix}

We define the functional response F as

F (%) = \frac{100 (O^{2} O^{4} + O_{A}^{2} O^{4} + f (C^{2} O^{4} + O^{2} C^{4} + C_{A}^{2} O^{4} + O_{A}^{2} C^{4}))}{[R_{T}]}

With 0 < f < 1.

We consider that OO states produce 100% F, CC states produce 0% F and OC states produce 0 < F(%)<100.

F (%) = \frac{100 (1 + \frac{[A]}{K_{1}} + f (X_{1} + X_{2} + (Y_{1} + Y_{2}) \frac{[A]}{K_{1}}))}{1 + X_{1} + X_{2} + X_{1} X_{3} + (1 + Y_{1} + Y_{2} + Y_{1} Y_{3}) \frac{[A]}{K_{1}}}

Equation 2 can be rearranged as the empirical Equation 3

F (%) = 100 \frac{a + b [A]}{c + d [A]}

With

\begin{array}{ll} a = & 1 + f (X_{1} + X_{2}) \\ b = & \frac{1}{K_{1}} (1 + f (Y_{1} + Y_{2})) \\ c = & 1 + X_{1} + X_{2} + X_{1} X_{3} \\ d = & \frac{1}{K_{1}} (1 + Y_{1} + Y_{2} + Y_{1} Y_{3}) \end{array}

If we divide the numerator and denominator of Equation 3 by d we have

F (%) = 100 \frac{a_{1} + a_{2} [A]}{a_{3} + [A]}

It can be shown that Equation 4 can be written as the typical Hill equation with a Hill coefficient of one (Equation 5).

F (%) = B o t t o m + \frac{T o p - B o t t o m}{1 + 10^{x - x_{50}}}

With $x = l o g [A], B o t t o m = 100 a_{2}, T o p = 100 \frac{a_{1}}{a_{3}}, x_{50} = l o g a_{3}$

If we retake the mechanistic constants that define a₁, a₂, and a₃ we can see that

The basal response (when [A] = 0) is defined by $T o p = 100 \frac{a_{1}}{a_{3}} = 100 \frac{a}{c} = 100 \frac{1 + f (X_{1} + X_{2})}{1 + X_{1} + X_{2} + X_{1} X_{3}}$ . Consistently with basal definition there is no constant related with the agonist A.
The minimum response, that is the asymptotic response as [A] increases, is defined by $Bottom = {100a}_{2} = 100 \frac{b}{d} = 100 \frac{1 + f (Y_{1} + Y_{2})}{1 + Y_{1} + Y_{2} + Y_{1} Y_{3}}$ . Bottom determines the efficacy of the ligand. . Consequently, the dissociation constant for binding is not present. A full agonist mGlu2 in the heterodimeric context would be one with a high Y₃, which leads to the formation of C²_AC⁴, that is both protomers are closed. Obviously, if we perform the concentration-response curve of the mGlu2 agonist in the presence of an mGlu4 agonist the closing of the mGlu4 subunit is facilitated, which affects both Top and Bottom.
The location of the curve along the X=log[A] axis is defined by $x_{50} = \log a_{3} = \log \frac{c}{d} = \log \frac{1 + X_{1} + X_{2} + X_{1} X_{3}}{\frac{1}{K_{1}} (1 + Y_{1} + Y_{2} + Y_{1} Y_{3})}$ . Consistently with potency definition, values related with efficacy (Y constants) and affinity (K₁constant) are present.

Finally, from the slope parameter point we conclude that the proposed mechanistic model with LY35 binding exclusively to mGlu2 and L-AP4 binding exclusively to mGlu4 produce Hill curves with Hill coefficients of one.

Note

The model can be used for the function of an agonist mGlu4 in the presence of fixed concentrations of an mGlu2 agonist.

Data analysis

Experimental data curves were fitted with Hill equations with nH = 1 and with nH allowed to be different from 1. To assess whether nH is statistically different from one different tests can be done.

Data are fitted with the Hill equation $F(%) = Bottom + \frac{Top - Bottom}{1 + 10^{nH(x - x_{50})}}$

Appendix 1—table 1

Hill equation parameters resulting of fitting curve data in Appendix 1—figure 1.

https://doi.org/10.7554/eLife.25233.019

L-AP4 conc (n)	Top Mean ± SEM	Bottom Mean ± SEM	x₅₀ Mean ± SEM	nH Mean ± SEM
0 (3)	102.72 ± 6.89	39.40 ± 2.63	−7.14 ± 0.11	0.77 ± 0.13
−10 (3)	100.29 ± 1.07	47.56 ± 1.29	−7.40 ± 0.09	0.65 ± 0.04
−9 (3)	104.33 ± 2.04	39.05 ± 1.57	−7.10 ± 0.14	0.68 ± 0.10
−8 (3)	88.64 ± 2.48	29.41 ± 1.71	−7.82 ± 0.11	0.67 ± 0.04
−7 (3)	82.23 ± 2.52	29.34 ± 0.62	−8.79 ± 0.17	0.93 ± 0.09
−6.5 (2)	61.65 ± 3.27	24.35 ± 0.74	−8.21 ± 0.05	1.40 ± 0.03
−6 (2)	54.20 ± 0.48	27.81 ± 3.12	−8.38 ± 0.23	1.32 ± 0.11
−5 (2)	46.14 ± 3.95	20.47 ± 3.17	−8.44 ± 0.20	1.21 ± 0.48

Appendix 1—figure 1

Download asset Open asset

Concentration-effect curves of LY35 in the presence of L-AP4.
https://doi.org/10.7554/eLife.25233.020

Addition of mGlu4 agonist displaces the curves downwards because closing of the mGlu4 protomer decreases F. Interestingly an apparent influence on the slope of the curves is observed: The Hill coefficient is less than one when L.AP4 is absent or at low concentration and increases to one at higher L-AP4 concentration.

Statistical analysis of the Hill coefficient

We analyze whether nH is statistically different from one in two ways.

1.We fit each of the curves with equation $F(%) = Bottom + \frac{Top - Bottom}{1 + 10^{nH(x - x_{50})}}$ and calculate the confidence interval at 95% for m parameter.

Appendix 1—table 2

Hill equation parameters resulting of fitting curve data in Appendix 1—figure 1.

https://doi.org/10.7554/eLife.25233.021

L-AP4 conc (n)	nH Mean ± SEM	Confidence interval of nH parameter (95%)
0 (3)	0.77 ± 0.13	(0.21, 1.33)
−10 (3)	0.65 ± 0.04	(0.46, 0.83)
−9 (3)	0.68 ± 0.10	(0.26, 1.09)
−8 (3)	0.67 ± 0.04	(0.49, 0.85)
−7 (3)	0.93 ± 0.09	(0.56, 1.30)
−6.5 (2)	1.40 ± 0.03	(1.04, 1.75)
−6 (2)	1.32 ± 0.11	(−0.12, 2.76)
−5 (2)	1.21 ± 0.48	(−4.94, 7.37)

Curves with [L-AP4]≤10⁻⁷ present nH <1 in average with some of them (10⁻¹⁰, 10⁻⁸) reaching statistical significance (the confidence interval of the nH parameter is below 1) and one of them (10⁻⁹) very close to be statistically significant.

2.We fit the collection of curves (3 or 2) for each L-AP4 concentration with two equations $F(%) = Bottom + \frac{Top - Bottom}{1 + 10^{nH(x - x_{50})}}$ and $F (%) = B o t t o m + \frac{T o p - B o t t o m}{1 + 10^{x - x_{50}}}$ and analyze with an F-test of the sum of squares errors whether the model including the nH parameter provides a better fit than that in which nH is not present.

Appendix 1—table 3

Statistical comparison of goodness of fit including the slope parameter (nH) or not to curves displayed in Appendix 1—figure 1.

https://doi.org/10.7554/eLife.25233.022

[LAP4]	SS1	df1	SS2	df2	F-value	p-value
0	1430.9	21	1330.1	20	1.515675513	0.232557
10**(−10)	541.3	17	415.4	16	4.849301878	0.042671
10**(−9)	852.8	21	727.1	20	3.457571173	0.077738
10**(−8)	560.5	21	455.6	20	4.604916594	0.044334
10**(−7)	273.4	21	273.1	20	0.021969974	0.883651
10**(−6.5)	220.5	12	207.8	11	0.672281039	0.429666
10**(−6)	206.5	13	206.5	12	0	1
10**(−5)	252.3	13	252	12	0.014285714	0.906839

Results are consistent with the previous analysis, the model including the slope parameter improves significantly the fitting for [LAP4] = 10⁻¹⁰ and 10⁻⁸ and close to significance for [LAP4] = 10⁻⁹.

Comments

The mechanistic model depicted in Appendix 1—figure 2 yields an empirical Hill equation with a Hill coefficient of 1. Experimental data suggest that the binding of LY35 to the heterodimer produces curves with Hill coefficient lower than one at low [L-AP4] and curves with Hill coefficient not different from one at high [L-AP4]. Thus, apparently, there is a contradiction between the mechanistic model and those results with the slope parameter lower than one.

Appendix 1—figure 2

Download asset Open asset

Model 1.

A heterodimeric mGlu2/4 model in which an mGlu2 agonist binds exclusively the mGlu2 protomer. The binding of an mGlu4 ligand to mGlu4 alters the constants of the model.

https://doi.org/10.7554/eLife.25233.023

A slope parameter lower than one could be explained assuming that LY35 binds at both mGlu2 and mGlu4 subunits in the heterodimer with crosstalk between them. Addition of mGlu4 agonist L-AP4 precludes the binding of LY35 to mGlu4 subunit and converts the heterodimeric receptor in a monomeric receptor for LY35.

Extending the model

Model 2

To account for concentration-effect curves with a Hill coefficient different from one the model displayed in Appendix 1—figure 2 was extended by allowing the possibility that LY35 could bind the mGlu4 protomer in addition to the mGlu2 one (Appendix 1—figure 3).

Appendix 1—figure 3

Download asset Open asset

Model 2.

A heterodimeric mGlu2/4 model in which an mGlu2 agonist binds both the mGlu2 and the mGlu4 protomers. The additional binding of an mGlu4 ligand alters the constants of the model.

https://doi.org/10.7554/eLife.25233.024

With

\begin{matrix} K_{1} = \frac{[O^{2} O^{4}] [A]}{[O_{A}^{2} O^{4}]}; K_{2} = \frac{[O^{2} O^{4}] [A]}{[O^{2} O_{A}^{4}]}; K_{3} = \frac{[O_{A}^{2} O^{4}] [A]}{[O_{A}^{2} O_{A}^{4}]}; K_{4} = \frac{[O^{2} O_{A}^{4}] [A]}{[O_{A}^{2} O_{A}^{4}]}; \\ X_{1} = \frac{[C^{2} O^{4}]}{[O^{2} O^{4}]}; X_{2} = \frac{[O^{2} C^{4}]}{[O^{2} O^{4}]}; X_{3} = \frac{[C^{2} C^{4}]}{[C^{2} O^{4}]}; X_{4} = \frac{[C^{2} C^{4}]}{[O^{2} C^{4}]}; \\ Y_{1} = \frac{[C_{A}^{2} O^{4}]}{[O_{A}^{2} O^{4}]}; Y_{2} = \frac{[O_{A}^{2} C^{4}]}{[O_{A}^{2} O^{4}]}; Y_{3} = \frac{[C_{A}^{2} C^{4}]}{[C_{A}^{2} O^{4}]}; Y_{4} = \frac{[C_{A}^{2} C^{4}]}{[O_{A}^{2} C^{4}]}; \\ Y_{5} = \frac{[C^{2} O_{A}^{4}]}{[O^{2} O_{A}^{4}]}; Y_{6} = \frac{[O^{2} C_{A}^{4}]}{[O^{2} O_{A}^{4}]}; Y_{7} = \frac{[O^{2} C_{A}^{4}]}{[O^{2} O_{A}^{4}]}; Y_{8} = \frac{[O^{2} C_{A}^{4}]}{[O^{2} O_{A}^{4}]}; \\ Y_{9} = \frac{[C_{A}^{2} O_{A}^{4}]}{[O_{A}^{2} O_{A}^{4}]}; Y_{10} = \frac{[C_{A}^{2} O_{A}^{4}]}{[O_{A}^{2} O_{A}^{4}]}; Y_{11} = \frac{[C_{A}^{2} O_{A}^{4}]}{[O_{A}^{2} O_{A}^{4}]}; Y_{12} = \frac{[C_{A}^{2} O_{A}^{4}]}{[O_{A}^{2} O_{A}^{4}]}; \end{matrix}

We define the functional response F as

F(%) = \frac{100 (O^{2} O^{4} + O^{2}_{A} O^{4} + O^{2} O^{4}_{A} + O^{2}_{A} O^{4}_{A} + f (\begin{matrix} C^{2} O^{4} + O^{2} C^{4} + C^{2}_{A} O^{4} + O^{2}_{A} C^{4} \\ + C^{2} O^{4}_{A} + O^{2} C^{4}_{A} + C^{2}_{A} O^{4}_{A} + O^{2}_{A} C^{4}_{A} + \end{matrix}))}{[R_{T}]}

With 0 < f < 1.

We consider that OO states produce 100% F, CC states produce 0% F and OC states produce 0 < F(%)<100.

F(%) = \frac{100 (\begin{matrix} 1 + f (X_{1} + X_{2}) + [A] (\frac{1}{K_{1}} + \frac{1}{K_{2}} + f (\frac{Y_{1} + Y_{2}}{K_{1}} + \frac{Y_{5} + Y_{6}}{K_{2}})) \\ + {[A]}^{2} \frac{1}{K_{1} K_{3}} (1 + {fY}_{9} (1 + Y_{10})) \end{matrix})}{\begin{matrix} 1 + X_{1} + X_{2} + X_{1} X_{3} + [A] (\frac{1}{K_{1}} + \frac{1}{K_{2}} + \frac{Y_{1} + Y_{2}}{K_{1}} + \frac{Y_{5} + Y_{6}}{K_{2}}) \\ + {[A]}^{2} \frac{1}{K_{1} K_{3}} (1 + Y_{9} (1 + Y_{10})) \end{matrix}}

Equation 7 can be rearranged as the empirical Equation 8

F(%) = 100 \frac{c_{1} + c_{2} [A] + c_{3} {[A]}^{2}}{c_{4} + c_{5} [A] + c_{6} {[A]}^{2}}

With

\begin{array}{ll} c_{1} = & 1 + f (X_{1} + X_{2}) \\ c_{2} = & \frac{1}{K_{1}} + \frac{1}{K_{2}} + f (\frac{Y_{1} + Y_{2}}{K_{1}} + \frac{Y_{5} + Y_{6}}{K_{2}}) \\ c_{3} = & \frac{1}{K_{1} K_{3}} (1 + f (Y_{9} + Y_{10})) \\ c_{4} = & 1 + X_{1} + X_{2} + X_{1} X_{3} \\ c_{5} = & \frac{1}{K_{1}} + \frac{1}{K_{2}} + \frac{Y_{1} + Y_{2} + Y_{1} Y_{3}}{K_{1}} + \frac{Y_{5} + Y_{6} + Y_{5} Y_{7}}{K_{2}} \\ c_{6} = & \frac{1}{K_{1} K_{3}} (1 + Y_{9} + Y_{10} + Y_{9} Y_{11}) \end{array}

The empirical ci coefficients reflect, when expressed as combination of mechanistic constants, what we would expect from the comparison between Equation 8 and Appendix 1—figure 3. That is, c1 and c4 include equilibrium constants related with free receptor species only. Analogously, c2 and c5 include equilibrium constants related with receptor species with only one bound agonist whereas c3 and c6 include equilibrium constants related with receptor species with two bound agonists.

Of note, the relationships between the empirical coefficients determine the shape of the concentration-effect curves (see below for shape quantification).

If we divide the numerator and denominator of Equation 8 by c₆ we have Equation 9.

F(\% ) = 100 \frac{a_{1} + a_{2} [A] + a_{3} {[A]}^{2}}{a_{4} + a_{5} [A] + {[A]}^{2}}

With a_i=c_i/c₆ for i = 1 to 6.

Equation 9 is an empirical equation for a ligand that binds to two receptor sites. This equation has been previously derived from some mechanistic models involving two receptor binding sites (see ^1,2 as examples and reviews ^3,4).

The shape of the concentration-effect curve determined by Equation 9 can be quantitatively characterized by the following geometric determinants.

The basal response (when [A] = 0) is defined by
$Top = 100 \frac{a_{1}}{a_{4}} = 100 \frac{c_{1}}{c_{4}} = 100 \frac{1 + f (X_{1} + X_{2})}{1 + X_{1} + X_{2} + X_{1} X_{3}}$
Consistently with basal definition there is no constant related with the agonist A.

Because 0<f < 1, Top will be always lower than 100. Furthermore, as higher is X3 (the propensity to form CC states) lower is Top.

The minimum response, that is the asymptotic response as [A] increases, is defined by
$Bottom = {100a}_{3} = 100 \frac{c_{3}}{c_{6}} = 100 \frac{1 + f (Y_{9} + Y_{10})}{1 + Y_{9} + Y_{10} + Y_{9} Y_{11}}$
Bottom determines the efficacy of the ligand. Considering the mechanistic constants included in Bottom definition, it follows that a3 <1. A ligand is a full agonist if a₃< < 1 and a partial agonist if a₃ <1.

In agreement with Bottom defined as efficacy, the dissociation constants for binding are not present in its mechanistic expression. A full agonist mGlu2 in the heterodimeric context would be one with a high Y₁₁, which leads to the formation of C²_AC⁴_A, that is, both protomers are closed. Obviously, if we perform the concentration-response curve of the mGlu2 agonist in the presence of an mGlu4 agonist the closing of the mGlu4 subunit is facilitated, which affects both Top and Bottom.

The location of the curve along the X=log[A] axis (X₅₀=log[A₅₀] or mid-point) defines the potency of the ligand and is defined by
$X_{50} = log (\frac{- b \pm \sqrt{b^{2} - 4ac}}{2a})$
Consistently with potency definition, values related with efficacy (Y constants) and affinity (K_iconstants) are present.

Where

a = a_{1} - a_{3} a_{4}; a = a_{3} a_{4} a_{5} - 2 a_{2} a_{4} + a_{1} a_{5}; a n d c = - a_{4} (a_{1} (a_{1} - a_{3} a_{4}))

•Quantification of cooperativity by the calculation of the Hill coefficient can be done by making use of the definition of the Hill coefficient at the mid-point ( $n_{H_{50}}$ ) for a given y(x) function.⁵
$n_{H_{50}} = \frac{4 {(\frac{dy}{dx})}_{x_{50}}}{aln10}$
With y=F(%), x = log[A]; a, the Bottom; ln, the natural logarithm; and d/dx, the derivative operator as expressed in Equation 14.
$\frac{dy}{dx} = \frac{100 (- (a_{1} + a_{2} 10^{x} + a_{3} 10^{2x}) (a_{5} 10^{x} + 2 \cdot 10^{2x}) + ({2a}_{3} 10^{2x} + a_{2} 10^{x}) (a_{4} + a_{5} 10^{x} + 10^{2x})) ln10}{{(a_{4} + a_{5} 10^{x} + 10^{2x})}^{2}}$

The value of the Hill coefficient as obtained from Equation 13, with empirical coefficients (a₁ to a₅), which in turn are defined in terms of mechanistic equilibrium constants, may provide a mechanistic interpretation to the Hill coefficient obtained by fitting with the empirical Hill equation including the slope parameter.

The ratio c₁/c₄ determines the Top asymptote (basal response) and the ratio c₃/c₆ determines the Bottom asymptote (efficacy). The ratio c₂/c₅ determines the sensitivity of the measured effect to agonist concentration. Considering the mechanistic constants included in c₂ and c₅ definition, it follows that c₂ <c₅. The induction constants that appear in c2 and c5 expressions are Y₁, Y₂, Y₃, Y₅, Y₆ and Y₇, which are those constants affecting receptor species with only one molecule of mGlu2 agonist present. We see that Y₃ and Y₇ are present in c₅ but not in c_2; then, the values of these constants may modulate the c₂/c₅ ratio. Because of the closure of the thermodynamic cycles included in the model displayed in Fig. S3, Y₄ and Y₈ ( $Y_{4} = \frac{Y_{1}}{Y_{2}} Y_{3}, Y_{8} = \frac{Y_{5}}{Y_{6}} Y_{7}$ ) can be used instead of Y₃ and Y₇, respectively. The pair (Y₃, Y₇) or the pair (Y₄, Y₈) measure how the closure of one protomer favors the closure of the other thus it can be considered as a measure of functional cooperativity. Thus, we can conclude that the functional cooperativity between the two protomers affects the sensitivity of the measured effect and be the cause of some of the flat curves observed.

Equation 9 contains five parameters and is difficult to fit to curves that do not display a clear biphasic shape. However it may be used for modeling different pharmacological conditions by assigning particular values to the parameters.

Simulation of pharmacological conditions under the mechanistic models

Appendix 1—figure 4 illustrates how Model two can explain the flat curve observed for LY35 mGlu2 agonist. We assume that the closed-closed state is not achieved by proposing the induction constants Y₃ = Y₇ = Y₁₁=10⁻⁶. These constants make LY35 to behave as a partial agonist with a bottom value of 51%. We assume that there is negative binding cooperativity and the ligand binds the mGlu4 protomer after occupying first the mGlu2 binding site (K₁ = 10⁻⁶; K₂ = 10³; K₃ = 10⁻³). This leads to a curve with two components, one related with the binding of the first molecule to the heterodimer and another one related with the binding of the second molecule. The induction constant for the closure of the mGlu2 subunit is greater in the doubly- than in the singly-bound heterodimer (Y₉ >Y₁). A f-value of 0.5 was used for the functional closed-open state.

Appendix 1—figure 4

Download asset Open asset

Theoretical concentration-effect curve for particular values of the mechanistic constants included in Model 2.

The theoretical concentration-effect data are the following.

https://doi.org/10.7554/eLife.25233.025

Fitting with the Hill equation $F(%) = Bottom + \frac{Top - Bottom}{1 + 10^{nH(x - x_{50})}}$ yielded the following parameters

Appendix 1—table 4

Data extracted from concentration-effect curve of Appendix 1—figure 4.

https://doi.org/10.7554/eLife.25233.026

Log[LY354740]	Effect (%)
−12.00	100.00
−11.00	100.00
−10.00	99.99
−9.00	99.90
−8.00	99.03
−7.00	92.28
−6.00	74.40
−5.00	63.50
−4.00	54.23
−3.00	50.96

Appendix 1—figure 5

Download asset Open asset

Curve data included in Appendix 1—table 4 (solid points) and the theoretical curve by using the Hill equation parameters of Appendix 1—table 5 (curve line).
https://doi.org/10.7554/eLife.25233.027

The Hill coefficient nH is ~0.6 in agreement with experimental data.

The following graph includes the theoretical data from Model two and the curve produced by using the Hill equation fitted parameters.

Appendix—table 5

Parameter values by fitting curve data in Appendix 1—table 4 with the Hill equation.

https://doi.org/10.7554/eLife.25233.028

Parameter	Estimate	Approx std error	Approximate 95% Confidence Limits
Bottom	50.9874	1.5333	47.2355	54.7393
Top	100.5	0.7699	98.6649	102.4
x50	−5.9487	0.0920	−6.1739	−5.7234
nH	0.5993	0.0669	0.4356	0.7631

References

1. Beurrier C
2. Lopez S
3. Révy D
4. Selvam C
5. Goudet C
6. Lhérondel M
7. Gubellini P
8. Kerkerian-LeGoff L
9. Acher F
10. Pin JP
11. Amalric M
(2009) Electrophysiological and behavioral evidence that modulation of metabotropic glutamate receptor 4 with a new agonist reverses experimental parkinsonism
The FASEB Journal 23:3619–3628.

https://doi.org/10.1096/fj.09-131789
- PubMed
- Google Scholar
1. Bouvier M
2. Hébert TE
(2014) CrossTalk proposal: weighing the evidence for class A GPCR dimers, the evidence favours dimers
The Journal of Physiology 592:2439–2441.

https://doi.org/10.1113/jphysiol.2014.272252
- PubMed
- Google Scholar
1. Bradley SR
2. Standaert DG
3. Rhodes KJ
4. Rees HD
5. Testa CM
6. Levey AI
7. Conn PJ
(1999) Immunohistochemical localization of subtype 4a metabotropic glutamate receptors in the rat and mouse basal ganglia
The Journal of Comparative Neurology 407:33–46.

https://doi.org/10.1002/(SICI)1096-9861(19990428)407:1<33::AID-CNE3>3.0.CO;2-G
- PubMed
- Google Scholar
1. Brock C
2. Oueslati N
3. Soler S
4. Boudier L
5. Rondard P
6. Pin JP
(2007) Activation of a dimeric metabotropic glutamate receptor by intersubunit rearrangement
Journal of Biological Chemistry 282:33000–33008.

https://doi.org/10.1074/jbc.M702542200
- PubMed
- Google Scholar
1. Calebiro D
2. Rieken F
3. Wagner J
4. Sungkaworn T
5. Zabel U
6. Borzi A
7. Cocucci E
8. Zürn A
9. Lohse MJ
(2013) Single-molecule analysis of fluorescently labeled G-protein-coupled receptors reveals complexes with distinct dynamics and organization
PNAS 110:743–748.

https://doi.org/10.1073/pnas.1205798110
- PubMed
- Google Scholar
1. Conn PJ
2. Pin JP
(1997) Pharmacology and functions of metabotropic glutamate receptors
Annual Review of Pharmacology and Toxicology 37:205–237.

https://doi.org/10.1146/annurev.pharmtox.37.1.205
- PubMed
- Google Scholar
(2005) Metabotropic glutamate receptors in the basal ganglia motor circuit
Nature Reviews Neuroscience 6:787–798.

https://doi.org/10.1038/nrn1763
- PubMed
- Google Scholar
(2014) Opportunities and challenges in the discovery of allosteric modulators of GPCRs for treating CNS disorders
Nature Reviews Drug Discovery 13:692–708.

https://doi.org/10.1038/nrd4308
- PubMed
- Google Scholar
(2002) Distribution and synaptic localisation of the metabotropic glutamate receptor 4 (mGluR4) in the rodent CNS
Neuroscience 110:403–420.

https://doi.org/10.1016/S0306-4522(01)00591-7
- PubMed
- Google Scholar
1. Doumazane E
2. Scholler P
3. Zwier JM
4. Trinquet E
5. Rondard P
6. Pin JP
(2011) A new approach to analyze cell surface protein complexes reveals specific heterodimeric metabotropic glutamate receptors
The FASEB Journal 25:66–77.

https://doi.org/10.1096/fj.10-163147
- PubMed
- Google Scholar
1. Doumazane E
2. Scholler P
3. Fabre L
4. Zwier JM
5. Trinquet E
6. Pin JP
7. Rondard P
(2013) Illuminating the activation mechanisms and allosteric properties of metabotropic glutamate receptors
PNAS 110:E1416–E1425.

https://doi.org/10.1073/pnas.1215615110
- PubMed
- Google Scholar
1. Faklaris O
2. Cottet M
3. Falco A
4. Villier B
5. Laget M
6. Zwier JM
7. Trinquet E
8. Mouillac B
9. Pin JP
10. Durroux T
(2015) Multicolor time-resolved förster resonance energy transfer microscopy reveals the impact of GPCR oligomerization on internalization processes
The FASEB Journal 29:2235–2246.

https://doi.org/10.1096/fj.14-260059
- PubMed
- Google Scholar
1. Feige JN
2. Sage D
3. Wahli W
4. Desvergne B
5. Gelman L
(2005) PixFRET, an ImageJ plug-in for FRET calculation that can accommodate variations in spectral bleed-throughs
Microscopy Research and Technique 68:51–58.

https://doi.org/10.1002/jemt.20215
- PubMed
- Google Scholar
1. Ferraguti F
2. Klausberger T
3. Cobden P
4. Baude A
5. Roberts JD
6. Szucs P
7. Kinoshita A
8. Shigemoto R
9. Somogyi P
10. Dalezios Y
(2005) Metabotropic glutamate receptor 8-expressing nerve terminals target subsets of GABAergic neurons in the Hippocampus
Journal of Neuroscience 25:10520–10536.

https://doi.org/10.1523/JNEUROSCI.2547-05.2005
- PubMed
- Google Scholar
1. Ferraguti F
2. Shigemoto R
(2006) Metabotropic glutamate receptors
Cell and Tissue Research 326:483–504.

https://doi.org/10.1007/s00441-006-0266-5
- PubMed
- Google Scholar
1. Ferré S
2. Casadó V
3. Devi LA
4. Filizola M
5. Jockers R
6. Lohse MJ
7. Milligan G
8. Pin JP
9. Guitart X
(2014) G protein-coupled receptor oligomerization revisited: functional and pharmacological perspectives
Pharmacological Reviews 66:413–434.

https://doi.org/10.1124/pr.113.008052
- PubMed
- Google Scholar
(2008) An engineered protein tag for multiprotein labeling in living cells
Chemistry & Biology 15:128–136.

https://doi.org/10.1016/j.chembiol.2008.01.007
- PubMed
- Google Scholar
1. Giraldo J
(2003) Empirical models and Hill coefficients
Trends in Pharmacological Sciences 24:63–65.

https://doi.org/10.1016/S0165-6147(02)00048-2
- PubMed
- Google Scholar
1. Giraldo J
(2013) Modeling cooperativity effects in dimeric G protein-coupled receptors
Progress in Molecular Biology and Translational Science 115:349–373.

https://doi.org/10.1016/B978-0-12-394587-7.00008-7
- PubMed
- Google Scholar
1. Goudet C
2. Vilar B
3. Courtiol T
4. Deltheil T
5. Bessiron T
6. Brabet I
7. Oueslati N
8. Rigault D
9. Bertrand HO
10. McLean H
11. Daniel H
12. Amalric M
13. Acher F
14. Pin JP
(2012) A novel selective metabotropic glutamate receptor 4 agonist reveals new possibilities for developing subtype selective ligands with therapeutic potential
The FASEB Journal 26:1682–1693.

https://doi.org/10.1096/fj.11-195941
- PubMed
- Google Scholar
1. Gu G
2. Lorrain DS
3. Wei H
4. Cole RL
5. Zhang X
6. Daggett LP
7. Schaffhauser HJ
8. Bristow LJ
9. Lechner SM
(2008) Distribution of metabotropic glutamate 2 and 3 receptors in the rat forebrain: implication in emotional responses and central disinhibition
Brain Research 1197:47–62.

https://doi.org/10.1016/j.brainres.2007.12.057
- PubMed
- Google Scholar
1. Hlavackova V
2. Goudet C
3. Kniazeff J
4. Zikova A
5. Maurel D
6. Vol C
7. Trojanova J
8. Prézeau L
9. Pin JP
10. Blahos J
(2005) Evidence for a single heptahelical domain being turned on upon activation of a dimeric GPCR
The EMBO Journal 24:499–509.

https://doi.org/10.1038/sj.emboj.7600557
- PubMed
- Google Scholar
1. Hlavackova V
2. Zabel U
3. Frankova D
4. Bätz J
5. Hoffmann C
6. Prezeau L
7. Pin JP
8. Blahos J
9. Lohse MJ
(2012) Sequential inter- and intrasubunit rearrangements during activation of dimeric metabotropic glutamate receptor 1
Science Signaling 5:ra59.

https://doi.org/10.1126/scisignal.2002720
- PubMed
- Google Scholar
1. Huang S
2. Cao J
3. Jiang M
4. Labesse G
5. Liu J
6. Pin JP
7. Rondard P
(2011) Interdomain movements in metabotropic glutamate receptor activation
PNAS 108:15480–15485.

https://doi.org/10.1073/pnas.1107775108
- PubMed
- Google Scholar
(2001)
Chemical induction of mGluR5- and protein synthesis--dependent long-term depression in hippocampal area CA1

Journal of Neurophysiology 86:321–325.
- PubMed
- Google Scholar
1. Iacovelli L
2. Molinaro G
3. Battaglia G
4. Motolese M
5. Di Menna L
6. Alfiero M
7. Blahos J
8. Matrisciano F
9. Corsi M
10. Corti C
11. Bruno V
12. De Blasi A
13. Nicoletti F
(2009) Regulation of group II metabotropic glutamate receptors by G protein-coupled receptor kinases: mglu2 receptors are resistant to homologous desensitization
Molecular Pharmacology 75:991–1003.

https://doi.org/10.1124/mol.108.052316
- PubMed
- Google Scholar
1. Johnson MP
2. Barda D
3. Britton TC
4. Emkey R
5. Hornback WJ
6. Jagdmann GE
7. McKinzie DL
8. Nisenbaum ES
9. Tizzano JP
10. Schoepp DD
(2005) Metabotropic glutamate 2 receptor potentiators: receptor modulation, frequency-dependent synaptic activity, and efficacy in preclinical anxiety and psychosis model(s)
Psychopharmacology 179:271–283.

https://doi.org/10.1007/s00213-004-2099-9
- PubMed
- Google Scholar
1. Juillerat A
2. Gronemeyer T
3. Keppler A
4. Gendreizig S
5. Pick H
6. Vogel H
7. Johnsson K
(2003) Directed evolution of O6-alkylguanine-DNA alkyltransferase for efficient labeling of fusion proteins with small molecules in vivo
Chemistry & Biology 10:313–317.

https://doi.org/10.1016/S1074-5521(03)00068-1
- PubMed
- Google Scholar
1. Kammermeier PJ
(2012) Functional and pharmacological characteristics of metabotropic glutamate receptors 2/4 heterodimers
Molecular Pharmacology 82:438–447.

https://doi.org/10.1124/mol.112.078501
- PubMed
- Google Scholar
1. Kniazeff J
2. Bessis AS
3. Maurel D
4. Ansanay H
5. Prézeau L
6. Pin JP
(2004) Closed state of both binding domains of homodimeric mGlu receptors is required for full activity
Nature Structural & Molecular Biology 11:706–713.

https://doi.org/10.1038/nsmb794
- PubMed
- Google Scholar
1. Kunishima N
2. Shimada Y
3. Tsuji Y
4. Sato T
5. Yamamoto M
6. Kumasaka T
7. Nakanishi S
8. Jingami H
9. Morikawa K
(2000) Structural basis of glutamate recognition by a dimeric metabotropic glutamate receptor
Nature 407:971–977.

https://doi.org/10.1038/35039564
- PubMed
- Google Scholar
1. Lambert NA
2. Javitch JA
(2014) CrossTalk opposing view: weighing the evidence for class A GPCR dimers, the jury is still out
The Journal of Physiology 592:2443–2445.

https://doi.org/10.1113/jphysiol.2014.272997
- PubMed
- Google Scholar
1. Levitz J
2. Habrian C
3. Bharill S
4. Fu Z
5. Vafabakhsh R
6. Isacoff EY
(2016) Mechanism of assembly and cooperativity of homomeric and heteromeric metabotropic glutamate receptors
Neuron 92:143–159.

https://doi.org/10.1016/j.neuron.2016.08.036
- PubMed
- Google Scholar
(1999) GABAB receptors - the first 7tm heterodimers
Trends in Pharmacological Sciences 20:396–399.

https://doi.org/10.1016/S0165-6147(99)01383-8
- PubMed
- Google Scholar
1. Maurel D
2. Comps-Agrar L
3. Brock C
4. Rives ML
5. Bourrier E
6. Ayoub MA
7. Bazin H
8. Tinel N
9. Durroux T
10. Prézeau L
11. Trinquet E
12. Pin JP
(2008) Cell-surface protein-protein interaction analysis with time-resolved FRET and snap-tag technologies: application to GPCR oligomerization
Nature Methods 5:561–567.

https://doi.org/10.1038/nmeth.1213
- PubMed
- Google Scholar
1. Neki A
2. Ohishi H
3. Kaneko T
4. Shigemoto R
5. Nakanishi S
6. Mizuno N
(1996) Pre- and postsynaptic localization of a metabotropic glutamate receptor, mGluR2, in the rat brain: an immunohistochemical study with a monoclonal antibody
Neuroscience Letters 202:197–200.

https://doi.org/10.1016/0304-3940(95)12248-6
- PubMed
- Google Scholar
1. Niswender CM
2. Conn PJ
(2010) Metabotropic glutamate receptors: physiology, pharmacology, and disease
Annual Review of Pharmacology and Toxicology 50:295–322.

https://doi.org/10.1146/annurev.pharmtox.011008.145533
- PubMed
- Google Scholar
1. Niswender CM
2. Jones CK
3. Lin X
4. Bubser M
5. Thompson Gray A
6. Blobaum AL
7. Engers DW
8. Rodriguez AL
9. Loch MT
10. Daniels JS
11. Lindsley CW
12. Hopkins CR
13. Javitch JA
14. Conn PJ
(2016) Development and Antiparkinsonian activity of VU0418506, a selective positive allosteric modulator of Metabotropic Glutamate receptor 4 homomers without activity at mGlu2/4 heteromers
ACS Chemical Neuroscience 7:1201–1211.

https://doi.org/10.1021/acschemneuro.6b00036
- PubMed
- Google Scholar
(1993) Distribution of the messenger RNA for a metabotropic glutamate receptor, mGluR2, in the central nervous system of the rat
Neuroscience 53:1009–1018.

https://doi.org/10.1016/0306-4522(93)90485-X
- PubMed
- Google Scholar
(2016) Group 1 metabotropic glutamate receptors 1 and 5 form a protein complex in mouse Hippocampus and cortex
PROTEOMICS 16:2698–2705.

https://doi.org/10.1002/pmic.201500400
- PubMed
- Google Scholar
(1996)
Impaired cerebellar synaptic plasticity and motor performance in mice lacking the mGluR4 subtype of metabotropic glutamate receptor

Journal of Neuroscience 16:6364–6373.
- PubMed
- Google Scholar
1. Pin JP
2. Neubig R
3. Bouvier M
4. Devi L
5. Filizola M
6. Javitch JA
7. Lohse MJ
8. Milligan G
9. Palczewski K
10. Parmentier M
11. Spedding M
(2007) International union of basic and clinical pharmacology. LXVII. recommendations for the recognition and nomenclature of G protein-coupled receptor heteromultimers
Pharmacological Reviews 59:5–13.

https://doi.org/10.1124/pr.59.1.5
- PubMed
- Google Scholar
1. Pin JP
2. Bettler B
(2016) Organization and functions of mGlu and GABAB receptor complexes
Nature 540:60–68.

https://doi.org/10.1038/nature20566
- PubMed
- Google Scholar
1. Pitsch J
2. Schoch S
3. Gueler N
4. Flor PJ
5. van der Putten H
6. Becker AJ
(2007) Functional role of mGluR1 and mGluR4 in pilocarpine-induced temporal lobe epilepsy
Neurobiology of Disease 26:623–633.

https://doi.org/10.1016/j.nbd.2007.03.003
- PubMed
- Google Scholar
1. Roche D
2. Gil D
3. Giraldo J
(2014) Mathematical modeling of G protein-coupled receptor function: what can we learn from empirical and mechanistic models?
Advances in Experimental Medicine and Biology 796:159–181.

https://doi.org/10.1007/978-94-007-7423-0_8
- PubMed
- Google Scholar
(1996) Metabotropic glutamate receptor 5 is a disulfide-linked dimer
Journal of Biological Chemistry 271:28612–28616.

https://doi.org/10.1074/jbc.271.45.28612
- PubMed
- Google Scholar
1. Rovira X
2. Roche D
3. Serra J
4. Kniazeff J
5. Pin JP
6. Giraldo J
(2008a) Modeling the binding and function of metabotropic glutamate receptors
Journal of Pharmacology and Experimental Therapeutics 325:443–456.

https://doi.org/10.1124/jpet.107.133967
- PubMed
- Google Scholar
(2010) The asymmetric/symmetric activation of GPCR dimers as a possible mechanistic rationale for multiple signalling pathways
Trends in Pharmacological Sciences 31:15–21.

https://doi.org/10.1016/j.tips.2009.10.008
- PubMed
- Google Scholar
(2012) NIH image to ImageJ: 25 years of image analysis
Nature Methods 9:671–675.

https://doi.org/10.1038/nmeth.2089
- PubMed
- Google Scholar
1. Scholler P
2. Moreno-Delgado D
3. Lecat-Guillet N
4. Doumazane E
5. Monnier C
6. Charrier-Savournin F
7. Fabre L
8. Chouvet C
9. Soldevila S
10. Lamarque L
11. Donsimoni G
12. Roux T
13. Zwier JM
14. Trinquet E
15. Rondard P
16. Pin JP
(2017) HTS-compatible FRET-based conformational sensors clarify membrane receptor activation
Nature Chemical Biology 13:372–380.

https://doi.org/10.1038/nchembio.2286
- PubMed
- Google Scholar
1. Shigemoto R
2. Kinoshita A
3. Wada E
4. Nomura S
5. Ohishi H
6. Takada M
7. Flor PJ
8. Neki A
9. Abe T
10. Nakanishi S
11. Mizuno N
(1997)
Differential presynaptic localization of metabotropic glutamate receptor subtypes in the rat Hippocampus

Journal of Neuroscience 17:7503–7522.
- PubMed
- Google Scholar
(2000) Dominant phenotypes produced by the HD mutation in STHdh(Q111) striatal cells
Human Molecular Genetics 9:2799–2809.

https://doi.org/10.1093/hmg/9.19.2799
- PubMed
- Google Scholar
1. Volk LJ
2. Daly CA
3. Huber KM
(2006) Differential roles for group 1 mGluR subtypes in induction and expression of chemically induced hippocampal long-term depression
Journal of Neurophysiology 95:2427–2438.

https://doi.org/10.1152/jn.00383.2005
- PubMed
- Google Scholar
1. Whorton MR
2. Bokoch MP
3. Rasmussen SG
4. Huang B
5. Zare RN
6. Kobilka B
7. Sunahara RK
(2007) A monomeric G protein-coupled receptor isolated in a high-density lipoprotein particle efficiently activates its G protein
PNAS 104:7682–7687.

https://doi.org/10.1073/pnas.0611448104
- PubMed
- Google Scholar
1. Wright RA
2. Johnson BG
3. Zhang C
4. Salhoff C
5. Kingston AE
6. Calligaro DO
7. Monn JA
8. Schoepp DD
9. Marek GJ
(2013) CNS distribution of metabotropic glutamate 2 and 3 receptors: transgenic mice and [³H]LY459477 autoradiography
Neuropharmacology 66:89–98.

https://doi.org/10.1016/j.neuropharm.2012.01.019
- PubMed
- Google Scholar
(2014) Selective actions of novel allosteric modulators reveal functional heteromers of metabotropic glutamate receptors in the CNS
Journal of Neuroscience 34:79–94.

https://doi.org/10.1523/JNEUROSCI.1129-13.2014
- PubMed
- Google Scholar
1. Zhao GQ
2. Zhang Y
3. Hoon MA
4. Chandrashekar J
5. Erlenbach I
6. Ryba NJ
7. Zuker CS
(2003) The receptors for mammalian sweet and umami taste
Cell 115:255–266.

https://doi.org/10.1016/S0092-8674(03)00844-4
- PubMed
- Google Scholar

Article and author information

Author details

David Moreno Delgado

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
DMD, Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing—original draft

Competing interests
No competing interests declared.

"This ORCID iD identifies the author of this article:" 0000-0003-0070-2566
Thor C Møller

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
TCM, Formal analysis, Investigation, Methodology, performed the TR-FRET microscopy experiments in cultured neurons, analyzed the results, prepared the figures and text related to these experiments

Competing interests
No competing interests declared.

"This ORCID iD identifies the author of this article:" 0000-0002-2256-577X
Jeanne Ster

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
JS, Conceptualization, Formal analysis, Investigation, Methodology, Writing—review and editing, performed the electrophysiological experiments and wrote the corresponding part of the ms

Competing interests
No competing interests declared.
Jesús Giraldo
1. Laboratory of Molecular Neuropharmacology and Bioinformatics, Institut de Neurociències and Unitat de Bioestadística, Universitat Autònoma de Barcelona, Bellaterra, Spain
2. Network Biomedical Research Center on Mental Health, Madrid, Spain
Contribution
JG, Conceptualization, Formal analysis, Funding acquisition, Investigation, supervised the mathematical modeling, and wrote the appendix

Competing interests
No competing interests declared.

"This ORCID iD identifies the author of this article:" 0000-0001-7082-4695
Damien Maurel

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
DM, Conceptualization, Supervision, Methodology, developped the XCelligence method in the laboratory, validating the assay with the proper controls, and trained D Moreno in setting up the assay for his experiments

Competing interests
No competing interests declared.
Xavier Rovira

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
XR, Investigation, Methodology

Competing interests
No competing interests declared.

"This ORCID iD identifies the author of this article:" 0000-0002-9764-9927
Pauline Scholler

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
PS, Supervision, Investigation, Methodology, Set up and validated the heterodimer selective sensor assay and trained D Moreno in using this assay. Performed the preliminary experiments at the origin of this project

Competing interests
No competing interests declared.
Jurriaan M Zwier

CisBio Bioassays, Codolet, France

Contribution
JMZ, Formal analysis, Methodology, Writing—review and editing, Essential role in the development of the TR-FRET microscope, and supervise the development of the sensor assay. Provided his knowledge in FRET

Competing interests
JMZ: Dr Zwier is working at CisBio Bioassays, a company selling HTRF tools

"This ORCID iD identifies the author of this article:" 0000-0001-9017-5856
Julie Perroy

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
JP, Conceptualization, Supervision, Writing—review and editing, suêrvised the electrophysiological experiments

Competing interests
No competing interests declared.
Thierry Durroux

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
TD, Supervision, Methodology, Developed the TR-FRET microscope and trained T Moeller in using it. Participated in the analysis of the TR-FRET microscopy data

Competing interests
No competing interests declared.
Eric Trinquet

CisBio Bioassays, Codolet, France

Contribution
ET, Conceptualization, Resources, Supervision, Methodology, Supervise the Cisbio team, and provided a lot of chemical tools essential to the project. Partcipated in the development of the TR-FRET sensor development

Competing interests
ET: Head of R&D at CisBio Bioassays
Laurent Prezeau

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
LP, Formal analysis, Funding acquisition, Visualization, Methodology, Supervise the arpege facility and initiated the Xcelligence technology in the laboratory. Raised the funding necessary for the acquisition of this technology

Competing interests
No competing interests declared.
Philippe Rondard

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
PR, Conceptualization, Formal analysis, Supervision, Writing—review and editing, supervise the preliminary experiments, and the TR-FRET sensor development

Competing interests
No competing interests declared.

"This ORCID iD identifies the author of this article:" 0000-0003-1134-2738
Jean-Philippe Pin

Institut de Génomique Fonctionnelle (IGF), CNRS, INSERM, Univ. Montpellier, Montpellier, France

Contribution
J-PP, Conceptualization, Supervision, Funding acquisition, Methodology, Writing—review and editing

For correspondence
jppin@igf.cnrs.fr

Competing interests
No competing interests declared.

"This ORCID iD identifies the author of this article:" 0000-0002-1423-345X

Funding

Seventh Framework Programme (627227)

Thor C Møller

European Commission (ERA-NET Neuron PCIN-2013-018-C03-02)

Jesús Giraldo

European Commission (ERA-NET Neuron SAF2014-58396-R)

Jesús Giraldo

Agence Nationale de la Recherche (ANR-09-BIOT-018)

Jurriaan M Zwier
Eric Trinquet
Philippe Rondard
Jean-Philippe Pin

Fondation pour la Recherche Médicale (DEQ20130326522)

Jean-Philippe Pin

CisBio Bioassays (Laboratoire Collaboratif)

Jean-Philippe Pin

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was made possible thanks to the Arpege facility from BioCampus (IGF, Montpellier, France). The authors wish to thank Drs A Llebaria and X Gomez (Barcelona, Spain), and F Acher (Paris, France) for providing VU0415374 and LSP4-2022, respectively; Dr Sílvia Ginés (University of Barcelona, Spain) for providing the STHdh^Q7 cells; and Dr Bernhard Bettler (Pharmacentrum, Basel, Switzerland) for the synapsin-1 promoter. The authors also thank Michaël Mathieu for the delta opioid receptor experiments. This work was supported by the Centre National de la Recherche Scientifique, the Institut National de la Recherche Médicale, the University of Montpellier, Cisbio Bioassays and by the Fondation pour la Recherche Médicale (grant no. DEQ20130326522) to JPP, and by the Agence Nationale de la Recherche (ANR-09-BIOT-018) to ET, JMZ, PR and JPP. JG was supported by the Ministerio de Economía y Competitividad (ERA-NET NEURON PCIN-2013–018 C03-02 and SAF2014-58396-R). This project has received funding from the European Union’s Seventh Framework Program for research, technological development and demonstration under grant agreement no 627227 (TCM)

Copyright

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.