Fine-tuning the spike: role of the nature and topology of the glycan shield in the structure and dynamics of the SARS-CoV-2 S

The N545S and K717N substitution at the N-glycosylation sites of the S2 subunit of avian infectious bronchitis virus can significantly enhance viral pathogenicity

The S2 subunit of infectious bronchitis virus (IBV) is a heavily glycosylated protein that can impact various characteristics of the virus. It is currently known that N-glycosylation modifications are predominantly located on the S2 subunit. However, the exact role of their N-glycosylation modification remains undisclosed. To elucidate the function of these N-glycosylation sites, we identified 14 common sites distributed on the S2 subunit of the 5 genotypes of IBV in present study. Subsequently, we selected 7 sites to generate mutants and assessed their impact on viral virulence, replication ability, and antigenicity. Our finding revealed that only 2 substitutions, N545S and K717N, increased the viral replication titer and antigenicity, and ultimately the pathogenicity in chicks. To delve into the mechanisms underlying this increased pathogenicity, we discovered that K717N can change the structure of antigenic epitopes. The N545S substitution not only influenced antigenic epitope structure, but also enhanced the ability of the virus to enter CEKs during the early stages of viral replication. These results suggest that the enhanced viral pathogenicity associated with N545S and K717N substitutions is multifaceted, with acceleration of the viral membrane fusion process and alterations in epitope structure representing crucial factors in the capability of N-glycosylation modifications to boost viral virulence. These insights provide valuable guidance for the efficient development of live attenuated vaccines.

Role of N343 glycosylation on the SARS-CoV-2 S RBD structure and co-receptor binding across variants of concern

Glycosylation of the SARS-CoV-2 spike (S) protein represents a key target for viral evolution because it affects both viral evasion and fitness. Successful variations in the glycan shield are difficult to achieve though, as protein glycosylation is also critical to folding and structural stability. Within this framework, the identification of glycosylation sites that are structurally dispensable can provide insight into the evolutionary mechanisms of the shield and inform immune surveillance. In this work, we show through over 45 μs of cumulative sampling from conventional and enhanced molecular dynamics (MD) simulations, how the structure of the immunodominant S receptor binding domain (RBD) is regulated by N-glycosylation at N343 and how this glycan's structural role changes from WHu-1, alpha (B.1.1.7), and beta (B.1.351), to the delta (B.1.617.2), and omicron (BA.1 and BA.2.86) variants. More specifically, we find that the amphipathic nature of the N-glycan is instrumental to preserve the structural integrity of the RBD hydrophobic core and that loss of glycosylation at N343 triggers a specific and consistent conformational change. We show how this change allosterically regulates the conformation of the receptor binding motif (RBM) in the WHu-1, alpha, and beta RBDs, but not in the delta and omicron variants, due to mutations that reinforce the RBD architecture. In support of these findings, we show that the binding of the RBD to monosialylated ganglioside co-receptors is highly dependent on N343 glycosylation in the WHu-1, but not in the delta RBD, and that affinity changes significantly across VoCs. Ultimately, the molecular and functional insight we provide in this work reinforces our understanding of the role of glycosylation in protein structure and function and it also allows us to identify the structural constraints within which the glycosylation site at N343 can become a hotspot for mutations in the SARS-CoV-2 S glycan shield.

Decoding glycosylation potential from protein structure across human glycoproteins with a multi-view recurrent neural network

Glycosylation is described as a non-templated biosynthesis. Yet, the template-free premise is antithetical to the observation that different N-glycans are consistently placed at specific sites. It has been proposed that glycosite-proximal protein structures could constrain glycosylation and explain the observed microheterogeneity. Using site-specific glycosylation data, we trained a hybrid neural network to parse glycosites (recurrent neural network) and match them to feasible N-glycosylation events (graph neural network). From glycosite-flanking sequences, the algorithm predicts most human N-glycosylation events documented in the GlyConnect database and proposed structures corresponding to observed monosaccharide composition of the glycans at these sites. The algorithm also recapitulated glycosylation in Enhanced Aromatic Sequons, SARS-CoV-2 spike, and IgG3 variants, thus demonstrating the ability of the algorithm to predict both glycan structure and abundance. Thus, protein structure constrains glycosylation, and the neural network enables predictive in silico glycosylation of uncharacterized or novel protein sequences and genetic variants.

Glycan-directed SARS-CoV-2 inhibition by leek extract and lectins with insights into the mode-of-action of Concanavalin A

Four years after its outbreak, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains a global challenge for human health. At its surface, SARS-CoV-2 features numerous extensively glycosylated spike proteins. This glycan coat supports virion docking and entry into host cells and at the same time renders the virus less susceptible to neutralizing antibodies. Given the high genetic plasticity of SARS-CoV-2 and the rapid emergence of immune escape variants, targeting the glycan shield by carbohydrate-binding agents emerges as a promising strategy. However, the potential of carbohydrate-targeting reagents as viral inhibitors remains underexplored. Here, we tested seven plant-derived carbohydrate-binding proteins, called lectins, and one crude plant extract for their antiviral activity against SARS-CoV-2 in two types of human lung cells: A549 cells ectopically expressing the ACE2 receptor and Calu-3 cells. We identified three lectins and an Allium porrum (leek) extract inhibiting SARS-CoV-2 infection in both cell systems with selectivity indices (SI) ranging between >2 and >299. Amongst these, the lectin Concanavalin A (Con A) exerted the most potent and broad activity against a panel of SARS-CoV-2 variants. We used multiplex super-resolution microscopy to address lectin interactions with SARS-CoV-2 and its host cells. Notably, we discovered that Con A not only binds to SARS-CoV-2 virions and their host cells, but also causes SARS-CoV-2 aggregation. Thus, Con A exerts a dual mode-of-action comprising both, antiviral and virucidal, mechanisms. These results establish Con A and other plant lectins as candidates for COVID-19 prevention and basis for further drug development.

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.

Molecular insights into the adaptive evolution of SARS-CoV-2 spike protein

The COVID-19 pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has substantially damaged the global economy and human health. The spike (S) protein of coronaviruses plays a pivotal role in viral entry by binding to host cell receptors. Additionally, it acts as the primary target for neutralizing antibodies in those infected and is the central focus for currently utilized or researched vaccines. During the virus's adaptation to the human host, the S protein of SARS-CoV-2 has undergone significant evolution. As the COVID-19 pandemic has unfolded, new mutations have arisen and vanished, giving rise to distinctive amino acid profiles within variant of concern strains of SARS-CoV-2. Notably, many of these changes in the S protein have been positively selected, leading to substantial alterations in viral characteristics, such as heightened transmissibility and immune evasion capabilities. This review aims to provide an overview of our current understanding of the structural implications associated with key amino acid changes in the S protein of SARS-CoV-2. These research findings shed light on the intricate and dynamic nature of viral evolution, underscoring the importance of continuous monitoring and analysis of viral genomes. Through these molecular-level investigations, we can attain deeper insights into the virus's adaptive evolution, offering valuable guidance for designing vaccines and developing antiviral drugs to combat the ever-evolving viral threats.

N-Glycosylation as a Modulator of Protein Conformation and Assembly in Disease

Glycosylation, a prevalent post-translational modification, plays a pivotal role in regulating intricate cellular processes by covalently attaching glycans to macromolecules. Dysregulated glycosylation is linked to a spectrum of diseases, encompassing cancer, neurodegenerative disorders, congenital disorders, infections, and inflammation. This review delves into the intricate interplay between glycosylation and protein conformation, with a specific focus on the profound impact of N-glycans on the selection of distinct protein conformations characterized by distinct interactomes-namely, protein assemblies-under normal and pathological conditions across various diseases. We begin by examining the spike protein of the SARS virus, illustrating how N-glycans regulate the infectivity of pathogenic agents. Subsequently, we utilize the prion protein and the chaperone glucose-regulated protein 94 as examples, exploring instances where N-glycosylation transforms physiological protein structures into disease-associated forms. Unraveling these connections provides valuable insights into potential therapeutic avenues and a deeper comprehension of the molecular intricacies that underlie disease conditions. This exploration of glycosylation's influence on protein conformation effectively bridges the gap between the glycome and disease, offering a comprehensive perspective on the therapeutic implications of targeting conformational mutants and their pathologic assemblies in various diseases. The goal is to unravel the nuances of these post-translational modifications, shedding light on how they contribute to the intricate interplay between protein conformation, assembly, and disease.

Many Roles of Carbohydrates: A Computational Spotlight on the Coronavirus S Protein Binding

Glycosylation is one of the post-translational modifications with more than 50% of human proteins being glycosylated. The exact nature and chemical composition of glycans are inaccessible to X-ray or cryo-electron microscopy imaging techniques. Therefore, computational modeling studies and molecular dynamics must be used as a "computational microscope". The spike (S) protein of SARS-CoV-2 is heavily glycosylated, and a few glycans play a more functional role "beyond shielding". In this mini-review, we discuss computational investigations of the roles of specific S-protein and ACE2 glycans in the overall ACE2-S protein binding. We highlight different functions of specific glycans demonstrated in myriad computational models and simulations in the context of the SARS-CoV-2 virus binding to the receptor. We also discuss interactions between glycocalyx and the S protein, which may be utilized to design prophylactic polysaccharide-based therapeutics targeting the S protein. In addition, we underline the recent emergence of coronavirus variants and their impact on the S protein and its glycans.

Stabilization of the Metastable Pre-Fusion Conformation of the SARS-CoV-2 Spike Glycoprotein through N-Linked Glycosylation of the S2 Subunit

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the novel coronavirus responsible for the coronavirus disease 2019 (COVID-19) pandemic, represents a serious threat to public health. The spike (S) glycoprotein of SARS-CoV-2 mediates viral entry into host cells and is heavily glycosylated. In this study, we systemically analyzed the roles of 22 putative N-linked glycans in SARS-CoV-2 S protein expression, membrane fusion, viral entry, and stability. Using the α-glycosidase inhibitors castanospermine and NB-DNJ, we confirmed that disruption of N-linked glycosylation blocked the maturation of the S protein, leading to the impairment of S protein-mediated membrane fusion. Single-amino-acid substitution of each of the 22 N-linked glycosylation sites with glutamine revealed that 9 out of the 22 N-linked glycosylation sites were critical for S protein folding and maturation. Thus, substitution at these sites resulted in reduced S protein-mediated cell-cell fusion and viral entry. Notably, the N1074Q mutation markedly affected S protein stability and induced significant receptor-independent syncytium (RIS) formation in HEK293T/hACE2-KO cells. Additionally, the removal of the furin cleavage site partially compensated for the instability induced by the N1074Q mutation. Although the corresponding mutation in the SARS-CoV S protein (N1056Q) did not induce RIS in HEK293T cells, the N669Q and N1080Q mutants exhibited increased fusogenic activity and did induce syncytium formation in HEK293T cells. Therefore, N-glycans on the SARS-CoV and SARS-CoV-2 S2 subunits are highly important for maintaining the pre-fusion state of the S protein. This study revealed the critical roles of N-glycans in S protein maturation and stability, information that has implications for the design of vaccines and antiviral strategies.

Broad receptor tropism and immunogenicity of a clade 3 sarbecovirus

Although Rhinolophus bats harbor diverse clade 3 sarbecoviruses, the structural determinants of receptor tropism along with the antigenicity of their spike (S) glycoproteins remain uncharacterized. Here, we show that the African Rhinolophus bat clade 3 sarbecovirus PRD-0038 S has a broad angiotensin-converting enzyme 2 (ACE2) usage and that receptor-binding domain (RBD) mutations further expand receptor promiscuity and enable human ACE2 utilization. We determine a cryo-EM structure of the PRD-0038 RBD bound to Rhinolophus alcyone ACE2, explaining receptor tropism and highlighting differences with SARS-CoV-1 and SARS-CoV-2. Characterization of PRD-0038 S using cryo-EM and monoclonal antibody reactivity reveals its distinct antigenicity relative to SARS-CoV-2 and identifies PRD-0038 cross-neutralizing antibodies for pandemic preparedness. PRD-0038 S vaccination elicits greater titers of antibodies cross-reacting with vaccine-mismatched clade 2 and clade 1a sarbecoviruses compared with SARS-CoV-2 S due to broader antigenic targeting, motivating the inclusion of clade 3 antigens in next-generation vaccines for enhanced resilience to viral evolution.

Investigating the role of glycans in Omicron sub-lineages XBB.1.5 and XBB.1.16 binding to host receptor using molecular dynamics and binding free energy calculations

Omicron derived lineages viz. BA.2, BA.3, BA.4 BA.5, BF.7 and XBBs show prominence with improved immune escape, transmissibility, infectivity, and pathogenicity in general. Sub-variants, XBB.1.5 and XBB.1.16 have shown rapid spread, with mutations embedded throughout the viral genome, including the spike protein. Changing atomic landscapes in spike contributes significantly to modulate host pathogen interactions and infections thereof. In the present work, we computationally analyzed the binding affinities of spike receptor binding domains (RBDs) of XBB.1.5 and XBB.1.16 towards human angiotensin-converting enzyme 2 (hACE2) compared to Omicron. We have employed simulations and binding energy estimation of molecular complexes of spike-hACE2 to assess the interplay of interaction pattern and effect of mutations if any in the binding mode of the RBDs of these novel mutants. We calculated the binding free energy (BFE) of the RBD of the Omicron, XBB.1.5 and XBB.1.16 spike protein to hACE2. We showed that XBB.1.5 and XBB.1.16 can bind to human cells more strongly than Omicron due to the increased charge of the RBD, which enhances the electrostatic interactions with negatively charged hACE2. The per-residue decompositions further show that the Asp339His, Asp405Asn and Asn460Lys mutations in the XBBs RBD play a crucial role in enhancing the electrostatic interactions, by acquiring positively charged residues, thereby influencing the formation/loss of interfacial bonds and thus strongly affecting the spike RBD-hACE2 binding affinity. Simulation results also indicate less interference of heterogeneous glycans of XBB.1.5 spike RBD towards binding to hACE2. Moreover, despite having less interaction at the three interfacial contacts between XBB S RBD and hACE2 compared to Omicron, variants XBB.1.5 and XBB.1.16 had higher total binding free energies (ΔGbind) than Omicron due to the contribution of non-interfacial residues to the free energy, providing insight into the increased binding affinity of XBB1.5 and XBB.1.16. Furthermore, the presence of large positively charged surface patches in the XBBs act as drivers of electrostatic interactions, thus support the possibility of a higher binding affinity to hACE2.

Comparative Analysis of Conformational Dynamics and Systematic Characterization of Cryptic Pockets in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 Spike Complexes with the ACE2 Host Receptor: Confluence of Binding and Structural Plasticity in Mediating Networks of Conserved Allosteric Sites

In the current study, we explore coarse-grained simulations and atomistic molecular dynamics together with binding energetics scanning and cryptic pocket detection in a comparative examination of conformational landscapes and systematic characterization of allosteric binding sites in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 spike full-length trimer complexes with the host receptor ACE2. Microsecond simulations, Markov state models and mutational scanning of binding energies of the SARS-CoV-2 BA.2 and BA.2.75 receptor binding domain complexes revealed the increased thermodynamic stabilization of the BA.2.75 variant and significant dynamic differences between these Omicron variants. Molecular simulations of the SARS-CoV-2 Omicron spike full-length trimer complexes with the ACE2 receptor complemented atomistic studies and enabled an in-depth analysis of mutational and binding effects on conformational dynamic and functional adaptability of the Omicron variants. Despite considerable structural similarities, Omicron variants BA.2, BA.2.75 and XBB.1 can induce unique conformational dynamic signatures and specific distributions of the conformational states. Using conformational ensembles of the SARS-CoV-2 Omicron spike trimer complexes with ACE2, we conducted a comprehensive cryptic pocket screening to examine the role of Omicron mutations and ACE2 binding on the distribution and functional mechanisms of the emerging allosteric binding sites. This analysis captured all experimentally known allosteric sites and discovered networks of inter-connected and functionally relevant allosteric sites that are governed by variant-sensitive conformational adaptability of the SARS-CoV-2 spike structures. The results detailed how ACE2 binding and Omicron mutations in the BA.2, BA.2.75 and XBB.1 spike complexes modulate the distribution of conserved and druggable allosteric pockets harboring functionally important regions. The results are significant for understanding the functional roles of druggable cryptic pockets that can be used for allostery-mediated therapeutic intervention targeting conformational states of the Omicron variants.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

Broad receptor tropism and immunogenicity of a clade 3 sarbecovirus

Although Rhinolophus bats harbor diverse clade 3 sarbecoviruses, the structural determinants of receptor tropism along with the antigenicity of their spike (S) glycoproteins remain uncharacterized. Here, we show that the African Rinolophus bat clade 3 sarbecovirus PRD-0038 S has a broad ACE2 usage and that RBD mutations further expand receptor promiscuity and enable human ACE2 utilization. We determined a cryoEM structure of the PRD-0038 RBD bound to R. alcyone ACE2, explaining receptor tropism and highlighting differences with SARS-CoV-1 and SARS-CoV-2. Characterization of PRD-0038 S using cryoEM and monoclonal antibody reactivity revealed its distinct antigenicity relative to SARS-CoV-2 and identified PRD-0038 cross-neutralizing antibodies for pandemic preparedness. PRD-0038 S vaccination elicited greater titers of antibodies cross-reacting with vaccine-mismatched clade 2 and clade 1a sarbecoviruses compared to SARS-CoV-2 S due to broader antigenic targeting, motivating the inclusion of clade 3 antigens in next-generation vaccines for enhanced resilience to viral evolution.

Development of an integrated and inferenceable RDF database of glycan, pathogen and disease resources

Glycans are known to play extremely important roles in infections by viruses and pathogens. In fact, the SARS-CoV-2 virus has been shown to have evolved due to a single change in glycosylation. However, data resources on glycans, pathogens and diseases are not well organized. To accurately obtain such information from these various resources, we have constructed a foundation for discovering glycan and virus interaction data using Semantic Web technologies to be able to semantically integrate such heterogeneous data. Here, we created an ontology to encapsulate the semantics of virus-glycan interactions, and used Resource Description Framework (RDF) to represent the data we obtained from non-RDF related databases and data associated with literature. These databases include PubChem, SugarBind, and PSICQUIC, which made it possible to refer to other RDF resources such as UniProt and GlyTouCan. We made these data publicly available as open data and provided a service that allows anyone to freely perform searches using SPARQL. In addition, the RDF resources created in this study are available at the GlyCosmos Portal.

How aberrant N-glycosylation can alter protein functionality and ligand binding: An atomistic view

Protein-assembly defects due to an enrichment of aberrant conformational protein variants are emerging as a new frontier in therapeutics design. Understanding the structural elements that rewire the conformational dynamics of proteins and pathologically perturb functionally oriented ensembles is important for inhibitor development. Chaperones are hub proteins for the assembly of multiprotein complexes and an enrichment of aberrant conformers can affect the cellular proteome, and in turn, phenotypes. Here, we integrate computational and experimental tools to investigte how N-glycosylation of specific residues in glucose-regulated protein 94 (GRP94) modulates internal dynamics and alters the conformational fitness of regions fundamental for the interaction with ATP and synthetic ligands and impacts substructures important for the recognition of interacting proteins. N-glycosylation plays an active role in modulating the energy landscape of GRP94, and we provide support for leveraging the knowledge on distinct glycosylation variants to design molecules targeting GRP94 disease-associated conformational states and assemblies.

An innovative strategy to investigate microbial protein modifications in a reliable fast and sensitive way: A therapy oriented proof of concept based on UV-C irradiation of SARS-CoV-2 spike protein

The characterization of modifications of microbial proteins is of primary importance to dissect pathogen lifecycle mechanisms and could be useful in identifying therapeutic targets. Attempts to solve this issue yielded only partial and non-exhaustive results. We developed a multidisciplinary approach by coupling in vitro infection assay, mass spectrometry (MS), protein 3D modelling, and surface plasma resonance (SPR). As a proof of concept, the effect of low UV-C (273 nm) irradiation on SARS-CoV-2 spike (S) protein was investigated. Following UV-C exposure, MS analysis identified, among other modifications, the disruption of a disulphide bond within the conserved S2 subunit of S protein. Computational analyses revealed that this bond breakage associates with an allosteric effect resulting in the generation of a closed conformation with a reduced ability to bind the ACE2 receptor. The UV-C-induced reduced affinity of S protein for ACE2 was further confirmed by SPR analyses and in vitro infection assays. This comprehensive approach pinpoints the S2 domain of S protein as a potential therapeutic target to prevent SARS-CoV-2 infection. Notably, this workflow could be used to screen a wide variety of microbial protein domains, resulting in a precise molecular fingerprint and providing new insights to adequately address future epidemics.

Mechanism and evolution of human ACE2 binding by SARS-CoV-2 spike

Spike glycoprotein of SARS-CoV-2 mediates viral entry into host cells by facilitating virus attachment and membrane fusion. ACE2 is the main receptor of SARS-CoV-2 and its interaction with spike has shaped the virus' emergence from an animal reservoir and subsequent evolution in the human host. Many structural studies on the spike:ACE2 interaction have provided insights into mechanisms driving viral evolution during the on-going pandemic. This review describes the molecular basis of spike binding to ACE2, outlines mechanisms that have optimised this interaction during viral evolution, and suggests directions for future research.

Heme binding to the SARS-CoV-2 spike glycoprotein

The target for humoral immunity, SARS-CoV-2 spike glycoprotein, has become the focus of vaccine research and development. Previous work demonstrated that the N-terminal domain (NTD) of SARS-CoV-2 spike binds biliverdin-a product of heme catabolism-causing a strong allosteric effect on the activity of a subset of neutralizing antibodies. Herein, we show that the spike glycoprotein is also able to bind heme (KD = 0.5 ± 0.2 μM). Molecular modeling indicated that the heme group fits well within the same pocket on the SARS-CoV-2 spike NTD. Lined by aromatic and hydrophobic residues (W104, V126, I129, F192, F194, I203, and L226), the pocket provides a suitable environment to stabilize the hydrophobic heme. Mutagenesis of N121 has a substantive effect on heme binding (KD = 3000 ± 220 μM), confirming the pocket as a major heme binding location of the viral glycoprotein. Coupled oxidation experiments in the presence of ascorbate indicated that the SARS-CoV-2 glycoprotein can catalyze the slow conversion of heme to biliverdin. The heme trapping and oxidation activities of the spike may allow the virus to reduce levels of free heme during infection to facilitate evasion of the adaptive and innate immunity.

Evolution of Sequence and Structure of SARS-CoV-2 Spike Protein: A Dynamic Perspective

Novel coronavirus (SARS-CoV-2) enters its host cell through a surface spike protein. The viral spike protein has undergone several modifications/mutations at the genomic level, through which it modulated its structure-function and passed through several variants of concern. Recent advances in high-resolution structure determination and multiscale imaging techniques, cost-effective next-generation sequencing, and development of new computational methods (including information theory, statistical methods, machine learning, and many other artificial intelligence-based techniques) have hugely contributed to the characterization of sequence, structure, function of spike proteins, and its different variants to understand viral pathogenesis, evolutions, and transmission. Laying on the foundation of the sequence-structure-function paradigm, this review summarizes not only the important findings on structure/function but also the structural dynamics of different spike components, highlighting the effects of mutations on them. As dynamic fluctuations of three-dimensional spike structure often provide important clues for functional modulation, quantifying time-dependent fluctuations of mutational events over spike structure and its genetic/amino acidic sequence helps identify alarming functional transitions having implications for enhanced fusogenicity and pathogenicity of the virus. Although these dynamic events are more difficult to capture than quantifying a static, average property, this review encompasses those challenging aspects of characterizing the evolutionary dynamics of spike sequence and structure and their implications for functions.

Respiratory viruses interacting with cells: the importance of electrostatics

The COVID-19 pandemic has rekindled interest in the molecular mechanisms involved in the early steps of infection of cells by viruses. Compared to SARS-CoV-1 which only caused a relatively small albeit deadly outbreak, SARS-CoV-2 has led to fulminant spread and a full-scale pandemic characterized by efficient virus transmission worldwide within a very short time. Moreover, the mutations the virus acquired over the many months of virus transmission, particularly those seen in the Omicron variant, have turned out to result in an even more transmissible virus. Here, we focus on the early events of virus infection of cells. We review evidence that the first decisive step in this process is the electrostatic interaction of the spike protein with heparan sulfate chains present on the surface of target cells: Patches of cationic amino acids located on the surface of the spike protein can interact intimately with the negatively charged heparan sulfate chains, which results in the binding of the virion to the cell surface. In a second step, the specific interaction of the receptor binding domain (RBD) within the spike with the angiotensin-converting enzyme 2 (ACE2) receptor leads to the uptake of bound virions into the cell. We show that these events can be expressed as a semi-quantitative model by calculating the surface potential of different spike proteins using the Adaptive Poison-Boltzmann-Solver (APBS). This software allows visualization of the positive surface potential caused by the cationic patches, which increased markedly from the original Wuhan strain of SARS-CoV-2 to the Omicron variant. The surface potential thus enhanced leads to a much stronger binding of the Omicron variant as compared to the original wild-type virus. At the same time, data taken from the literature demonstrate that the interaction of the RBD of the spike protein with the ACE2 receptor remains constant within the limits of error. Finally, we briefly digress to other viruses and show the usefulness of these electrostatic processes and calculations for cell-virus interactions more generally.

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.

How helpful were molecular dynamics simulations in shaping our understanding of SARS-CoV-2 spike protein dynamics?

The SARS-CoV-2 spike protein (S) represents an important viral component that is required for successful viral infection in humans owing to its essential role in recognition of and entry to host cells. The spike is also an appealing target for drug designers who develop vaccines and antivirals. This article is important as it summarizes how molecular simulations successfully shaped our understanding of spike conformational behavior and its role in viral infection. MD simulations found that the higher affinity of SARS-CoV-2-S to ACE2 is linked to its unique residues that add extra electrostatic and van der Waal interactions in comparison to the SARS-CoV S. This illustrates the spread potential of the pandemic SARS-CoV-2 relative to the epidemic SARS-CoV. Different mutations at the S-ACE2 interface, which is believed to increase the transmission of the new variants, affected the behavior and binding interactions in different simulations. The contributions of glycans to the opening of S were revealed via simulations. The immune evasion of S was linked to the spatial distribution of glycans. This help the virus to escape the immune system recognition. This article is important as it summarizes how molecular simulations successfully shaped our understanding of spike conformational behavior and its role in viral infection. This will pave the way to us preparing for the next pandemic as the computational tools are tailored to help fight new challenges.

Lectin Fingerprinting Distinguishes Antibody Neutralization in SARS-CoV-2

Enveloped viruses co-opt host glycosylation pathways to decorate their surface proteins. As viruses evolve, emerging strains can modify their glycosylation patterns to influence host interactions and subvert immune recognition. Still, changes in viral glycosylation or their impact on antibody protection cannot be predicted from genomic sequences alone. Using the highly glycosylated SARS-CoV-2 Spike protein as a model system, we present a lectin fingerprinting method that rapidly reports on changes in variant glycosylation state, which are linked to antibody neutralization. In the presence of antibodies or convalescent and vaccinated patient sera, unique lectin fingerprints emerge that distinguish neutralizing versus non-neutralizing antibodies. This information could not be inferred from direct binding interactions between antibodies and the Spike receptor-binding domain (RBD) binding data alone. Comparative glycoproteomics of the Spike RBD of wild-type (Wuhan-Hu-1) and Delta (B.1.617.2) variants reveal O-glycosylation differences as a key determinant of immune recognition differences. These data underscore the interplay between viral glycosylation and immune recognition and reveal lectin fingerprinting to be a rapid, sensitive, and high-throughput assay to distinguish the neutralization potential of antibodies that target critical viral glycoproteins.

The diversity of the glycan shield of sarbecoviruses related to SARS-CoV-2

Animal reservoirs of sarbecoviruses represent a significant risk of emergent pandemics, as evidenced by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic. Vaccines remain successful at limiting severe disease and death, but the potential for further coronavirus zoonosis motivates the search for pan-coronavirus vaccines. This necessitates a better understanding of the glycan shields of coronaviruses, which can occlude potential antibody epitopes on spike glycoproteins. Here, we compare the structure of 12 sarbecovirus glycan shields. Of the 22 N-linked glycan attachment sites present on SARS-CoV-2, 15 are shared by all 12 sarbecoviruses. However, there are significant differences in the processing state at glycan sites in the N-terminal domain, such as N165. Conversely, glycosylation sites in the S2 domain are highly conserved and contain a low abundance of oligomannose-type glycans, suggesting a low glycan shield density. The S2 domain may therefore provide a more attractive target for immunogen design efforts aiming to generate a pan-coronavirus antibody response.

Glycosaminoglycans: What Remains To Be Deciphered?

Glycosaminoglycans (GAGs) are complex polysaccharides exhibiting a vast structural diversity and fulfilling various functions mediated by thousands of interactions in the extracellular matrix, at the cell surface, and within the cells where they have been detected in the nucleus. It is known that the chemical groups attached to GAGs and GAG conformations comprise "glycocodes" that are not yet fully deciphered. The molecular context also matters for GAG structures and functions, and the influence of the structure and functions of the proteoglycan core proteins on sulfated GAGs and vice versa warrants further investigation. The lack of dedicated bioinformatic tools for mining GAG data sets contributes to a partial characterization of the structural and functional landscape and interactions of GAGs. These pending issues will benefit from the development of new approaches reviewed here, namely (i) the synthesis of GAG oligosaccharides to build large and diverse GAG libraries, (ii) GAG analysis and sequencing by mass spectrometry (e.g., ion mobility-mass spectrometry), gas-phase infrared spectroscopy, recognition tunnelling nanopores, and molecular modeling to identify bioactive GAG sequences, biophysical methods to investigate binding interfaces, and to expand our knowledge and understanding of glycocodes governing GAG molecular recognition, and (iii) artificial intelligence for in-depth investigation of GAGomic data sets and their integration with proteomics.

O-Linked Sialoglycans Modulate the Proteolysis of SARS-CoV-2 Spike and Likely Contribute to the Mutational Trajectory in Variants of Concern

The emergence of a polybasic cleavage motif for the protease furin in SARS-CoV-2 spike has been established as a major factor for human viral transmission. The region N-terminal to that motif is extensively mutated in variants of concern (VOCs). Besides furin, spikes from these variants appear to rely on other proteases for maturation, including TMPRSS2. Glycans near the cleavage site have raised questions about proteolytic processing and the consequences of variant-borne mutations. Here, we identify that sialic acid-containing O-linked glycans on Thr678 of SARS-CoV-2 spike influence furin and TMPRSS2 cleavage and posit O-linked glycosylation as a likely driving force for the emergence of VOC mutations. We provide direct evidence that the glycosyltransferase GalNAc-T1 primes glycosylation at Thr678 in the living cell, an event that is suppressed by mutations in the VOCs Alpha, Delta, and Omicron. We found that the sole incorporation of N-acetylgalactosamine did not impact furin activity in synthetic O-glycopeptides, but the presence of sialic acid reduced the furin rate by up to 65%. Similarly, O-glycosylation with a sialylated trisaccharide had a negative impact on TMPRSS2 cleavage. With a chemistry-centered approach, we substantiate O-glycosylation as a major determinant of spike maturation and propose disruption of O-glycosylation as a substantial driving force for VOC evolution.

Variations within the Glycan Shield of SARS-CoV-2 Impact Viral Spike Dynamics

The emergence of SARS-CoV-2 variants alters the efficacy of existing immunity, whether arisen naturally or through vaccination. Understanding the structure of the viral spike assists in determining the impact of mutations on the antigenic surface. One class of mutation impacts glycosylation attachment sites, which have the capacity to influence the antigenic structure beyond the immediate site of attachment. Here, we compare the site-specific glycosylation of recombinant viral spike mimetics of B.1.351 (Beta), P.1 (Gamma), B.1.617.2 (Delta), B.1.1.529 (Omicron). The P.1 strain exhibits two additional N-linked glycan sites compared to the other variants analyzed and we investigate the impact of these glycans by molecular dynamics. The acquired N188 site is shown to exhibit very limited glycan maturation, consistent with limited enzyme accessibility. Structural modeling and molecular dynamics reveal that N188 is located within a cavity by the receptor binding domain, which influences the dynamics of these attachment domains. These observations suggest a mechanism whereby mutations affecting viral glycosylation sites have a structural impact across the protein surface.

Pineapple Lectin AcmJRL Binds SARS-CoV-2 Spike Protein in a Carbohydrate-Dependent Fashion

The highly glycosylated spike protein of SARS-CoV-2 is essential for infection and constitutes a prime target for antiviral agents and vaccines. The pineapple-derived jacalin-related lectin AcmJRL is present in the medication bromelain in significant quantities and has previously been described to bind mannosides. Here, we performed a large ligand screening of AcmJRL by glycan array analysis, quantified the interaction with carbohydrates and validated high-mannose glycans as preferred ligands. Because the SARS-CoV-2 spike protein was previously reported to carry a high proportion of high-mannose N-glycans, we tested the binding of AcmJRL to the recombinantly produced extraviral domain of spike protein. We could demonstrate that AcmJRL binds the spike protein with a low-micromolar KD in a carbohydrate-dependent fashion.

Involvement of sialoglycans in SARS-COV-2 infection: Opportunities and challenges for glyco-based inhibitors

Viral infections have been the causes of global pandemics, including the ongoing coronavirus disease 2019, which prompted the investigation into the infection mechanisms to find treatment and aid the vaccine design. Betacoronaviruses use spike glycoprotein on their surface to bind to host receptors, aiding their host attachment and cell fusion. Protein-glycan interaction has been implicated in the viral entry mechanism of many viruses and has recently been shown in SARS-CoV-2. Here, we reviewed the current knowledge on protein-glycan interactions that facilitate SARS-CoV-2 host entry, with special interest in sialoglycans present on both the virions and host cell surfaces. We also analyze how such information provides opportunities and challenges in glyco-based inhibitors.

SARS-CoV-2 spike opening dynamics and energetics reveal the individual roles of glycans and their collective impact

The trimeric spike (S) glycoprotein, which protrudes from the SARS-CoV-2 viral envelope, binds to human ACE2, initiated by at least one protomer's receptor binding domain (RBD) switching from a "down" (closed) to an "up" (open) state. Here, we used large-scale molecular dynamics simulations and two-dimensional replica exchange umbrella sampling calculations with more than a thousand windows and an aggregate total of 160 μs of simulation to investigate this transition with and without glycans. We find that the glycosylated spike has a higher barrier to opening and also energetically favors the down state over the up state. Analysis of the S-protein opening pathway reveals that glycans at N165 and N122 interfere with hydrogen bonds between the RBD and the N-terminal domain in the up state, while glycans at N165 and N343 can stabilize both the down and up states. Finally, we estimate how epitope exposure for several known antibodies changes along the opening path. We find that the BD-368-2 antibody's epitope is continuously exposed, explaining its high efficacy.

The Spike Mutants Website: A Worldwide Used Resource against SARS-CoV-2

A large number of SARS-CoV-2 mutations in a short period of time has driven scientific research related to vaccines, new drugs, and antibodies to combat the new variants of the virus. Herein, we present a web portal containing the structural information, the tridimensional coordinates, and the molecular dynamics trajectories of the SARS-CoV-2 spike protein and its main variants. The Spike Mutants website can serve as a rapid online tool for investigating the impact of novel mutations on virus fitness. Taking into account the high variability of SARS-CoV-2, this application can help the scientific community when prioritizing molecules for experimental assays, thus, accelerating the identification of promising drug candidates for COVID-19 treatment. Below we describe the main features of the platform and illustrate the possible applications for speeding up the drug discovery process and hypothesize new effective strategies to overcome the recurrent mutations in SARS-CoV-2 genome.

Development of a Next-Generation Vaccine Platform for Porcine Epidemic Diarrhea Virus Using a Reverse Genetics System

For the past three decades, the porcine epidemic diarrhea virus (PEDV) has remained an enormous threat to the South Korean swine industry. The scarcity of an effective method for manipulating viral genomes has impeded research progress in PEDV biology and vaccinology. Here, we report the development of reverse genetics systems using two novel infectious full-length cDNA clones of a Korean highly pathogenic-G2b strain, KNU-141112, and its live attenuated vaccine strain, S DEL5/ORF3, in a bacterial artificial chromosome (BAC) under the control of a eukaryotic promoter. Direct transfection of cells with each recombinant BAC clone induced cytopathic effects and produced infectious progeny. The reconstituted viruses, icKNU-141112 and icS DEL5/ORF3, harboring genetic markers, displayed phenotypic and genotypic properties identical to their respective parental viruses. Using the DNA-launched KNU-141112 infectious cDNA clone as a backbone, two types of recombinant viruses were generated. First, we edited the open reading frame 3 (ORF3) gene, as cell-adapted strains lose full-length ORF3, and replaced this region with an enhanced green fluorescent protein (EGFP) gene to generate icPEDV-EGFP. This mutant virus presented parental virus-like growth kinetics and stably retained robust EGFP expression, indicating that ORF3 is dispensable for PEDV replication in cell culture and is a tolerant location for exogeneous gene acceptance. However, the plaque size and syncytia phenotypes of ORF3-null icPEDV-EGFP were larger than those of icKNU-141112 but similar to ORF3-null icS DEL5/ORF3, suggesting a potential role of ORF3 in PEDV cytopathology. Second, we substituted the spike (S) gene with a heterologous S protein, designated S51, from a variant of interest (VOI), which was the most genetically and phylogenetically distant from KNU-141112. The infectious recombinant VOI, named icPEDV-S51, could be recovered, and the rescued virus showed indistinguishable growth characteristics compared to icKNU-141112. Virus cross-neutralization and structural analyses revealed antigenic differences in S between icKNU-141112 and icPEDV-S51, suggesting that genetic and conformational changes mapped within the neutralizing epitopes of S51 could impair the neutralization capacity and cause considerable immune evasion. Collectively, while the established molecular clones afford convenient, versatile platforms for PEDV genome manipulation, allowing for corroborating the molecular basis of viral replication and pathogenesis, they also provide key infrastructural frameworks for developing new vaccines and coronaviral vectors.

Distinct core glycan and O-glycoform utilization of SARS-CoV-2 Omicron variant Spike protein RBD revealed by top-down mass spectrometry

The SARS-CoV-2 Omicron (B.1.1.529) variant possesses numerous spike (S) mutations particularly in the S receptor-binding domain (S-RBD) that significantly improve transmissibility and evasion of neutralizing antibodies. But exactly how the mutations in the Omicron variant enhance viral escape from immunological protection remains to be understood. The S-RBD remains the principal target for neutralizing antibodies and therapeutics, thus new structural insights into the Omicron S-RBD and characterization of the post-translational glycosylation changes can inform rational design of vaccines and therapeutics. Here we report the molecular variations and O-glycoform changes of the Omicron S-RBD variant as compared to wild-type (WA1/2020) and Delta (B.1.617.2) variants using high-resolution top-down mass spectrometry (MS). A novel O-glycosite (Thr376) unique to the Omicron variant is identified. Moreover, we have directly quantified the Core 1 and Core 2 O-glycan structures and characterized the O-glycoform structural heterogeneity of the three variants. Our findings reveal high resolution detail of Omicron O-glycoforms and their utilization to provide direct molecular evidence of proteoform alterations in the Omicron variant which could shed light on how this variant escapes immunological protection.

Molecular dynamics studies reveal structural and functional features of the SARS-CoV-2 spike protein

The SARS-CoV-2 virus is responsible for the COVID-19 pandemic the world experience since 2019. The protein responsible for the first steps of cell invasion, the spike protein, has probably received the most attention in light of its central role during infection. Computational approaches are among the tools employed by the scientific community in the enormous effort to study this new affliction. One of these methods, namely molecular dynamics (MD), has been used to characterize the function of the spike protein at the atomic level and unveil its structural features from a dynamic perspective. In this review, we focus on these main findings, including spike protein flexibility, rare S protein conformational changes, cryptic epitopes, the role of glycans, drug repurposing, and the effect of spike protein variants.

Omicron Binding Mode: Contact Analysis and Dynamics of the Omicron Receptor-Binding Domain in Complex with ACE2

On 26 November 2021, the WHO classified the Omicron variant of the SARS-CoV-2 virus (B.1.1.529 lineage) as a variant of concern (VOC) (COVID-19 Variant Data, Department of Health, 2022). The Omicron variant contains as many as 26 unique mutations of effects not yet determined (Venkatakrishnan, A., Open Science Framework, 2021). Out of its total of 34 Spike protein mutations, 15 are located on the receptor-binding domain (S-RBD) (Stanford Coronavirus Antiviral & Resistance Database, 2022) that directly contacts the angiotensin-converting enzyme 2 (ACE2) host receptor and is also a primary target for antibodies. Here, we studied the binding mode of the S-RBD domain of the Spike protein carrying the Omicron mutations and the globular domain of human ACE2 using molecular dynamics (MD) simulations. We identified new and key Omicron-specific interactions such as R493 (of mutation Q493R), which forms salt bridges both with E35 and D38 of ACE2, Y501 (N501Y), which forms an edge-to-face aromatic interaction with Y41, and Y505 (Y505H), which makes an H-bond with E37 and K353. The glycan chains of ACE2 also bind differently in the WT and Omicron variants in response to different charge distributions on the surface of Spike proteins. However, while the Omicron mutations considerably improve the overall electrostatic fit of the two interfaces, the total number of specific and favorable interactions between the two does not increase. The dynamics of the complexes are highly affected too, making the Omicron S-RBD:ACE2 complex more rigid; the two main interaction sites, Patches I and II, isolated in the WT complex, become connected in the Omicron complex through the alternating interaction of R493 and R498 with E35 and D38.

Structural Evolution of Delta (B.1.617.2) and Omicron (BA.1) Spike Glycoproteins

The vast amount of epidemiologic and genomic data that were gathered as a global response to the COVID-19 pandemic that was caused by SARS-CoV-2 offer a unique opportunity to shed light on the structural evolution of coronaviruses and in particular on the spike (S) glycoprotein, which mediates virus entry into the host cell by binding to the human ACE2 receptor. Herein, we carry out an investigation into the dynamic properties of the S glycoprotein, focusing on the much more transmissible Delta and Omicron variants. Notwithstanding the great number of mutations that have accumulated, particularly in the Omicron S glycoprotein, our data clearly showed the conservation of some structural and dynamic elements, such as the global motion of the receptor binding domain (RBD). However, our studies also revealed structural and dynamic alterations that were concentrated in the aa 627–635 region, on a small region of the receptor binding motif (aa 483–485), and the so-called “fusion-peptide proximal region”. In particular, these last two S regions are known to be involved in the human receptor ACE2 recognition and membrane fusion. Our structural evidence, therefore, is likely involved in the observed different transmissibility of these S mutants. Finally, we highlighted the role of glycans in the increased RBD flexibility of the monomer in the up conformation of Omicron.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Point mutations in SARS-CoV-2 variants induce long-range dynamical perturbations in neutralizing antibodies

Monoclonal antibodies are emerging as a viable treatment for the coronavirus disease 19 (COVID-19). However, newly evolved variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can reduce the efficacy of currently available antibodies and can diminish vaccine-induced immunity. Here, we demonstrate that the microscopic dynamics of neutralizing monoclonal antibodies can be profoundly modified by the mutations present in the spike proteins of the SARS-COV-2 variants currently circulating in the world population. The dynamical perturbations within the antibody structure, which alter the thermodynamics of antigen recognition, are diverse and can depend both on the nature of the antibody and on the spatial location of the spike mutation. The correlation between the motion of the antibody and that of the spike receptor binding domain (RBD) can also be changed, modulating binding affinity. Using protein-graph-connectivity networks, we delineated the mutant-induced modifications in the information-flow along allosteric pathway throughout the antibody. Changes in the collective dynamics were spatially distributed both locally and across long-range distances within the antibody. On the receptor side, we identified an anchor-like structural element that prevents the detachment of the antibodies; individual mutations there can significantly affect the antibody binding propensity. Our study provides insight into how virus neutralization by monoclonal antibodies can be impacted by local mutations in the epitope via a change in dynamics. This realization adds a new layer of sophistication to the efforts for rational design of monoclonal antibodies against new variants of SARS-CoV2, taking the allostery in the antibody into consideration.

Mutations in the new variants of SARS-CoV-2 spike protein modulates the dynamics of the neutralizing antibodies. Capturing such modulations from MD simulations and graph network model identifies the role of mutations in facilitating immune evasion.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

Computational Analysis of the SARS-CoV-2 RBD-ACE2-Binding Process Based on MD and the 3D-RISM Theory

The binding process of angiotensin-converting enzyme 2 (ACE2) to the receptor-binding domain (RBD) of the severe acute respiratory syndrome-like coronavirus 2 spike protein was investigated using molecular dynamics simulation and the three-dimensional reference interaction-site model theory. The results suggested that the protein-binding process consists of a protein-protein approaching step, followed by a local structural rearrangement step. In the approaching step, the interprotein interaction energy decreased as the proteins approached each other, whereas the solvation free energy increased. As the proteins approached, the glycan of ACE2 first established a hydrogen bond with the RBD. Thereafter, the number of interprotein hydrogen bonds increased rapidly. The solvation free energy increased because of the desolvation of the protein as it approached its partner. The spatial distribution function of the solvent revealed the presence of hydrogen bonds bridged by water molecules on the RBD-ACE2 interface. Finally, principal component analysis revealed that ACE2 showed a pronounced conformational change, whereas there was no significant change in RBD.

Principles of SARS-CoV-2 Glycosylation

The structure and post-translational processing of the SARS-CoV-2 spike glycoprotein (S) is intimately associated with the function of the virus and of sterilising vaccines. The surface of the S protein is extensively modified by glycans, and their biosynthesis is driven by both the wider cellular context, and importantly, the underlining protein structure and local glycan density. Comparison of virally derived S protein with both recombinantly derived and adenovirally induced proteins, reveal hotspots of protein-directed glycosylation that drive conserved glycosylation motifs. Molecular dynamics simulations revealed that, while the S surface is extensively shielded by N-glycans, it presents regions vulnerable to neutralising antibodies. Furthermore, glycans have been shown to influence the accessibility of the receptor binding domain and the binding to the cellular receptor. The emerging picture is one of unifying, principles of S protein glycosylation and an intimate role of glycosylation in immunogen structure and efficacy.

Distinct Core Glycan and O-Glycoform Utilization of SARS-CoV-2 Omicron Variant Spike Protein RBD Revealed by Top-Down Mass Spectrometry

The SARS-CoV-2 Omicron (B.1.1.529) variant possesses numerous spike (S) mutations particularly in the S receptor-binding domain (S-RBD) that significantly improve transmissibility and evasion of neutralizing antibodies. But exactly how the mutations in the Omicron variant enhance viral escape from immunological protection remains to be understood. The S-RBD remains the principal target for neutralizing antibodies and therapeutics, thus new structural insights into the Omicron S-RBD and characterization of the post-translational glycosylation changes can inform rational design of vaccines and therapeutics. Here we report the molecular variations and O-glycoform changes of the Omicron S-RBD variant as compared to wild-type (WA1/2020) and Delta (B.1.617.2) variants using high-resolution top-down mass spectrometry (MS). A novel O-glycosite (Thr376) unique to the Omicron variant is identified. Moreover, we have directly quantified the Core 1 and Core 2 O-glycan structures and characterized the O-glycoform structural heterogeneity of the three variants. Our findings reveal high resolution detail of Omicron O-glycoforms and their utilization to provide direct molecular evidence of proteoform alterations in the Omicron variant which could shed light on how this variant escapes immunological protection.

Analysis of Viral Spike Protein N-Glycosylation Using Ultraviolet Photodissociation Mass Spectrometry

Characterization of protein glycosylation by tandem mass spectrometry remains challenging owing to the vast diversity of oligosaccharides bound to proteins, the variation in monosaccharide linkage patterns, and the lability of the linkage between the glycan and protein. Here, we have adapted an HCD-triggered-ultraviolet photodissociation (UVPD) approach for the simultaneous localization of glycosites and full characterization of both glycan compositions and intersaccharide linkages, the latter provided by extensive cross-ring cleavages enabled by UVPD. The method is applied to study glycan compositions based on analysis of glycopeptides from proteolytic digestion of recombinant human coronaviruse spike proteins from SARS-CoV-2 and HKU1. UVPD reveals unique intersaccharide linkage information and is leveraged to localize N-linked glycoforms with confidence.

GlycoGrip: Cell Surface-Inspired Universal Sensor for Betacoronaviruses

Inspired by the role of cell-surface glycoproteins as coreceptors for pathogens, we report the development of GlycoGrip: a glycopolymer-based lateral flow assay for detecting SARS-CoV-2 and its variants. GlycoGrip utilizes glycopolymers for primary capture and antispike antibodies labeled with gold nanoparticles for signal-generating detection. A lock-step integration between experiment and computation has enabled efficient optimization of GlycoGrip test strips which can selectively, sensitively, and rapidly detect SARS-CoV-2 and its variants in biofluids. Employing the power of the glycocalyx in a diagnostic assay has distinct advantages over conventional immunoassays as glycopolymers can bind to antigens in a multivalent capacity and are highly adaptable for mutated strains. As new variants of SARS-CoV-2 are identified, GlycoGrip will serve as a highly reconfigurable biosensor for their detection. Additionally, via extensive ensemble-based docking simulations which incorporate protein and glycan motion, we have elucidated important clues as to how heparan sulfate and other glycocalyx components may bind the spike glycoprotein during SARS-CoV-2 host-cell infection. GlycoGrip is a promising and generalizable alternative to costly, labor-intensive RT-PCR, and we envision it will be broadly useful, including for rural or low-income populations that are historically undertested and under-reported in infection statistics.

Structural basis for continued antibody evasion by the SARS-CoV-2 receptor binding domain

Many studies have examined the impact of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants on neutralizing antibody activity after they have become dominant strains. Here, we evaluate the consequences of further viral evolution. We demonstrate mechanisms through which the SARS-CoV-2 receptor binding domain (RBD) can tolerate large numbers of simultaneous antibody escape mutations and show that pseudotypes containing up to seven mutations, as opposed to the one to three found in previously studied variants of concern, are more resistant to neutralization by therapeutic antibodies and serum from vaccine recipients. We identify an antibody that binds the RBD core to neutralize pseudotypes for all tested variants but show that the RBD can acquire an N-linked glycan to escape neutralization. Our findings portend continued emergence of escape variants as SARS-CoV-2 adapts to humans.