SARS-CoV-2 evolution in the Omicron era

Temperature-dependent Spike-ACE2 interaction of Omicron subvariants is associated with viral transmission

The continued evolution of severe acute respiratory syndrome 2 (SARS-CoV-2) requires persistent monitoring of its subvariants. Omicron subvariants are responsible for the vast majority of SARS-CoV-2 infections worldwide, with XBB and BA.2.86 sublineages representing more than 90% of circulating strains as of January 2024. To better understand parameters involved in viral transmission, we characterized the functional properties of Spike glycoproteins from BA.2.75, CH.1.1, DV.7.1, BA.4/5, BQ.1.1, XBB, XBB.1, XBB.1.16, XBB.1.5, FD.1.1, EG.5.1, HK.3, BA.2.86 and JN.1. We tested their capacity to evade plasma-mediated recognition and neutralization, binding to angiotensin-converting enzyme 2 (ACE2), their susceptibility to cold inactivation, Spike processing, as well as the impact of temperature on Spike-ACE2 interaction. We found that compared to the early wild-type (D614G) strain, most Omicron subvariants' Spike glycoproteins evolved to escape recognition and neutralization by plasma from individuals who received a fifth dose of bivalent (BA.1 or BA.4/5) mRNA vaccine and improve ACE2 binding, particularly at low temperatures. Moreover, BA.2.86 had the best affinity for ACE2 at all temperatures tested. We found that Omicron subvariants' Spike processing is associated with their susceptibility to cold inactivation. Intriguingly, we found that Spike-ACE2 binding at low temperature was significantly associated with growth rates of Omicron subvariants in humans. Overall, we report that Spikes from newly emerged Omicron subvariants are relatively more stable and resistant to plasma-mediated neutralization, present improved affinity for ACE2 which is associated, particularly at low temperatures, with their growth rates.IMPORTANCEThe persistent evolution of SARS-CoV-2 gave rise to a wide range of variants harboring new mutations in their Spike glycoproteins. Several factors have been associated with viral transmission and fitness such as plasma-neutralization escape and ACE2 interaction. To better understand whether additional factors could be of importance in SARS-CoV-2 variants' transmission, we characterize the functional properties of Spike glycoproteins from several Omicron subvariants. We found that the Spike glycoprotein of Omicron subvariants presents an improved escape from plasma-mediated recognition and neutralization, Spike processing, and ACE2 binding which was further improved at low temperature. Intriguingly, Spike-ACE2 interaction at low temperature is strongly associated with viral growth rate, as such, low temperatures could represent another parameter affecting viral transmission.

Antigenic cartography using variant-specific hamster sera reveals substantial antigenic variation among Omicron subvariants

Severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) has developed substantial antigenic variability. As the majority of the population now has pre-existing immunity due to infection or vaccination, the use of experimentally generated animal immune sera can be valuable for measuring antigenic differences between virus variants. Here, we immunized Syrian hamsters by two successive infections with one of nine SARS-CoV-2 variants. Their sera were titrated against 16 SARS-CoV-2 variants, and the resulting titers were visualized using antigenic cartography. The antigenic map shows a condensed cluster containing all pre-Omicron variants (D614G, Alpha, Delta, Beta, Mu, and an engineered B.1+E484K variant) and considerably more diversity among a selected panel of Omicron subvariants (BA.1, BA.2, BA.4/BA.5, the BA.5 descendants BF.7 and BQ.1.18, the BA.2.75 descendant BN.1.3.1, the BA.2-derived recombinants XBB.2 and EG.5.1, and the BA.2.86 descendant JN.1). Some Omicron subvariants were as antigenically distinct from each other as the wildtype is from the Omicron BA.1 variant. Compared to titers measured in human sera, titers in hamster sera are of higher magnitude, show less fold change, and result in a more compact antigenic map topology. The results highlight the potential of sera from hamsters for the continued antigenic characterization of SARS-CoV-2.

Multi-compartmental diversification of neutralizing antibody lineages dissected in SARS-CoV-2 spike-immunized macaques

The continued evolution of SARS-CoV-2 underscores the need to understand qualitative aspects of the humoral immune response elicited by spike immunization. Here, we combine monoclonal antibody (mAb) isolation with deep B cell receptor (BCR) repertoire sequencing of rhesus macaques immunized with prefusion-stabilized spike glycoprotein. Longitudinal tracing of spike-sorted B cell lineages in multiple immune compartments demonstrates increasing somatic hypermutation and broad dissemination of vaccine-elicited B cells in draining and non-draining lymphoid compartments, including the bone marrow, spleen and, most notably, periaortic lymph nodes. Phylogenetic analysis of spike-specific monoclonal antibody lineages identified through deep repertoire sequencing delineates extensive intra-clonal diversification that shaped neutralizing activity. Structural analysis of the spike in complex with a broadly neutralizing mAb provides a molecular basis for the observed differences in neutralization breadth between clonally related antibodies. Our findings highlight that immunization leads to extensive intra-clonal B cell evolution where members of the same lineage can both retain the original epitope specificity and evolve to recognize additional spike variants not previously encountered.

Combining genomic data and infection estimates to characterize the complex dynamics of SARS-CoV-2 Omicron variants in the US

Omicron surged as a variant of concern in late 2021. Several distinct Omicron variants appeared and overtook each other. We combined variant frequencies and infection estimates from a nowcasting model for each US state to estimate variant-specific infections, attack rates, and effective reproduction numbers (Rt). BA.1 rapidly emerged, and we estimate that it infected 47.7% of the US population before it was replaced by BA.2. We estimate that BA.5 infected 35.7% of the US population, persisting in circulation for nearly 6 months. Other variants-BA.2, BA.4, and XBB-together infected 30.7% of the US population. We found a positive correlation between the state-level BA.1 attack rate and social vulnerability and a negative correlation between the BA.1 and BA.2 attack rates. Our findings illustrate the complex interplay between viral evolution, population susceptibility, and social factors during the Omicron emergence in the US.

A potent, broadly neutralizing human monoclonal antibody that efficiently protects hACE2-transgenic mice from infection with the Wuhan, BA.5, and XBB.1.5 SARS-CoV-2 variants

The COVID-19 pandemic has uncovered the high genetic variability of the SARS-CoV-2 virus and its ability to evade the immune responses that were induced by earlier viral variants. Only a few monoclonal antibodies that have been reported to date are capable of neutralizing a broad spectrum of SARS-CoV-2 variants. Here, we report the isolation of a new broadly neutralizing human monoclonal antibody, iC1. The antibody was identified through sorting the SARS-CoV-1 RBD-stained individual B cells that were isolated from the blood of a vaccinated donor following a breakthrough infection. In vitro, iC1 potently neutralizes pseudoviruses expressing a wide range of SARS-CoV-2 Spike variants, including those of the XBB sublineage. In an hACE2-transgenic mouse model, iC1 provided effective protection against the Wuhan strain of the virus as well as the BA.5 and XBB.1.5 variants. Therefore, iC1 can be considered as a potential component of the broadly neutralizing antibody cocktails resisting the SARS-CoV-2 mutation escape.

Broad sarbecovirus neutralization by combined memory B cell antibodies to ancestral SARS-CoV-2

Antibodies play a pivotal role in protecting from SARS-CoV-2 infection, but their efficacy is challenged by the continuous emergence of viral variants. In this study, we describe two broadly neutralizing antibodies cloned from the memory B cells of a single convalescent individual after infection with ancestral SARS-CoV-2. Cv2.3194, a resilient class 1 anti-RBD antibody, remains active against Omicron sub-variants up to BA.2.86. Cv2.3132, a near pan-Sarbecovirus neutralizer, targets the heptad repeat 2 membrane proximal region. When combined, Cv2.3194 and Cv2.3132 form a complementary SARS-CoV-2 neutralizing antibody cocktail exhibiting a local dose-dependent synergy. Thus, remarkably robust neutralizing memory B cell antibodies elicited in response to ancestral SARS-CoV-2 infection can withstand viral evolution and immune escape. The cooperative effect of such antibody combination may confer a certain level of protection against the latest SARS-CoV-2 variants.

Spike deep mutational scanning helps predict success of SARS-CoV-2 clades

SARS-CoV-2 variants acquire mutations in the spike protein that promote immune evasion1 and affect other properties that contribute to viral fitness, such as ACE2 receptor binding and cell entry2,3. Knowledge of how mutations affect these spike phenotypes can provide insight into the current and potential future evolution of the virus. Here we use pseudovirus deep mutational scanning4 to measure how more than 9,000 mutations across the full XBB.1.5 and BA.2 spikes affect ACE2 binding, cell entry or escape from human sera. We find that mutations outside the receptor-binding domain (RBD) have meaningfully affected ACE2 binding during SARS-CoV-2 evolution. We also measure how mutations to the XBB.1.5 spike affect neutralization by serum from individuals who recently had SARS-CoV-2 infections. The strongest serum escape mutations are in the RBD at sites 357, 420, 440, 456 and 473; however, the antigenic effects of these mutations vary across individuals. We also identify strong escape mutations outside the RBD; however, many of them decrease ACE2 binding, suggesting they act by modulating RBD conformation. Notably, the growth rates of human SARS-CoV-2 clades can be explained in substantial part by the measured effects of mutations on spike phenotypes, suggesting our data could enable better prediction of viral evolution.

Efficient overexpression and purification of severe acute respiratory syndrome coronavirus 2 nucleocapsid proteins in Escherichia coli

The fundamental biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) nucleocapsid protein (Ncap), its use in diagnostic assays and its potential application as a vaccine component have received considerable attention since the outbreak of the Covid19 pandemic in late 2019. Here we report the scalable expression and purification of soluble, immunologically active, SARS-CoV-2 Ncap in Escherichia coli. Codon-optimised synthetic genes encoding the original Ncap sequence and four common variants with an N-terminal 6His affinity tag (sequence MHHHHHHG) were cloned into an inducible expression vector carrying a regulated bacteriophage T5 synthetic promoter controlled by lac operator binding sites. The constructs were used to express Ncap proteins and protocols developed which allow efficient production of purified Ncap with yields of over 200 mg per litre of culture media. These proteins were deployed in ELISA assays to allow comparison of their responses to human sera. Our results suggest that there was no detectable difference between the 6His-tagged and untagged original Ncap proteins but there may be a slight loss of sensitivity of sera to other Ncap isolates.

PHDtools: A platform for pathogen detection and multi-dimensional genetic signatures decoding to realize pathogen genomics data analyses online

OBJECTIVES

Facing the emerging diseases, rapid identification of the pathogen and multi-dimensional characterization of the genomic features at both isolate-level and population-level through high-throughput sequencing data can provide invaluable information to guide the development of antiviral agents and strategies. However, a user-friendly program is in urgent need for clinical laboratories without bioinformatics background to decode the complex big genomics data.

METHODS

In this study, we developed an interactive online platform named PHDtools with a total of 15 functions to analyze metagenomics data to identify the potential pathogen and decode multi-dimensional genetic signatures including intra-/inter-host variations and lineage-level variations. The platform was applied to analyze the meta-genomic data of the samples collected from the 172 imported COVID-19 cases.

RESULTS

According to the analytical results of mNGS, 27 patients were found to have the co-infections of SARS-CoV-2 with either influenza virus (n = 9) or human picobirnavirus (n = 19). Enough coverages of all the assembled SARS-CoV-2 genomes provided the sub-lineages of Omicron variant, and the number of mutations in the non-structural genes and M gene was increased, as well as the intra-host variations occurred in E and M gene were under positive selection (Ka/Ks > 1). These findings of increased or changed mutations in the SARS-CoV-2 genome characterized the current adaptive evolution patterns of Omicron sub-lineages, and revealed the evolution speed of these sub-lineages might increase.

CONCLUSIONS

Consequently, the application of PHDtools has proved that this platform is accurate, user-friendly and convenient for clinical users who are deficient in bioinformatics, and the clinical monitor of SARS-CoV-2 genomes by PHDtools also highlighted the potential evolution features of current SARS-CoV-2 and indicated that the development of anti-SARS-CoV-2 agents and new-designed vaccines should incorporate the gene variations other than S gene.

Characterizing fitness and immune escape of SARS-CoV-2 EG.5 sublineage using elderly serum and nasal organoid

SARS-CoV-2 Omicron variant has evolved into sublineages. Here, we compared the neutralization susceptibility and viral fitness of EG.5.1 and XBB.1.9.1. Serum neutralization antibody titer against EG.5.1 was 1.71-fold lower than that for XBB.1.9.1. However, there was no significant difference in virus replication between EG.5.1 and XBB.1.9.1 in human nasal organoids and TMPRSS2/ACE2 over-expressing A549 cells. No significant difference was observed in competitive fitness and cytokine/chemokine response between EG.5.1 and XBB.1.9.1. Both EG.5.1 and XBB.1.9.1 replicated more robustly in the nasal organoid from a younger adult than that from an older adult. Our findings suggest that enhanced immune escape contributes to the dominance of EG.5.1 over earlier sublineages. The combination of population serum susceptibility testing and viral fitness evaluation with nasal organoids may hold promise in risk assessment of upcoming variants. Utilization of serum specimens and nasal organoid derived from older adults provides a targeted risk assessment for this vulnerable population.

Evidence of SARS-CoV-2 convergent evolution in immunosuppressed patients treated with antiviral therapies

BACKGROUND

The factors contributing to the accelerated convergent evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are not fully understood. Unraveling the contribution of viral replication in immunocompromised patients is important for the early detection of novel mutations and developing approaches to limit COVID-19.

METHODS

We deep sequenced SARS-CoV-2 RNA from 192 patients (64% hospitalized, 39% immunosuppressed) and compared the viral genetic diversity within the patient groups of different immunity and hospitalization status. Serial sampling of 14 patients was evaluated for viral evolution in response to antiviral treatments.

RESULTS

We identified hospitalized and immunosuppressed patients with significantly higher levels of viral genetic diversity and variability. Further evaluation of serial samples revealed accumulated mutations associated with escape from neutralizing antibodies in a subset of the immunosuppressed patients treated with antiviral therapies. Interestingly, the accumulated viral mutations that arose in this early Omicron wave, which were not common in the patient viral lineages, represent convergent mutations that are prevalent in the later Omicron sublineages, including the XBB, BA.2.86.1 and its descendent JN sublineages.

CONCLUSIONS

Our results illustrate the importance of identifying convergent mutations generated during antiviral therapy in immunosuppressed patients, as they may contribute to the future evolutionary landscape of SARS-CoV-2. Our study also provides evidence of a correlation between SARS-CoV-2 convergent mutations and specific antiviral treatments. Evaluating high-confidence genomes from distinct waves in the pandemic with detailed patient metadata allows for discerning of convergent mutations that contribute to the ongoing evolution of SARS-CoV-2.

SARS-CoV-2 strains bearing Omicron BA.1 spike replicate in C57BL/6 mice

Introduction

SARS-CoV-2, the cause of the COVID pandemic, is an RNA virus with a high propensity to mutate. Successive virus variants, including variants of concern (VOC), have emerged with increased transmission or immune escape. The original pandemic virus and early variants replicated poorly, if at all, in mice at least partly due to a mismatch between the receptor binding domain on the viral spike protein and the murine angiotensin converting enzyme 2 (ACE2). Omicron VOC emerged in late 2021 harboring > 50 new mutations, 35 of them in the spike protein. This variant resulted in a very large wave of infections, even in the face of prior immunity, albeit being inherently less severe than earlier variants. Reflecting the lower severity reported in humans, Omicron displayed attenuated infection in hamsters and also in the K18-hACE2 mouse model. K18-hACE2 mice express both the human ACE2 as well as the endogenous mouse ACE2.

Methods

Here we infected hACE2 knock-in mice that express only human ACE2 and no murine ACE2, or C57BL/6 wildtype mice with SARS-CoV-2 D614G (first-wave isolate), Delta or Omicron BA.1 variants and assessed infectivity and downstream innate immune responses.

Results

While replication of SARS-CoV-2 Omicron was lower in the lungs of hACE2 knock-in mice compared with SARS-CoV-2 D614G and VOC Delta, it replicated more efficiently than the earlier variants in C57BL/6 wildtype mice. This opens the opportunity to test the effect of host genetics on SARS-CoV-2 infections in wildtype mice. As a proof of principle, we tested Omicron infection in mice lacking expression of the interferon-alpha receptor-1 (IFNAR1). In these mice we found that loss of type I IFN receptor signaling resulted in higher viral loads in the lungs were detected. Finally, using a chimeric virus of first wave SARS-CoV-2 harboring the Omicron spike protein, we show that Omicron spike increase infection of C57BL/6 wildtype mice, but non-spike genes of Omicron confer attenuation of viral replication.

Discussion

Since this chimeric virus efficiently infected C57BL/6 wildtype mice, and replicated in their lungs, our findings illustrate a pathway for genetic mapping of virushost interactions during SARS-CoV-2 infection.

SARS-CoV-2 Omicron: Viral Evolution, Immune Evasion, and Alternative Durable Therapeutic Strategies

Since the SARS-CoV-2 Omicron virus has gained dominance worldwide, its continual evolution with unpredictable mutations and patterns has revoked all authorized immunotherapeutics. Rapid viral evolution has also necessitated several rounds of vaccine updates in order to provide adequate immune protection. It remains imperative to understand how Omicron evolves into different subvariants and causes immune escape as this could help reevaluate the current intervention strategies mostly implemented in the clinics as emergency measures to counter the pandemic and, importantly, develop new solutions. Here, we provide a review focusing on the major events of Omicron viral evolution, including the features of spike mutation that lead to immune evasion against monoclonal antibody (mAb) therapy and vaccination, and suggest alternative durable options such as the ACE2-based experimental therapies superior to mAbs to address this unprecedented evolution of Omicron virus. In addition, this type of unique ACE2-based virus-trapping molecules can counter all zoonotic SARS coronaviruses, either from unknown animal hosts or from established wild-life reservoirs of SARS-CoV-2, and even seasonal alpha coronavirus NL63 that depends on human ACE2 for infection.

Severity of Omicron Subvariants and Vaccine Impact in Catalonia, Spain

In the current COVID-19 landscape dominated by Omicron subvariants, understanding the timing and efficacy of vaccination against emergent lineages is crucial for planning future vaccination campaigns, yet detailed studies stratified by subvariant, vaccination timing, and age groups are scarce. This retrospective study analyzed COVID-19 cases from December 2021 to January 2023 in Catalonia, Spain, focusing on vulnerable populations affected by variants BA.1, BA.2, BA.5, and BQ.1 and including two national booster campaigns. Our database includes detailed information such as dates of diagnosis, hospitalization and death, last vaccination, and cause of death, among others. We evaluated the impact of vaccination on disease severity by age, variant, and vaccination status, finding that recent vaccination significantly mitigated severity across all Omicron subvariants, although efficacy waned six months post-vaccination, except for BQ.1, which showed more stable levels. Unvaccinated individuals had higher hospitalization and mortality rates. Our results highlight the importance of periodic vaccination to reduce severe outcomes, which are influenced by variant and vaccination timing. Although the seasonality of COVID-19 is uncertain, our analysis suggests the potential benefit of annual vaccination in populations >60 years old, probably in early fall, if COVID-19 eventually exhibits a major peak similar to other respiratory viruses.

Surfing the Waves of SARS-CoV-2: Analysis of Viral Genome Variants Using an NGS Survey in Verona, Italy

The availability of new technologies for deep sequencing, including next-generation sequencing (NGS), allows for the detection of viral genome variations. The epidemiological determination of SARS-CoV-2 viral genome changes during the pandemic waves displayed the genome evolution and subsequent onset of variants over time. These variants were often associated with a different impact on viral transmission and disease severity. We investigated, in a retrospective study, the trend of SARS-CoV-2-positive samples collected from the start of the Italian pandemic (January 2020) to June 2023. In addition, viral RNAs extracted from 938 nasopharyngeal swab samples were analyzed using NGS between February 2022 and June 2023. Sequences were analyzed with bioinformatic tools to identify lineages and mutations and for phylogenetic studies. Six pandemic waves were detected. In our samples, we predominantly detected BA.2, BQ.1, BA.5.1, BA.5.2, and, more recently, XBB.1 and its subvariants. The data describe the SARS-CoV-2 genome evolution involved in viral interactions with the host and the dynamics of specific genome mutations and deletions.

Blood Count and Renal Functionality Assessments in the Emergency Section Disclose Morbidity and Mortality in Omicron COVID-19 Patients: A Retrospective Study

Background: SARS-CoV-2 is the coronavirus responsible for the COVID-19 pandemic. Even though we are no longer in a pandemic situation, people are still getting infected, some of them need hospitalization and a few of them die. Methods: We conducted a retrospective study including 445 patients who accessed the Emergency Section of Policlinico Umberto I, Rome, Italy, where they had routine blood exams. In this study, we focused on the complete blood count, serum creatinine and azotemia. The data were analyzed using ANOVA, Spearman correlation and ROC analyses. They were divided into four groups based on their clinical outcomes: (1) the emergency group (patients who had mild forms and were quickly discharged); (2) the hospital ward group (patients who were admitted to the emergency section and were then hospitalized in a COVID-19 ward); (3) the intensive care unit (ICU) group (patients who required intensive assistance after the admission in the emergency section); (4) the deceased group (patients who had a fatal outcome after admission to the emergency section). Results: We found significant changes for creatinine, azotemia, hematocrit, mean corpuscular hemoglobin concentration, basophils, monocytes, red blood cell distribution width, hemoglobin, hematocrit and red blood cell numbers using ANOVA according to their clinical outcomes, particularly for the deceased group. Also, we found linear correlations of clinical outcomes with eosinophils, hemoglobin, hematocrit, mean corpuscular hemoglobin concentration, lymphocyte, neutrophil, platelet and red blood cell number and red blood cell distribution width. Conclusions: This study discloses an early association between "classical" routine blood biomarkers and the severity of clinical outcomes in Omicron patients.

Emerging variants develop total escape from potent monoclonal antibodies induced by BA.4/5 infection

The rapid evolution of SARS-CoV-2 is driven in part by a need to evade the antibody response in the face of high levels of immunity. Here, we isolate spike (S) binding monoclonal antibodies (mAbs) from vaccinees who suffered vaccine break-through infections with Omicron sub lineages BA.4 or BA.5. Twenty eight potent antibodies are isolated and characterised functionally, and in some cases structurally. Since the emergence of BA.4/5, SARS-CoV-2 has continued to accrue mutations in the S protein, to understand this we characterize neutralization of a large panel of variants and demonstrate a steady attrition of neutralization by the panel of BA.4/5 mAbs culminating in total loss of function with recent XBB.1.5.70 variants containing the so-called 'FLip' mutations at positions 455 and 456. Interestingly, activity of some mAbs is regained on the recently reported variant BA.2.86.

Ensemble-Based Mutational Profiling and Network Analysis of the SARS-CoV-2 Spike Omicron XBB Lineages for Interactions with the ACE2 Receptor and Antibodies: Cooperation of Binding Hotspots in Mediating Epistatic Couplings Underlies Binding Mechanism and Immune Escape

In this study, we performed a computational study of binding mechanisms for the SARS-CoV-2 spike Omicron XBB lineages with the host cell receptor ACE2 and a panel of diverse class one antibodies. The central objective of this investigation was to examine the molecular factors underlying epistatic couplings among convergent evolution hotspots that enable optimal balancing of ACE2 binding and antibody evasion for Omicron variants BA.1, BA2, BA.3, BA.4/BA.5, BQ.1.1, XBB.1, XBB.1.5, and XBB.1.5 + L455F/F456L. By combining evolutionary analysis, molecular dynamics simulations, and ensemble-based mutational scanning of spike protein residues in complexes with ACE2, we identified structural stability and binding affinity hotspots that are consistent with the results of biochemical studies. In agreement with the results of deep mutational scanning experiments, our quantitative analysis correctly reproduced strong and variant-specific epistatic effects in the XBB.1.5 and BA.2 variants. It was shown that Y453W and F456L mutations can enhance ACE2 binding when coupled with Q493 in XBB.1.5, while these mutations become destabilized when coupled with the R493 position in the BA.2 variant. The results provided a molecular rationale of the epistatic mechanism in Omicron variants, showing a central role of the Q493/R493 hotspot in modulating epistatic couplings between convergent mutational sites L455F and F456L in XBB lineages. The results of mutational scanning and binding analysis of the Omicron XBB spike variants with ACE2 receptors and a panel of class one antibodies provide a quantitative rationale for the experimental evidence that epistatic interactions of the physically proximal binding hotspots Y501, R498, Q493, L455F, and F456L can determine strong ACE2 binding, while convergent mutational sites F456L and F486P are instrumental in mediating broad antibody resistance. The study supports a mechanism in which the impact on ACE2 binding affinity is mediated through a small group of universal binding hotspots, while the effect of immune evasion could be more variant-dependent and modulated by convergent mutational sites in the conformationally adaptable spike regions.

AlphaFold2-Enabled Atomistic Modeling of Structure, Conformational Ensembles, and Binding Energetics of the SARS-CoV-2 Omicron BA.2.86 Spike Protein with ACE2 Host Receptor and Antibodies: Compensatory Functional Effects of Binding Hotspots in Modulating Mechanisms of Receptor Binding and Immune Escape

The latest wave of SARS-CoV-2 Omicron variants displayed a growth advantage and increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with atomistic simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both the structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that the AlphaFold2-predicted structural ensemble of the BA.2.86 spike protein complex with ACE2 can accurately capture the main conformational states of the Omicron variant. Complementary to AlphaFold2 structural predictions, microsecond molecular dynamics simulations reveal the details of the conformational landscape and produced equilibrium ensembles of the BA.2.86 structures that are used to perform mutational scanning of spike residues and characterize structural stability and binding energy hotspots. The ensemble-based mutational profiling of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 revealed a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 convergent mutational hotspots R403K, F486P, and R493Q. To examine the immune evasion properties of BA.2.86 in atomistic detail, we performed structure-based mutational profiling of the spike protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against the BA.2.86 variant. The results revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have evolved to outcompete other Omicron subvariants by improving immune evasion while preserving binding affinity with ACE2 via through a compensatory effect of R493Q and F486P convergent mutational hotspots. This study demonstrated that an integrative approach combining AlphaFold2 predictions with complementary atomistic molecular dynamics simulations and robust ensemble-based mutational profiling of spike residues can enable accurate and comprehensive characterization of structure, dynamics, and binding mechanisms of newly emerging Omicron variants.

Divergent pathogenetic outcomes in BALB/c mice following Omicron subvariant infection

Following the emergence of B.1.1.529 Omicron, the SARS-CoV-2 virus evolved into a significant number of sublineage variants that possessed numerous mutations throughout the genome, but particularly within the spike glycoprotein (S) gene. For example, the BQ.1.1 and the XBB.1 and XBB.1.5 subvariants contained 34 and 41 mutations in S, respectively. However, these variants elicited largely replication only or mild disease phenotypes in mice. To better model pathogenic outcomes and measure countermeasure performance, we developed mouse adapted versions (BQ.1.1 MA; XBB.1 MA; XBB.1.5 MA) that reflect more pathogenic acute phase pulmonary disease symptoms of SARS-CoV-2, as well as derivative strains expressing nano-luciferase (nLuc) in place of ORF7 (BQ.1.1 nLuc; XBB.1 nLuc; XBB.1.5 nLuc). Amongst the mouse adapted (MA) viruses, a wide range of disease outcomes were observed including mortality, weight loss, lung dysfunction, and tissue viral loads in the lung and nasal turbinates. Intriguingly, XBB.1 MA and XBB.1.5 MA strains, which contained identical mutations throughout except at position F486S/P in S, exhibited divergent disease outcomes in mice (Ao et al., 2023). XBB.1.5 MA infection was associated with significant weight loss and ∼45 % mortality across two independent studies, while XBB.1 MA infected animals suffered from mild weight loss and only 10 % mortality across the same two independent studies. Additionally, the development and use of nanoluciferase expressing strains provided moderate throughput for live virus neutralization assays. The availability of small animal models for the assessment of Omicron VOC disease potential will enable refined capacity to evaluate the efficacy of on market and pre-clinical therapeutics and interventions.

Post-pandemic memory T cell response to SARS-CoV-2 is durable, broadly targeted, and cross-reactive to the hypermutated BA.2.86 variant

Ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) evolution has given rise to recombinant Omicron lineages that dominate globally (XBB.1), as well as the emergence of hypermutated variants (BA.2.86). In this context, durable and cross-reactive T cell immune memory is critical for continued protection against severe COVID-19. We examined T cell responses to SARS-CoV-2 approximately 1.5 years since Omicron first emerged. We describe sustained CD4+ and CD8+ spike-specific T cell memory responses in healthcare workers in South Africa (n = 39) who were vaccinated and experienced at least one SARS-CoV-2 infection. Spike-specific T cells are highly cross-reactive with all Omicron variants tested, including BA.2.86. Abundant nucleocapsid and membrane-specific T cells are detectable in most participants. The bulk of SARS-CoV-2-specific T cell responses have an early-differentiated phenotype, explaining their persistent nature. Overall, hybrid immunity leads to the accumulation of spike and non-spike T cells evident 3.5 years after the start of the pandemic, with preserved recognition of highly mutated SARS-CoV-2 variants.

RBD mutations at the residues K417, E484, N501 reduced immunoreactivity with antisera from vaccinated and COVID-19 recovered patients

Introduction

It is unclear whether induced spike protein-specific antibodies due to infections with SARS-CoV-2 or to the prototypic Wuhan isolate-based vaccination can immune-react with the emerging variants of SARS-CoV-2.

Aim/objectives

The main objective of the study was to measure the immunoreactivity of induced antibodies postvaccination with Covishield™ (ChAdOx1 nCoV-19 coronavirus vaccines) or infections with SARS-CoV-2 by using selected peptides of the spike protein of wild type and variants of SARS-CoV-2.

Methodology

Thirty patients who had recovered from SARS-CoV-2 infections and 30 individuals vaccinated with both doses of Covishield™ were recruited for the study. Venous blood samples (5 mL) were collected at a single time point from patients within 3-4 weeks of recovery from SARS-CoV-2 infections or receiving both doses of Covishield™ vaccines. The serum levels of total immunoglobulin were measured in both study groups. A total of 12 peptides of 10 to 24 amino acids length spanning to the receptor-binding domain (RBD) of wild type of SARS-CoV-2 and their variants were synthesized. The serum levels of immune-reactive antibodies were measured using these peptides.

Results

The serum levels of total antibodies were found to be significantly (p<0.001) higher in the vaccinated individuals as compared to COVID-19 recovered patients. Our study reported that the mutations in the RBD at the residues K417, E484, and N501 have been associated with reduced immunoreactivity with anti-sera of vaccinated people and COVID-19 recovered patients.

Conclusion

The amino acid substitutions at the RBD of SARS-CoV-2 have been associated with a higher potential to escape the humoral immune response.

Humoral immune escape by current SARS-CoV-2 variants BA.2.86 and JN.1, December 2023

Variant BA.2.86 and its descendant, JN.1, of SARS-CoV-2 are rising in incidence across Europe and globally. We isolated recent JN.1, BA.2.86, EG.5, XBB.1.5 and earlier variants. We tested live virus neutralisation of sera taken in September 2023 from vaccinated and exposed healthy persons (n = 39). We found clear neutralisation escape against recent variants but no specific pronounced escape for BA.2.86 or JN.1. Neutralisation escape corresponds to recent variant predominance but may not be causative of the recent upsurge in JN.1 incidence.

AlphaFold2-Enabled Atomistic Modeling of Epistatic Binding Mechanisms for the SARS-CoV-2 Spike Omicron XBB.1.5, EG.5 and FLip Variants: Convergent Evolution Hotspots Cooperate to Control Stability and Conformational Adaptability in Balancing ACE2 Binding and Antibody Resistance

In this study, we combined AI-based atomistic structural modeling and microsecond molecular simulations of the SARS-CoV-2 Spike complexes with the host receptor ACE2 for XBB.1.5+L455F, XBB.1.5+F456L(EG.5) and XBB.1.5+L455F/F456L (FLip) lineages to examine the mechanisms underlying the role of convergent evolution hotspots in balancing ACE2 binding and antibody evasion. Using the ensemble-based mutational scanning of the spike protein residues and physics-based rigorous computations of binding affinities, we identified binding energy hotspots and characterized molecular basis underlying epistatic couplings between convergent mutational hotspots. Consistent with the experiments, the results revealed the mediating role of Q493 hotspot in synchronization of epistatic couplings between L455F and F456L mutations providing a quantitative insight into the mechanism underlying differences between XBB lineages. Mutational profiling is combined with network-based model of epistatic couplings showing that the Q493, L455 and F456 sites mediate stable communities at the binding interface with ACE2 and can serve as stable mediators of non-additive couplings. Structure-based mutational analysis of Spike protein binding with the class 1 antibodies quantified the critical role of F456L and F486P mutations in eliciting strong immune evasion response. The results of this analysis support a mechanism in which the emergence of EG.5 and FLip variants may have been dictated by leveraging strong epistatic effects between several convergent revolutionary hotspots that provide synergy between the improved ACE2 binding and broad neutralization resistance. This interpretation is consistent with the notion that functionally balanced substitutions which simultaneously optimize immune evasion and high ACE2 affinity may continue to emerge through lineages with beneficial pair or triplet combinations of RBD mutations involving mediators of epistatic couplings and sites in highly adaptable RBD regions.

Accurate Characterization of Conformational Ensembles and Binding Mechanisms of the SARS-CoV-2 Omicron BA.2 and BA.2.86 Spike Protein with the Host Receptor and Distinct Classes of Antibodies Using AlphaFold2-Augmented Integrative Computational Modeling

The latest wave SARS-CoV-2 Omicron variants displayed a growth advantage and the increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with all-atom MD simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that AlphaFold2-predicted conformational ensemble of the BA.2.86 spike protein complex can accurately capture the main dynamics signatures obtained from microscond molecular dynamics simulations. The ensemble-based dynamic mutational scanning of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 dissected the role of the BA.2 and BA.2.86 backgrounds in modulating binding free energy changes revealing a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 mutational sites R403K, F486P and R493Q. To examine immune evasion properties of BA.2.86 in atomistic detail, we performed large scale structure-based mutational profiling of the S protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against BA.2.86 variant. The results quantified specific function of the BA.2.86 mutations to ensure broad resistance against different classes of RBD antibodies. This study revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have primarily evolved to improve immune escape while modulating binding affinity with ACE2 through cooperative effect of R403K, F486P and R493Q mutations. The study supports a hypothesis that the impact of the increased ACE2 binding affinity on viral fitness is more universal and is mediated through cross-talk between convergent mutational hotspots, while the effect of immune evasion could be more variant-dependent.

Full-spike deep mutational scanning helps predict the evolutionary success of SARS-CoV-2 clades

SARS-CoV-2 variants acquire mutations in spike that promote immune evasion and impact other properties that contribute to viral fitness such as ACE2 receptor binding and cell entry. Knowledge of how mutations affect these spike phenotypes can provide insight into the current and potential future evolution of the virus. Here we use pseudovirus deep mutational scanning to measure how >9,000 mutations across the full XBB.1.5 and BA.2 spikes affect ACE2 binding, cell entry, or escape from human sera. We find that mutations outside the receptor-binding domain (RBD) have meaningfully impacted ACE2 binding during SARS-CoV-2 evolution. We also measure how mutations to the XBB.1.5 spike affect neutralization by serum from individuals who recently had SARS-CoV-2 infections. The strongest serum escape mutations are in the RBD at sites 357, 420, 440, 456, and 473-however, the antigenic impacts of these mutations vary across individuals. We also identify strong escape mutations outside the RBD; however many of them decrease ACE2 binding, suggesting they act by modulating RBD conformation. Notably, the growth rates of human SARS-CoV-2 clades can be explained in substantial part by the measured effects of mutations on spike phenotypes, suggesting our data could enable better prediction of viral evolution.