Distinct shifts in site-specific glycosylation pattern of SARS-CoV-2 spike proteins associated with arising mutations in the D614G and Alpha variants

Enhanced Surface Accessibility of SARS-CoV-2 Omicron Spike Protein Due to an Altered Glycosylation Profile

SARS-CoV-2 spike (S) proteins undergo extensive glycosylation, aiding in proper folding, enhancing stability, and evading host immune surveillance. In this study, we used mass spectrometric analysis to elucidate the N-glycosylation characteristics and disulfide bonding of recombinant spike proteins derived from the SARS-CoV-2 Omicron variant (B.1.1.529) in comparison with the D614G spike variant. Furthermore, we conducted microsecond-long molecular dynamics simulations on spike proteins to resolve how the different N-glycans impact spike conformational sampling in the two variants. Our findings reveal that the Omicron spike protein maintains an overall resemblance to the D614G spike variant in terms of site-specific glycan processing and disulfide bond formation. Nonetheless, alterations in glycans were observed at certain N-glycosylation sites. These changes, in synergy with mutations within the Omicron spike protein, result in increased surface accessibility of the macromolecule, including the ectodomain, receptor-binding domain, and N-terminal domain. Additionally, mutagenesis and pull-down assays reveal the role of glycosylation of a specific sequon (N149); furthermore, the correlation of MD simulation and HDX-MS identified several high-dynamic areas of the spike proteins. These insights contribute to our understanding of the interplay between structure and function, thereby advancing effective vaccination and therapeutic strategies.

The role of N-glycosylation in spike antigenicity for the SARS-CoV-2 gamma variant

The emergence of SARS-CoV-2 variants alters the efficacy of existing immunity towards the viral spike protein, whether acquired from infection or vaccination. Mutations that impact N-glycosylation of spike may be particularly important in influencing antigenicity, but their consequences are difficult to predict. Here, we compare the glycosylation profiles and antigenicity of recombinant viral spike of ancestral Wu-1 and the Gamma strain, which has two additional N-glycosylation sites due to amino acid substitutions in the N-terminal domain (NTD). We found that a mutation at residue 20 from threonine to asparagine within the NTD caused the loss of NTD-specific antibody COVA2-17 binding. Glycan site-occupancy analyses revealed that the mutation resulted in N-glycosylation switching to the new sequon at N20 from the native N17 site. Site-specific glycosylation profiles demonstrated distinct glycoform differences between Wu-1, Gamma, and selected NTD variant spike proteins, but these did not affect antibody binding. Finally, we evaluated the specificity of spike proteins against convalescent COVID-19 sera and found reduced cross-reactivity against some mutants, but not Gamma spike compared to Wuhan spike. Our results illustrate the impact of viral divergence on spike glycosylation and SARS-CoV-2 antibody binding profiles.

Rapid simulation of glycoprotein structures by grafting and steric exclusion of glycan conformer libraries

Most membrane proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, shielding potentially large fractions of protein surface. High glycan conformational freedom hinders complete structural elucidation of glycoproteins. Computer simulations may be used to model glycosylated proteins but require hundreds of thousands of computing hours on supercomputers, thus limiting routine use. Here, we describe GlycoSHIELD, a reductionist method that can be implemented on personal computers to graft realistic ensembles of glycan conformers onto static protein structures in minutes. Using molecular dynamics simulation, small-angle X-ray scattering, cryoelectron microscopy, and mass spectrometry, we show that this open-access toolkit provides enhanced models of glycoprotein structures. Focusing on N-cadherin, human coronavirus spike proteins, and gamma-aminobutyric acid receptors, we show that GlycoSHIELD can shed light on the impact of glycans on the conformation and activity of complex glycoproteins.

Proteomics-based mass spectrometry profiling of SARS-CoV-2 infection from human nasopharyngeal samples

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the on-going global pandemic of coronavirus disease 2019 (COVID-19) that continues to pose a significant threat to public health worldwide. SARS-CoV-2 encodes four structural proteins namely membrane, nucleocapsid, spike, and envelope proteins that play essential roles in viral entry, fusion, and attachment to the host cell. Extensively glycosylated spike protein efficiently binds to the host angiotensin-converting enzyme 2 initiating viral entry and pathogenesis. Reverse transcriptase polymerase chain reaction on nasopharyngeal swab is the preferred method of sample collection and viral detection because it is a rapid, specific, and high-throughput technique. Alternate strategies such as proteomics and glycoproteomics-based mass spectrometry enable a more detailed and holistic view of the viral proteins and host-pathogen interactions and help in detection of potential disease markers. In this review, we highlight the use of mass spectrometry methods to profile the SARS-CoV-2 proteome from clinical nasopharyngeal swab samples. We also highlight the necessity for a comprehensive glycoproteomics mapping of SARS-CoV-2 from biological complex matrices to identify potential COVID-19 markers.

An N-glycopeptide MS/MS data analysis workflow leveraging two complementary glycoproteomic software tools for more confident identification and assignments

Complete coverage of all N-glycosylation sites on the SARS-CoV2 spike protein would require the use of multiple proteases in addition to trypsin. Subsequent identification of the resulting glycopeptides by searching against database often introduces assignment errors due to similar mass differences between different permutations of amino acids and glycosyl residues. By manually interpreting the individual MS2 spectra, we report here the common sources of errors in assignment, especially those introduced by the use of chymotrypsin. We show that by applying a stringent threshold of acceptance, erroneous assignment by the commonly used Byonic software can be controlled within 15%, which can be reduced further if only those also confidently identified by a different search engine, pGlyco3, were considered. A representative site-specific N-glycosylation pattern could be constructed based on quantifying only the overlapping subset of N-glycopeptides identified at higher confidence. Applying the two complimentary glycoproteomic software in a concerted data analysis workflow, we found and confirmed that glycosylation at several sites of an unstable Omicron spike protein differed significantly from those of the stable trimeric product of the parental D614G variant.

The Effect of Select SARS-CoV-2 N-Linked Glycan and Variant of Concern Spike Protein Mutations on C-Type Lectin-Receptor-Mediated Infection

The SARS-CoV-2 virion has shown remarkable resilience, capable of mutating to escape immune detection and re-establishing infectious capabilities despite new vaccine rollouts. Therefore, there is a critical need to identify relatively immutable epitopes on the SARS-CoV-2 virion that are resistant to future mutations the virus may accumulate. While hACE2 has been identified as the receptor that mediates SARS-CoV-2 susceptibility, it is only modestly expressed in lung tissue. C-type lectin receptors like DC-SIGN can act as attachment sites to enhance SARS-CoV-2 infection of cells with moderate or low hACE2 expression. We developed an easy-to-implement assay system that allows for the testing of SARS-CoV-2 trans-infection. Using our assay, we assessed how SARS-CoV-2 Spike S1-domain glycans and spike proteins from different strains affected the ability of pseudotyped lentivirions to undergo DC-SIGN-mediated trans-infection. Through our experiments with seven glycan point mutants, two glycan cluster mutants and four strains of SARS-CoV-2 spike, we found that glycans N17 and N122 appear to have significant roles in maintaining COVID-19's infectious capabilities. We further found that the virus cannot retain infectivity upon the loss of multiple glycosylation sites, and that Omicron BA.2 pseudovirions may have an increased ability to bind to other non-lectin receptor proteins on the surface of cells. Taken together, our work opens the door to the development of new therapeutics that can target overlooked epitopes of the SARS-CoV-2 virion to prevent C-type lectin-receptor-mediated trans-infection in lung tissue.

Polyvalent Nano-Lectin Potently Neutralizes SARS-CoV-2 by Targeting Glycans on the Viral Spike Protein

Mutations in spike (S) protein epitopes allow SARS-CoV-2 variants to evade antibody responses induced by infection and/or vaccination. In contrast, mutations in glycosylation sites across SARS-CoV-2 variants are very rare, making glycans a potential robust target for developing antivirals. However, this target has not been adequately exploited for SARS-CoV-2, mostly due to intrinsically weak monovalent protein-glycan interactions. We hypothesize that polyvalent nano-lectins with flexibly linked carbohydrate recognition domains (CRDs) can adjust their relative positions and bind multivalently to S protein glycans, potentially exerting potent antiviral activity. Herein, we displayed the CRDs of DC-SIGN, a dendritic cell lectin known to bind to diverse viruses, polyvalently onto 13 nm gold nanoparticles (named G13-CRD). G13-CRD bound strongly and specifically to target glycan-coated quantum dots with sub-nM K_d. Moreover, G13-CRD neutralized particles pseudotyped with the S proteins of Wuhan Hu-1, B.1, Delta variant and Omicron subvariant BA.1 with low nM EC₅₀. In contrast, natural tetrameric DC-SIGN and its G13 conjugate were ineffective. Further, G13-CRD potently inhibited authentic SARS-CoV-2 B.1 and BA.1, with <10 pM and <10 nM EC₅₀, respectively. These results identify G13-CRD as the 1st polyvalent nano-lectin with broad activity against SARS-CoV-2 variants that merits further exploration as a novel approach to antiviral therapy.

Analysis of the N-glycosylation profiles of the spike proteins from the Alpha, Beta, Gamma, and Delta variants of SARS-CoV-2

N-Glycosylation plays an important role in the structure and function of membrane and secreted proteins. Viral proteins used in cell entry are often extensively glycosylated to assist in protein folding, provide stability, and shield the virus from immune recognition by its host (described as a "glycan shield"). The SARS-CoV-2 spike protein (S) is a prime example, having 22 potential sites of N-glycosylation per protein protomer, as predicted from the primary sequence. In this report, we conducted mass spectrometric analysis of the N-glycosylation profiles of recombinant spike proteins derived from four common SARS-CoV-2 variants classified as Variant of Concern, including Alpha, Beta, Gamma, and Delta along with D614G variant spike as a control. Our data reveal that the amino acid substitutions and deletions between variants impact the abundance and type of glycans on glycosylation sites of the spike protein. Some of the N-glycosylation sequons in S show differences between SARS-CoV-2 variants in the distribution of glycan forms. In comparison with our previously reported site-specific glycan analysis on the S-D614G and its ancestral protein, glycan types on later variants showed high similarity on the site-specific glycan content to S-D614G. Additionally, we applied multiple digestion methods on each sample, and confirmed the results for individual glycosylation sites from different experiment conditions to improve the identification and quantification of glycopeptides. Detailed site-specific glycan analysis of a wide variety of SARS-CoV-2 variants provides useful information toward the understanding of the role of protein glycosylation on viral protein structure and function and development of effective vaccines and therapeutics.

Site specific N- and O-glycosylation mapping of the spike proteins of SARS-CoV-2 variants of concern

The glycosylation on the spike (S) protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, modulates the viral infection by altering conformational dynamics, receptor interaction and host immune responses. Several variants of concern (VOCs) of SARS-CoV-2 have evolved during the pandemic, and crucial mutations on the S protein of the virus have led to increased transmissibility and immune escape. In this study, we compare the site-specific glycosylation and overall glycomic profiles of the wild type Wuhan-Hu-1 strain (WT) S protein and five VOCs of SARS-CoV-2: Alpha, Beta, Gamma, Delta and Omicron. Interestingly, both N- and O-glycosylation sites on the S protein are highly conserved among the spike mutant variants, particularly at the sites on the receptor-binding domain (RBD). The conservation of glycosylation sites is noteworthy, as over 2 million SARS-CoV-2 S protein sequences have been reported with various amino acid mutations. Our detailed profiling of the glycosylation at each of the individual sites of the S protein across the variants revealed intriguing possible association of glycosylation pattern on the variants and their previously reported infectivity. While the sites are conserved, we observed changes in the N- and O-glycosylation profile across the variants. The newly emerged variants, which showed higher resistance to neutralizing antibodies and vaccines, displayed a decrease in the overall abundance of complex-type glycans with both fucosylation and sialylation and an increase in the oligomannose-type glycans across the sites. Among the variants, the glycosylation sites with significant changes in glycan profile were observed at both the N-terminal domain and RBD of S protein, with Omicron showing the highest deviation. The increase in oligomannose-type happens sequentially from Alpha through Delta. Interestingly, Omicron does not contain more oligomannose-type glycans compared to Delta but does contain more compared to the WT and other VOCs. O-glycosylation at the RBD showed lower occupancy in the VOCs in comparison to the WT. Our study on the sites and pattern of glycosylation on the SARS-CoV-2 S proteins across the VOCs may help to understand how the virus evolved to trick the host immune system. Our study also highlights how the SARS-CoV-2 virus has conserved both N- and O- glycosylation sites on the S protein of the most successful variants even after undergoing extensive mutations, suggesting a correlation between infectivity/ transmissibility and glycosylation.

Characterization and Function of Glycans on the Spike Proteins of SARS-CoV-2 Variants of Concern

SARS-CoV-2 variants of concern (VOCs) pose a great challenge to viral prevention and treatment owing to spike (S) protein mutations, which enhance their infectivity and capacity for immune evasion. However, whether these S protein mutations affect glycosylation patterns and thereby influence infectivity and immunogenicity remains unclear. In this study, four VOC S proteins-S-Alpha, S-Beta, S-Delta, and S-Omicron-were expressed and purified. Lectin microarrays were performed to characterize their glycosylation patterns. Several glycans were differentially expressed among the four VOC S proteins. Furthermore, the functional examination of glycans differentially expressed on S-Omicron revealed a higher expression of fucose-containing glycans, which modestly increased the binding of S-Omicron to angiotensin converting enzyme 2 (ACE2). A higher abundance of sialic acid and galactose-containing glycan was observed on S-Omicron, which significantly reduced its sensitivity against broad S protein-neutralizing antibodies. These findings contribute to the further understanding of SARS-CoV-2 infection mechanisms and provide novel glycan targets for emerging and future variants of SARS-CoV-2. IMPORTANCE Though glycosylation sites of SARS-CoV-2 S protein remain highly conserved, we confirmed that mutations in the Spike gene affect the S protein glycan expression pattern in different variants. More importantly, we found that glycans were differentially expressed on the S protein of the Omicron variant, enabling different forms of receptor binding and neutralization resistance. This study improves our understanding of SARS-CoV-2 glycomics and glycobiology and provides novel therapeutic and preventive strategies for SARS-CoV-2 VOCs.

Site Specific N- and O-glycosylation mapping of the Spike Proteins of SARS-CoV-2 Variants of Concern.. [PMID: 36415454 PMCID: PMC9681045 DOI: 10.21203/rs.3.rs-2188138/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The glycosylation on the spike (S) protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, modulates the viral infection by altering conformational dynamics, receptor interaction and host immune responses. Several variants of concern (VOCs) of SARS-CoV-2 have evolved during the pandemic, and crucial mutations on the S protein of the virus led to increased transmissibility and immune escape. In this study, we compare the site-specific glycosylation and overall glycomic profile of the wild type Wuhan-Hu-1 strain (WT) S protein and five VOCs of SARS-CoV-2: Alpha, Beta, Gamma, Delta and Omicron. Interestingly, both N- and O-glycosylation sites on the S protein are highly conserved among the spike mutant variants, particularly at the sites on the receptor-binding domain (RBD). The conservation of glycosylation sites is noteworthy, as over 2 million SARS-CoV-2 S protein sequences have been reported with various amino acid mutations. Our detailed profiling of the glycosylation at each of the individual sites of the S protein across the variants revealed intriguing possible association of glycosylation pattern on the variants and their previously reported infectivity. While the sites are conserved, we observed changes in the N- and O-glycosylation profile across the variants. The newly emerged variants, which showed higher resistance to neutralizing antibodies and vaccines, displayed a decrease in the overall abundance of complex-type glycans with both fucosylation and sialylation and an increase in the oligomannose-type glycans across the sites. Among the variants, the glycosylation sites with significant changes in glycan profile were observed at both the N-terminal domain (NTD) and RBD of S protein, with Omicron showing the highest deviation. The increase in oligomannose-type happens sequentially from Alpha through Delta. Interestingly, Omicron does not contain more oligomannose-type glycans compared to Delta but does contain more compared to the WT and other VOCs. O-glycosylation at the RBD showed lower occupancy in the VOCs in comparison to the WT. Our study on the sites and pattern of glycosylation on the SARS-CoV-2 S proteins across the VOCs may help to understand how the virus evolved to trick the host immune system. Our study also highlights how the SARS-CoV-2 virus has conserved both N- and O- glycosylation sites on the S protein of the most successful variants even after undergoing extensive mutations, suggesting a correlation between infectivity/transmissibility and glycosylation.

Core fucose-specific Pholiota squarrosa lectin (PhoSL) as a potent broad-spectrum inhibitor of SARS-CoV-2 infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein (S protein) is highly N-glycosylated, and a "glycan shield" is formed to limit the access of other molecules; however, a small open area coincides with the interface to the host's receptor and also neutralising antibodies. Most of the variants of concern have mutations in this area, which could reduce the efficacy of existing antibodies. In contrast, N-glycosylation sites are relatively invariant, and some are essential for infection. Here, we observed that the S proteins of the ancestral (Wuhan) and Omicron strains bind with Pholiota squarrosa lectin (PhoSL), a 40-amino-acid chemically synthesised peptide specific to core-fucosylated N-glycans. The affinities were at a low nanomolar level, which were ~ 1000-fold stronger than those between PhoSL and the core-fucosylated N-glycans at the micromolar level. We demonstrated that PhoSL inhibited infection by both strains at similar submicromolar levels, suggesting its broad-spectrum effect on SARS-CoV-2 variants. Cryogenic electron microscopy revealed that PhoSL caused an aggregation of the S protein, which was likely due to the multivalence of both the trimeric PhoSL and S protein. This characteristic is likely relevant to the inhibitory mechanism. Structural modelling of the PhoSL-S protein complex indicated that PhoSL was in contact with the amino acids of the S protein, which explains the enhanced affinity with S protein and also indicates the significant potential for developing specific binders by the engineering of PhoSL.

In situ structure and dynamics of an alphacoronavirus spike protein by cryo-ET and cryo-EM

Porcine epidemic diarrhea (PED) is a highly contagious swine disease caused by porcine epidemic diarrhea virus (PEDV). PED causes enteric disorders with an exceptionally high fatality in neonates, bringing substantial economic losses in the pork industry. The trimeric spike (S) glycoprotein of PEDV is responsible for virus-host recognition, membrane fusion, and is the main target for vaccine development and antigenic analysis. The atomic structures of the recombinant PEDV S proteins of two different strains have been reported, but they reveal distinct N-terminal domain 0 (D0) architectures that may correspond to different functional states. The existence of the D0 is a unique feature of alphacoronavirus. Here we combined cryo-electron tomography (cryo-ET) and cryo-electron microscopy (cryo-EM) to demonstrate in situ the asynchronous S protein D0 motions on intact viral particles of a highly virulent PEDV Pintung 52 strain. We further determined the cryo-EM structure of the recombinant S protein derived from a porcine cell line, which revealed additional domain motions likely associated with receptor binding. By integrating mass spectrometry and cryo-EM, we delineated the complex compositions and spatial distribution of the PEDV S protein N-glycans, and demonstrated the functional role of a key N-glycan in modulating the D0 conformation.

Hsu and co-workers integrate cryo-electron tomography, cryo-electron microscopy and mass spectrometry to reveal the structural polymorphism of a pig coronavirus spike protein within intact viral particles, and how glycosylation modulates the conformational changes pertinent to host recognition.

N-glycosylation profiles of the SARS-CoV-2 spike D614G mutant and its ancestral protein characterized by advanced mass spectrometry

N-glycosylation plays an important role in the structure and function of membrane and secreted proteins. The spike protein on the surface of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, is heavily glycosylated and the major target for developing vaccines, therapeutic drugs and diagnostic tests. The first major SARS-CoV-2 variant carries a D614G substitution in the spike (S-D614G) that has been associated with altered conformation, enhanced ACE2 binding, and increased infectivity and transmission. In this report, we used mass spectrometry techniques to characterize and compare the N-glycosylation of the wild type (S-614D) or variant (S-614G) SARS-CoV-2 spike glycoproteins prepared under identical conditions. The data showed that half of the N-glycosylation sequons changed their distribution of glycans in the S-614G variant. The S-614G variant showed a decrease in the relative abundance of complex-type glycans (up to 45%) and an increase in oligomannose glycans (up to 33%) on all altered sequons. These changes led to a reduction in the overall complexity of the total N-glycosylation profile. All the glycosylation sites with altered patterns were in the spike head while the glycosylation of three sites in the stalk remained unchanged between S-614G and S-614D proteins.