Structural diversity of the SARS-CoV-2 Omicron spike

Multimodal, broadly neutralizing antibodies against SARS-CoV-2 identified by high-throughput native pairing of BCRs from bulk B cells

TruAB Discovery is an approach that integrates cellular immunology, high-throughput immunosequencing, bioinformatics, and computational biology in order to discover naturally occurring human antibodies for prophylactic or therapeutic use. We adapted our previously described pairSEQ technology to pair B cell receptor heavy and light chains of SARS-CoV-2 spike protein-binding antibodies derived from enriched antigen-specific memory B cells and bulk antibody-secreting cells. We identified approximately 60,000 productive, in-frame, paired antibody sequences, from which 2,093 antibodies were selected for functional evaluation based on abundance, isotype and patterns of somatic hypermutation. The exceptionally diverse antibodies included RBD-binders with broad neutralizing activity against SARS-CoV-2 variants, and S2-binders with broad specificity against betacoronaviruses and the ability to block membrane fusion. A subset of these RBD- and S2-binding antibodies demonstrated robust protection against challenge in hamster and mouse models. This high-throughput approach can accelerate discovery of diverse, multifunctional antibodies against any target of interest.

Population immunity predicts evolutionary trajectories of SARS-CoV-2

The large-scale evolution of the SARS-CoV-2 virus has been marked by rapid turnover of genetic clades. New variants show intrinsic changes, notably increased transmissibility, and antigenic changes that reduce cross-immunity induced by previous infections or vaccinations. How this functional variation shapes global evolution has remained unclear. Here, we establish a predictive fitness model for SARS-CoV-2 that integrates antigenic and intrinsic selection. The model is informed by tracking of time-resolved sequence data, epidemiological records, and cross-neutralization data of viral variants. Our inference shows that immune pressure, including contributions of vaccinations and previous infections, has become the dominant force driving the recent evolution of SARS-CoV-2. The fitness model can serve continued surveillance in two ways. First, it successfully predicts the short-term evolution of circulating strains and flags emerging variants likely to displace the previously predominant variant. Second, it predicts likely antigenic profiles of successful escape variants prior to their emergence.

Structural understanding of SARS-CoV-2 virus entry to host cells

Coronavirus disease 2019 (COVID-19), caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a major global health concern associated with millions of fatalities worldwide. Mutant variants of the virus have further exacerbated COVID-19 mortality and infection rates, emphasizing the urgent need for effective preventive strategies. Understanding the viral infection mechanism is crucial for developing therapeutics and vaccines. The entry of SARS-CoV-2 into host cells is a key step in the infection pathway and has been targeted for drug development. Despite numerous reviews of COVID-19 and the virus, there is a lack of comprehensive reviews focusing on the structural aspects of viral entry. In this review, we analyze structural changes in Spike proteins during the entry process, dividing the entry process into prebinding, receptor binding, proteolytic cleavage, and membrane fusion steps. By understanding the atomic-scale details of viral entry, we can better target the entry step for intervention strategies. We also examine the impacts of mutations in Spike proteins, including the Omicron variant, on viral entry. Structural information provides insights into the effects of mutations and can guide the development of therapeutics and vaccines. Finally, we discuss available structure-based approaches for the development of therapeutics and vaccines. Overall, this review provides a detailed analysis of the structural aspects of SARS-CoV-2 viral entry, highlighting its significance in the development of therapeutics and vaccines against COVID-19. Therefore, our review emphasizes the importance of structural information in combating SARS-CoV-2 infection.

Effects of N-glycan modifications on spike expression, virus infectivity, and neutralization sensitivity in ancestral compared to Omicron SARS-CoV-2 variants

The SARS-CoV-2 spike glycoprotein has 22 potential N-linked glycosylation sites per monomer that are highly conserved among diverse variants, but how individual glycans affect virus entry and neutralization of Omicron variants has not been extensively characterized. Here we compared the effects of specific glycan deletions or modifications in the Omicron BA.1 and D614G spikes on spike expression, processing, and incorporation into pseudoviruses, as well as on virus infectivity and neutralization by therapeutic antibodies. We found that loss of potential glycans at spike residues N717 and N801 each conferred a loss of pseudovirus infectivity for Omicron but not for D614G or Delta variants. This decrease in infectivity correlated with decreased spike processing and incorporation into Omicron pseudoviruses. Oligomannose-enriched Omicron pseudoviruses generated in GnTI- cells or in the presence of kifunensine were non-infectious, whereas D614G or Delta pseudoviruses generated under similar conditions remained infectious. Similarly, growth of live (authentic) SARS-CoV-2 in the presence of kifunensine resulted in a greater reduction of titers for the BA.1.1 variant than Delta or D614G variants relative to their respective, untreated controls. Finally, we found that loss of some N-glycans, including N343 and N234, increased the maximum percent neutralization by the class 3 S309 monoclonal antibody against D614G but not BA.1 variants, while these glycan deletions altered the neutralization potency of the class 1 COV2-2196 and Etesevimab monoclonal antibodies without affecting maximum percent neutralization. The maximum neutralization by some antibodies also varied with the glycan composition, with oligomannose-enriched pseudoviruses conferring the highest percent neutralization. These results highlight differences in the interactions between glycans and residues among SARS-CoV-2 variants that can affect spike expression, virus infectivity, and susceptibility of variants to antibody neutralization.

Molecular dynamics simulations of the spike trimeric ectodomain of the SARS-CoV-2 Omicron variant: structural relationships with infectivity, evasion to immune system and transmissibility

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron is currently the most prevalent SARS-CoV-2 variant worldwide. Herein, we calculated molecular dynamics simulations of the trimeric spikeWT and SpikeBA.1 for 300 ns. Our results show that SpikeBA.1 has more conformational flexibility than SpikeWT. Our principal component analysis (PCA) allowed us to observe a broader spectrum of different conformations for SpikeBA.1, mainly at N-terminal domain (NTD) and receptor-binding domain (RBD). Such increased flexibility could contribute to decreased neutralizing antibody recognition of this variant. Our molecular dynamics data show that the RBDBA.1 easily visits an up-conformational state and the prevalent D614G mutation is pivotal to explain molecular dynamics results for this variant because to lost hydrogen bonding interactions between the residue pairs K854SC/D614SC, Y837MC/D614MC, K835SC/D614SC, T859SC/D614SC. In addition, SpikeBA.1 residues near the furin cleavage site are more flexible than in SpikeWT, probably due to P681H and D614G substitutions. Finally, dynamical cross-correlation matrix (DCCM) analysis reveals that D614G and P681H may allosterically affect the cleavage site S1/S2. Conversely, S2' site may be influenced by residues located between NTD and RBD of a neighboring protomer of the SpikeWT. Such communication may be lost in SpikeBA.1, explaining the changes of the cell tropism in the viral infection. In addition, the movements of the NTDWT and NTDBA.1 may modulate the RBD conformation through allosteric effects. Taken together, our results explain how the structural aspects may explain the observed gains in infectivity, immune system evasion and transmissibility of the Omicron variant.Communicated by Ramaswamy H. Sarma.

Association of Social Vulnerability and COVID-19 Mortality Rates in Texas between 15 March 2020, and 21 July 2022: An Ecological Analysis

BACKGROUND

Despite the key role of social vulnerability such as economic disadvantage in health outcomes, research is limited on the impact of social vulnerabilities on COVID-19-related deaths, especially at the state and county level in the USA.

METHODS

We conducted a cross-sectional ecologic analysis of COVID-19 mortality by the county-level Minority Health Social Vulnerability Index (MH SVI) and each of its components in Texas. Negative binomial regression (NBR) analyses were used to estimate the association between the composite MH SVI (and its components) and COVID-19 mortality.

RESULTS

A 0.1-unit increase in the overall MH SVI (IRR, 1.27; 95% CI, 1.04-1.55; p = 0.017) was associated with a 27% increase in the COVID-19 mortality rate. Among the MH SVI component measures, only low socioeconomic status (IRR, 1.55; 95% CI, 1.28-1.89; p = 0.001) and higher household composition (e.g., proportion of older population per county) and disability scores (IRR, 1.47; 95% CI, 1.29-1.68; p < 0.001) were positively associated with COVID-19 mortality rates.

CONCLUSIONS

This study provides further evidence of disparities in COVID-19 mortality by social vulnerability and can inform decisions on the allocation of social resources and services as a strategy for reducing COVID-19 mortality rates and similar pandemics in the future.

Molecular insights and optimization strategies for the competitive binding of engineered ACE2 proteins: a multiple replica molecular dynamics study

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) continues to spread globally, and rapid viral evolution and the emergence of new variants pose challenges to pandemic control. During infection, the spike protein of SARS-CoV-2 interacts with the human ACE2 protein via its receptor binding domain (RBD), and it is known that engineered forms of ACE2 can compete with wild-type (WT) ACE2 for binding to inhibit infection. Here, we conducted multiple replica molecular dynamics (MRMD) simulations to study the mechanisms of the engineered ACE2 variants 3N39 and 3N94 and provide directions for optimization. Our findings reveal that engineered ACE2 is notably more efficacious in systems that show weaker binding to WT ACE2 (i.e., WT and BA.1 RBD), but also faces immune escape as the virus evolves. Moreover, by modifying residue types near the binding interface, engineered ACE2 alters the electrostatic potential distribution and reconfigures the hydrogen bonding network, which results in modified binding to the RBD. However, this structural rearrangement does not occur in all RBD variants. In addition, we identified potentially engineerable beneficial residues and potentially engineerable detrimental residues in both ACE2 and RBD. Functional conservation can thus enable the optimization of these residues and improve the binding competitiveness of engineered ACE2, which therefore provides additional immune escape prevention. Finally, we conclude that these findings have implications for understanding the mechanisms responsible for engineered ACE2 and can help us to develop engineered ACE2 proteins that show superior performance.

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.

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.

De Novo Human Angiotensin-Converting Enzyme 2 Decoy NL-CVX1 Protects Mice From Severe Disease After Severe Acute Respiratory Syndrome Coronavirus 2 Infection

The emergence of novel variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) underscores the need to investigate alternative approaches to prevent infection and treat patients with coronavirus disease 2019. Here, we report the preclinical efficacy of NL-CVX1, a de novo decoy that blocks virus entry into cells by binding with nanomolar affinity and high specificity to the receptor-binding domain of the SARS-CoV-2 spike protein. Using a transgenic mouse model of SARS-CoV-2 infection, we showed that a single prophylactic intranasal dose of NL-CVX1 conferred complete protection from severe disease following SARS-CoV-2 infection. Multiple therapeutic administrations of NL-CVX1 also protected mice from succumbing to infection. Finally, we showed that infected mice treated with NL-CVX1 developed both anti-SARS-CoV-2 antibodies and memory T cells and were protected against reinfection a month after treatment. Overall, these observations suggest NL-CVX1 is a promising therapeutic candidate for preventing and treating severe SARS-CoV-2 infections.

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.

The Effect of Select SARS-CoV-2 N-Linked Glycan and Variant of Concern Spike Protein Mutations on C-Type Lectin-Receptor-Mediated Infection

The SARS-CoV-2 virion has shown remarkable resilience, capable of mutating to escape immune detection and re-establishing infectious capabilities despite new vaccine rollouts. Therefore, there is a critical need to identify relatively immutable epitopes on the SARS-CoV-2 virion that are resistant to future mutations the virus may accumulate. While hACE2 has been identified as the receptor that mediates SARS-CoV-2 susceptibility, it is only modestly expressed in lung tissue. C-type lectin receptors like DC-SIGN can act as attachment sites to enhance SARS-CoV-2 infection of cells with moderate or low hACE2 expression. We developed an easy-to-implement assay system that allows for the testing of SARS-CoV-2 trans-infection. Using our assay, we assessed how SARS-CoV-2 Spike S1-domain glycans and spike proteins from different strains affected the ability of pseudotyped lentivirions to undergo DC-SIGN-mediated trans-infection. Through our experiments with seven glycan point mutants, two glycan cluster mutants and four strains of SARS-CoV-2 spike, we found that glycans N17 and N122 appear to have significant roles in maintaining COVID-19's infectious capabilities. We further found that the virus cannot retain infectivity upon the loss of multiple glycosylation sites, and that Omicron BA.2 pseudovirions may have an increased ability to bind to other non-lectin receptor proteins on the surface of cells. Taken together, our work opens the door to the development of new therapeutics that can target overlooked epitopes of the SARS-CoV-2 virion to prevent C-type lectin-receptor-mediated trans-infection in lung tissue.

Disulfide stabilization reveals conserved dynamic features between SARS-CoV-1 and SARS-CoV-2 spikes

SARS-CoV-2 spike protein (S) is structurally dynamic and has been observed by cryo-EM to adopt a variety of prefusion conformations that can be categorized as locked, closed, and open. S-trimers adopting locked conformations are tightly packed featuring structural elements incompatible with RBD in the "up" position. For SARS-CoV-2 S, it has been shown that the locked conformations are transient under neutral pH. Probably because of their transience, locked conformations remain largely uncharacterized for SARS-CoV-1 S. In this study, we introduced x1, x2, and x3 disulfides into SARS-CoV-1 S. Some of these disulfides have been shown to preserve rare locked conformations when introduced to SARS-CoV-2 S. Introduction of these disulfides allowed us to image a variety of locked and other rare conformations for SARS-CoV-1 S by cryo-EM. We identified bound cofactors and structural features that are associated with SARS-CoV-1 S locked conformations. We compare newly determined structures with other available spike structures of SARS-related CoVs to identify conserved features and discuss their possible functions.

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.

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.

Genomic surveillance of genes encoding the SARS-CoV-2 spike protein to monitor for emerging variants on Jeju Island, Republic of Korea

Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has been fueled by new variants emerging from circulating strains. Here, we report results from a genomic surveillance study of SARS-CoV-2 on Jeju Island, Republic of Korea, from February 2021 to September 2022.

Methods

A total of 3,585 SARS-CoV-2 positive samples were analyzed by Sanger sequencing of the gene encoding the spike protein before performing phylogenetic analyses.

Results

We found that the Alpha variant (B.1.1.7) was dominant in May 2021 before being replaced by the Delta variant (B.1.617.2) in July 2021, which was dominant until December 2021 before being replaced by the Omicron variant. Mutations in the spike protein, including N440K and G446S, have been proposed to contribute to immune evasion, accelerating the spread of Omicron variants.

Discussion

Our results from Juju Island, Republic of Korea, are consistent with and contribute to global surveillance efforts crucial for identifying new variants of concern of SARS-CoV-2 and for monitoring the transmission dynamics and characteristics of known strains.

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.

Oxalic acid blocked the binding of spike protein from SARS-CoV-2 Delta (B.1.617.2) and Omicron (B.1.1.529) variants to human angiotensin-converting enzymes 2

An epidemic of Corona Virus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is spreading worldwide. Moreover, the emergence of SARS-CoV-2 variants of concern, such as Delta and Omicron, has seriously challenged the application of current therapeutics including vaccination and drugs. Relying on interaction of spike protein with receptor angiotensin-converting enzymes 2 (ACE2), SARS-CoV-2 successfully invades to the host cells, which indicates a strategy that identification of small-molecular compounds to block the entry is of great significance for COVID-19 prevention. Our study evaluated the potential efficacy of natural compound oxalic acid (OA) as an inhibitory agent against SARS-CoV-2 invasion, particular on the interaction of the receptor binding domain (RBD) of Delta and Omicron variants to ACE2. By employing a competitive binding assay in vitro, OA significantly blocked the binding of RBDs from Delta B.1.617.2 and Omicron B.1.1.529 to ACE2, but has no effect on the wide-type SARS-CoV-2 strain. Furthermore, OA inhibited the entries of Delta and Omicron pseudovirus into ACE2 high expressing-HEK293T cells. By surface plasmon resonance (SPR) assay, the direct bindings of OA to RBD and ACE2 were analyzed and OA had both affinities with RBDs of B.1.617.2 and B.1.1.529 and with ACE2. Molecular docking predicted the binding sites on the RBD-ACE2 complex and it showed similar binding abilities to both complex of variant Delta or Omicron RBD and ACE2. In conclusion, we provided a promising novel small-molecule compound OA as an antiviral candidate by blocking the cellular entries of SARS-CoV-2 variants.

The SARS-CoV-2 spike protein structure: a symmetry tale on distortion trail

A preliminary step in the SARS-CoV-2 human infection process is a conformational change of the receptor binding domain (RBD) of its spike protein, characterized by a significant loss of symmetry. During this process, the residues which later on bind to the human angiotensin converting enzyme 2 (ACE2) receptor, become exposed at the surface of the protein. Symmetry analysis of a data set of 33 protein structures from experimental measurements and 32 structures from molecular dynamics simulation, show that the initial state carries clear indications on the structure of the final state, with respect to the local distortion along the sequence. This surprising finding implies that this type of analysis predicts the mechanism of change. We further show that the level of local distortion at the initial state increases with variant's transmissibility, for the wild type (WT) along with past and present variants of concern (WT ∼ alpha < beta < delta < Omicron BA.1), in accordance with the trend of their evolutionary path. In other words, the initial structure of the variant which is most infectious is also the most distorted, making its path to the final state shorter. It has been claimed that the RBD migration of the spike protein is allosterically controlled. Our analysis provides a quantitative support to a major theorem in this respect - that information about an allosteric process is encoded in the structure itself, suggesting that the path of local distortion is related to an allosteric information network.

Omicron Spike Protein Is Vulnerable to Reduction

SARS-CoV-2 virus spike (S) protein is an envelope protein responsible for binding to the ACE2 receptor, driving subsequent entry into host cells. The existence of multiple disulfide bonds in the S protein makes it potentially susceptible to reductive cleavage. Using a tri-part split luciferase-based binding assay, we evaluated the impacts of chemical reduction on S proteins from different virus variants and found that those from the Omicron family are highly vulnerable to reduction. Through manipulation of different Omicron mutations, we found that alterations in the receptor binding module (RBM) are the major determinants of this vulnerability. Specifically we discovered that Omicron mutations facilitate the cleavage of C480-C488 and C379-C432 disulfides, which consequently impairs binding activity and protein stability. The vulnerability of Omicron S proteins suggests a mechanism that can be harnessed to treat specific SARS-CoV-2 strains.

Site of vulnerability on SARS-CoV-2 spike induces broadly protective antibody against antigenically distinct Omicron subvariants

The rapid evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variants has emphasized the need to identify antibodies with broad neutralizing capabilities to inform future monoclonal therapies and vaccination strategies. Herein, we identified S728-1157, a broadly neutralizing antibody (bnAb) targeting the receptor-binding site (RBS) that was derived from an individual previously infected with WT SARS-CoV-2 prior to the spread of variants of concern (VOCs). S728-1157 demonstrated broad cross-neutralization of all dominant variants, including D614G, Beta, Delta, Kappa, Mu, and Omicron (BA.1/BA.2/BA.2.75/BA.4/BA.5/BL.1/XBB). Furthermore, S728-1157 protected hamsters against in vivo challenges with WT, Delta, and BA.1 viruses. Structural analysis showed that this antibody targets a class 1/RBS-A epitope in the receptor binding domain via multiple hydrophobic and polar interactions with its heavy chain complementarity determining region 3 (CDR-H3), in addition to common motifs in CDR-H1/CDR-H2 of class 1/RBS-A antibodies. Importantly, this epitope was more readily accessible in the open and prefusion state, or in the hexaproline (6P)-stabilized spike constructs, as compared with diproline (2P) constructs. Overall, S728-1157 demonstrates broad therapeutic potential and may inform target-driven vaccine designs against future SARS-CoV-2 variants.

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.

Identification of a conserved S2 epitope present on spike proteins from all highly pathogenic coronaviruses

To address the ongoing SARS-CoV-2 pandemic and prepare for future coronavirus outbreaks, understanding the protective potential of epitopes conserved across SARS-CoV-2 variants and coronavirus lineages is essential. We describe a highly conserved, conformational S2 domain epitope present only in the prefusion core of β-coronaviruses: SARS-CoV-2 S2 apex residues 980-1006 in the flexible hinge. Antibody RAY53 binds the native hinge in MERS-CoV and SARS-CoV-2 spikes on the surface of mammalian cells and mediates antibody-dependent cellular phagocytosis and cytotoxicity against SARS-CoV-2 spike in vitro. Hinge epitope mutations that ablate antibody binding compromise pseudovirus infectivity, but changes elsewhere that affect spike opening dynamics, including those found in Omicron BA.1, occlude the epitope and may evade pre-existing serum antibodies targeting the S2 core. This work defines a third class of S2 antibody while providing insights into the potency and limitations of S2 core epitope targeting.

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.

Natural selection differences detected in key protein domains between non-pathogenic and pathogenic feline coronavirus phenotypes

Feline coronaviruses (FCoVs) commonly cause mild enteric infections in felines worldwide (termed feline enteric coronavirus [FECV]), with around 12 per cent developing into deadly feline infectious peritonitis (FIP; feline infectious peritonitis virus [FIPV]). Genomic differences between FECV and FIPV have been reported, yet the putative genotypic basis of the highly pathogenic phenotype remains unclear. Here, we used state-of-the-art molecular evolutionary genetic statistical techniques to identify and compare differences in natural selection pressure between FECV and FIPV sequences, as well as to identify FIPV- and FECV-specific signals of positive selection. We analyzed full-length FCoV protein coding genes thought to contain mutations associated with FIPV (Spike, ORF3abc, and ORF7ab). We identified two sites exhibiting differences in natural selection pressure between FECV and FIPV: one within the S1/S2 furin cleavage site (FCS) and the other within the fusion domain of Spike. We also found fifteen sites subject to positive selection associated with FIPV within Spike, eleven of which have not previously been suggested as possibly relevant to FIP development. These sites fall within Spike protein subdomains that participate in host cell receptor interaction, immune evasion, tropism shifts, host cellular entry, and viral escape. There were fourteen sites (twelve novel sites) within Spike under positive selection associated with the FECV phenotype, almost exclusively within the S1/S2 FCS and adjacent to C domain, along with a signal of relaxed selection in FIPV relative to FECV, suggesting that furin cleavage functionality may not be needed for FIPV. Positive selection inferred in ORF7b was associated with the FECV phenotype and included twenty-four positively selected sites, while ORF7b had signals of relaxed selection in FIPV. We found evidence of positive selection in ORF3c in FCoV-wide analyses, but no specific association with the FIPV or FECV phenotype. We hypothesize that some combination of mutations in FECV may contribute to FIP development, and that it is unlikely to be one singular 'switch' mutational event. This work expands our understanding of the complexities of FIP development and provides insights into how evolutionary forces may alter pathogenesis in coronavirus genomes.

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.

Heterologous boost with mRNA vaccines against SARS-CoV-2 Delta/Omicron variants following an inactivated whole-virus vaccine

The coronavirus SARS-CoV-2 has mutated quickly and caused significant global damage. This study characterizes two mRNA vaccines ZSVG-02 (Delta) and ZSVG-02-O (Omicron BA.1), and associating heterologous prime-boost strategy following the prime of a most widely administrated inactivated whole-virus vaccine (BBIBP-CorV). The ZSVG-02-O induces neutralizing antibodies that effectively cross-react with Omicron subvariants. In naïve animals, ZSVG-02 or ZSVG-02-O induce humoral responses skewed to the vaccine's targeting strains, but cellular immune responses cross-react to all variants of concern (VOCs) tested. Following heterologous prime-boost regimes, animals present comparable neutralizing antibody levels and superior protection against Delta and Omicron BA.1variants. Single-boost only generated ancestral and omicron dual-responsive antibodies, probably by "recall" and "reshape" the prime immunity. New Omicron-specific antibody populations, however, appeared only following the second boost with ZSVG-02-O. Overall, our results support a heterologous boost with ZSVG-02-O, providing the best protection against current VOCs in inactivated virus vaccine-primed populations.

Variants of SARS-CoV-2: Influences on the Vaccines' Effectiveness and Possible Strategies to Overcome Their Consequences

The immune response elicited by the current COVID-19 vaccinations declines with time, especially among the immunocompromised population. Furthermore, the emergence of novel SARS-CoV-2 variants, particularly the Omicron variant, has raised serious concerns about the efficacy of currently available vaccines in protecting the most vulnerable people. Several studies have reported that vaccinated people get breakthrough infections amid COVID-19 cases. So far, five variants of concern (VOCs) have been reported, resulting in successive waves of infection. These variants have shown a variable amount of resistance towards the neutralising antibodies (nAbs) elicited either through natural infection or the vaccination. The spike (S) protein, membrane (M) protein, and envelope (E) protein on the viral surface envelope and the N-nucleocapsid protein in the core of the ribonucleoprotein are the major structural vaccine target proteins against COVID-19. Among these targets, S Protein has been extensively exploited to generate effective vaccines against COVID-19. Hence, amid the emergence of novel variants of SARS-CoV-2, we have discussed their impact on currently available vaccines. We have also discussed the potential roles of S Protein in the development of novel vaccination approaches to contain the negative consequences of the variants' emergence and acquisition of mutations in the S Protein of SARS-CoV-2. Moreover, the implications of SARS-CoV-2's structural proteins were also discussed in terms of their variable potential to elicit an effective amount of immune response.

Variations within the Glycan Shield of SARS-CoV-2 Impact Viral Spike Dynamics

The emergence of SARS-CoV-2 variants alters the efficacy of existing immunity, whether arisen naturally or through vaccination. Understanding the structure of the viral spike assists in determining the impact of mutations on the antigenic surface. One class of mutation impacts glycosylation attachment sites, which have the capacity to influence the antigenic structure beyond the immediate site of attachment. Here, we compare the site-specific glycosylation of recombinant viral spike mimetics of B.1.351 (Beta), P.1 (Gamma), B.1.617.2 (Delta), B.1.1.529 (Omicron). The P.1 strain exhibits two additional N-linked glycan sites compared to the other variants analyzed and we investigate the impact of these glycans by molecular dynamics. The acquired N188 site is shown to exhibit very limited glycan maturation, consistent with limited enzyme accessibility. Structural modeling and molecular dynamics reveal that N188 is located within a cavity by the receptor binding domain, which influences the dynamics of these attachment domains. These observations suggest a mechanism whereby mutations affecting viral glycosylation sites have a structural impact across the protein surface.

Computational analysis of the sequence-structure relation in SARS-CoV-2 spike protein using protein contact networks

The structure of proteins impacts directly on the function they perform. Mutations in the primary sequence can provoke structural changes with consequent modification of functional properties. SARS-CoV-2 proteins have been extensively studied during the pandemic. This wide dataset, related to sequence and structure, has enabled joint sequence-structure analysis. In this work, we focus on the SARS-CoV-2 S (Spike) protein and the relations between sequence mutations and structure variations, in order to shed light on the structural changes stemming from the position of mutated amino acid residues in three different SARS-CoV-2 strains. We propose the use of protein contact network (PCN) formalism to: (i) obtain a global metric space and compare various molecular entities, (ii) give a structural explanation of the observed phenotype, and (iii) provide context dependent descriptors of single mutations. PCNs have been used to compare sequence and structure of the Alpha, Delta, and Omicron SARS-CoV-2 variants, and we found that omicron has a unique mutational pattern leading to different structural consequences from mutations of other strains. The non-random distribution of changes in network centrality along the chain has allowed to shed light on the structural (and functional) consequences of mutations.

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.

Distinct Conformations of SARS-CoV-2 Omicron Spike Protein and Its Interaction with ACE2 and Antibody

Since November 2021, Omicron has been the dominant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant that causes the coronavirus disease 2019 (COVID-19) and has continuously impacted human health. Omicron sublineages are still increasing and cause increased transmission and infection rates. The additional 15 mutations on the receptor binding domain (RBD) of Omicron spike proteins change the protein conformation, enabling the Omicron variant to evade neutralizing antibodies. For this reason, many efforts have been made to design new antigenic variants to induce effective antibodies in SARS-CoV-2 vaccine development. However, understanding the different states of Omicron spike proteins with and without external molecules has not yet been addressed. In this review, we analyze the structures of the spike protein in the presence and absence of angiotensin-converting enzyme 2 (ACE2) and antibodies. Compared to previously determined structures for the wildtype spike protein and other variants such as alpha, beta, delta, and gamma, the Omicron spike protein adopts a partially open form. The open-form spike protein with one RBD up is dominant, followed by the open-form spike protein with two RBD up, and the closed-form spike protein with the RBD down. It is suggested that the competition between antibodies and ACE2 induces interactions between adjacent RBDs of the spike protein, which lead to a partially open form of the Omicron spike protein. The comprehensive structural information of Omicron spike proteins could be helpful for the efficient design of vaccines against the Omicron variant.

Ancestral Lineage of SARS-CoV-2 Is More Stable in Human Biological Fluids than Alpha, Beta, and Omicron Variants of Concern

SARS-CoV-2 is a zoonotic virus first identified in 2019, and has quickly spread worldwide. The virus is primarily transmitted through respiratory droplets from infected persons; however, the virus-laden excretions can contaminate surfaces which can serve as a potential source of infection. Since the beginning of the pandemic, SARS-CoV-2 has continued to evolve and accumulate mutations throughout its genome leading to the emergence of variants of concern (VOCs) which exhibit increased fitness, transmissibility, and/or virulence. However, the stability of SARS-CoV-2 VOCs in biological fluids has not been thoroughly investigated. The aim of this study was to determine and compare the stability of different SARS-CoV-2 strains in human biological fluids. Here, we demonstrate that the ancestral strain of the Wuhan-like lineage A was more stable than the Alpha VOC B.1.1.7, and the Beta VOC B.1.351 strains in human liquid nasal mucus and sputum. In contrast, there was no difference in stability among the three strains in dried biological fluids. Furthermore, we also show that the Omicron VOC B.1.1.529 strain was less stable than the ancestral Wuhan-like strain in liquid nasal mucus. These studies provide insight into the effect of the molecular evolution of SARS-CoV-2 on environmental virus stability, which is important information for the development of countermeasures against SARS-CoV-2. IMPORTANCE Genetic evolution of SARS-CoV-2 leads to the continuous emergence of novel virus variants, posing a significant concern to global public health. Five of these variants have been classified to date into variants of concern (VOCs); Alpha, Beta, Gamma, Delta, and Omicron. Previous studies investigated the stability of SARS-CoV-2 under various conditions, but there is a gap of knowledge on the survival of SARS-CoV-2 VOCs in human biological fluids which are clinically relevant. Here, we present evidence that Alpha, Beta, and Omicron VOCs were less stable than the ancestral Wuhan-like strain in human biological fluids. Our findings highlight the potential risk of contaminated human biological fluids in SARS-CoV-2 transmission and contribute to the development of countermeasures against SARS-CoV-2.

SARS-CoV-2 Omicron Subvariants Balance Host Cell Membrane, Receptor, and Antibody Docking via an Overlapping Target Site

Variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are emerging rapidly and offer surfaces that are optimized for recognition of host cell membranes while also evading antibodies arising from vaccinations and previous infections. Host cell infection is a multi-step process in which spike heads engage lipid bilayers and one or more angiotensin-converting enzyme 2 (ACE-2) receptors. Here, the membrane binding surfaces of Omicron subvariants are compared using cryo-electron microscopy (cEM) structures of spike trimers from BA.2, BA.2.12.1, BA.2.13, BA.2.75, BA.3, BA.4, and BA.5 viruses. Despite significant differences around mutated sites, they all maintain strong membrane binding propensities that first appeared in BA.1. Both their closed and open states retain elevated membrane docking capacities, although the presence of more closed than open states diminishes opportunities to bind receptors while enhancing membrane engagement. The electrostatic dipoles are generally conserved. However, the BA.2.75 spike dipole is compromised, and its ACE-2 affinity is increased, and BA.3 exhibits the opposite pattern. We propose that balancing the functional imperatives of a stable, readily cleavable spike that engages both lipid bilayers and receptors while avoiding host defenses underlies betacoronavirus evolution. This provides predictive criteria for rationalizing future pandemic waves and COVID-19 transmissibility while illuminating critical sites and strategies for simultaneously combating multiple variants.

Smart healthcare: A prospective future medical approach for COVID-19

COVID-19 has greatly affected human life for over 3 years. In this review, we focus on smart healthcare solutions that address major requirements for coping with the COVID-19 pandemic, including (1) the continuous monitoring of severe acute respiratory syndrome coronavirus 2, (2) patient stratification with distinct short-term outcomes (eg, mild or severe diseases) and long-term outcomes (eg, long COVID), and (3) adherence to medication and treatments for patients with COVID-19. Smart healthcare often utilizes medical artificial intelligence (AI) and cloud computing and integrates cutting-edge biological and optoelectronic techniques. These are valuable technologies for addressing the unmet needs in the management of COVID. By leveraging deep learning/machine learning capabilities and big data, medical AI can perform precise prognosis predictions and provide reliable suggestions for physicians' decision-making. Through the assistance of the Internet of Medical Things, which encompasses wearable devices, smartphone apps, internet-based drug delivery systems, and telemedicine technologies, the status of mild cases can be continuously monitored and medications provided at home without the need for hospital care. In cases that develop into severe cases, emergency feedback can be provided through the hospital for rapid treatment. Smart healthcare can possibly prevent the development of severe COVID-19 cases and therefore lower the burden on intensive care units.

Long-term passaging of pseudo-typed SARS-CoV-2 reveals the breadth of monoclonal and bispecific antibody cocktails

The continuous emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants poses challenges to the effectiveness of neutralizing antibodies. Rational design of antibody cocktails is a realizable approach addressing viral immune evasion. However, evaluating the breadth of antibody cocktails is essential for understanding the development potential. Here, based on a replication competent vesicular stomatitis virus model that incorporates the spike of SARS-CoV-2 (VSV-SARS-CoV-2), we evaluated the breadth of a number of antibody cocktails consisting of monoclonal antibodies and bispecific antibodies by long-term passaging the virus in the presence of the cocktails. Results from over two-month passaging of the virus showed that 9E12 + 10D4 + 2G1 and 7B9-9D11 + 2G1 from these cocktails were highly resistant to random mutation, and there was no breakthrough after 30 rounds of passaging. As a control, antibody REGN10933 was broken through in the third passage. Next generation sequencing was performed and several critical mutations related to viral evasion were identified. These mutations caused a decrease in neutralization efficiency, but the reduced replication rate and ACE2 susceptibility of the mutant virus suggested that they might not have the potential to become epidemic strains. The 9E12 + 10D4 + 2G1 and 7B9-9D11 + 2G1 cocktails that picked from the VSV-SARS-CoV-2 system efficiently neutralized all current variants of concern and variants of interest including the most recent variants Delta and Omicron, as well as SARS-CoV-1. Our results highlight the feasibility of using the VSV-SARS-CoV-2 system to develop SARS-CoV-2 antibody cocktails and provide a reference for the clinical selection of therapeutic strategies to address the mutational escape of SARS-CoV-2.

Network analysis uncovers the communication structure of SARS-CoV-2 spike protein identifying sites for immunogen design

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has triggered myriad efforts to understand the structure and dynamics of this complex pathogen. The spike glycoprotein of SARS-CoV-2 is a significant target for immunogens as it is the means by which the virus enters human cells, while simultaneously sporting mutations responsible for immune escape. These functional and escape processes are regulated by complex molecular-level interactions. Our study presents quantitative insights on domain and residue contributions to allosteric communication, immune evasion, and local- and global-level control of functions through the derivation of a weighted graph representation from all-atom MD simulations. Focusing on the ancestral form and the D614G-variant, we provide evidence of the utility of our approach by guiding the selection of a mutation that alters the spike's stability. Taken together, the network approach serves as a valuable tool to evaluate communication "hot-spots" in proteins to guide design of stable immunogens.

Long-term variations and potency of neutralizing antibodies against Omicron subvariants after CoronaVac-inactivated booster: A 7-month follow-up study

The long-term protective efficacy of neutralizing antibodies (Nabs) against Omicron subvariants after inactivated booster vaccines remains elusive. During the follow-up study, 54 healthy volunteers aged 20-31 years received inactivated CoronaVac booster vaccinations and were monitored for 221 days. The dynamic efficacy and durability of Nab against Omicron subvariants BA.1, BA.2, BA.2.12.2, and BA4/5 were assessed using a pseudotyped virus neutralization assay at up to nine time points post immunization. The antibody response against Omicron subvariants was substantially weaker than D614G, with BA.4/5 being the least responsive. The geometric mean titer (GMT) of Nab against Omicron subvariants BA.1, BA.2, BA.2.12.1, and BA.4/5 was 2.2-, 1.7-, 1.8-, and 2.2-fold lower than that against D614G (p_s  < 0.0001). The gap in Nab response between Omicron subvariants was pronounced during the 2 weeks-2 months following booster vaccination (p_s  < 0.05). Seven months post booster, the antibody potency against D614G was maintained at 100% (50% for Nab titers ≥ 100 50% inhibitory dilution [EC₅₀ ]), whereas at 77.3% for BA.1, 90.9% for BA.2, 86.4% for BA.2.12.1, and 86.4% for BA.4/5 (almost 20% for Nab titers ≥ 100 EC₅₀ ). Despite the inevitable immune escape, Omicron subvariants maintained sustained and measurable antibody potency post-booster vaccination during long-term monitoring, which could help optimize immunization strategies.

Force-tuned avidity of spike variant-ACE2 interactions viewed on the single-molecule level

Recent waves of COVID-19 correlate with the emergence of the Delta and the Omicron variant. We report that the Spike trimer acts as a highly dynamic molecular caliper, thereby forming up to three tight bonds through its RBDs with ACE2 expressed on the cell surface. The Spike of both Delta and Omicron (B.1.1.529) Variant enhance and markedly prolong viral attachment to the host cell receptor ACE2, as opposed to the early Wuhan-1 isolate. Delta Spike shows rapid binding of all three Spike RBDs to three different ACE2 molecules with considerably increased bond lifetime when compared to the reference strain, thereby significantly amplifying avidity. Intriguingly, Omicron (B.1.1.529) Spike displays less multivalent bindings to ACE2 molecules, yet with a ten time longer bond lifetime than Delta. Delta and Omicron (B.1.1.529) Spike variants enhance and prolong viral attachment to the host, which likely not only increases the rate of viral uptake, but also enhances the resistance of the variants against host-cell detachment by shear forces such as airflow, mucus or blood flow. We uncover distinct binding mechanisms and strategies at single-molecule resolution, employed by circulating SARS-CoV-2 variants to enhance infectivity and viral transmission.

The SARS-CoV-2 spike S375F mutation characterizes the Omicron BA.1 variant

Recent studies have revealed the unique virological characteristics of Omicron, particularly those of its spike protein, such as less cleavage efficacy in cells, reduced ACE2 binding affinity, and poor fusogenicity. However, it remains unclear which mutation(s) determine these three virological characteristics of Omicron spike. Here, we show that these characteristics of the Omicron spike protein are determined by its receptor-binding domain. Of interest, molecular phylogenetic analysis revealed that acquisition of the spike S375F mutation was closely associated with the explosive spread of Omicron in the human population. We further elucidated that the F375 residue forms an interprotomer pi-pi interaction with the H505 residue of another protomer in the spike trimer, conferring the attenuated cleavage efficiency and fusogenicity of Omicron spike. Our data shed light on the evolutionary events underlying the emergence of Omicron at the molecular level.

SARSNTdb database: Factors affecting SARS-CoV-2 sequence conservation

SARSNTdb offers a curated, nucleotide-centric database for users of varying levels of SARS-CoV-2 knowledge. Its user-friendly interface enables querying coding regions and coordinate intervals to find out the various functional and selective constraints that act upon the corresponding nucleotides and amino acids. Users can easily obtain information about viral genes and proteins, functional domains, repeats, secondary structure formation, intragenomic interactions, and mutation prevalence. Currently, many databases are focused on the phylogeny and amino acid substitutions, mainly in the spike protein. We took a novel, more nucleotide-focused approach as RNA does more than just code for proteins and many insights can be gleaned from its study. For example, RNA-targeted drug therapies for SARS-CoV-2 are currently being developed and it is essential to understand the features only visible at that level. This database enables the user to identify regions that are more prone to forming secondary structures that drugs can target. SARSNTdb also provides illustrative mutation data from a subset of ~25,000 patient samples with a reliable read coverage across the whole genome (from different locations and time points in the pandemic. Finally, the database allows for comparing SARS-CoV-2 and SARS-CoV domains and sequences. SARSNTdb can serve the research community by being a curated repository for information that gives a jump start to analyze a mutation’s effect far beyond just determining synonymous/non-synonymous substitutions in protein sequences.

Novel chimeric proteins mimicking SARS-CoV-2 spike epitopes with broad inhibitory activity

SARS-CoV-2 spike (S) protein mediates virus attachment to the cells and fusion between viral and cell membranes. Membrane fusion is driven by mutual interaction between the highly conserved heptad-repeat regions 1 and 2 (HR1 and HR2) of the S2 subunit of the spike. For this reason, these S2 regions are interesting therapeutic targets for COVID-19. Although HR1 and HR2 have been described as transiently exposed during the fusion process, no significant antibody responses against these S2 regions have been reported. Here we designed chimeric proteins that imitate highly stable HR1 helical trimers and strongly bind to HR2. The proteins have broad inhibitory activity against WT B.1 and BA.1 viruses. Sera from COVID-19 convalescent donors showed significant levels of reactive antibodies (IgG and IgA) against the HR1 mimetic proteins, whereas these antibody responses were absent in sera from uninfected donors. Moreover, both inhibitory activity and antigenicity of the proteins correlate positively with their structural stability but not with the number of amino acid changes in their HR1 sequences, indicating a conformational and conserved nature of the involved epitopes. Our results reveal previously undetected spike epitopes that may guide the design of new robust COVID-19 vaccines and therapies.