Cytosolic N-GlcNAc proteins are formed by the action of endo-\xce\xb2-N-acetylglucosaminidase
NGLY1 is a widely conserved eukaryotic cytosolic deglycosylating enzyme involved in the endoplasmic reticulum-associated degradation (ERAD) process, which eliminates misfolded proteins through retrograde translocation and proteasomal degradation. A human genetic disorder called NGLY1-deficiency has been reported, indicating the functional importance of NGLY1 in humans. Evidence suggests that Ngly1-KO is embryonic lethal in mice, while additional deletion of the Engase gene, encoding another cytosolic deglycosylating enzyme (endo-\xce\xb2-N-acetylglucosaminidase; ENGase), partially rescued lethality. Upon compromised Ngly1 activity, ENGase-mediated deglycosylation of misfolded glycoproteins may cause excess formation of N-GlcNAc proteins in the cytosol, leading to detrimental effects in the mice. Whether endogenous N-GlcNAc proteins are really formed in Ngly1-KO cells/animals or not remains unclarified. Here, comprehensive identification of O- and N-GlcNAc proteins was carried out using purified cytosol from wild type, Ngly1-KO, Engase-KO, and Ngly1/Engase double KO mouse embryonic fibroblasts. It was revealed that while there is no dramatic change in the level of O-GlcNAc proteins among cells examined, there was a vast increase of N-GlcNAc proteins in Ngly1-KO cells upon proteasome inhibition. Importantly, few N-GlcNAc proteins were observed in Engase-KO or Ngly1/Engase double-KO cells, clearly indicating that the cytosolic ENGase is responsible for the formation of N-GlcNAc proteins. The excess formation of N-GlcNAc proteins may at least in part account for the pathogenesis of NGLY1-deficiency.
INTRODUCTION<!>Cell Culture Procedures<!>Cytosol Enrichment<!>Glycopeptide Enrichment<!>Ubiquitin-Peptide Enrichment<!>Mass Spectrometry<!>O-GlcNAc Network Analysis<!>Formation of N-GlcNAc Proteins in Mouse Embryonic Fibroblasts<!>O-GlcNAcylation of ubiquitin-associated proteins<!>General Discussion<!>
<p>Endoplasmic reticulum (ER)1-associated degradation (ERAD) constitutes one of the protein homeostasis/quality control systems for newly synthesized proteins in the ER [1]. It has been proposed that N-glycosylated proteins have dedicated ERAD machinery [2]. Misfolded N-glycoproteins are recognized by the ERAD system and are retrotranslocated into the cytosol, where the proteasome plays a critical role in their degradation. During degradation, N-glycans are released by cytosolic peptide:N-glycanase (PNGase) (Fig. 1A) [3].</p><p>Cytosolic PNGase was first identified and reported in mammalian cells [4,5]. The gene encoding the cytosolic PNGase (Ngly1 in mammals [6]) is well conserved in eukaryotes [7]. Recently, patients bearing mutations in the NGLY1 gene have been described and the clinical disorder termed NGLY1-deficiency [8–10]. Patients with NGLY1-deficiency display systemic symptoms, including developmental delay, movement disorders, peripheral neuropathy, hypo/alacrima, and hypotonia [9,10], clearly indicating the functional importance of NGLY1.</p><p>Previously it was revealed that Ngly1-KO in C57BL/6 mice is embryonic/perinatal lethal [11]. Unexpectedly, the embryonic lethality of Ngly1-KO mice was partially rescued by the additional deletion of the Engase gene [11]. Engase encodes a cytosolic endo-β-N-acetylglucosaminidase (ENGase), another cytosolic N-glycan deglycosylating enzyme [12]. ENGase is involved in the catabolism of free N-glycans generated by multiple pathways (Fig. 1A) [13]. This observation led to the hypothesis that ENGase may be a drug target candidate for NGLY1-deficiency. We previously proposed the "N-GlcNAc hypothesis" as a potential mechanism for how ENGase can partially rescue the defects caused by Ngly1 deletion [11,14]. In this proposed process, ENGase cleaves glycans from cytosolic misfolded glycoproteins in a stochastic fashion, and N-GlcNAc modified proteins would be generated in the cytosol. In fact, degradation of a model ERAD substrate (RTAΔm derived from ricin toxin A chain) was delayed in Ngly1-KO mouse embryonic fibroblasts (MEF) cells, and the substrate was still deglycosylated by ENGase, forming N-GlcNAc RTAΔm [14]. Moreover, the N-GlcNAc-modified RTAΔm was found to form detergent-insoluble aggregates, implying that some N-GlcNAc proteins may be prone to form proteasome-resistant aggregates. These observations suggest that upon compromised NGLY1 activity, excess formation of N-GlcNAc proteins is expected, and the formation of excessive levels of N-GlcNAc proteins could lead to detrimental effects in part accounting for the pathophysiology of NGLY1-deficiency. The N-GlcNAc hypothesis also proposes that the presence of N-GlcNAc proteins may affect pre-existing O-GlcNAc mediated intracellular signaling.</p><p>While N-GlcNAc-containing peptides/proteins have been identified in nature from various sources [15–20], how NGLY1 and ENGase are involved in their formation remains unclear. In this study, a comprehensive analysis to detect and identify cytosolic N-GlcNAc and O-GlcNAc glycopeptides in wild type, Ngly1-KO, Engase-KO and Ngly1/Engase double KO MEF cells was performed to provide insight into the roles of Ngly1 and ENGase in the formation of N-GlcNAc proteins.</p><!><p>MEF cells from all genetic backgrounds were prepared and cultivated as described previously [14].</p><!><p>Two million MEF cells were seeded on a 150 cm2 flask 48 h before extraction. For MG132 treatment, the MEF cells were cultured with medium containing 5 μM of MG132 for 3 h. Detergent extraction was used to separate cytosolic contents [21]. Briefly, MEF cells were washed 3 times with 10 mL of PBS. Cytosol was released with 1 mL of permeabilization buffer (25 mM KHEPES (pH 7.2), 110 mM potassium acetate, 2.5 mM magnesium acetate, 1 mM ethyleneglycoltetraacetic acid (EGTA), 1 mM dithiothreitol (DTT), 1 mM phenylmethylsulfonyl fluoride (PMSF), 100 μM O-(2-Acetamido-2-deoxy-D-glucopyranosylidenamino) N-phenylcarbamate (PUGNAc), 0.015% (w/v) digitonin) while slowly rocked at 4°C for 5 min. The collected buffer was then centrifuged at 7,500×g for 10 min at 4°C. The supernatant was used as the cytosolic fraction for glycoproteomic analysis.</p><!><p>Glycopeptide enrichment by lectin weak affinity chromatography (LWAC) was described previously [16,19]. Briefly, 6-11 mg of protein from MEF cytosol were denatured in 8 M urea, reduced with 2.5 mM DTT for 1 hour at 25°C, and cysteines were alkylated with 5 mM iodoacetic acid at 25°C in the dark for 45 mins. Samples were diluted to 2 M urea with 25 mM KHEPES (pH 7.2) and digested overnight at 37°C with 5% (w/w) sequencing-grade trypsin (Promega). Digested samples were acidified with formic acid (Sigma-Aldrich) and subsequently desalted using a 360 mg C18 Sep-Pak SPE cartridge (Waters). Desalted samples were dried using a SpeedVac concentrator (Thermo). Desalted cytosolic peptides were resuspended in LWAC buffer (100 mM Tris-HCl (pH 7.5), 150 mM NaCl, 10 mm MgCl2, 10 mm CaCl2, and 5% acetonitrile). Glycopeptides were enriched with a POROS- wheat germ agglutinin (WGA) column, collected, and desalted inline using a Luna 10μ C18 column (Phenomenex). A total of three rounds of LWAC enrichment were performed. LWAC-enriched glycopeptides were subsequently fractionated by high-pH reversed-phase liquid chromatography.</p><!><p>Peptides from LWAC flowthrough were desalted as above. Ubiquitin modified peptides were enriched using the PTMScan ubiquitin remnant motif (K;-ε-GG) kit (Cell Signaling Technology) according to the protocol from Udeshi et al. [22].</p><!><p>Glycopeptide enriched samples were analyzed on an Orbitrap Velos (Thermo Fisher Scientific) equipped with a NanoAcquity UPLC (Waters). Peptides were fractionated on a 15 cm × 75 μm ID 3 μm C18 EASY-Spray column using a linear gradient from 2 to 30% solvent B over 60 min. Survey mass measurements were performed using the Orbitrap, scanning from m/z 350–2000. The three most abundant multiply charged ions were selected for higher energy collisional dissociation (HCD) and electron transfer dissociation (ETD) analysis. The trigger intensity was set to 2000. Supplemental activation was enabled. The ETD fragments were measured in the linear trap; HCD fragments were measured in the Orbitrap. Each sample was injected twice; the first analysis selected only 2+ precursor ions, and in the second analysis, 2+ precursor ions were excluded. Peaklists were extracted using Proteome Discoverer 1.4. The data were searched against the SwissProt Mus musculus and Bos taurus databases (22941 entries, downloaded January 1, 2017) (and concatenated with a randomized sequence for each entry) using Protein Prospector (version 5.23.0). Cleavage specificity was set as tryptic, allowing for two missed cleavages. Carboxymethylation of Cys was set as a constant modification. The required mass accuracy was 10 ppm for precursor ions and fragments ion mass accuracy was 0.6 Da for ETD files or 30 ppm for HCD files. Variable modifications are in Supplemental Table 8. Three modifications per peptide were permitted. Modified peptides were identified with a protein and peptide false discovery rate of 1%. HexNAc modified peptides were manually verified and O-GlcNAc and O-GalNAc modifications were differentiated based on protein localization and HexNAc oxonium ion fragment ratios [23]. Results can be viewed in MS Viewer (ETD search key: beypz6zf20, HCD search key: mriemlhilv) [24]. When comparing across cell conditions, peptide spectral matches (PSM) were normalized using total spectra.</p><p>Ubiquitin enriched samples were analyzed on a Q-Exactive Plus (Thermo Fisher Scientific) equipped with a NanoAcquity UPLC (Waters). Peptides were fractionated as above using a gradient from 2-25% B over 78 mins, then to 37% B over 4 mins, and to 40% B over 3 mins before washing and equilibrating the column. Survey mass measurements were scanned from m/z 350-1500. The ten most abundant multiply charged ions were selected with dynamic exclusion for HCD fragmentation with 25 normalized collision energy. Protein prospector search parameters were as above also allowing for GlyGly on uncleaved lysine as a variable modification. Modified peptides were identified with a protein and peptide false discovery rate of 1%. Results can be viewed in MS Viewer (search key: yxtxwkxdg0).</p><p>Raw mass spectrometry data has been deposited to MassIVE with the dataset identifier: MSV000084330.</p><!><p>Network analysis was performed in Cytoscape v3.7.1 with stringApp v1.4.2 and Omics Visualizer v1.1.2. A confidence score cutoff of 0.9 was set for the network analysis and the functional enrichment analysis used an FDR value cutoff of 0.05.title.</p><!><p>Cytosolic proteins were isolated from MEF cells [11] and digested with trypsin. Tryptic glycopeptides were enriched using WGA-LWAC and analyzed the by mass spectrometry (Fig. 1B). HCD and ETD fragmentation techniques were used for glycopeptide identification to ensure spectra identified as glycopeptides do in fact contain a glycan and to properly assign the type and site of the glycan modification.</p><p>The vast majority (90%) of assigned glycopeptide spectral matches (gPSMs) in wild type (WT) MEF cells were O-GlcNAc-modified (Supplemental Tables 1–4). However, glycopeptides from non-nucleocytoplasmic proteins were also identified and included O-GalNAc and N-linked glycans. The pool of N-linked glycans included N-GlcNAc modified peptides. The identification of N-GlcNAc modified peptides was used to assess ENGase activity in MEF cells.</p><p>The number of N-GlcNAc spectral matches in the cytosol of both WT and Ngly1-KO (NKO) MEF cells was similar, i.e. 1-2% of identified gPSMs (Fig. 2A, Supplemental Table 1)Over. Over 200 N-GlcNAc-containing peptides from almost 100 different proteins were identified in the NKO cells upon treatment with the proteosomal inhibitor. Ninety percent of these identified sites are annotated in Uniprot as known N-glycosylation sites. Figure 2B highlights an example of an N-GlcNAc peptide from thrombospondin-1, a secreted glycoprotein that mediates cell-cell and cell-matrix interactions and also plays a role in ER stress response [25]. The GlcNAc modification is present on Asn-1067, a known N-glycosylation site [26]. The c3 fragment ion confidently places the GlcNAc on Asn-1067, and the HCD fragmentation pattern (inset) is consistent with GlcNAc rather than GalNAc [23]. To determine if these peptides could originate from ERAD substrates, the LWAC flowthrough was analyzed for ubiquitin modifications. Tryptic digestion of ubiquitinated proteins results in Lys residues modified with two glycines (K-ε-GG). Peptides bearing this modification can be enriched with an antibody directed against the K-ε-GG motif [27]. More than 40 percent of the proteins with identified N-GlcNAc moieties in the MG132 treated NKO sample were found to also contain K-ε-GG moieties, establishing that these proteins represent in vivo ERAD substrates (Supplemental Table 5). Figure 2C illustrates an example of a thrombospondin-1 peptide containing a K-ε-GG modification. Three additional ubiquitination sites were identified on thrombospondin-1 from the MG132 treated NKO cytosol (Supplemental Table 6).</p><p>The presence of N-GlcNAc moieties in the Ngly1-KO samples suggests that the cytoplasmic ENGase processes the protein-bound N-linked glycan structures of ERAD substrates in the cytosol. To verify this, cytosolic glycopeptides from Engase-KO (EKO) or Ngly1/Engase double KO (DKO) MEF cells were enriched and identified. Only a few N-GlcNAc peptides were identified in EKO or DKO cells irrespective of MG132 treatment (Fig. 2A, Supplemental Table 1). These data provide strong evidence that the cytosolic ENGase generates N-GlcNAc moieties.</p><p>It was interesting to note that, in the absence of proteasome inhibitor, there was little difference in the number of detected N-GlcNAc proteins between Ngly1-KO and WT MEF cells. This result indicates that N-GlcNAc formation by ENGase can occur under normal conditions. Upon treatment with proteasome inhibitor, however, there were fewer N-GlcNAc proteins identified in WT MEF, suggesting that the level of cytosolic N-GlcNAc proteins may be decreased in the presence of proteasome inhibitor. If N-GlcNAc modifications lead to protein aggregation, it is possible that cytosol enrichment with digitonin is not sensitive to the detection of those aggregated proteins. Cytosolic aggregates are reported to be segregated to mitochondria for degradation as a part of the cytosolic proteostasis machinery [28] and it is possible that mitochondria may be involved in the clearance of excess N-GlcNAc proteins. Mitochondrial function appears compromised in Ngly1-KO cells [29,30], and Ngly1-KO cells/animals are sensitive to proteasome inhibitors [31–34]. Those combination effects may be the reason for differential response towards proteasome inhibitor treatment between NKO and WT cells.</p><!><p>LWAC-based glycopeptide enrichment was originally developed to enrich O-GlcNAc modified peptides [15,16,35]. Unsurprisingly, nearly 90% of cytosolic-glycopeptide spectral matches across all cell conditions were O-GlcNAc-containing peptides (Fig. 3A). Across all samples, over 3000 unique O-GlcNAc modified peptides were identified from 466 proteins and nearly 700 unambiguous sites of O-GlcNAcylation were assigned. When normalized to total spectra, no large differences were apparent between cell conditions (Fig. 3B).</p><p>Network and enrichment analysis of the O-GlcNAc dataset identified clustered networks involved in many cellular processes including transcriptional regulation, mRNA processing, nuclear and intracellular trafficking, and clathrinmediated endocytosis (Fig. S1, Supplemental Table 7). One such network contains proteins involved in the protein ubiquitination conjugation pathway (Fig. 3C). Nearly 40 percent of the identified O-GlcNAc modified proteins are related to ubiquitin conjugation (Uniprot keywords KW-0832 and KW-0833). Many were identified regardless of MG132 treatment or genetic background (Supplemental Table 7). Within this subset of proteins, at least one ERAD-associated protein, F-box only protein 6 (FBXO6), was shown to be O-GlcNAcylated. FBXO6 is part of an SCF complex (Skp, Cullin, F-box containing complex) which recognizes unfolded glycoproteins and targets them for degradation [36,37]. Although we only found O-GlcNAc-modified FBXO6 in EKO and DKO samples (Supplemental Table 4), further experimentation and analysis should be undertaken to determine the cellular mechanism of this modification and how it relates to NGLY1 expression.</p><!><p>Recent studies suggest NGLY1 is involved not only in the degradation of misfolded glycoproteins, but also activation of a transcription factor, NFE2L1 [31,32]. NFE2L1 is known to be involved in various cellular stress responses, including the proteasome recovery pathway upon inhibition of proteasome activity [38]. Accordingly, NGLY1-KO cells/animals exhibited sensitivity towards proteasome inhibitors [31–34]. It has been shown that NFE2L1-dependent NGLY1 function is independent from ENGase, as DKO MEF are as sensitive as NKO MEF to proteasome inhibitors [39]. While the contribution of N-GlcNAc protein accumulation to the pathophysiology of NGLY1-deficiency remains unclear, it is important to note that we can conclude that almost all the cytosolic N-GlcNAc proteins are formed by the action of ENGase.</p><p>Part of the N-GlcNAc hypothesis proposes that the presence of N-GlcNAc proteins may affect pre-existing O-GlcNAc mediated signaling. Our analysis showed little change in the number of O-GlcNAc proteins. It remains to be investigated whether the amount or localization of a subset of specific O-GlcNAc proteins may be altered. Our analysis also showed little difference in the level of ubiquitinated proteins between WT and NKO cells, with or without proteasome inhibitortreatment. This result is consistent with the recent observation that increase in the level of polyubiquitinated proteins are not observed in Ngly1-KO liver tissue [39] and ubiquitin-positive aggregates are only seen in limited tissue in Ngly1-KO rat [40]. We can therefore safely assume that accumulation of ubiquitinated proteins may be distinct among different tissues and cells. Further studies will help clarify how the accumulation of N-GlcNAc proteins impact the pathophysiology of NGLY1 deficiency.</p><!><p>1. String network and enrichment analysis of O-GlcNAcylated proteins across all MEF cell lines and treatments. Many cellular processes including transcriptional regulation, mRNA processing (pre-mRNA processing, mRNA decay, and splicing), nuclear and intracellular trafficking, and clathrin mediated endocytosis are clustered within these proteins.</p>