Correlation of MCHR1 and MRAP2 co-expression
We developed a custom script to interrogate the co-expression correlation of MCHR1 and MRAP1 or MRAP2 in 1698 human or mice bulk RNA-seq samples of central nervous system (including 1000 human whole brain samples, 435 mice whole brain samples, 202 human hypothalamic samples and 61 mouse hypothalamic samples) from the Genotype-Tissue Expression (GTEx) project database and all RNA-seq and ChIP-seq sample and signature search (ARCHS4) database. As shown in Fig. 1A-B, a few number of cells co-expressed MCHR1 and MRAP1 in bulk RNA-seq datasets, which showed a low correlation (mean_cor = 0.2) between MCHR1 and MRAP1. While the number of co-expressing MCHR1 and MRAP2 is comparatively larger, which revealed high correlation (mean_cor = 0.6) between MCHR1 and MRAP2. In order to better explore the effect of MCHR1 on the phenotype of MRAP2 null mice and in view of the accuracy and specificity of single-cell RNA sequencing, we next analyzed the co-expression correlation of MCHR1 and MRAP2 in 28,320 cells of 20 mice from 4 published single-cell RNA-seq datasets of Gene Expression Omnibus (GEO) database. We found that MCHR1 and MRAP2 have lots of number of co-expressed cells in single-cell sequencing data as well (Fig. 1C). The average correlation between MCHR1 and MRAP2 in the 4 datasets is relatively large (mean_cor = 0.5).
MCHR1 and MRAP2 both play a critical role in energy homeostasis. Therefore, we analyzed the cell number ratio change of MCHR1 and MRAP2 under different metabolic condition to explore their metabolic relevance. The results showed that MCHR1 and MRAP2 had consistent cell ratio shift in HFD group compared to control, which were both significantly up-regulated (Fig. 1D). Furthermore, we found that MRAP1, MRAP2 and MCHR1 can be enriched in some energy metabolism pathways. Among them, MRAP1 and MRAP2 are both related to the binding of neuropeptides and GPCR signaling pathways. More importantly, GO(Gene Ontology) functional analysis exhibits that MCHR1 and MRAP2 both attributed to regulation of feeding behavior (Fig. 1E). These results all preliminarily indicate that MCHR1 and MRAP2 may have a synergistic effect on regulating energy metabolism in central nervous system.
Tissue-distribution and interaction of MCHR1 with MRAP1 or MRAP2
To characterize mRNA expression profiles of MCHR1, MRAP1 and MRAP2 in vivo, Our RT-PCR analysis was extended to 14 tissues from mice. MRAP1 and MRAP2 show universal expression in different tissues compared to MCHR1 (Fig. 2A). However, the expression of MCHR1, MRAP1 and MRAP2 show consistency in some organizations, which are all highly expressed in the brain, cerebellum and eye. These results suggest that the MCHR1, MRAP1 and MRAP2 may function together in vivo.
This finding provides evidence for co-expression of these proteins, to further determine whether MRAP2 could interact with MCHR1 protein, we performed co-immunoprecipitation (Co-IP) assays in HEK293T cells transfected with 3HA-MCHR1 and 2Flag-MRAP1 or 2Flag-MRAP2. As shown in Fig. 2B, MCHR1 did not co-purify with MRAP1 whereas co-purified with MRAP2 (Fig. 2C). Besides, we found that neither MCHR1 or MRAP proteins were observable when the purification was executed with protein A + G agarose beads only (without immunoprecipitation antibody) (Fig. 2B-C: the first two lanes in IP group). Western blot on cell lysates transfected with MCHR1, MRAP1 or MRAP2 alone was also performed to affirm that the bands observed for those proteins matched the ones detected in the IP group(Fig. 2B-C: the right half of lysate side).
Next, to further assess the ability of MCHR1 to form a functional protein complex with MRAP2 but not MRAP1 in live cells, we performed bimolecular fluorescence complementation (BiFC) assays. To achieve this goal, we generated MCHR1 constructs fused to YFP-F1 fragment in the C-terminus, while the C-terminal of MRAP1 or MRAP2 was fused to the complementary YFP-F2 fragment and Flag tag (Fig. 2D,2H). As we expected, the expression of MRAP1-Flag-F2 with MCHR1-F1 alone did not exhibit any fluorescence (Fig. 2E-G). However, YFP fluorescence was detected when co-transfected of MCHR1-F1 and MRAP2-Flag-F2 (Fig. 2I-K), suggesting that MRAP2 but not MRAP1 could form a functional complex with MCHR1, because the emission of a fluorescent YFP molecule requires the close proximity of YFP-F1 and YFP-F2 fragments.
Modulation of the surface translocation of MCHR1 by MRAP2 proteins
The trafficking of several GPCRs has been shown to be modulated by MRAP2 (Chaly et al., 2016; Chan et al., 2009; Rouault et al., 2017a; Srisai et al., 2017). To test whether MRAP2 alters the membrane translocation of MCHR1, we performed an Enzyme-linked Immunosorbent Assay (ELISA) to measure the cell surface receptors expression quantitatively. The 3HA tag was added to the N-terminus of MCHR1 and the construct was then expressed with or without MRAP2 in HEK293T cells. To determine the effect of MRAP2 on the expression of MCHR1, 3HA-MCHR1 and either MRAP2 or RAMP3 (control) at different receptor to accessory ratios (from 1:0 to 1:9) were transfected to cells. MCHR1 surface expression was detected at OD 450 nm after the addition of the HA antibody and substrate with or without cell lysis blocking. Cell numbers were measured at OD 595 nm and normalized by cell normalization stain, Janus Green. We showed that MRAP2 significant decreased the surface expression of MCHR1 compared to the control group (Fig. 3A). In general, as the ratio increased, the difference between the MRAP2 group and the control group increased (Fig. 3B). These results indicate that more than two molecules of MRAP2 interact with MCHR1, which was consistent with the results of OX1R studies. Our results showed that, at a 1:9 ratio, while MRAP2 decreased the maximal surface density of MCHR1 by > 60%, the total expression of MCHR1 was only reduced by 20% when MRAP2 existed compared to the control group (Fig. 3C). GFP-MCHR1 alone, or with MRAP2 or RAMP3 was transfected to cells to further visualize the impact of MRAP2 on MCHR1 location. We can see from the images that expression GFP-MCHR1 only or with RAMP3, GFP fluorescence was mainly present at the plasma membrane (Fig. 3D: the left and right panel), whereas GFP-MCHR1 was mostly retained in the intracellular compartment in the presence of MRAP2, which is consistent with the ELISA results. (Fig. 3D: the middle panel).
Repression of the surface expression of MCHR1 by specific regions of MRAP2
These preferential localizations suggested that MRAP2 was involved in MCHR1 trafficking to the plasma membrane. But the functional domain of MRAP2 on MCHR1 trafficking needed to be further explored. Several studies have shown that specific regions of MRAP2 are not required for the regulation of different GPCRs (Rouault et al., 2017a; Rouault et al., 2020b). To find the regions of MRAP2 that were necessary for its inhibition on MCHR1 trafficking, a series of MRAP2 mutants were generated, in which the N-terminal fragments were deleted and the C-terminal fragments were truncated (Fig. 3F). The regions of mouse MRAP2 that were truncated were same as the human MRAP2 in a previous study (Rouault et al., 2017a), since we conducted an amino acid sequence comparison in mice and human MRAP2, and found that they are high conservative, especially their TM area are completely consistent. (Fig. 3E). Using ELISA experiments, we explored the effect of each MRAP2 mutant on the surface expression of MCHR1. Our results showed that the deletion of the whole C-tail caused a loss of MRAP2 function on MCHR1 surface expression (Fig. 3G).
To explore whether the effect of these MRAP2 mutants on receptor membrane transport was caused by the disappearance of their interaction, we performed Co-IP experiments to verify the interplay between each MRAP2 mutant with MCHR1. The results showed that each MRAP2 mutant still exhibited interactions with MCHR1 (Fig. 4). Here, TM dimerization domain was retained in all MRAP2 constructs, with reference to Chen et al (Chen et al., 2020). This is important because these authors identify the TM as the dimerization domain. By retaining it the MRAP2 protein is allowed to dimerize.
Co-localization of MRAP2 mutants with MCHR1
Next, we also performed co-immunofluorescence to investigate the role of the identified functional regions of MRAP2 on MCHR1 trafficking. Different MRAP2 mutants tagged with 2Flag (WT, Δ4–11, Δ12–25, Δ26–35, Δ36–45, Δct, t97, t118 and t143) were co-transfected with 3HA-MCHR1 into cells. The results showed that the nonfunctional MRAP2 mutants were all expressed (Fig. 5), meaning that the loss of activity of these mutants was not due to lack of expression. Moreover, the coIF results confirmed again that each MRAP2 mutant had interactions with MCHR1.
The influence of MCHR1 signaling by the functional domains of MRAP2
It was previously shown that MRAP2 inhibited the signaling of OX1R and PKR1 (Rouault et al., 2017a). Based on this, we detected whether the regions of MRAP2 required for the inhibition of MCHR1 were same. To achieve this goal, we transfected HEK293T cells with MCHR1 in the presence of WT MRAP2 or MRAP2 mutants at a 1:9 ratio and measured the Ca2 + influx of MCHR1 via the CRE-luciferase reporter assay. The results showed that MRAP2 strongly decreased the efficacy of MCH (Fig. 6A-H: the blue curve). Moreover, MRAP2-Δ12–45 retained a significant inhibitory action on MCHR1 signaling (Fig. 6B-D: the black curve), while the inhibition of MRAP2 on MCHR1 signaling was almost completely reversed when transfected with MRAP2-Δ4–11 or the C-terminal truncations of MRAP2 (Fig. 6A, E-H: the black curve). Moreover, as we can see in Table1, the sensitivity of MCHR1 to MCH was reduced to varying degrees by adding MRAP2 or MRAP2 truncated constructs. Especially with the addition of mRAP2 Δ12–25 and Δ36–45 mutants, the EC50 value changed by an order of magnitude, from 10-8M to more than 10-7M.
Table 1
Statistical analysis of MCHR1 in the presence of different MRAP2 mutants in response to MCH.
Data statistics of Fig. 6A-H | LogEC50 | P value comparison |
vs. receptor alone | vs. MRAP2 WT |
MCHR1 alone | -8.149 ± 0.10 | - | < 0.0001 |
MCHR1:MRAP2 WT | -7.061 ± 0.40 | < 0.0001 | - |
MCHR1:MRAP2 Δ4–11 | -7.242 ± 0.19 | < 0.0001 | < 0.0001 |
MCHR1:MRAP2 Δ12–25 | -6.765 ± 1.15 | < 0.0001 | 0.9803 |
MCHR1:MRAP2 Δ26–35 | -7.551 ± 0.31 | < 0.0001 | < 0.0001 |
MCHR1:MRAP2 Δ36–45 | -6.643 ± 0.43 | < 0.0001 | 0.8704 |
MCHR1:MRAP2 Δct | -7.515 ± 0.54 | < 0.0001 | < 0.0001 |
MCHR1:MRAP2 t97 | -7.338 ± 0.68 | < 0.0001 | < 0.0001 |
MCHR1:MRAP2 t118 | -7.593 ± 0.63 | < 0.0001 | < 0.0001 |
MCHR1:MRAP2 t143 | -7.599 ± 0.46 | < 0.0001 | < 0.0001 |
Values were expressed as the mean ± S.E.M. of at least three independent experiments. Two-way ANOVA with Tukey post-test was applied in the statistical analysis. |
In addition, we also performed competitive inhibition experiments (Fig. 6I-P). We added different concentrations of MCHR1 antagonist (SNAP-94847) in the presence of EC80 MCH (10 − 7 M) to observe the effect of MRAP2 mutants on MCHR1. We found that, consistent with the results observed in Fig. 6A-H, the inhibition of MCHR1 signaling relied on the 4–11 amino acids and the C-terminus of MRAP2 (Fig. 6I, M-P). Except for the addition of Mrap2 Δ 26–35 and T143, IC50 of MCHR1 has a slight decrease, other variants all increase the sensitivity of MCHR1 to the inverse agonist, but in general IC50 values were not significantly different (Table 2). In short, these results indicated that the trafficking and signaling of MCHR1 were both inhibited by MRAP2, suggesting that different MRAP2 regions were involved in regulating distinct GPCR (Fig. 6Q).
Table 2
Statistical analysis of MCHR1 in the presence of different MRAP2 mutants in response to the antagonist of MCHR1.
Data statistics of Fig. 6I-P | LogIC50 | P value comparison |
vs. receptor alone | vs. MRAP2 WT |
MCHR1 alone | -8.584 ± 0.53 | - | < 0.0001 |
MCHR1:MRAP2 WT | -7.752 ± 0.68 | < 0.0001 | - |
MCHR1:MRAP2 Δ4–11 | -7.977 ± 0.57 | < 0.0001 | < 0.0001 |
MCHR1:MRAP2 Δ12–25 | -8.451 ± 0.89 | < 0.0001 | 0.2763 |
MCHR1:MRAP2 Δ26–35 | -9.278 ± 0.57 | < 0.0001 | 0.6027 |
MCHR1:MRAP2 Δ36–45 | -8.408 ± 0.68 | < 0.0001 | 0.2858 |
MCHR1:MRAP2 Δct | -7.641 ± 0.94 | 0.0002 | < 0.0001 |
MCHR1:MRAP2 t97 | -8.251 ± 0.69 | 0.0002 | < 0.0001 |
MCHR1:MRAP2 t118 | -8.314 ± 0.69 | < 0.0001 | 0.4248 |
MCHR1:MRAP2 t143 | -8.683 ± 0.51 | < 0.0001 | 0.0875 |
Values were expressed as the mean ± S.E.M. of at least three independent experiments. Two-way ANOVA with Tukey post-test was applied in the statistical analysis. |