Analysis of the SARS-CoV-2 spike protein glycan shield: implications for immune recognition

Differences between Omicron SARS-CoV-2 RBD and other variants in their ability to interact with cell receptors and monoclonal antibodies

SARS-CoV-2 remains a health threat with the continuous emergence of new variants. This work aims to expand the knowledge about the SARS-CoV-2 receptor-binding domain (RBD) interactions with cell receptors and monoclonal antibodies (mAbs). By using constant-pH Monte Carlo simulations, the free energy of interactions between the RBD from different variants and several partners (Angiotensin-Converting Enzyme-2 (ACE2) polymorphisms and various mAbs) were predicted. Computed RBD-ACE2-binding affinities were higher for two ACE2 polymorphisms (rs142984500 and rs4646116) typically found in Europeans which indicates a genetic susceptibility. This is amplified for Omicron (BA.1) and its sublineages BA.2 and BA.3. The antibody landscape was computationally investigated with the largest set of mAbs so far in the literature. From the 32 studied binders, groups of mAbs were identified from weak to strong binding affinities (e.g. S2K146). These mAbs with strong binding capacity and especially their combination are amenable to experimentation and clinical trials because of their high predicted binding affinities and possible neutralization potential for current known virus mutations and a universal coronavirus.Communicated by Ramaswamy H. Sarma.

Immune Response and Molecular Mechanisms of Cardiovascular Adverse Effects of Spike Proteins from SARS-CoV-2 and mRNA Vaccines

The SARS-CoV-2 (severe acute respiratory syndrome coronavirus responsible for the COVID-19 disease) uses the Spike proteins of its envelope for infecting target cells expressing on the membrane the angiotensin converting enzyme 2 (ACE2) enzyme that acts as a receptor. To control the pandemic, genetically engineered vaccines have been designed for inducing neutralizing antibodies against the Spike proteins. These vaccines do not act like traditional protein-based vaccines, as they deliver the message in the form of mRNA or DNA to host cells that then produce and expose the Spike protein on the membrane (from which it can be shed in soluble form) to alert the immune system. Mass vaccination has brought to light various adverse effects associated with these genetically based vaccines, mainly affecting the circulatory and cardiovascular system. ACE2 is present as membrane-bound on several cell types, including the mucosa of the upper respiratory and of the gastrointestinal tracts, the endothelium, the platelets, and in soluble form in the plasma. The ACE2 enzyme converts the vasoconstrictor angiotensin II into peptides with vasodilator properties. Here we review the pathways for immunization and the molecular mechanisms through which the Spike protein, either from SARS-CoV-2 or encoded by the mRNA-based vaccines, interferes with the Renin-Angiotensin-System governed by ACE2, thus altering the homeostasis of the circulation and of the cardiovascular system. Understanding the molecular interactions of the Spike protein with ACE2 and the consequent impact on cardiovascular system homeostasis will direct the diagnosis and therapy of the vaccine-related adverse effects and provide information for development of a personalized vaccination that considers pathophysiological conditions predisposing to such adverse events.

An insight into the mechanisms of COVID-19, SARS-CoV2 infection severity concerning β-cell survival and cardiovascular conditions in diabetic patients

A significantly high percentage of hospitalized COVID-19 patients with diabetes mellitus (DM) had severe conditions and were admitted to ICU. In this review, we have delineated the plausible molecular mechanisms that could explain why there are increased clinical complications in patients with DM that become critically ill when infected with SARS-CoV2. RNA viruses have been classically implicated in manifestation of new onset diabetes. SARS-CoV2 infection through cytokine storm leads to elevated levels of pro-inflammatory cytokines creating an imbalance in the functioning of T helper cells affecting multiple organs. Inflammation and Th1/Th2 cell imbalance along with Th17 have been associated with DM, which can exacerbate SARS-CoV2 infection severity. ACE-2-Ang-(1-7)-Mas axis positively modulates β-cell and cardiac tissue function and survival. However, ACE-2 receptors dock SARS-CoV2, which internalize and deplete ACE-2 and activate Renin-angiotensin system (RAS) pathway. This induces inflammation promoting insulin resistance that has positive effect on RAS pathway, causes β-cell dysfunction, promotes inflammation and increases the risk of cardiovascular complications. Further, hyperglycemic state could upregulate ACE-2 receptors for viral infection thereby increasing the severity of the diabetic condition. SARS-CoV2 infection in diabetic patients with heart conditions are linked to worse outcomes. SARS-CoV2 can directly affect cardiac tissue or inflammatory response during diabetic condition and worsen the underlying heart conditions.

Capturing a Crucial ‘Disorder-to-Order Transition’ at the Heart of the Coronavirus Molecular Pathology—Triggered by Highly Persistent, Interchangeable Salt-Bridges

The COVID-19 origin debate has greatly been influenced by genome comparison studies of late, revealing the emergence of the Furin-like cleavage site at the S1/S2 junction of the SARS-CoV-2 Spike (FLCSSpike) containing its 681PRRAR685 motif, absent in other related respiratory viruses. Being the rate-limiting (i.e., the slowest) step, the host Furin cleavage is instrumental in the abrupt increase in transmissibility in COVID-19, compared to earlier onsets of respiratory viral diseases. In such a context, the current paper entraps a ‘disorder-to-order transition’ of the FLCSSpike (concomitant to an entropy arrest) upon binding to Furin. The interaction clearly seems to be optimized for a more efficient proteolytic cleavage in SARS-CoV-2. The study further shows the formation of dynamically interchangeable and persistent networks of salt-bridges at the Spike–Furin interface in SARS-CoV-2 involving the three arginines (R682, R683, R685) of the FLCSSpike with several anionic residues (E230, E236, D259, D264, D306) coming from Furin, strategically distributed around its catalytic triad. Multiplicity and structural degeneracy of plausible salt-bridge network archetypes seem to be the other key characteristic features of the Spike–Furin binding in SARS-CoV-2, allowing the system to breathe—a trademark of protein disorder transitions. Interestingly, with respect to the homologous interaction in SARS-CoV (2002/2003) taken as a baseline, the Spike–Furin binding events, generally, in the coronavirus lineage, seems to have preference for ionic bond formation, even with a lesser number of cationic residues at their potentially polybasic FLCSSpike patches. The interaction energies are suggestive of characteristic metastabilities attributed to Spike–Furin interactions, generally to the coronavirus lineage, which appears to be favorable for proteolytic cleavages targeted at flexible protein loops. The current findings not only offer novel mechanistic insights into the coronavirus molecular pathology and evolution, but also add substantially to the existing theories of proteolytic cleavages.

Site-Specific Glycosylation Patterns of the SARS-CoV-2 Spike Protein Derived From Recombinant Protein and Viral WA1 and D614G Strains

The SARS-CoV-2 spike protein is heavily glycosylated, having 22 predicted N-glycosylation sites per monomer. It is also O-glycosylated, although the number of O-glycosites is less defined. Recent studies show that spike protein glycans play critical roles in viral entry and infection. The spike monomer has two subdomains, S1 and S2, and a receptor-binding domain (RBD) within the S1 domain. In this study, we have characterized the site-specific glycosylation patterns of the HEK293 recombinant spike RBD and S1 domains as well as the intact spike derived from the whole virus produced in Vero cells. The Vero cell-derived spike from the WA1 strain and a D614G variant was analyzed. All spike proteins, S1, and RBDs were analyzed using hydrophilic interaction chromatography (HILIC) and LC-MS/MS on an Orbitrap Eclipse Tribrid mass spectrometer. N-glycans identified in HEK293-derived S1 were structurally diverse. Those found in the HEK293-derived RBD were highly similar to those in HEK293 S1 where N-glycosites were shared. Comparison of the whole cell-derived WA1 and D614G spike proteins revealed that N-glycosites local to the mutation site appeared to be more readily detected, hinting that these sites are more exposed to glycosylation machinery. Moreover, recombinant HEK293-derived S1 was occupied almost completely with complex glycan, while both WA1 and D614G derived from the Vero E6 cell whole virus were predominantly high-mannose glycans. This stands in stark contrast to glycosylation patterns seen in both CHO- and HEK cell-derived recombinant S1, S2, and the whole spike previously reported. Concerning O-glycosylation, our analyses revealed that HEK293 recombinant proteins possessed a range of O-glycosites with compositions consistent with Core type 1 and 2 glycans. The O-glycosites shared between the S1 and RBD constructs, sites T323 and T523, were occupied by a similar range of Core 1 and 2 type O-glycans. Overall, this study reveals that the sample nature and cell substrate used for production of these proteins can have a dramatic impact on the glycosylation profile. SARS-CoV-2 spike glycans are associated with host ACE2 receptor interaction efficiency. Therefore, understanding such differences will serve to better understand these host–pathogen interactions and inform the choice of cell substrates to suite downstream investigations.

Landscape and selection of vaccine epitopes in SARS-CoV-2

BACKGROUND

Early in the pandemic, we designed a SARS-CoV-2 peptide vaccine containing epitope regions optimized for concurrent B cell, CD4+ T cell, and CD8+ T cell stimulation. The rationale for this design was to drive both humoral and cellular immunity with high specificity while avoiding undesired effects such as antibody-dependent enhancement (ADE).

METHODS

We explored the set of computationally predicted SARS-CoV-2 HLA-I and HLA-II ligands, examining protein source, concurrent human/murine coverage, and population coverage. Beyond MHC affinity, T cell vaccine candidates were further refined by predicted immunogenicity, sequence conservation, source protein abundance, and coverage of high frequency HLA alleles. B cell epitope regions were chosen from linear epitope mapping studies of convalescent patient serum, followed by filtering for surface accessibility, sequence conservation, spatial localization near functional domains of the spike glycoprotein, and avoidance of glycosylation sites.

RESULTS

From 58 initial candidates, three B cell epitope regions were identified. From 3730 (MHC-I) and 5045 (MHC-II) candidate ligands, 292 CD8+ and 284 CD4+ T cell epitopes were identified. By combining these B cell and T cell analyses, as well as a manufacturability heuristic, we proposed a set of 22 SARS-CoV-2 vaccine peptides for use in subsequent murine studies. We curated a dataset of ~ 1000 observed T cell epitopes from convalescent COVID-19 patients across eight studies, showing 8/15 recurrent epitope regions to overlap with at least one of our candidate peptides. Of the 22 candidate vaccine peptides, 16 (n = 10 T cell epitope optimized; n = 6 B cell epitope optimized) were manually selected to decrease their degree of sequence overlap and then synthesized. The immunogenicity of the synthesized vaccine peptides was validated using ELISpot and ELISA following murine vaccination. Strong T cell responses were observed in 7/10 T cell epitope optimized peptides following vaccination. Humoral responses were deficient, likely due to the unrestricted conformational space inhabited by linear vaccine peptides.

CONCLUSIONS

Overall, we find our selection process and vaccine formulation to be appropriate for identifying T cell epitopes and eliciting T cell responses against those epitopes. Further studies are needed to optimize prediction and induction of B cell responses, as well as study the protective capacity of predicted T and B cell epitopes.

Plausible blockers of Spike RBD in SARS-CoV2-molecular design and underlying interaction dynamics from high-level structural descriptors

Abstract

COVID-19 is characterized by an unprecedented abrupt increase in the viral transmission rate (SARS-CoV-2) relative to its pandemic evolutionary ancestor, SARS-CoV (2003). The complex molecular cascade of events related to the viral pathogenicity is triggered by the Spike protein upon interacting with the ACE2 receptor on human lung cells through its receptor binding domain (RBD_Spike). One potential therapeutic strategy to combat COVID-19 could thus be limiting the infection by blocking this key interaction. In this current study, we adopt a protein design approach to predict and propose non-virulent structural mimics of the RBD_Spike which can potentially serve as its competitive inhibitors in binding to ACE2. The RBD_Spike is an independently foldable protein domain, resilient to conformational changes upon mutations and therefore an attractive target for strategic re-design. Interestingly, in spite of displaying an optimal shape fit between their interacting surfaces (attributed to a consequently high mutual affinity), the RBD_Spike–ACE2 interaction appears to have a quasi-stable character due to a poor electrostatic match at their interface. Structural analyses of homologous protein complexes reveal that the ACE2 binding site of RBD_Spike has an unusually high degree of solvent-exposed hydrophobic residues, attributed to key evolutionary changes, making it inherently “reaction-prone.” The designed mimics aimed to block the viral entry by occupying the available binding sites on ACE2, are tested to have signatures of stable high-affinity binding with ACE2 (cross-validated by appropriate free energy estimates), overriding the native quasi-stable feature. The results show the apt of directly adapting natural examples in rational protein design, wherein, homology-based threading coupled with strategic “hydrophobic ↔ polar” mutations serve as a potential breakthrough.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s00894-021-04779-0.

Potent Molecular Feature-based Neutralizing Monoclonal Antibodies as Promising Therapeutics Against SARS-CoV-2 Infection

The 2019-2020 winter was marked by the emergence of a new coronavirus (SARS-CoV-2) related disease (COVID-19), which started in Wuhan, China. Its high human-to-human transmission ability led to a worldwide spread within few weeks and has caused substantial human loss. Mechanical antiviral control approach, drug repositioning, and use of COVID-19 convalescent plasmas (CPs) were the first line strategies utilized to mitigate the viral spread, yet insufficient. The urgent need to contain this deadly pandemic has led searchers and pharmaceutical companies to develop vaccines. However, not all vaccines manufactured are safe. Besides, an alternative and effective treatment option for such an infectious disease would include pure anti-viral neutralizing monoclonal antibodies (NmAbs), which can block the virus at specific molecular targets from entering cells by inhibiting virus-cell structural complex formation, with more safety and efficiency than the CP. Indeed, there is a lot of molecular evidence about the protector effect and the use of molecular feature-based NmAbs as promising therapeutics to contain COVID-19. Thus, from the scientific publication database screening, we here retrieved antibody-related papers and summarized the repertory of characterized NmAbs against SARS-CoV-2, their molecular neutralization mechanisms, and their immunotherapeutic pros and cons. About 500 anti-SARS-CoV-2 NmAbs, characterized through competitive binding assays and neutralization efficacy, were reported at the writing time (January 2021). All NmAbs bind respectively to SARS-CoV-2 S and exhibit high molecular neutralizing effects against wild-type and/or pseudotyped virus. Overall, we defined six NmAb groups blocking SARS-CoV-2 through different molecular neutralization mechanisms, from which five potential neutralization sites on SARS-CoV-2 S protein are described. Therefore, more efforts are needed to develop NmAbs-based cocktails to mitigate COVID-19.

Scrutinizing Coronaviruses Using Publicly Available Bioinformatic Tools: The Viral Structural Proteins as a Case Study

Since early 2020, the world suffers from a new beta-coronavirus, called SARS-CoV-2, that has devastating effects globally due to its associated disease, Covid-19. Until today, Covid-19, which not only causes life-threatening lung infections but also impairs various other organs and tissues, has killed hundreds of thousands of people and caused irreparable damage to many others. Since the very onset of the pandemic, huge efforts were made worldwide to fully understand this virus and numerous studies were, and still are, published. Many of these deal with structural analyses of the viral spike glycoprotein and with vaccine development, antibodies and antiviral molecules or immunomodulators that are assumed to become essential tools in the struggle against the virus. This paper summarizes knowledge on the properties of the four structural proteins (spike protein S, membrane protein M, envelope protein E and nucleocapsid protein N) of the SARS-CoV-2 virus and its relatives, SARS-CoV and MERS-CoV, that emerged few years earlier. Moreover, attention is paid to ways to analyze such proteins using freely available bioinformatic tools and, more importantly, to bring these proteins alive by looking at them on a computer/laptop screen with the easy-to-use but highly performant and interactive molecular graphics program DeepView. It is hoped that this paper will stimulate non-bioinformaticians and non-specialists in structural biology to scrutinize these and other macromolecules and as such will contribute to establishing procedures to fight these and maybe other forthcoming viruses.

The Role of the SARS-CoV-2 S-Protein Glycosylation in the Interaction of SARS-CoV-2/ACE2 and Immunological Responses

The current pandemic is caused by the coronavirus disease 2019 (COVID-19), which is, in turn, induced by a novel coronavirus (SARS-CoV-2) that triggers an acute respiratory disease. In recent years, the emergence of SARS-CoV-2 is the third highly pathogenic event and large-scale epidemic affecting the human population. It follows the severe acute respiratory syndrome coronavirus (SARS-CoV) in 2003 and the Middle East respiratory syndrome coronavirus (MERS-CoV) in 2012. This novel SARS-CoV-2 employs the angiotensin-converting enzyme 2 (ACE2) receptor, like SARS-CoV, and spreads principally in the respiratory tract. The viral spike (S) protein of coronaviruses facilities the attachment to the cellular receptor, entrance, and membrane fusion. The S protein is a glycoprotein and is critical to elicit an immune response. Glycosylation is a biologically significant post-translational modification in virus surface proteins. These glycans play important roles in the viral life cycle, structure, immune evasion, and cell infection. However, it is necessary to search for new information about viral behavior and immunological host's response after SARS-CoV-2 infection. The present review discusses the implications of the CoV-2 S protein glycosylation in the SARS-CoV-2/ACE2 interaction and the immunological response. Elucidation of the glycan repertoire on the spike protein can propel research for the development of an appropriate vaccine.

The Importance of Glycosylation in COVID-19 Infection

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is currently one of the major health problems worldwide. SARS-CoV-2 survival and virulence are shown to be impacted by glycans, covalently attached to proteins in a process of glycosylation, making glycans an area of interest in SARS-CoV-2 biology and COVID-19 infection. The SARS-CoV-2 uses its highly glycosylated spike (S) glycoproteins to bind to the cell surface receptor angiotensin-converting enzyme 2 (ACE2) glycoprotein and facilitate host cell entry. Viral glycosylation has wide-ranging roles in viral pathobiology, including mediating protein folding and stability, immune evasion, host receptor attachment, and cell entry. Modification of SARS-CoV-2 envelope membrane with glycans is important in host immune recognition and interaction between S and ACE2 glycoproteins. On the other hand, immunoglobulin G, a key molecule in immune response, shows a distinct glycosylation profile in COVID-19 infection and with increased disease severity. Hence, further studies on the role of glycosylation in SARS-CoV-2 infectivity and COVID-19 infection are needed for its successful prevention and treatment. This chapter focuses on recent findings on the importance of glycosylation in COVID-19 infection.

Pathophysiology of SARS-CoV-2 in Lung of Diabetic Patients

Novel coronavirus disease (COVID-19) is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Its impact on patients with comorbidities is clearly related to fatality cases, and diabetes has been linked to one of the most important causes of severity and mortality in SARS-CoV-2 infected patients. Substantial research progress has been made on COVID-19 therapeutics; however, effective treatments remain unsatisfactory. This unmet clinical need is robustly associated with the complexity of pathophysiological mechanisms described for COVID-19. Several key lung pathophysiological mechanisms promoted by SARS-CoV-2 have driven the response in normoglycemic and hyperglycemic subjects. There is sufficient evidence that glucose metabolism pathways in the lung are closely tied to bacterial proliferation, inflammation, oxidative stress, and pro-thrombotic responses, which lead to severe clinical outcomes. It is also likely that SARS-CoV-2 proliferation is affected by glucose metabolism of type I and type II cells. This review summarizes the current understanding of pathophysiology of SARS-CoV-2 in the lung of diabetic patients and highlights the changes in clinical outcomes of COVID-19 in normoglycemic and hyperglycemic conditions.

Molecular Basis of SARS-CoV-2 Infection and Rational Design of Potential Antiviral Agents: Modeling and Simulation Approaches

The emergence in late 2019 of the coronavirus SARS-CoV-2 has resulted in the breakthrough of the COVID-19 pandemic that is presently affecting a growing number of countries. The development of the pandemic has also prompted an unprecedented effort of the scientific community to understand the molecular bases of the virus infection and to propose rational drug design strategies able to alleviate the serious COVID-19 morbidity. In this context, a strong synergy between the structural biophysics and molecular modeling and simulation communities has emerged, resolving at the atomistic level the crucial protein apparatus of the virus and revealing the dynamic aspects of key viral processes. In this Review, we focus on how in silico studies have contributed to the understanding of the SARS-CoV-2 infection mechanism and the proposal of novel and original agents to inhibit the viral key functioning. This Review deals with the SARS-CoV-2 spike protein, including the mode of action that this structural protein uses to entry human cells, as well as with nonstructural viral proteins, focusing the attention on the most studied proteases and also proposing alternative mechanisms involving some of its domains, such as the SARS unique domain. We demonstrate that molecular modeling and simulation represent an effective approach to gather information on key biological processes and thus guide rational molecular design strategies.

Beyond Shielding: The Roles of Glycans in the SARS-CoV-2 Spike Protein

The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 28,000,000 infections and 900,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viral fusion proteins, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of the glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans and on the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift toward the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of the SARS-CoV-2 S protein, which may be exploited in the therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.

Genomic Diversity and Hotspot Mutations in 30,983 SARS-CoV-2 Genomes: Moving Toward a Universal Vaccine for the "Confined Virus"?

The COVID-19 pandemic has been ongoing since its onset in late November 2019 in Wuhan, China. Understanding and monitoring the genetic evolution of the virus, its geographical characteristics, and its stability are particularly important for controlling the spread of the disease and especially for the development of a universal vaccine covering all circulating strains. From this perspective, we analyzed 30,983 complete SARS-CoV-2 genomes from 79 countries located in the six continents and collected from 24 December 2019, to 13 May 2020, according to the GISAID database. Our analysis revealed the presence of 3206 variant sites, with a uniform distribution of mutation types in different geographic areas. Remarkably, a low frequency of recurrent mutations has been observed; only 169 mutations (5.27%) had a prevalence greater than 1% of genomes. Nevertheless, fourteen non-synonymous hotspot mutations (>10%) have been identified at different locations along the viral genome; eight in ORF1ab polyprotein (in nsp2, nsp3, transmembrane domain, RdRp, helicase, exonuclease, and endoribonuclease), three in nucleocapsid protein, and one in each of three proteins: Spike, ORF3a, and ORF8. Moreover, 36 non-synonymous mutations were identified in the receptor-binding domain (RBD) of the spike protein with a low prevalence (<1%) across all genomes, of which only four could potentially enhance the binding of the SARS-CoV-2 spike protein to the human ACE2 receptor. These results along with intra-genomic divergence of SARS-CoV-2 could indicate that unlike the influenza virus or HIV viruses, SARS-CoV-2 has a low mutation rate which makes the development of an effective global vaccine very likely.

The Answer Lies in the Energy: How Simple Atomistic Molecular Dynamics Simulations May Hold the Key to Epitope Prediction on the Fully Glycosylated SARS-CoV-2 Spike Protein

SARS-CoV-2 is a health threat with dire socioeconomical consequences. As the crucial mediator of infection, the viral glycosylated spike protein (S) has attracted the most attention and is at the center of efforts to develop therapeutics and diagnostics. Herein, we use an original decomposition approach to identify energetically uncoupled substructures as antibody binding sites on the fully glycosylated S. Crucially, all that is required are unbiased MD simulations; no prior knowledge of binding properties or ad hoc parameter combinations is needed. Our results are validated by experimentally confirmed structures of S in complex with anti- or nanobodies. We identify poorly coupled subdomains that are poised to host (several) epitopes and potentially involved in large functional conformational transitions. Moreover, we detect two distinct behaviors for glycans: those with stronger energetic coupling are structurally relevant and protect underlying peptidic epitopes, and those with weaker coupling could themselves be prone to antibody recognition.

Competitive SARS-CoV-2 Serology Reveals Most Antibodies Targeting the Spike Receptor-Binding Domain Compete for ACE2 Binding

With the emergence and continued spread of the SARS-CoV-2 virus, and of the associated disease, coronavirus disease 2019 (COVID-19), there is an urgent need for improved understanding of how the body mounts an immune response to the virus. Here, we developed a competitive SARS-CoV-2 serological assay that can simultaneously determine whether an individual has developed antibodies against the SARS-CoV-2 Spike protein receptor-binding domain (RBD) and measure the proportion of these antibodies that block interaction with the human angiotensin-converting enzyme 2 (ACE2) required for viral entry. Using this assay and 144 SARS-CoV-2 patient serum samples, we found that a majority of anti-RBD antibodies compete for ACE2 binding. These results not only highlight the need to design vaccines to generate such blocking antibodies but also demonstrate the utility of this assay to rapidly screen patient sera for potentially neutralizing antibodies.

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to spread around the world, there is an urgent need for new assay formats to characterize the humoral response to infection. Here, we present an efficient, competitive serological assay that can simultaneously determine an individual’s seroreactivity against the SARS-CoV-2 Spike protein and determine the proportion of anti-Spike antibodies that block interaction with the human angiotensin-converting enzyme 2 (ACE2) required for viral entry. In this approach based on the use of enzyme-linked immunosorbent assays (ELISA), we present natively folded viral Spike protein receptor-binding domain (RBD)-containing antigens via avidin-biotin interactions. Sera are then competed with soluble ACE2-Fc, or with a higher-affinity variant thereof, to determine the proportion of ACE2 blocking anti-RBD antibodies. Assessment of sera from 144 SARS-CoV-2 patients ultimately revealed that a remarkably consistent and high proportion of antibodies in the anti-RBD pool targeted the epitope responsible for ACE2 engagement (83% ± 11%; 50% to 107% signal inhibition in our largest cohort), further underscoring the importance of tailoring vaccines to promote the development of such antibodies.

IMPORTANCE With the emergence and continued spread of the SARS-CoV-2 virus, and of the associated disease, coronavirus disease 2019 (COVID-19), there is an urgent need for improved understanding of how the body mounts an immune response to the virus. Here, we developed a competitive SARS-CoV-2 serological assay that can simultaneously determine whether an individual has developed antibodies against the SARS-CoV-2 Spike protein receptor-binding domain (RBD) and measure the proportion of these antibodies that block interaction with the human angiotensin-converting enzyme 2 (ACE2) required for viral entry. Using this assay and 144 SARS-CoV-2 patient serum samples, we found that a majority of anti-RBD antibodies compete for ACE2 binding. These results not only highlight the need to design vaccines to generate such blocking antibodies but also demonstrate the utility of this assay to rapidly screen patient sera for potentially neutralizing antibodies.

Beyond Shielding: The Roles of Glycans in SARS-CoV-2 Spike Protein

The ongoing COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has resulted in more than 15,000,000 infections and 600,000 deaths worldwide to date. Antibody development efforts mainly revolve around the extensively glycosylated SARS-CoV-2 spike (S) protein, which mediates the host cell entry by binding to the angiotensin-converting enzyme 2 (ACE2). Similar to many other viruses, the SARS-CoV-2 spike utilizes a glycan shield to thwart the host immune response. Here, we built a full-length model of glycosylated SARS-CoV-2 S protein, both in the open and closed states, augmenting the available structural and biological data. Multiple microsecond-long, all-atom molecular dynamics simulations were used to provide an atomistic perspective on the roles of glycans, and the protein structure and dynamics. We reveal an essential structural role of N-glycans at sites N165 and N234 in modulating the conformational dynamics of the spike's receptor binding domain (RBD), which is responsible for ACE2 recognition. This finding is corroborated by biolayer interferometry experiments, which show that deletion of these glycans through N165A and N234A mutations significantly reduces binding to ACE2 as a result of the RBD conformational shift towards the "down" state. Additionally, end-to-end accessibility analyses outline a complete overview of the vulnerabilities of the glycan shield of SARS-CoV-2 S protein, which may be exploited by therapeutic efforts targeting this molecular machine. Overall, this work presents hitherto unseen functional and structural insights into the SARS-CoV-2 S protein and its glycan coat, providing a strategy to control the conformational plasticity of the RBD that could be harnessed for vaccine development.

Variations in SARS-CoV-2 Spike Protein Cell Epitopes and Glycosylation Profiles During Global Transmission Course of COVID-19

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has led to an outbreak of a pandemic worldwide. For better understanding the viral spike (S) protein variations and its potential effects on the interaction with the host immune system and also in vaccine development, the cell epitopes, glycosylation profile and their changes during the global transmission course were characterized and compared with SARS-CoV for their glycosylation profile. We analyzed totally 7,813 sequences screened from 8,897 whole genome sequences on GISAID database up to April 26, and 18 S protein amino acid variations with relatively high frequency (≥10^-3) were identified. A total of 228 sequences of variants had multiple variations, of note, most of them harboring the D614G mutation. Among the predicted 69 linear B cell epitopes, 175 discontinuous B cell epitopes and 41 cytotoxic T lymphocyte epitopes in the viral S protein, we found that the protein structure and its potential function of some sites changed, such as the linear epitope length shortened and discontinuous epitope disappeared of G476S. In addition, we detected 9 predicted N-glycosylation sites and 3 O-glycosylation sites unique to SARS-CoV-2, but no evidently observed variation of the glycan sites so far. Our findings provided an important snapshot of temporal and geographical distributions on SARS-CoV-2 S protein cell epitopes and glycosylation sites, which would be an essential basis for the selection of vaccine candidates.

Clinical and Analytical Performance of an Automated Serological Test That Identifies S1/S2-Neutralizing IgG in COVID-19 Patients Semiquantitatively

In the coronavirus (CoV) disease 2019 (COVID-19) pandemic, highly selective serological testing is essential to define exposure to severe acute respiratory syndrome CoV 2 (SARS-CoV-2). Many tests have been developed, yet with variable speeds to first results, and are of unknown quality, particularly when considering the prediction of neutralizing capacity. The LIAISON SARS-CoV-2 S1/S2 IgG assay was designed to measure antibodies against the SARS-CoV-2 native S1/S2 proteins in a standardized automated chemiluminescence assay.

In the coronavirus (CoV) disease 2019 (COVID-19) pandemic, highly selective serological testing is essential to define exposure to severe acute respiratory syndrome CoV 2 (SARS-CoV-2). Many tests have been developed, yet with variable speeds to first results, and are of unknown quality, particularly when considering the prediction of neutralizing capacity. The LIAISON SARS-CoV-2 S1/S2 IgG assay was designed to measure antibodies against the SARS-CoV-2 native S1/S2 proteins in a standardized automated chemiluminescence assay. The clinical and analytical performances of the test were validated in an observational study using residual samples (>1,500) with a positive or negative COVID-19 diagnosis. The LIAISON SARS-CoV-2 S1/S2 IgG assay proved to be highly selective and specific and offered semiquantitative measures of serum or plasma levels of anti-S1/S2 IgG with neutralizing activity. The assay’s diagnostic sensitivities were 91.3% and 95.7% at >5 or ≥15 days from diagnosis, respectively, and 100% when assessed against a neutralizing assay. The assay’s specificity ranged between 97% and 98.5%. The average imprecision of the assay was a <5% coefficient of variation. Assay performance at 2 different cutoffs was evaluated to optimize predictive values. The automated LIAISON SARS-CoV-2 S1/S2 IgG assay brings efficient, sensitive, specific, and precise serological testing to the laboratory, with the capacity to test large amounts of samples per day; first results are available within 35 min, with a throughput of 170 tests/hour. The semiquantitative results provided by the test also associate with the presence of neutralizing antibodies and may provide a useful tool for the large-scale screening of convalescent-phase plasma for safe therapeutic use.

Developing a Fully Glycosylated Full-Length SARS-CoV-2 Spike Protein Model in a Viral Membrane

This technical study describes all-atom modeling and simulation of a fully glycosylated full-length SARS-CoV-2 spike (S) protein in a viral membrane. First, starting from PDB: 6VSB and 6VXX, full-length S protein structures were modeled using template-based modeling, de-novo protein structure prediction, and loop modeling techniques in GALAXY modeling suite. Then, using the recently determined most occupied glycoforms, 22 N-glycans and 1 O-glycan of each monomer were modeled using Glycan Reader & Modeler in CHARMM-GUI. These fully glycosylated full-length S protein model structures were assessed and further refined against the low-resolution data in their respective experimental maps using ISOLDE. We then used CHARMM-GUI Membrane Builder to place the S proteins in a viral membrane and performed all-atom molecular dynamics simulations. All structures are available in CHARMM-GUI COVID-19 Archive (http://www.charmm-gui.org/docs/archive/covid19) so that researchers can use these models to carry out innovative and novel modeling and simulation research for the prevention and treatment of COVID-19.

Mapping the SARS-CoV-2 spike glycoprotein-derived peptidome presented by HLA class II on dendritic cells

Understanding and eliciting protective immune responses to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an urgent priority. To facilitate these objectives, we have profiled the repertoire of human leukocyte antigen class II (HLA-II)-bound peptides presented by HLA-DR diverse monocyte-derived dendritic cells pulsed with SARS-CoV-2 spike (S) protein. We identify 209 unique HLA-II-bound peptide sequences, many forming nested sets, which map to sites throughout S including glycosylated regions. Comparison of the glycosylation profile of the S protein to that of the HLA-II-bound S peptides revealed substantial trimming of glycan residues on the latter, likely introduced during antigen processing. Our data also highlight the receptor-binding motif in S1 as a HLA-DR-binding peptide-rich region. Results from this study have application in vaccine design, and will aid analysis of CD4+ T cell responses in infected individuals and vaccine recipients.

Epitope-based chimeric peptide vaccine design against S, M and E proteins of SARS-CoV-2, the etiologic agent of COVID-19 pandemic: an in silico approach

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiologic agent of the ongoing pandemic of coronavirus disease 2019 (COVID-19), a public health emergency of international concerns declared by the World Health Organization (WHO). An immuno-informatics approach along with comparative genomics was applied to design a multi-epitope-based peptide vaccine against SARS-CoV-2 combining the antigenic epitopes of the S, M, and E proteins. The tertiary structure was predicted, refined and validated using advanced bioinformatics tools. The candidate vaccine showed an average of ≥90.0% world population coverage for different ethnic groups. Molecular docking and dynamics simulation of the chimeric vaccine with the immune receptors (TLR3 and TLR4) predicted efficient binding. Immune simulation predicted significant primary immune response with increased IgM and secondary immune response with high levels of both IgG1 and IgG2. It also increased the proliferation of T-helper cells and cytotoxic T-cells along with the increased IFN-γ and IL-2 cytokines. The codon optimization and mRNA secondary structure prediction revealed that the chimera is suitable for high-level expression and cloning. Overall, the constructed recombinant chimeric vaccine candidate demonstrated significant potential and can be considered for clinical validation to fight against this global threat, COVID-19.

Developing a Fully-glycosylated Full-length SARS-CoV-2 Spike Protein Model in a Viral Membrane

This technical study describes all-atom modeling and simulation of a fully-glycosylated full-length SARS-CoV-2 spike (S) protein in a viral membrane. First, starting from PDB:6VSB and 6VXX, full-length S protein structures were modeled using template-based modeling, de-novo protein structure prediction, and loop modeling techniques in GALAXY modeling suite. Then, using the recently-determined most occupied glycoforms, 22 N-glycans and 1 O-glycan of each monomer were modeled using Glycan Reader & Modeler in CHARMM-GUI. These fully-glycosylated full-length S protein model structures were assessed and further refined against the low-resolution data in their respective experimental maps using ISOLDE. We then used CHARMM-GUI Membrane Builder to place the S proteins in a viral membrane and performed all-atom molecular dynamics simulations. All structures are available in CHARMM-GUI COVID-19 Archive (http://www.charmm-gui.org/docs/archive/covid19), so researchers can use these models to carry out innovative and novel modeling and simulation research for the prevention and treatment of COVID-19.

Landscape and Selection of Vaccine Epitopes in SARS-CoV-2

There is an urgent need for a vaccine with efficacy against SARS-CoV-2. We hypothesize that peptide vaccines containing epitope regions optimized for concurrent B cell, CD4+ T cell, and CD8+ T cell stimulation would drive both humoral and cellular immunity with high specificity, potentially avoiding undesired effects such as antibody-dependent enhancement (ADE). Additionally, such vaccines can be rapidly manufactured in a distributed manner. In this study, we combine computational prediction of T cell epitopes, recently published B cell epitope mapping studies, and epitope accessibility to select candidate peptide vaccines for SARS-CoV-2. We begin with an exploration of the space of possible T cell epitopes in SARS-CoV-2 with interrogation of predicted HLA-I and HLA-II ligands, overlap between predicted ligands, protein source, as well as concurrent human/murine coverage. Beyond MHC affinity, T cell vaccine candidates were further refined by predicted immunogenicity, viral source protein abundance, sequence conservation, coverage of high frequency HLA alleles and co-localization of CD4+ and CD8+ T cell epitopes. B cell epitope regions were chosen from linear epitope mapping studies of convalescent patient serum, followed by filtering to select regions with surface accessibility, high sequence conservation, spatial localization near functional domains of the spike glycoprotein, and avoidance of glycosylation sites. From 58 initial candidates, three B cell epitope regions were identified. By combining these B cell and T cell analyses, as well as a manufacturability heuristic, we propose a set of SARS-CoV-2 vaccine peptides for use in subsequent murine studies and clinical trials.