Membrane fusion and immune evasion by the spike protein of SARS-CoV-2 Delta variant

Overview of immunological & virological factors driving the evolution & global spread of SARS-CoV-2 variants

The SARS-CoV-2, a highly infectious positive strand RNA virus first identified in December 2019, has produced multiple genetic variants that have rapidly and sequentially spread worldwide during the coronavirus disease 2019 (COVID-19) pandemic. Genetic changes in SARS-CoV-2 for greater infectivity, replication and transmission were selected during the early stages of the pandemic. More recently, after widespread infection and vaccination, SARS-CoV-2 variants that evade antigen-specific adaptive immunity, have begun to be selected. This article provides an overview of the molecular immunological and virological factors underlying the origin and global spread of important SARS-CoV-2 variant lineages.

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.

Neutralization against SARS-CoV-2 Delta/Omicron variants and B cell response after inactivated vaccination among COVID-19 convalescents

Emerging SARS-CoV-2 variants have made COVID-19 convalescents susceptible to re-infection and have raised concern about the efficacy of inactivated vaccination in neutralization against emerging variants and antigen-specific B cell response. To this end, a study on a long-term cohort of 208 participants who have recovered from COVID-19 was conducted, and the participants were followed up at 3.3 (Visit 1), 9.2 (Visit 2), and 18.5 (Visit 3) months after SARS-CoV-2 infection. They were classified into three groups (no-vaccination (n = 54), one-dose (n = 62), and two-dose (n = 92) groups) on the basis of the administration of inactivated vaccination. The neutralizing antibody (NAb) titers against the wild-type virus continued to decrease in the no-vaccination group, but they rose significantly in the one-dose and two-dose groups, with the highest NAb titers being observed in the two-dose group at Visit 3. The NAb titers against the Delta variant for the no-vaccination, one-dose, and two-dose groups decreased by 3.3, 1.9, and 2.3 folds relative to the wild-type virus, respectively, and those against the Omicron variant decreased by 7.0, 4.0, and 3.8 folds, respectively. Similarly, the responses of SARS-CoV-2 RBD-specific B cells and memory B cells were boosted by the second vaccine dose. Results showed that the convalescents benefited from the administration of the inactivated vaccine (one or two doses), which enhanced neutralization against highly mutated SARS-CoV-2 variants and memory B cell responses. Two doses of inactivated vaccine among COVID-19 convalescents are therefore recommended for the prevention of the COVID-19 pandemic, and vaccination guidelines and policies need to be updated.

Free Energy Perturbation Calculations of Mutation Effects on SARS-CoV-2 RBD::ACE2 Binding Affinity

The strength of binding between human angiotensin converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of viral spike protein plays a role in the transmissibility of the SARS-CoV-2 virus. In this study we focus on a subset of RBD mutations that have been frequently observed in infected individuals and probe binding affinity changes to ACE2 using surface plasmon resonance (SPR) measurements and free energy perturbation (FEP) calculations. Our SPR results are largely in accord with previous studies but discrepancies do arise due to differences in experimental methods and to protocol differences even when a single method is used. Overall, we find that FEP performance is superior to that of other computational approaches examined as determined by agreement with experiment and, in particular, by its ability to identify stabilizing mutations. Moreover, the calculations successfully predict the observed cooperative stabilization of binding by the Q498R N501Y double mutant present in Omicron variants and offer a physical explanation for the underlying mechanism. Overall, our results suggest that despite the significant computational cost, FEP calculations may offer an effective strategy to understand the effects of interfacial mutations on protein-protein binding affinities and, hence, in a variety of practical applications such as the optimization of neutralizing antibodies.

Structural effects of spike protein D614G mutation in SARS-CoV-2

A single mutation from aspartate to glycine at position 614 has dominated all circulating variants of the severe acute respiratory syndrome coronavirus 2. D614G mutation induces structural changes in the spike (S) protein that strengthen the virus infectivity. Here, we use molecular dynamics simulations to dissect the effects of mutation and 630-loop rigidification on S-protein structure. The introduction of the mutation orders the 630-loop structure and thereby induces global structural changes toward the cryoelectron microscopy structure of the D614G S-protein. The ordered 630-loop weakens local interactions between the 614th residue and others in contrast to disordered structures in the wild-type protein. The mutation allosterically alters global interactions between receptor-binding domains, forming an asymmetric and mobile down conformation and facilitating transitions toward up conformation. The loss of salt bridge between D614 and K854 upon the mutation generally stabilizes S-protein protomer, including the fusion peptide proximal region that mediates membrane fusion. Understanding the molecular basis of D614G mutation is crucial as it dominates in all variants of concern, including Delta and Omicron.

SARS-CoV-2 Spike Protein Is Capable of Inducing Cell-Cell Fusions Independent from Its Receptor ACE2 and This Activity Can Be Impaired by Furin Inhibitors or a Subset of Monoclonal Antibodies

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which was responsible for the COVID-19 pandemic, efficiently spreads cell-to-cell through mechanisms facilitated by its membrane glycoprotein spike. We established a dual split protein (DSP) assay based on the complementation of GFP and luciferase to quantify the fusogenic activity of the SARS-CoV-2 spike protein. We provide several lines of evidence that the spike protein of SARS-CoV-2, but not SARS-CoV-1, induced cell-cell fusion even in the absence of its receptor, angiotensin-converting enzyme 2 (ACE2). This poorly described ACE2-independent cell fusion activity of the spike protein was strictly dependent on the proteasomal cleavage of the spike by furin while TMPRSS2 was dispensable. Previous and current variants of concern (VOCs) differed significantly in their fusogenicity. The Delta spike was extremely potent compared to Alpha, Beta, Gamma and Kappa, while the Omicron spike was almost devoid of receptor-independent fusion activity. Nonetheless, for all analyzed variants, cell fusion was dependent on furin cleavage and could be pharmacologically inhibited with CMK. Mapping studies revealed that amino acids 652-1273 conferred the ACE2-independent fusion activity of the spike. Unexpectedly, residues proximal to the furin cleavage site were not of major relevance, whereas residue 655 critically regulated fusion. Finally, we found that the spike's fusion activity in the absence of ACE2 could be inhibited by antibodies directed against its N-terminal domain (NTD) but not by antibodies targeting its receptor-binding domain (RBD). In conclusion, our BSL-1-compatible DSP assay allowed us to screen for inhibitors or antibodies that interfere with the spike's fusogenic activity and may therefore contribute to both rational vaccine design and development of novel treatment options against SARS-CoV-2.

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

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

The anti-SARS-CoV-2 effect and mechanism of Chiehyuan herbal oral protection solution

The Chiehyuan herbal oral protection solution (GB-2) is a herbal mixture commonly utilized in Taiwan for combating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as per traditional Chinese medicine practices. This study assessed the clinical impact of GB-2 through prospective clinical trials. With twice-daily use for a week, GB-2 was shown to diminish the expression of angiotensin-converting enzyme 2 (ACE2) in oral mucosal cells. Moreover, after two weeks of use, it could reduce transmembrane protease, serine 2 (TMRPSS2) expression in these cells. Additionally, in vitro experiments demonstrated that GB-2 lessened the entry efficiency of the Omicron, L452R-D614G, T478K-D614G, and L452R-T478K-D614G variants of the SARS-CoV-2 pseudotyped lentivirus. It also impeded the interaction between ACE2 and the receptor-binding domain (RBD) presenting N501Y-K417N-E484A-G339D-Q493R-G496S-Q498R and L452R-T478K mutations. Glycyrrhizic acid, a major compound in GB-2, also hindered the entry of the Omicron variant (BA.1) of the SARS-CoV-2 pseudotyped lentivirus by obstructing the binding between ACE2 and the RBD presenting the N501Y-K417N-E484A-G339D-Q493R-G496S-Q498R mutation. To sum up, these findings suggest that GB-2 can decrease ACE2 and TMPRSS2 expression in oral mucosal cells. Both glycyrrhizic acid and GB-2 were found to reduce the entry efficiency of the Omicron variant (BA.1) of the SARS-CoV-2 pseudotyped lentivirus and block the binding between ACE2 and the RBD with the N501Y-K417N-E484A-G339D-Q493R-G496S-Q498R mutation. This evidence implies that GB-2 might be a potential candidate for further study as a preventative measure against SARS-CoV-2 infection.

Differences between Omicron SARS-CoV-2 RBD and other variants in their ability to interact with cell receptors and monoclonal antibodies

SARS-CoV-2 remains a health threat with the continuous emergence of new variants. This work aims to expand the knowledge about the SARS-CoV-2 receptor-binding domain (RBD) interactions with cell receptors and monoclonal antibodies (mAbs). By using constant-pH Monte Carlo simulations, the free energy of interactions between the RBD from different variants and several partners (Angiotensin-Converting Enzyme-2 (ACE2) polymorphisms and various mAbs) were predicted. Computed RBD-ACE2-binding affinities were higher for two ACE2 polymorphisms (rs142984500 and rs4646116) typically found in Europeans which indicates a genetic susceptibility. This is amplified for Omicron (BA.1) and its sublineages BA.2 and BA.3. The antibody landscape was computationally investigated with the largest set of mAbs so far in the literature. From the 32 studied binders, groups of mAbs were identified from weak to strong binding affinities (e.g. S2K146). These mAbs with strong binding capacity and especially their combination are amenable to experimentation and clinical trials because of their high predicted binding affinities and possible neutralization potential for current known virus mutations and a universal coronavirus.Communicated by Ramaswamy H. Sarma.

Harringtonine: A more effective antagonist for Omicron variant

Fusion with host cell membrane is the main mechanism of infection of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here, we propose that a new strategy to screen small-molecule antagonists blocking SARS-CoV-2 membrane fusion. Using cell membrane chromatography (CMC), we found that harringtonine (HT) simultaneously targeted SARS-CoV-2 S protein and host cell surface TMPRSS2 expressed by the host cell, and subsequently confirmed that HT can inhibit membrane fusion. HT effectively blocked SARS-CoV-2 original strain entry with the IC50 of 0.217 μM, while the IC50 in delta variant decreased to 0.101 μM, the IC50 in Omicron BA.1 variant was 0.042 μM. Due to high transmissibility and immune escape, Omicron subvariant BA.5 has become the dominant strain of the SARS-CoV-2 virus and led to escalating COVID-19 cases, however, against BA.5, HT showed a surprising effectiveness. The IC50 in Omicron BA.5 was even lower than 0.0019 μM. The above results revealed the effect of HT on Omicron is very significant. In summary, we characterize HT as a small-molecule antagonist by direct targeting on the Spike protein and TMPRSS2.

Kaempferol inhibits SARS-CoV-2 invasion by impairing heptad repeats-mediated viral fusion

BACKGROUND

The continuous evolution of SARS-CoV-2 has underscored the development of broad-spectrum prophylaxis. Antivirals targeting the membrane fusion process represent promising paradigms. Kaempferol (Kae), an ubiquitous plant flavonol, has been shown efficacy against various enveloped viruses. However, its potential in anti-SARS-CoV-2 invasion remains obscure.

PURPOSE

To evaluate capabilities and mechanisms of Kae in preventing SARS-CoV-2 invasion.

METHODS

To avoid interference of viral replication, virus-like particles (VLPs) constructed with luciferase reporter were applied. To investigate the antiviral potency of Kae, human induced pluripotent stem cells (hiPSC)-derived alveolar epithelial cells type II (AECII) and human ACE2 (hACE2) transgenic mice were utilized as in vitro and in vivo models, respectively. Using dual split protein (DSP) assays, inhibitory activities of Kae in viral fusion were determined in Alpha, Delta and Omicron variants of SARS-CoV-2, as well as in SARS-CoV and MERS-CoV. To further reveal molecular determinants of Kae in restricting viral fusion, synthetic peptides corresponding to the conserved heptad repeat (HR) 1 and 2, involved in viral fusion, and the mutant form of HR2 were explored by circular dichroism and native polyacrylamide gel electrophoresis.

RESULTS

Kae inhibited SARS-CoV-2 invasion both in vitro and in vivo, which was mainly attributed to its suppressive effects on viral fusion, but not endocytosis, two pathways that mediate viral invasion. In accordance with the proposed model of anti-fusion prophylaxis, Kae functioned as a pan-inhibitor of viral fusion, including three emerged highly pathogenic coronaviruses, and the currently circulating Omicron BQ.1.1 and XBB.1 variants of SARS-CoV-2. Consistent with the typical target of viral fusion inhibitors, Kae interacted with HR regions of SARS-CoV-2 S2 subunits. Distinct from previous inhibitory fusion peptides which prevent the formation of six-helix bundle (6-HB) by competitively interacting with HRs, Kae deformed HR1 and directly reacted with lysine residues within HR2 region, the latter of which was considered critical for the preservation of stabilized S2 during SARS-CoV-2 invasion.

CONCLUSIONS

Kae prevents SARS-CoV-2 infection by blocking membrane fusion and possesses a broad-spectrum anti-fusion ability. These findings provide valuable insights into potential benefits of Kae-containing botanical products as a complementary prophylaxis, especially during the waves of breakthrough infections and re-infections.

Changes within the P681 residue of spike dictate cell fusion and syncytia formation of Delta and Omicron variants of SARS-CoV-2 with no effects on neutralization or infectivity

The rapid spread and dominance of the Omicron SARS-CoV-2 lineages have posed severe health challenges worldwide. While extensive research on the role of the Receptor Binding Domain (RBD) in promoting viral infectivity and vaccine sensitivity has been well documented, the functional significance of the ₆₈₁PRRAR/SV₆₈₇ polybasic motif of the viral spike is less clear. In this work, we monitored the infectivity levels and neutralization potential of the wild-type human coronavirus 2019 (hCoV-19), Delta, and Omicron SARS-CoV-2 pseudoviruses against sera samples drawn four months post administration of a third dose of the BNT162b2 mRNA vaccine. Our findings show that in comparison to hCoV-19 and Delta SARS-CoV-2, Omicron lineages BA.1 and BA.2 exhibit enhanced infectivity and a sharp decline in their sensitivity to vaccine-induced neutralizing antibodies. Interestingly, P681 mutations within the viral spike do not play a role in the neutralization potential or infectivity of SARS Cov-2 pseudoviruses carrying mutations in this position. The P681 residue however, dictates the ability of the spike protein to promote fusion and syncytia formation between infected cells. While spike from hCoV-19 (P681) and Omicron (H681) promote only modest cell fusion and formation of syncytia between cells that express the spike-protein, Delta spike (R681) displays enhanced fusogenic activity and promotes syncytia formation. Additional analysis shows that a single P681R mutation within the hCoV-19 spike, or H681R within the Omicron spike, restores fusion potential to similar levels observed for the Delta R681 spike. Conversely, R681P point mutation within the spike of Delta pseudovirus abolishes efficient fusion and syncytia formation. Our investigation also demonstrates that spike proteins from hCoV-19 and Delta SARS-CoV-2 are efficiently incorporated into viral particles relative to the spike of Omicron lineages. We conclude that the third dose of the Pfizer-BNT162b2 provides appreciable protection against the newly emerged Omicron sub-lineages. However, the neutralization sensitivity of these new variants is diminished relative to that of the hCoV-19 or Delta SARS-CoV-2. We further show that the P681 residue within spike dictates cell fusion and syncytia formation with no effects on the infectivity of the specific viral variant and on its sensitivity to vaccine-mediated neutralization.

Differences in airborne stability of SARS-CoV-2 variants of concern is impacted by alkalinity of surrogates of respiratory aerosol

The mechanistic factors hypothesized to be key drivers for the loss of infectivity of viruses in the aerosol phase often remain speculative. Using a next-generation bioaerosol technology, we report measurements of the aero-stability of several SARS-CoV-2 variants of concern in aerosol droplets of well-defined size and composition at high (90%) and low (40%) relative humidity (RH) upwards of 40 min. When compared with the ancestral virus, the infectivity of the Delta variant displayed different decay profiles. At low RH, a loss of viral infectivity of approximately 55% was observed over the initial 5 s for both variants. Regardless of RH and variant, greater than 95% of the viral infectivity was lost after 40 min of being aerosolized. Aero-stability of the variants correlate with their sensitivities to alkaline pH. Removal of all acidic vapours dramatically increased the rate of infectivity decay, with 90% loss after 2 min, while the addition of nitric acid vapour improved aero-stability. Similar aero-stability in droplets of artificial saliva and growth medium was observed. A model to predict loss of viral infectivity is proposed: at high RH, the high pH of exhaled aerosol drives viral infectivity loss; at low RH, high salt content limits the loss of viral infectivity.

Licorice-saponin A3 is a broad-spectrum inhibitor for COVID-19 by targeting viral spike and anti-inflammation

Currently, human health due to corona virus disease 2019 (COVID-19) pandemic has been seriously threatened. The coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) protein plays a crucial role in virus transmission and several S-based therapeutic approaches have been approved for the treatment of COVID-19. However, the efficacy is compromised by the SARS-CoV-2 evolvement and mutation. Here we report the SARS-CoV-2 S protein receptor-binding domain (RBD) inhibitor licorice-saponin A3 (A3) could widely inhibit RBD of SARS-CoV-2 variants, including Beta, Delta, and Omicron BA.1, XBB and BQ1.1. Furthermore, A3 could potently inhibit SARS-CoV-2 Omicron virus in Vero E6 cells, with EC₅₀ of 1.016 μM. The mechanism was related with binding with Y453 of RBD determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) analysis combined with quantum mechanics/molecular mechanics (QM/MM) simulations. Interestingly, phosphoproteomics analysis and multi fluorescent immunohistochemistry (mIHC) respectively indicated that A3 also inhibits host inflammation by directly modulating the JNK and p38 MAPK pathways and rebalancing the corresponding immune dysregulation. This work supports A3 as a promising broad-spectrum small molecule drug candidate for COVID-19.

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

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

For Table of Contents Use Only

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

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

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

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

S:D614G and S:H655Y are gateway mutations that act epistatically to promote SARS-CoV-2 variant fitness

SARS-CoV-2 variants bearing complex combinations of mutations that confer increased transmissibility, COVID-19 severity, and immune escape, were first detected after S:D614G had gone to fixation, and likely originated during persistent infection of immunocompromised hosts. To test the hypothesis that S:D614G facilitated emergence of such variants, S:D614G was reverted to the ancestral sequence in the context of sequential Spike sequences from an immunocompromised individual, and within each of the major SARS-CoV-2 variants of concern. In all cases, infectivity of the S:D614G revertants was severely compromised. The infectivity of atypical SARS-CoV-2 lineages that propagated in the absence of S:D614G was found to be dependent upon either S:Q613H or S:H655Y. Notably, Gamma and Omicron variants possess both S:D614G and S:H655Y, each of which contributed to infectivity of these variants. Among sarbecoviruses, S:Q613H, S:D614G, and S:H655Y are only detected in SARS-CoV-2, which is also distinguished by a polybasic S1/S2 cleavage site. Genetic and biochemical experiments here showed that S:Q613H, S:D614G, and S:H655Y each stabilize Spike on virions, and that they are dispensable in the absence of S1/S2 cleavage, consistent with selection of these mutations by the S1/S2 cleavage site. CryoEM revealed that either S:D614G or S:H655Y shift the Spike receptor binding domain (RBD) towards the open conformation required for ACE2-binding and therefore on pathway for infection. Consistent with this, an smFRET reporter for RBD conformation showed that both S:D614G and S:H655Y spontaneously adopt the conformation that ACE2 induces in the parental Spike. Data from these orthogonal experiments demonstrate that S:D614G and S:H655Y are convergent adaptations to the polybasic S1/S2 cleavage site which stabilize S1 on the virion in the open RBD conformation and act epistatically to promote the fitness of variants bearing complex combinations of clinically significant mutations.

An intranasal influenza virus-vectored vaccine prevents SARS-CoV-2 replication in respiratory tissues of mice and hamsters

Current available vaccines for COVID-19 are effective in reducing severe diseases and deaths caused by SARS-CoV-2 infection but less optimal in preventing infection. Next-generation vaccines which are able to induce mucosal immunity in the upper respiratory to prevent or reduce infections caused by highly transmissible variants of SARS-CoV-2 are urgently needed. We have developed an intranasal vaccine candidate based on a live attenuated influenza virus (LAIV) with a deleted NS1 gene that encodes cell surface expression of the receptor-binding-domain (RBD) of the SARS-CoV-2 spike protein, designated DelNS1-RBD4N-DAF. Immune responses and protection against virus challenge following intranasal administration of DelNS1-RBD4N-DAF vaccines were analyzed in mice and compared with intramuscular injection of the BioNTech BNT162b2 mRNA vaccine in hamsters. DelNS1-RBD4N-DAF LAIVs induced high levels of neutralizing antibodies against various SARS-CoV-2 variants in mice and hamsters and stimulated robust T cell responses in mice. Notably, vaccination with DelNS1-RBD4N-DAF LAIVs, but not BNT162b2 mRNA, prevented replication of SARS-CoV-2 variants, including Delta and Omicron BA.2, in the respiratory tissues of animals. The DelNS1-RBD4N-DAF LAIV system warrants further evaluation in humans for the control of SARS-CoV-2 transmission and, more significantly, for creating dual function vaccines against both influenza and COVID-19 for use in annual vaccination strategies.

SARS-CoV-2 before and after Omicron: two different viruses and two different diseases?

For the first time in the history of medicine, it has been possible to describe-after a spillover-the evolution of a new human virus spreading in a non-immune population. This allowed not only to observe the subsequent emersion of variants endowed with features providing the virus with an evolutionary advantage, but also the shift of the pathways of virus replication and the acquisition of immunoevasive features. These characteristics had a remarkable influence on the diffusion of the SARS-CoV-2 and on the clinical presentation and prognosis of COVID-19, aspects that are described and commented in this review.

In SARS-CoV-2 delta variants, Spike-P681R and D950N promote membrane fusion, Spike-P681R enhances spike cleavage, but neither substitution affects pathogenicity in hamsters

BACKGROUND

The SARS-CoV-2 delta (B.1.617.2 lineage) variant was first identified at the end of 2020 and possessed two unique amino acid substitutions in its spike protein: S-P681R, at the S1/S2 cleavage site, and S-D950N, in the HR1 of the S2 subunit. However, the roles of these substitutions in virus phenotypes have not been fully characterized.

METHODS

We used reverse genetics to generate Wuhan-D614G viruses with these substitutions and delta viruses lacking these substitutions and explored how these changes affected their viral characteristics in vitro and in vivo.

FINDINGS

S-P681R enhanced spike cleavage and membrane fusion, whereas S-D950N slightly promoted membrane fusion. Although S-681R reduced the virus replicative ability especially in VeroE6 cells, neither substitution affected virus replication in Calu-3 cells and hamsters. The pathogenicity of all recombinant viruses tested in hamsters was slightly but not significantly affected.

INTERPRETATION

Our observations suggest that the S-P681R and S-D950N substitutions alone do not increase virus pathogenicity, despite of their enhancement of spike cleavage or fusogenicity.

FUNDING

A full list of funding bodies that contributed to this study can be found under Acknowledgments.

Dual mechanism: Epigenetic inhibitor apabetalone reduces SARS-CoV-2 Delta and Omicron variant spike binding and attenuates SARS-CoV-2 RNA induced inflammation

The SARS-CoV-2 virus initiates infection via interactions between the viral spike protein and the ACE2 receptors on host cells. Variants of concern have mutations in the spike protein that enhance ACE2 binding affinity, leading to increased virulence and transmission. Viral RNAs released after entry into host cells trigger interferon-I (IFN-I) mediated inflammatory responses for viral clearance and resolution of infection. However, overreactive host IFN-I responses and pro-inflammatory signals drive COVID-19 pathophysiology and disease severity during acute infection. These immune abnormalities also lead to the development of post-COVID syndrome if persistent. Novel therapeutics are urgently required to reduce short- and long-term pathologic consequences associated with SARS-CoV-2 infection. Apabetalone, an inhibitor of epigenetic regulators of the BET protein family, is a candidate for COVID-19 treatment via a dual mechanism of action. In vitro, apabetalone downregulates ACE2 gene expression to limit SARS-CoV-2 entry and propagation. In pre-clinical models and patients treated for cardiovascular disease, apabetalone inhibits expression of inflammatory mediators involved in the pathologic cytokine storm (CS) stimulated by various cytokines. Here we show apabetalone treatment of human lung epithelial cells reduces binding of viral spike protein regardless of mutations found in the highly contagious Delta variant and heavily mutated Omicron. Additionally, we demonstrate that apabetalone counters expression of pro-inflammatory factors with roles in CS and IFN-I signaling in lung cells stimulated with SARS-CoV-2 RNA. Our results support clinical evaluation of apabetalone to treat COVID-19 and post-COVID syndrome regardless of the SARS-CoV-2 variant.

Comprehensive deep mutational scanning reveals the pH induced stability and binding differences between SARS-CoV-2 spike RBD and human ACE2

The SARS-CoV-2 spike (S) glycoprotein with its mobile receptor-binding domain (RBD), binds to the human ACE2 receptor and thus facilitates virus entry through low-pH-endosomal pathways. The high degree of SARS-CoV-2 mutability has raised concern among scientists and medical professionals because it created doubt about the effectiveness of drugs and vaccinations designed specifically for COVID-19. In this study, we used computational saturation mutagenesis approach, including structure-based free energy calculations to analyse the effects of the missense mutations on the SARS-CoV-2 S-RBD stability and the S-RBD binding affinity with ACE2 at three different pH (pH 4.5, pH 6.5, and pH 7.4). A total of 3705 mutations in the S-RBD protein were analyzed, and we discovered that most of these mutations destabilize the RBD protein. Specifically, residues G404, G431, G447, A475, and G526 were important for RBD protein stability. In addition, RBD residues Y449, Y489, Y495, Q498, and N487 were critical for the RBD-ACE2 interaction. Next, we found that the distribution of the mean stability changes and mean binding energy changes of RBD due to mutations at both serological and endosomal pH correlated well, indicating the similar effects of mutations. Overall, this computational analysis is useful for understanding the effects of missense mutations in SARS-CoV-2 pathogenesis at different pH.Communicated by Ramaswamy H. Sarma.

Heparan sulfate promotes ACE2 super-cluster assembly to enhance SARS-CoV-2-associated syncytium formation

The mechanism of syncytium formation, caused by spike-induced cell-cell fusion in severe COVID-19, is largely unclear. Here we combine chemical genetics with 4D confocal imaging to establish the cell surface heparan sulfate (HS) as a critical host factor exploited by SARS-CoV-2 to enhance spike’s fusogenic activity. HS binds spike to facilitate ACE2 clustering, generating synapse-like cell-cell contacts to promote fusion pore formation. ACE2 clustering, and thus, syncytium formation is significantly mitigated by chemical or genetic elimination of cell surface HS, while in a cell-free system consisting of purified HS, spike, and lipid-anchored ACE2, HS directly induces ACE2 clustering. Importantly, the interaction of HS with spike allosterically enables a conserved ACE2 linker in receptor clustering, which concentrates spike at the fusion site to overcome fusion-associated activity loss. This fusion-boosting mechanism can be effectively targeted by an investigational HS-binding drug, which reduces syncytium formation in vitro and viral infection in mice.

Research progress in spike mutations of SARS-CoV-2 variants and vaccine development

The coronavirus disease 2019 (COVID-19) pandemic can hardly end with the emergence of different variants over time. In the past 2 years, several variants of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), such as the Delta and Omicron variants, have emerged with higher transmissibility, immune evasion and drug resistance, leading to higher morbidity and mortality in the population. The prevalent variants of concern (VOCs) share several mutations on the spike that can affect virus characteristics, including transmissibility, antigenicity, and immune evasion. Increasing evidence has demonstrated that the neutralization capacity of sera from COVID-19 convalescent or vaccinated individuals is decreased against SARS-CoV-2 variants. Moreover, the vaccine effectiveness of current COVID-19 vaccines against SARS-CoV-2 VOCs is not as high as that against wild-type SARS-CoV-2. Therefore, more attention might be paid to how the mutations impact vaccine effectiveness. In this review, we summarized the current studies on the mutations of the SARS-CoV-2 spike, particularly of the receptor binding domain, to elaborate on how the mutations impact the infectivity, transmissibility and immune evasion of the virus. The effects of mutations in the SARS-CoV-2 spike on the current therapeutics were highlighted, and potential strategies for future vaccine development were suggested.

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

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

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

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

The impact of SARS-CoV-2 3CL protease mutations on nirmatrelvir inhibitory efficiency. Computational insights into potential resistance mechanisms

The use of antiviral drugs can promote the appearance of mutations in the target protein that increase the resistance of the virus to the treatment. This is also the case of nirmatrelvir, a covalent inhibitor of the 3CL protease, or main protease, of SARS-CoV-2. In this work we show how the by-residue decomposition of noncovalent interactions established between the drug and the enzyme, in combination with an analysis of naturally occurring mutations, can be used to detect potential mutations in the 3CL protease conferring resistance to nirmatrelvir. We also investigate the consequences of these mutations on the reaction mechanism to form the covalent enzyme-inhibitor complex using QM/MM methods. In particular, we show that the E166V variant of the protease displays smaller binding affinity to nirmatrelvir and larger activation free energy for the formation of the covalent complex, both factors contributing to the observed resistance to the treatment with this drug. The conclusions derived from our work can be used to anticipate the consequences of the introduction of nirmatrelvir in the fitness landscape of the virus and to design new inhibitors adapted to some of the possible resistance mechanisms.

Development of highly effective LCB1-based lipopeptides targeting the spike receptor-binding motif of SARS-CoV-2

LCB1 is a computationally designed 56-mer miniprotein targeting the spike (S) receptor-binding motif of SARS-CoV- 2 with high potent activity (Science, 2020; Cell host microbe, 2021); however, recent studies have demonstrated that emerging SARS-CoV-2 variants are highly resistant to LCB1's inhibition. In this study, we first identified a truncated peptide termed LCB1v8, which maintained the high antiviral potency. Then, a group of lipopeptides were generated by modifying LCB1v8 with diverse lipids, and of two lipopeptides, the C-terminally stearicacid-conjugtaed LCB1v17 and cholesterol-conjugated LCB1v18, were highly effective in inhibiting both S protein-pseudovirus and authentic SARS-CoV-2 infections. We further showed that LCB1-based inhibitors had similar α-helicity and thermostability in structure and bound to the target-mimic RBD protein with high affinity, and the lipopeptides exhibited greatly enhanced binding with the viral and cellular membranes, improved inhibitory activities against emerging SARS-CoV-2 variants. Moreover, LCB1v18 was validated with high preventive and therapeutic efficacies in K18-hACE2 transgenic mice against lethal SARS-CoV-2 challenge. In conclusion, our studies have provided important information for understanding the structure and activity relationship (SAR) of LCB1 inhibitor and would guide the future development of novel antivirals.

Drug Cocktail Formulation via Circuit Design

Electronic circuits intuitively visualize and quantitatively simulate biological systems with nonlinear differential equations that exhibit complicated dynamics. Drug cocktail therapies are a powerful tool against diseases that exhibit such dynamics. We show that just six key states, which are represented in a feedback circuit, enable drug-cocktail formulation: 1) healthy cell number; 2) infected cell number; 3) extracellular pathogen number; 4) intracellular pathogenic molecule number; 5) innate immune system strength; and 6) adaptive immune system strength. To enable drug cocktail formulation, the model represents the effects of the drugs in the circuit. For example, a nonlinear feedback circuit model fits measured clinical data, represents cytokine storm and adaptive autoimmune behavior, and accounts for age, sex, and variant effects for SARS-CoV-2 with few free parameters. The latter circuit model provided three quantitative insights on the optimal timing and dosage of drug components in a cocktail: 1) antipathogenic drugs should be given early in the infection, but immunosuppressant timing involves a tradeoff between controlling pathogen load and mitigating inflammation; 2) both within and across-class combinations of drugs have synergistic effects; 3) if they are administered sufficiently early in the infection, anti-pathogenic drugs are more effective at mitigating autoimmune behavior than immunosuppressant drugs.

Distinct in vitro and in vivo neutralization profiles of monoclonal antibodies elicited by the receptor binding domain of the ancestral SARS-CoV-2

Broadly neutralizing antibodies against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants are sought to curb coronavirus disease 2019 (COVID-19) infections. Here we produced and characterized a set of mouse monoclonal antibodies (mAbs) specific for the ancestral SARS-CoV-2 receptor binding domain (RBD). Two of them, 17A7 and 17B10, were highly potent in microneutralization assay with 50% inhibitory concentration (IC50 ) ≤135 ng/mL against infectious SARS-CoV-2 variants, including G614, Alpha, Beta, Gamma, Delta, Epsilon, Zeta, Kappa, Lambda, B.1.1.298, B.1.222, B.1.5, and R.1. Both mAbs (especially 17A7) also exhibited strong in vivo efficacy in protecting K18-hACE2 transgenic mice from the lethal infection with G614, Alpha, Beta, Gamma, and Delta viruses. Structural analysis indicated that 17A7 and 17B10 target the tip of the receptor binding motif in the RBD-up conformation. A third RBD-reactive mAb (3A6) although escaped by Beta and Gamma, was highly effective in cross-neutralizing Delta and Omicron BA.1 variants in vitro and in vivo. In competition experiments, antibodies targeting epitopes similar to these 3 mAbs were rarely enriched in human COVID-19 convalescent sera or postvaccination sera. These results are helpful to inform new antibody/vaccine design and these mAbs can be useful tools for characterizing SARS-CoV-2 variants and elicited antibody responses.

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.

Profiling humoral immunity after mixing and matching COVID-19 vaccines using SARS-CoV-2 variant protein microarrays

In November 2022, 68% of the population received at least one dose of COVID-19 vaccines. Due to the ongoing mutations, especially for the variants of concern (VOCs), it is important to monitor the humoral immune responses after different vaccination strategies. In this study, we developed a SARS-CoV-2 variant protein microarray that contained the spike proteins from the VOCs, e.g., alpha, beta, gamma, delta, and omicron, to quantify the binding antibody and surrogate neutralizing antibody. Plasmas were collected after two doses of matching AZD1222 (AZx2), two doses of matching mRNA-1273 (Mx2), or mixing AZD1222 and mRNA-1273 (AZ+M). The results showed a significant decrease of surrogate neutralizing antibodies against the receptor-binding domain in all VOCs in AZx2 and Mx2 but not AZ+M. A similar but minor reduction pattern of surrogate neutralizing antibodies against the extracellular domain was observed. While Mx2 exhibited a higher surrogate neutralizing level against all VOCs compared to AZx2, AZ+M showed an even higher surrogate neutralizing level in gamma and omicron compared to Mx2. It is worth noting that the binding antibody displayed a low correlation to the surrogate neutralizing antibody (R-square 0.130-0.382). This study delivers insights into humoral immunities, SARS-CoV-2 mutations, and mixing and matching vaccine strategies, which may provide a more effective vaccine strategy especially in preventing omicron.

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid-amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

Characterization of a Vesicular Stomatitis Virus-Vectored Recombinant Virus Bearing Spike Protein of SARS-CoV-2 Delta Variant

The frequent emergence of SARS-CoV-2 variants thwarts the prophylactic and therapeutic countermeasures confronting COVID-19. Among them, the Delta variant attracts widespread attention due to its high pathogenicity and fatality rate compared with other variants. However, with the emergence of new variants, studies on Delta variants have been gradually weakened and ignored. In this study, a replication-competent recombinant virus carrying the S protein of the SARS-CoV-2 Delta variant was established based on the vesicular stomatitis virus (VSV), which presented a safe alternative model for studying the Delta variant. The recombinant virus showed a replication advantage in Vero E6 cells, and the viral titers reach 10^7.3 TCID₅₀/mL at 36 h post-inoculation. In the VSV-vectored recombinant platform, the spike proteins of the Delta variant mediated higher fusion activity and syncytium formation than the wild-type strain. Notably, the recombinant virus was avirulent in BALB/c mice, Syrian hamsters, 3-day ICR suckling mice, and IFNAR/GR^-/- mice. It induced protective neutralizing antibodies in rodents, and protected the Syrian hamsters against the SARS-CoV-2 Delta variant infection. Meanwhile, the eGFP reporter of recombinant virus enabled the visual assay of neutralizing antibodies. Therefore, the recombinant virus could be a safe and convenient surrogate tool for authentic SARS-CoV-2. This efficient and reliable model has significant potential for research on viral-host interactions, epidemiological investigation of serum-neutralizing antibodies, and vaccine development.

Chiral-Selective Antigen-Presentation by Supramolecular Chiral Polymer Micelles

Chirality is ubiquitous in biological systems, which is closely related to biological functions, life processes, and even the pathogenesis of diseases. However, the interface between the chirality of synthetic materials and organisms, particularly the immune system, remains poorly understood. Here, supramolecular chiral polymer micelles (SCPMs) are prepared by complexing antigenic proteins with chiral amino acid-modified polyethyleneimine. The introduction of chirality not only reduces the toxicity of cationic polymer, but also benefits cell uptake and antigen presentation. Especially, D-chirality presents the lowest cytotoxicity, while promoting the highest expression level of costimulatory molecules on dendritic cells compared to L-chirality and achirality. The superiority of D-chirality to stimulate dendritic cell maturation is supported by immunization with D-SCPMs, which achieves significant antigen-specific proliferation of T cells in the spleen, lymph nodes, and tumor of mice. Chirality-mediated antigen processing and presentation are demonstrated by D-SCPMs self-assembled from chiral alkaline histidine or neutral phenylalanine modified polyethyleneimine and tumor associated ovalbumin or severe acute respiratory syndrome coronavirus 2 spike 1 antigenic protein. Immunoactivation enabled by D-chirality opens a window to prepare potent nanotherapeutics for disease prevention and treatment.

Delta variant (B.1.617.2) of SARS-CoV-2: current understanding of infection, transmission, immune escape, and mutational landscape

The Delta variant is one of the alarming variants of the SARS-CoV-2 virus that have been immensely detrimental and a significant cause of the prolonged pandemic (B.1.617.2). During the SARS-CoV-2 pandemic from December 2020 to October 2021, the Delta variant showed global dominance, and afterwards, the Omicron variant showed global dominance. Delta shows high infectivity rate which accounted for nearly 70% of the cases after December 2020. This review discusses the additional attributes that make the Delta variant so infectious and transmissible. The study also focuses on the significant mutations, namely the L452R and T478K present on the receptor-binding domain of spike (S)-glycoprotein, which confers specific alterations to the Delta variant. Considerably, we have also highlighted other notable factors such as the immune escape, infectivity and re-infectivity, vaccine escape, Ro number, S-glycoprotein stability, cleavage pattern, and its binding affinity with the host cell receptor protein. We have also emphasized clinical manifestations, symptomatology, morbidity, and mortality for the Delta variant compared with other significant SARS-CoV-2 variants. This review will help the researchers to get an elucidative view of the Delta variant to adopt some practical strategies to minimize the escalating spread of the SARS-CoV-2 Delta variant.

Structural analysis of receptor engagement and antigenic drift within the BA.2 spike protein

The BA.2 sub-lineage of the Omicron (B.1.1.529) severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant rapidly supplanted the original BA.1 sub-lineage in early 2022. Both lineages threatened the efficacy of vaccine-elicited antibodies and acquired increased binding to several mammalian ACE2 receptors. Cryoelectron microscopy (cryo-EM) analysis of the BA.2 spike (S) glycoprotein in complex with mouse ACE2 (mACE2) identifies BA.1- and BA.2-mutated residues Q493R, N501Y, and Y505H as complementing non-conserved residues between human and mouse ACE2, rationalizing the enhanced S protein-mACE2 interaction for Omicron variants. Cryo-EM structures of the BA.2 S-human ACE2 complex and of the extensively mutated BA.2 amino-terminal domain (NTD) reveal a dramatic reorganization of the highly antigenic N1 loop into a β-strand, providing an explanation for decreased binding of the BA.2 S protein to antibodies isolated from BA.1-convalescent patients. Our analysis reveals structural mechanisms underlying the antigenic drift in the rapidly evolving Omicron variant landscape.

Defending against SARS-CoV-2: The T cell perspective

SARS-CoV-2-specific T cell response has been proven essential for viral clearance, COVID-19 outcome and long-term memory. Impaired early T cell-driven immunity leads to a severe form of the disease associated with lymphopenia, hyperinflammation and imbalanced humoral response. Analyses of acute SARS-CoV-2 infection have revealed that mild COVID-19 course is characterized by an early induction of specific T cells within the first 7 days of symptoms, coordinately followed by antibody production for an effective control of viral infection. In contrast, patients who do not develop an early specific cellular response and initiate a humoral immune response with subsequent production of high levels of antibodies, develop severe symptoms. Yet, delayed and persistent bystander CD8+ T cell activation has been also reported in hospitalized patients and could be a driver of lung pathology. Literature supports that long-term maintenance of T cell response appears more stable than antibody titters. Up to date, virus-specific T cell memory has been detected 22 months post-symptom onset, with a predominant IL-2 memory response compared to IFN-γ. Furthermore, T cell responses are conserved against the emerging variants of concern (VoCs) while these variants are mostly able to evade humoral responses. This could be partly explained by the high HLA polymorphism whereby the viral epitope repertoire recognized could differ among individuals, greatly decreasing the likelihood of immune escape. Current COVID-19-vaccination has been shown to elicit Th1-driven spike-specific T cell response, as does natural infection, which provides substantial protection against severe COVID-19 and death. In addition, mucosal vaccination has been reported to induce strong adaptive responses both locally and systemically and to protect against VoCs in animal models. The optimization of vaccine formulations by including a variety of viral regions, innovative adjuvants or diverse administration routes could result in a desirable enhanced cellular response and memory, and help to prevent breakthrough infections. In summary, the increasing evidence highlights the relevance of monitoring SARS-CoV-2-specific cellular immune response, and not only antibody levels, as a correlate for protection after infection and/or vaccination. Moreover, it may help to better identify target populations that could benefit most from booster doses and to personalize vaccination strategies.

Guedes J, Nogueira FC, Silva JL, Castilho LR, de Oliveira GA. Conformational stability of SARS-CoV-2 glycoprotein spike variants

The severe acute respiratory syndrome spread worldwide, causing a pandemic. SARS-CoV-2 mutations have arisen in the spike, a glycoprotein at the viral envelope and an antigenic candidate for vaccines against COVID-19. Here, we present comparative data of the glycosylated full-length ancestral and D614G spike together with three other transmissible strains classified by the World Health Organization as variants of concern: beta, gamma, and delta. By showing that D614G has less hydrophobic surface exposure and trimer persistence, we place D614G with features that support a model of temporary fitness advantage for virus spillover. Furthermore, during the SARS-CoV-2 adaptation, the spike accumulates alterations leading to less structural stability for some variants. The decreased trimer stability of the ancestral and gamma and the presence of D614G uncoupled conformations mean higher ACE-2 affinities compared to the beta and delta strains. Mapping the energetics and flexibility of variants is necessary to improve vaccine development.

The Relationship between the Transmission of Different SARS-CoV-2 Strains and Air Quality: A Case Study in China

Coronavirus Disease 2019 (COVID-19) has been a global public health concern for almost three years, and the transmission characteristics vary among different virus variants. Previous studies have investigated the relationship between air pollutants and COVID-19 infection caused by the original strain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, it is unclear whether individuals might be more susceptible to COVID-19 due to exposure to air pollutants, with the SARS-CoV-2 mutating faster and faster. This study aimed to explore the relationship between air pollutants and COVID-19 infection caused by three major SARS-CoV-2 strains (the original strain, Delta variant, and Omicron variant) in China. A generalized additive model was applied to investigate the associations of COVID-19 infection with six air pollutants (PM_2.5, PM₁₀, SO₂, CO, NO₂, and O₃). A positive correlation might be indicated between air pollutants (PM_2.5, PM₁₀, and NO₂) and confirmed cases of COVID-19 caused by different SARS-CoV-2 strains. It also suggested that the mutant variants appear to be more closely associated with air pollutants than the original strain. This study could provide valuable insight into control strategies that limit the concentration of air pollutants at lower levels and would better control the spread of COVID-19 even as the virus continues to mutate.

Rapid transmission and tight bottlenecks constrain the evolution of highly transmissible SARS-CoV-2 variants

Transmission bottlenecks limit the spread of novel mutations and reduce the efficiency of selection along a transmission chain. While increased force of infection, receptor binding, or immune evasion may influence bottleneck size, the relationship between transmissibility and the transmission bottleneck is unclear. Here we compare the transmission bottleneck of non-VOC SARS-CoV-2 lineages to those of Alpha, Delta, and Omicron. We sequenced viruses from 168 individuals in 65 households. Most virus populations had 0-1 single nucleotide variants (iSNV). From 64 transmission pairs with detectable iSNV, we identify a per clade bottleneck of 1 (95% CI 1-1) for Alpha, Delta, and Omicron and 2 (95% CI 2-2) for non-VOC. These tight bottlenecks reflect the low diversity at the time of transmission, which may be more pronounced in rapidly transmissible variants. Tight bottlenecks will limit the development of highly mutated VOC in transmission chains, adding to the evidence that selection over prolonged infections may drive their evolution.

The Delta variant wave in Tunisia: Genetic diversity, spatio-temporal distribution and evidence of the spread of a divergent AY.122 sub-lineage

Introduction

The Delta variant posed an increased risk to global public health and rapidly replaced the pre-existent variants worldwide. In this study, the genetic diversity and the spatio-temporal dynamics of 662 SARS-CoV2 genomes obtained during the Delta wave across Tunisia were investigated.

Methods

Viral whole genome and partial S-segment sequencing was performed using Illumina and Sanger platforms, respectively and lineage assignemnt was assessed using Pangolin version 1.2.4 and scorpio version 3.4.X. Phylogenetic and phylogeographic analyses were achieved using IQ-Tree and Beast programs.

Results

The age distribution of the infected cases showed a large peak between 25 to 50 years. Twelve Delta sub-lineages were detected nation-wide with AY.122 being the predominant variant representing 94.6% of sequences. AY.122 sequences were highly related and shared the amino-acid change ORF1a:A498V, the synonymous mutations 2746T>C, 3037C>T, 8986C>T, 11332A>G in ORF1a and 23683C>T in the S gene with respect to the Wuhan reference genome (NC_045512.2). Spatio-temporal analysis indicates that the larger cities of Nabeul, Tunis and Kairouan constituted epicenters for the AY.122 sub-lineage and subsequent dispersion to the rest of the country.

Discussion

This study adds more knowledge about the Delta variant and sub-variants distribution worldwide by documenting genomic and epidemiological data from Tunisia, a North African region. Such results may be helpful to the understanding of future COVID-19 waves and variants.

Structural characterization of SARS-CoV-2 dimeric ORF9b reveals potential fold-switching trigger mechanism

The constant emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants indicates the evolution and adaptation of the virus. Enhanced innate immune evasion through increased expression of viral antagonist proteins, including ORF9b, contributes to the improved transmission of the Alpha variant; hence, more attention should be paid to these viral proteins. ORF9b is an accessory protein that suppresses innate immunity via a monomer conformation by binding to Tom70. Here, we solved the dimeric structure of SARS-CoV-2 ORF9b with a long hydrophobic tunnel containing a lipid molecule that is crucial for the dimeric conformation and determined the specific lipid ligands as monoglycerides by conducting a liquid chromatography with tandem mass spectrometry analysis, suggesting an important role in the viral life cycle. Notably, a long intertwined loop accessible for host factor binding was observed in the structure. Eight phosphorylated residues in ORF9b were identified, and residues S50 and S53 were found to contribute to the stabilization of dimeric ORF9b. Additionally, we proposed a model of multifunctional ORF9b with a distinct conformation, suggesting that ORF9b is a fold-switching protein, while both lipids and phosphorylation contribute to the switching. Specifically, the ORF9b monomer interacts with Tom70 to suppress the innate immune response, whereas the ORF9b dimer binds to the membrane involving mature virion assembly. Our results provide a better understanding of the multiple functions of ORF9b.

Towards an optimal monoclonal antibody with higher binding affinity to the receptor-binding domain of SARS-CoV-2 spike proteins from different variants

A highly efficient and robust multiple scales in silico protocol, consisting of atomistic Molecular Dynamics (MD), coarse-grain (CG) MD, and constant-pH CG Monte Carlo (MC), has been developed and used to study the binding affinities of selected antigen-binding fragments of the monoclonal antibody (mAbs) CR3022 and several of its here optimized versions against 11 SARS-CoV-2 variants including the wild type. Totally 235,000 mAbs structures were initially generated using the RosettaAntibodyDesign software, resulting in top 10 scored CR3022-like-RBD complexes with critical mutations and compared to the native one, all having the potential to block virus-host cell interaction. Of these 10 finalists, two candidates were further identified in the CG simulations to be the best against all SARS-CoV-2 variants. Surprisingly, all 10 candidates and the native CR3022 exhibited a higher affinity for the Omicron variant despite its highest number of mutations. The multiscale protocol gives us a powerful rational tool to design efficient mAbs. The electrostatic interactions play a crucial role and appear to be controlling the affinity and complex building. Studied mAbs carrying a more negative total net charge show a higher affinity. Structural determinants could be identified in atomistic simulations and their roles are discussed in detail to further hint at a strategy for designing the best RBD binder. Although the SARS-CoV-2 was specifically targeted in this work, our approach is generally suitable for many diseases and viral and bacterial pathogens, leukemia, cancer, multiple sclerosis, rheumatoid, arthritis, lupus, and more.

Rational Development of Hypervalent Glycan Shield-Binding Nanoparticles with Broad-Spectrum Inhibition against Fatal Viruses Including SARS-CoV-2 Variants

Infectious virus diseases, particularly coronavirus disease 2019, have posed a severe threat to public health, whereas the developed therapeutic and prophylactic strategies are seriously challenged by viral evolution and mutation. Therefore, broad-spectrum inhibitors of viruses are highly demanded. Herein, an unprecedented antiviral strategy is reported, targeting the viral glycan shields with hypervalent mannose-binding nanoparticles. The nanoparticles exhibit a unique double-punch mechanism, being capable of not only blocking the virus-receptor interaction but also inducing viral aggregation, thereby allowing for inhibiting the virus entry and facilitating the phagocytosis of viruses. The nanoparticles exhibit potent and broad-spectrum antiviral efficacy to multiple pseudoviruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its major variants (D614G, N501Y, N439K, Δ69-70, Delta, and Omicron; lentiviruses expressing only the spike proteins), as well as other vital viruses (human immunodeficiency virus 1 and Lassa virus), with apparent EC50 values around the 10^-9  m level. Significantly, the broad-spectrum inhibition of authentic viruses of both wild-type SARS-CoV-2 and Delta variants is confirmed. Therefore, this hypervalent glycan-shield targeting strategy opens new access to broad-spectrum viral inhibition.

Broadly neutralizing and protective nanobodies against SARS-CoV-2 Omicron subvariants BA.1, BA.2, and BA.4/5 and diverse sarbecoviruses

As SARS-CoV-2 Omicron and other variants of concern (VOCs) continue spreading worldwide, development of antibodies and vaccines to confer broad and protective activity is a global priority. Here, we report on the identification of a special group of nanobodies from immunized alpaca with potency against diverse VOCs including Omicron subvariants BA.1, BA.2 and BA.4/5, SARS-CoV-1, and major sarbecoviruses. Crystal structure analysis of one representative nanobody, 3-2A2-4, discovers a highly conserved epitope located between the cryptic and the outer face of the receptor binding domain (RBD), distinctive from the receptor ACE2 binding site. Cryo-EM and biochemical evaluation reveal that 3-2A2-4 interferes structural alteration of RBD required for ACE2 binding. Passive delivery of 3-2A2-4 protects K18-hACE2 mice from infection of authentic SARS-CoV-2 Delta and Omicron. Identification of these unique nanobodies will inform the development of next generation antibody therapies and design of pan-sarbecovirus vaccines.

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.

Enhanced virulence and waning vaccine-elicited antibodies account for breakthrough infections caused by SARS-CoV-2 delta and beyond

Here we interrogate the factors responsible for SARS-CoV-2 breakthrough infections in a K18-hACE2 transgenic mouse model. We show that Delta and the closely related Kappa variant cause viral pneumonia and severe lung lesions in K18-hACE2 mice. Human COVID-19 mRNA post-vaccination sera after the 2nd dose are significantly less efficient in neutralizing Delta/Kappa than early 614G virus in vitro and in vivo. By 5 months post-vaccination, ≥50% of donors lack detectable neutralizing antibodies against Delta and Kappa and all mice receiving 5-month post-vaccination sera die after the lethal challenges. Although a 3rd vaccine dose can boost antibody neutralization against Delta in vitro and in vivo, the mean log neutralization titers against the latest Omicron subvariants are 1/3-1/2 of those against the original 614D virus. Our results suggest that enhanced virulence, greater immune evasion, and waning of vaccine-elicited protection account for SARS-CoV-2 variants caused breakthrough infections.