Structural remodeling of SARS-CoV-2 spike protein glycans reveals the regulatory roles in receptor-binding affinity

Structural and Functional Glycosylation of the Abdala COVID-19 Vaccine

Abdala is a COVID-19 vaccine produced in Pichia pastoris and is based on the receptor-binding domain (RBD) of the SARS-CoV-2 spike. Abdala is currently approved for use in multiple countries with clinical trials confirming its safety and efficacy in preventing severe illness and death. Although P. pastoris is used as an expression system for protein-based vaccines, yeast glycosylation remains largely uncharacterised across immunogens. Here, we characterise N-glycan structures and their site of attachment on Abdala and show how yeast-specific glycosylation decreases binding to the ACE2 receptor and a receptor-binding motif (RBM) targeting antibody compared to the equivalent mammalian-derived RBD. Reduced receptor and antibody binding is attributed to changes in conformational dynamics resulting from N-glycosylation. These data highlight the critical importance of glycosylation in vaccine design and demonstrate how individual glycans can influence host interactions and immune recognition via protein structural dynamics.

The accomplices: Heparan sulfates and N-glycans foster SARS-CoV-2 spike:ACE2 receptor binding and virus priming

Although it is well established that the SARS-CoV-2 spike glycoprotein binds to the host cell ACE2 receptor to initiate infection, far less is known about the tissue tropism and host cell susceptibility to the virus. Differential expression across different cell types of heparan sulfate (HS) proteoglycans, with variably sulfated glycosaminoglycans (GAGs), and their synergistic interactions with host and viral N-glycans may contribute to tissue tropism and host cell susceptibility. Nevertheless, their contribution remains unclear since HS and N-glycans evade experimental characterization. We, therefore, carried out microsecond-long all-atom molecular dynamics simulations, followed by random acceleration molecular dynamics simulations, of the fully glycosylated spike:ACE2 complex with and without highly sulfated GAG chains bound. By considering the model GAGs as surrogates for the highly sulfated HS expressed in lung cells, we identified key cell entry mechanisms of spike SARS-CoV-2. We find that HS promotes structural and energetic stabilization of the active conformation of the spike receptor-binding domain (RBD) and reorientation of ACE2 toward the N-terminal domain in the same spike subunit as the RBD. Spike and ACE2 N-glycans exert synergistic effects, promoting better packing, strengthening the protein:protein interaction, and prolonging the residence time of the complex. ACE2 and HS binding trigger rearrangement of the S2' functional protease cleavage site through allosteric interdomain communication. These results thus show that HS has a multifaceted role in facilitating SARS-CoV-2 infection, and they provide a mechanistic basis for the development of GAG derivatives with anti-SARS-CoV-2 potential.

Water-glycan interactions drive the SARS-CoV-2 spike dynamics: insights into glycan-gate control and camouflage mechanisms

To develop therapeutic strategies against COVID-19, we introduce a high-resolution all-atom polarizable model capturing many-body effects of protein, glycan, solvent, and membrane components in SARS-CoV-2 spike protein open and closed states. Employing μs-long molecular dynamics simulations powered by high-performance cloud-computing and unsupervised density-driven adaptive sampling, we investigated the differences in bulk-solvent-glycan and protein-solvent-glycan interfaces between these states. We unraveled a sophisticated solvent-glycan polarization interaction network involving the N165/N343 glycan-gate patterns that provide structural support for the open state and identified key water molecules that could potentially be targeted to destabilize this configuration. In the closed state, the reduced solvent polarization diminishes the overall N165/N343 dipoles, yet internal interactions and a reorganized sugar coat stabilize this state. Despite variations, our glycan-solvent accessibility analysis reveals the glycan shield capability to conserve constant interactions with the solvent, effectively camouflaging the virus from immune detection in both states. The presented insights advance our comprehension of viral pathogenesis at an atomic level, offering potential to combat COVID-19.

Inclusion of deuterated glycopeptides provides increased sequence coverage in hydrogen/deuterium exchange mass spectrometry analysis of SARS-CoV-2 spike glycoprotein

RATIONALE

Hydrogen/deuterium exchange mass spectrometry (HDX-MS) can provide precise analysis of a protein's conformational dynamics across varied states, such as heat-denatured versus native protein structures, localizing regions that are specifically affected by such conditional changes. Maximizing protein sequence coverage provides high confidence that regions of interest were located by HDX-MS, but one challenge for complete sequence coverage is N-glycosylation sites. The deuteration of peptides post-translationally modified by asparagine-bound glycans (glycopeptides) has not always been identified in previous reports of HDX-MS analyses, causing significant sequence coverage gaps in heavily glycosylated proteins and uncertainty in structural dynamics in many regions throughout a glycoprotein.

METHODS

We detected deuterated glycopeptides with a Tribrid Orbitrap Eclipse mass spectrometer performing data-dependent acquisition. An MS scan was used to identify precursor ions; if high-energy collision-induced dissociation MS/MS of the precursor indicated oxonium ions diagnostic for complex glycans, then electron transfer low-energy collision-induced dissociation MS/MS scans of the precursor identified the modified asparagine residue and the glycan's mass. As in traditional HDX-MS, the identified glycopeptides were then analyzed at the MS level in samples labeled with D2 O.

RESULTS

We report HDX-MS analysis of the SARS-CoV-2 spike protein ectodomain in its trimeric prefusion form, which has 22 predicted N-glycosylation sites per monomer, with and without heat treatment. We identified glycopeptides and calculated their average isotopic mass shifts from deuteration. Inclusion of the deuterated glycopeptides increased sequence coverage of spike ectodomain from 76% to 84%, demonstrated that glycopeptides had been deuterated, and improved confidence in results localizing structural rearrangements.

CONCLUSION

Inclusion of deuterated glycopeptides improves the analysis of the conformational dynamics of glycoproteins such as viral surface antigens and cellular receptors.

Expression and purification of the receptor-binding domain of SARS-CoV-2 spike protein in mammalian cells for immunological assays

The receptor-binding domain (RBD) of the spike glycoprotein of SARS-CoV-2 virus mediates the interaction with the host cell and is required for virus internalization. It is, therefore, the primary target of neutralizing antibodies. The receptor-binding domain soon became the major target for COVID-19 research and the development of diagnostic tools and new-generation vaccines. Here, we provide a detailed protocol for high-yield expression and one-step affinity purification of recombinant RBD from transiently transfected Expi293F cells. Expi293F mammalian cells can be grown to extremely high densities in a specially formulated serum-free medium in suspension cultures, which makes them an excellent tool for secreted protein production. The highly purified RBD is glycosylated, structurally intact, and forms homomeric complexes. With this quick and easy method, we are able to produce large quantities of RBD (80 mg·L-1 culture) that we have successfully used in immunological assays to examine antibody titers and seroconversion after mRNA-based vaccination of mice.

When inflammatory stressors dramatically change, disease phenotypes may transform between autoimmune hematopoietic failure and myeloid neoplasms

Aplastic anemia (AA) and hypoplastic myelodysplastic syndrome are paradigms of autoimmune hematopoietic failure (AHF). Myelodysplastic syndrome and acute myeloid leukemia are unequivocal myeloid neoplasms (MNs). Currently, AA is also known to be a clonal hematological disease. Genetic aberrations typically observed in MNs are detected in approximately one-third of AA patients. In AA patients harboring MN-related genetic aberrations, a poor response to immunosuppressive therapy (IST) and an increased risk of transformation to MNs occurring either naturally or after IST are predicted. Approximately 10%-15% of patients with severe AA transform the disease phenotype to MNs following IST, and in some patients, leukemic transformation emerges during or shortly after IST. Phenotypic transformations between AHF and MNs can occur reciprocally. A fraction of advanced MN patients experience an aplastic crisis during which leukemic blasts are repressed. The switch that shapes the disease phenotype is a change in the strength of extramedullary inflammation. Both AHF and MNs have an immune-active bone marrow (BM) environment (BME). In AHF patients, an inflamed BME can be evoked by infiltrated immune cells targeting neoplastic molecules, which contributes to the BM-specific autoimmune impairment. Autoimmune responses in AHF may represent an antileukemic mechanism, and inflammatory stressors strengthen antileukemic immunity, at least in a significant proportion of patients who have MN-related genetic aberrations. During active inflammatory episodes, normal and leukemic hematopoieses are suppressed, which leads to the occurrence of aplastic cytopenia and leukemic cell regression. The successful treatment of underlying infections mitigates inflammatory stress-related antileukemic activities and promotes the penetration of leukemic hematopoiesis. The effect of IST is similar to that of treating underlying infections. Investigating inflammatory stress-powered antileukemic immunity is highly important in theoretical studies and clinical practice, especially given the wide application of immune-activating agents and immune checkpoint inhibitors in the treatment of hematological neoplasms.

Effects of the glycosylation of the receptor binding domain (RBD dimer)-based Covid-19 vaccine (ZF2001) on its humoral immunogenicity and immunoreactivity

The SARS-CoV-2 spike protein receptor-binding domain (RBD), which is a key target for the development of SARS-CoV-2 neutralizing antibodies and vaccines, mediates the binding of the host receptor angiotensin-converting enzyme 2 (ACE2). However, the high heterogeneity of RBD glycoforms may lead to an incomplete neutralization effect and impact the immunogenicity of RBD-based vaccines (Ye et al., 2021). Here, our data suggested that the glycosylation significantly affected the humoral immunogenicity and immunoreactivity of the RBD-dimer-based Covid-19 vaccine (ZF2001) (Yang et al., 2021). Several deglycosylated types of ZF2001 (with sialic acid removed (ZF2001-ΔSA), sialic acid & O-glycans removed (ZF2001-ΔSA&O), N-glycans removed (ZF2001-ΔN), N- & O-glycans removed (ZF2001-ΔN&O)) were obtained by treatment with glycosidases. The binding affinity between deglycosylated types of ZF2001 and ACE2 was slightly weakened and that between deglycosylated types of ZF2001 and several monoclonal antibodies (mAbs) were also changed compared with ZF2001. The results of pseudovirus neutralization assay and binding affinity assay of all ZF2001 types revealed that the antigens with complex glycosylation had better humoral immunogenicity and immunoreactivity. Molecular dynamics simulation indicated that the more complex glycosylation of RBD corresponded to more hydrogen bonds formed between helper T-cell epitopes of RBD and major histocompatibility complex II (MHC-II). In summary, these results demonstrated that the glycosylation of RBD affects antigen presentation, humoral immunogenicity and immunoreactivity, which may be an important consideration for vaccine design and production technology.

Glycosyltransferases as versatile tools to study the biology of glycans

All cells are decorated with complex carbohydrate structures called glycans that serve as ligands for glycan-binding proteins (GBPs) to mediate a wide range of biological processes. Understanding the specific functions of glycans is key to advancing an understanding of human health and disease. However, the lack of convenient and accessible tools to study glycan-based interactions has been a defining challenge in glycobiology. Thus, the development of chemical and biochemical strategies to address these limitations has been a rapidly growing area of research. In this review, we describe the use of glycosyltransferases (GTs) as versatile tools to facilitate a greater understanding of the biological roles of glycans. We highlight key examples of how GTs have streamlined the preparation of well-defined complex glycan structures through chemoenzymatic synthesis, with an emphasis on synthetic strategies allowing for site- and branch-specific display of glyco-epitopes. We also describe how GTs have facilitated expansion of glyco-engineering strategies, on both glycoproteins and cell surfaces. Coupled with advancements in bioorthogonal chemistry, GTs have enabled selective glyco-epitope editing of glycoproteins and cells, selective glycan subclass labeling, and the introduction of novel biomolecule functionalities onto cells, including defined oligosaccharides, antibodies, and other proteins. Collectively, these approaches have contributed great insight into the fundamental biological roles of glycans and are enabling their application in drug development and cellular therapies, leaving the field poised for rapid expansion.

Impact of mutations on the plant-based production of recombinant SARS-CoV-2 RBDs

Subunit vaccines based on recombinant viral antigens are valuable interventions to fight existing and evolving viruses and can be produced at large-scale in plant-based expression systems. The recombinant viral antigens are often derived from glycosylated envelope proteins of the virus and glycosylation plays an important role for the immunogenicity by shielding protein epitopes. The receptor-binding domain (RBD) of the SARS-CoV-2 spike is a principal target for vaccine development and has been produced in plants, but the yields of recombinant RBD variants were low and the role of the N-glycosylation in RBD from different SARS-CoV-2 variants of concern is less studied. Here, we investigated the expression and glycosylation of six different RBD variants transiently expressed in leaves of Nicotiana benthamiana. All of the purified RBD variants were functional in terms of receptor binding and displayed almost full N-glycan occupancy at both glycosylation sites with predominately complex N-glycans. Despite the high structural sequence conservation of the RBD variants, we detected a variation in yield which can be attributed to lower expression and differences in unintentional proteolytic processing of the C-terminal polyhistidine tag used for purification. Glycoengineering towards a human-type complex N-glycan profile with core α1,6-fucose, showed that the reactivity of the neutralizing antibody S309 differs depending on the N-glycan profile and the RBD variant.

Docking and Molecular Dynamics Simulations Clarify Binding Sites for Interactions of Novel Marine Sulfated Glycans with SARS-CoV-2 Spike Glycoprotein

The entry of SARS-CoV-2 into the host cell is mediated by its S-glycoprotein (SGP). Sulfated glycans bind to the SGP receptor-binding domain (RBD), which forms a ternary complex with its receptor angiotensin converting enzyme 2. Here, we have conducted a thorough and systematic computational study of the binding of four oligosaccharide building blocks from novel marine sulfated glycans (isolated from Pentacta pygmaea and Isostichopus badionotus) to the non-glycosylated and glycosylated RBD. Blind docking studies using three docking programs identified five potential cryptic binding sites. Extensive site-targeted docking and molecular dynamics simulations using two force fields confirmed only two binding sites (Sites 1 and 5) for these novel, highly charged sulfated glycans, which were also confirmed by previously published reports. This work showed the structural features and key interactions driving ligand binding. A previous study predicted Site 2 to be a potential binding site, which was not observed here. The use of several molecular modeling approaches gave a comprehensive assessment. The detailed comparative study utilizing multiple modeling approaches is the first of its kind for novel glycan-SGP interaction characterization. This study provided insights into the key structural features of these novel glycans as they are considered for development as potential therapeutics.