Structure and Mutations of SARS-CoV-2 Spike Protein: A Focused Overview

Changes in total charge on spike protein of SARS-CoV-2 in emerging lineages

Motivation

Charged amino acid residues on the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been shown to influence its binding to different cell surface receptors, its non-specific electrostatic interactions with the environment, and its structural stability and conformation. It is therefore important to obtain a good understanding of amino acid mutations that affect the total charge on the spike protein which have arisen across different SARS-CoV-2 lineages during the course of the virus' evolution.

Results

We analyse the change in the number of ionizable amino acids and the corresponding total charge on the spike proteins of almost 2200 SARS-CoV-2 lineages that have emerged over the span of the pandemic. Our results show that the previously observed trend toward an increase in the positive charge on the spike protein of SARS-CoV-2 variants of concern has essentially stopped with the emergence of the early omicron variants. Furthermore, recently emerged lineages show a greater diversity in terms of their composition of ionizable amino acids. We also demonstrate that the patterns of change in the number of ionizable amino acids on the spike protein are characteristic of related lineages within the broader clade division of the SARS-CoV-2 phylogenetic tree. Due to the ubiquity of electrostatic interactions in the biological environment, our findings are relevant for a broad range of studies dealing with the structural stability of SARS-CoV-2 and its interactions with the environment.

Availability and implementation

The data underlying the article are available in the Supplementary material.

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

RATIONALE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Hydroxypropyl-Beta Cyclodextrin Barrier Prevents Respiratory Viral Infections: A Preclinical Study

The emergence and mutation of pathogenic viruses have been occurring at an unprecedented rate in recent decades. The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has developed into a global public health crisis due to extensive viral transmission. In situ RNA mapping has revealed angiotensin-converting enzyme 2 (ACE2) expression to be highest in the nose and lower in the lung, pointing to nasal susceptibility as a predominant route for infection and the cause of subsequent pulmonary effects. By blocking viral attachment and entry at the nasal airway using a cyclodextrin-based formulation, a preventative therapy can be developed to reduce viral infection at the site of entry. Here, we assess the safety and antiviral efficacy of cyclodextrin-based formulations. From these studies, hydroxypropyl beta-cyclodextrin (HPBCD) and hydroxypropyl gamma-cyclodextrin (HPGCD) were then further evaluated for antiviral effects using SARS-CoV-2 pseudotypes. Efficacy findings were confirmed with SARS-CoV-2 Delta variant infection of Calu-3 cells and using a K18-hACE2 murine model. Intranasal pre-treatment with HPBCD-based formulations reduced viral load and inflammatory signaling in the lung. In vitro efficacy studies were further conducted using lentiviruses, murine hepatitis virus (MHV), and influenza A virus subtype H1N1. These findings suggest HPBCD may be used as an agnostic barrier against transmissible pathogens, including but not limited to SARS-CoV-2.

Predicting virus Fitness: Towards a structure-based computational model

Predicting the impact of new emerging virus mutations is of major interest in surveillance and for understanding the evolutionary forces of the pathogens. The SARS-CoV-2 surface spike-protein (S-protein) binds to human ACE2 receptors as a critical step in host cell infection. At the same time, S-protein binding to human antibodies neutralizes the virus and prevents interaction with ACE2. Here we combine these two binding properties in a simple virus fitness model, using structure-based computation of all possible mutation effects averaged over 10 ACE2 complexes and 10 antibody complexes of the S-protein (∼380,000 computed mutations), and validated the approach against diverse experimental binding/escape data of ACE2 and antibodies. The ACE2-antibody selectivity change caused by mutation (i.e., the differential change in binding to ACE2 vs. immunity-inducing antibodies) is proposed to be a key metric of fitness model, enabling systematic error cancelation when evaluated. In this model, new mutations become fixated if they increase the selective binding to ACE2 relative to circulating antibodies, assuming that both are present in the host in a competitive binding situation. We use this model to categorize viral mutations that may best reach ACE2 before being captured by antibodies. Our model may aid the understanding of variant-specific vaccines and molecular mechanisms of viral evolution in the context of a human host.

Novel Polymyxin-Inspired Peptidomimetics Targeting the SARS-CoV-2 Spike:hACE2 Interface

Though the bulk of the COVID-19 pandemic is behind, the search for effective and safe anti-SARS-CoV-2 drugs continues to be relevant. A highly pursued approach for antiviral drug development involves targeting the viral spike (S) protein of SARS-CoV-2 to prevent its attachment to the cellular receptor ACE2. Here, we exploited the core structure of polymyxin B, a naturally occurring antibiotic, to design and synthesize unprecedented peptidomimetics (PMs), intended to target contemporarily two defined, non-overlapping regions of the S receptor-binding domain (RBD). Monomers 1, 2, and 8, and heterodimers 7 and 10 bound to the S-RBD with micromolar affinity in cell-free surface plasmon resonance assays (K_D ranging from 2.31 μM to 2.78 μM for dimers and 8.56 μM to 10.12 μM for monomers). Although the PMs were not able to fully protect cell cultures from infection with authentic live SARS-CoV-2, dimer 10 exerted a minimal but detectable inhibition of SARS-CoV-2 entry in U87.ACE2⁺ and A549.ACE2.TMPRSS2⁺ cells. These results validated a previous modeling study and provided the first proof-of-feasibility of using medium-sized heterodimeric PMs for targeting the S-RBD. Thus, heterodimers 7 and 10 may serve as a lead for the development of optimized compounds, which are structurally related to polymyxin, with improved S-RBD affinity and anti-SARS-CoV-2 potential.

Topological data analysis of protein structure and inter/intra-molecular interaction changes attributable to amino acid mutations

The presence of some amino acid mutations in the amino acid sequence that determines a protein's structure can significantly affect that 3D structure and its biological function. However, the effects upon structural and functional changes differ for each displaced amino acid, and it is very difficult to predict these changes in advance. Although computer simulations are very effective at predicting conformational changes, they struggle to determine whether the amino acid mutation of interest induces sufficient conformational changes, unless the researcher is a specialist in molecular structure calculations. Therefore, we created a framework that efficiently utilizes molecular dynamics and persistent homology methods to identify amino acid mutations that induce structural changes. We show that this framework can be used not only to predict conformational changes produced by amino acid mutations but also to extract groups of mutations that significantly alter similar molecular interactions, by capturing the resultant protein-protein interaction changes.

Abiotic Synthetic Antibody Inhibitor with Broad-Spectrum Neutralization and Antiviral Efficacy against Escaping SARS-CoV-2 Variants

The rapid emergence and spread of vaccine/antibody-escaping variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posed serious challenges to our efforts in combating corona virus disease 2019 (COVID-19) pandemic. A potent and broad-spectrum neutralizing reagent against these escaping mutants is extremely important for the development of strategies for the prevention and treatment of SARS-CoV-2 infection. We herein report an abiotic synthetic antibody inhibitor as a potential anti-SARS-CoV-2 therapeutic agent. The inhibitor, Aphe-NP14, was selected from a synthetic hydrogel polymer nanoparticle library created by incorporating monomers with functionalities complementary to key residues of the SARS-CoV-2 spike glycoprotein receptor binding domain (RBD) involved in human angiotensin-converting enzyme 2 (ACE2) binding. It has high capacity, fast adsorption kinetics, strong affinity, and broad specificity in biologically relevant conditions to both the wild type and the current variants of concern, including Beta, Delta, and Omicron spike RBD. The Aphe-NP14 uptake of spike RBD results in strong blockage of spike RBD-ACE2 interaction and thus potent neutralization efficacy against these escaping spike protein variant pseudotyped viruses. It also inhibits live SARS-CoV-2 virus recognition, entry, replication, and infection in vitro and in vivo. The Aphe-NP14 intranasal administration is found to be safe due to its low in vitro and in vivo toxicity. These results establish a potential application of abiotic synthetic antibody inhibitors in the prevention and treatment of the infection of emerging or possibly future SARS-CoV-2 variants.

Exploring the Potential of Broadly Neutralizing Antibodies for Treating SARS-CoV-2 Variants of Global Concern in 2023: A Comprehensive Clinical Review

In the aftermath of the coronavirus disease 2019 (COVID-19) pandemic, the world is still seeing outbreaks of COVID-19 infections as of 2023, especially in populations that have been adequately vaccinated. This situation across the globe gives rise to important questions regarding the efficacy of current treatments and the real rate of mutations in the COVID-19 virus itself which can make the currently available treatments and vaccines obsolete. We have tried to answer a few of those questions and put forth some new questions of our own. Our efforts in this paper were directed towards understanding the utilization of broadly neutralizing antibodies as a treatment for COVID-19 infection with a particular focus on the Omicron variant and other newer variants. We gathered our data from three major databases: PubMed, Google Scholar, and Cochrane Central Register of Controlled Trials (CENTRAL). We have screened 7070 studies from inception till March 5, 2023, and gathered 63 articles that were relevant to the topic of interest. Based on the existing medical literature regarding the topic of interest and also based on our own personal and clinical experience treating COVID-19 patients across the multiple waves in the United States and India since the beginning of the pandemic, we have concluded that broad neutralizing antibodies could be an effective option for treatment and prophylaxis for current and future outbreaks of COVID-19 including the Omicron variant and newer variants. Further research, including clinical trials, is required to tailor optimal dosages, prevent adverse reactions/side effects, and develop treatment strategies.

Antitarget, Anti-SARS-CoV-2 Leads, Drugs, and the Drug Discovery-Genetics Alliance Perspective

The most advanced antiviral molecules addressing major SARS-CoV-2 targets (Main protease, Spike protein, and RNA polymerase), compared with proteins of other human pathogenic coronaviruses, may have a short-lasting clinical efficacy. Accumulating knowledge on the mechanisms underlying the target structural basis, its mutational progression, and the related biological significance to virus replication allows envisaging the development of better-targeted therapies in the context of COVID-19 epidemic and future coronavirus outbreaks. The identification of evolutionary patterns based solely on sequence information analysis for those targets can provide meaningful insights into the molecular basis of host-pathogen interactions and adaptation, leading to drug resistance phenomena. Herein, we will explore how the study of observed and predicted mutations may offer valuable suggestions for the application of the so-called "synthetic lethal" strategy to SARS-CoV-2 Main protease and Spike protein. The synergy between genetics evidence and drug discovery may prioritize the development of novel long-lasting antiviral agents.

Evolution of increased positive charge on the SARS-CoV-2 spike protein may be adaptation to human transmission

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve and infect individuals. The exterior surface of the SARS-CoV-2 virion is dominated by the spike protein, and the current work examined spike protein biochemical features that have changed during the 3 years in which SARS-CoV-2 has infected humans. Our analysis identified a striking change in spike protein charge, from -8.3 in the original Lineage A and B viruses to -1.26 in most of the current Omicron viruses. We conclude that in addition to immune selection pressure, the evolution of SARS-CoV-2 has also altered viral spike protein biochemical properties, which may influence virion survival and promote transmission. Future vaccine and therapeutic development should also exploit and target these biochemical properties.

Modelling SARS-CoV-2 spike-protein mutation effects on ACE2 binding

The binding affinity of the SARS-CoV-2 spike (S)-protein to the human membrane protein ACE2 is critical for virus function. Computational structure-based screening of new S-protein mutations for ACE2 binding lends promise to rationalize virus function directly from protein structure and ideally aid early detection of potentially concerning variants. We used a computational protocol based on cryo-electron microscopy structures of the S-protein to estimate the change in ACE2-affinity due to S-protein mutation (ΔΔG_bind) in good trend agreement with experimental ACE2 affinities. We then expanded predictions to all possible S-protein mutations in 21 different S-protein-ACE2 complexes (400,000 ΔΔG_bind data points in total), using mutation group comparisons to reduce systematic errors. The results suggest that mutations that have arisen in major variants as a group maintain ACE2 affinity significantly more than random mutations in the total protein, at the interface, and at evolvable sites. Omicron mutations as a group had a modest change in binding affinity compared to mutations in other major variants. The single-mutation effects seem consistent with ACE2 binding being optimized and maintained in omicron, despite increased importance of other selection pressures (antigenic drift), however, epistasis, glycosylation and in vivo conditions will modulate these effects. Computational prediction of SARS-CoV-2 evolution remains far from achieved, but the feasibility of large-scale computation is substantially aided by using many structures and mutation groups rather than single mutation effects, which are very uncertain. Our results demonstrate substantial challenges but indicate ways forward to improve the quality of computer models for assessing SARS-CoV-2 mutation effects.

Binding of human ACE2 and RBD of Omicron enhanced by unique interaction patterns among SARS-CoV-2 variants of concern

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing COVID-19, has continued to mutate and spread worldwide despite global vaccination efforts. In particular, the Omicron variant, first identified in South Africa in late November 2021, has become the dominant strain worldwide. Compared to the original strain identified in Wuhan, Omicron features 50 genetic mutations, with 15 mutations in the receptor-binding domain (RBD) of the spike protein, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor for viral entry. However, it is not completely understood how these mutations alter the interaction and binding strength between the Omicron RBD and ACE2. In this study, we used a combined steered molecular dynamics (SMD) simulation and experimental microscale thermophoresis (MST) approach to quantify the interaction between Omicron RBD and ACE2. We report that the Omicron brings an enhanced RBD-ACE2 interface through N501Y, Q498R, and T478K mutations; the changes further lead to unique interaction patterns, reminiscing the features of previously dominated variants, Alpha (N501Y) and Delta (L452R and T478K). Among the Q493K and Q493R, we report that Q493R shows stronger binding to ACE2 than Q493K due to increased interactions. Our MST data confirmed that the Omicron mutations in RBD are associated with a five-fold higher binding affinity to ACE2 compared to the RBD of the original strain. In conclusion, our results could help explain the Omicron variant's prevalence in human populations, as higher interaction forces or affinity for ACE2 likely promote greater viral binding and internalization, leading to increased infectivity.

Stability and expression of SARS-CoV-2 spike-protein mutations

Protein fold stability likely plays a role in SARS-CoV-2 S-protein evolution, together with ACE2 binding and antibody evasion. While few thermodynamic stability data are available for S-protein mutants, many systematic experimental data exist for their expression. In this paper, we explore whether such expression levels relate to the thermodynamic stability of the mutants. We studied mutation-induced SARS-CoV-2 S-protein fold stability, as computed by three very distinct methods and eight different protein structures to account for method- and structure-dependencies. For all methods and structures used (24 comparisons), computed stability changes correlate significantly (99% confidence level) with experimental yeast expression from the literature, such that higher expression is associated with relatively higher fold stability. Also significant, albeit weaker, correlations were seen between stability and ACE2 binding effects. The effect of thermodynamic fold stability may be direct or a correlate of amino acid or site properties, notably the solvent exposure of the site. Correlation between computed stability and experimental expression and ACE2 binding suggests that functional properties of the SARS-CoV-2 S-protein mutant space are largely determined by a few simple features, due to underlying correlations. Our study lends promise to the development of computational tools that may ideally aid in understanding and predicting SARS-CoV-2 S-protein evolution.

Membrane attachment and fusion of HIV-1, influenza A, and SARS-CoV-2: resolving the mechanisms with biophysical methods

Attachment to and fusion with cell membranes are two major steps in the replication cycle of many human viruses. We focus on these steps for three enveloped viruses, i.e., HIV-1, IAVs, and SARS-CoV-2. Viral spike proteins drive the membrane attachment and fusion of these viruses. Dynamic interactions between the spike proteins and membrane receptors trigger their specific attachment to the plasma membrane of host cells. A single virion on cell membranes can engage in binding with multiple receptors of the same or different types. Such dynamic and multivalent binding of these viruses result in an optimal attachment strength which in turn leads to their cellular entry and membrane fusion. The latter process is driven by conformational changes of the spike proteins which are also class I fusion proteins, providing the energetics of membrane tethering, bending, and fusion. These viruses exploit cellular and membrane factors in regulating the conformation changes and membrane processes. Herein, we describe the major structural and functional features of spike proteins of the enveloped viruses including highlights on their structural dynamics. The review delves into some of the case studies in the literature discussing the findings on multivalent binding, membrane hemifusion, and fusion of these viruses. The focus is on applications of biophysical tools with an emphasis on single-particle methods for evaluating mechanisms of these processes at the molecular and cellular levels.

Structural heterogeneity and precision of implications drawn from cryo-electron microscopy structures: SARS-CoV-2 spike-protein mutations as a test case

Protein structures may be used to draw functional implications at the residue level, but how sensitive are these implications to the exact structure used? Calculation of the effects of SARS-CoV-2 S-protein mutations based on experimental cryo-electron microscopy structures have been abundant during the pandemic. To understand the precision of such estimates, we studied three distinct methods to estimate stability changes for all possible mutations in 23 different S-protein structures (3.69 million ΔΔG values in total) and explored how random and systematic errors can be remedied by structure-averaged mutation group comparisons. We show that computational estimates have low precision, due to method and structure heterogeneity making results for single mutations uninformative. However, structure-averaged differences in mean effects for groups of substitutions can yield significant results. Illustrating this protocol, functionally important natural mutations, despite individual variations, average to a smaller stability impact compared to other possible mutations, independent of conformational state (open, closed). In summary, we document substantial issues with precision in structure-based protein modeling and recommend sensitivity tests to quantify these effects, but also suggest partial solutions to the problem in the form of structure-averaged “ensemble” estimates for groups of residues when multiple structures are available.

In silico analysis of mutations near S1/S2 cleavage site in SARS-CoV-2 spike protein reveals increased propensity of glycosylation in Omicron strain

Cleavage of the severe respiratory syndrome coronavirus-2 (SARS-CoV-2) spike protein has been demonstrated to contribute to viral-cell fusion and syncytia formation. Studies have shown that variants of concern (VOC) and variants of interest (VOI) show differing membrane fusion capacity. Mutations near cleavage motifs, such as the S1/S2 and S2' sites, may alter interactions with host proteases and, thus, the potential for fusion. The biochemical basis for the differences in interactions with host proteases for the VOC/VOI spike proteins has not yet been explored. Using sequence and structure-based bioinformatics, mutations near the VOC/VOI spike protein cleavage sites were inspected for their structural effects. All mutations found at the S1/S2 sites were predicted to increase affinity to the furin protease but not TMPRSS2. Mutations at the spike residue P681 in several strains, such P681R in the Delta strain, resulted in the disruption of a proline-directed kinase phosphorylation motif at the S1/S2 site, which may lessen the impact of phosphorylation for these variants. However, the unique N679K mutation in the Omicron strain was found to increase the propensity for O-linked glycosylation at the S1/S2 cleavage site, which may prevent recognition by proteases. Such glycosylation in the Omicron strain may hinder entry at the cell surface and, thus, decrease syncytia formation and induce cell entry through the endocytic pathway as has been shown in previous studies. Further experimental work is needed to confirm the effect of mutations and posttranslational modifications on SARS-CoV-2 spike protein cleavage sites.

Computational Analysis of Short Linear Motifs in the Spike Protein of SARS-CoV-2 Variants Provides Possible Clues into the Immune Hijack and Evasion Mechanisms of Omicron Variant

Short linear motifs (SLiMs) are short linear sequences that can mediate protein–protein interaction. Mimicking eukaryotic SLiMs to compete with extra- or intracellular binding partners, or to sequester host proteins is the crucial strategy of viruses to pervert the host system. Evolved proteins in viruses facilitate minimal protein–protein interactions that significantly affect intracellular signaling networks. Unfortunately, very little information about SARS-CoV-2 SLiMs is known, especially across SARS-CoV-2 variants. Through the ELM database-based sequence analysis of spike proteins from all the major SARS-CoV-2 variants, we identified four overriding SLiMs in the SARS-CoV-2 Omicron variant, namely, LIG_TRFH_1, LIG_REV1ctd_RIR_1, LIG_CaM_NSCaTE_8, and MOD_LATS_1. These SLiMs are highly likely to interfere with various immune functions, interact with host intracellular proteins, regulate cellular pathways, and lubricate viral infection and transmission. These cellular interactions possibly serve as potential therapeutic targets for these variants, and this approach can be further exploited to combat emerging SARS-CoV-2 variants.

Review of concerned SARS-CoV-2 variants like Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529), as well as novel methods for reducing and inactivating SARS-CoV-2 mutants in wastewater treatment facilities

The coronavirus known as COVID-19, which causes pandemics, is causing a global epidemic at a critical stage today. Furthermore, novel mutations in the SARS-CoV-2 spike protein have been discovered in an entirely new strain, impacting the clinical and epidemiological features of COVID-19. Variants of these viruses can increase the transmission in wastewater, lead to reinfection, and reduce immunity provided by monoclonal antibodies and vaccinations. According to the research, a large quantity of viral RNA was discovered in wastewater, suggesting that wastewater can be a crucial source of epidemiological data and health hazards. The purpose of this paper is to introduce a few basic concepts regarding wastewater surveillance as a starting point for comprehending COVID-19′s epidemiological aspects. Next, the observation of Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529) in wastewater is discussed in detail. Secondly, the essential information for the initial, primary, and final treating sewage in SARS-CoV-2 is introduced. Following that, a thorough examination is provided to highlight the newly developed methods for eradicating SARS-CoV-2 using a combination of solar water disinfection (SODIS) and ultraviolet radiation A (UVA (315-400 nm)), ultraviolet radiation B (UVB (280-315 nm)), and ultraviolet radiation C (UVC (100-280 nm)) processes. SARS-CoV-2 eradication requires high temperatures (above 56°C) and UVC. However, SODIS technologies are based on UVA and operate at cooler temperatures (less than 45°C). Hence, it is not appropriate for sewage treatment (or water consumption) to be conducted using SODIS methods in the current pandemic. Finally, SARS-CoV-2 may be discovered in sewage utilizing the wastewater-based epidemiology (WBE) monitoring method.

Omicron Spike Protein Has a Positive Electrostatic Surface That Promotes ACE2 Recognition and Antibody Escape

High transmissibility is a hallmark of the Omicron variant of SARS-CoV-2. Understanding the molecular determinants of Omicron’s transmissibility will impact development of intervention strategies. Here we map the electrostatic potential surface of the Spike protein to show that major SARS-CoV-2 variants have accumulated positive charges in solvent-exposed regions of the Spike protein, especially its ACE2-binding interface. Significantly, the Omicron Spike-ACE2 complex has complementary electrostatic surfaces. In contrast, interfaces between Omicron and neutralizing antibodies tend to have similar positively charged surfaces. Structural modeling demonstrates that the electrostatic property of Omicron’s Spike receptor binding domain (S RBD) plays a role in enhancing ACE2 recognition and destabilizing Spike-antibody complexes. Specifically, the Omicron S RBD has favorable electrostatic interaction energy with ACE2 that is 3-5 times greater than the Delta variant over a range of 20 Å, implying efficient recognition of host receptors. Computed binding affinities of six representative S RBD-antibody complexes show that Omicron can escape most antibodies targeting the ACE2-binding region of S RBD. Interestingly, a straightforward assessment of the electrostatic surfaces of 18 neutralizing antibodies correctly predicted the Omicron escape status of 80% of cases. Collectively, our structural analysis implies that Omicron S RBD interaction interfaces have been optimized to simultaneously promote access to human ACE2 receptors and evade antibodies. These findings suggest that electrostatic interactions are a major contributing factor for increased Omicron transmissibility relative to other variants.

Changing Dynamics of SARS-CoV-2: A Global Challenge

Since November 2019, SARS-CoV-2 has been a matter of global concern due to its rapid spread, the millions of deaths it caused, and repeated waves of infections. One after another, many variants of this novel virus have come into existence due to its constant mutability, specifically in the spike glycoprotein region. The tally for variants of concern (VOCs), which already include Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1) and Delta (B.1.617.2), has increased to five with the latest appearance of Omicron (B.1.1.529). In our study, we examine the effect of the transmissibility and infectious potential of the virus due to various mutations of SARS-CoV-2, especially in the receptor-binding domain (RBD). We discuss the role of genome sequencing in tracing all the mutations and the importance of the R value (reproductive number) to understand the virus spread. We also review the effectiveness of the available vaccines on the variants of concern, as the rapid spread of the newly emergent Omicron variant has raised doubts about the usefulness of the current vaccines. The use of a mixed vaccination strategy has proved to be effective, yet the newer variants, such as Omicron, demand booster doses for the population. Multivalent immunogens could be considered as the plausible solution for conferring protection against potential new mutants of the virus in the future.

An Insight Based on Computational Analysis of the Interaction between the Receptor-Binding Domain of the Omicron Variants and Human Angiotensin-Converting Enzyme 2

Simple Summary

The Omicron variant has recently been divided into BA.1, BA.2, and BA.3 subvariants. In the present study, we focused on comparing the interaction between the receptor-binding domain (RBD) of BA.1 and BA.2 spike proteins with human angiotensin-converting enzyme 2 (hACE2) using a computational approach. The RBD BA.2 was modeled after the BA.1. The results from molecular docking and molecular dynamics studies showed that RBD BA.2 has a higher and more stable affinity for hACE2 compared to RBD BA.1.

Abstract

Concerns have been raised about the high number of mutations in the spike protein of the new emergence of the highly transmissible Omicron variant (B.1.1529 lineage) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This variant’s extraordinary ability to evade antibodies would significantly impair the current vaccination program. This present study aimed to computationally analyze the interaction between the receptor-binding domain (RBD) in the spike protein of Omicron variants and human angiotensin-converting enzyme 2 (hACE2). The docking results indicated that Omicron BA.2 has exceptionally strong interactions with hACE2 in comparison to Omicron BA.1, Delta, and wild-type, as indicated by various parameters such as salt bridge, hydrogen bond, and non-bonded interactions. The results of the molecular dynamics simulation study corroborate these findings, indicating that Omicron BA.2 has a strong and stable interaction with hACE2. This study provides insight into the development of an effective intervention against this variant.

Binding of Human ACE2 and RBD of Omicron Enhanced by Unique Interaction Patterns Among SARS-CoV-2 Variants of Concern

The 2019 coronavirus disease (COVID-19) pandemic has had devastating impacts on our global health. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing COVID-19, has continued to mutate and spread worldwide despite global vaccination efforts. In particular, the Omicron variant, first identified in South Africa in late November 2021, has now overtaken the Delta variant and become the dominant strain worldwide. Compared to the original strain identified in Wuhan, Omicron features 50 genetic mutations, with 15 mutations in the receptor-binding domain (RBD) of the spike protein, which binds to the human angiotensin-converting enzyme 2 (ACE2) receptor for viral entry. However, it is not completely understood how these mutations alter the interaction and binding strength between the Omicron RBD and ACE2. In this study, we used a combined steered molecular dynamics (SMD) simulation and experimental microscale thermophoresis (MST) approach to quantify the interaction between Omicron RBD and ACE2. We report that the Omicron brings an enhanced RBD-ACE2 interface through N501Y, Q493K/R, and T478K mutations; the changes further lead to unique interaction patterns, reminiscing the features of previously dominated variants, Alpha (N501Y) and Delta (L452R and T478K). Our MST data confirmed that the Omicron mutations in RBD are associated with a five-fold higher binding affinity to ACE2 compared to the RBD of the original strain. In conclusion, our result could help explain the Omicron variant’s prevalence in human populations, as higher interaction forces or affinity for ACE2 likely promote greater viral binding and internalization, leading to increased infectivity.

TOC GRAPHIC