Systematic Evaluation of Fragmentation Methods for Unlabeled and Isobaric Mass Tag-Labeled O-Glycopeptides

Quantitative proteome-wide O-glycoproteomics analysis with FragPipe

Identification of O-glycopeptides from tandem mass spectrometry data is complicated by the near complete dissociation of O-glycans from the peptide during collisional activation and by the combinatorial explosion of possible glycoforms when glycans are retained intact in electron-based activation. The recent O-Pair search method provides an elegant solution to these problems, using a collisional activation scan to identify the peptide sequence and total glycan mass, and a follow-up electron-based activation scan to localize the glycosite(s) using a graph-based algorithm in a reduced search space. Our previous O-glycoproteomics methods with MSFragger-Glyco allowed for extremely fast and sensitive identification of O-glycopeptides from collisional activation data but had limited support for site localization of glycans and quantification of glycopeptides. Here, we report an improved pipeline for O-glycoproteomics analysis that provides proteome-wide, site-specific, quantitative results by incorporating the O-Pair method as a module within FragPipe. In addition to improved search speed and sensitivity, we add flexible options for oxonium ion-based filtering of glycans and support for a variety of MS acquisition methods and provide a comparison between all software tools currently capable of O-glycosite localization in proteome-wide searches.

Hydrophilic Interaction Chromatography Coupled to Ultraviolet Photodissociation Affords Identification, Localization, and Relative Quantitation of Glycans on Intact Glycoproteins

Protein glycosylation is implicated in a wide array of diseases, yet glycoprotein analysis remains elusive owing to the extreme heterogeneity of glycans, including microheterogeneity of some of the glycosites (amino acid residues). Various mass spectrometry (MS) strategies have proven tremendously successful for localizing and identifying glycans, typically utilizing a bottom-up workflow in which glycoproteins are digested to create glycopeptides to facilitate analysis. An emerging alternative is top-down MS that aims to characterize intact glycoproteins to allow precise identification and localization of glycans. The most comprehensive characterization of intact glycoproteins requires integration of a suitable separation method and high performance tandem mass spectrometry to provide both protein sequence information and glycosite localization. Here, we couple ultraviolet photodissociation and hydrophilic interaction chromatography with high resolution mass spectrometry to advance the characterization of intact glycoproteins ranging from 15 to 34 kDa, offering site localization of glycans, providing sequence coverages up to 93%, and affording relative quantitation of individual glycoforms.

Targeting host O-linked glycan biosynthesis affects Ebola virus replication efficiency and reveals differential GalNAc-T acceptor site preferences on the Ebola virus glycoprotein

Ebola virus glycoprotein (EBOV GP) is one of the most heavily O-glycosylated viral glycoproteins, yet we still lack a fundamental understanding of the structure of its large O-glycosylated mucin-like domain and to what degree the host O-glycosylation capacity influences EBOV replication. Using tandem mass spectrometry, we identified 47 O-glycosites on EBOV GP and found similar glycosylation signatures on virus-like particle- and cell lysate-derived GP. Furthermore, we performed quantitative differential O-glycoproteomics on proteins produced in wild-type HEK293 cells and cell lines ablated for the three key initiators of O-linked glycosylation, GalNAc-T1, -T2, and -T3. The data show that 12 out of the 47 O-glycosylated sites were regulated, predominantly by GalNAc-T1. Using the glycoengineered cell lines for authentic EBOV propagation, we demonstrate the importance of O-linked glycan initiation and elongation for the production of viral particles and the titers of progeny virus. The mapped O-glycan positions and structures allowed to generate molecular dynamics simulations probing the largely unknown spatial arrangements of the mucin-like domain. The data highlight targeting GALNT1 or C1GALT1C1 as a possible way to modulate O-glycan density on EBOV GP for novel vaccine designs and tailored intervention approaches.IMPORTANCEEbola virus glycoprotein acquires its extensive glycan shield in the host cell, where it is decorated with N-linked glycans and mucin-type O-linked glycans. The latter is initiated by a family of polypeptide GalNAc-transferases that have different preferences for optimal peptide substrates resulting in a spectrum of both very selective and redundant substrates for each isoform. In this work, we map the exact locations of O-glycans on Ebola virus glycoprotein and identify subsets of sites preferentially initiated by one of the three key isoforms of GalNAc-Ts, demonstrating that each enzyme contributes to the glycan shield integrity. We further show that altering host O-glycosylation capacity has detrimental effects on Ebola virus replication, with both isoform-specific initiation and elongation playing a role. The combined structural and functional data highlight glycoengineered cell lines as useful tools for investigating molecular mechanisms imposed by specific glycans and for steering the immune responses in future vaccine designs.

ETD-Based Proteomic Profiling Improves Arginine Methylation Identification and Reveals Novel PRMT5 Substrates

Protein arginine methylations are important post-translational modifications (PTMs) in eukaryotes, regulating many biological processes. However, traditional collision-based mass spectrometry methods inevitably cause neutral losses of methylarginines, preventing the deep mining of biologically important sites. Herein we developed an optimized mass spectrometry workflow based on electron-transfer dissociation (ETD) with supplemental activation for proteomic profiling of arginine methylation in human cells. Using symmetric dimethylarginine (sDMA) as an example, we show that the ETD-based optimized workflow significantly improved the identification and site localization of sDMA. Quantitative proteomics identified 138 novel sDMA sites as potential PRMT5 substrates in HeLa cells. Further biochemical studies on SERBP1, a newly identified PRMT5 substrate, confirmed the coexistence of sDMA and asymmetric dimethylarginine in the central RGG/RG motif, and loss of either methylation caused increased the recruitment of SERBP1 to stress granules under oxidative stress. Overall, our optimized workflow not only enabled the identification and localization of extensive, nonoverlapping sDMA sites in human cells but also revealed novel PRMT5 substrates whose sDMA may play potentially important biological functions.

Variations in O-Glycosylation Patterns Influence Viral Pathogenicity, Infectivity, and Transmissibility in SARS-CoV-2 Variants

The highly glycosylated S protein plays a vital role in host cell invasion, making it the principal target for vaccine development. Differences in mutations observed on the spike (S) protein of SARS-CoV-2 variants may result in distinct glycosylation patterns, thus influencing immunological evasion, infectivity, and transmissibility. The glycans can mask key epitopes on the S1 protein and alter its structural conformation, allowing the virus to escape the immune system. Therefore, we comprehensively characterize O-glycosylation in eleven variants of SARS-CoV-2 S1 subunits to understand the differences observed in the biology of the variants. In-depth characterization was performed with a double digestion strategy and an efficient LC-MS/MS approach. We observed that O-glycosylation is highly conserved across all variants in the region between the NTD and RBD, whereas other domains and regions exhibit variation in O-glycosylation. Notably, omicron has the highest number of O-glycosylation sites on the S1 subunit. Also, omicron has the highest level of sialylation in the RBD and RBM functional motifs. Our findings may shed light on how differences in O-glycosylation impact viral pathogenicity in variants of SARS-CoV-2 and facilitate the development of a robust vaccine with high protective efficacy against the variants of concern.

In-Depth Mass Spectrometry Analysis Reveals the Plasma Proteomic and N-Glycoproteomic Impact of an Amish-Enriched Cardioprotective Variant in B4GALT1

B4GALT1 encodes β-1,4-galactosyltransferase 1, an enzyme that plays a major role in glycan synthesis in the Golgi apparatus by catalyzing the addition of terminal galactose. Studies increasingly suggest that B4GALT1 may be involved in the regulation of lipid metabolism pathways. Recently, we discovered a single-site missense variant Asn352Ser (N352S) in the functional domain of B4GALT1 in an Amish population, which decreases the level of LDL-cholesterol (LDL-c) as well as the protein levels of ApoB, fibrinogen, and IgG in the blood. To systematically evaluate the effects of this missense variant on protein glycosylation, expression, and secretion, we developed a nano-LC-MS/MS-based platform combined with TMT-labeling for in-depth quantitative proteomic and glycoproteomic analyses in the plasma of individuals homozygous for the B4GALT1 missense variant N352S versus non-carriers (n = 5 per genotype). A total of 488 secreted proteins in the plasma were identified and quantified, 34 of which showed significant fold changes in protein levels between N352S homozygotes and non-carriers. We determined N-glycosylation profiles from 370 glycosylation sites in 151 glycoproteins and identified ten proteins most significantly associated with decreased galactosylation and sialyation in B4GALT1 N352S homozygotes. These results further support that B4GALT1 N352S alters the glycosylation profiles of a variety of critical target proteins, thus governing the functions of these proteins in multiple pathways, such as those involved in lipid metabolism, coagulation, and the immune response.

Global mapping of GalNAc-T isoform-specificities and O-glycosylation site-occupancy in a tissue-forming human cell line

Mucin-type-O-glycosylation on proteins is integrally involved in human health and disease and is coordinated by an enzyme family of 20 N-acetylgalactosaminyltransferases (GalNAc-Ts). Detailed knowledge on the biological effects of site-specific O-glycosylation is limited due to lack of information on specific glycosylation enzyme activities and O-glycosylation site-occupancies. Here we present a systematic analysis of the isoform-specific targets of all GalNAc-Ts expressed within a tissue-forming human skin cell line, and demonstrate biologically significant effects of O-glycan initiation on epithelial formation. We find over 300 unique glycosylation sites across a diverse set of proteins specifically regulated by one of the GalNAc-T isoforms, consistent with their impact on the tissue phenotypes. Notably, we discover a high variability in the O-glycosylation site-occupancy of 70 glycosylated regions of secreted proteins. These findings revisit the relevance of individual O-glycosylation sites in the proteome, and provide an approach to establish which sites drive biological functions.

Mass spectrometric analysis of chondroitin sulfate-linked peptides

Chondroitin sulfate proteoglycans (CSPGs) are extracellular matrix components composed of linear glycosaminoglycan (GAG) side chains attached to a core protein. CSPGs play a vital role in neurodevelopment, signal transduction, cellular proliferation and differentiation and tumor metastasis through interaction with growth factors and signaling proteins. These pleiotropic functions of proteoglycans are regulated spatiotemporally by the GAG chains attached to the core protein. There are over 70 chondroitin sulfate-linked proteoglycans reported in cells, cerebrospinal fluid and urine. A core glycan linker of 3-6 monosaccharides attached to specific serine residues can be extended by 20-200 disaccharide repeating units making intact CSPGs very large and impractical to analyze. The current paradigm of CSPG analysis involves digesting the GAG chains by chondroitinase enzymes and analyzing either the protein part, the disaccharide repeats, or both by mass spectrometry. This method, however, provides no information about the site of attachment or the composition of linker oligosaccharides and the degree of sulfation and/or phosphorylation. Further, the analysis by mass spectrometry and subsequent identification of novel CSPGs is hampered by technical challenges in their isolation, less optimal ionization and data analysis. Unknown identity of the linker oligosaccharide also makes it more difficult to identify the glycan composition using database searching approaches. Following chondroitinase digestion of long GAG chains linked to tryptic peptides, we identified intact GAG-linked peptides in clinically relevant samples including plasma, urine and dermal fibroblasts. These intact glycopeptides including their core linker glycans were identified by mass spectrometry using optimized stepped higher energy collision dissociation and electron-transfer/higher energy collision dissociation combined with hybrid database search/de novo glycan composition search. We identified 25 CSPGs including three novel CSPGs that have not been described earlier. Our findings demonstrate the utility of combining enrichment strategies and optimized high-resolution mass spectrometry analysis including alternative fragmentation methods for the characterization of CSPGs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42485-022-00092-3.