Title: Impact of Missense Mutations on Spike Protein Stability and Binding Affinity in the Omicron Variant
The global effort to combat the COVID-19 pandemic faces ongoing uncertainty with the emergence of Variants of Concern featuring numerous mutations on the Spike (S) protein. In particular, the Omicron Variant is distinguished by 32 mutations, including 10 within its receptor-binding domain (RBD). These mutations significantly impact viral infectivity and the efficacy of vaccines and antibodies currently in use for therapeutic purposes. In our study, we employed structure-based computational saturation mutagenesis approaches to predict the effects of Omicron missense mutations on RBD stability and binding affinity, comparing them to the original Wuhan-Hu-1 strain. Our results predict that mutations such as G431W and P507W induce the most substantial destabilizations in the Wuhan-Hu-1-S/Omicron-S RBD. Notably, we postulate that mutations in the Omicron-S exhibit a higher percentage of enhancing binding affinity compared to Wuhan-S. We found that the mutations at residue positions G447, Y449, F456, F486, and S496 led to significant changes in binding affinity. In summary, our findings may shed light on the widespread prevalence of Omicron mutations in human populations. The Omicron mutations that potentially enhance their affinity for human receptors may facilitate increased viral binding and internalization in infected cells, thereby enhancing infectivity. This informs the development of new neutralizing antibodies capable of targeting Omicron’s immune-evading mutations, potentially aiding in the ongoing battle against the COVID-19 pandemic.  more » « less
Award ID(s):
2000296
PAR ID:
10595219
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
MDPI
Date Published:
Journal Name:
Viruses
Volume:
16
Issue:
7
ISSN:
1999-4915
Page Range / eLocation ID:
1150
Subject(s) / Keyword(s):
SARS-CoV-2 COVID-19 Spike protein Omicron Variant computational saturation mutagenesis
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Liu, Shan-Lu (Ed.)
    ABSTRACT The human angiotensin-converting enzyme 2 (hACE2) is the primary receptor for the entry of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Some human alleles of ACE2 exhibit an improved affinity for the SARS-CoV-2 Spike protein. However, the impact of ACE2 polymorphisms on SARS-CoV-2 infection remains unclear. Our previous study predicted that G431 and S514 in the receptor-binding domain (RBD) of SARS-CoV-2 S1 domain are important for S protein stability, and that S protein residues G496 and F497 and ACE2 residues D355 and Y41 are critical for the RBD-ACE2 interaction. In this study, we explored the potential of hACE2-derived neutralizing peptides as a therapeutic strategy against SARS-CoV-2 and investigated how ACE2 polymorphisms affect RBD-ACE2 binding affinity. We applied computational saturation mutagenesis to systematically screen the binding affinity changes among all possible ACE2 missense mutations within the ACE2-Wuhan-S1 complex. Mutations at ACE2 residues D355 and Y41 were predicted to weaken binding affinity, whereas those at N330 and D30 enhanced it. We identified six ACE2 regions (19–49, 65–102, 320–333, 348–359, 378–395, and 552–563) to be vital for ACE2-RBD interaction. We synthesized peptides corresponding to these six regions and tested them using a pseudotyped viral particle system and dot blot assay. Three peptides were confirmed to bind with the S protein, and four exhibited inhibitory effects. We aligned ACE2-Wuhan-S1 and ACE2-Omicron-S1 complexes, conducted correlation analysis, and observed similar binding patterns, suggesting that these peptides also have the potential to neutralize Omicron strains.IMPORTANCESARS-CoV-2 continues its global spread. In this research, we identified six regions within ACE2 that are vital for interaction with the viral S receptor-binding domain and have the potential to neutralize SARS-CoV-2 infection. Among the six peptides derived from ACE2, three were confirmed to bind with the S protein of the Wuhan strain, and four exhibited inhibitory effects on the Wuhan strain SARS-CoV-2. We also found ACE2 residues D355 and Y41 as weakening affinity, and N330 and D30 as enhancing it. We also aligned this complex with the ACE2-Omicron-S1 complex, performed correlation analyses, and compared their patterns of stability changes upon mutations and obtained similar results, indicating that these peptides may also be effective against Omicron variants. These results provide insight into the role of ACE2 polymorphism in viral entry and suggest that hACE2-derived peptides may offer a promising therapeutic strategy against SARS-CoV-2, demonstrating strong consistency between our computational predictions and experimental outcomes. 
    more » « less
  2. Abstract The glycosylation on the spike (S) protein of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus that causes COVID-19, modulates the viral infection by altering conformational dynamics, receptor interaction and host immune responses. Several variants of concern (VOCs) of SARS-CoV-2 have evolved during the pandemic, and crucial mutations on the S protein of the virus have led to increased transmissibility and immune escape. In this study, we compare the site-specific glycosylation and overall glycomic profiles of the wild type Wuhan-Hu-1 strain (WT) S protein and five VOCs of SARS-CoV-2: Alpha, Beta, Gamma, Delta and Omicron. Interestingly, both N- and O-glycosylation sites on the S protein are highly conserved among the spike mutant variants, particularly at the sites on the receptor-binding domain (RBD). The conservation of glycosylation sites is noteworthy, as over 2 million SARS-CoV-2 S protein sequences have been reported with various amino acid mutations. Our detailed profiling of the glycosylation at each of the individual sites of the S protein across the variants revealed intriguing possible association of glycosylation pattern on the variants and their previously reported infectivity. While the sites are conserved, we observed changes in the N- and O-glycosylation profile across the variants. The newly emerged variants, which showed higher resistance to neutralizing antibodies and vaccines, displayed a decrease in the overall abundance of complex-type glycans with both fucosylation and sialylation and an increase in the oligomannose-type glycans across the sites. Among the variants, the glycosylation sites with significant changes in glycan profile were observed at both theN-terminal domain and RBD of S protein, with Omicron showing the highest deviation. The increase in oligomannose-type happens sequentially from Alpha through Delta. Interestingly, Omicron does not contain more oligomannose-type glycans compared to Delta but does contain more compared to the WT and other VOCs. O-glycosylation at the RBD showed lower occupancy in the VOCs in comparison to the WT. Our study on the sites and pattern of glycosylation on the SARS-CoV-2 S proteins across the VOCs may help to understand how the virus evolved to trick the host immune system. Our study also highlights how the SARS-CoV-2 virus has conserved bothN- andO- glycosylation sites on the S protein of the most successful variants even after undergoing extensive mutations, suggesting a correlation between infectivity/ transmissibility and glycosylation. 
    more » « less
  3. Abstract The Omicron BA.1 variant emerged in late 2021 and quickly spread across the world. Compared to the earlier SARS-CoV-2 variants, BA.1 has many mutations, some of which are known to enable antibody escape. Many of these antibody-escape mutations individually decrease the spike receptor-binding domain (RBD) affinity for ACE2, but BA.1 still binds ACE2 with high affinity. The fitness and evolution of the BA.1 lineage is therefore driven by the combined effects of numerous mutations. Here, we systematically map the epistatic interactions between the 15 mutations in the RBD of BA.1 relative to the Wuhan Hu-1 strain. Specifically, we measure the ACE2 affinity of all possible combinations of these 15 mutations (2 15  = 32,768 genotypes), spanning all possible evolutionary intermediates from the ancestral Wuhan Hu-1 strain to BA.1. We find that immune escape mutations in BA.1 individually reduce ACE2 affinity but are compensated by epistatic interactions with other affinity-enhancing mutations, including Q498R and N501Y. Thus, the ability of BA.1 to evade immunity while maintaining ACE2 affinity is contingent on acquiring multiple interacting mutations. Our results implicate compensatory epistasis as a key factor driving substantial evolutionary change for SARS-CoV-2 and are consistent with Omicron BA.1 arising from a chronic infection. 
    more » « less
  4. null (Ed.)
    Antibody therapeutics and vaccines are among our last resort to end the raging COVID-19 pandemic. They, however, are prone to over 5000 mutations on the spike (S) protein uncovered by a Mutation Tracker based on over 200 000 genome isolates. It is imperative to understand how mutations will impact vaccines and antibodies in development. In this work, we first study the mechanism, frequency, and ratio of mutations on the S protein which is the common target of most COVID-19 vaccines and antibody therapies. Additionally, we build a library of 56 antibody structures and analyze their 2D and 3D characteristics. Moreover, we predict the mutation-induced binding free energy (BFE) changes for the complexes of S protein and antibodies or ACE2. By integrating genetics, biophysics, deep learning, and algebraic topology, we reveal that most of the 462 mutations on the receptor-binding domain (RBD) will weaken the binding of S protein and antibodies and disrupt the efficacy and reliability of antibody therapies and vaccines. A list of 31 antibody disrupting mutants is identified, while many other disruptive mutations are detailed as well. We also unveil that about 65% of the existing RBD mutations, including those variants recently found in the United Kingdom (UK) and South Africa, will strengthen the binding between the S protein and human angiotensin-converting enzyme 2 (ACE2), resulting in more infectious COVID-19 variants. We discover the disparity between the extreme values of RBD mutation-induced BFE strengthening and weakening of the bindings with antibodies and angiotensin-converting enzyme 2 (ACE2), suggesting that SARS-CoV-2 is at an advanced stage of evolution for human infection, while the human immune system is able to produce optimized antibodies. This discovery, unfortunately, implies the vulnerability of current vaccines and antibody drugs to new mutations. Our predictions were validated by comparison with more than 1400 deep mutations on the S protein RBD. Our results show the urgent need to develop new mutation-resistant vaccines and antibodies and to prepare for seasonal vaccinations. 
    more » « less
  5. Li, Jinyan (Ed.)
    SARS-CoV-2 variants with enhanced transmissibility represent a serious threat to global health. Here we report machine learning models that can predict the impact of receptor-binding domain (RBD) mutations on receptor (ACE2) affinity, which is linked to infectivity, and escape from human serum antibodies, which is linked to viral neutralization. Importantly, the models predict many of the known impacts of RBD mutations in current and former Variants of Concern on receptor affinity and antibody escape as well as novel sets of mutations that strongly modulate both properties. Moreover, these models reveal key opposing impacts of RBD mutations on transmissibility, as many sets of RBD mutations predicted to increase antibody escape are also predicted to reduce receptor affinity and vice versa. These models, when used in concert, capture the complex impacts of SARS-CoV-2 mutations on properties linked to transmissibility and are expected to improve the development of next-generation vaccines and biotherapeutics. 
    more » « less