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

Protein Complex Heterogeneity and Topology Revealed by Electron Capture Charge Reduction and Surface Induced Dissociation

We illustrate the utility of native mass spectrometry (nMS) combined with a fast, tunable gas-phase charge reduction, electron capture charge reduction (ECCR), for the characterization of protein complex topology and glycoprotein heterogeneity. ECCR efficiently reduces the charge states of tetradecameric GroEL, illustrating Orbitrap m/z measurements to greater than 100,000 m/z. For pentameric C-reactive protein and tetradecameric GroEL, our novel device combining ECCR with surface induced dissociation (SID) reduces the charge states and yields more topologically informative fragmentation. This is the first demonstration that ECCR yields more native-like SID fragmentation. ECCR also significantly improved mass and glycan heterogeneity measurements of heavily glycosylated SARS-CoV-2 spike protein trimer and thyroglobulin dimer. Protein glycosylation is important for structural and functional properties and plays essential roles in many biological processes. The immense heterogeneity in glycosylation sites and glycan structure poses significant analytical challenges that hinder a mechanistic understanding of the biological role of glycosylation. Without ECCR, average mass determination of glycoprotein complexes is available only through charge detection mass spectrometry or mass photometry. With narrow m/z selection windows followed by ECCR, multiple glycoform m/z values are apparent, providing quick global glycoform profiling and providing a future path for glycan localization on individual intact glycoforms.

The Malaria Box molecules: a source for targeting the RBD and NTD cryptic pocket of the spike glycoprotein in SARS-CoV-2

CONTEXT

SARS-CoV-2, responsible for COVID-19, has led to over 500 million infections and more than 6 million deaths globally. There have been limited effective treatments available. The study aims to find a drug that can prevent the virus from entering host cells by targeting specific sites on the virus's spike protein.

METHOD

We examined 13,397 compounds from the Malaria Box library against two specific sites on the spike protein: the receptor-binding domain (RBD) and a predicted cryptic pocket. Using virtual screening, molecular docking, molecular dynamics, and MMPBSA techniques, they evaluated the stability of two compounds. TCMDC-124223 showed high stability and binding energy in the RBD, while TCMDC-133766 had better binding energy in the cryptic pocket. The study also identified that the interacting residues are conserved, which is crucial for addressing various virus variants. The findings provide insights into the potential of small molecules as drugs against the spike protein.

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 IRE1α-XBP1 arm of the unfolded protein response is a host factor activated in SARS-CoV-2 infection

SARS-CoV-2 infection can cause severe pneumonia, wherein exacerbated inflammation plays a major role. This is reminiscent of the process commonly termed cytokine storm, a condition dependent on a disproportionated production of cytokines. This state involves the activation of the innate immune response by viral patterns and coincides with the biosynthesis of the biomass required for viral replication, which may overwhelm the capacity of the endoplasmic reticulum and drive the unfolded protein response (UPR). The UPR is a signal transduction pathway composed of three branches that is initiated by a set of sensors: inositol-requiring protein 1 (IRE1), protein kinase RNA-like ER kinase (PERK), and activating transcription factor 6 (ATF6). These sensors control adaptive processes, including the transcriptional regulation of proinflammatory cytokines. Based on this background, the role of the UPR in SARS-CoV-2 replication and the ensuing inflammatory response was investigated using in vivo and in vitro models of infection. Mice and Syrian hamsters infected with SARS-CoV-2 showed a sole activation of the Ire1α-Xbp1 arm of the UPR associated with a robust production of proinflammatory cytokines. Human lung epithelial cells showed the dependence of viral replication on the expression of UPR-target proteins branching on the IRE1α-XBP1 arm and to a lower extent on the PERK route. Likewise, activation of the IRE1α-XBP1 branch by Spike (S) proteins from different variants of concern was a uniform finding. These results show that the IRE1α-XBP1 system enhances viral replication and cytokine expression and may represent a potential therapeutic target in SARS-CoV-2 severe pneumonia.

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.

Evolving spike-protein N-glycosylation in SARS-CoV-2 variants

Since >3 years, SARS-CoV-2 has plunged humans into a colossal pandemic. Henceforth, multiple waves of infection have swept through the human population, led by variants that were able to partially evade acquired immunity. The co-evolution of SARS-CoV-2 variants with human immunity provides an excellent opportunity to study the interaction between viral pathogens and their human hosts. The heavily N-glycosylated spike-protein of SARS-CoV-2 plays a pivotal role in initiating infection and is the target for host immune-response, both of which are impacted by host-installed N-glycans. Using highly-sensitive DeGlyPHER approach, we compared the N-glycan landscape on spikes of the SARS-CoV-2 Wuhan-Hu-1 strain to seven WHO-defined variants of concern/interest, using recombinantly expressed, soluble spike-protein trimers, sharing same stabilizing-mutations. We found that N-glycan processing is conserved at most sites. However, in multiple variants, processing of N-glycans from high mannose- to complex-type is reduced at sites N165, N343 and N616, implicated in spike-protein function.

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.

Guedes J, Nogueira FC, Silva JL, Castilho LR, de Oliveira GA. Conformational stability of SARS-CoV-2 glycoprotein spike variants

The severe acute respiratory syndrome spread worldwide, causing a pandemic. SARS-CoV-2 mutations have arisen in the spike, a glycoprotein at the viral envelope and an antigenic candidate for vaccines against COVID-19. Here, we present comparative data of the glycosylated full-length ancestral and D614G spike together with three other transmissible strains classified by the World Health Organization as variants of concern: beta, gamma, and delta. By showing that D614G has less hydrophobic surface exposure and trimer persistence, we place D614G with features that support a model of temporary fitness advantage for virus spillover. Furthermore, during the SARS-CoV-2 adaptation, the spike accumulates alterations leading to less structural stability for some variants. The decreased trimer stability of the ancestral and gamma and the presence of D614G uncoupled conformations mean higher ACE-2 affinities compared to the beta and delta strains. Mapping the energetics and flexibility of variants is necessary to improve vaccine development.

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.

Molecular evolution of human coronavirus-NL63, -229E, -HKU1 and -OC43 in hospitalized children in China

Human coronaviruses (HCoVs) HCoV-NL63, HCoV-229E, HCoV-HKU1 and HCoV-OC43 have been circulated in the human population worldwide, and they are associated with a broad range of respiratory diseases with varying severity. However, there are neither effective therapeutic drugs nor licensed vaccines available for the treatment and prevention of infections by the four HCoVs. In this study, we collected nasopharyngeal aspirates of children hospitalized for respiratory tract infection in China during 2014-2018 and conducted next-generation sequencing. Sequences of four HCoVs were then selected for an in-depth analysis. Genome sequences of 2 HCoV-NL63, 8 HCoV-229E, 2 HCoV-HKU1, and 6 HCoV-OC43 were obtained. Based on the full-length S gene, a strong temporal signal was found in HCoV-229E and the molecular evolutionary rate was 6 × 10^-4 substitutions/site/year. Based on the maximum-likelihood (ML) phylogenetic tree of complete S gene, we designated H78 as a new sub-genotype C2 of HCoV-HKU1, and the obtained P43 sequence was grouped into the reported novel genotype K of HCoV-OC43 circulating in Guangzhou, China. Based on the complete genome, potential recombination events were found to occur as two phenomena, namely intraspecies and interspecies. Moreover, we observed two amino acid substitutions in the S1 subunit of obtained HCoV-NL63 (G534V) and HCoV-HKU1 (H512R), while residues 534 and 512 are important for the binding of angiotensin-converting enzyme 2 and neutralizing antibodies, respectively. Our findings might provide a clue for the molecular evolution of the four HCoVs and help in the early diagnosis, treatment and prevention of broad-spectrum HCoV infection.

Investigation of the N-Glycosylation of the SARS-CoV-2 S Protein Contained in VLPs Produced in Nicotiana benthamiana

The emergence of the SARS-CoV-2 coronavirus pandemic in China in late 2019 led to the fast development of efficient therapeutics. Of the major structural proteins encoded by the SARS-CoV-2 genome, the SPIKE (S) protein has attracted considerable research interest because of the central role it plays in virus entry into host cells. Therefore, to date, most immunization strategies aim at inducing neutralizing antibodies against the surface viral S protein. The SARS-CoV-2 S protein is heavily glycosylated with 22 predicted N-glycosylation consensus sites as well as numerous mucin-type O-glycosylation sites. As a consequence, O- and N-glycosylations of this viral protein have received particular attention. Glycans N-linked to the S protein are mainly exposed at the surface and form a shield-masking specific epitope to escape the virus antigenic recognition. In this work, the N-glycosylation status of the S protein within virus-like particles (VLPs) produced in Nicotiana benthamiana (N. benthamiana) was investigated using a glycoproteomic approach. We show that 20 among the 22 predicted N-glycosylation sites are dominated by complex plant N-glycans and one carries oligomannoses. This suggests that the SARS-CoV-2 S protein produced in N. benthamiana adopts an overall 3D structure similar to that of recombinant homologues produced in mammalian cells.

Analysis of memory B cells identifies conserved neutralizing epitopes on the N-terminal domain of variant SARS-Cov-2 spike proteins

SARS-CoV-2 infection or vaccination produces neutralizing antibody responses that contribute to better clinical outcomes. The receptor-binding domain (RBD) and the N-terminal domain (NTD) of the spike trimer (S) constitute the two major neutralizing targets for antibodies. Here, we use NTD-specific probes to capture anti-NTD memory B cells in a longitudinal cohort of infected individuals, some of whom were vaccinated. We found 6 complementation groups of neutralizing antibodies. 58% targeted epitopes outside the NTD supersite, 58% neutralized either Gamma or Omicron, and 14% were broad neutralizers that also neutralized Omicron. Structural characterization revealed that broadly active antibodies targeted three epitopes outside the NTD supersite including a class that recognized both the NTD and SD2 domain. Rapid recruitment of memory B cells producing these antibodies into the plasma cell compartment upon re-infection likely contributes to the relatively benign course of subsequent infections with SARS-CoV-2 variants, including Omicron.

Conserved Neutralizing Epitopes on the N-Terminal Domain of Variant SARS-CoV-2 Spike Proteins

SARS-CoV-2 infection or vaccination produces neutralizing antibody responses that contribute to better clinical outcomes. The receptor binding domain (RBD) and the N-terminal domain (NTD) of the spike trimer (S) constitute the two major neutralizing targets for the antibody system. Neutralizing antibodies targeting the RBD bind to several different sites on this domain. In contrast, most neutralizing antibodies to NTD characterized to date bind to a single supersite, however these antibodies were obtained by methods that were not NTD specific. Here we use NTD specific probes to focus on anti-NTD memory B cells in a cohort of pre-omicron infected individuals some of which were also vaccinated. Of 275 NTD binding antibodies tested 103 neutralized at least one of three tested strains: Wuhan-Hu-1, Gamma, or PMS20, a synthetic variant which is extensively mutated in the NTD supersite. Among the 43 neutralizing antibodies that were further characterized, we found 6 complementation groups based on competition binding experiments. 58% targeted epitopes outside the NTD supersite, and 58% neutralized either Gamma or Omicron, but only 14% were broad neutralizers. Three of the broad neutralizers were characterized structurally. C1520 and C1791 recognize epitopes on opposite faces of the NTD with a distinct binding pose relative to previously described antibodies allowing for greater potency and cross-reactivity with 7 different variants including Beta, Delta, Gamma and Omicron. Antibody C1717 represents a previously uncharacterized class of NTD-directed antibodies that recognizes the viral membrane proximal side of the NTD and SD2 domain, leading to cross-neutralization of Beta, Gamma and Omicron. We conclude SARS-CoV-2 infection and/or Wuhan-Hu-1 mRNA vaccination produces a diverse collection of memory B cells that produce anti-NTD antibodies some of which can neutralize variants of concern. Rapid recruitment of these cells into the antibody secreting plasma cell compartment upon re-infection likely contributes to the relatively benign course of subsequent infections with SARS-CoV-2 variants including omicron.