Convergent evolution of SARS-CoV-2 XBB lineages on receptor-binding domain 455-456 synergistically enhances antibody evasion and ACE2 binding

Antigenicity assessment of SARS-CoV-2 saltation variant BA.2.87.1

The recent emergence of a SARS-CoV-2 saltation variant, BA.2.87.1, which features 65 spike mutations relative to BA.2, has attracted worldwide attention. In this study, we elucidate the antigenic characteristics and immune evasion capability of BA.2.87.1. Our findings reveal that BA.2.87.1 is more susceptible to XBB-induced humoral immunity compared to JN.1. Notably, BA.2.87.1 lacks critical escaping mutations in the receptor binding domain (RBD) thus allowing various classes of neutralizing antibodies (NAbs) that were escaped by XBB or BA.2.86 subvariants to neutralize BA.2.87.1, although the deletions in the N-terminal domain (NTD), specifically 15-23del and 136-146del, compensate for the resistance to humoral immunity. Interestingly, several neutralizing antibody drugs have been found to restore their efficacy against BA.2.87.1, including SA58, REGN-10933 and COV2-2196. Hence, our results suggest that BA.2.87.1 may not become widespread until it acquires multiple RBD mutations to achieve sufficient immune evasion comparable to that of JN.1.

Structural basis for the evolution and antibody evasion of SARS-CoV-2 BA.2.86 and JN.1 subvariants

The Omicron subvariants of SARS-CoV-2, especially for BA.2.86 and JN.1, have rapidly spread across multiple countries, posing a significant threat in the ongoing COVID-19 pandemic. Distinguished by 34 additional mutations on the Spike (S) protein compared to its BA.2 predecessor, the implications of BA.2.86 and its evolved descendant, JN.1 with additional L455S mutation in receptor-binding domains (RBDs), are of paramount concern. In this work, we systematically examine the neutralization susceptibilities of SARS-CoV-2 Omicron subvariants and reveal the enhanced antibody evasion of BA.2.86 and JN.1. We also determine the cryo-EM structures of the trimeric S proteins from BA.2.86 and JN.1 in complex with the host receptor ACE2, respectively. The mutations within the RBDs of BA.2.86 and JN.1 induce a remodeling of the interaction network between the RBD and ACE2. The L455S mutation of JN.1 further induces a notable shift of the RBD-ACE2 interface, suggesting the notably reduced binding affinity of JN.1 than BA.2.86. An analysis of the broadly neutralizing antibodies possessing core neutralizing epitopes reveals the antibody evasion mechanism underlying the evolution of Omicron BA.2.86 subvariant. In general, we construct a landscape of evolution in virus-receptor of the circulating Omicron subvariants.

Neutralization escape, infectivity, and membrane fusion of JN.1-derived SARS-CoV-2 SLip, FLiRT, and KP.2 variants

We investigate JN.1-derived subvariants SLip, FLiRT, and KP.2 for neutralization by antibodies in vaccinated individuals, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected patients, or class III monoclonal antibody S309. Compared to JN.1, SLip, KP.2, and especially FLiRT exhibit increased resistance to bivalent-vaccinated and BA.2.86/JN.1-wave convalescent human sera. XBB.1.5 monovalent-vaccinated hamster sera robustly neutralize FLiRT and KP.2 but have reduced efficiency for SLip. All subvariants are resistant to S309 and show decreased infectivity, cell-cell fusion, and spike processing relative to JN.1. Modeling reveals that L455S and F456L in SLip reduce spike binding for ACE2, while R346T in FLiRT and KP.2 strengthens it. These three mutations, alongside D339H, alter key epitopes in spike, likely explaining the reduced sensitivity of these subvariants to neutralization. Our findings highlight the increased neutralization resistance of JN.1 subvariants and suggest that future vaccine formulations should consider the JN.1 spike as an immunogen, although the current XBB.1.5 monovalent vaccine could still offer adequate protection.

Molecular epidemiology and population immunity of SARS-CoV-2 in Guangdong (2022-2023) following a pivotal shift in the pandemic

The SARS-CoV-2 Omicron variant sparked the largest wave of infections worldwide. Mainland China eased its strict COVID-19 measures in late 2022 and experienced two nationwide Omicron waves in 2023. Here, we investigated lineage distribution and virus evolution in Guangdong, China, 2022-2023 by comparing 5813 local viral genomes with the datasets from other regions of China and worldwide. Additionally, we conducted three large-scale serological surveys involving 1696 participants to measure their immune response to the BA.5 and XBB.1.9 before and after the corresponding waves. Our findings revealed the Omicron variants, mainly the BA.5.2.48 lineage, causing infections in over 90% of individuals across different age groups within a month. This rapid spread led to the establishment of widespread immunity, limiting the virus's ability to further adaptive mutation and dissemination. While similar immune responses to BA.5 were observed across all age groups after the initial wave, children aged 3 to 11 developed a stronger cross immune response to the XBB.1.9 strain, possibly explaining their lower infection rates in the following XBB.1 wave. Reinfection with Omicron XBB.1 variant triggered a more potent neutralizing immune response among older adults. These findings highlight the impact of age-specific immune responses on viral spread in potential future waves.

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.

SARS-CoV-2 JN.1 variant evasion of IGHV3-53/3-66 B cell germlines

The severe acute respiratory syndrome coronavirus 2 variant JN.1 recently emerged as the dominant variant despite having only one amino acid change on the spike (S) protein receptor binding domain (RBD) compared with the ancestral BA.2.86, which never represented more than 5% of global variants. To define at the molecular level the JN.1 ability to spread globally, we interrogated a panel of 899 neutralizing human monoclonal antibodies. Our data show that the single leucine-455-to-serine mutation in the JN.1 spike protein RBD unleashed the global spread of JN.1, likely occurring by elimination of more than 70% of the neutralizing antibodies mediated by IGHV3-53/3-66 germlines. However, the resilience of class 3 antibodies with low neutralization potency but strong Fc functions may explain the absence of JN.1 severe disease.

Subsequent Waves of Convergent Evolution in SARS-CoV-2 Genes and Proteins

Beginning in 2022, following widespread infection and vaccination among the global population, the SARS-CoV-2 virus mainly evolved to evade immunity derived from vaccines and past infections. This review covers the convergent evolution of structural, nonstructural, and accessory proteins in SARS-CoV-2, with a specific look at common mutations found in long-lasting infections that hint at the virus potentially reverting to an enteric sarbecovirus type.

An update on the anti-spike monoclonal antibody pipeline for SARS-CoV-2

BACKGROUND

Anti-spike monoclonal antibodies represent one of the most tolerable prophylaxis and therapies for COVID-19 in frail and immunocompromised patients. Unfortunately, viral evolution in Omicron has led all of them to failure.

OBJECTIVES

We review here the current pipeline of anti-spike mAb's, discussing in detail the most promising candidates.

SOURCES

We scanned PubMed, ClinicalTrials.gov and manufacturers' press releases for clinical studies on anti-spike monoclonal antibodies.

CONTENT

We present state-of-art data clinical progress for AstraZeneca's AZD3152, Invivyd's VYD222, Regeneron's REGN-17092 and Aerium Therapeutics' AER-800.

IMPLICATIONS

The anti-spike monoclonal antibody clinical pipeline is currently limited to few agents (most being single antibodies) with unknown efficacy against the dominant JN.1 sublineage. The field of antibody-based therapies requires boosting by both manufacturers and institutions.

SARS-CoV-2 Omicron XBB lineage spike structures, conformations, antigenicity, and receptor recognition

A recombinant lineage of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant, named XBB, appeared in late 2022 and evolved descendants that successively swept local and global populations. XBB lineage members were noted for their improved immune evasion and transmissibility. Here, we determine cryoelectron microscopy (cryo-EM) structures of XBB.1.5, XBB.1.16, EG.5, and EG.5.1 spike (S) ectodomains to reveal reinforced 3-receptor binding domain (RBD)-down receptor-inaccessible closed states mediated by interprotomer RBD interactions previously observed in BA.1 and BA.2. Improved XBB.1.5 and XBB.1.16 RBD stability compensated for stability loss caused by early Omicron mutations, while the F456L substitution reduced EG.5 RBD stability. S1 subunit mutations had long-range impacts on conformation and epitope presentation in the S2 subunit. Our results reveal continued S protein evolution via simultaneous optimization of multiple parameters, including stability, receptor binding, and immune evasion, and the dramatic effects of relatively few residue substitutions in altering the S protein conformational landscape.

Deep mutational scanning of SARS-CoV-2 Omicron BA.2.86 and epistatic emergence of the KP.3 variant

Deep mutational scanning experiments aid in the surveillance and forecasting of viral evolution by providing prospective measurements of mutational effects on viral traits, but epistatic shifts in the impacts of mutations can hinder viral forecasting when measurements were made in outdated strain backgrounds. Here, we report measurements of the impact of all single amino acid mutations on ACE2-binding affinity and protein folding and expression in the SARS-CoV-2 Omicron BA.2.86 spike receptor-binding domain (RBD). As with other SARS-CoV-2 variants, we find a plastic and evolvable basis for receptor binding, with many mutations at the ACE2 interface maintaining or even improving ACE2-binding affinity. Despite its large genetic divergence, mutational effects in BA.2.86 have not diverged greatly from those measured in its Omicron BA.2 ancestor. However, we do identify strong positive epistasis among subsequent mutations that have accrued in BA.2.86 descendants. Specifically, the Q493E mutation that decreased ACE2-binding affinity in all previous SARS-CoV-2 backgrounds is reversed in sign to enhance human ACE2-binding affinity when coupled with L455S and F456L in the currently emerging KP.3 variant. Our results point to a modest degree of epistatic drift in mutational effects during recent SARS-CoV-2 evolution but highlight how these small epistatic shifts can have important consequences for the emergence of new SARS-CoV-2 variants.

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.

Atomistic Prediction of Structures, Conformational Ensembles and Binding Energetics for the SARS-CoV-2 Spike JN.1, KP.2 and KP.3 Variants Using AlphaFold2 and Molecular Dynamics Simulations: Mutational Profiling and Binding Free Energy Analysis Reveal Epistatic Hotspots of the ACE2 Affinity and Immune Escape

The most recent wave of SARS-CoV-2 Omicron variants descending from BA.2 and BA.2.86 exhibited improved viral growth and fitness due to convergent evolution of functional hotspots. These hotspots operate in tandem to optimize both receptor binding for effective infection and immune evasion efficiency, thereby maintaining overall viral fitness. The lack of molecular details on structure, dynamics and binding energetics of the latest FLiRT and FLuQE variants with the ACE2 receptor and antibodies provides a considerable challenge that is explored in this study. We combined AlphaFold2-based atomistic predictions of structures and conformational ensembles of the SARS-CoV-2 Spike complexes with the host receptor ACE2 for the most dominant Omicron variants JN.1, KP.1, KP.2 and KP.3 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 computations of binding affinities, we identified binding energy hotspots and characterized molecular basis underlying epistatic couplings between convergent mutational hotspots. The results suggested that the existence of epistatic interactions between convergent mutational sites at L455, F456, Q493 positions that enable to protect and restore ACE2 binding affinity while conferring beneficial immune escape. To examine immune escape mechanisms, we performed structure-based mutational profiling of the spike protein binding with several classes of antibodies that displayed impaired neutralization against BA.2.86, JN.1, KP.2 and KP.3. The results confirmed the experimental data that JN.1, KP.2 and KP.3 harboring the L455S and F456L mutations can significantly impair the neutralizing activity of class-1 monoclonal antibodies, while the epistatic effects mediated by F456L can facilitate the subsequent convergence of Q493E changes to rescue ACE2 binding. Structural and energetic analysis provided a rationale to the experimental results showing that BD55-5840 and BD55-5514 antibodies that bind to different binding epitopes can retain neutralizing efficacy against all examined variants BA.2.86, JN.1, KP.2 and KP.3. The results support the notion that evolution of Omicron variants may favor emergence of lineages with beneficial combinations of mutations involving mediators of epistatic couplings that control balance of high ACE2 affinity and immune evasion.

Spike N354 glycosylation augments SARS-CoV-2 fitness for human adaptation through structural plasticity

Selective pressures have given rise to a number of SARS-CoV-2 variants during the prolonged course of the COVID-19 pandemic. Recently evolved variants differ from ancestors in additional glycosylation within the spike protein receptor-binding domain (RBD). Details of how the acquisition of glycosylation impacts viral fitness and human adaptation are not clearly understood. Here, we dissected the role of N354-linked glycosylation, acquired by BA.2.86 sub-lineages, as a RBD conformational control element in attenuating viral infectivity. The reduced infectivity is recovered in the presence of heparin sulfate, which targets the 'N354 pocket' to ease restrictions of conformational transition resulting in a 'RBD-up' state, thereby conferring an adjustable infectivity. Furthermore, N354 glycosylation improved spike cleavage and cell-cell fusion, and in particular escaped one subset of ADCC antibodies. Together with reduced immunogenicity in hybrid immunity background, these indicate a single spike amino acid glycosylation event provides selective advantage in humans through multiple mechanisms.

Structural Basis for Evasion of New SARS-CoV-2 Variants from the Potent Virus-Neutralizing Nanobody Targeting the S-Protein Receptor-Binding Domain

COVID-19 has caused millions of deaths and many times more infections worldwide, emphasizing the unpreparedness of the global health system in the face of new infections and the key role for vaccines and therapeutics, including virus-neutralizing antibodies, in prevention and containment of the disease. Continuous evolution of the SARS-CoV-2 coronavirus has been causing its new variants to evade the action of the immune system, which highlighted the importance of detailed knowledge of the epitopes of already selected potent virus-neutralizing antibodies. A single-chain antibody ("nanobody") targeting the SARS-CoV-2 receptor-binding domain (RBD), clone P2C5, had exhibited robust virus-neutralizing activity against all SARS-CoV-2 variants and, being a major component of the anti-COVID-19 formulation "GamCoviMab", had successfully passed Phase I of clinical trials. However, after the emergence of the Delta and XBB variants, a decrease in the neutralizing activity of this nanobody was observed. Here we report on the successful crystal structure determination of the RBD:P2C5 complex at 3.1 Å, which revealed the intricate protein-protein interface, sterically occluding full ACE2 receptor binding by the P2C5-neutralized RBD. Moreover, the structure revealed the developed RBD:P2C5 interface centered around residues Leu452 and Phe490, thereby explaining the evasion of the Delta or Omicron XBB, but not Omicron B.1.1.529 variant, as a result of the single L452R or F490S mutations, respectively, from the action of P2C5. The structure obtained is expected to foster nanobody engineering in order to rescue neutralization activity and will facilitate epitope mapping for other neutralizing nanobodies by competition assays.

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.

Exploring conformational landscapes and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variant complexes with the ACE2 receptor using AlphaFold2-based structural ensembles and molecular dynamics simulations

In this study, we combined AlphaFold-based approaches for atomistic modeling of multiple protein states and microsecond molecular simulations to accurately characterize conformational ensembles evolution and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variants BA.1, BA.2, BA.2.75, BA.3, BA.4/BA.5 and BQ.1.1. We employed and validated several different adaptations of the AlphaFold methodology for modeling of conformational ensembles including the introduced randomized full sequence scanning for manipulation of sequence variations to systematically explore conformational dynamics of Omicron spike protein complexes with the ACE2 receptor. Microsecond atomistic molecular dynamics (MD) simulations provide a detailed characterization of the conformational landscapes and thermodynamic stability of the Omicron variant complexes. By integrating the predictions of conformational ensembles from different AlphaFold adaptations and applying statistical confidence metrics we can expand characterization of the conformational ensembles and identify functional protein conformations that determine the equilibrium dynamics for the Omicron spike complexes with the ACE2. Conformational ensembles of the Omicron RBD-ACE2 complexes obtained using AlphaFold-based approaches for modeling protein states and MD simulations are employed for accurate comparative prediction of the binding energetics revealing an excellent agreement with the experimental data. In particular, the results demonstrated that AlphaFold-generated extended conformational ensembles can produce accurate binding energies for the Omicron RBD-ACE2 complexes. The results of this study suggested complementarities and potential synergies between AlphaFold predictions of protein conformational ensembles and MD simulations showing that integrating information from both methods can potentially yield a more adequate characterization of the conformational landscapes for the Omicron RBD-ACE2 complexes. This study provides insights in the interplay between conformational dynamics and binding, showing that evolution of Omicron variants through acquisition of convergent mutational sites may leverage conformational adaptability and dynamic couplings between key binding energy hotspots to optimize ACE2 binding affinity and enable immune evasion.

'Grasshopper sign': the novel imaging of post-COVID-19 myelopathy with delayed longitudinal white matter abnormalities

Introduction

Recently, there have been a few reports of atypical post-coronavirus disease 2019 (COVID-19) myelopathy manifesting tract-specific lesions similar to those due to vitamin B12 deficiency. However, the precise characteristics of imaging or clinical course remain not well understood.

Methods

A retrospective analysis of the clinical and imaging characteristics of four patients who were referred to our hospital with a unique post-COVID-19 myelopathy was performed.

Results

Four-to-six weeks following COVID-19 infection in the summer of 2023, four middle-aged men developed paraparesis, hypo/dysesthesia and bladder/bowel disturbance, suggesting myelopathy. Although spinal MRI showed no abnormalities in the early stages, tract-specific longitudinal lesions along the dorsal and lateral columns became apparent as the symptoms progressed. Owing to the lack of MRI findings at the early stage, all cases were challenging to diagnose. However, the patients remained partially responsive to aggressive immunosuppressive therapies, even in the advanced stage.

Discussion

We termed these tract-specific longitudinal lesions in the presented case series 'Grasshopper sign' because brain coronal and spine axial MRI findings looked like a grasshopper's antennae and face. Early identification of the characteristic MRI abnormality could allow for early intervention using intensive immunosuppressive therapy, which could improve patient outcomes.

The presence of broadly neutralizing anti-SARS-CoV-2 RBD antibodies elicited by primary series and booster dose of COVID-19 vaccine

Antibody-mediated immunity plays a key role in protection against SARS-CoV-2. We characterized B-cell-derived anti-SARS-CoV-2 RBD antibody repertoires from vaccinated and infected individuals and elucidate the mechanism of action of broadly neutralizing antibodies and dissect antibodies at the epitope level. The breadth and clonality of anti-RBD B cell response varies among individuals. The majority of neutralizing antibody clones lose or exhibit reduced activities against Beta, Delta, and Omicron variants. Nevertheless, a portion of anti-RBD antibody clones that develops after a primary series or booster dose of COVID-19 vaccination exhibit broad neutralization against emerging Omicron BA.2, BA.4, BA.5, BQ.1.1, XBB.1.5 and XBB.1.16 variants. These broadly neutralizing antibodies share genetic features including a conserved usage of the IGHV3-53 and 3-9 genes and recognize three clustered epitopes of the RBD, including epitopes that partially overlap the classically defined set identified early in the pandemic. The Fab-RBD crystal and Fab-Spike complex structures corroborate the epitope grouping of antibodies and reveal the detailed binding mode of broadly neutralizing antibodies. Structure-guided mutagenesis improves binding and neutralization potency of antibody with Omicron variants via a single amino-substitution. Together, these results provide an immunological basis for partial protection against severe COVID-19 by the ancestral strain-based vaccine and indicate guidance for next generation monoclonal antibody development and vaccine design.

FLip mutations (L455F + F456L) in newly emerging VOI, JN.1: Its antibody and immune escape

During the COVID-19 pandemic, one of the biggest challenges was the continuous evolution of SARS-CoV-2 through various mutations. This has resulted in the emergence of several variants and subvariants. The escape mutations are reported as significant mutations in several variants and subvariants responsible for immune, antibody, and nAb escape. It has been reported that FLip mutations (L455F and F456L) in the spike RBD are responsible for immune evasion and antibody escape. Recently, WHO has included a new SARS-CoV-2 VOI, JN.1 lineage, a descendent of BA.2.86. The variant is reported from more than 41 countries, including France, the USA, Canada, the UK, Singapore, Sweden, and India. It contains FLip mutations in the spike protein in RBD (L455F and F456L). The risk assessment of the variant by WHO shows it has increased transmission, immune escape, and antibody escape due to the mutations. The article illustrated that FLip mutations in RBD (L455F and F456L) are responsible for augmented transmission and immune and antibody escape.

Characteristics of JN.1-derived SARS-CoV-2 subvariants SLip, FLiRT, and KP.2 in neutralization escape, infectivity and membrane fusion

SARS-CoV-2 variants derived from the immune evasive JN.1 are on the rise worldwide. Here, we investigated JN.1-derived subvariants SLip, FLiRT, and KP.2 for their ability to be neutralized by antibodies in bivalent-vaccinated human sera, XBB.1.5 monovalent-vaccinated hamster sera, sera from people infected during the BA.2.86/JN.1 wave, and class III monoclonal antibody (Mab) S309. We found that compared to parental JN.1, SLip and KP.2, and especially FLiRT, exhibit increased resistance to COVID-19 bivalent-vaccinated human sera and BA.2.86/JN.1-wave convalescent sera. Interestingly, antibodies in XBB.1.5 monovalent vaccinated hamster sera robustly neutralized FLiRT and KP.2 but had reduced efficiency for SLip. These JN.1 subvariants were resistant to neutralization by Mab S309. In addition, we investigated aspects of spike protein biology including infectivity, cell-cell fusion and processing, and found that these subvariants, especially SLip, had a decreased infectivity and membrane fusion relative to JN.1, correlating with decreased spike processing. Homology modeling revealed that L455S and F456L mutations in SLip reduced local hydrophobicity in the spike and hence its binding to ACE2. In contrast, the additional R346T mutation in FLiRT and KP.2 strengthened conformational support of the receptor-binding motif, thus counteracting the effects of L455S and F456L. These three mutations, alongside D339H, which is present in all JN.1 sublineages, alter the epitopes targeted by therapeutic Mabs, including class I and class III S309, explaining their reduced sensitivity to neutralization by sera and S309. Together, our findings provide insight into neutralization resistance of newly emerged JN.1 subvariants and suggest that future vaccine formulations should consider JN.1 spike as immunogen, although the current XBB.1.5 monovalent vaccine could still offer adequate protection.

AlphaFold2 Predictions of Conformational Ensembles and Atomistic Simulations of the SARS-CoV-2 Spike XBB Lineages Reveal Epistatic Couplings between Convergent Mutational Hotspots that Control ACE2 Affinity

In this study, we combined AlphaFold-based atomistic structural modeling, microsecond molecular simulations, mutational profiling, and network analysis to characterize binding mechanisms of the SARS-CoV-2 spike protein with the host receptor ACE2 for a series of Omicron XBB variants including XBB.1.5, XBB.1.5+L455F, XBB.1.5+F456L, and XBB.1.5+L455F+F456L. AlphaFold-based structural and dynamic modeling of SARS-CoV-2 Spike XBB lineages can accurately predict the experimental structures and characterize conformational ensembles of the spike protein complexes with the ACE2. Microsecond molecular dynamics simulations identified important differences in the conformational landscapes and equilibrium ensembles of the XBB variants, suggesting that combining AlphaFold predictions of multiple conformations with molecular dynamics simulations can provide a complementary approach for the characterization of functional protein states and binding mechanisms. Using the ensemble-based mutational profiling of protein residues and physics-based rigorous calculations of binding affinities, we identified binding energy hotspots and characterized the molecular basis underlying epistatic couplings between convergent mutational hotspots. Consistent with the experiments, the results revealed the mediating role of the Q493 hotspot in the synchronization of epistatic couplings between L455F and F456L mutations, providing a quantitative insight into the energetic determinants underlying binding differences between XBB lineages. We also proposed a network-based perturbation approach for mutational profiling of allosteric communications and uncovered the important relationships between allosteric centers mediating long-range communication and binding hotspots of epistatic couplings. The results of this study support a mechanism in which the binding mechanisms of the XBB variants may be determined by epistatic effects between convergent evolutionary hotspots that control ACE2 binding.

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.

Persistent immune imprinting occurs after vaccination with the COVID-19 XBB.1.5 mRNA booster in humans

Immune imprinting describes how the first exposure to a virus shapes immunological outcomes of subsequent exposures to antigenically related strains. Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) Omicron breakthrough infections and bivalent COVID-19 vaccination primarily recall cross-reactive memory B cells induced by prior Wuhan-Hu-1 spike mRNA vaccination rather than priming Omicron-specific naive B cells. These findings indicate that immune imprinting occurs after repeated Wuhan-Hu-1 spike exposures, but whether it can be overcome remains unclear. To understand the persistence of immune imprinting, we investigated memory and plasma antibody responses after administration of the updated XBB.1.5 COVID-19 mRNA vaccine booster. We showed that the XBB.1.5 booster elicited neutralizing antibody responses against current variants that were dominated by recall of pre-existing memory B cells previously induced by the Wuhan-Hu-1 spike. Therefore, immune imprinting persists after multiple exposures to Omicron spikes through vaccination and infection, including post XBB.1.5 booster vaccination, which will need to be considered to guide future vaccination.

Predicting Functional Conformational Ensembles and Binding Mechanisms of Convergent Evolution for SARS-CoV-2 Spike Omicron Variants Using AlphaFold2 Sequence Scanning Adaptations and Molecular Dynamics Simulations

In this study, we combined AlphaFold-based approaches for atomistic modeling of multiple protein states and microsecond molecular simulations to accurately characterize conformational ensembles and binding mechanisms of convergent evolution for the SARS-CoV-2 Spike Omicron variants BA.1, BA.2, BA.2.75, BA.3, BA.4/BA.5 and BQ.1.1. We employed and validated several different adaptations of the AlphaFold methodology for modeling of conformational ensembles including the introduced randomized full sequence scanning for manipulation of sequence variations to systematically explore conformational dynamics of Omicron Spike protein complexes with the ACE2 receptor. Microsecond atomistic molecular dynamic simulations provide a detailed characterization of the conformational landscapes and thermodynamic stability of the Omicron variant complexes. By integrating the predictions of conformational ensembles from different AlphaFold adaptations and applying statistical confidence metrics we can expand characterization of the conformational ensembles and identify functional protein conformations that determine the equilibrium dynamics for the Omicron Spike complexes with the ACE2. Conformational ensembles of the Omicron RBD-ACE2 complexes obtained using AlphaFold-based approaches for modeling protein states and molecular dynamics simulations are employed for accurate comparative prediction of the binding energetics revealing an excellent agreement with the experimental data. In particular, the results demonstrated that AlphaFold-generated extended conformational ensembles can produce accurate binding energies for the Omicron RBD-ACE2 complexes. The results of this study suggested complementarities and potential synergies between AlphaFold predictions of protein conformational ensembles and molecular dynamics simulations showing that integrating information from both methods can potentially yield a more adequate characterization of the conformational landscapes for the Omicron RBD-ACE2 complexes. This study provides insights in the interplay between conformational dynamics and binding, showing that evolution of Omicron variants through acquisition of convergent mutational sites may leverage conformational adaptability and dynamic couplings between key binding energy hotspots to optimize ACE2 binding affinity and enable immune evasion.

SARS-CoV-2 Omicron XBB lineage spike structures, conformations, antigenicity, and receptor recognition

A recombinant lineage of the SARS-CoV-2 Omicron variant, named XBB, appeared in late 2022 and evolved descendants that successively swept local and global populations. XBB lineage members were noted for their improved immune evasion and transmissibility. Here, we determine cryo-EM structures of XBB.1.5, XBB.1.16, EG.5 and EG.5.1 spike (S) ectodomains to reveal reinforced 3-RBD-down receptor inaccessible closed states mediated by interprotomer receptor binding domain (RBD) interactions previously observed in BA.1 and BA.2. Improved XBB.1.5 and XBB.1.16 RBD stability compensated for stability loss caused by early Omicron mutations, while the F456L substitution reduced EG.5 RBD stability. S1 subunit mutations had long-range impacts on conformation and epitope presentation in the S2 subunit. Our results reveal continued S protein evolution via simultaneous optimization of multiple parameters including stability, receptor binding and immune evasion, and the dramatic effects of relatively few residue substitutions in altering the S protein conformational landscape.

Spatiotemporal dynamics and epidemiological impact of SARS-CoV-2 XBB lineage dissemination in Brazil in 2023

The SARS-CoV-2 XBB is a group of highly immune-evasive lineages of the Omicron variant of concern that emerged by recombining BA.2-descendent lineages and spread worldwide during 2023. In this study, we combine SARS-CoV-2 genomic data (n = 11,065 sequences) with epidemiological data of severe acute respiratory infection (SARI) cases collected in Brazil between October 2022 and July 2023 to reconstruct the space-time dynamics and epidemiologic impact of XBB dissemination in the country. Our analyses revealed that the introduction and local emergence of lineages carrying convergent mutations within the Spike protein, especially F486P, F456L, and L455F, propelled the spread of XBB* lineages in Brazil. The average relative instantaneous reproduction numbers of XBB* + F486P, XBB* + F486P + F456L, and XBB* + F486P + F456L + L455F lineages in Brazil were estimated to be 1.24, 1.33, and 1.48 higher than that of other co-circulating lineages (mainly BQ.1*/BE*), respectively. Despite such a growth advantage, the dissemination of these XBB* lineages had a reduced impact on Brazil's epidemiological scenario concerning previous Omicron subvariants. The peak number of SARI cases from SARS-CoV-2 during the XBB wave was approximately 90%, 80%, and 70% lower than that observed during the previous BA.1*, BA.5*, and BQ.1* waves, respectively. These findings revealed the emergence of multiple XBB lineages with progressively increasing growth advantage, yet with relatively limited epidemiological impact in Brazil throughout 2023. The XBB* + F486P + F456L + L455F lineages stand out for their heightened transmissibility, warranting close monitoring in the months ahead.

IMPORTANCE

Brazil was one the most affected countries by the SARS-CoV-2 pandemic, with more than 700,000 deaths by mid-2023. This study reconstructs the dissemination of the virus in the country in the first half of 2023, a period characterized by the dissemination of descendants of XBB.1, a recombinant of Omicron BA.2 lineages evolved in late 2022. The analysis supports that XBB dissemination was marked by the continuous emergence of indigenous lineages bearing similar mutations in key sites of their Spike protein, a process followed by continuous increments in transmissibility, and without repercussions in the incidence of severe cases. Thus, the results suggest that the epidemiological impact of the spread of a SARS-CoV-2 variant is influenced by an intricate interplay of factors that extend beyond the virus's transmissibility alone. The study also underlines the need for SARS-CoV-2 genomic surveillance that allows the monitoring of its ever-shifting composition.

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

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