Blog

Defining the ligand-dependent proximatome of the sigma 1 receptor

Defining the ligand-dependent proximatome of the sigma 1 receptor

By Jing Zhao, Rajalakshmi Veeranan-Karmegam, Frederick C. Baker, Barbara A. Mysona, Pritha Bagchi, Yutao Liu, Sylvia B. Smith, Graydon B. Gonsalvez, and Kathryn E. Bollinger

Excerpt from the article published in Frontiers in Cell and Developmental Biology, 11,  07 June 2023
Sec. Membrane Traffic, DOI: 1045759.https://doi.org/10.3389/fcell.2023.1045759

Editor’s Highlights

  • Sigma 1 Receptor (S1R) functions as a “pluripotent modulator” of multiple signaling pathways and therefore affects a wide range of cellular activities including calcium homeostasis, ion channel regulation, and responses to ER and oxidative stress
  • The experiments used to support direct interactions between S1R and other proteins have been mainly limited to low-throughput, candidate-based methods. 
  • The promiscuous biotin ligase Apex2 fused to the C-terminus of S1R was used to label proximal proteins under native conditions and in a ligand dependent manner.
  •  The S1R-Apex proximatome contains over 200 novel proteins, many of which localize within the endomembrane system (ER, Golgi, secretory vesicles) and function within the cellular secretory apparatus.
  • These include proteins involved in cholesterol and lipid metabolism as well as matricellular and extracellular proteins.
  • The S1R-Apex proximatome was not significantly altered upon treating cells with either a classical antagonist or agonist. It should be noted.
  • The S1R proximatome might be different in different cell types and tissues. 

Abstract

Sigma 1 Receptor (S1R) is a therapeutic target for a wide spectrum of pathological conditions ranging from neurodegenerative diseases to cancer and COVID-19. S1R is ubiquitously expressed throughout the visceral organs, nervous, immune and cardiovascular systems. It is proposed to function as a ligand-dependent molecular chaperone that modulates multiple intracellular signaling pathways. The purpose of this study was to define the S1R proximatome under native conditions and upon binding to well-characterized ligands. This was accomplished by fusing the biotin ligase, Apex2, to the C terminus of S1R. Cells stably expressing S1R-Apex or a GFP-Apex control were used to map proximal proteins. Biotinylated proteins were labeled under native conditions and in a ligand dependent manner, then purified and identified using quantitative mass spectrometry. Under native conditions, S1R biotinylates over 200 novel proteins, many of which localize within the endomembrane system (endoplasmic reticulum, Golgi, secretory vesicles) and function within the secretory pathway. Under conditions of cellular exposure to either S1R agonist or antagonist, results show enrichment of proteins integral to secretion, extracellular matrix formation, and cholesterol biosynthesis. Notably, Proprotein Convertase Subtilisin/Kexin Type 9 (PCSK9) displays increased binding to S1R under conditions of treatment with Haloperidol, a well-known S1R antagonist; whereas Low density lipoprotein receptor (LDLR) binds more efficiently to S1R upon treatment with (+)-Pentazocine ((+)-PTZ), a classical S1R agonist. Furthermore, we demonstrate that the ligand bound state of S1R correlates with specific changes to the cellular secretome. Our results are consistent with the postulated role of S1R as an intracellular chaperone and further suggest important and novel functionalities related to secretion and cholesterol metabolism.

Introduction

The sigma 1 receptor (S1R) is a small (25kD), ubiquitously expressed, transmembrane protein that is localized within the endoplasmic reticulum and its mitochondria-associated membranes (Hayashi and Su, 2005Hayashi and Su, 2007). Studies also show that it can translocate to nuclear and plasma membranes under certain conditions (Su et al., 2010Mavlyutov et al., 2015aTsai et al., 2015). S1R has been the subject of intense pharmacologic analysis over the past several decades due to the therapeutic potential of its ligands. In many instances, S1R-mediated therapeutics are advancing from bench to bedside. For example, clinical trials targeting S1R for treatment of neurodegenerative diseases including Amyotrophic Lateral Sclerosis (ALS) and Huntingtons Disease are ongoing (Arnold, 2021). In addition, promising preclinical studies indicate that S1R may offer a treatment target for visual disorders including glaucoma, retinal degeneration and traumatic optic neuropathy (Mavlyutov et al., 2015bNguyen et al., 2015Wang et al., 2016Sambo et al., 2017Geva et al., 2021Li et al., 2021). Furthermore, in recent studies, S1R was shown to link the SARS-CoV2 replicase/transcriptase complex to the ER membrane by binding directly to nonstructural protein 6 (NSP 6) (Gordon et al., 2020). Therefore, S1R ligands may provide antiviral activity against severe acute respiratory syndrome CoV-2 (SARS-CoV-2) (Gordon et al., 2020Vela, 2020). Finally, there is interest in using S1R ligands for treating and imaging cancer (Happy et al., 2015Das et al., 2016).

Despite intense interest, the molecular mechanisms that underlie effects of S1R ligands are not well understood. In general, agonists for S1R show pro-survival effects whereas antagonists inhibit tumor cell proliferation and induce apoptosis (Maurice and Su, 2009). Studies indicate that S1R functions as a “pluripotent modulator” of multiple signaling pathways and therefore affects a wide range of cellular activities including calcium homeostasis, ion channel regulation, and responses to ER and oxidative stress (Hayashi and Su, 2003). The proposed general mechanism for S1R function is through protein-protein interactions (Hayashi and Su, 2007). In support of this paradigm, S1R has been reported to bind to at least 49 different proteins that generally show diverse structure and function (Schmidt and Kruse, 2019). However, published experiments used to support direct interactions between S1R and other proteins have been mainly limited to low-throughput, candidate-based methods. These include co-immunoprecipitation and resonance energy transfer experiments (Navarro et al., 2010Balasuriya et al., 2013Tsai et al., 2015). Previous studies have also utilized proximity ligation assays, but few have combined these evaluations with high-throughput proteomic analyses.

In this study, a cell line-based proximity biotin labeling assay was developed and combined with proteomic identification (Rhee et al., 2013). The promiscuous biotin ligase Apex2 fused to the C-terminus of S1R was used to label proximal proteins under native conditions and in a ligand dependent manner (Lam et al., 2015Hung et al., 2016). The biotinylated proteins were then purified and identified using quantitative mass spectrometry. The labeling radius of Apex2 is thought to be around 20 nm (Trinkle-Mulcahy, 2019). Thus, although this strategy does not necessarily identify direct protein interaction partners, biotinylated proteins are proximal to the bait protein and may be present in a complex along with it. We refer to the proteins identified in this manner as the S1R-Apex proximatome. Under native conditions, we find that S1R-Apex proximatome contains over 200 novel proteins, many of which localize within the endomembrane system (ER, Golgi, secretory vesicles) and function within the cellular secretory apparatus. In addition, under conditions of cellular exposure to either S1R agonist or antagonist we identify proximatome changes that highlight the molecular pathways critical to S1R-mediated ligand-dependent effects. These include proteins involved in cholesterol and lipid metabolism as well as matricellular and extracellular proteins. Thus our results offer insight into the cellular changes that accompany S1R ligand binding.

Results

Defining the proximatome of the Sigma1 receptor

In order to identify proteins central the ligand-dependent function of Sigma1 receptor (S1R), we employed a proximity biotinylation approach. With this strategy, the protein of interest is tagged with a promiscuous biotin ligase. Proteins that are present within 20 nm of the tagged protein become biotinylated in vivo and can be easily purified using streptavidin conjugated beads (Figure 1A) (Cho et al., 2020). One of the main advantages of this approach is that because the biotin-streptavidin interaction is extremely strong, the purification can be done using harsh wash buffers. This minimizes binding of non-specific proteins to beads, thus reducing the potential for false-positive hits (Samavarchi-Tehrani et al., 2020). For our studies, we chose to tag S1R at the C-terminus with Apex2 (Schmidt et al., 2016Mavlyutov et al., 2017Schmidt et al., 2018). We chose this approach because previous studies have shown that S1R-Apex is functional and produces the expected localization pattern (Mavlyutov et al., 2017). In addition, a structural analysis of S1R revealed that only a few N-terminal residues are exposed within the cytoplasm. The bulk of the protein, including the C-teminus, resides within the ER (Schmidt et al., 2016). Thus, tagging S1R on the C-terminus with Apex2 was the most logical approach. A V5 epitope was incorporated into the tag in order to enable Western blotting and immunofluorescence analysis using commercially available antibodies. A construct expressing GFP-Apex was used as the negative control.

FIGURE 1.
Generation and validation of S1R-Apex cells. 
(A)
 Schematic of the proximity biotinylation strategy (B) HeLa cells stably expressing S1R-Apex were fixed and processed for immunofluorescence using a V5 antibody (green). S1R-Apex localizes around the nuclear envelope and to ER tubules. (C) HeLa cells stably expressing S1R-Apex were fixed and processed for immunofluorescence using a V5 antibody (green). The cells were also incubated with Streptavidin-647 (red) to reveal the localization of biotinylated proteins and were counterstained using DAPI (cyan). Biotinylated proteins co-localize with S1R-Apex. (D) Control HeLa cells not expressing a biotin ligase were processed using Streptavidin-647 (red). The cells were counterstained with DAPI (cyan). Minimal Streptavidin signal is observed in these cells. The scale bar in B is 15 microns and the scale bar in C is 20 microns.

HeLa cells stably expressing S1R-Apex and GFP-Apex were generated using CRISPR-based incorporation of the constructs at the AAVS1 safe harbor locus using a published protocol (Oceguera-Yanez et al., 2016). Immunostaining of cells expressing S1R-Apex produced the expected localization pattern (Figure 1B). In order to validate the proximity biotinylation approach, cells expressing S1R-Apex were labeled using biotin phenol. Control cells that did not express a biotin ligase were similarly treated. The cells were then processed for immunofluorescence using a V5 antibody for detecting the Apex tagged protein, and Streptavidin conjugated Alexa647 for detecting biotinylated proteins. In contrast to control cells, robust biotinylation signal was observed in cells expressing S1R-Apex (Figures 1C, D). Furthermore, the biotinylation signal co-localized with the V5 signal for S1R-Apex (Figure 1C). The biotinylation pattern is consistent with the published localization of S1R to ER and nuclear membranes (Jiang et al., 2006Mavlyutov et al., 2015aMavlyutov et al., 2017). By contrast, cells expressing GFP-Apex displayed a nuclear and cytoplasmic biotinylation pattern (Supplementary Figure S1). Based on these results, we conclude that S1R-Apex is localizing as expected and biotinylating proximal proteins.

Next, cells expressing GFP-Apex and S1R-Apex were grown and labeled for proteomics analysis. Lysates were prepared from the labeled cells and the biotinylated proteins were purified using high-capacity Streptavidin agarose beads. After extensive wash steps, the bound proteins were eluted using Trypsin digestion and processed for mass-spectrometry. The entire experiment was done in triplicate. Proteins that were enriched at least two-fold in the S1R-Apex pellet in comparison to GFP-Apex and had a p-value of at least 0.05 were considered to be part of the S1R proximatome (Figure 2ASupplementary Table S1). Using these criteria, the S1R proximatome contains 233 specific proteins in HeLa cells. The top 60 candidates are shown in Table 1.

FIGURE 2.
The S1R proximatome. 
(A)
 A volcano plot depicting the S1R proximatome. A line demarcating 2-fold enrichment and a p-value of 0.05 is shown. The grey shaded box indicates proteins that are at least 2-fold enriched with S1R in comparison to the GFP-Apex control with a p-value of at least 0.05. (B) A Cellular Component GO analysis of S1R proximal proteins. (C) A Biological Process GO analysis of S1R proximal proteins. (D,E) Biotinylated proteins were purified from HEK293T cells stably expressing either GFP-Apex or S1R-Apex using streptavidin conjugated beads. Biotinylated proteins were also purified from HEK293T cells stably expressing S1R-Apex that were treated with either Haloperidol or (+)-PTZ. The bound proteins were eluted and analyzed by Western blotting using an antibody against PDI (D) or Calnexin (E). S1R-Apex expressed in HEK293T cells is able to biotinylate PDI and Calnexin. This proximal association is not altered by treatment with Haloperidol or (+)-PTZ.
TABLE 1.
A list of the top 60 S1R proximal proteins.

S1R is known to exist as a trimeric complex in vivo (Schmidt et al., 2016). Thus, peptides corresponding to S1R should be highly enriched in the S1R-Apex dataset. This was indeed the case (Figure 2A). In addition, the ER chaperone BiP (GRP78), a known S1R interacting partner, was also specifically enriched in the S1R-Apex pellet (Figure 2A) (Hayashi and Su, 2007). Previous studies have shown that S1R localizes to the endoplasmic reticulum (ER), to nuclear lamellae, to sites of ER-mitochondria contact and to the plasma membrane (Hayashi and Su, 2007Su et al., 2010Tsai et al., 2015Mavlyutov et al., 2017). Consistent with these studies, a “cellular component” Gene Ontology (GO) analysis indicates that the S1R proximatome is highly enriched for proteins that localize to the nuclear membrane, to the ER lumen and membrane, and to organelles involved in protein secretion. (Figure 2B). Previous studies have implicated a role for S1R in the ER stress response and as a potential molecular chaperone. In line with these findings, a “biological process” GO analysis is enriched for terms such as “protein targeting to ER” and “response to unfolded protein.” In addition, proteins with a role in N-linked glycosylation are highly enriched within the S1R proximatome (Figure 2C).

In order to validate these results, we generated HEK293T cells that stably express either GFP-Apex or S1R-Apex integrated at the AAVS1 safe harbor locus (Oceguera-Yanez et al., 2016). Two ER proteins that were highly enriched in the S1R proximatome were Protein Disulfide Isomerase (PDI) and the chaperone Calnexin. Consistent with the proteomics results, PDI and Calnexin were specifically biotinylated and precipitated in HEK293T cells expressing S1R-Apex (Figures 2D, E). Biotinylation of these proteins by S1R-Apex was unchanged upon treating cells with either Haloperidol, an S1R antagonist, or with (+)-PTZ, an S1R agonist (Figures 2D, E).

We next determined whether PDI and Calnexin were capable of co-precipitating with S1R. Using a standard co-immunoprecipitation protocol involving a lysis buffer containing 0.2% CHAPS, we were able to detect a small amount of co-precipitaitng Calnexin but not PDI (Supplementary Figure S1B and data not shown). Surprisingly, using this approach, we were not able to efficiently co-precipitate Bip along with S1R-Apex (Supplementary Figure S1B). Thus, although Calnexin, PDI and Bip were biotinylated by S1R-Apex, detecting these interactions via co-immunoprecipitation proved much more challenging. It is possible that the interactions are transient or weak and are thus not retained during the process of co-immunoprecipitation.

S1R-apex biotinylates components of the ER translocation machinery

For proteins that contain a signal sequence, import into the ER occurs co-translationally. Once the signal sequence has been translated, it is bound by the signal recognition particle (SRP). This complex is then docked onto the surface of the ER by binding to the SRP receptor. Next, the docked ribosome associates with the Sec61 translocation complex and its accessory proteins. This enables the translating peptide to be imported into the ER. Once in the ER, luminal chaperones such as BiP and Calnexin help fold the protein into its native conformation (Figure 3A) (Lang et al., 2017Linxweiler et al., 2017). Multi-pass transmembrane proteins are imported into the ER with the aid of the Nicalin-TMEM147-NOMO complex (McGilvray et al., 2020). Numerous proteins within the S1R proximatome function at various steps in the ER translocation process (Figure 3B). This finding was quite intriguing given the proposed function of S1R as a molecular chaperone. These results position S1R at the entry way into the ER, an ideal location for a molecular chaperone.

FIGURE 3.
S1R biotinylates components of the ER translocation complex. 
(A)
 A schematic of the components involved in translocation of proteins into the ER. (B) A list of ER translocation components identified in the S1R proximatome. (C) Biotinylated proteins were purified from HEK293T cells stably expressing either GFP-Apex or S1R-Apex that were untreated or treated with either Haloperidol or (+)-PTZ. The proteins were purified using streptavidin conjugated beads. The bound proteins were eluted and analyzed by Western blotting using an antibody against Sec61alpha. S1R-Apex expressed in HEK293T cells biotinylates Sec61alpha. (D) HEK293T cells were transiently transfected with a construct expressing either GFP-Apex or Sec61 beta-Apex. Biotinylated proteins were purified using streptavidin conjugated beads, the bound proteins were eluted and analyzed by Western blotting using an antibody against endogenous S1R. Endogenous S1R is biotinylated by Sec61 beta-Apex, a component of the ER translocation complex. (E)HEK293T cells stably expressing GFP or S1R-Apex were either untreated (DMSO) or treated with the drug Eeyarestatin. Eeyarestatin blocks the transport of proteins across the Sec61 translocation channel. Biotinylated proteins were purified using streptavidin conjugated beads and the bound proteins were eluted and analyzed by Western blotting using antibodies against Sec61 alpha and Calnexin. (F) The above experiment was repeated in triplicate and the binding of S1R-Apex with Sec61 alpha and Calnexin was quantified. Unpaired t tests were used for these analyses; ****p < 0.0001, ns = not significant. Blocking protein transport across the Sec61 channel disrupts the proximal association between S1R and Sec61 alpha.

In order to validate these findings, we focused our studies on Sec61alpha, the protein that forms the central translocation channel. Consistent with the proteomics result, Sec61alpha was specifically detected in the S1R-Apex pellet from HEK293T cells (Figure 3C). As with PDI and Calnexin, the biotinylation of Sec61alpha by S1R-Apex was unchanged upon treating cells with either Haloperidol or (+)-PTZ. As further validation of this result, we determined whether endogenous untagged S1R was also present in a complex with the ER translocation channel. For this experiment, HEK293T cells were transiently transfected with either GFP-Apex or Sec61beta-Apex. Sec61beta-Apex was chosen for this experiment because this construct has been used in a previous publication (Lee et al., 2016). In addition, we were concerned that tagging Sec61alpha with Apex might compromise its function and would negatively impact import of proteins into the ER. After binding and wash steps, the bound proteins were eluted and analyzed by Western blotting using an antibody against S1R. Consistent with our proteomics results which suggests proximity between S1R-Apex and the ER translocation channel, significantly more S1R was biotinylated by Sec61beta-Apex in comparison to GFP-Apex (Figure 3D).

We next determined whether the proximal association between S1R and the Sec61 complex was sensitive to active protein import into the ER. The drug Eeyarestatin is known to inhibit Sec61 alpha, and as a consequence, ER protein import is blocked (Wang et al., 2008Wang et al., 2010Gamayun et al., 2019). Interestingly treating cells with Eeyarestatin reduced the biotinylation of Sec61alpha by S1R-Apex (Figures 3E, F). The biotinylation of Calnexin by S1R-Apex was also somewhat reduced. However, these results were not statistically significant (Figures 3E, F). Based on these results, we conclude that S1R localizes promixal to the ER translocation machinery and likely does so in an import-dependent manner.

S1R-apex biotinylates components at the ER-Golgi intermediate compartment

Another category of proteins enriched in the S1R proximatome are proteins that localize to the ER-Golgi intermediate compartment (ERGIC) (Figure 2B). Proteins that localize at this site are involved in trafficking between the ER and Golgi compartments and often have a role in the secretory process (Appenzeller-Herzog and Hauri, 2006). The S1R proximatome contains several proteins that are secreted or localized on the plasma membrane such as Collagens, Integrins, and Low Density Lipoprotein Receptor (Supplementary Table S1). It is therefore possible that S1R plays an active role in the secretory process. In order to validate this finding, the binding experiment was repeated, and the pellets were analyzed by Western blotting using an antibody against Lman1, a marker protein of the ERGIC and a protein that was highly enriched in the S1R proteomics dataset (Supplementary Table S1) (Schindler et al., 1993). Lman1 was specifically biotinyated by S1R-Apex irrespective of the ligan bound state of receptor (Figure 4A). In addition, we observed significant co-localization between S1R and GFP-Lman1 (Figure 4B).

FIGURE 4.
S1R biotinylates the ER-Golgi intermediate compartment component Lman1. 
(A)
 Biotinylated proteins were purified using streptavidin conjugated beads from HEK293T cells stably expressing either GFP-Apex or S1R-Apex that were untreated or treated with either Haloperidol or (+)-PTZ. The bound proteins were eluted and analyzed by Western blotting using an antibody against Lman1. S1R-Apex expressed in HEK293T cells biotinylates the ERGIC component Lman1. (B) HeLa cells stably expressing S1R-Apex were transfected with a plasmid expressing GFP-Lamn1 (green). The cells were fixed and processed for immunofluorescence using a V5 antibody (red). S1R-Apex co-localizes with GFP-Lman1 in the area of the ER-Golgi intermediate compartment.

Ligand dependent proximatome of S1R

Having established the S1R proximatome under native conditions, we next wished to determine how the proximatome would change when cells were treated with Haloperidol, an S1R antagonist, or with (+)-PTZ, an agonist of S1R (de Costa et al., 1989Hayashi and Su, 2004Vavers et al., 2019). For this experiment, HeLa cells expressing S1R-Apex were either untreated, treated with 25uM Haloperidol for 24 h, or treated with 20uM (+)-PTZ for 24 h. The localization of S1R was relatively unchanged under these treatment conditions (Supplementary Figures S2A–C). The treated cells were labeled, the biotinylated proteins were purified using streptavidin agarose and the bound proteins were analyzed using mass-spectrometry. As before, the entire experiment was done in triplicate.

Somewhat surprisingly, treatment of cells with Haloperidol or (+)-PTZ did not change the biotinylation profile for the vast majority of proteins (Figures 5A, BSupplementary Tables S2, S3). For this analysis, we considered proteins that were enriched at least two-fold (1.0 fold in the log2 scale) in the drug treated sample in comparison to the control and with a p-value of at least 0.05 as being significantly changed. Using these criteria, eight proteins displayed increased biotinylation by S1R-Apex in the presence of Haloperidol (Figure 5A) and thirteen proteins displayed increased biotinylation by S1R-Apex in the presence of (+)-PTZ (Figure 5B). Interestingly, even within this small set, three proteins were shared between the two drug-treated samples, FDFT1 which encodes Squalene synthase, PCKS9 which encodes Proprotein Convertase Subtilisin/Kexin Type 9, and GARS which encodes glycyl-tRNA-synthetase. Although Haloperidol is considered an antagonist and (+)-PTZ is considered an agonist of S1R, a recent structural study found that both compounds bind within the same binding pocket of S1R (Schmidt et al., 2018). Thus, in addition to unique structural changes that might be induced upon ligand binding, there might also be some common changes that are promoted by Haloperidol and (+)-PTZ.

FIGURE 5.
The ligand-dependent S1R proximatome. 
(A)
 A volcano plot comparing the proximatome of S1R-Apex under untreated conditions versus treatment with 25uM Haloperidol for 24 h. A line demarcating 2-fold enrichment and a p-value of 0.05 is shown. The grey shaded box indicates proteins that show at least 2-fold greater binding to S1R under Haloperidol treatment conditions with a p-value of at least 0.05. A list of these proteins is shown below the volcano plot. (B) A volcano plot comparing the proximatome of S1R-Apex under untreated conditions versus treatment with 20uM (+)-PTZ for 24 h. The layout is similar to panel (A). The proteins in blue were shown to have an increased proximal association with S1R under both Haloperidol and (+)-PTZ treatment conditions. (C–E) Biotinylated proteins were purified from HeLa cells stably expressing either GFP-Apex or S1R-Apex that were untreated or treated with either Haloperidol or (+)-PTZ. The proteins were purified using streptavidin conjugated beads. The bound proteins were eluted and analyzed by Western blotting using antibodies against PCSK9 (C), LDLR (D) or Thrombospondin1 (E). Haloperidol treatment increased the proximal association between S1R and PCSK9, whereas (+)-PTZ treatment increases the proximal association between S1R LDLR and Thrombospondin1.

For follow-up studies, we decided to focus our efforts primarily on PCSK9 and LDLR. PCSK9 (Proprotein Convertase Subtilisin/Kexin Type 9) displays increased biotinylation by S1R-Apex under conditions of Haloperidol treatment and to a lesser extent under PTZ treatment conditions, whereas LDLR (Low density lipoprotein receptor) is biotinylated more efficiently to S1R-Apex upon treatment with (+)-PTZ. Disease-causing variants in PCSK9 and loss of function mutations in LDLR are associated with familial hypercholesterolemia (Brown and Goldstein, 1986Abifadel et al., 2003). PCSK9 is a secreted protein and extracellular PCSK9 binds to LDLR present on the cell surface. This results in endocytosis of LDLR and targeting of LDLR to the lysosome for degradation (Weinreich and Frishman, 2014). Increased turnover of LDLR results in high serum cholesterol levels (Cohen et al., 2006). Thus, both proteins play an important role in cholesterol metabolism. In addition, FDFT1, which encodes Squalene synthase, and CYP51A1, which encodes Lanosterol 14-alpha demethylase, are also a critical players in the cholesterol biosynthetic pathway (Shechter et al., 1992Lorbek et al., 2012).

Cells expressing either GFP-Apex, S1R-Apex (untreated), or S1R-Apex treated with either Haloperidol or (+)-PTZ were labeled and the binding reaction was performed as before. Bound proteins were eluted and analyzed by Western blotting using an antibody against PCSK9 (Figure 5C) or LDLR (Figure 5D). Consistent with the proteomics results, treatment of cells with Haloperidol increased the amount of PCSK9 that was labeled and precipitated by S1R-Apex (Figure 5C). (+)-PTZ treatment also resulted in an enrichment of PCKS9 in the S1R-Apex bound sample. Similarly, treatment with (+)-PTZ increased the amount of LDLR that was labeled and precipitated in the S1R-Apex sample (Figure 5D). Another protein that also displayed an increased biotinylation by S1R-Apex in the presence of (+)-PTZ was Thrombospondin1 (THBS1), a finding that was also validated by Western blotting (Figures 5B, E).

S1R regulates the cellular secretome

In order to follow up on these results, we examined the localization of PCSK9 in control cells or in cells treated with either Haloperidol or (+)-PTZ. PCSK9 was localized in small cytoplasmic foci. Although the pattern of localization was unchanged upon drug treatment, the level of PCKS9 was elevated upon treatment with Haloperidol and to a lesser extent with (+)-PTZ (Figures 6A–C). This result was confirmed using Western blot analysis (Figures 6H, J). Given the increased biotinylating of PCSK9 by S1R-Apex after drug treatment, we evaluated potential co-localizaiton between these two proteins. Partial co-localization was observed between PCSK9 and S1R-Apex under native conditions (Figures 6D, G). Moreover, the level of co-localizaiton was increased in cells treated with either Haloperidol or (+)-PTZ (Figures 6E–G).

FIGURE 6.
The effect of S1R ligand binding on PCSK9 and Thrombospondin1. 
(A–C)
 HeLa cells were untreated (A) or treated with either Haloperidol (B) or (+)-PTZ (C) for 24 h. The cells were then fixed and processed for immunofluorescence using an antibody against PCSK9. The cells were also counterstained using Alexa-633 conjugated Phalloidin to reveal the actin cytoskeleton (cyan). The signal for PCSK9 is depicted using a red to white range indicator, also known as a look-up-table (LUT). Low intensity pixels are shown in red and high intensity pixels are shown in white. The scale bar is 20 microns. (D–F) HeLa cells were untreated (D) or treated with either Haloperidol (E) or (+)-PTZ (F) for 24 h. The cells were then fixed and processed for immunofluorescence using antibodies against PCSK9 (cyan) and V5 (magenta). Individual channels as well as a merged image is shown. The scale bar is 20 microns. (G) Co-localization between S1R (V5) and PCKS9 was analyzed using Imaris 10. A Pearson’s coefficient is shown. Using this analysis method, perfect co-localization has a value of 1.0. S1R partially co-localizes with PCKS9 and the level of co-localizaiton is increased under conditions of ligand treatment. n = 12 images per condition. A one-way Anova was used for these analyses; ****p < 0.0001. (H) Lysates were prepared from HeLa cells that were either untreated or treated with either Haloperidol or (+)-PTZ. The lysates were analyzed by Western blotting using antibodies against PCKS9 and GAPDH. (I) HeLa cells were left untreated or were treated with either Haloperidol or (+)-PTZ for 24 h. For the last 4 h of this treatment, the media was replaced with serum free medium. The cells were then cultured for the final 4 h. The culture supernatant was collected, concentrated and analyzed by Western blotting using antibodies against either PCSK9 or Thrombospondin1. (J)The experiment in panel H was repeated in triplicate and the level of PCSK9 was quantified relative to the control untreated sample. A one-way Anova was used for these analyses; **p < 0.01, ns = not significant. Haloperidol treatment results in the upregulation of PCSK9 levels. (K,L)The experiment on panel I was repeated in triplicate and the level of secreted PCSK9 (K) and Thrombospondin1 (L) was quantified relative to the untreated sample. A one-way Anova was used for these analyses; ***p < 0.001, ns = not significant. Haloperidol treatment results in increased secretion of PCKS9. However, the sercretion of Thrombospondin1 is decreaed under these same conditions.

PCSK9 is a secreted protein. We therefore monitored the secretion of PCSK9 under conditions of Haloperidol treatment. Culture supernatants were collected and analyzed by Western blotting. In comparison to untreated and (+)-PTZ treated cells, culture supernatant from Haloperidol treated cells contained significantly more PCSK9 (Figures 6I, K). In order to determine whether Haloperidol treatment increased global secretion we examined culture supernatants using an antibody against Thrombospondin1. In contrast to PCSK9, Thrombospondin1 was slightly decreased in culture supernatants from Haloperidol treated cells (Figures 6I, L). These findings suggest that Haloperidol does not globally affect protein secretion and that the effect on PCSK9 is likely to be specific. Interestingly, under these conditions, (+)-PTZ treatment also reduced the secretion of Thrombospondin1 (Figures 6I, L).

We next compared the global secretome from control HeLa cells versus cells treated with either Haloperidol or (+)-PTZ. For this experiment, cells were cultured in serum containing media for 20 h, after which time the media was removed and replaced with serum free media for 4 h. As with previous experiments, cells were either untreated or treated with Haloperidol or (+)-PTZ for 24 h. Conditioned medium was collected, concentrated, and analyzed by mass spectrometry. A total of 1,087 proteins were detected in culture supernatant (Supplementary Table S4). The vast majority of these proteins were unchanged in their secretion upon treatment with either Haloperidol or (+)-PTZ. However, 15 proteins, including PCSK9, showed a greater than 2 fold increase in their secretion upon Haloperidol treatment (Figure 7ASupplementary Table S5). Only 1 protein was reduced in its secretion more than 2 fold upon Haloperidol treatment (Figure 7ASupplementary Table S5). Consistent with our prior observation via Western blot analysis (Figure 6I), we detected a statistically significant decrease in the level of Thrombospondin1 upon Haloperidol treatment (Figure 7A). However, the difference by mass spectrometry analysis of conditioned media was less than 2 fold (Supplementary Table S5). In (+)-PTZ treated samples, 9 proteins showed a greater than 2 fold increase in secretion (Figure 7BSupplementary Table S6). Interestingly, 3 of these proteins also displayed increased secretion under Haloperidol treatment conditions (Figures 7A, B, blue font). (+)-PTZ treatment resulted in a greater than 2 fold decrease in the secretion of 3 proteins (Figure 7BSupplementary Table S6). Collectively, these results suggest that the ligand bound state of S1R can specifically alter the cellular secretome.

FIGURE 7.
Global secretome analysis upon ligand treatment. 
(A)
 A volcano plot comparing the secretome between control cells versus cells treated with 25uM Haloperidol for 24 h. A line demarcating 2-fold enrichment (+1 and −1 on the x-axis) and a p-value of 0.05 is shown (1.3 on the y-axis). A list of proteins whose secretion is affected greater than two fold is shown below the volcano plot. (B) A volcano plot comparing the secretome between control cells versus cells treated with 20uM (+)-PTZ for 24 h. The layout is similar to panel (A). The proteins in blue were shown to be affected in their secretion under both treatment conditions. (C) A pie chart showing the total number of secreted proteins, the number of proteins found to be biotinylated by S1R-Apex, and the overlap between the two categories.

S1R, LDLR and cholesterol

We next turned our attention to LDLR. In control cells, LDLR was diffusely localized within the cytoplasm with occasional localization within small puncta (Figure 8A). Interestingly, Haloperidol treatment caused a dramatic relocalization of LDLR to large intracellular foci (Figure 8B). A similar, but milder phenotype was observed in (+)-PTZ treated cells (Figure 8C). Thus, although Haloperidol and (+)-PTZ do not affect the cellular level of LDLR, both drugs affect the intracellular localization of the receptor. Co-staining cells with antibodies recognizing LDLR and S1R-Apex indicated a partial co-localization between both proteins (Supplementary Figures S2D–G). However, unlike with PCKS9, the level of co-localizaton between S1R and LDLR was unchanged upon treatment of cells with either Haloperidol or (+)-PTZ (Supplementary Figures S2E–G).

FIGURE 8.
The effect of S1R ligand binding on LDLR, LDL uptake and Cholesterol. 
(A–C)
 HeLa cells were untreated (A) or treated with either Haloperidol (B) or (+)-PTZ (C) for 24 h. The cells were then fixed and processed for immunofluorescence using an antibody against LDLR. The cells were also counterstained using Alexa-633 conjugated Phalloidin to reveal the actin cytoskeleton (cyan). The signal for LDLR is depicted using a red to white look-up-table (LUT). Haloperidol treatment causes LDLR to accumulate within large intracellular puncta. A similar but milder phenotype was observed in (+)-PTZ treated cells. The scale bar is 20 microns. (D–F) HeLa cells were untreated (G) or treated with either Haloperidol (H) or (+)-PTZ (I) for 24 h. The cells were then incubated with fluorescently labeled LDL. After 2 h of incubation, the cells were fixed and counterstained with DAPI to reveal nuclei (cyan). The signal for LDL is shown using the red to white LUT. The scale bar is 20 microns. (G) The LDL assay was repeated in triplicate. The fluorescence intensity per cell was calculated for twenty frames for each treatment condition. This was used to determine the mean intensity for each treatment. The LDL uptake value for the Haloperidol and (+)-PTZ treated samples are reported relative to the control sample. A one-way Anova was used for these analyses; ****p < 0.0001, ns = not significant. The ability of cells to endocytose LDL is increased upon treatment with either Haloperidol or (+)-PTZ. (H) Cells were left untreated or were treated for 48 h with either Haloperidol or (+)-PTZ. Lipids were extracted from the cells and the cholesertol content was measured using the Amplex Red kit. The mean value for each condition was determined. The mean cholesterol content in the treated samples is reported relative to the mean of the untreated sample. In comparison to the control, Haloperidol treated cells contain a higher cholesterol content. Although a trend towards higher cholesterol content was also observed in the (+)-PTZ treated cells, this result was not statistically significant. A one-way Anova was used for these analyses. *p < 0.05, ns = not significant.

LDLR is present on the cell surface where it binds to extracellular LDL resulting in its endocytosis. Internalized LDL is trafficked through the endocytic pathway to the lysosome. Lysosomal degradation of LDL releases cholesterol, which can then be used by the cell (Goldstein and Brown, 1977Brown and Goldstein, 1997). We therefore examined the ability of cells that were treated with Haloperidol or (+)-PTZ to endocytose labelled LDL. In comparison to untreated cells, both (+)-PTZ and Haloperidol treatment resulted in increased LDL uptake (Figures 8G–J). This result was somewhat surprising for Haloperidol treated cells. Higher levels of PSCK9 are generally thought to result in lower LDLR levels and therefore reduced LDL uptake. Possible scenarios that might explain this unanticipated result are discussed below. We next measured the cholesterol content of cells treated with either (+)-PTZ or Haloperidol. In comparison to control, the cholesterol content was higher in cells treated with Haloperidol (Figure 8H). A similar trend of higher cholesterol content was also observed in (+)-PTZ treated cells. However, under the conditions of this experiment, the difference was not statistically significant (Figure 8H).

Collectively, our findings reveal that the S1R proximatome consists of over 200 proteins. Ligand dependent changes involve proteins that are involved in cholesterol metabolism as well as proteins that are secreted into the extracellular space. The therapeutic benefit of S1R activation likely involves these and other proximatome changes.

Discussion

We report here use of a cell line-based proximity labeling assay to probe the S1R proximatome under native conditions and in a ligand-dependent manner. This methodology, combined with quantitative mass spectrometry-based proteomics, identified over 200 S1R-proximal proteins. Proteins identified in this manner do not necessarily represent direct interaction partners. However, given the limited labeling radius of Apex2, it is likely that many of the identified proteins are present in a complex with S1R. However, additional experiments will be required to determine whether these proximal associations represent true biochemical interactions. The advantage of proximity biotin labeling is that weak interactions, which might otherwise be missed in a classical purification approach, might be identified using this strategy. Consistent with this notion, although we were able to easily biotinylate and purify Bip, Calnexin, and Protein Disulfide Isomerase (PDI) using S1R-Apex, these proteins did not efficiently co-immunoprecipitate with S1R-Apex (Supplementary Figure S1B and data not shown).

The recently published crystal structure of S1R shows that it consists of a single transmembrane segment with a short cytoplasmic tail and a large luminal ligand-binding domain (Schmidt et al., 2016). High resolution electron microscopy studies using an S1R-GFP-APEX2 fusion protein are consistent with luminal localization of the ligand-binding S1R C-terminus (Mavlyutov et al., 2017). Importantly, our experiments also involved tagging S1R with Apex on the C-terminus of the protein. Our results indicate that the S1R proximatome consists mostly of proteins residing within the ER lumen, the ER membrane, or within the secretory pathway. This provides functional evidence for a luminal localization of the ligand-binding domain of S1R, and is consistent with structural studies.

Previous studies indicate that S1R can influence the cellular function of many proteins (Su et al., 2016Schmidt and Kruse, 2019). A proposed mechanism for S1R’s multiple intracellular effects has been through direct protein interactions, and S1R has been termed a ligand-operated molecular chaperone (Hayashi and Su, 2007Su et al., 2010). In support of this model, experiments performed in Chinese Hamster Ovary (CHO) cells indicate that S1R forms a complex with the chaperone protein, BiP, which is known to play a central role in protein folding and quality control (Hayashi and Su, 2007). Studies also showed that activation of S1R by the agonist (+)PTZ, led to dissociation of S1R from BiP. Release from BiP was proposed to trigger S1R’s multiple interactions with client proteins (Hayashi and Su, 2003Hayashi and Su, 2007).

Our results (Figure 2) suggest a proximal association between S1R and BiP, consistent with previous studies. However, we did not observe a difference in the biotinylation status of Bip by S1R-Apex when cells were exposed to either Haloperidol or (+)-PTZ (Figure 5). In addition, in native condition HeLa cells (without antagonist or agonist treatment), S1R-Apex biotinylates not only Bip, but also hundreds of other proteins (Figure 2). Furthermore, the number of biotinylated proteins does not substantially increase with exposure to ligands (Figure 5). In contrast to previous studies, our results offer a comprehensive and quantitative view of the S1R proximatome, and future work should address how S1R-protein interactions vary dependent on cell type analyzed, ligand specificity and exposure time as well as cell culture conditions.

Our data do show a large, robust and reproducible S1R-Apex proximatome, consistent with the paradigm of S1R-mediated chaperone function. Our studies indicate a proximal association between S1R and the ER translocation complex. We have validated the biotinylation of Sec61alpha by S1R-Apex and show that this proximal association is sensitive to protein import into the ER. Overall, consistent with previous work, our studies support ER-localized chaperone functionalities as central to the mechanism of S1R’s pleotropic effects.

A recent study by Zhemkov and colleagues also reported on the S1R proximatome in native HeLa cells (Zhemkov et al., 2021a). Our results are broadly consistent with theirs. For instance, numerous proteins such as PDI, Calnexin, Integrin, and PCKS9 were identified as S1R proximal proteins in both studies. However, the proximatomes were not completely overlapping. The study by Zhemkov et al., used a Tetracyline inducible expression system for expressing S1R-Apex, whereas in this study S1R-Apex was expressed constitutively after integration at the AAVS1 safe harbor locus. In addition, the binding control used in the Zhemkov et al. study corresponeded to uninduced cells (i.e., no expression of S1R-Apex), whereas in the current study, cells expressing GFP-Apex were used as the control.

In addition to evaluating the native S1R proximatome, our study assessed proximatome changes in cells treated with S1R ligands. For these studies, we used two classic ligands with high affinity for S1R. Traditionally, these chemicals have been classified as S1R agonists or antagonists based on their ability to recapitulate the effects of genetic knockout (KO) or knockdown of S1R. For example, ligands that mimic S1R genetic KO are considered antagonists (Nguyen et al., 2015). S1R KO animals are viable and fertile, but show signs and symptoms consistent with exacerbation of age-associated neurodegenerative disorders (Sabino et al., 2009). The KO studies collectively highlight the innate prosurvival and neuroprotective properties of S1R activity. Therefore, S1R ligands that promote cellular survival and neuroprotection are generally considered agonists.

The two ligands in our studies, (+)-Pentazocine, (+)-PTZ, an agonist, and Haloperidol, an antagonist, were each used in the original structural analyses that characterized S1R ligand recognition (Schmidt et al., 2018). PTZ is an FDA-approved analgesic that exists as a racemic mixture of dextro-(+) and levo-(−) isomers. The levo-(−) isomer shows high affinity for opioid receptors and the dextro-(+) isomer (referred to in this study as (+)-PTZ) shows high affinity for S1R. Effects of (+)-PTZ are generally neuroprotective and pro-survival, and it is considered the prototypical S1R agonist (Walker et al., 1990). Haloperidol, the S1R antagonist used in this study, is a first generation (typical), FDA-approved, antipsychotic (Shader and Shader, 1994). It mitigates positive symptoms of schizophrenia by blocking dopamine D2 receptors in the brain. However, Haloperidol is not selective for D2 receptors, and also has high affinity for S1Rs. As such it has been used extensively in cellular, molecular and biochemical studies of ligand-dependent S1R effects (Mei and Pasternak, 2007Bai et al., 2017).

Our results provide a unique and important contribution to knowledge of S1R biology. Results are notable for enrichment of proteins integral to secretion and extracellular matrix formation as well as cholesterol biosynthesis and metabolism (Figure 5). Consistent with the latter result, previous studies have shown that S1R interacts with cholesterol and have implicated S1R in lipid metabolism and transport (Palmer et al., 2007Hayashi and Su, 2010Zhemkov et al., 2021b). However, the mechanisms that underlie these interactions are not well understood.

Our results indicate that under conditions of cellular exposure to Haloperidol, the proximal association between S1R-Apex and Proprotein convertase subtilisin/kexin type 9 (PCSK9) was increased. In addition, exposure to (+)-PTZ also increased the proximal association between S1R and PCSK9 and between S1R and the low density lipoprotein receptor (LDLR) (Figure 5). PCSK9 and LDLR are known to interact with each other and both are medically important because mutations in both genes are associated with familial hypercholesterolemia (Abifadel et al., 2003Abifadel et al., 2009aAbifadel et al., 2009bAbifadel et al., 2012). LDLR, situated on cellular membranes, binds to and mediates endocytosis of cholesterol-rich LDL from the extracellular space. The endocytosed LDL is eventually broken down and cholesterol is released for use by the cell. Interaction between PCSK9 and LDLR has been shown to lead to lysosome-mediated destruction of LDLR and therefore to decreased levels of LDLR at the cellular surface (Brown and Goldstein, 1986). In general, a decrease in cell surface LDLR levels increases concentration of systemically circulating LDL-cholesterol (Brown and Goldstein, 1986). Therefore, inhibition of PCSK9 is clinically useful for treatment of hypercholesterolemia, and inhibitors of PCSK9 have been FDA-approved as cholesterol-lowering therapeutics (Weinreich and Frishman, 2014).

Our findings indicate that Haloperidol treatment not only increases the proximal association between S1R-Apex and PCSK9, but also leads to increased levels of intracellular and extracellular (secreted) PCSK9 (Figures 67). Somewhat surprisingly, LDLR levels were unaffected by Haloperidol treatment, at least under the conditions of our experimental setup (Figure 5D). This finding was unexpected given the documented role of PCSK9 in mediating turnover of LDLR. However, the organ that is mostly responsible for producing and secreting PCSK9 is the liver. Secreted PCSK9 can then alter the level of LDLR in distal tissues, a scenario that is different from the cell culture setup used in our experiments. Another possibility is that treatment with Haloperidol for an extended duration might be required to observe a PCSK9-induced reduction in LDLR levels. Although Haloperidol treatment did not affect the level of LDLR, the intracellular localization of LDLR was dramatically altered. Instead of mostly diffuse cytoplasmic staining, LDLR was localized to large puncta in Haloperidol treated cells (Figure 8). At present, the mechanism that causes this relocalization of intracellular LDLR is unknown.

In addition to exploring effects of Haloperidol, a known S1R antagonist, on the S1R proximatome, we also evaluated the effects of the high affinity S1R agonist, (+)-Pentazocine. Most of the proximatome changes induced by (+)-PTZ correspond to secreted and extracellular matrix proteins such as Thrombospondin1 (THBS1), Collagen12 (COL12A1) and Fibulin 1 and 3 (EFEMP1 and FBN1 respectively) (Figure 5). For this set, we validated the biotinylation between S1R-Apex and Thrombospondin1 (Figure 5E). In addition, proteins involved in cholesterol metabolism and homeostasis such as squalene synthase (FDFT1), PCKS9, and LDLR were also enriched in this dataset. Consistent with our proteomics results, (+)-PTZ treatment resulted in increased biotinylation of LDLR by S1R-Apex (Figure 5D). The intracellular localization of LDLR was also altered upon (+)-PTZ treatment, but to a lesser degree than when cells were treated with Haloperidol (Figure 7).

One common theme between (+)-PTZ and Haloperidol treatment appears to involve cholesterol metabolism. Both ligands affect the localization of intracellular LDLR and treatment with both compounds promoted the cellular uptake of labeled LDL (Figure 7). Consistent with a role for S1R in cholesterol metabolism, treatment with Haloperidol increased the total cellular level of cholesterol (Figure 8H). A similar trend was observed in (+)-PTZ treated cells. However, the difference was not statistically significant. The effects of S1R ligand binding on cholesterol metabolism is medically important because medications with S1R affinity are already in widespread use (Gordon et al., 2020Ye et al., 2020Maurice, 2021Reis et al., 2022).

Another common theme that emerges from our studies is the involvement of S1R in modulation of protein secretion. Proteins destined for secretion are known to move from ER exit sites through the ER-Golgi intermediate compartment (ERGIC) (Appenzeller-Herzog and Hauri, 2006). As shown in Figure 2B, components of the ERGIC are enriched in the S1R proximatome. This finding was validated by Western blot analysis of pellets using an antibody against Lman1, a marker protein of the ERGIC, and a protein that was highly enriched in the S1R proximatome (Figure 4A and Supplementary Table S1) (Hauri et al., 2000Fu et al., 2019). Furthermore, a globlal secretome analysis identified a total of 1,087 proteins in HeLa cell conditioned media. Of these, 78 proteins were also detected in the S1R proximatome (Figure 7CSupplementary Table S7). Interestingly, the ligand-bound state of S1R correlated with increased or decreased secretion of a subset of these proteins (Figure 7). Our results are therefore consistent with previous studies that implicate S1R as an important regulator of the cellular secretome. For example, stimulation of S1R with (+)-PTZ leads to increased astrocytic release of brain-derived neurotrophic factor (BDNF), a neurotrophin that supports neuronal growth and survival (Mysona et al., 2018). In addition, S1R has been shown to regulate levels of several other proteins that are processed and transported through the cellular secretory pathway (Maher et al., 2018Lopez et al., 2019).

Collectively, our findings indicate that the S1R-Apex proximatome contains a large subset of proteins. Somewhat unexpectedly, the S1R-Apex proximatome was not significantly altered upon treating cells with either a classical antagonist or agonist. It should be noted, however, that our analysis only considered proteins showing a greater than two-fold, ligand-dependent biotinylation change as significant. In an organismal setting, proteins displaying more subtle changes upon ligand binding might nevertheless elicit physiological effects. Furthermore, the time scale of ligand incubation might affect the degree and nature of proximatome changes. Our proximatome experiments were performed after 24 h of ligand treatment. Additional time points might have revealed different changes. A final consideration is that the S1R proximatome might be different in different cell types and tissues. Despite these caveats, S1R agonists and antagonists are under consideration for an ever-widening spectrum of pathologies ranging from COVID-19 treatment to cancer diagnosis, chronic pain remedies, and neurodegenerative disease therapeutics (Happy et al., 2015Das et al., 2016Vavers et al., 2019Vela, 2020). Our results are therefore critically relevant to a broad range of translational outcomes.