Structural diversity of the SARS-CoV-2 Omicron spike

Conformational trajectory of the HIV-1 fusion peptide during CD4-induced envelope opening

The hydrophobic fusion peptide (FP), a critical component of the HIV-1 entry machinery, is located at the N terminus of the envelope (Env) gp41 subunit. The receptor-binding gp120 subunit of Env forms a heterodimer with gp41. The gp120/gp41 heterodimer assembles into a homotrimer, in which FP is accessible for antibody binding. Env conformational changes or "opening" that follow receptor binding result in FP relocating to a newly formed interprotomer pocket at the gp41-gp120 interface where it is sterically inaccessible to antibodies. The mechanistic steps connecting the entry-related transition of antibody accessible-to-inaccessible FP configurations remain unresolved. Here, using SOSIP-stabilized Env ectodomains, we visualize that the FP remains accessible for antibody binding despite substantial receptor-induced Env opening. We delineate stepwise Env opening from its closed state to a functional CD4-bound symmetrically open Env in which we show that FP was accessible for antibody binding. We define downstream re-organizations that lead to the formation of a gp120/gp41 cavity into which the FP buries to become inaccessible for antibody binding. These findings improve our understanding of HIV-1 entry and delineate the entry-related conformational trajectory of a key site of HIV vulnerability to neutralizing antibody.

A global collaboration for systematic analysis of broad-ranging antibodies against the SARS-CoV-2 spike protein

The Coronavirus Immunotherapeutic Consortium (CoVIC) conducted side-by-side comparisons of over 400 anti-SARS-CoV-2 spike therapeutic antibody candidates contributed by large and small companies as well as academic groups on multiple continents. Nine reference labs analyzed antibody features, including in vivo protection in a mouse model of infection, spike protein affinity, high-resolution epitope binning, ACE-2 binding blockage, structures, and neutralization of pseudovirus and authentic virus infection, to build a publicly accessible dataset in the database CoVIC-DB. High-throughput, high-resolution binning of CoVIC antibodies defines a broad and predictive landscape of antibody epitopes on the SARS-CoV-2 spike protein and identifies features associated with durable potency against multiple SARS-CoV-2 variants of concern and high in vivo efficacy. Results of the CoVIC studies provide a guide for selecting effective and durable antibody therapeutics and for immunogen design as well as providing a framework for rapid response to future viral disease outbreaks.

Structural and Energetic Insights into SARS-CoV-2 Evolution: Analysis of hACE2-RBD Binding in Wild-Type, Delta, and Omicron Subvariants

The evolution of SARS-CoV-2, particularly the emergence of Omicron variants, has raised questions regarding changes in its binding affinity to the human angiotensin-converting enzyme 2 receptor (hACE2). Understanding the impact of mutations on the interaction between the receptor-binding domain (RBD) of the spike protein and hACE2 is critical for evaluating viral transmissibility, immune evasion, and the efficacy of therapeutic strategies. Here, we used molecular dynamics (MD) simulations and binding energy calculations to investigate the structural and energetic differences between the hACE2- RBD complexes of wild-type (WT), Delta, and Omicron subvariants. Our results indicate that the Delta and the first Omicron variants showed the highest and the second-highest binding energy among the variants studied. Furthermore, while Omicron variants exhibit increased structural stability and altered electrostatic potential at the hACE2-RBD interface when compared to the ancestral WT, their binding strength to hACE2 does not consistently increase with viral evolution. Moreover, newer Omicron subvariants like JN.1 exhibit a bimodal conformational strategy, alternating between a high-affinity state for hACE2 and a low-affinity state, which could potentially facilitate immune evasion. These findings suggest that, in addition to enhanced hACE2 binding affinity, other factors, such as immune evasion and structural adaptability, shape SARS-CoV-2 evolution.

Exploring Diverse Binding Mechanisms of Broadly Neutralizing Antibodies S309, S304, CYFN-1006 and VIR-7229 Targeting SARS-CoV-2 Spike Omicron Variants: Integrative Computational Modeling Reveals Balance of Evolutionary and Dynamic Adaptability in Shaping Molecular Determinants of Immune Escape

Evolution of SARS-CoV-2 has led to the emergence of variants with increased immune evasion capabilities, posing significant challenges to antibody-based therapeutics and vaccines. The cross-neutralization activity of antibodies against Omicron variants is governed by a complex and delicate interplay of multiple energetic factors and interaction contributions. In this study, we conducted a comprehensive analysis of the interactions between the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein and four neutralizing antibodies S309, S304, CYFN1006, and VIR-7229. Using integrative computational modeling that combined all-atom molecular dynamics (MD) simulations, mutational scanning, and MM-GBSA binding free energy calculations, we elucidated the structural, energetic, and dynamic determinants of antibody binding. Our findings reveal distinct dynamic binding mechanisms and evolutionary adaptation driving broad neutralization effect of these antibodies. We show that S309 targets conserved residues near the ACE2 interface, leveraging synergistic van der Waals and electrostatic interactions, while S304 focuses on fewer but sensitive residues, making it more susceptible to escape mutations. The analysis of CYFN-1006.1 and CYFN-1006.2 antibody binding highlights broad epitope coverage with critical anchors at T345, K440, and T346, enhancing its efficacy against variants carrying the K356T mutation which caused escape from S309 binding. Our analysis of broadly potent VIR-7229 antibody binding to XBB.1.5 and EG.5 Omicron variants emphasized a large and structurally complex epitope, demonstrating certain adaptability and compensatory effects to F456L and L455S mutations. Mutational profiling identified key residues crucial for antibody binding, including T345, P337, and R346 for S309, and T385 and K386 for S304, underscoring their roles as evolutionary "weak spots" that balance viral fitness and immune evasion. The results of this energetic analysis demonstrate a good agreement between the predicted binding hotspots and critical mutations with respect to the latest experiments on average antibody escape scores. The results of this study dissect distinct energetic mechanisms of binding and importance of targeting conserved residues and diverse epitopes to counteract viral resistance. Broad-spectrum antibodies CYFN1006 and VIR-7229 maintain efficacy across multiple variants and achieve neutralization by targeting convergent evolution hotspots while enabling tolerance to mutations in these positions through structural adaptability and compensatory interactions at the binding interface. The results of this study underscore the diversity of binding mechanisms employed by different antibodies and molecular basis for high affinity and excellent neutralization activity of the latest generation of antibodies.

SARS-CoV-2 neutralizing antibody specificities differ dramatically between recently infected infants and immune-imprinted individuals

The immune response to viral infection is shaped by past exposures to related virus strains, a phenomenon known as imprinting. For severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), much of the population has been imprinted by a viral spike from an early strain, either through vaccination or infection during the early stages of the COVID-19 pandemic. As a consequence of this imprinting, infection with more recent SARS-CoV-2 strains primarily boosts cross-reactive antibodies elicited by the imprinting strain. Here we compare the neutralizing antibody specificities of imprinted individuals versus infants infected with a recent strain. Specifically, we use pseudovirus-based deep mutational scanning to measure how spike mutations affect neutralization by the serum antibodies of adults and children imprinted by the original vaccine versus infants with a primary infection by an XBB* variant. While the serum neutralizing activity of the imprinted individuals primarily targets the spike receptor-binding domain (RBD), the serum neutralizing activity of infants infected with only XBB* mostly targets the spike N-terminal domain. In these infants, secondary exposure to the XBB* spike via vaccination shifts more of the neutralizing activity toward the RBD, although the specific RBD sites targeted are different from imprinted adults. The dramatic differences in neutralization specificities among individuals with different exposure histories likely impact SARS-CoV-2 evolution.IMPORTANCEWe show that a person's exposure history to different SARS-CoV-2 strains strongly affects which regions on the viral spike that their neutralizing antibodies target. In particular, infants who have just been infected once with a recent viral strain make neutralizing antibodies that target different regions of the viral spike than adults or children who have been exposed to both older and more recent strains. This person-to-person heterogeneity means that the same viral mutation can have different impacts on the antibody immunity of different people.

Balancing stability and function: impact of the surface charge of SARS-CoV-2 Omicron spike protein

The ectodomain of the Omicron SARS-CoV-2 spike has an increased positive surface charge, favoring binding to the host cell surface, but may affect the stability of the ectodomain. Thermal stability studies identified two transitions associated with the flexibility of the receptor binding domain and the unfolding of the whole ectodomain, respectively. Despite destabilizing effects of some mutations, compensatory mutations maintain ECD stability and functional advantages thus supporting viral fitness.

Identification of Acanthopanax trifoliatus (L.) Merr as a Novel Potential Therapeutic Agent Against COVID-19 and Pharyngitis

Individuals infected with COVID-19 often experience the distressing discomfort of pharyngitis. Thus, it is crucial to develop novel drugs to improve therapeutic options. In this study, we investigated the interaction between bioactive compounds isolated from Acanthopanax trifoliatus (L.) Merr and proteins associated with COVID-19 and pharyngitis through in silico analysis. Several molecules demonstrated high affinities to multiple targets, indicating significant potential for alleviating pharyngitis and other COVID-19-related symptoms. Among them, rutin and isochlorogenic acid C, two major components in Acanthopanax trifoliatus (L.) Merr ethanol extracts, were further experimentally demonstrated to exhibit strong inhibitory effects against SARS-CoV-2 and to possess significant anti-inflammatory activities. Inhibition of over 50% in several key genes was observed, demonstrating the efficacy of in silico methods in identifying high-affinity target binders. Our findings provide a theoretical foundation for the development of Acanthopanax trifoliatus (L.) Merr as a novel multi-target therapeutic agent for both COVID-19 and pharyngitis.

Mutational Scanning and Binding Free Energy Computations of the SARS-CoV-2 Spike Complexes with Distinct Groups of Neutralizing Antibodies: Energetic Drivers of Convergent Evolution of Binding Affinity and Immune Escape Hotspots

The rapid evolution of SARS-CoV-2 has led to the emergence of variants with increased immune evasion capabilities, posing significant challenges to antibody-based therapeutics and vaccines. In this study, we conducted a comprehensive structural and energetic analysis of SARS-CoV-2 spike receptor-binding domain (RBD) complexes with neutralizing antibodies from four distinct groups (A-D), including group A LY-CoV016, group B AZD8895 and REGN10933, group C LY-CoV555, and group D antibodies AZD1061, REGN10987, and LY-CoV1404. Using coarse-grained simplified simulation models, rapid energy-based mutational scanning, and rigorous MM-GBSA binding free energy calculations, we elucidated the molecular mechanisms of antibody binding and escape mechanisms, identified key binding hotspots, and explored the evolutionary strategies employed by the virus to evade neutralization. The residue-based decomposition analysis revealed energetic mechanisms and thermodynamic factors underlying the effect of mutations on antibody binding. The results demonstrate excellent qualitative agreement between the predicted binding hotspots and the latest experiments on antibody escape. These findings provide valuable insights into the molecular determinants of antibody binding and viral escape, highlighting the importance of targeting conserved epitopes and leveraging combination therapies to mitigate the risk of immune evasion.

Quantitative Characterization and Prediction of the Binding Determinants and Immune Escape Hotspots for Groups of Broadly Neutralizing Antibodies Against Omicron Variants: Atomistic Modeling of the SARS-CoV-2 Spike Complexes with Antibodies

A growing body of experimental and computational studies suggests that the cross-neutralization antibody activity against Omicron variants may be driven by the balance and tradeoff between multiple energetic factors and interaction contributions of the evolving escape hotspots involved in antigenic drift and convergent evolution. However, the dynamic and energetic details quantifying the balance and contribution of these factors, particularly the balancing nature of specific interactions formed by antibodies with epitope residues, remain largely uncharacterized. In this study, we performed molecular dynamics simulations, an ensemble-based deep mutational scanning of SARS-CoV-2 spike residues, and binding free energy computations for two distinct groups of broadly neutralizing antibodies: the E1 group (BD55-3152, BD55-3546, and BD5-5840) and the F3 group (BD55-3372, BD55-4637, and BD55-5514). Using these approaches, we examined the energetic determinants by which broadly potent antibodies can largely evade immune resistance. Our analysis revealed the emergence of a small number of immune escape positions for E1 group antibodies that correspond to the R346 and K444 positions in which the strong van der Waals and interactions act synchronously, leading to the large binding contribution. According to our results, the E1 and F3 groups of Abs effectively exploit binding hotspot clusters of hydrophobic sites that are critical for spike functions along with the selective complementary targeting of positively charged sites that are important for ACE2 binding. Together with targeting conserved epitopes, these groups of antibodies can lead expand the breadth and resilience of neutralization to the antigenic shifts associated with viral evolution. The results of this study and the energetic analysis demonstrate excellent qualitative agreement between the predicted binding hotspots and critical mutations with respect to the latest experiments on average antibody escape scores. We argue that the E1 and F3 groups of antibodies targeting binding epitopes may leverage strong hydrophobic interactions with the binding epitope hotspots that are critical for the spike stability and ACE2 binding, while escape mutations tend to emerge in sites associated with synergistically strong hydrophobic and electrostatic interactions.

SARS-CoV-2 neutralizing antibody specificities differ dramatically between recently infected infants and immune-imprinted individuals

The immune response to viral infection is shaped by past exposures to related virus strains, a phenomenon known as imprinting. For SARS-CoV-2, much of the population has been imprinted by a viral spike from an early strain, either through vaccination or infection during the early stages of the COVID-19 pandemic. As a consequence of this imprinting, infection with more recent SARS-CoV-2 strains primarily boosts cross-reactive antibodies elicited by the imprinting strain. Here we compare the neutralizing antibody specificities of imprinted individuals versus infants infected with a recent strain. Specifically, we use pseudovirus-based deep mutational scanning to measure how spike mutations affect neutralization by the serum antibodies of adults and children imprinted by the original vaccine versus infants with a primary infection by a XBB* variant. While the serum neutralizing activity of the imprinted individuals primarily targets the spike receptor-binding domain (RBD), serum neutralizing activity of infants only infected with XBB* mostly targets the spike N-terminal domain (NTD). In these infants, secondary exposure to the XBB* spike via vaccination shifts more of the neutralizing activity towards the RBD, although the specific RBD sites targeted are different than for imprinted adults. The dramatic differences in neutralization specificities among individuals with different exposure histories likely impact SARS-CoV-2 evolution.

Design of SARS-CoV-2 RBD Immunogens to Focus Immune Responses Towards Conserved Coronavirus Epitopes

SARS-CoV-2 continues to evolve, with new variants emerging that evade pre-existing immunity and limit the efficacy of existing vaccines. One approach towards developing superior, variant-proof vaccines is to engineer immunogens that preferentially elicit antibodies with broad cross-reactivity against SARS-CoV-2 and its variants by targeting conserved epitopes on spike. The inner and outer faces of the Receptor Binding Domain (RBD) are two such conserved regions targeted by antibodies that recognize diverse human and animal coronaviruses. To promote the elicitation of such antibodies by vaccination, we engineered "resurfaced" RBD immunogens that contained mutations at exposed RBD residues outside the target epitopes. In the context of pre-existing immunity, these vaccine candidates aim to disfavor the elicitation of strain-specific antibodies against the immunodominant Receptor Binding Motif (RBM) while boosting the induction of inner and outer face antibodies. The engineered resurfaced RBD immunogens were stable, lacked binding to monoclonal antibodies with limited breadth, and maintained strong interactions with target broadly neutralizing antibodies. When used as vaccines, they limited humoral responses against the RBM as intended. Multimerization on nanoparticles further increased the immunogenicity of the resurfaced RBDs immunogens, thus supporting resurfacing as a promising immunogen design approach to rationally shift natural immune responses to develop more protective vaccines.

Phenotypic evolution of SARS-CoV-2 spike during the COVID-19 pandemic

SARS-CoV-2 variants are mainly defined by mutations in their spike. It is therefore critical to understand how the evolutionary trajectories of spike affect virus phenotypes. So far, it has been challenging to comprehensively compare the many spikes that emerged during the pandemic in a single experimental platform. Here we generated a panel of recombinant viruses carrying different spike proteins from 27 variants circulating between 2020 and 2024 in the same genomic background. We then assessed several of their phenotypic traits both in vitro and in vivo. We found distinct phenotypic trajectories of spike among and between variants circulating before and after the emergence of Omicron variants. Spike of post-Omicron variants maintained enhanced tropism for the nasal epithelium and large airways but displayed, over time, several phenotypic traits typical of the pre-Omicron variants. Hence, spike with phenotypic features of both pre- and post-Omicron variants may continue to emerge in the future.

Structural Immunology of SARS-CoV-2

The SARS-CoV-2 spike (S) protein has undergone significant evolution, enhancing both receptor binding and immune evasion. In this review, we summarize ongoing efforts to develop antibodies targeting various epitopes of the S protein, focusing on their neutralization potency, breadth, and escape mechanisms. Antibodies targeting the receptor-binding site (RBS) typically exhibit high neutralizing potency but are frequently evaded by mutations in SARS-CoV-2 variants. In contrast, antibodies targeting conserved regions, such as the S2 stem helix and fusion peptide, exhibit broader reactivity but generally lower neutralization potency. However, several broadly neutralizing antibodies have demonstrated exceptional efficacy against emerging variants, including the latest omicron subvariants, underscoring the potential of targeting vulnerable sites such as RBS-A and RBS-D/CR3022. We also highlight public classes of antibodies targeting different sites on the S protein. The vulnerable sites targeted by public antibodies present opportunities for germline-targeting vaccine strategies. Overall, developing escape-resistant, potent antibodies and broadly effective vaccines remains crucial for combating future variants. This review emphasizes the importance of identifying key epitopes and utilizing antibody affinity maturation to inform future therapeutic and vaccine design.

The evolution of vaccine hesitancy through the COVID-19 pandemic: A semi-structured interview study on booster and bivalent doses

We sought in-depth understanding on the evolution of factors influencing COVID-19 booster dose and bivalent vaccine hesitancy in a longitudinal semi-structured interview-based qualitative study. Serial interviews were conducted between July 25th and September 1st, 2022 (Phase I: univalent booster dose availability), and between November 21st, 2022 and January 11th, 2023 (Phase II: bivalent vaccine availability). Adults (≥18 years) in Canada who had received an initial primary series and had not received a COVID-19 booster dose were eligible for Phase I, and subsequently invited to participate in Phase II. Twenty-two of twenty-three (96%) participants completed interviews for both phases (45 interviews). Nearly half of participants identified as a woman (n = 11), the median age was 37 years (interquartile range: 32-48), and most participants were employed full-time (n = 12); no participant reported needing to vaccinate (with a primary series) for their workplace. No participant reported having received a COVID-19 booster dose at the time of their interview in Phase II. Three themes relating to the development of hesitancy toward continued vaccination against COVID-19 were identified: 1) effectiveness (frequency concerns; infection despite vaccination); 2) necessity (less threatening, low urgency, alternate protective measures); and 3) information (need for data, contradiction and confusion, lack of trust, decreased motivation). The data from interviews with individuals who had not received a COVID-19 booster dose or bivalent vaccine despite having received a primary series of COVID-19 vaccines highlights actionable targets to address vaccine hesitancy and improve public health literacy.

Structure and dynamics of the interaction of Delta and Omicron BA.1 SARS-CoV-2 variants with REGN10987 Fab reveal mechanism of antibody action

Study of mechanisms by which antibodies recognize different viral strains is necessary for the development of new drugs and vaccines to treat COVID-19 and other infections. Here, we report 2.5 Å cryo-EM structure of the SARS-CoV-2 Delta trimeric S-protein in complex with Fab of the recombinant analog of REGN10987 neutralizing antibody. S-protein adopts "two RBD-down and one RBD-up" conformation. Fab interacts with RBDs in both conformations, blocking the recognition of angiotensin converting enzyme-2. Three-dimensional variability analysis reveals high mobility of the RBD/Fab regions. Interaction of REGN10987 with Wuhan, Delta, Omicron BA.1, and mutated variants of RBDs is analyzed by microscale thermophoresis, molecular dynamics simulations, and ΔG calculations with umbrella sampling and one-dimensional potential of mean force. Variability in molecular dynamics trajectories results in a large scatter of calculated ΔG values, but Boltzmann weighting provides an acceptable correlation with experiment. REGN10987 evasion of the Omicron variant is found to be due to the additive effect of the N440K and G446S mutations located at the RBD/Fab binding interface with a small effect of Q498R mutation. Our study explains the influence of known-to-date SARS-CoV-2 RBD mutations on REGN10987 recognition and highlights the importance of dynamics data beyond the static structure of the RBD/Fab complex.

Quantitative Characterization and Prediction of the Binding Determinants and Immune Escape Hotspots for Groups of Broadly Neutralizing Antibodies Against Omicron Variants: Atomistic Modeling of the SARS-CoV-2 Spike Complexes with Antibodies

The growing body of experimental and computational studies suggested that the cross-neutralization antibody activity against Omicron variants may be driven by balance and tradeoff of multiple energetic factors and interaction contributions of the evolving escape hotspots involved in antigenic drift and convergent evolution. However, the dynamic and energetic details quantifying the balance and contribution of these factors, particularly the balancing nature of specific interactions formed by antibodies with the epitope residues remain scarcely characterized. In this study, we performed molecular dynamics simulations, ensemble-based deep mutational scanning of SARS-CoV-2 spike residues and binding free energy computations for two distinct groups of broadly neutralizing antibodies : E1 group (BD55-3152, BD55-3546 and BD5-5840) and F3 group (BD55-3372, BD55-4637 and BD55-5514). Using these approaches, we examine the energetic determinants by which broadly potent antibodies can largely evade immune resistance. Our analysis revealed the emergence of a small number of immune escape positions for E1 group antibodies that correspond to R346 and K444 positions in which the strong van der Waals and interactions act synchronously leading to the large binding contribution. According to our results, E1 and F3 groups of Abs effectively exploit binding hotspot clusters of hydrophobic sites critical for spike functions along with selective complementary targeting of positively charged sites that are important for ACE2 binding. Together with targeting conserved epitopes, these groups of antibodies can lead to the expanded neutralization breadth and resilience to antigenic shift associated with viral evolution. The results of this study and the energetic analysis demonstrate excellent qualitative agreement between the predicted binding hotspots and critical mutations with respect to the latest experiments on average antibody escape scores. We argue that E1 and F3 groups of antibodies targeting binding epitopes may leverage strong hydrophobic interactions with the binding epitope hotspots critical for the spike stability and ACE2 binding, while escape mutations tend to emerge in sites associated with synergistically strong hydrophobic and electrostatic interactions.

Interfacial subregions of SARS-CoV-2 spike RBD to hACE2 affect intermolecular affinity by their distinct roles played in association and dissociation kinetics

SARS-CoV-2's rapid global transmission depends on spike RBD's strong affinity to hACE2. In the context of binding hot spots well defined, the work investigated how interfacial subregions of SARS-CoV-2 spike RBD to hACE2 affect intermolecular affinity and their potential distinct roles involved in association and dissociation kinetics due to their local structural characteristics. Three spatially consecutive subregions of SARS-CoV-2 RBD were structurally partitioned across RBD's receptor binding motif (RBM). Their impacts on binding affinity and kinetics were differentiated through a comprehensive SPR measurement of hACE2 binding by chimeric swap mutants of respective subdomains from SARS-CoV-2 VOCs & phylogenetically close sarbecoviruses, and further compared with those of included single mutations across RBM and around the RBD core. The data supports that the intermediate interfacial subregion of RBD involving key residue at 417 is the rate-limiting effector of association kinetics and the subregion encompassing residues at 501/498/449 is the key binding energy contributor dictating dissociation kinetics, both of which relate to SARS-CoV-2's adaptive mutational evolution and host tropism closely. The kinetic data and structural analysis of local mutations' impact on spike RBD's binding and thermal stability provide a new perspective in evaluating SARS-CoV-2 evolution and other sarbecoviruses' evolvable binding to hACE2. The inherent binding mode offers direct clues of valid epitope in designing new antibodies that the coronavirus can't elude.

SARS-CoV-2 Omicron variations reveal mechanisms controlling cell entry dynamics and antibody neutralization

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) is adapting to continuous presence in humans. Transitions to endemic infection patterns are associated with changes in the spike (S) proteins that direct virus-cell entry. These changes generate antigenic drift and thereby allow virus maintenance in the face of prevalent human antiviral antibodies. These changes also fine tune virus-cell entry dynamics in ways that optimize transmission and infection into human cells. Focusing on the latter aspect, we evaluated the effects of several S protein substitutions on virus-cell membrane fusion, an essential final step in enveloped virus-cell entry. Membrane fusion is executed by integral-membrane "S2" domains, yet we found that substitutions in peripheral "S1" domains altered late-stage fusion dynamics, consistent with S1-S2 heterodimers cooperating throughout cell entry. A specific H655Y change in S1 stabilized a fusion-intermediate S protein conformation and thereby delayed membrane fusion. The H655Y change also sensitized viruses to neutralization by S2-targeting fusion-inhibitory peptides and stem-helix antibodies. The antibodies did not interfere with early fusion-activating steps; rather they targeted the latest stages of S2-directed membrane fusion in a novel neutralization mechanism. These findings demonstrate that single amino acid substitutions in the S proteins both reset viral entry-fusion kinetics and increase sensitivity to antibody neutralization. The results exemplify how selective forces driving SARS-CoV-2 fitness and antibody evasion operate together to shape SARS-CoV-2 evolution.

Conformational dynamics of SARS-CoV-2 Omicron spike trimers during fusion activation at single molecule resolution

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron entry involves spike (S) glycoprotein-mediated fusion of viral and late endosomal membranes. Here, using single-molecule Förster resonance energy transfer (sm-FRET) imaging and biochemical measurements, we directly visualized conformational changes of individual spike trimers on the surface of SARS-CoV-2 Omicron pseudovirions during fusion activation. We observed that the S2 domain of the Omicron spike is a dynamic fusion machine. S2 reversibly interchanges between the pre-fusion conformation and two previously undescribed intermediate conformations. Acidic pH shifts the conformational equilibrium of S2 toward an intermediate conformation and promotes the membrane hemi-fusion reaction. Moreover, we captured conformational reversibility in the S2 domain, which suggests that spike can protect itself from pre-triggering. Furthermore, we determined that Ca2+ directly promotes the S2 conformational change from an intermediate conformation to post-fusion conformation. In the presence of a target membrane, low pH and Ca2+ stimulate the irreversible transition to S2 post-fusion state and promote membrane fusion.

Molecular and structural insights into SARS-CoV-2 evolution: from BA.2 to XBB subvariants

Due to the incessant emergence of various SARS-CoV-2 variants with enhanced fitness in the human population, controlling the COVID-19 pandemic has been challenging. Understanding how the virus enhances its fitness during a pandemic could offer valuable insights for more effective control of viral epidemics. In this manuscript, we review the evolution of SARS-CoV-2 from early 2022 to the end of 2023-from Omicron BA.2 to XBB descendants. Focusing on viral evolution during this period, we provide concrete examples that SARS-CoV-2 has increased its fitness by enhancing several functions of the spike (S) protein, including its binding affinity to the ACE2 receptor and its ability to evade humoral immunity. Furthermore, we explore how specific mutations modify these functions of the S protein through structural alterations. This review provides evolutionary, molecular, and structural insights into how SARS-CoV-2 has increased its fitness and repeatedly caused epidemic surges during the pandemic.

Conformational trajectory of the HIV-1 fusion peptide during CD4-induced envelope opening

The hydrophobic fusion peptide (FP), a critical component of the HIV-1 entry machinery, is located at the N terminal stretch of the envelope (Env) gp41 subunit 1-3 . The receptor-binding gp120 subunit of Env forms a heterodimer with gp41 and assembles into a trimer, in which FP is accessible for antibody binding 3 . Env conformational changes or "opening" that follow receptor binding result in FP relocating to a newly formed interprotomer pocket at the gp41-gp120 interface where it is sterically inaccessible to antibody 4 . The mechanistic steps connecting the entry-related transition of antibody accessible-to-inaccessible FP configurations remain unresolved. Here, using SOSIP-stabilized Env ectodomains 5 , we visualized atomic-level details of a functional entry intermediate, where partially open Env was bound to receptor CD4, co-receptor mimetic antibody 17b, and FP-targeting antibody VRC34.01, demonstrating that FP remains antibody accessible despite substantial receptor-induced Env opening. We determined a series of structures delineating stepwise opening of Env from its closed state to a newly resolved intermediate and defining downstream re-organizations of the gp120-gp41 interface that ultimately resulted in FP burial in an antibody-inaccessible configuration. Our studies improve our understanding of HIV-1 entry and provide information on entry-related conformation reorganization of a key site of HIV vulnerability to neutralizing antibody.

Prying the lid open: Atomic-level insights on sialoglycan-TMPRSS2 coordination in HKU1 entry

The pre-fusion coronavirus HKU1 spike binds host sialoglycans and proteinaceous receptor TMPRSS2 for cell entry. In this issue of Cell, three papers by Fernández et al., McCallum et al., and Wang et al. provide structural information on HKU1 spike interactions with host receptors, providing insights into its multi-step opening.

The S2 subunit of spike encodes diverse targets for functional antibody responses to SARS-CoV-2

The SARS-CoV-2 virus responsible for the COVID-19 global pandemic has exhibited a striking capacity for viral evolution that drives continued evasion from vaccine and infection-induced immune responses. Mutations in the receptor binding domain of the S1 subunit of the spike glycoprotein have led to considerable escape from antibody responses, reducing the efficacy of vaccines and monoclonal antibody (mAb) therapies. Therefore, there is a need to interrogate more constrained regions of spike, such as the S2 subdomain. Here, we present a collection of S2 mAbs from two SARS-CoV-2 convalescent individuals that target multiple regions in S2, including regions outside of those commonly reported. One of the S2 mAbs, C20.119, which bound to a highly conserved epitope in the fusion peptide, was able to broadly neutralize across SARS-CoV-2 variants, SARS-CoV-1, and closely related zoonotic sarbecoviruses. The majority of the mAbs were non-neutralizing; however, many of them could mediate antibody-dependent cellular cytotoxicity (ADCC) at levels similar to the S1-targeting mAb S309 that was previously authorized for treatment of SARS-CoV-2 infections. Several of the mAbs with ADCC function also bound to spike trimers from other human coronaviruses (HCoVs), such as MERS-CoV and HCoV-HKU1. Our findings suggest S2 mAbs can target diverse epitopes in S2, including functional mAbs with HCoV and sarbecovirus breadth that likely target functionally constrained regions of spike. These mAbs could be developed for potential future pandemics, while also providing insight into ideal epitopes for eliciting a broad HCoV response.

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.

Cryo-electron microscopy in the study of virus entry and infection

Viruses have been responsible for many epidemics and pandemics that have impacted human life globally. The COVID-19 pandemic highlighted both our vulnerability to viral outbreaks, as well as the mobilization of the scientific community to come together to combat the unprecedented threat to humanity. Cryo-electron microscopy (cryo-EM) played a central role in our understanding of SARS-CoV-2 during the pandemic and continues to inform about this evolving pathogen. Cryo-EM with its two popular imaging modalities, single particle analysis (SPA) and cryo-electron tomography (cryo-ET), has contributed immensely to understanding the structure of viruses and interactions that define their life cycles and pathogenicity. Here, we review how cryo-EM has informed our understanding of three distinct viruses, of which two - HIV-1 and SARS-CoV-2 infect humans, and the third, bacteriophages, infect bacteria. For HIV-1 and SARS-CoV-2 our focus is on the surface glycoproteins that are responsible for mediating host receptor binding, and host and cell membrane fusion, while for bacteriophages, we review their structure, capsid maturation, attachment to the bacterial cell surface and infection initiation mechanism.

Assessing pH-Dependent Conformational Changes in the Fusion Peptide Proximal Region of the SARS-CoV-2 Spike Glycoprotein

One of the entry mechanisms of the SARS-CoV-2 coronavirus into host cells involves endosomal acidification. It has been proposed that under acidic conditions, the fusion peptide proximal region (FPPR) of the SARS-CoV-2 spike glycoprotein acts as a pH-dependent switch, modulating immune response accessibility by influencing the positioning of the receptor binding domain (RBD). This would provide indirect coupling of RBD opening to the environmental pH. Here, we explored this possible pH-dependent conformational equilibrium of the FPPR within the SARS-CoV-2 spike glycoprotein. We analyzed hundreds of experimentally determined spike structures from the Protein Data Bank and carried out pH-replica exchange molecular dynamics to explore the extent to which the FPPR conformation depends on pH and the positioning of the RBD. A meta-analysis of experimental structures identified alternate conformations of the FPPR among structures in which this flexible regions was resolved. However, the results did not support a correlation between the FPPR conformation and either RBD position or the reported pH of the cryo-EM experiment. We calculated pKa values for titratable side chains in the FPPR region using PDB structures, but these pKa values showed large differences between alternate PDB structures that otherwise adopt the same FPPR conformation type. This hampers the comparison of pKa values in different FPPR conformations to rationalize a pH-dependent conformational change. We supplemented these PDB-based analyses with all-atom simulations and used constant-pH replica exchange molecular dynamics to estimate pKa values in the context of flexibility and explicit water. The resulting titration curves show good reproducibility between simulations, but they also suggest that the titration curves of the different FPPR conformations are the same within the error bars. In summary, we were unable to find evidence supporting the previously published hypothesis of an FPPR pH-dependent equilibrium: neither from existing experimental data nor from constant-pH MD simulations. The study underscores the complexity of the spike system and opens avenues for further exploration into the interplay between pH and SARS-CoV-2 viral entry mechanisms.

Micro-second time-resolved X-ray single-molecule internal motions of SARS-CoV-2 spike variants

Single-molecule intramolecular dynamics were successfully measured for three variants of SARS-CoV-2 spike protein, alpha: B.1.1.7, delta: B.1.617, and omicron: B.1.1.529, with a time resolution of 100 μs using X-rays. The results were then compared with respect to the magnitude and directions of motions for the three variants. The largest 3-D intramolecular movement was observed for the omicron variant irrespective of ACE2 receptor binding. A more detailed analysis of the intramolecular motions revealed that the distribution state of intramolecular motion for the three variants was completely different with and without ACE2 receptor binding. The molecular dynamics for the trimeric spike protein of the omicron variant increased when ACE2 binding occurred. At that time, the diffusion constant increased from 71.0 [mrad2/ms] to 91.1 [mrad2/ms].

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.

APOBEC3-related mutations in the spike protein-encoding region facilitate SARS-CoV-2 evolution

SARS-CoV-2 evolves gradually to cause COVID-19 epidemic. One of driving forces of SARS-CoV-2 evolution might be activation of apolipoprotein B mRNA editing catalytic subunit-like protein 3 (APOBEC3) by inflammatory factors. Here, we aimed to elucidate the effect of the APOBEC3-related viral mutations on the infectivity and immune evasion of SARS-CoV-2. The APOBEC3-related C > U mutations ranked as the second most common mutation types in the SARS-CoV-2 genome. mRNA expression of APOBEC3A (A3A), APOBEC3B (A3B), and APOBEC3G (A3G) in peripheral blood cells increased with disease severity. A3B, a critical member of the APOBEC3 family, was significantly upregulated in both severe and moderate COVID-19 patients and positively associated with neutrophil proportion and COVID-19 severity. We identified USP18 protein, a key molecule centralizing the protein-protein interaction network of key APOBEC3 proteins. Furthermore, mRNA expression of USP18 was significantly correlated to ACE2 and TMPRSS2 expression in the tissue of upper airways. Knockdown of USP18 mRNA significantly decreased A3B expression. Ectopic expression of A3B gene increased SARS-CoV-2 infectivity. C > U mutations at S371F, S373L, and S375F significantly conferred with the immune escape of SARS-CoV-2. Thus, APOBEC3, whose expression are upregulated by inflammatory factors, might promote SARS-CoV-2 evolution and spread via upregulating USP18 level and facilitating the immune escape. A3B and USP18 might be therapeutic targets for interfering with SARS-CoV-2 evolution.

In Silico Analyses Indicate a Lower Potency for Dimerization of TLR4/MD-2 as the Reason for the Lower Pathogenicity of Omicron Compared to Wild-Type Virus and Earlier SARS-CoV-2 Variants

The SARS-CoV-2 Omicron variants have replaced all earlier variants, due to increased infectivity and effective evasion from infection- and vaccination-induced neutralizing antibodies. Compared to earlier variants of concern (VoCs), the Omicron variants show high TMPRSS2-independent replication in the upper airway organs, but lower replication in the lungs and lower mortality rates. The shift in cellular tropism and towards lower pathogenicity of Omicron was hypothesized to correlate with a lower toll-like receptor (TLR) activation, although the underlying molecular mechanisms remained undefined. In silico analyses presented here indicate that the Omicron spike protein has a lower potency to induce dimerization of TLR4/MD-2 compared to wild type virus despite a comparable binding activity to TLR4. A model illustrating the molecular consequences of the different potencies of the Omicron spike protein vs. wild-type spike protein for TLR4 activation is presented. Further analyses indicate a clear tendency for decreasing TLR4 dimerization potential during SARS-CoV-2 evolution via Alpha to Gamma to Delta to Omicron variants.

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.

A chimeric adenovirus-vectored vaccine based on Beta spike and Delta RBD confers a broad-spectrum neutralization against Omicron-included SARS-CoV-2 variants

Urgent research into innovative severe acute respiratory coronavirus-2 (SARS-CoV-2) vaccines that may successfully prevent various emerging emerged variants, particularly the Omicron variant and its subvariants, is necessary. Here, we designed a chimeric adenovirus-vectored vaccine named Ad5-Beta/Delta. This vaccine was created by incorporating the receptor-binding domain from the Delta variant, which has the L452R and T478K mutations, into the complete spike protein of the Beta variant. Both intramuscular (IM) and intranasal (IN) vaccination with Ad5-Beta/Deta vaccine induced robust broad-spectrum neutralization against Omicron BA.5-included variants. IN immunization with Ad5-Beta/Delta vaccine exhibited superior mucosal immunity, manifested by higher secretory IgA antibodies and more tissue-resident memory T cells (TRM) in respiratory tract. The combination of IM and IN delivery of the Ad5-Beta/Delta vaccine was capable of synergically eliciting stronger systemic and mucosal immune responses. Furthermore, the Ad5-Beta/Delta vaccination demonstrated more effective boosting implications after two dosages of mRNA or subunit recombinant protein vaccine, indicating its capacity for utilization as a booster shot in the heterologous vaccination. These outcomes quantified Ad5-Beta/Delta vaccine as a favorable vaccine can provide protective immunity versus SARS-CoV-2 pre-Omicron variants of concern and BA.5-included Omicron subvariants.

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.

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

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

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.

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

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

Fusing Peptide Epitopes for Advanced Multiplex Serological Testing for SARS-CoV-2 Antibody Detection

The tragic COVID-19 pandemic, which has seen a total of 655 million cases worldwide and a death toll of over 6.6 million seems finally tailing off. Even so, new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continue to arise, the severity of which cannot be predicted in advance. This is concerning for the maintenance and stability of public health, since immune evasion and increased transmissibility may arise. Therefore, it is crucial to continue monitoring antibody responses to SARS-CoV-2 in the general population. As a complement to polymerase chain reaction tests, multiplex immunoassays are elegant tools that use individual protein or peptide antigens simultaneously to provide a high level of sensitivity and specificity. To further improve these aspects of SARS-CoV-2 antibody detection, as well as accuracy, we have developed an advanced serological peptide-based multiplex assay using antigen-fused peptide epitopes derived from both the spike and the nucleocapsid proteins. The significance of the epitopes selected for antibody detection has been verified by in silico molecular docking simulations between the peptide epitopes and reported SARS-CoV-2 antibodies. Peptides can be more easily and quickly modified and synthesized than full length proteins and can, therefore, be used in a more cost-effective manner. Three different fusion-epitope peptides (FEPs) were synthesized and tested by enzyme-linked immunosorbent assay (ELISA). A total of 145 blood serum samples were used, compromising 110 COVID-19 serum samples from COVID-19 patients and 35 negative control serum samples taken from COVID-19-free individuals before the outbreak. Interestingly, our data demonstrate that the sensitivity, specificity, and accuracy of the results for the FEP antigens are higher than for single peptide epitopes or mixtures of single peptide epitopes. Our FEP concept can be applied to different multiplex immunoassays testing not only for SARS-CoV-2 but also for various other pathogens. A significantly improved peptide-based serological assay may support the development of commercial point-of-care tests, such as lateral-flow-assays.

Oral administration of a recombinant modified RBD antigen of SARS-CoV-2 as a possible immunostimulant for the care of COVID-19

BACKGROUND

Developing effective vaccines against SARS-CoV-2 that consider manufacturing limitations, equitable access, and acceptance is necessary for developing platforms to produce antigens that can be efficiently presented for generating neutralizing antibodies and as a model for new vaccines.

RESULTS

This work presents the development of an applicable technology through the oral administration of the SARS-CoV-2 RBD antigen fused with a peptide to improve its antigenic presentation. We focused on the development and production of the recombinant receptor binding domain (RBD) produced in E. coli modified with the addition of amino acids extension designed to improve antigen presentation. The production was carried out in shake flask and bioreactor cultures, obtaining around 200 mg/L of the antigen. The peptide-fused RBD and peptide-free RBD proteins were characterized and compared using SDS-PAGE gel, high-performance chromatography, and circular dichroism. The peptide-fused RBD was formulated in an oil-in-water emulsion for oral mice immunization. The peptide-fused RBD, compared to RBD, induced robust IgG production in mice, capable of recognizing the recombinant RBD in Enzyme-linked immunosorbent assays. In addition, the peptide-fused RBD generated neutralizing antibodies in the sera of the dosed mice. The formulation showed no reactive episodes and no changes in temperature or vomiting.

CONCLUSIONS

Our study demonstrated the effectiveness of the designed peptide added to the RBD to improve antigen immunostimulation by oral administration.

Recurrent SARS-CoV-2 mutations at Spike D796 evade antibodies from pre-Omicron convalescent and vaccinated subjects

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages of the Omicron variant rapidly became dominant in early 2022 and frequently cause human infections despite vaccination or prior infection with other variants. In addition to antibody-evading mutations in the receptor-binding domain, Omicron features amino acid mutations elsewhere in the Spike protein; however, their effects generally remain ill defined. The Spike D796Y substitution is present in all Omicron sub-variants and occurs at the same site as a mutation (D796H) selected during viral evolution in a chronically infected patient. Here, we map antibody reactivity to a linear epitope in the Spike protein overlapping position 796. We show that antibodies binding this region arise in pre-Omicron SARS-CoV-2 convalescent and vaccinated subjects but that both D796Y and D796H abrogate their binding. These results suggest that D796Y contributes to the fitness of Omicron in hosts with pre-existing immunity to other variants of SARS-CoV-2 by evading antibodies targeting this site.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved substantially through the coronavirus disease 2019 (COVID-19) pandemic: understanding the drivers and consequences of this evolution is essential for projecting the course of the pandemic and developing new countermeasures. Here, we study the immunological effects of a particular mutation present in the Spike protein of all Omicron strains and find that it prevents the efficient binding of a class of antibodies raised by pre-Omicron vaccination and infection. These findings reveal a novel consequence of a poorly understood Omicron mutation and shed light on the drivers and effects of SARS-CoV-2 evolution.

Combined in vitro/in silico approaches, molecular dynamics simulations and safety assessment of the multifunctional properties of thymol and carvacrol: A comparative insight

Bioactive compounds derived from medicinal plants have acquired immense attentiveness in drug discovery and development. The present study investigated in vitro and predicted in silico the antibacterial, antifungal, and antiviral properties of thymol and carvacrol, and assessed their safety. The performed microbiological assays against Pseudomonas aeruginosa, Escherichia coli, Salmonella enterica Typhimurium revealed that the minimal inhibitory concentration values ranged from (0.078 to 0.312 mg/mL) and the minimal fungicidal concentration against Candida albicans was 0.625 mg/mL. Molecular docking simulations, stipulated that these compounds could inhibit bacterial replication and transcription functions by targeting DNA and RNA polymerases receptors with docking scores varying between (-5.1 to -6.9 kcal/mol). Studied hydroxylated monoterpenes could hinder C. albicans growth by impeding lanosterol 14α-demethylase enzyme and showed a (ΔG=-6.2 and -6.3 kcal/mol). Computational studies revealed that thymol and carvacrol could target the SARS-Cov-2 spike protein of the Omicron variant RBD domain. Molecular dynamics simulations disclosed that these compounds have a stable dynamic behavior over 100 ns as compared to remdesivir. Chemo-computational toxicity prediction using Protox II webserver indicated that thymol and carvacrol could be safely and effectively used as drug candidates to tackle bacterial, fungal, and viral infections as compared to chemical medication.

Caveats of chimpanzee ChAdOx1 adenovirus-vectored vaccines to boost anti-SARS-CoV-2 protective immunity in mice

Several COVID-19 vaccines use adenovirus vectors to deliver the SARS-CoV-2 spike (S) protein. Immunization with these vaccines promotes immunity against the S protein, but against also the adenovirus itself. This could interfere with the entry of the vaccine into the cell, reducing its efficacy. Herein, we evaluate the efficiency of an adenovirus-vectored vaccine (chimpanzee ChAdOx1 adenovirus, AZD1222) in boosting the specific immunity compared to that induced by a recombinant receptor-binding domain (RBD)-based vaccine without viral vector. Mice immunized with the AZD1222 human vaccine were given a booster 6 months later, with either the homologous vaccine or a recombinant vaccine based on RBD of the delta variant, which was prevalent at the start of this study. A significant increase in anti-RBD antibody levels was observed in rRBD-boosted mice (31-61%) compared to those receiving two doses of AZD1222 (0%). Significantly higher rates of PepMix™- or RBD-elicited proliferation were also observed in IFNγ-producing CD4 and CD8 cells from mice boosted with one or two doses of RBD, respectively. The lower efficiency of the ChAdOx1-S vaccine in boosting specific immunity could be the result of a pre-existing anti-vector immunity, induced by increased levels of anti-adenovirus antibodies found both in mice and humans. Taken together, these results point to the importance of avoiding the recurrent use of the same adenovirus vector in individuals with immunity and memory against them. It also illustrates the disadvantages of ChAdOx1 adenovirus-vectored vaccine with respect to recombinant protein vaccines, which can be used without restriction in vaccine-booster programs. KEY POINTS: • ChAdOx1 adenovirus vaccine (AZD1222) may not be effective in boosting anti-SARS-CoV-2 immunity • A recombinant RBD protein vaccine is effective in boosting anti-SARS-CoV-2 immunity in mice • Antibodies elicited by the rRBD-delta vaccine persisted for up to 3 months in mice.

Novel Omicron Variants Enhance Anchored Recognition of TMEM106B: A New Pathway for SARS-CoV-2 Cellular Invasion

The recent discovery that TMEM106B serves as a receptor mediating ACE2-independent SARS-CoV-2 entry into cells deserves attention, especially in the background of the frequent emergence of mutant strains. Here, the structure-dynamic features of this novel pathway are dissected deeply. Our investigation revealed that the large loop (RBD@471-491) could anchor TMEM106B, which was then firmly locked by another loop (RBD@444-451). The novel and widely disseminated Omicron variants (BA.2.86/EG.5.1) affect the anchoring recognition of proteins, with BA.2.86 being more likely to impact cells with limited or undetectable ACE2 expression. The large loop of the EG.5.1 variant captures TMEM106B poorly due to impaired electrostatic complementarity. Furthermore, we emphasize that antibody design against these two loops could enhance the protection of ACE2 low-expressing cells according to the alanine scanning mutagenesis of multiple antibodies. We hope this study will provide a novel perspective for the prevention and treatment against this new viral invasion pathway.

Structural and biochemical rationale for Beta variant protein booster vaccine broad cross-neutralization of SARS-CoV-2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), responsible for the COVID-19 pandemic, uses a surface expressed trimeric spike glycoprotein for cell entry. This trimer is the primary target for neutralizing antibodies making it a key candidate for vaccine development. During the global pandemic circulating variants of concern (VOC) caused several waves of infection, severe disease, and death. The reduced efficacy of the ancestral trimer-based vaccines against emerging VOC led to the need for booster vaccines. Here we present a detailed characterization of the Sanofi Beta trimer, utilizing cryo-EM for structural elucidation. We investigate the conformational dynamics and stabilizing features using orthogonal SPR, SEC, nanoDSF, and HDX-MS techniques to better understand how this antigen elicits superior broad neutralizing antibodies as a variant booster vaccine. This structural analysis confirms the Beta trimer preference for canonical quaternary structure with two RBD in the up position and the reversible equilibrium between the canonical spike and open trimer conformations. Moreover, this report provides a better understanding of structural differences between spike antigens contributing to differential vaccine efficacy.

Highly potent dual-targeting angiotensin-converting enzyme 2 (ACE2) and Neuropilin-1 (NRP1) peptides: A promising broad-spectrum therapeutic strategy against SARS-CoV-2 infection

The efficacy of approved vaccines has been diminishing due to the increasing advent of SARS-CoV-2 variants with diverse mutations that favor sneak entry. Nonetheless, these variants recognize the conservative host receptors angiotensin-converting enzyme 2 (ACE2) and neuropilin-1 (NRP1) for entry, rendering the dual blockade of ACE2 and NRP1 an advantageous pan-inhibition strategy. Here, we identified a highly potent dual-targeting peptide AP-1 using structure-based virtual screening protocol. AP-1 had nanoscale binding affinities for ACE2 (Kd = 6.1 ± 0.2 nM) and NRP1 (Kd = 13.4 ± 1.2 nM) and approximately 102- and 8-fold stronger than positive inhibitors S471-503 and NMTP-5, respectively. Further evidence in pseudovirus cell infection and cytotoxicity assays demonstrated that AP-1 exhibited remarkable entry inhibition of variants of concern (VOCs) of SARS-CoV-2 without impairing host cell viability. Together, our findings suggest that AP-1 with dual-targeting ACE2/NRP1 efficacy could be a promising broad-spectrum agent for treating SARS-CoV-2 emerging VOCs.

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.

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.

Mutations in S2 subunit of SARS-CoV-2 Omicron spike strongly influence its conformation, fusogenicity, and neutralization sensitivity

IMPORTANCE

The Omicron subvariants have substantially evaded host-neutralizing antibodies and adopted an endosomal route of entry. The virus has acquired several mutations in the receptor binding domain and N-terminal domain of S1 subunit, but remarkably, also incorporated mutations in S2 which are fixed in Omicron sub-lineage. Here, we found that the mutations in the S2 subunit affect the structural and biological properties such as neutralization escape, entry route, fusogenicity, and protease requirement. In vivo, these mutations may have significant roles in tropism and replication. A detailed understanding of the effects of S2 mutations on Spike function, immune evasion, and viral entry would inform the vaccine design, as well as therapeutic interventions aiming to block the essential proteases for virus entry. Thus, our study has identified the crucial role of S2 mutations in stabilizing the Omicron spike and modulating neutralization resistance to antibodies targeting the S1 subunit.

Impact of mutations defining SARS-CoV-2 Omicron subvariants BA.2.12.1 and BA.4/5 on Spike function and neutralization

Additional mutations in the viral Spike protein helped the BA.2.12.1 and BA.4/5 SARS-CoV-2 Omicron subvariants to outcompete the parental BA.2 subvariant. Here, we determined the functional impact of mutations that newly emerged in the BA.2.12.1 (L452Q, S704L) and BA.4/5 (Δ69-70, L452R, F486V, R493Q) Spike proteins. Our results show that mutation of L452Q/R or F486V typically increases and R493Q or S704L impair BA.2 Spike-mediated infection. In combination, changes of Δ69-70, L452R, and F486V contribute to the higher infectiousness and fusogenicity of the BA.4/5 Spike. L452R/Q and F486V in Spike are mainly responsible for reduced sensitivity to neutralizing antibodies. However, the combined mutations are required for full infectivity, reduced TMPRSS2 dependency, and immune escape of BA.4/5 Spike. Thus, it is the specific combination of mutations in BA.4/5 Spike that allows increased functionality and immune evasion, which helps to explain the temporary dominance and increased pathogenicity of these Omicron subvariants.