A thermostable, closed SARS-CoV-2 spike protein trimer

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

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

Humoral and cellular response of two different vaccines against SARS-CoV-2 in a group of healthcare workers: An observational study

On March 13, 2021, Tunisia started a widespread immunization program against SARS-CoV-2 utilizing different vaccinations that had been given emergency approval. Herein, we followed prospectively a cohort of participant who received COVID-19 vaccine (Pfizer BioNTech and Sputnik-Gameleya V). The goal of this follow-up was to define the humoral and cellular immunological profile after immunization by assessing neutralizing antibodies and IFN- γ release. 26 vaccinated health care workers by Pfizer BioNTech (n=12) and Sputnik-Gameleya V (n=14) were enrolled from June to December 2021 in Military hospital of Tunis. All consenting participants were sampled for peripheral blood after three weeks of vaccination. The humoral response was investigated by the titer of anti-SARS-CoV-2 immunoglobulin G (IgG) antibodies to S1 protein. The CD4 and CD8 T cell responses were evaluated by the QuantiFERON® SARS-CoV-2 (Qiagen® Basel, Switzerland). Regardless the type of vaccine, the assessment of humoral and cellular response following vaccination showed a strong involvement of the later with expression of IFN-γ as compared to antibodies secretion. Moreover, we showed that people with past SARS-CoV-2 infection developed high levels of antibodies than those who are not previously infected. However, no significant difference was detected concerning interferon gamma (IFN-γ) expression by CD4 and CD8 T cells in health care worker (HCW) previously infection or not with COVID-19 infection. Analysis of immune response according to the type of vaccine, we found that Pfizer BioNTech induced high level of humoral response (91.66%) followed by Sputnik-Gameleya V (64.28%). However, adenovirus vaccine gave a better cellular response (57.14%) than mRNA vaccine (41.66%). Regarding the immune response following vaccine doses, we revealed a significant increase of neutralizing antibodies and IFN-γ release by T cells in patients fully vaccinated as compared to those who have received just one vaccine. Collectively, our data revealed a similar immune response between Pfizer BioNTech and Sputnik-Gameleya V vaccine with a slight increase of humoral response by mRNA vaccine and cellular response by adenovirus vaccine. It's evident that past SARS-CoV-2 infection was a factor that contributed to the vaccination's increased immunogenicity. However, the administration of full doses of vaccines (Pfizer BioNTech or Sputnik-Gameleya V) induces better humoral and cellular responses detectable even more than three months following vaccination.

Exploring Binding Pockets in the Conformational States of the SARS-CoV-2 Spike Trimers for the Screening of Allosteric Inhibitors Using Molecular Simulations and Ensemble-Based Ligand Docking

Understanding mechanisms of allosteric regulation remains elusive for the SARS-CoV-2 spike protein, despite the increasing interest and effort in discovering allosteric inhibitors of the viral activity and interactions with the host receptor ACE2. The challenges of discovering allosteric modulators of the SARS-CoV-2 spike proteins are associated with the diversity of cryptic allosteric sites and complex molecular mechanisms that can be employed by allosteric ligands, including the alteration of the conformational equilibrium of spike protein and preferential stabilization of specific functional states. In the current study, we combine conformational dynamics analysis of distinct forms of the full-length spike protein trimers and machine-learning-based binding pocket detection with the ensemble-based ligand docking and binding free energy analysis to characterize the potential allosteric binding sites and determine structural and energetic determinants of allosteric inhibition for a series of experimentally validated allosteric molecules. The results demonstrate a good agreement between computational and experimental binding affinities, providing support to the predicted binding modes and suggesting key interactions formed by the allosteric ligands to elicit the experimentally observed inhibition. We establish structural and energetic determinants of allosteric binding for the experimentally known allosteric molecules, indicating a potential mechanism of allosteric modulation by targeting the hinges of the inter-protomer movements and blocking conformational changes between the closed and open spike trimer forms. The results of this study demonstrate that combining ensemble-based ligand docking with conformational states of spike protein and rigorous binding energy analysis enables robust characterization of the ligand binding modes, the identification of allosteric binding hotspots, and the prediction of binding affinities for validated allosteric modulators, which is consistent with the experimental data. This study suggested that the conformational adaptability of the protein allosteric sites and the diversity of ligand bound conformations are both in play to enable efficient targeting of allosteric binding sites and interfere with the conformational changes.

Single-Molecule Investigation of the Binding Interface Stability of SARS-CoV-2 Variants with ACE2

The SARS-CoV-2 pandemic spurred numerous research endeavors to comprehend the virus and mitigate its global severity. Understanding the binding interface between the virus and human receptors is pivotal to these efforts and paramount to curbing infection and transmission. Here we employ atomic force microscopy and steered molecular dynamics simulation to explore SARS-CoV-2 receptor binding domain (RBD) variants and angiotensin-converting enzyme 2 (ACE2), examining the impact of mutations at key residues upon binding affinity. Our results show that the Omicron and Delta variants possess strengthened binding affinity in comparison to the Mu variant. Further, using sera from individuals either vaccinated or with acquired immunity following Delta strain infection, we assess the impact of immunity upon variant RBD/ACE2 complex formation. Single-molecule force spectroscopy analysis suggests that vaccination before infection may provide stronger protection across variants. These results underscore the need to monitor antigenic changes in order to continue developing innovative and effective SARS-CoV-2 abrogation strategies.

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

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

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

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

Characterization of a neutralizing antibody that recognizes a loop region adjacent to the receptor-binding interface of the SARS-CoV-2 spike receptor-binding domain

Although the global crisis caused by the coronavirus disease 2019 (COVID-19) pandemic is over, the global epidemic of the disease continues. Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the cause of COVID-19, initiates infection via the binding of the receptor-binding domain (RBD) of its spike protein to the human angiotensin-converting enzyme II (ACE2) receptor, and this interaction has been the primary target for the development of COVID-19 therapeutics. Here, we identified neutralizing antibodies against SARS-CoV-2 by screening mouse monoclonal antibodies and characterized an antibody, CSW1-1805, that targets a narrow region at the RBD ridge of the spike protein. CSW1-1805 neutralized several variants in vitro and completely protected mice from SARS-CoV-2 infection. Cryo-EM and biochemical analyses revealed that this antibody recognizes the loop region adjacent to the ACE2-binding interface with the RBD in both a receptor-inaccessible "down" state and a receptor-accessible "up" state and could stabilize the RBD conformation in the up-state. CSW1-1805 also showed different binding orientations and complementarity determining region properties compared to other RBD ridge-targeting antibodies with similar binding epitopes. It is important to continuously characterize neutralizing antibodies to address new variants that continue to emerge. Our characterization of this antibody that recognizes the RBD ridge of the spike protein will aid in the development of future neutralizing antibodies.IMPORTANCESARS-CoV-2 cell entry is initiated by the interaction of the viral spike protein with the host cell receptor. Therefore, mechanistic findings regarding receptor recognition by the spike protein help uncover the molecular mechanism of SARS-CoV-2 infection and guide neutralizing antibody development. Here, we characterized a SARS-CoV-2 neutralizing antibody that recognizes an epitope, a loop region adjacent to the receptor-binding interface, that may be involved in the conformational transition of the receptor-binding domain (RBD) of the spike protein from a receptor-inaccessible "down" state into a receptor-accessible "up" state, and also stabilizes the RBD in the up-state. Our mechanistic findings provide new insights into SARS-CoV-2 receptor recognition and guidance for neutralizing antibody development.

Stability of gut microbiome after COVID-19 vaccination in healthy and immuno-compromised individuals

Bidirectional interactions between the immune system and the gut microbiota are key contributors to various physiological functions. Immune-associated diseases such as cancer and autoimmunity, and efficacy of immunomodulatory therapies, have been linked to microbiome variation. Although COVID-19 infection has been shown to cause microbial dysbiosis, it remains understudied whether the inflammatory response associated with vaccination also impacts the microbiota. Here, we investigate the temporal impact of COVID-19 vaccination on the gut microbiome in healthy and immuno-compromised individuals; the latter included patients with primary immunodeficiency and cancer patients on immunomodulating therapies. We find that the gut microbiome remained remarkably stable post-vaccination irrespective of diverse immune status, vaccine response, and microbial composition spanned by the cohort. The stability is evident at all evaluated levels including diversity, phylum, species, and functional capacity. Our results indicate the resilience of the gut microbiome to host immune changes triggered by COVID-19 vaccination and suggest minimal, if any, impact on microbiome-mediated processes. These findings encourage vaccine acceptance, particularly when contrasted with the significant microbiome shifts observed during COVID-19 infection.

Immunization with V987H-stabilized Spike glycoprotein protects K18-hACE2 mice and golden Syrian hamsters upon SARS-CoV-2 infection

Safe and effective severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccines are crucial to fight against the coronavirus disease 2019 pandemic. Most vaccines are based on a mutated version of the Spike glycoprotein [K986P/V987P (S-2P)] with improved stability, yield and immunogenicity. However, S-2P is still produced at low levels. Here, we describe the V987H mutation that increases by two-fold the production of the recombinant Spike and the exposure of the receptor binding domain (RBD). S-V987H immunogenicity is similar to S-2P in mice and golden Syrian hamsters (GSH), and superior to a monomeric RBD. S-V987H immunization confer full protection against severe disease in K18-hACE2 mice and GSH upon SARS-CoV-2 challenge (D614G or B.1.351 variants). Furthermore, S-V987H immunized K18-hACE2 mice show a faster tissue viral clearance than RBD- or S-2P-vaccinated animals challenged with D614G, B.1.351 or Omicron BQ1.1 variants. Thus, S-V987H protein might be considered for future SARS-CoV-2 vaccines development.

Exploring cell-free assays for COVID-19 serosurvey

Serosurveys to monitor immunity toward COVID-19 in the population are primarily performed using an ELISA to screen samples for SARS-CoV-2 antibodies, followed by confirmation by a virus neutralization test, which is considered the Gold Standard. However, virus neutralization test may not be feasible for some laboratories because of the requirement for specific facilities and trained personnel. In an attempt to address this limitation, we evaluated three cell-free methods as potential alternatives for assessing SARS-CoV-2 seroprevalence in human population from plasma. We report the establishment of two inhibition ELISAs designed to detect anti-Spike RBD IgG antibodies and a microsphere quantitative suspension array technology assay, based on the Luminex xMAP platform, to measure the presence of antibodies against various SARS-CoV-2 antigens, including anti-RBD. These methods were also compared to a commercial chemiluminescent immunoassay designed for anti-RBD antibodies detection and to the combined ELISA + virus neutralization test strategy. These cell-free assays performed equally to estimate the percentage of positive and negative samples and could be used to determine the prevalence of SARS-CoV-2 antibodies in human population, at least in cohort with high-expected prevalence, without the use of seroneutralization assay.

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

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

High-throughput screening of genetic and cellular drivers of syncytium formation induced by the spike protein of SARS-CoV-2

Mapping mutations and discovering cellular determinants that cause the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to induce infected cells to form syncytia would facilitate the development of strategies for blocking the formation of such cell-cell fusion. Here we describe high-throughput screening methods based on droplet microfluidics and the size-exclusion selection of syncytia, coupled with large-scale mutagenesis and genome-wide knockout screening via clustered regularly interspaced short palindromic repeats (CRISPR), for the large-scale identification of determinants of cell-cell fusion. We used the methods to perform deep mutational scans in spike-presenting cells to pinpoint mutable syncytium-enhancing substitutions in two regions of the spike protein (the fusion peptide proximal region and the furin-cleavage site). We also used a genome-wide CRISPR screen in cells expressing the receptor angiotensin-converting enzyme 2 to identify inhibitors of clathrin-mediated endocytosis that impede syncytium formation, which we validated in hamsters infected with SARS-CoV-2. Finding genetic and cellular determinants of the formation of syncytia may reveal insights into the physiological and pathological consequences of cell-cell fusion.

Pseudo-Luciferase Activity of the SARS-CoV-2 Spike Protein for Cypridina Luciferin

Enzymatic reactions that involve a luminescent substrate (luciferin) and enzyme (luciferase) from luminous organisms enable a luminescence detection of target proteins and cells with high specificity, albeit that conventional assay design requires a prelabeling of target molecules with luciferase. Here, we report a luciferase-independent luminescence assay in which the target protein directly catalyzes the oxidative luminescence reaction of luciferin. The SARS-CoV-2 antigen (spike) protein catalyzes the light emission of Cypridina luciferin, whereas no such catalytic function was observed for salivary proteins. This selective luminescence reaction is due to the enzymatic recognition of the 3-(1-guanidino)propyl group in luciferin at the interfaces between the units of the spike protein, allowing a specific detection of the spike protein in human saliva without sample pretreatment. This method offers a novel platform to detect virus antigens simply and rapidly without genetic manipulation or antibodies.

Spontaneous, persistent, T cell-dependent IFN-γ release in patients who progress to Long Covid

After acute infection with severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), a proportion of patients experience persistent symptoms beyond 12 weeks, termed Long Covid. Understanding the mechanisms that cause this debilitating disease and identifying biomarkers for diagnostic, therapeutic, and monitoring purposes are urgently required. We detected persistently high levels of interferon-γ (IFN-γ) from peripheral blood mononuclear cells of patients with Long Covid using highly sensitive FluoroSpot assays. This IFN-γ release was seen in the absence of ex vivo peptide stimulation and remains persistently elevated in patients with Long Covid, unlike the resolution seen in patients recovering from acute SARS-CoV-2 infection. The IFN-γ release was CD8+ T cell-mediated and dependent on antigen presentation by CD14+ cells. Longitudinal follow-up of our study cohort showed that symptom improvement and resolution correlated with a decrease in IFN-γ production to baseline levels. Our study highlights a potential mechanism underlying Long Covid, enabling the search for biomarkers and therapeutics in patients with Long Covid.

Microfluidic antibody profiling after repeated SARS-CoV-2 vaccination links antibody affinity and concentration to impaired immunity and variant escape in patients on anti-CD20 therapy

Background

Patients with autoimmune/inflammatory conditions on anti-CD20 therapies, such as rituximab, have suboptimal humoral responses to vaccination and are vulnerable to poorer clinical outcomes following SARS-CoV-2 infection. We aimed to examine how the fundamental parameters of antibody responses, namely, affinity and concentration, shape the quality of humoral immunity after vaccination in these patients.

Methods

We performed in-depth antibody characterisation in sera collected 4 to 6 weeks after each of three vaccine doses to wild-type (WT) SARS-CoV-2 in rituximab-treated primary vasculitis patients (n = 14) using Luminex and pseudovirus neutralisation assays, whereas we used a novel microfluidic-based immunoassay to quantify polyclonal antibody affinity and concentration against both WT and Omicron (B.1.1.529) variants. We performed comparative antibody profiling at equivalent timepoints in healthy individuals after three antigenic exposures to WT SARS-CoV-2 (one infection and two vaccinations; n = 15) and in convalescent patients after WT SARS-CoV-2 infection (n = 30).

Results

Rituximab-treated patients had lower antibody levels and neutralisation titres against both WT and Omicron SARS-CoV-2 variants compared to healthy individuals. Neutralisation capacity was weaker against Omicron versus WT both in rituximab-treated patients and in healthy individuals. In the rituximab cohort, this was driven by lower antibody affinity against Omicron versus WT [median (range) KD: 21.6 (9.7-38.8) nM vs. 4.6 (2.3-44.8) nM, p = 0.0004]. By contrast, healthy individuals with hybrid immunity produced a broader antibody response, a subset of which recognised Omicron with higher affinity than antibodies in rituximab-treated patients [median (range) KD: 1.05 (0.45-1.84) nM vs. 20.25 (13.2-38.8) nM, p = 0.0002], underpinning the stronger serum neutralisation capacity against Omicron in the former group. Rituximab-treated patients had similar anti-WT antibody levels and neutralisation titres to unvaccinated convalescent individuals, despite two more exposures to SARS-CoV-2 antigen. Temporal profiling of the antibody response showed evidence of affinity maturation in healthy convalescent patients after a single SARS-CoV-2 infection, which was not observed in rituximab-treated patients, despite repeated vaccination.

Discussion

Our results enrich previous observations of impaired humoral immune responses to SARS-CoV-2 in rituximab-treated patients and highlight the significance of quantitative assessment of serum antibody affinity and concentration in monitoring anti-viral immunity, viral escape, and the evolution of the humoral response.

Computational Investigation of Selected Spike Protein Mutations in SARS-CoV-2: Delta, Omicron, and Some Circulating Subvariants

Among the multiple SARS-CoV-2 variants recently reported, the Delta variant has generated the most perilous and widespread effects. Another variant, Omicron, has been identified specifically for its high transmissibility. Omicron contains numerous spike (S) protein mutations and numbers much larger than those of its predecessor variants. In this report, the author has discussed some essential structural aspects and time-based structure changes of a selected set of spike protein mutations within the Delta and Omicron variants. The expected impact of multiple point mutations within the spike protein's receptor-binding domain (RBD) and S1 of these variants are examined. Additionally, the RBDs of the more recently emerged subvariants BA.4, BA.5, and BA.2.12.1 are discussed. Within the latter group, BA.5 represents the most prevalent form of SARS-CoV-2 globally until recently. This computational work also briefly explores the temporal mutation profile for the currently circulating variants of interest (VOIs), variants under monitoring (VUMs), and variants being monitored (VBMs) including XBB.1.5, BQ.1, BA.2.75, CH.1.1, XBB, XBF, EG.5 (or Eris), and BA.2.86 (or Pirola). It is expected that these structural data can facilitate the tasks of identifying drug targets and neutralizing antibodies for the evolving variants/subvariants of SARS-CoV-2.

AlphaFold2-Enabled Atomistic Modeling of Epistatic Binding Mechanisms for the SARS-CoV-2 Spike Omicron XBB.1.5, EG.5 and FLip Variants: Convergent Evolution Hotspots Cooperate to Control Stability and Conformational Adaptability in Balancing ACE2 Binding and Antibody Resistance

In this study, we combined AI-based atomistic structural modeling and microsecond molecular simulations of the SARS-CoV-2 Spike complexes with the host receptor ACE2 for XBB.1.5+L455F, XBB.1.5+F456L(EG.5) and XBB.1.5+L455F/F456L (FLip) lineages to examine the mechanisms underlying the role of convergent evolution hotspots in balancing ACE2 binding and antibody evasion. Using the ensemble-based mutational scanning of the spike protein residues and physics-based rigorous computations of binding affinities, we identified binding energy hotspots and characterized molecular basis underlying epistatic couplings between convergent mutational hotspots. Consistent with the experiments, the results revealed the mediating role of Q493 hotspot in synchronization of epistatic couplings between L455F and F456L mutations providing a quantitative insight into the mechanism underlying differences between XBB lineages. Mutational profiling is combined with network-based model of epistatic couplings showing that the Q493, L455 and F456 sites mediate stable communities at the binding interface with ACE2 and can serve as stable mediators of non-additive couplings. Structure-based mutational analysis of Spike protein binding with the class 1 antibodies quantified the critical role of F456L and F486P mutations in eliciting strong immune evasion response. The results of this analysis support a mechanism in which the emergence of EG.5 and FLip variants may have been dictated by leveraging strong epistatic effects between several convergent revolutionary hotspots that provide synergy between the improved ACE2 binding and broad neutralization resistance. This interpretation is consistent with the notion that functionally balanced substitutions which simultaneously optimize immune evasion and high ACE2 affinity may continue to emerge through lineages with beneficial pair or triplet combinations of RBD mutations involving mediators of epistatic couplings and sites in highly adaptable RBD regions.

Novel Spike-stabilized trimers with improved production protect K18-hACE2 mice and golden Syrian hamsters from the highly pathogenic SARS-CoV-2 Beta variant

Most COVID-19 vaccines are based on the SARS-CoV-2 Spike glycoprotein (S) or their subunits. However, S shows some structural instability that limits its immunogenicity and production, hampering the development of recombinant S-based vaccines. The introduction of the K986P and V987P (S-2P) mutations increases the production and immunogenicity of the recombinant S trimer, suggesting that these two parameters are related. Nevertheless, S-2P still shows some molecular instability and it is produced with low yield. Here we described a novel set of mutations identified by molecular modeling and located in the S2 region of the S-2P that increase its production up to five-fold. Besides their immunogenicity, the efficacy of two representative S-2P-based mutants, S-29 and S-21, protecting from a heterologous SARS-CoV-2 Beta variant challenge was assayed in K18-hACE2 mice (an animal model of severe SARS-CoV-2 disease) and golden Syrian hamsters (GSH) (a moderate disease model). S-21 induced higher level of WH1 and Delta variants neutralizing antibodies than S-2P in K18-hACE2 mice three days after challenge. Viral load in nasal turbinate and oropharyngeal samples were reduced in S-21 and S-29 vaccinated mice. Despite that, only the S-29 protein protected 100% of K18-hACE2 mice from severe disease. When GSH were analyzed, all immunized animals were protected from disease development irrespectively of the immunogen they received. Therefore, the higher yield of S-29, as well as its improved immunogenicity and efficacy protecting from the highly pathogenic SARS-CoV-2 Beta variant, pinpoint the S-29 mutant as an alternative to the S-2P protein for future SARS-CoV-2 vaccine development.

Analysis of the Tropism of SARS-CoV-2 Based on the Host Interactome of the Spike Protein

The β-coronavirus SARS-CoV-2 causes severe acute respiratory syndrome (COVID-19) in humans. It enters and infects epithelial airway cells upon binding of the receptor binding domain (RBD) of the virus entry protein spike to the host receptor protein Angiotensin Converting Enzyme 2 (ACE2). Here, we used coimmunoprecipitation coupled with bottom-up mass spectrometry to identify host proteins that engaged with the spike protein in human bronchial epithelial cells (16HBEo-). We found that the spike protein bound to extracellular laminin and thrombospondin and endoplasmatic reticulum (ER)-resident DJB11 and FBX2 proteins. The ER-resident proteins UGGT1, CALX, HSP7A, and GRP78/BiP bound preferentially to the original Wuhan D614 over the mutated G614 spike protein in the more rapidly spreading Alpha SARS-CoV-2 strain. The increase in protein binding to the D614 spike might be explained by higher accessibility of cryptic sites in "RDB open" and "S2 only" D614 spike protein conformations and may enable SARS-CoV-2 to infect additional, ACE2-negative cell types. Moreover, a novel proteome-based cell type set enrichment analysis (pCtSEA) found that host factors like laminin might render additional cell types such as macrophages and epithelial cells in the nephron permissive to SARS-CoV-2 infection.

Discovering covalent inhibitors of protein-protein interactions from trillions of sulfur(VI) fluoride exchange-modified oligonucleotides

Molecules that covalently engage target proteins are widely used as activity-based probes and covalent drugs. The performance of these covalent inhibitors is, however, often compromised by the paradox of efficacy and risk, which demands a balance between reactivity and selectivity. The challenge is more evident when targeting protein-protein interactions owing to their low ligandability and undefined reactivity. Here we report sulfur(VI) fluoride exchange (SuFEx) in vitro selection, a general platform for high-throughput discovery of covalent inhibitors from trillions of SuFEx-modified oligonucleotides. With SuFEx in vitro selection, we identified covalent inhibitors that cross-link distinct residues of the SARS-CoV-2 spike protein at its protein-protein interaction interface with the human angiotensin-converting enzyme 2. A separate suite of covalent inhibitors was isolated for the human complement C5 protein. In both cases, we observed a clear disconnection between binding affinity and cross-linking reactivity, indicating that direct search for the aimed reactivity-as enabled by SuFEx in vitro selection-is vital for discovering covalent inhibitors of high selectivity and potency.

Lys417 acts as a molecular switch that regulates the conformation of SARS-CoV-2 spike protein

SARS-CoV-2 spike protein plays a key role in mediating viral entry and inducing host immune responses. It can adopt either an open or closed conformation based on the position of its receptor-binding domain (RBD). It is yet unclear what causes these conformational changes or how they influence the spike's functions. Here, we show that Lys417 in the RBD plays dual roles in the spike's structure: it stabilizes the closed conformation of the trimeric spike by mediating inter-spike-subunit interactions; it also directly interacts with ACE2 receptor. Hence, a K417V mutation has opposing effects on the spike's function: it opens up the spike for better ACE2 binding while weakening the RBD's direct binding to ACE2. The net outcomes of this mutation are to allow the spike to bind ACE2 with higher probability and mediate viral entry more efficiently, but become more exposed to neutralizing antibodies. Given that residue 417 has been a viral mutational hotspot, SARS-CoV-2 may have been evolving to strike a balance between infection potency and immune evasion, contributing to its pandemic spread.

Accurate Characterization of Conformational Ensembles and Binding Mechanisms of the SARS-CoV-2 Omicron BA.2 and BA.2.86 Spike Protein with the Host Receptor and Distinct Classes of Antibodies Using AlphaFold2-Augmented Integrative Computational Modeling

The latest wave SARS-CoV-2 Omicron variants displayed a growth advantage and the increased viral fitness through convergent evolution of functional hotspots that work synchronously to balance fitness requirements for productive receptor binding and efficient immune evasion. In this study, we combined AlphaFold2-based structural modeling approaches with all-atom MD simulations and mutational profiling of binding energetics and stability for prediction and comprehensive analysis of the structure, dynamics, and binding of the SARS-CoV-2 Omicron BA.2.86 spike variant with ACE2 host receptor and distinct classes of antibodies. We adapted several AlphaFold2 approaches to predict both structure and conformational ensembles of the Omicron BA.2.86 spike protein in the complex with the host receptor. The results showed that AlphaFold2-predicted conformational ensemble of the BA.2.86 spike protein complex can accurately capture the main dynamics signatures obtained from microscond molecular dynamics simulations. The ensemble-based dynamic mutational scanning of the receptor binding domain residues in the BA.2 and BA.2.86 spike complexes with ACE2 dissected the role of the BA.2 and BA.2.86 backgrounds in modulating binding free energy changes revealing a group of conserved hydrophobic hotspots and critical variant-specific contributions of the BA.2.86 mutational sites R403K, F486P and R493Q. To examine immune evasion properties of BA.2.86 in atomistic detail, we performed large scale structure-based mutational profiling of the S protein binding interfaces with distinct classes of antibodies that displayed significantly reduced neutralization against BA.2.86 variant. The results quantified specific function of the BA.2.86 mutations to ensure broad resistance against different classes of RBD antibodies. This study revealed the molecular basis of compensatory functional effects of the binding hotspots, showing that BA.2.86 lineage may have primarily evolved to improve immune escape while modulating binding affinity with ACE2 through cooperative effect of R403K, F486P and R493Q mutations. The study supports a hypothesis that the impact of the increased ACE2 binding affinity on viral fitness is more universal and is mediated through cross-talk between convergent mutational hotspots, while the effect of immune evasion could be more variant-dependent.

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.

Design, structure and plasma binding of ancestral β-CoV scaffold antigens

We report the application of ancestral sequence reconstruction on coronavirus spike protein, resulting in stable and highly soluble ancestral scaffold antigens (AnSAs). The AnSAs interact with plasma of patients recovered from COVID-19 but do not bind to the human angiotensin-converting enzyme 2 (ACE2) receptor. Cryo-EM analysis of the AnSAs yield high resolution structures (2.6-2.8 Å) indicating a closed pre-fusion conformation in which all three receptor-binding domains (RBDs) are facing downwards. The structures reveal an intricate hydrogen-bonding network mediated by well-resolved loops, both within and across monomers, tethering the N-terminal domain and RBD together. We show that AnSA-5 can induce and boost a broad-spectrum immune response against the wild-type RBD as well as circulating variants of concern in an immune organoid model derived from tonsils. Finally, we highlight how AnSAs are potent scaffolds by replacing the ancestral RBD with the wild-type sequence, which restores ACE2 binding and increases the interaction with convalescent plasma.

Homologous Pairs of Low and High Temperature Originating Proteins Spanning the Known Prokaryotic Universe

Stability of proteins at high temperature has been a topic of interest for many years, as this attribute is favourable for applications ranging from therapeutics to industrial chemical manufacturing. Our current understanding and methods for designing high-temperature stability into target proteins are inadequate. To drive innovation in this space, we have curated a large dataset, learn2thermDB, of protein-temperature examples, totalling 24 million instances, and paired proteins across temperatures based on homology, yielding 69 million protein pairs - orders of magnitude larger than the current largest. This important step of pairing allows for study of high-temperature stability in a sequence-dependent manner in the big data era. The data pipeline is parameterized and open, allowing it to be tuned by downstream users. We further show that the data contains signal for deep learning. This data offers a new doorway towards thermal stability design models.

Performance of SARS COV-2 IgG Anti-N as an Independent Marker of Exposure to SARS COV-2 in an Unvaccinated West African Population

Determination of previous SARS-COV-2 infection is hampered by the absence of a standardized test. The marker used to assess previous exposure is IgG antibody to the nucleocapsid (IgG anti-N), although it is known to wane quickly from peripheral blood. The accuracies of seven antibody tests (virus neutralization test, IgG anti-N, IgG anti-spike [anti-S], IgG anti-receptor binding domain [anti-RBD], IgG anti-N + anti-RBD, IgG anti-N + anti-S, and IgG anti-S + anti-RBD), either singly or in combination, were evaluated on 502 cryopreserved serum samples collected before the COVID-19 vaccination rollout in Kumasi, Ghana. The accuracy of each index test was measured using a composite reference standard based on a combination of neutralization test and IgG anti-N antibody tests. According to the composite reference, 262 participants were previously exposed; the most sensitive test was the virus neutralization test, with 95.4% sensitivity (95% CI: 93.6-97.3), followed by 79.0% for IgG anti-N + anti-S (95% CI: 76.3-83.3). The most specific tests were virus neutralization and IgG anti-N, both with 100% specificity. Viral neutralization and IgG anti-N + anti-S were the overall most accurate tests, with specificity/sensitivity of 100/95.2% and 79.0/92.1%, respectively. Our findings indicate that IgG anti-N alone is an inadequate marker of prior exposure to SARS COV-2 in this population. Virus neutralization assay appears to be the most accurate assay in discerning prior infection. A combination of IgG anti-N and IgG anti-S is also accurate and suited for assessment of SARS COV-2 exposure in low-resource settings.

Antibody response to SARS-CoV-2 vaccines in patients with relapsing multiple sclerosis treated with evobrutinib: A Bruton's tyrosine kinase inhibitor

BACKGROUND

Evobrutinib is an oral, central nervous system (CNS)-penetrant and highly selective covalent Bruton's tyrosine kinase inhibitor under clinical development for patients with relapsing multiple sclerosis (RMS).

OBJECTIVE

To investigate the effect of evobrutinib on immune responses in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccinated patients with RMS from a Phase II trial (NCT02975349).

METHODS

A post hoc analysis of patients with RMS who received evobrutinib 75 mg twice daily and SARS-CoV-2 vaccines during the open-label extension (n = 45) was conducted. Immunoglobulin (Ig)G anti-S1/S2-specific SARS-CoV-2 antibodies were measured using an indirect chemiluminescence immunoassay.

RESULTS

In the vaccinated subgroup, mean/minimum evobrutinib exposure pre-vaccination was 105.2/88.7 weeks. In total, 43 of 45 patients developed/increased S1/S2 IgG antibody levels post-vaccination; one patient's antibody response remained negative post-vaccination and the other had antibody levels above the upper limit of detection, both pre- and post-vaccination. Most patients (n = 36/45), regardless of pre-vaccination serostatus, had a 10-100-fold increase of antibody levels pre- to post-vaccination. Antibody levels post-booster were higher versus post-vaccination.

CONCLUSION

These results suggest evobrutinib, an investigational drug with therapeutic potential for patients with RMS, acts as an immunomodulator, that is, it inhibits aberrant immune cell responses in patients with RMS, while responsiveness to foreign de novo and recall antigens is maintained.

Cryo-EM and cryo-ET of the spike, virion, and antibody neutralization of SARS-CoV-2 and VOCs

Since the outbreak of the COVID-19 pandemic, cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET) have been demonstrated to be powerful and efficient tools in structural studies of distinct conformational states of the trimeric spike protein of SARS-CoV-2 and the VOCs as well as the intact virion. Cryo-EM has also contributed greatly to revealing the molecular basis of receptor recognition and antibody neutralization of the S trimer. Additionally, it has provided structural insights into the enhanced transformation and immune evasion of the VOCs, thus facilitating antiviral antibody and drug discovery. In this review, we summarize the contributions of cryo-EM and cryo-ET in revealing the structures of SARS-CoV-2 S trimer and intact virion and the mechanisms of receptor binding and antibody neutralization. We also highlight their prospective utilities in the development of vaccines and future therapeutics against emerging SARS-CoV-2 variants and other epidemic viruses.

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.

Direct blue 53, a biological dye, inhibits SARS-CoV-2 infection by blocking ACE2 and spike interaction in vitro and in vivo

COVID-19 is a global health problem caused by SARS-CoV-2, which has led to over 600 million infections and 6 million deaths. Developing novel antiviral drugs is of pivotal importance to slow down the epidemic swiftly. In this study, we identified five azo compounds as effective antiviral drugs to SARS-CoV-2, and mechanism study revealed their targets for impeding viral particles' ability to bind to host receptors. Direct Blue 53, which displayed the strongest inhibitory impact, inhibited five mutant strains at micromole. In vitro, mechanism study demonstrated Direct Blue 53 inhibited viral infection through interaction with the spike of SARS-CoV-2. And 25 mg/kg/d compound treatment showed 50% or 60% survival protection against lethal Delta or Omicron BA.2 infection in vivo. Taken together, our results demonstrate that azo compounds with dimethyl-biphenyl-diyl-bis(azo)bis structure may be promising anti-SARS-CoV-2 drug candidates, which provide practicable therapies with the aid of structural optimizations and further research.

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.

Atypical B cells and impaired SARS-CoV-2 neutralization following heterologous vaccination in the elderly

Suboptimal responses to a primary vaccination course have been reported in the elderly, but there is little information regarding the impact of age on responses to booster third doses. Here, we show that individuals 70 years or older (median age 73, range 70-75) who received a primary two-dose schedule with AZD1222 and booster third dose with mRNA vaccine achieve significantly lower neutralizing antibody responses against SARS-CoV-2 spike pseudotyped virus compared with those younger than 70 (median age 66, range 54-69) at 1 month post booster. Impaired neutralization potency and breadth post third dose in the elderly is associated with circulating "atypical" spike-specific B cells expressing CD11c and FCRL5. However, when considering individuals who received three doses of mRNA vaccine, we did not observe differences in neutralization or enrichment in atypical B cells. This work highlights the finding that AdV and mRNA COVID-19 vaccine formats differentially instruct the memory B cell response.

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.

Omicron Coronavirus: pH-Dependent Electrostatic Potential and Energy of Association of Spike Protein to ACE2 Receptor

The association of the S-protein of the SARS-CoV-2 beta coronavirus to ACE2 receptors of the human epithelial cells determines its contagiousness and pathogenicity. We computed the pH-dependent electric potential on the surface of the interacting globular proteins and pH-dependent Gibbs free energy at the association of the wild-type strain and the omicron variant. The calculated isoelectric points of the ACE2 receptor (pI 5.4) and the S-protein in trimeric form (pI 7.3, wild type), (pI 7.8, omicron variant), experimentally verified by isoelectric focusing, show that at pH 6-7, the S1-ACE2 association is conditioned by electrostatic attraction of the oppositely charged receptor and viral protein. The comparison of the local electrostatic potentials of the omicron variant and the wild-type strain shows that the point mutations alter the electrostatic potential in a relatively small area on the surface of the receptor-binding domain (RBD) of the S1 subunit. The appearance of seven charge-changing point mutations in RBD (equivalent to three additional positive charges) leads to a stronger S1-ACE2 association at pH 5.5 (typical for the respiratory tract) and a weaker one at pH 7.4 (characteristic of the blood plasma); this reveals the reason for the higher contagiousness but lower pathogenicity of the omicron variant in comparison to the wild-type strain.

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.

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.

Heme binding to the SARS-CoV-2 spike glycoprotein

The target for humoral immunity, SARS-CoV-2 spike glycoprotein, has become the focus of vaccine research and development. Previous work demonstrated that the N-terminal domain (NTD) of SARS-CoV-2 spike binds biliverdin-a product of heme catabolism-causing a strong allosteric effect on the activity of a subset of neutralizing antibodies. Herein, we show that the spike glycoprotein is also able to bind heme (KD = 0.5 ± 0.2 μM). Molecular modeling indicated that the heme group fits well within the same pocket on the SARS-CoV-2 spike NTD. Lined by aromatic and hydrophobic residues (W104, V126, I129, F192, F194, I203, and L226), the pocket provides a suitable environment to stabilize the hydrophobic heme. Mutagenesis of N121 has a substantive effect on heme binding (KD = 3000 ± 220 μM), confirming the pocket as a major heme binding location of the viral glycoprotein. Coupled oxidation experiments in the presence of ascorbate indicated that the SARS-CoV-2 glycoprotein can catalyze the slow conversion of heme to biliverdin. The heme trapping and oxidation activities of the spike may allow the virus to reduce levels of free heme during infection to facilitate evasion of the adaptive and innate immunity.

An ultrasensitive ELISA to assay femtomolar level SARS-CoV-2 antigen based on specific peptide and tyramine signal amplification

The SARS-CoV-2 spreading rapidly has aroused catastrophic public healthcare issues and economy crisis worldwide. It plays predominant role to rapidly and accurately diagnose the virus for effective prevention and treatment. As an abundant transmembrane protein, spike protein (SP) is one of the most valuable antigenic biomarkers for diagnosis of COVID-19. Herein a phage expression of WNLDLSQWLPPM peptide specific to SARS-CoV-2 SP was screened. Molecular docking revealed that the isolated peptide binds to major antigenic epitope locating at S2 subunit with hydrogen bonding. Taking the specific peptide as antigen sensing probe and tyramine signal amplification (TSA), an ultrasensitive "peptide-antigen-antibody" ELISA (p-ELISA) was explored, by which the limit of detection (LOD) was 14 fM and 2.8 fM SARS-CoV-2 SP antigen for first TSA and secondary TSA, respectively. Compared with the LOD by the p-ELISA by direct mode, the sensitivity with 2nd TSA enhanced 100 times. Further, the proposed p-ELISA method can detect SARS-CoV-2 pseudoviruses down to 10 and 3 TCID₅₀/mL spiked in healthy nasal swab sample with 1st TSA and 2nd TSA, separately. Thus, the proposed p-ELISA method with TSA is expected to be a promising ultrasensitive tool for rapidly detecting SARS-CoV-2 antigen to help control the infectious disease.

Age-associated B cells predict impaired humoral immunity after COVID-19 vaccination in patients receiving immune checkpoint blockade

Age-associated B cells (ABC) accumulate with age and in individuals with different immunological disorders, including cancer patients treated with immune checkpoint blockade and those with inborn errors of immunity. Here, we investigate whether ABCs from different conditions are similar and how they impact the longitudinal level of the COVID-19 vaccine response. Single-cell RNA sequencing indicates that ABCs with distinct aetiologies have common transcriptional profiles and can be categorised according to their expression of immune genes, such as the autoimmune regulator (AIRE). Furthermore, higher baseline ABC frequency correlates with decreased levels of antigen-specific memory B cells and reduced neutralising capacity against SARS-CoV-2. ABCs express high levels of the inhibitory FcγRIIB receptor and are distinctive in their ability to bind immune complexes, which could contribute to diminish vaccine responses either directly, or indirectly via enhanced clearance of immune complexed-antigen. Expansion of ABCs may, therefore, serve as a biomarker identifying individuals at risk of suboptimal responses to vaccination.

Multivalent bicyclic peptides are an effective antiviral modality that can potently inhibit SARS-CoV-2

COVID-19 has stimulated the rapid development of new antibody and small molecule therapeutics to inhibit SARS-CoV-2 infection. Here we describe a third antiviral modality that combines the drug-like advantages of both. Bicycles are entropically constrained peptides stabilized by a central chemical scaffold into a bi-cyclic structure. Rapid screening of diverse bacteriophage libraries against SARS-CoV-2 Spike yielded unique Bicycle binders across the entire protein. Exploiting Bicycles' inherent chemical combinability, we converted early micromolar hits into nanomolar viral inhibitors through simple multimerization. We also show how combining Bicycles against different epitopes into a single biparatopic agent allows Spike from diverse variants of concern (VoC) to be targeted (Alpha, Beta, Delta and Omicron). Finally, we demonstrate in both male hACE2-transgenic mice and Syrian golden hamsters that both multimerized and biparatopic Bicycles reduce viraemia and prevent host inflammation. These results introduce Bicycles as a potential antiviral modality to tackle new and rapidly evolving viruses.

Systematic conformation-to-phenotype mapping via limited deep sequencing of proteins

Non-native conformations drive protein-misfolding diseases, complicate bioengineering efforts, and fuel molecular evolution. No current experimental technique is well suited for elucidating them and their phenotypic effects. Especially intractable are the transient conformations populated by intrinsically disordered proteins. We describe an approach to systematically discover, stabilize, and purify native and non-native conformations, generated in vitro or in vivo, and directly link conformations to molecular, organismal, or evolutionary phenotypes. This approach involves high-throughput disulfide scanning (HTDS) of the entire protein. To reveal which disulfides trap which chromatographically resolvable conformers, we devised a deep-sequencing method for double-Cys variant libraries of proteins that precisely and simultaneously locates both Cys residues within each polypeptide. HTDS of the abundant E. coli periplasmic chaperone HdeA revealed distinct classes of disordered hydrophobic conformers with variable cytotoxicity depending on where the backbone was cross-linked. HTDS can bridge conformational and phenotypic landscapes for many proteins that function in disulfide-permissive environments.

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.

Accelerated waning of the humoral response to COVID-19 vaccines in obesity

Obesity is associated with an increased risk of severe Coronavirus Disease 2019 (COVID-19) infection and mortality. COVID-19 vaccines reduce the risk of serious COVID-19 outcomes; however, their effectiveness in people with obesity is incompletely understood. We studied the relationship among body mass index (BMI), hospitalization and mortality due to COVID-19 among 3.6 million people in Scotland using the Early Pandemic Evaluation and Enhanced Surveillance of COVID-19 (EAVE II) surveillance platform. We found that vaccinated individuals with severe obesity (BMI > 40 kg/m2) were 76% more likely to experience hospitalization or death from COVID-19 (adjusted rate ratio of 1.76 (95% confidence interval (CI), 1.60-1.94). We also conducted a prospective longitudinal study of a cohort of 28 individuals with severe obesity compared to 41 control individuals with normal BMI (BMI 18.5-24.9 kg/m2). We found that 55% of individuals with severe obesity had unquantifiable titers of neutralizing antibody against authentic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus compared to 12% of individuals with normal BMI (P = 0.0003) 6 months after their second vaccine dose. Furthermore, we observed that, for individuals with severe obesity, at any given anti-spike and anti-receptor-binding domain (RBD) antibody level, neutralizing capacity was lower than that of individuals with a normal BMI. Neutralizing capacity was restored by a third dose of vaccine but again declined more rapidly in people with severe obesity. We demonstrate that waning of COVID-19 vaccine-induced humoral immunity is accelerated in individuals with severe obesity. As obesity is associated with increased hospitalization and mortality from breakthrough infections, our findings have implications for vaccine prioritization policies.

Factor Xa cleaves SARS-CoV-2 spike protein to block viral entry and infection

Serine proteases (SP), including furin, trypsin, and TMPRSS2 cleave the SARS-CoV-2 spike (S) protein, enabling the virus to enter cells. Here, we show that factor (F) Xa, an SP involved in blood coagulation, is upregulated in COVID-19 patients. In contrast to other SPs, FXa exerts antiviral activity. Mechanistically, FXa cleaves S protein, preventing its binding to ACE2, and thus blocking viral entry and infection. However, FXa is less effective against variants carrying the D614G mutation common in all pandemic variants. The anticoagulant rivaroxaban, a direct FXa inhibitor, inhibits FXa-mediated S protein cleavage and facilitates viral entry, whereas the indirect FXa inhibitor fondaparinux does not. In the lethal SARS-CoV-2 K18-hACE2 model, FXa prolongs survival yet its combination with rivaroxaban but not fondaparinux abrogates that protection. These results identify both a previously unknown function for FXa and an associated antiviral host defense mechanism against SARS-CoV-2 and suggest caution in considering direct FXa inhibitors for preventing or treating thrombotic complications in COVID-19 patients.

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.

Intranasal trimeric sherpabody inhibits SARS-CoV-2 including recent immunoevasive Omicron subvariants

The emergence of increasingly immunoevasive SARS-CoV-2 variants emphasizes the need for prophylactic strategies to complement vaccination in fighting the COVID-19 pandemic. Intranasal administration of neutralizing antibodies has shown encouraging protective potential but there remains a need for SARS-CoV-2 blocking agents that are less vulnerable to mutational viral variation and more economical to produce in large scale. Here we describe TriSb92, a highly manufacturable and stable trimeric antibody-mimetic sherpabody targeted against a conserved region of the viral spike glycoprotein. TriSb92 potently neutralizes SARS-CoV-2, including the latest Omicron variants like BF.7, XBB, and BQ.1.1. In female Balb/c mice intranasal administration of just 5 or 50 micrograms of TriSb92 as early as 8 h before but also 4 h after SARS-CoV-2 challenge can protect from infection. Cryo-EM and biochemical studies reveal triggering of a conformational shift in the spike trimer as the inhibitory mechanism of TriSb92. The potency and robust biochemical properties of TriSb92 together with its resistance against viral sequence evolution suggest that TriSb92 could be useful as a nasal spray for protecting susceptible individuals from SARS-CoV-2 infection.

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.

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.