Cryo-EM structures of SARS-CoV-2 Omicron BA.2 spike

Biophysical principles predict fitness of SARS-CoV-2 variants

SARS-CoV-2 employs its spike protein's receptor binding domain (RBD) to enter host cells. The RBD is constantly subjected to immune responses, while requiring efficient binding to host cell receptors for successful infection. However, our understanding of how RBD's biophysical properties contribute to SARS-CoV-2's epidemiological fitness remains largely incomplete. Through a comprehensive approach, comprising large-scale sequence analysis of SARS-CoV-2 variants and the identification of a fitness function based on binding thermodynamics, we unravel the relationship between the biophysical properties of RBD variants and their contribution to viral fitness. We developed a biophysical model that uses statistical mechanics to map the molecular phenotype space, characterized by dissociation constants of RBD to ACE2, LY-CoV016, LY-CoV555, REGN10987, and S309, onto an epistatic fitness landscape. We validate our findings through experimentally measured and machine learning (ML) estimated binding affinities, coupled with infectivity data derived from population-level sequencing. Our analysis reveals that this model effectively predicts the fitness of novel RBD variants and can account for the epistatic interactions among mutations, including explaining the later reversal of Q493R. Our study sheds light on the impact of specific mutations on viral fitness and delivers a tool for predicting the future epidemiological trajectory of previously unseen or emerging low-frequency variants. These insights offer not only greater understanding of viral evolution but also potentially aid in guiding public health decisions in the battle against COVID-19 and future pandemics.

Epitopes of an antibody that neutralizes a wide range of SARS-CoV-2 variants in a conserved subdomain 1 of the spike protein

The evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continued, enabling the virus to escape from host immunity by changing its spike antigen, while biased toward the receptor-binding domain and N-terminal domain. Here, we isolated a novel pan-SARS-CoV-2 neutralizing antibody (which we named MO11) for even the recent dominators XBB.1.16 and EG.5.1, from a convalescent patient who had received three doses of an original mRNA COVID-19 vaccination. A cryo-electron microscopy analysis of the spike-MO11 complex at 2.3 Å atomic resolution revealed that it recognizes a conserved epitope hidden behind a glycan shield at N331 on subdomain 1 (SD1), holding both the N- and C-terminal segments comprising SD1. Our identification of MO11 unveiled the functional importance of SD1 for the spike's function, and we discuss the potential availability of a novel common epitope among the SARS-CoV-2 variants.IMPORTANCENovel severe acute respiratory syndrome coronavirus 2 variants with immune evasion ability are still repeatedly emerging, nonetheless, a part of immunity developed in responding to the antigen of earlier variants retains efficacy against recent variants irrespective of the numerous mutations. In exploration for the broadly effective antibodies, we identified a cross-neutralizing antibody, named MO11, from the B cells of the convalescent patient. MO11 targets a novel epitope in subdomain 1 (SD1) and was effective against all emerging variants including XBB.1.16 and EG.5.1. The neutralizing activity covering from D614G to EG.5.1 variants was explained by the conservation of the epitope, and it revealed the importance of the subdomain on regulating the function of the antigen for viral infection. Demonstrated identification of the neutralizing antibody that recognizes a conserved epitope implies basal contribution of such group of antibodies for prophylaxis against COVID-19.

The D Gene in CDR H3 Determines a Public Class of Human Antibodies to SARS-CoV-2

Public antibody responses have been found against many infectious agents. Structural convergence of public antibodies is usually determined by immunoglobulin V genes. Recently, a human antibody public class against SARS-CoV-2 was reported, where the D gene (IGHD3-22) encodes a common YYDxxG motif in heavy-chain complementarity-determining region 3 (CDR H3), which determines specificity for the receptor-binding domain (RBD). In this review, we discuss the isolation, structural characterization, and genetic analyses of this class of antibodies, which have been isolated from various cohorts of COVID-19 convalescents and vaccinees. All eleven YYDxxG antibodies with available structures target the SARS-CoV-2 RBD in a similar binding mode, where the CDR H3 dominates the interaction with antigen. The antibodies target a conserved site on the RBD that does not overlap with the receptor-binding site, but their particular angle of approach results in direct steric hindrance to receptor binding, which enables both neutralization potency and breadth. We also review the properties of CDR H3-dominant antibodies that target other human viruses. Overall, unlike most public antibodies, which are identified by their V gene usage, this newly discovered public class of YYDxxG antibodies is dominated by a D-gene-encoded motif and uncovers further opportunities for germline-targeting vaccine design.

Genome sequencing of severe acute respiratory syndrome COVID-19 EF.1 variant strain obtained from a Moroccan patient

Here, we report the identification and coding-complete genome sequence of a severe acute respiratory syndrome COVID-19 (SARS-CoV-2) strain obtained from a Moroccan patient. The detected strain EF.1 belongs to the BQ1.1 subvariant of the BA.5 Omicron variant.

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.

Altered receptor binding, antibody evasion and retention of T cell recognition by the SARS-CoV-2 XBB.1.5 spike protein

The XBB.1.5 variant of SARS-CoV-2 has rapidly achieved global dominance and exhibits a high growth advantage over previous variants. Preliminary reports suggest that the success of XBB.1.5 stems from mutations within its spike glycoprotein, causing immune evasion and enhanced receptor binding. We present receptor binding studies that demonstrate retention of binding contacts with the human ACE2 receptor and a striking decrease in binding to mouse ACE2 due to the revertant R493Q mutation. Despite extensive evasion of antibody binding, we highlight a region on the XBB.1.5 spike protein receptor binding domain (RBD) that is recognized by serum antibodies from a donor with hybrid immunity, collected prior to the emergence of the XBB.1.5 variant. T cell assays reveal high frequencies of XBB.1.5 spike-specific CD4+ and CD8+ T cells amongst donors with hybrid immunity, with the CD4+ T cells skewed towards a Th1 cell phenotype and having attenuated effector cytokine secretion as compared to ancestral spike protein-specific cells. Thus, while the XBB.1.5 variant has retained efficient human receptor binding and gained antigenic alterations, it remains susceptible to recognition by T cells induced via vaccination and previous infection.

Evolution of the SARS-CoV-2 Omicron spike

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern, first identified in November 2021, rapidly spread worldwide and diversified into several subvariants. The Omicron spike (S) protein accumulated an unprecedented number of sequence changes relative to previous variants. In this review, we discuss how Omicron S protein structural features modulate host cell receptor binding, virus entry, and immune evasion and highlight how these structural features differentiate Omicron from previous variants. We also examine how key structural properties track across the still-evolving Omicron subvariants and the importance of continuing surveillance of the S protein sequence evolution over time.

Broadly neutralizing antibody induction by non-stabilized SARS-CoV-2 Spike mRNA vaccination in nonhuman primates

Immunization with mRNA or viral vectors encoding spike with diproline substitutions (S-2P) has provided protective immunity against severe COVID-19 disease. How immunization with Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) spike elicits neutralizing antibodies (nAbs) against difficult-to-neutralize variants of concern (VOCs) remains an area of great interest. Here, we compare immunization of macaques with mRNA vaccines expressing ancestral spike either including or lacking diproline substitutions, and show the diproline substitutions were not required for protection against SARS-CoV-2 challenge or induction of broadly neutralizing B cell lineages. One group of nAbs elicited by the ancestral spike lacking diproline substitutions targeted the outer face of the receptor binding domain (RBD), neutralized all tested SARS-CoV-2 VOCs including Omicron XBB.1.5, but lacked cross-Sarbecovirus neutralization. Structural analysis showed that the macaque broad SARS-CoV-2 VOC nAbs bound to the same epitope as a human broad SARS-CoV-2 VOC nAb, DH1193. Vaccine-induced antibodies that targeted the RBD inner face neutralized multiple Sarbecoviruses, protected mice from bat CoV RsSHC014 challenge, but lacked Omicron variant neutralization. Thus, ancestral SARS-CoV-2 spike lacking proline substitutions encoded by nucleoside-modified mRNA can induce B cell lineages binding to distinct RBD sites that either broadly neutralize animal and human Sarbecoviruses or recent Omicron VOCs.

Broad receptor tropism and immunogenicity of a clade 3 sarbecovirus

Although Rhinolophus bats harbor diverse clade 3 sarbecoviruses, the structural determinants of receptor tropism along with the antigenicity of their spike (S) glycoproteins remain uncharacterized. Here, we show that the African Rhinolophus bat clade 3 sarbecovirus PRD-0038 S has a broad angiotensin-converting enzyme 2 (ACE2) usage and that receptor-binding domain (RBD) mutations further expand receptor promiscuity and enable human ACE2 utilization. We determine a cryo-EM structure of the PRD-0038 RBD bound to Rhinolophus alcyone ACE2, explaining receptor tropism and highlighting differences with SARS-CoV-1 and SARS-CoV-2. Characterization of PRD-0038 S using cryo-EM and monoclonal antibody reactivity reveals its distinct antigenicity relative to SARS-CoV-2 and identifies PRD-0038 cross-neutralizing antibodies for pandemic preparedness. PRD-0038 S vaccination elicits greater titers of antibodies cross-reacting with vaccine-mismatched clade 2 and clade 1a sarbecoviruses compared with SARS-CoV-2 S due to broader antigenic targeting, motivating the inclusion of clade 3 antigens in next-generation vaccines for enhanced resilience to viral evolution.

Structural insights into broadly neutralizing antibodies elicited by hybrid immunity against SARS-CoV-2

ABSTRACTIncreasing spread by SARS-CoV-2 Omicron variants challenges existing vaccines and broadly reactive neutralizing antibodies (bNAbs) against COVID-19. Here we determine the diversity, potency, breadth and structural insights of bNAbs derived from memory B cells of BNT162b2-vaccinee after homogeneous Omicron BA.1 breakthrough infection. The infection activates diverse memory B cell clonotypes for generating potent class I/II and III bNAbs with new epitopes mapped to the receptor-binding domain (RBD). The top eight bNAbs neutralize wildtype and BA.1 potently but display divergent IgH/IgL sequences and neuralization profiles against other variants of concern (VOCs). Two of them (P2D9 and P3E6) belonging to class III NAbs display comparable potency against BA.4/BA.5, although structural analysis reveals distinct modes of action. P3E6 neutralizes all variants tested through a unique bivalent interaction with two RBDs. Our findings provide new insights into hybrid immunity on BNT162b2-induced diverse memory B cells in response to Omicron breakthrough infection for generating diverse bNAbs with distinct structural basis.

A potent and broad-spectrum neutralizing nanobody for SARS-CoV-2 viruses, including all major Omicron strains

SARS-CoV-2 viruses are highly transmissible and immune evasive. It is critical to develop broad-spectrum prophylactic and therapeutic antibodies for potential future pandemics. Here, we used the phage display method to discover nanobodies (Nbs) for neutralizing SARS-CoV-2 viruses especially Omicron strains. The leading nanobody (Nb), namely, Nb4, with excellent physicochemical properties, can neutralize Delta and Omicron subtypes, including BA.1, BA.1.1 (BA.1 + R346K), BA.2, BA.5, BQ.1, and XBB.1. The crystal structure of Nb4 in complex with the receptor-binding domain (RBD) of BA.1 Spike protein reveals that Nb4 interacts with an epitope on the RBD overlapping with the receptor-binding motif, and thus competes with angiotensin-converting enzyme 2 (ACE2) binding. Nb4 is expected to be effective for neutralizing most recent Omicron variants, since the epitopes are evolutionarily conserved among them. Indeed, trivalent Nb4 interacts with the XBB1.5 Spike protein with low nM affinity and competes for ACE2 binding. Prophylactic and therapeutic experiments in mice indicated that Nb4 could reduce the Omicron virus loads in the lung. In particular, in prophylactic experiments, intranasal administration of multivalent Nb4 completely protected mice from Omicron infection. Taken together, these results demonstrated that Nb4 could serve as a potent and broad-spectrum prophylactic and therapeutic Nb for COVID-19.

Cooperative and structural relationships of the trimeric Spike with infectivity and antibody escape of the strains Delta (B.1.617.2) and Omicron (BA.2, BA.5, and BQ.1)

Herein, we conducted simulations of trimeric Spike from several SARS-CoV-2 variants of concern (Delta and Omicron sub-variants BA.2, BA.5, and BQ.1) and investigated the mechanisms by which specific mutations confer resistance to neutralizing antibodies. We observed that the mutations primarily affect the cooperation between protein domains within and between protomers. The substitutions K417N and L452R expand hydrogen bonding interactions, reducing their interaction with neutralizing antibodies. By interacting with nearby residues, the K444T and N460K mutations in the SpikeBQ.1 variant potentially reduces solvent exposure, thereby promoting resistance to antibodies. We also examined the impact of D614G, P681R, and P681H substitutions on Spike protein structure that may be related to infectivity. The D614G substitution influences communication between a glycine residue and neighboring domains, affecting the transition between up- and -down RBD states. The P681R mutation, found in the Delta variant, enhances correlations between protein subunits, while the P681H mutation in Omicron sub-variants weakens long-range interactions that may be associated with reduced fusogenicity. Using a multiple linear regression model, we established a connection between inter-protomer communication and loss of sensitivity to neutralizing antibodies. Our findings underscore the importance of structural communication between protein domains and provide insights into potential mechanisms of immune evasion by SARS-CoV-2. Overall, this study deepens our understanding of how specific mutations impact SARS-CoV-2 infectivity and shed light on how the virus evades the immune system.

Engineered immunogens to elicit antibodies against conserved coronavirus epitopes

Immune responses to SARS-CoV-2 primarily target the receptor binding domain of the spike protein, which continually mutates to escape acquired immunity. Other regions in the spike S2 subunit, such as the stem helix and the segment encompassing residues 815-823 adjacent to the fusion peptide, are highly conserved across sarbecoviruses and are recognized by broadly reactive antibodies, providing hope that vaccines targeting these epitopes could offer protection against both current and emergent viruses. Here we employ computational modeling to design scaffolded immunogens that display the spike 815-823 peptide and the stem helix epitopes without the distracting and immunodominant receptor binding domain. These engineered proteins bind with high affinity and specificity to the mature and germline versions of previously identified broadly protective human antibodies. Epitope scaffolds interact with both sera and isolated monoclonal antibodies with broadly reactivity from individuals with pre-existing SARS-CoV-2 immunity. When used as immunogens, epitope scaffolds elicit sera with broad betacoronavirus reactivity and protect as "boosts" against live virus challenge in mice, illustrating their potential as components of a future pancoronavirus vaccine.

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.

Infection and Vaccine Induced Spike Antibody Responses Against SARS-CoV-2 Variants of Concern in COVID-19-Naïve Children and Adults

Although a more efficient adaptive humoral immune response has been proposed to underlie the usually favorable outcome of pediatric COVID-19, the breadth of viral and vaccine cross-reactivity toward the ever-mutating Spike protein among variants of concern (VOCs) has not yet been compared between children and adults. We assessed antibodies to conformational Spike in COVID-19-naïve children and adults vaccinated by BNT162b2 and ChAdOx1, and naturally infected with SARS-CoV-2 Early Clade, Delta, and Omicron. Sera were analyzed against Spike including naturally occurring VOCs Alpha, Beta, Gamma, Delta, and Omicron BA.1, BA.2, BA.5, BQ.1.1, BA2.75.2, and XBB.1, and variants of interest Epsilon, Kappa, Eta, D.2, and artificial mutant Spikes. There was no notable difference between breadth and longevity of antibody against VOCs in children and adults. Vaccinated individuals displayed similar immunoreactivity profiles across variants compared with naturally infected individuals. Delta-infected patients had an enhanced cross-reactivity toward Delta and earlier VOCs compared to patients infected by Early Clade SARS-CoV-2. Although Omicron BA.1, BA.2, BA.5, BQ.1.1, BA2.75.2, and XBB.1 antibody titers were generated after Omicron infection, cross-reactive binding against Omicron subvariants was reduced across all infection, immunization, and age groups. Some mutations, such as 498R and 501Y, epistatically combined to enhance cross-reactive binding, but could not fully compensate for antibody-evasive mutations within the Omicron subvariants tested. Our results reveal important molecular features central to the generation of high antibody titers and broad immunoreactivity that should be considered in future vaccine design and global serosurveillance in the context of limited vaccine boosters available to the pediatric population.

Comparative Analysis of Conformational Dynamics and Systematic Characterization of Cryptic Pockets in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 Spike Complexes with the ACE2 Host Receptor: Confluence of Binding and Structural Plasticity in Mediating Networks of Conserved Allosteric Sites

In the current study, we explore coarse-grained simulations and atomistic molecular dynamics together with binding energetics scanning and cryptic pocket detection in a comparative examination of conformational landscapes and systematic characterization of allosteric binding sites in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 spike full-length trimer complexes with the host receptor ACE2. Microsecond simulations, Markov state models and mutational scanning of binding energies of the SARS-CoV-2 BA.2 and BA.2.75 receptor binding domain complexes revealed the increased thermodynamic stabilization of the BA.2.75 variant and significant dynamic differences between these Omicron variants. Molecular simulations of the SARS-CoV-2 Omicron spike full-length trimer complexes with the ACE2 receptor complemented atomistic studies and enabled an in-depth analysis of mutational and binding effects on conformational dynamic and functional adaptability of the Omicron variants. Despite considerable structural similarities, Omicron variants BA.2, BA.2.75 and XBB.1 can induce unique conformational dynamic signatures and specific distributions of the conformational states. Using conformational ensembles of the SARS-CoV-2 Omicron spike trimer complexes with ACE2, we conducted a comprehensive cryptic pocket screening to examine the role of Omicron mutations and ACE2 binding on the distribution and functional mechanisms of the emerging allosteric binding sites. This analysis captured all experimentally known allosteric sites and discovered networks of inter-connected and functionally relevant allosteric sites that are governed by variant-sensitive conformational adaptability of the SARS-CoV-2 spike structures. The results detailed how ACE2 binding and Omicron mutations in the BA.2, BA.2.75 and XBB.1 spike complexes modulate the distribution of conserved and druggable allosteric pockets harboring functionally important regions. The results are significant for understanding the functional roles of druggable cryptic pockets that can be used for allostery-mediated therapeutic intervention targeting conformational states of the Omicron variants.

Impact of mutations on the plant-based production of recombinant SARS-CoV-2 RBDs

Subunit vaccines based on recombinant viral antigens are valuable interventions to fight existing and evolving viruses and can be produced at large-scale in plant-based expression systems. The recombinant viral antigens are often derived from glycosylated envelope proteins of the virus and glycosylation plays an important role for the immunogenicity by shielding protein epitopes. The receptor-binding domain (RBD) of the SARS-CoV-2 spike is a principal target for vaccine development and has been produced in plants, but the yields of recombinant RBD variants were low and the role of the N-glycosylation in RBD from different SARS-CoV-2 variants of concern is less studied. Here, we investigated the expression and glycosylation of six different RBD variants transiently expressed in leaves of Nicotiana benthamiana. All of the purified RBD variants were functional in terms of receptor binding and displayed almost full N-glycan occupancy at both glycosylation sites with predominately complex N-glycans. Despite the high structural sequence conservation of the RBD variants, we detected a variation in yield which can be attributed to lower expression and differences in unintentional proteolytic processing of the C-terminal polyhistidine tag used for purification. Glycoengineering towards a human-type complex N-glycan profile with core α1,6-fucose, showed that the reactivity of the neutralizing antibody S309 differs depending on the N-glycan profile and the RBD variant.

Exploring Conformational Landscapes and Cryptic Binding Pockets in Distinct Functional States of the SARS-CoV-2 Omicron BA.1 and BA.2 Trimers: Mutation-Induced Modulation of Protein Dynamics and Network-Guided Prediction of Variant-Specific Allosteric Binding Sites

A significant body of experimental structures of SARS-CoV-2 spike trimers for the BA.1 and BA.2 variants revealed a considerable plasticity of the spike protein and the emergence of druggable binding pockets. Understanding the interplay of conformational dynamics changes induced by the Omicron variants and the identification of cryptic dynamic binding pockets in the S protein is of paramount importance as exploring broad-spectrum antiviral agents to combat the emerging variants is imperative. In the current study, we explore conformational landscapes and characterize the universe of binding pockets in multiple open and closed functional spike states of the BA.1 and BA.2 Omicron variants. By using a combination of atomistic simulations, a dynamics network analysis, and an allostery-guided network screening of binding pockets in the conformational ensembles of the BA.1 and BA.2 spike conformations, we identified all experimentally known allosteric sites and discovered significant variant-specific differences in the distribution of binding sites in the BA.1 and BA.2 trimers. This study provided a structural characterization of the predicted cryptic pockets and captured the experimentally known allosteric sites, revealing the critical role of conformational plasticity in modulating the distribution and cross-talk between functional binding sites. We found that mutational and dynamic changes in the BA.1 variant can induce the remodeling and stabilization of a known druggable pocket in the N-terminal domain, while this pocket is drastically altered and may no longer be available for ligand binding in the BA.2 variant. Our results predicted the experimentally known allosteric site in the receptor-binding domain that remains stable and ranks as the most favorable site in the conformational ensembles of the BA.2 variant but could become fragmented and less probable in BA.1 conformations. We also uncovered several cryptic pockets formed at the inter-domain and inter-protomer interface, including functional regions of the S2 subunit and stem helix region, which are consistent with the known role of pocket residues in modulating conformational transitions and antibody recognition. The results of this study are particularly significant for understanding the dynamic and network features of the universe of available binding pockets in spike proteins, as well as the effects of the Omicron-variant-specific modulation of preferential druggable pockets. The exploration of predicted druggable sites can present a new and previously underappreciated opportunity for therapeutic interventions for Omicron variants through the conformation-selective and variant-specific targeting of functional sites involved in allosteric changes.

Interaction of SARS-CoV-2 with host cells and antibodies: experiment and simulation

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the devastating global COVID-19 pandemic announced by WHO in March 2020. Through unprecedented scientific effort, several vaccines, drugs and antibodies have been developed, saving millions of lives, but the fight against COVID-19 continues as immune escape variants of concern such as Delta and Omicron emerge. To develop more effective treatments and to elucidate the side effects caused by vaccines and therapeutic agents, a deeper understanding of the molecular interactions of SARS-CoV-2 with them and human cells is required. With special interest in computational approaches, we will focus on the structure of SARS-CoV-2 and the interaction of its spike protein with human angiotensin-converting enzyme-2 (ACE2) as a prime entry point of the virus into host cells. In addition, other possible viral receptors will be considered. The fusion of viral and human membranes and the interaction of the spike protein with antibodies and nanobodies will be discussed, as well as the effect of SARS-CoV-2 on protein synthesis in host cells.

Engineered Immunogens to Elicit Antibodies Against Conserved Coronavirus Epitopes

Immune responses to SARS-CoV-2 primarily target the receptor binding domain of the spike protein, which continually mutates to escape acquired immunity. Other regions in the spike S2 subunit, such as the stem helix and the segment encompassing residues 815-823 adjacent to the fusion peptide, are highly conserved across sarbecoviruses and are recognized by broadly reactive antibodies, providing hope that vaccines targeting these epitopes could offer protection against both current and emergent viruses. Here we employed computational modeling to design scaffolded immunogens that display the spike 815-823 peptide and the stem helix epitopes without the distracting and immunodominant RBD. These engineered proteins bound with high affinity and specificity to the mature and germline versions of previously identified broadly protective human antibodies. Epitope scaffolds interacted with both sera and isolated monoclonal antibodies with broadly reactivity from individuals with pre-existing SARS-CoV-2 immunity. When used as immunogens, epitope scaffolds elicited sera with broad betacoronavirus reactivity and protected as "boosts" against live virus challenge in mice, illustrating their potential as components of a future pancoronavirus vaccine.

Broad receptor tropism and immunogenicity of a clade 3 sarbecovirus

Although Rhinolophus bats harbor diverse clade 3 sarbecoviruses, the structural determinants of receptor tropism along with the antigenicity of their spike (S) glycoproteins remain uncharacterized. Here, we show that the African Rinolophus bat clade 3 sarbecovirus PRD-0038 S has a broad ACE2 usage and that RBD mutations further expand receptor promiscuity and enable human ACE2 utilization. We determined a cryoEM structure of the PRD-0038 RBD bound to R. alcyone ACE2, explaining receptor tropism and highlighting differences with SARS-CoV-1 and SARS-CoV-2. Characterization of PRD-0038 S using cryoEM and monoclonal antibody reactivity revealed its distinct antigenicity relative to SARS-CoV-2 and identified PRD-0038 cross-neutralizing antibodies for pandemic preparedness. PRD-0038 S vaccination elicited greater titers of antibodies cross-reacting with vaccine-mismatched clade 2 and clade 1a sarbecoviruses compared to SARS-CoV-2 S due to broader antigenic targeting, motivating the inclusion of clade 3 antigens in next-generation vaccines for enhanced resilience to viral evolution.

Examining Functional Linkages Between Conformational Dynamics, Protein Stability and Evolution of Cryptic Binding Pockets in the SARS-CoV-2 Omicron Spike Complexes with the ACE2 Host Receptor: Recombinant Omicron Variants Mediate Variability of Conserved Allosteric Sites and Binding Epitopes

In the current study, we explore coarse-grained simulations and atomistic molecular dynamics together with binding energetics scanning and cryptic pocket detection in a comparative examination of conformational landscapes and systematic characterization of allosteric binding sites in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB.1 spike full-length trimer complexes with the host receptor ACE2. Microsecond simulations, Markov state models and mutational scanning of binding energies of the SARS-CoV-2 BA.2 and BA.2.75 receptor binding domain complexes revealed the increased thermodynamic stabilization of the BA.2.75 variant and significant dynamic differences between these Omicron variants. Molecular simulations of the SARS-CoV-2 Omicron spike full length trimer complexes with the ACE2 receptor complemented atomistic studies and enabled an in-depth analysis of mutational and binding effects on conformational dynamic and functional adaptability of the Omicron variants. Despite considerable structural similarities, Omicron variants BA.2, BA.2.75 and XBB.1 can induce unique conformational dynamic signatures and specific distributions of the conformational states. Using conformational ensembles of the SARS-CoV-2 Omicron spike trimer complexes with ACE2, we conducted a comprehensive cryptic pocket screening to examine the role of Omicron mutations and ACE2 binding on the distribution and functional mechanisms of the emerging allosteric binding sites. This analysis captured all experimentally known allosteric sites and discovered networks of inter-connected and functionally relevant allosteric sites that are governed by variant-sensitive conformational adaptability of the SARS-CoV-2 spike structures. The results detailed how ACE2 binding and Omicron mutations in the BA.2, BA.2.75 and XBB.1 spike complexes modulate the distribution of conserved and druggable allosteric pockets harboring functionally important regions. The results of are significant for understanding functional roles of druggable cryptic pockets that can be used for allostery-mediated therapeutic intervention targeting conformational states of the Omicron variants.

Structural Framework for Analysis of CD4+ T-Cell Epitope Dominance in Viral Fusion Proteins

Antigen conformation shapes CD4+ T-cell specificity through mechanisms of antigen processing, and the consequences for immunity may rival those from conformational effects on antibody specificity. CD4+ T cells initiate and control immunity to pathogens and cancer and are at least partly responsible for immunopathology associated with infection, autoimmunity, and allergy. The primary trigger for CD4+ T-cell maturation is the presentation of an epitope peptide in the MHC class II antigen-presenting protein (MHCII), most commonly on an activated dendritic cell, and then the T-cell responses are recalled by subsequent presentations of the epitope peptide by the same or other antigen-presenting cells. Peptide presentation depends on the proteolytic fragmentation of the antigen in an endosomal/lysosomal compartment and concomitant loading of the fragments into the MHCII, a multistep mechanism called antigen processing and presentation. Although the role of peptide affinity for MHCII has been well studied, the role of proteolytic fragmentation has received less attention. In this Perspective, we will briefly summarize evidence that antigen resistance to unfolding and proteolytic fragmentation shapes the specificity of the CD4+ T-cell response to selected viral envelope proteins, identify several remarkable examples in which the immunodominant CD4+ epitopes most likely depend on the interaction of processing machinery with antigen conformation, and outline how knowledge of antigen conformation can inform future efforts to design vaccines.

Neutralization, effector function and immune imprinting of Omicron variants

Currently circulating SARS-CoV-2 variants have acquired convergent mutations at hot spots in the receptor-binding domain1 (RBD) of the spike protein. The effects of these mutations on viral infection and transmission and the efficacy of vaccines and therapies remains poorly understood. Here we demonstrate that recently emerged BQ.1.1 and XBB.1.5 variants bind host ACE2 with high affinity and promote membrane fusion more efficiently than earlier Omicron variants. Structures of the BQ.1.1, XBB.1 and BN.1 RBDs bound to the fragment antigen-binding region of the S309 antibody (the parent antibody for sotrovimab) and human ACE2 explain the preservation of antibody binding through conformational selection, altered ACE2 recognition and immune evasion. We show that sotrovimab binds avidly to all Omicron variants, promotes Fc-dependent effector functions and protects mice challenged with BQ.1.1 and hamsters challenged with XBB.1.5. Vaccine-elicited human plasma antibodies cross-react with and trigger effector functions against current Omicron variants, despite a reduced neutralizing activity, suggesting a mechanism of protection against disease, exemplified by S309. Cross-reactive RBD-directed human memory B cells remained dominant even after two exposures to Omicron spikes, underscoring the role of persistent immune imprinting.

Adaptive variations in SARS-CoV-2 spike proteins: effects on distinct virus-cell entry stages

Evolved SARS-CoV-2 variants of concern (VOCs) spread through human populations in succession. Major virus variations are in the entry-facilitating viral spike (S) proteins; Omicron VOCs have 29-40 S mutations relative to ancestral D614G viruses. The impacts of this Omicron divergence on S protein structure, antigenicity, cell entry pathways, and pathogenicity have been extensively evaluated, yet gaps remain in correlating specific alterations with S protein functions. In this study, we compared the functions of ancestral D614G and Omicron VOCs using cell-free assays that can reveal differences in several distinct steps of the S-directed virus entry process. Relative to ancestral D614G, Omicron BA.1 S proteins were hypersensitized to receptor activation, to conversion into intermediate conformational states, and to membrane fusion-activating proteases. We identified mutations conferring these changes in S protein character by evaluating domain-exchanged D614G/Omicron recombinants in the cell-free assays. Each of the three functional alterations was mapped to specific S protein domains, with the recombinants providing insights on inter-domain interactions that fine-tune S-directed virus entry. Our results provide a structure-function atlas of the S protein variations that may promote the transmissibility and infectivity of current and future SARS-CoV-2 VOCs. IMPORTANCE Continuous SARS-CoV-2 adaptations generate increasingly transmissible variants. These succeeding variants show ever-increasing evasion of suppressive antibodies and host factors, as well as increasing invasion of susceptible host cells. Here, we evaluated the adaptations enhancing invasion. We used reductionist cell-free assays to compare the entry steps of ancestral (D614G) and Omicron (BA.1) variants. Relative to D614G, Omicron entry was distinguished by heightened responsiveness to entry-facilitating receptors and proteases and by enhanced formation of intermediate states that execute virus-cell membrane fusion. We found that these Omicron-specific characteristics arose from mutations in specific S protein domains and subdomains. The results reveal the inter-domain networks controlling S protein dynamics and efficiencies of entry steps, and they offer insights on the evolution of SARS-CoV-2 variants that arise and ultimately dominate infections worldwide.

Structural evolution of SARS-CoV-2 omicron in human receptor recognition

Understanding the evolutionary strategies of the SARS-CoV-2 omicron variant is crucial for comprehending the COVID-19 pandemic and preventing future coronavirus pandemics. In this study, we determined the crystal structures of the receptor-binding domains (RBDs) from currently circulating omicron subvariants XBB.1 and XBB.1.5 (also the emerging XBB.1.9.1), each complexed with human ACE2. We studied how individual RBD residues evolved structurally in omicron subvariants, specifically how they adapted to human ACE2. Our findings revealed that residues 493 and 496, which exhibited good human ACE2 adaptation in pre-omicron variants, evolved to poor adaptation in early omicron subvariants (but with good adaption to mouse ACE2) and then reverted to good adaptation in recent omicron subvariants. This result is consistent with the hypothesis that non-human animals facilitated the evolution of early omicron subvariants. Additionally, residue 486, which exhibited good human ACE2 adaptation in early omicron subvariants, evolved to poor adaptation in later omicron subvariants and then returned to good adaptation in recent omicron subvariants. This result is consistent with the hypothesis that immune evasion facilitated the evolution of later omicron subvariants. Thus, our study suggests that both non-human animals and immune evasion may have contributed to driving omicron evolution at different stages of the pandemic. IMPORTANCE The sudden emergence and continued evolution of the SARS-CoV-2 omicron variant have left many mysteries unanswered, such as the origin of early omicron subvariants and the factors driving omicron evolution. To address these questions, we studied the crystal structures of human ACE2-bound receptor-binding domains (RBDs) from omicron subvariants XBB.1 and XBB.1.5 (XBB.1.9.1). Our in-depth structural analysis sheds light on how specific RBD mutations adapt to either human or mouse ACE2 and suggests non-human animals and immune evasion may have influenced omicron evolution during different stages of the pandemic. These findings provide valuable insights into the mechanisms underlying omicron evolution, deepen our understanding of the COVID-19 pandemic, and have significant implications for preventing future coronavirus pandemics.

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variant Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, microsecond molecular dynamics simulations, and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the functional conformational states and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant, which can be contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of the conformational states. The results suggested that variant-specific changes of the conformational mobility in the functional interfacial loops of the receptor-binding domain in the SARS-CoV-2 spike protein can be fine-tuned through crosstalk between convergent mutations which could provide an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulation of conformational plasticity and regulation of allosteric communications. This study also revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions.

Probing conformational landscapes of binding and allostery in the SARS-CoV-2 omicron variant complexes using microsecond atomistic simulations and perturbation-based profiling approaches: hidden role of omicron mutations as modulators of allosteric signaling and epistatic relationships

In this study, we systematically examine the conformational dynamics, binding and allosteric communications in the Omicron BA.1, BA.2, BA.3 and BA.4/BA.5 spike protein complexes with the ACE2 host receptor using molecular dynamics simulations and perturbation-based network profiling approaches. Microsecond atomistic simulations provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the BA.2 variant which can be contrasted with the BA.4/BA.5 variants inducing a significant mobility of the complexes. Using the dynamics-based mutational scanning of spike residues, we identified structural stability and binding affinity hotspots in the Omicron complexes. Perturbation response scanning and network-based mutational profiling approaches probed the effect of the Omicron mutations on allosteric interactions and communications in the complexes. The results of this analysis revealed specific roles of Omicron mutations as conformationally plastic and evolutionary adaptable modulators of binding and allostery which are coupled to the major regulatory positions through interaction networks. Through perturbation network scanning of allosteric residue potentials in the Omicron variant complexes performed in the background of the original strain, we characterized regions of epistatic couplings that are centered around the binding affinity hotspots N501Y and Q498R. Our results dissected the vital role of these epistatic centers in regulating protein stability, efficient ACE2 binding and allostery which allows for accumulation of multiple Omicron immune escape mutations at other sites. Through integrative computational approaches, this study provides a systematic analysis of the effects of Omicron mutations on thermodynamics, binding and allosteric signaling in the complexes with ACE2 receptor.

SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid

The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, which is also found in some other coronaviruses, e.g. SARS-CoV and MERS-CoV. The occupancy of the FA site by linoleic acid (LA) reduces infectivity by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the allosteric responses of spike variants to LA removal. D-NEMD simulations show that the FA site is coupled to other functional regions of the protein, e.g. the receptor-binding motif (RBM), N-terminal domain (NTD), furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to these functional regions. The comparison between the wild-type spike and four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their responses to LA removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the RBM and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most different variant, exhibiting significant differences in the RBM, NTD, V622-L629, and furin cleavage site. These differences in the allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of LA on SARS-CoV-2 variants, including emerging variants, is warranted.

Jointly modeling deep mutational scans identifies shifted mutational effects among SARS-CoV-2 spike homologs

Deep mutational scanning (DMS) is a high-throughput experimental technique that measures the effects of thousands of mutations to a protein. These experiments can be performed on multiple homologs of a protein or on the same protein selected under multiple conditions. It is often of biological interest to identify mutations with shifted effects across homologs or conditions. However, it is challenging to determine if observed shifts arise from biological signal or experimental noise. Here, we describe a method for jointly inferring mutational effects across multiple DMS experiments while also identifying mutations that have shifted in their effects among experiments. A key aspect of our method is to regularize the inferred shifts, so that they are nonzero only when strongly supported by the data. We apply this method to DMS experiments that measure how mutations to spike proteins from SARS-CoV-2 variants (Delta, Omicron BA.1, and Omicron BA.2) affect cell entry. Most mutational effects are conserved between these spike homologs, but a fraction have markedly shifted. We experimentally validate a subset of the mutations inferred to have shifted effects, and confirm differences of > 1,000-fold in the impact of the same mutation on spike-mediated viral infection across spikes from different SARS-CoV-2 variants. Overall, our work establishes a general approach for comparing sets of DMS experiments to identify biologically important shifts in mutational effects.

Mechanism and evolution of human ACE2 binding by SARS-CoV-2 spike

Spike glycoprotein of SARS-CoV-2 mediates viral entry into host cells by facilitating virus attachment and membrane fusion. ACE2 is the main receptor of SARS-CoV-2 and its interaction with spike has shaped the virus' emergence from an animal reservoir and subsequent evolution in the human host. Many structural studies on the spike:ACE2 interaction have provided insights into mechanisms driving viral evolution during the on-going pandemic. This review describes the molecular basis of spike binding to ACE2, outlines mechanisms that have optimised this interaction during viral evolution, and suggests directions for future research.

Potent pan huACE2-dependent sarbecovirus neutralizing monoclonal antibodies isolated from a BNT162b2-vaccinated SARS survivor

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern such as Omicron hampered efforts in controlling the ongoing coronavirus disease 2019 pandemic due to their ability to escape neutralizing antibodies induced by vaccination or prior infection, highlighting the need to develop broad-spectrum vaccines and therapeutics. Most human monoclonal antibodies (mAbs) reported to date have not demonstrated true pan-sarbecovirus neutralizing breadth especially against animal sarbecoviruses. Here, we report the isolation and characterization of highly potent mAbs targeting the receptor binding domain (RBD) of huACE2-dependent sarbecovirus from a SARS-CoV survivor vaccinated with BNT162b2. Among the six mAbs identified, one (E7) showed better huACE2-dependent sarbecovirus neutralizing potency and breadth than any other mAbs reported to date. Mutagenesis and cryo-electron microscopy studies indicate that these mAbs have a unique RBD contact footprint and that E7 binds to a quaternary structure-dependent epitope.

Spike-protein proteolytic antibodies in COVID-19 convalescent plasma contribute to SARS-CoV-2 neutralization

Understanding the mechanisms of antibody-mediated neutralization of SARS-CoV-2 is critical in combating the COVID-19 pandemic. Based on previous reports of antibody catalysis, we investigated the proteolysis of spike (S) by antibodies in COVID-19 convalescent plasma (CCP) and its contribution to viral neutralization. Quenched fluorescent peptides were designed based on S epitopes to sensitively detect antibody-mediated proteolysis. We observed epitope cleavage by CCP from different donors which persisted when plasma was heat-treated or when IgG was isolated from plasma. Further, purified CCP antibodies proteolyzed recombinant S domains, as well as authentic viral S. Cleavage of S variants suggests CCP antibody-mediated proteolysis is a durable phenomenon despite antigenic drift. We differentiated viral neutralization occurring via direct interference with receptor binding from that occurring by antibody-mediated proteolysis, demonstrating that antibody catalysis enhanced neutralization. These results suggest that antibody-catalyzed damage of S is an immunologically relevant function of neutralizing antibodies against SARS-CoV-2.

Structural and functional characteristics of the SARS-CoV-2 Omicron subvariant BA.2 spike protein

The Omicron subvariant BA.2 has become the dominant circulating strain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in many countries. Here, we have characterized structural, functional and antigenic properties of the full-length BA.2 spike (S) protein and compared replication of the authentic virus in cell culture and an animal model with previously prevalent variants. BA.2 S can fuse membranes slightly more efficiently than Omicron BA.1, but still less efficiently than other previous variants. Both BA.1 and BA.2 viruses replicated substantially faster in animal lungs than the early G614 (B.1) strain in the absence of pre-existing immunity, possibly explaining the increased transmissibility despite their functionally compromised spikes. As in BA.1, mutations in the BA.2 S remodel its antigenic surfaces, leading to strong resistance to neutralizing antibodies. These results suggest that both immune evasion and replicative advantage may contribute to the heightened transmissibility of the Omicron subvariants.

Structural Basis for the Enhanced Infectivity and Immune Evasion of Omicron Subvariants

The Omicron variants of SARS-CoV-2 have emerged as the dominant strains worldwide, causing the COVID-19 pandemic. Each Omicron subvariant contains at least 30 mutations on the spike protein (S protein) compared to the original wild-type (WT) strain. Here we report the cryo-EM structures of the trimeric S proteins from the BA.1, BA.2, BA.3, and BA.4/BA.5 subvariants, with BA.4 and BA.5 sharing the same S protein mutations, each in complex with the surface receptor ACE2. All three receptor-binding domains of the S protein from BA.2 and BA.4/BA.5 are "up", while the BA.1 S protein has two "up" and one "down". The BA.3 S protein displays increased heterogeneity, with the majority in the all "up" RBD state. The different conformations preferences of the S protein are consistent with their varied transmissibility. By analyzing the position of the glycan modification on Asn343, which is located at the S309 epitopes, we have uncovered the underlying immune evasion mechanism of the Omicron subvariants. Our findings provide a molecular basis of high infectivity and immune evasion of Omicron subvariants, thereby offering insights into potential therapeutic interventions against SARS-CoV-2 variants.

Markov State Models and Perturbation-Based Approaches Reveal Distinct Dynamic Signatures and Hidden Allosteric Pockets in the Emerging SARS-Cov-2 Spike Omicron Variants Complexes with the Host Receptor: The Interplay of Dynamics and Convergent Evolution Modulates Allostery and Functional Mechanisms

The new generation of SARS-CoV-2 Omicron variants displayed a significant growth advantage and the increased viral fitness by acquiring convergent mutations, suggesting that the immune pressure can promote convergent evolution leading to the sudden acceleration of SARS-CoV-2 evolution. In the current study, we combined structural modeling, extensive microsecond MD simulations and Markov state models to characterize conformational landscapes and identify specific dynamic signatures of the SARS-CoV-2 spike complexes with the host receptor ACE2 for the recently emerged highly transmissible XBB.1, XBB.1.5, BQ.1, and BQ.1.1 Omicron variants. Microsecond simulations and Markovian modeling provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the XBB.1.5 subvariant which is contrasted to more dynamic BQ.1 and BQ.1.1 subvariants. Despite considerable structural similarities, Omicron mutations can induce unique dynamic signatures and specific distributions of conformational states. The results suggested that variant-specific changes of conformational mobility in the functional interfacial loops of the spike receptor binding domain can be fine-tuned through cross-talk between convergent mutations thereby providing an evolutionary path for modulation of immune escape. By combining atomistic simulations and Markovian modeling analysis with perturbation-based approaches, we determined important complementary roles of convergent mutation sites as effectors and receivers of allosteric signaling involved in modulating conformational plasticity at the binding interface and regulating allosteric responses. This study also characterized the dynamics-induced evolution of allosteric pockets in the Omicron complexes that revealed hidden allosteric pockets and suggested that convergent mutation sites could control evolution and distribution of allosteric pockets through modulation of conformational plasticity in the flexible adaptable regions. Through integrative computational approaches, this investigation provides a systematic analysis and comparison of the effects of Omicron subvariants on conformational dynamics and allosteric signaling in the complexes with the ACE2 receptor.

For Table of Contents Use Only

Virological characteristics of the SARS-CoV-2 XBB variant derived from recombination of two Omicron subvariants

In late 2022, SARS-CoV-2 Omicron subvariants have become highly diversified, and XBB is spreading rapidly around the world. Our phylogenetic analyses suggested that XBB emerged through the recombination of two cocirculating BA.2 lineages, BJ.1 and BM.1.1.1 (a progeny of BA.2.75), during the summer of 2022. XBB.1 is the variant most profoundly resistant to BA.2/5 breakthrough infection sera to date and is more fusogenic than BA.2.75. The recombination breakpoint is located in the receptor-binding domain of spike, and each region of the recombinant spike confers immune evasion and increases fusogenicity. We further provide the structural basis for the interaction between XBB.1 spike and human ACE2. Finally, the intrinsic pathogenicity of XBB.1 in male hamsters is comparable to or even lower than that of BA.2.75. Our multiscale investigation provides evidence suggesting that XBB is the first observed SARS-CoV-2 variant to increase its fitness through recombination rather than substitutions.

Balancing Functional Tradeoffs between Protein Stability and ACE2 Binding in the SARS-CoV-2 Omicron BA.2, BA.2.75 and XBB Lineages: Dynamics-Based Network Models Reveal Epistatic Effects Modulating Compensatory Dynamic and Energetic Changes

Evolutionary and functional studies suggested that the emergence of the Omicron variants can be determined by multiple fitness trade-offs including the immune escape, binding affinity for ACE2, conformational plasticity, protein stability and allosteric modulation. In this study, we systematically characterize conformational dynamics, structural stability and binding affinities of the SARS-CoV-2 Spike Omicron complexes with the host receptor ACE2 for BA.2, BA.2.75, XBB.1 and XBB.1.5 variants. We combined multiscale molecular simulations and dynamic analysis of allosteric interactions together with the ensemble-based mutational scanning of the protein residues and network modeling of epistatic interactions. This multifaceted computational study characterized molecular mechanisms and identified energetic hotspots that can mediate the predicted increased stability and the enhanced binding affinity of the BA.2.75 and XBB.1.5 complexes. The results suggested a mechanism driven by the stability hotspots and a spatially localized group of the Omicron binding affinity centers, while allowing for functionally beneficial neutral Omicron mutations in other binding interface positions. A network-based community model for the analysis of epistatic contributions in the Omicron complexes is proposed revealing the key role of the binding hotspots R498 and Y501 in mediating community-based epistatic couplings with other Omicron sites and allowing for compensatory dynamics and binding energetic changes. The results also showed that mutations in the convergent evolutionary hotspot F486 can modulate not only local interactions but also rewire the global network of local communities in this region allowing the F486P mutation to restore both the stability and binding affinity of the XBB.1.5 variant which may explain the growth advantages over the XBB.1 variant. The results of this study are consistent with a broad range of functional studies rationalizing functional roles of the Omicron mutation sites that form a coordinated network of hotspots enabling a balance of multiple fitness tradeoffs and shaping up a complex functional landscape of virus transmissibility.

Probing Conformational Landscapes of Binding and Allostery in the SARS-CoV-2 Omicron Variant Complexes Using Microsecond Atomistic Simulations and Perturbation-Based Profiling Approaches: Hidden Role of Omicron Mutations as Modulators of Allosteric Signaling and Epistatic Relationships

In this study, we systematically examine the conformational dynamics, binding and allosteric communications in the Omicron BA.1, BA.2, BA.3 and BA.4/BA.5 complexes with the ACE2 host receptor using molecular dynamics simulations and perturbation-based network profiling approaches. Microsecond atomistic simulations provided a detailed characterization of the conformational landscapes and revealed the increased thermodynamic stabilization of the BA.2 variant which is contrasted with the BA.4/BA.5 variants inducing a significant mobility of the complexes. Using ensemble-based mutational scanning of binding interactions, we identified binding affinity and structural stability hotspots in the Omicron complexes. Perturbation response scanning and network-based mutational profiling approaches probed the effect of the Omicron variants on allosteric communications. The results of this analysis revealed specific roles of Omicron mutations as "plastic and evolutionary adaptable" modulators of binding and allostery which are coupled to the major regulatory positions through interaction networks. Through perturbation network scanning of allosteric residue potentials in the Omicron variant complexes, which is performed in the background of the original strain, we identified that the key Omicron binding affinity hotspots N501Y and Q498R could mediate allosteric interactions and epistatic couplings. Our results suggested that the synergistic role of these hotspots in controlling stability, binding and allostery can enable for compensatory balance of fitness tradeoffs with conformationally and evolutionary adaptable immune-escape Omicron mutations. Through integrative computational approaches, this study provides a systematic analysis of the effects of Omicron mutations on thermodynamics, binding and allosteric signaling in the complexes with ACE2 receptor. The findings support a mechanism in which Omicron mutations can evolve to balance thermodynamic stability and conformational adaptability in order to ensure proper tradeoff between stability, binding and immune escape.

Enhanced stability of the SARS CoV-2 spike glycoprotein following modification of an alanine cavity in the protein core

The spike (S) glycoprotein of SARS CoV-2 is the target of neutralizing antibodies (NAbs) that are crucial for vaccine effectiveness. The S1 subunit binds ACE2 while the S2 subunit mediates virus-cell membrane fusion. S2 is a class I fusion glycoprotein subunit and contains a central coiled coil that acts as a scaffold for the conformational changes associated with fusion function. The coiled coil of S2 is unusual in that the 3-4 repeat of inward-facing positions are mostly occupied by polar residues that mediate few inter-helical contacts in the prefusion trimer. We examined how insertion of bulkier hydrophobic residues (Val, Leu, Ile, Phe) to fill a cavity next to Ala1016 and Ala1020 in the 3-4 repeat affects the stability and antigenicity of S trimers. Substitution of Ala1016 with bulkier hydrophobic residues in the context of a prefusion-stabilized S trimer, S2P-FHA, was associated with increased thermal stability. S glycoprotein membrane fusion function was retained with Ala1016/Ala1020 cavity-filling mutations associated with improved recombinant S2P-FHA thermostability, however 2 mutants, A1016L and A1016V/A1020I, lacked ability to mediate entry of S-HIV-1 pseudoparticles into 293-ACE2 cells. When assessed as immunogens, two thermostable S2P-FHA mutants derived from the ancestral isolate, A1016L (16L) and A1016V/A1020I (VI) elicited neutralizing antibody with 50%-inhibitory dilutions (ID50s) in the range 2,700-5,110 for ancestral and Delta-derived viruses, and 210-1,744 for Omicron BA.1. The antigens elicited antibody specificities directed to the receptor-binding domain (RBD), N-terminal domain (NTD), fusion peptide and stem region of S2. The VI mutation enabled the production of intrinsically stable Omicron BA.1 and Omicron BA.4/5 S2P-FHA-like ectodomain oligomers in the absence of an external trimerization motif (T4 foldon), thus representing an alternative approach for stabilizing oligomeric S glycoprotein vaccines.

Coarse-Grained Molecular Simulations and Ensemble-Based Mutational Profiling of Protein Stability in the Different Functional Forms of the SARS-CoV-2 Spike Trimers: Balancing Stability and Adaptability in BA.1, BA.2 and BA.2.75 Variants

Evolutionary and functional studies have suggested that the emergence of Omicron variants can be determined by multiple fitness tradeoffs including immune escape, binding affinity, conformational plasticity, protein stability, and allosteric modulation. In this study, we embarked on a systematic comparative analysis of the conformational dynamics, electrostatics, protein stability, and allostery in the different functional states of spike trimers for BA.1, BA.2, and BA.2.75 variants. Using efficient and accurate coarse-grained simulations and atomistic reconstruction of the ensembles, we examined the conformational dynamics of the spike trimers that agree with the recent functional studies, suggesting that BA.2.75 trimers are the most stable among these variants. A systematic mutational scanning of the inter-protomer interfaces in the spike trimers revealed a group of conserved structural stability hotspots that play a key role in the modulation of functional dynamics and are also involved in the inter-protomer couplings through local contacts and interaction networks with the Omicron mutational sites. The results of mutational scanning provided evidence that BA.2.75 trimers are more stable than BA.2 and comparable in stability to the BA.1 variant. Using dynamic network modeling of the S Omicron BA.1, BA.2, and BA.2.75 trimers, we showed that the key network mediators of allosteric interactions are associated with the major stability hotspots that are interconnected along potential communication pathways. The network analysis of the BA.1, BA.2, and BA.2.75 trimers suggested that the increased thermodynamic stability of the BA.2.75 variant may be linked with the organization and modularity of the residue interaction network that allows for allosteric communications between structural stability hotspots and Omicron mutational sites. This study provided a plausible rationale for a mechanism in which Omicron mutations may evolve by targeting vulnerable sites of conformational adaptability to elicit immune escape while maintaining their control on balancing protein stability and functional fitness through robust allosteric communications with the stability hotspots.

Genomic surveillance of SARS-CoV-2 strains circulating in Iran during six waves of the pandemic

Background

SARS-CoV-2 genomic surveillance is necessary for the detection, monitoring, and evaluation of virus variants, which can have increased transmissibility, disease severity, or other adverse effects. We sequenced 330 SARS-CoV-2 genomes during the sixth wave of the COVID pandemic in Iran and compared them with five previous waves, for identifying SARS-CoV-2 variants, the genomic behavior of the virus, and understanding its characteristics.

Methods

After viral RNA extraction from clinical samples collected during the COVID-19 pandemic, next generation sequencing was performed using the Nextseq and Nanopore platforms. The sequencing data were analyzed and compared with reference sequences.

Results

In Iran during the first wave, V and L clades were detected. The second wave was recognized by G, GH, and GR clades. Circulating clades during the third wave were GH and GR. In the fourth wave, GRY (alpha variant), GK (delta variant), and one GH clade (beta variant) were detected. All viruses in the fifth wave were in GK clade (delta variant). In the sixth wave, Omicron variant (GRA clade) was circulating.

Conclusions

Genome sequencing, a key strategy in genomic surveillance systems, helps to detect and monitor the prevalence of SARS-CoV-2 variants, monitor the viral evolution of SARS-CoV-2, identify new variants for disease prevention, control, and treatment, and also provide information for and conduct public health measures in this area. With this system, Iran could be ready for surveillance of other respiratory virus diseases besides influenza and SARS-CoV-2.

Conformational Behavior of SARS-Cov-2 Spike Protein Variants: Evolutionary Jumps in Sequence Reverberate in Structural Dynamic Differences

SARS-CoV-2 has evolved rapidly in the first 3 years of pandemic diffusion. The initial evolution of the virus appeared to proceed through big jumps in sequence changes rather than through the stepwise accumulation of point mutations on already established variants. Here, we examine whether this nonlinear mutational process reverberates in variations of the conformational dynamics of the SARS-CoV-2 Spike protein (S-protein), the first point of contact between the virus and the human host. We run extensive microsecond-scale molecular dynamics simulations of seven distinct variants of the protein in their fully glycosylated state and set out to elucidate possible links between the mutational spectrum of the S-protein and the structural dynamics of the respective variant, at global and local levels. The results reveal that mutation-dependent structural and dynamic modulations mostly consist of increased coordinated motions in variants that acquire stability and in an increased internal flexibility in variants that are less stable. Importantly, a limited number of functionally important substructures (the receptor binding domain, in particular) share the same time of movements in all variants, indicating efficient preorganization for functional regions dedicated to host interactions. Our results support a model in which the internal dynamics of the S-proteins from different strains varies in a way that reflects the observed random and non-stepwise jumps in sequence evolution, while conserving the functionally oriented traits of conformational dynamics necessary to support productive interactions with host receptors.

Structural dynamics in the evolution of SARS-CoV-2 spike glycoprotein

SARS-CoV-2 spike glycoprotein mediates receptor binding and subsequent membrane fusion. It exists in a range of conformations, including a closed state unable to bind the ACE2 receptor, and an open state that does so but displays more exposed antigenic surface. Spikes of variants of concern (VOCs) acquired amino acid changes linked to increased virulence and immune evasion. Here, using HDX-MS, we identified changes in spike dynamics that we associate with the transition from closed to open conformations, to ACE2 binding, and to specific mutations in VOCs. We show that the RBD-associated subdomain plays a role in spike opening, whereas the NTD acts as a hotspot of conformational divergence of VOC spikes driving immune evasion. Alpha, beta and delta spikes assume predominantly open conformations and ACE2 binding increases the dynamics of their core helices, priming spikes for fusion. Conversely, substitutions in omicron spike lead to predominantly closed conformations, presumably enabling it to escape antibodies. At the same time, its core helices show characteristics of being pre-primed for fusion even in the absence of ACE2. These data inform on SARS-CoV-2 evolution and omicron variant emergence.

Molecular Characterization of AZD7442 (Tixagevimab-Cilgavimab) Neutralization of SARS-CoV-2 Omicron Subvariants

Therapeutic anti-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) monoclonal antibodies (MAbs) provide immunosuppressed and vulnerable populations with prophylactic and treatment interventions against coronavirus disease 2019 (COVID-19). AZD7442 (tixagevimab-cilgavimab) is a combination of extended-half-life neutralizing MAbs that bind to distinct epitopes on the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. The Omicron variant of concern carries mutations at >35 positions in the spike protein and has undergone further genetic diversification since its emergence in November 2021. Here, we characterize the in vitro neutralization activity of AZD7442 toward major viral subvariants circulating worldwide during the first 9 months of the Omicron wave. BA.2 and its derived subvariants showed the highest susceptibility to AZD7442, while BA.1 and BA.1.1 showed a lower susceptibility. BA.4/BA.5 had a susceptibility level intermediate between BA.1 and BA.2. Mutagenesis of parental Omicron subvariant spike proteins was performed to establish a molecular model to describe the underlying determinants of neutralization by AZD7442 and its component MAbs. The concurrent mutation of residues at positions 446 and 493, located in the tixagevimab and cilgavimab binding sites, was sufficient to enhance in vitro susceptibility of BA.1 to AZD7442 and its component MAbs to levels similar to the Wuhan-Hu-1+D614G virus. AZD7442 maintained neutralization activity against all Omicron subvariants tested up to and including BA.5. The evolving nature of the SARS-CoV-2 pandemic warrants continuing real-time molecular surveillance and assessment of in vitro activity of MAbs used in prophylaxis against and the treatment of COVID-19. IMPORTANCE MAbs are key therapeutic options for COVID-19 prophylaxis and treatment in immunosuppressed and vulnerable populations. Due to the emergence of SARS-CoV-2 variants, including Omicron, it is vital to ensure that neutralization is maintained for MAb-based interventions. We studied the in vitro neutralization of AZD7442 (tixagevimab-cilgavimab), a cocktail of two long-acting MAbs targeting the SARS-CoV-2 spike protein, toward Omicron subvariants circulating from November 2021 to July 2022. AZD7442 neutralized major Omicron subvariants up to and including BA.5. The mechanism of action responsible for the lower in vitro susceptibility of BA.1 to AZD7442 was investigated using in vitro mutagenesis and molecular modeling. A combination of mutations at two spike protein positions, namely, 446 and 493, was sufficient to enhance BA.1 susceptibility to AZD7442 to levels similar to the Wuhan-Hu-1+D614G ancestral virus. The evolving nature of the SARS-CoV-2 pandemic warrants continuing real-time global molecular surveillance and mechanistic studies of therapeutic MAbs for COVID-19.

Functional HIV-1/HCV cross-reactive antibodies isolated from a chronically co-infected donor

Despite prolific efforts to characterize the antibody response to human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV) mono-infections, the response to chronic co-infection with these two ever-evolving viruses is poorly understood. Here, we investigate the antibody repertoire of a chronically HIV-1/HCV co-infected individual using linking B cell receptor to antigen specificity through sequencing (LIBRA-seq). We identify five HIV-1/HCV cross-reactive antibodies demonstrating binding and functional cross-reactivity between HIV-1 and HCV envelope glycoproteins. All five antibodies show exceptional HCV neutralization breadth and effector functions against both HIV-1 and HCV. One antibody, mAb688, also cross-reacts with influenza and coronaviruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We examine the development of these antibodies using next-generation sequencing analysis and lineage tracing and find that somatic hypermutation established and enhanced this reactivity. These antibodies provide a potential future direction for therapeutic and vaccine development against current and emerging infectious diseases. More broadly, chronic co-infection represents a complex immunological challenge that can provide insights into the fundamental rules that underly antibody-antigen specificity.

Therapeutic and vaccine-induced cross-reactive antibodies with effector function against emerging Omicron variants

Currently circulating SARS-CoV-2 variants acquired convergent mutations at receptor-binding domain (RBD) hot spots. Their impact on viral infection, transmission, and efficacy of vaccines and therapeutics remains poorly understood. Here, we demonstrate that recently emerged BQ.1.1. and XBB.1 variants bind ACE2 with high affinity and promote membrane fusion more efficiently than earlier Omicron variants. Structures of the BQ.1.1 and XBB.1 RBDs bound to human ACE2 and S309 Fab (sotrovimab parent) explain the altered ACE2 recognition and preserved antibody binding through conformational selection. We show that sotrovimab binds avidly to all Omicron variants, promotes Fc-dependent effector functions and protects mice challenged with BQ.1.1, the variant displaying the greatest loss of neutralization. Moreover, in several donors vaccine-elicited plasma antibodies cross-react with and trigger effector functions against Omicron variants despite reduced neutralizing activity. Cross-reactive RBD-directed human memory B cells remained dominant even after two exposures to Omicron spikes, underscoring persistent immune imprinting. Our findings suggest that this previously overlooked class of cross-reactive antibodies, exemplified by S309, may contribute to protection against disease caused by emerging variants through elicitation of effector functions.

Omicron mutations increase interdomain interactions and reduce epitope exposure in the SARS-CoV-2 spike

Omicron BA.1 is a highly infectious variant of SARS-CoV-2 that carries more than thirty mutations on the spike protein in comparison to the Wuhan wild type (WT). Some of the Omicron mutations, located on the receptor-binding domain (RBD), are exposed to the surrounding solvent and are known to help evade immunity. However, the impact of buried mutations on the RBD conformations and on the mechanics of the spike opening is less evident. Here, we use all-atom molecular dynamics (MD) simulations with metadynamics to characterize the thermodynamic RBD-opening ensemble, identifying significant differences between WT and Omicron. Specifically, the Omicron mutations S371L, S373P, and S375F make more RBD interdomain contacts during the spike's opening. Moreover, Omicron takes longer to reach the transition state than WT. It stabilizes up-state conformations with fewer RBD epitopes exposed to the solvent, potentially favoring immune or antibody evasion.