Structures of the Omicron Spike trimer with ACE2 and an anti-Omicron antibody

Mechanistic study of the transmission pattern of the SARS-CoV-2 omicron variant

The omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) characterized by 30 mutations in its spike protein, has rapidly spread worldwide since November 2021, significantly exacerbating the ongoing COVID-19 pandemic. In order to investigate the relationship between these mutations and the variant's high transmissibility, we conducted a systematic analysis of the mutational effect on spike-angiotensin-converting enzyme-2 (ACE2) interactions and explored the structural/energy correlation of key mutations, utilizing a reliable coarse-grained model. Our study extended beyond the receptor-binding domain (RBD) of spike trimer through comprehensive modeling of the full-length spike trimer rather than just the RBD. Our free-energy calculation revealed that the enhanced binding affinity between the spike protein and the ACE2 receptor is correlated with the increased structural stability of the isolated spike protein, thus explaining the omicron variant's heightened transmissibility. The conclusion was supported by our experimental analyses involving the expression and purification of the full-length spike trimer. Furthermore, the energy decomposition analysis established those electrostatic interactions make major contributions to this effect. We categorized the mutations into four groups and established an analytical framework that can be employed in studying future mutations. Additionally, our calculations rationalized the reduced affinity of the omicron variant towards most available therapeutic neutralizing antibodies, when compared with the wild type. By providing concrete experimental data and offering a solid explanation, this study contributes to a better understanding of the relationship between theories and observations and lays the foundation for future investigations.

Molecular Mechanism of Interaction between DNA Aptamer and Receptor-Binding Domain of Severe Acute Respiratory Syndrome Coronavirus 2 Variants Revealed by Steered Molecular Dynamics Simulations

The ongoing SARS-CoV-2 pandemic has underscored the urgent need for versatile and rapidly deployable antiviral strategies. While vaccines have been pivotal in controlling the spread of the virus, the emergence of new variants continues to pose significant challenges to global health. Here, our study focuses on a novel approach to antiviral therapy using DNA aptamers, short oligonucleotides with high specificity and affinity for their targets, as potential inhibitors against the spike protein of SARS-CoV-2 variants Omicron and JN.1. Our research utilizes steered molecular dynamics (SMD) simulations to elucidate the binding mechanisms of a specifically designed DNA aptamer, AM032-4, to the receptor-binding domain (RBD) of the aforementioned variants. The simulations reveal detailed molecular insights into the aptamer-RBD interaction, demonstrating the aptamer's potential to maintain effective binding in the face of rapid viral evolution. Our work not only demonstrates the dynamic interaction between aptamer-RBD for possible antiviral therapy but also introduces a computational method to study aptamer-protein interactions.

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.

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.

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.

A panel of multivalent nanobodies broadly neutralizing Omicron subvariants and recombinant

The emerging Omicron subvariants have a remarkable ability to spread and escape nearly all current monoclonal antibody (mAb) treatments. Although the virulence of SARS-CoV-2 has now diminished, it remains a significant threat to public health due to its high transmissibility and susceptibility to mutation. Therefore, it is urgent to develop broad-acting and potent therapeutics targeting current and emerging Omicron variants. Here, we identified a panel of Omicron BA.1 spike receptor-binding domain (RBD)-targeted nanobodies (Nbs) from a naive alpaca VHH library. This panel of Nbs exhibited high binding affinity to the spike RBD of wild-type, Alpha B.1.1.7, Beta B.1.351, Delta plus, Omicron BA.1, and BA.2. Through multivalent Nb construction, we obtained a subpanel of ultrapotent neutralizing Nbs against Omicron BA.1, BA.2, BF.7 and even emerging XBB.1.5, and XBB.1.16 pseudoviruses. Protein structure prediction and docking analysis showed that Nb trimer 2F2E5 targets two independent RBD epitopes, thus minimizing viral escape. Taken together, we obtained a panel of broad and ultrapotent neutralizing Nbs against Omicron BA.1, Omicron BA.2, BF.7, XBB.1.5, and XBB.1.16. These multivalent Nbs hold great promise for the treatment against SARS-CoV-2 infection and could possess a superwide neutralizing breadth against novel omicron mutants or recombinants.

The Diverse Nature of the Molecular Interactions That Govern the COV-2 Variants' Cell Receptor Affinity Ranking and Its Experimental Variability

A critical determinant of infectivity and virulence of the most infectious and or lethal variants of concern (VOCs): Wild Type, Delta and Omicron is related to the binding interactions between the receptor-binding domain of the spike and its host receptor, the initial step in cell infection. It is of the utmost importance to understand how mutations of a viral strain, especially those that are in the viral spike, affect the resulting infectivity of the emerging VOC, knowledge that could help us understand the variant virulence and inform the therapies applied or the vaccines developed. For this sake, we have applied a battery of computational protocols of increasing complexity to the calculation of the spike binding affinity for three variants of concern to the ACE2 cell receptor. The results clearly illustrate that the attachment of the spikes of the Delta and Omicron variants to the receptor originates through different molecular interaction mechanisms. All our protocols unanimously predict that the Delta variant has the highest receptor-binding affinity, while the Omicron variant displays a substantial variability in the binding affinity of the spike that relates to the structural plasticity of the Omicron spike-receptor complex. We suggest that the latter result could explain (at least in part) the variability of the in vitro binding results for this VOC and has led us to suggest a reason for the lower virulence of the Omicron variant as compared to earlier strains. Several hypotheses have been developed around this subject.

Transcriptional regulation of SARS-CoV-2 receptor ACE2 by SP1

Angiotensin-converting enzyme 2 (ACE2) is a major cell entry receptor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The induction of ACE2 expression may serve as a strategy by SARS-CoV-2 to facilitate its propagation. However, the regulatory mechanisms of ACE2 expression after viral infection remain largely unknown. Using 45 different luciferase reporters, the transcription factors SP1 and HNF4α were found to positively and negatively regulate ACE2 expression, respectively, at the transcriptional level in human lung epithelial cells (HPAEpiCs). SARS-CoV-2 infection increased the transcriptional activity of SP1 while inhibiting that of HNF4α. The PI3K/AKT signaling pathway, activated by SARS-CoV-2 infection, served as a crucial regulatory node, inducing ACE2 expression by enhancing SP1 phosphorylation-a marker of its activity-and reducing the nuclear localization of HNF4α. However, colchicine treatment inhibited the PI3K/AKT signaling pathway, thereby suppressing ACE2 expression. In Syrian hamsters (Mesocricetus auratus) infected with SARS-CoV-2, inhibition of SP1 by either mithramycin A or colchicine resulted in reduced viral replication and tissue injury. In summary, our study uncovers a novel function of SP1 in the regulation of ACE2 expression and identifies SP1 as a potential target to reduce SARS-CoV-2 infection.

Leveraging Artificial Intelligence to Expedite Antibody Design and Enhance Antibody-Antigen Interactions

This perspective sheds light on the transformative impact of recent computational advancements in the field of protein therapeutics, with a particular focus on the design and development of antibodies. Cutting-edge computational methods have revolutionized our understanding of protein-protein interactions (PPIs), enhancing the efficacy of protein therapeutics in preclinical and clinical settings. Central to these advancements is the application of machine learning and deep learning, which offers unprecedented insights into the intricate mechanisms of PPIs and facilitates precise control over protein functions. Despite these advancements, the complex structural nuances of antibodies pose ongoing challenges in their design and optimization. Our review provides a comprehensive exploration of the latest deep learning approaches, including language models and diffusion techniques, and their role in surmounting these challenges. We also present a critical analysis of these methods, offering insights to drive further progress in this rapidly evolving field. The paper includes practical recommendations for the application of these computational techniques, supplemented with independent benchmark studies. These studies focus on key performance metrics such as accuracy and the ease of program execution, providing a valuable resource for researchers engaged in antibody design and development. Through this detailed perspective, we aim to contribute to the advancement of antibody design, equipping researchers with the tools and knowledge to navigate the complexities of this field.

Computational design and engineering of self-assembling multivalent microproteins with therapeutic potential against SARS-CoV-2

Multivalent drugs targeting homo-oligomeric viral surface proteins, such as the SARS-CoV-2 trimeric spike (S) protein, have the potential to elicit more potent and broad-spectrum therapeutic responses than monovalent drugs by synergistically engaging multiple binding sites on viral targets. However, rational design and engineering of nanoscale multivalent protein drugs are still lacking. Here, we developed a computational approach to engineer self-assembling trivalent microproteins that simultaneously bind to the three receptor binding domains (RBDs) of the S protein. This approach involves four steps: structure-guided linker design, molecular simulation evaluation of self-assembly, experimental validation of self-assembly state, and functional testing. Using this approach, we first designed trivalent constructs of the microprotein miniACE2 (MP) with different trimerization scaffolds and linkers, and found that one of the constructs (MP-5ff) showed high trimerization efficiency, good conformational homogeneity, and strong antiviral neutralizing activity. With its trimerization unit (5ff), we then engineered a trivalent nanobody (Tr67) that exhibited potent and broad neutralizing activity against the dominant Omicron variants, including XBB.1 and XBB.1.5. Cryo-EM complex structure confirmed that Tr67 stably binds to all three RBDs of the Omicron S protein in a synergistic form, locking them in the "3-RBD-up" conformation that could block human receptor (ACE2) binding and potentially facilitate immune clearance. Therefore, our approach provides an effective strategy for engineering potent protein drugs against SARS-CoV-2 and other deadly coronaviruses.

Deciphering the free energy landscapes of SARS-CoV-2 wild type and Omicron variant interacting with human ACE2

The binding of the receptor binding domain (RBD) of the SARS-CoV-2 spike protein to the host cell receptor angiotensin-converting enzyme 2 (ACE2) is the first step in human viral infection. Therefore, understanding the mechanism of interaction between RBD and ACE2 at the molecular level is critical for the prevention of COVID-19, as more variants of concern, such as Omicron, appear. Recently, atomic force microscopy has been applied to characterize the free energy landscape of the RBD-ACE2 complex, including estimation of the distance between the transition state and the bound state, xu. Here, using a coarse-grained model and replica-exchange umbrella sampling, we studied the free energy landscape of both the wild type and Omicron subvariants BA.1 and XBB.1.5 interacting with ACE2. In agreement with experiment, we find that the wild type and Omicron subvariants have similar xu values, but Omicron binds ACE2 more strongly than the wild type, having a lower dissociation constant KD.

Rapid Profiling of the Glycosylation Effects on the Binding of SARS-CoV-2 Spike Protein to Angiotensin-Converting Enzyme 2 Using MALDI-MS with High Mass Detection

The spike protein receptor-binding domain (RBD) of SARS-CoV-2 binds directly to angiotensin-converting enzyme 2 (ACE2), mediating the host cell entry of SARS-CoV-2. Both spike protein and ACE2 are highly glycosylated, which can regulate the binding. Here, we utilized high-mass MALDI-MS with chemical cross-linking for profiling the glycosylation effects on the binding between RBD and ACE2. Overall, it was found that ACE2 glycosylation affects the binding more strongly than does RBD glycosylation. The binding affinity was improved after desialylation or partial deglycosylation (N690) of ACE2, while it decreased after degalactosylation. ACE2 can form dimers in solution, which bind more tightly to the RBD than the ACE2 monomers. The ACE2 dimerization and the binding of RBD to dimeric ACE2 can also be improved by the desialylation or deglycosylation of ACE2. Partial deglycosylation of ACE2 increased the dimerization of ACE2 and the binding affinity of RBD and ACE2 by more than a factor of 2, suggesting its high potential for neutralizing SARS-CoV-2. The method described in the work provided a simple way to analyze the protein-protein interaction without sample purification. It can be widely used for rapid profiling of glycosylation effects on protein-protein interaction for glycosylation-related diseases and the study of multiple interactions between protein and protein aggregates in a single system.

A flexible, image-based, high-throughput platform encompassing in-depth cell profiling to identify broad-spectrum coronavirus antivirals with limited off-target effects

The recent pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) posed a major threat to global health. Although the World Health Organization ended the public health emergency status, antiviral drugs are needed to address new variants of SARS-CoV-2 and future pandemics. To identify novel broad-spectrum coronavirus drugs, we developed a high-content imaging platform compatible with high-throughput screening. The platform is broadly applicable as it can be adapted to include various cell types, viruses, antibodies, and dyes. We demonstrated that the antiviral activity of compounds against SARS-CoV-2 variants (Omicron BA.5 and Omicron XBB.1.5), SARS-CoV, and human coronavirus 229E could easily be assessed. The inclusion of cellular dyes and immunostaining in combination with in-depth image analysis enabled us to identify compounds that induced undesirable phenotypes in host cells, such as changes in cell morphology or in lysosomal activity. With the platform, we screened ∼900K compounds and triaged hits, thereby identifying potential candidate compounds carrying broad-spectrum activity with limited off-target effects. The flexibility and early-stage identification of compounds with limited host cell effects provided by this high-content imaging platform can facilitate coronavirus drug discovery. We anticipate that its rapid deployability and fast turnaround can also be applied to combat future pandemics.

Marine natural products and human immunity: novel biomedical resources for anti-infection of SARS-CoV-2 and related cardiovascular disease

Marine natural products (MNPs) and marine organisms include sea urchin, sea squirts or ascidians, sea cucumbers, sea snake, sponge, soft coral, marine algae, and microalgae. As vital biomedical resources for the discovery of marine drugs, bioactive molecules, and agents, these MNPs have bioactive potentials of antioxidant, anti-infection, anti-inflammatory, anticoagulant, anti-diabetic effects, cancer treatment, and improvement of human immunity. This article reviews the role of MNPs on anti-infection of coronavirus, SARS-CoV-2 and its major variants (such as Delta and Omicron) as well as tuberculosis, H. Pylori, and HIV infection, and as promising biomedical resources for infection related cardiovascular disease (irCVD), diabetes, and cancer. The anti-inflammatory mechanisms of current MNPs against SARS-CoV-2 infection are also discussed. Since the use of other chemical agents for COVID-19 treatment are associated with some adverse effects in cardiovascular system, MNPs have more therapeutic advantages. Herein, it's time to protect this ecosystem for better sustainable development in the new era of ocean economy. As huge, novel and promising biomedical resources for anti-infection of SARS-CoV-2 and irCVD, the novel potential mechanisms of MNPs may be through multiple targets and pathways regulating human immunity and inhibiting inflammation. In conclusion, MNPs are worthy of translational research for further clinical application.

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

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

Identification of a broad sarbecovirus neutralizing antibody targeting a conserved epitope on the receptor-binding domain

Omicron, as the emerging variant with enhanced vaccine tolerance, has sharply disrupted most therapeutic antibodies. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) belongs to the subgenus Sarbecovirus, members of which share high sequence similarity. Herein, we report one sarbecovirus antibody, 5817, which has broad-spectrum neutralization capacity against SARS-CoV-2 variants of concern (VOCs) and SARS-CoV, as well as related bat and pangolin viruses. 5817 can hardly compete with six classes of receptor-binding-domain-targeted antibodies grouped by structural classifications. No obvious impairment in the potency is detected against SARS-CoV-2 Omicron and subvariants. The cryoelectron microscopy (cryo-EM) structure of neutralizing antibody 5817 in complex with Omicron spike reveals a highly conserved epitope, only existing at the receptor-binding domain (RBD) open state. Prophylactic and therapeutic administration of 5817 potently protects mice from SARS-CoV-2 Beta, Delta, Omicron, and SARS-CoV infection. This study reveals a highly conserved cryptic epitope targeted by a broad sarbecovirus neutralizing antibody, which would be beneficial to meet the potential threat of pre-emergent SARS-CoV-2 VOCs.

Generation of SARS-CoV-2 spike receptor binding domain mutants and functional screening for immune evaders using a novel lentivirus-based system

The emergence of rapid and continuous mutations of severe acute respiratory syndrome 2 (SARS-CoV-2) spike glycoprotein that increased with the Omicron variant points out the necessity to anticipate such mutations for conceiving specific and adaptable therapies to avoid another pandemic. The crucial target for the antibody treatment and vaccine design is the receptor binding domain (RBD) of the SARS-CoV-2 spike. It is also the site where the virus has shown its high ability to mutate and consequently escape immune response. We developed a robust and simple method for generating a large number of functional SARS-CoV-2 spike RBD mutants by error-prone PCR and a novel nonreplicative lentivirus-based system. We prepared anti-RBD wild type (WT) polyclonal antibodies and used them to screen and select for mutant libraries that escape inhibition of virion entry into recipient cells expressing human angiotensin-converting enzyme 2 and transmembrane serine protease 2. We isolated, cloned, and sequenced six mutants totally bearing nine mutation sites. Eight mutations were found in successive WT variants, including Omicron and other recombinants, whereas one is novel. These results, together with the detailed functional analyses of two mutants provided the proof of concept for our approach.

Molecular interactions between perlecan LG3 and the SARS-CoV-2 spike protein receptor binding domain

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused a global health crisis with significant clinical morbidity and mortality. While angiotensin-converting enzyme 2 (ACE2) is the primary receptor for viral entry, other cell surface and extracellular matrix proteins may also bind to the viral receptor binding domain (RBD) within the SARS-CoV-2 spike protein. Recent studies have implicated heparan sulfate proteoglycans, specifically perlecan LG3, in facilitating SARS-CoV-2 binding to ACE2. However, the role of perlecan LG3 in SARS-CoV-2 pathophysiology is not well understood. In this study, we investigated the binding interactions between the SARS-CoV-2 spike protein RBD and perlecan LG3 through molecular modeling simulations and surface plasmon resonance (SPR) experiments. Our results indicate stable binding between LG3 and SARS-CoV-2 spike protein RBD, which may potentially enhance RBD-ACE2 interactions. These findings shed light on the role of perlecan LG3 in SARS-CoV-2 infection and provide insight into SARS-CoV-2 pathophysiology and potential therapeutic strategy for COVID-19.

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.

Absolute binding free energies of the antiviral peptide ATN-161 with protein targets of SARS-CoV-2

The interactions of the antiviral pentapeptide ATN-161 with the closed and open conformations of the α5β1 integrin, the SARS-CoV-2 major protease, and the omicron variant spike protein complexed with hACE2 were studied using molecular docking and molecular dynamics simulation. Molecular docking was performed to obtain ATN-161 binding poses with these studied protein targets. Subsequently, molecular dynamics simulations were performed to verify the ligand stability at the binding site of each protein target. Pulling simulations, umbrella sampling, and weighted histogram analysis method were used to obtain the potential of mean force of each system and calculate the Gibbs free energy of binding for the ATN-161 peptide in each binding site of these protein targets. The results showed that ATN-161 binds to α5β1 integrin in its active and inactive form, binds weakly to the omicron variant spike protein complexed with hACE2, and strongly binds to the main protease target.Communicated by Ramaswamy H. Sarma.

Long-term immune response to Omicron-specific mRNA vaccination in mice, hamsters, and nonhuman primates

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron and its subvariants (such as BQ.1, XBB and the latest variants, including XBB.1.16, EG.5, and BA.2.86), as the dominant variants, currently account for almost all new infections in the world due to their high transmissibility and immune escape ability. Omicron-specific mRNA vaccines showed great potential to protect against Omicron infections. However, whether the vaccine could provide long-term protection is unknown. Toward this goal, we evaluated the immunogenicity of a preclinical Omicron (BA.1)-specific mRNA vaccine (SOmicron-6P) in different animal models. SOmicron-6P induced the highest levels of antibody titers at 1-2 weeks in different animals after the second dose. Even 9 months after the immunization, we observed modest neutralizing activity against Omicron subvariants in macaques. In addition, immunological memory cells can be rapidly reactivated upon stimulation. SOmicron-6P at concentrations higher than 10 μg effectively protected hamsters from BA.1 challenge 253 days after the first immunization, which could be attributed to the reactivation of immune systems. In addition, the toxicity tests conducted in rats revealed a highly favorable biosafety profile for SOmicron-6P, even at high dosages. Our data suggest that the Omicron-specific mRNA vaccine is highly effective and safe in animal models and provides long-term immunologic protection against SARS-CoV-2 Omicron infections.

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

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

Exploring the binding capacity of lactic acid bacteria derived bacteriocins against RBD of SARS-CoV-2 Omicron variant by molecular simulations

The changes in the SARS-CoV-2 genome have resulted in the emergence of new variants. Some of the variants have been classified as variants of concern (VOC). These strains have higher transmission rate and improved fitness. One of the prevalent were the Omicron variant. Unlike previous VOCs, the Omicron possesses fifteen mutations on the spike protein's receptor binding domain (RBD). The modifications of spike protein's key amino acid residues facilitate the virus' binding capability against ACE2, resulting in an increase in the infectiousness of Omicron variant. Consequently, investigating the prevention and treatment of the Omicron variant is crucial. In the present study, we aim to explore the binding capacity of twenty-two bacteriocins derived from Lactic Acid Bacteria (LAB) against the Omicron variant by using protein-peptidedocking and molecular dynamics (MD) simulations. The Omicron variant RBD was prepared by introducing fifteen mutations using PyMol. The protein-peptide complexes were obtained using HADDOCK v2.4 docking webserver. Top scoring complexes obtained from HADDOCK webserver were retrieved and submitted to the PRODIGY server for the prediction of binding energies. RBD-bacteriocin complexes were subjected to MD simulations. We discovered promising peptide-based therapeutic candidates for the inhibition of Omicron variant for example Salivaricin B, Pediocin PA 1, Plantaricin W, Lactococcin mmfii and Enterocin A. The lead bacteriocins, except Enterocin A, are biosynthesized by food-grade lactic acid bacteria. Our study puts forth a preliminary information regarding potential utilization of food-grade LAB-derived bacteriocins, particularly Salivaricin B and Pediocin PA 1, for Covid-19 treatment and prophylaxis.Communicated by Ramaswamy H. Sarma.

Protection of the receptor binding domain (RBD) dimer against SARS-CoV-2 and its variants

IMPORTANCE

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants achieved immune escape and became less virulent and easily transmissible through rapid mutation in the spike protein, thus the efficacy of vaccines on the market or in development continues to be challenged. Updating the vaccine, exploring compromise vaccination strategies, and evaluating the efficacy of candidate vaccines for the emerging variants in a timely manner are important to combat complex and volatile SARS-CoV-2. This study reports that vaccines prepared from the dimeric receptor-binding domain (RBD) recombinant protein, which can be quickly produced using a mature and stable process platform, had both good immunogenicity and protection in vivo and could completely protect rodents from lethal challenge by SARS-CoV-2 and its variants, including the emerging Omicron XBB.1.16, highlighting the value of dimeric recombinant vaccines in the post-COVID-19 era.

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.

Antibodies to S2 domain of SARS-CoV-2 spike protein in Moderna mRNA vaccinated subjects sustain antibody-dependent NK cell-mediated cell cytotoxicity against Omicron BA.1

Vaccination with the primary two-dose series of SARS-CoV-2 mRNA protects against infection with the ancestral strain, and limits the presentation of severe disease after re-infection by multiple variants of concern (VOC), including Omicron, despite the lack of a strong neutralizing response to these variants. We compared antibody responses in serum samples collected from mRNA-1273 (Moderna) vaccinated subjects to identify mechanisms of immune escape and cross-protection. Using pseudovirus constructs containing domain-specific amino acid changes representative of Omicron BA.1, combined with domain competition and RBD-antibody depletion, we showed that RBD antibodies were primarily responsible for virus neutralization and variant escape. Antibodies to NTD played a less significant role in antibody neutralization but acted along with RBD to enhance neutralization. S2 of Omicron BA.1 had no impact on neutralization escape, suggesting it is a less critical domain for antibody neutralization; however, it was as capable as S1 at eliciting IgG3 responses and NK-cell mediated, antibody-dependent cell cytotoxicity (ADCC). Antibody neutralization and ADCC activities to RBD, NTD, and S1 were all prone to BA.1 escape. In contrast, ADCC activities to S2 resisted BA.1 escape. In conclusion, S2 antibodies showed potent ADCC function and resisted Omicron BA.1 escape, suggesting that S2 contributes to cross-protection against Omicron BA.1. In line with its conserved nature, S2 may hold promise as a vaccine target against future variants of SARS-CoV-2.

Evolving antibody evasion and receptor affinity of the Omicron BA.2.75 sublineage of SARS-CoV-2

SARS-CoV-2 Omicron BA.2.75 has diversified into multiple subvariants with additional spike mutations and several are expanding in prevalence, particularly CH.1.1 and BN.1. Here, we investigated the viral receptor affinities and neutralization evasion properties of major BA.2.75 subvariants actively circulating in different regions worldwide. We found two distinct evolutionary pathways and three newly identified mutations that shaped the virological features of these subvariants. One phenotypic group exhibited a discernible decrease in viral receptor affinities, but a noteworthy increase in resistance to antibody neutralization, as exemplified by CH.1.1, which is apparently as resistant as XBB.1.5. In contrast, a second group demonstrated a substantial increase in viral receptor affinity but only a moderate increase in antibody evasion, as exemplified by BN.1. We also observed that all prevalent SARS-CoV-2 variants in the circulation presently, except for BN.1, exhibit profound levels of antibody evasion, suggesting this is the dominant determinant of virus transmissibility today.

Antibody evasion associated with the RBD significant mutations in several emerging SARS-CoV-2 variants and its subvariants

Since the origin of the wild strain of SARS-CoV-2, several variants have emerged, which were designated as VOC, VOI, and VUM from time to time. The Omicron variant is noted as the recent VOC. After the origin of the Omicron variant on November 2021, several subvariants of Omicron have originated subsequently, like BA.1/2, BA.2.75/2.75.2, BA.4/5, BF.7, BQ.1/1.1, XBB.1/1.5, etc. which are circulated throughout the globe. Scientists reported that antibody escape is a common phenomenon observed in all the previous VOCs, VOIs, including Omicron and its subvariants. The mutations in the NTD (N-terminal domain) and RBD (Receptor-binding domain) of the spike of these variants and subvariants are responsible for antibody escape. At the same time, it has been noted that spike RBD mutations have been increasing in the last few months. This review illustrates significant RBD mutations namely R346T, K417N/T, L452R, N460K E484A/K/Q, and N501Y found in the previous emerging SARS-CoV-2 variants, including Omicron and its subvariants in high frequency and their role in antibody evasion and immune evasion. The review also describes the different classes of nAb responsible for antibody escape in SARS-CoV-2 variants and the molecular perspective of the mutation in nAb escape. It will help the future researchers to develop efficient vaccines which can finally prevent the pandemic.

Defining neutralization and allostery by antibodies against COVID-19 variants

The changing landscape of SARS-CoV-2 Spike protein is linked to the emergence of variants, immune-escape and reduced efficacy of the existing repertoire of anti-viral antibodies. The functional activity of neutralizing antibodies is linked to their quaternary changes occurring as a result of antibody-Spike trimer interactions. Here, we reveal the conformational dynamics and allosteric perturbations linked to binding of novel human antibodies and the viral Spike protein. We identified epitope hotspots, and associated changes in Spike dynamics that distinguish weak, moderate and strong neutralizing antibodies. We show the impact of mutations in Wuhan-Hu-1, Delta, and Omicron variants on differences in the antibody-induced conformational changes in Spike and illustrate how these render certain antibodies ineffective. Antibodies with similar binding affinities may induce destabilizing or stabilizing allosteric effects on Spike, with implications for neutralization efficacy. Our results provide mechanistic insights into the functional modes and synergistic behavior of human antibodies against COVID-19 and may assist in designing effective antiviral strategies.

Structure adaptation in Omicron SARS-CoV-2/hACE2: Biophysical origins of evolutionary driving forces

Since its emergence, the COVID-19 threat has been sustained by a series of transmission waves initiated by new variants of the SARS-CoV-2 virus. Some of these arise with higher transmissivity and/or increased disease severity. Here, we use molecular dynamics simulations to examine the modulation of the fundamental interactions between the receptor binding domain (RBD) of the spike glycoprotein and the host cell receptor (human angiotensin-converting enzyme 2 [hACE2]) arising from Omicron variant mutations (BA.1 and BA.2) relative to the original wild-type strain. Our key findings are that glycans play a vital role at the RBD···hACE2 interface for the Omicrons, and the interplay between glycans and sequence mutations leads to enhanced binding. We find significant structural differences in the complexes, which overall bring the spike protein and its receptor into closer proximity. These are consistent with and attributed to the higher positive charge on the RBD conferred by BA.1 and BA.2 mutations relative to the wild-type. However, further differences between subvariants BA.1 and BA.2 (which have equivalent RBD charges) are also evident: mutations reduce interdomain interactions between the up chain and its clockwise neighbor chain in particular for the latter, resulting in enhanced flexibility for BA.2. Consequently, we see occurrence of additional close contacts in one replica of BA.2, which include binding to hACE2 by a second RBD in addition to the up chain. Although this motif is not seen in BA.1, we find that the Omicrons can directly/indirectly bind a down-RBD to hACE2 through glycans: the role of the glycan on N90 of hACE2 switches from inhibiting to facilitating the binding to Omicron spike protein via glycan-protein lateral interactions. These structural and electrostatic differences offer further insight into the mechanisms by which viral mutations modulate host cell binding and provide a biophysical basis for evolutionary driving forces.

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.

Designing peptides predicted to bind to the omicron variant better than ACE2 via computational protein design and molecular dynamics

Brought about by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), coronavirus disease (COVID-19) pandemic has resulted in large numbers of worldwide deaths and cases. Several SARS-CoV-2 variants have evolved, and Omicron (B.1.1.529) was one of the important variants of concern. It gets inside human cells by using its S1 subunit's receptor-binding domain (SARS-CoV-2-RBD) to bind to Angiotensin-converting enzyme 2 receptor's peptidase domain (ACE2-PD). Using peptides to inhibit binding interactions (BIs) between ACE2-PD and SARS-CoV-2-RBD is one of promising COVID-19 therapies. Employing computational protein design (CPD) as well as molecular dynamics (MD), this study used ACE2-PD's α1 helix to generate novel 25-mer peptide binders (SPB25) of Omicron RBD that have predicted binding affinities (ΔGbind (MM‑GBSA)) better than ACE2 by increasing favorable BIs between SPB25 and the conserved residues of RBD. Results from MD and the MM-GBSA method identified two best designed peptides (SPB25T7L/K11A and SPB25T7L/K11L with ΔGbind (MM‑GBSA) of -92.4 ± 0.4 and -95.7 ± 0.5 kcal/mol, respectively) that have better ΔGbind (MM‑GBSA) to Omicron RBD than ACE2 (-87.9 ± 0.5 kcal/mol) and SPB25 (-71.6 ± 0.5 kcal/mol). Additionally, they were predicted to have slightly higher stabilities, based on their percent helicities in water, than SBP1 (the experimentally proven inhibitor of SARS-CoV-2-RBD). Our two best designed SPB25s are promising candidates as omicron variant inhibitors.

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.

At the crossroads of epidemiology and biology: Bridging the gap between SARS-CoV-2 viral strain properties and epidemic wave characteristics

The COVID-19 pandemic has given rise to numerous articles from different scientific fields (epidemiology, virology, immunology, airflow physics …) without any effort to link these different insights. In this review, we aim to establish relationships between epidemiological data and the characteristics of the virus strain responsible for the epidemic wave concerned. We have carried out this study on the Wuhan, Alpha, Delta and Omicron strains allowing us to illustrate the evolution of the relationships we have highlighted according to these different viral strains. We addressed the following questions. 1) How can the mean infectious dose (one quantum, by definition in epidemiology) be measured and expressed as an amount of viral RNA molecules (in genome units, GU) or as a number of replicative viral particles (in plaque-forming units, PFU)? 2) How many infectious quanta are exhaled by an infected person per unit of time? 3) How many infectious quanta are exhaled, on average, integrated over the whole contagious period? 4) How do these quantities relate to the epidemic reproduction rate R as measured in epidemiology, and to the viral load, as measured by molecular biological methods? 5) How has the infectious dose evolved with the different strains of SARS-CoV-2? We make use of state-of-the-art modelling, reviewed and explained in the appendix of the article (Supplemental Information, SI), to answer these questions using data from the literature in both epidemiology and virology. We have considered the modification of these relationships according to the vaccination status of the population.

BCEDB: a linear B-cell epitopes database for SARS-CoV-2

The 2019 Novel Coronavirus (SARS-CoV-2) has infected millions of people worldwide and caused millions of deaths. The virus has gone numerous mutations to replicate faster, which can overwhelm the immune system of the host. Linear B-cell epitopes are becoming promising in prevention of various deadly infectious diseases, breaking the general idea of their low immunogenicity and partial protection. However, there is still no public repository to host the linear B-cell epitopes for facilitating the development vaccines against SARS-CoV-2. Therefore, we developed BCEDB, a linear B-cell epitopes database specifically designed for hosting, exploring and visualizing linear B-cell epitopes and their features. The database provides a comprehensive repository of computationally predicted linear B-cell epitopes from Spike protein; a systematic annotation of epitopes including sequence, antigenicity score, genomic locations of epitopes, mutations in different virus lineages, mutation sites on the 3D structure of Spike protein and a genome browser to visualize them in an interactive manner. It represents a valuable resource for peptide-based vaccine development. Database URL: http://www.oncoimmunobank.cn/bcedbindex.

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.

Computational redesign of Beta-27 Fab with substantially better predicted binding affinity to the SARS-CoV-2 Omicron variant than human ACE2 receptor

During the COVID-19 pandemic, SARS-CoV-2 has caused large numbers of morbidity and mortality, and the Omicron variant (B.1.1.529) was an important variant of concern. To enter human cells, the receptor-binding domain (RBD) of the S1 subunit of SARS-CoV-2 (SARS-CoV-2-RBD) binds to the peptidase domain (PD) of Angiotensin-converting enzyme 2 (ACE2) receptor. Disrupting the binding interactions between SARS-CoV-2-RBD and ACE2-PD using neutralizing antibodies is an effective COVID-19 therapeutic solution. Previous study found that Beta-27 Fab, which was obtained by digesting the full IgG antibodies that were isolated from a patient infected with SARS-CoV-2 Beta variant, can neutralize Victoria, Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), and Delta (B.1.617.2) variants. This study employed computational protein design and molecular dynamics (MD) to investigate and enhance the binding affinity of Beta-27 Fab to SARS-CoV-2-RBD Omicron variant. MD results show that five best designed Beta-27 Fabs (Beta-27-D01 Fab, Beta-27-D03 Fab, Beta-27-D06 Fab, Beta-27-D09 Fab and Beta-27-D10 Fab) were predicted to bind to Omicron RBD in the area, where ACE2 binds, with significantly better binding affinities than Beta-27 Fab and ACE2. Their enhanced binding affinities are mostly caused by increased binding interactions of CDR L2 and L3. They are promising candidates that could potentially be employed to disrupt the binding between ACE2 and Omicron RBD.

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

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

Bench-to-bedside: Innovation of small molecule anti-SARS-CoV-2 drugs in China

The ongoing COVID-19 pandemic has resulted in millions of deaths globally, highlighting the need to develop potent prophylactic and therapeutic strategies against SARS-CoV-2. Small molecule inhibitors (remdesivir, Paxlovid, and molnupiravir) are essential complements to vaccines and play important roles in clinical treatment of SARS-CoV-2. Many advances have been made in development of anti-SARS-CoV-2 inhibitors in China, but progress in discovery and characterization of pharmacological activity, antiviral mechanisms, and clinical efficacy are limited. We review development of small molecule anti-SARS-CoV-2 drugs (azvudine [approved by the NMPA of China on July 25, 2022], VV116 [approved by the NMPA of China on January 29, 2023], FB2001, WPV01, pentarlandir, and cepharanthine) in China and summarize their pharmacological activity, potential mechanisms of action, clinical trials and use, and important milestones in their discovery. The role of structural biology in drug development is also reviewed. Future studies should focus on development of diverse second-generation inhibitors with excellent oral bioavailability, superior plasma half-life, increased antiviral activity against SARS-CoV-2 and its variants, high target specificity, minimal side effects, reduced drug-drug interactions, and improved lung histopathology.

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.

Diverse SARS-CoV-2 aptamers overcome variant antigenic shift

SARS-CoV-2 mutates rapidly as evidenced by the emergence of Omicron which causes changes in the recognition epitopes of most current neutralizing antibodies and immune evasion. Although aptamers are potential neutralizing agents for SARS-CoV-2 due to their unique molecular properties, it is difficult to compare their performances as assay conditions vary greatly, and their activity levels against variants remain unknown. Here, we evaluated the performances of 14 SARS-CoV-2 aptamers and provided a comprehensive analysis them, which we expect will improve the development of aptamer tools for SARS-CoV-2 diagnostics and therapeutics.

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.

A Comprehensive Technology Platform for the Rapid Discovery of Peptide Inhibitors against SARS-CoV-2 Pseudovirus Infection

We developed and validated a technology platform for designing and testing peptides inhibiting the infectivity of SARS-CoV-2 spike protein-based pseudoviruses. This platform integrates target evaluation, in silico inhibitor design, peptide synthesis, and efficacy screening. We generated a cyclic peptide library derived from the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein and the angiotensin-converting enzyme 2 (ACE2) receptor. The cell-free validation process by ELISA competition assays and Surface Plasmon Resonance (SPR) studies revealed that the cyclic peptide c9_05, but not its linear form, binds well to ACE2. Moreover, it effectively inhibited the transduction in HEK293, stably expressing the human ACE2 receptor of pseudovirus particles displaying the SARS-CoV-2 spike in the Wuhan or UK variants. However, the inhibitory efficacy of c9_05 was negligible against the Omicron variant, and it failed to impede the entry of pseudoviruses carrying the B.1.351 (South African) spike. These variants contain three or more mutations known to increase affinity to ACE2. This suggests further refinement is needed for potential SARS-CoV-2 inhibition. Our study hints at a promising approach to develop inhibitors targeting viral infectivity receptors, including SARS-CoV-2's. This platform also promises swift identification and evaluation of inhibitors for other emergent viruses.

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

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

The effects of amino acid substitution of spike protein and genomic recombination on the evolution of SARS-CoV-2

Over three years' pandemic of 2019 novel coronavirus disease (COVID-19), multiple variants and novel subvariants have emerged successively, outcompeted earlier variants and become predominant. The sequential emergence of variants reflects the evolutionary process of mutation-selection-adaption of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Amino acid substitution/insertion/deletion in the spike protein causes altered viral antigenicity, transmissibility, and pathogenicity of SARS-CoV-2. Early in the pandemic, D614G mutation conferred virus with advantages over previous variants and increased transmissibility, and it also laid a conservative background for subsequent substantial mutations. The role of genomic recombination in the evolution of SARS-CoV-2 raised increasing concern with the occurrence of novel recombinants such as Deltacron, XBB.1.5, XBB.1.9.1, and XBB.1.16 in the late phase of pandemic. Co-circulation of different variants and co-infection in immunocompromised patients accelerate the emergence of recombinants. Surveillance for SARS-CoV-2 genomic variations, particularly spike protein mutation and recombination, is essential to identify ongoing changes in the viral genome and antigenic epitopes and thus leads to the development of new vaccine strategies and interventions.

S373P Mutation Stabilizes the Receptor-Binding Domain of the Spike Protein in Omicron and Promotes Binding

A cluster of several newly occurring mutations on Omicron is found at the β-core region of the spike protein's receptor-binding domain (RBD), where mutation rarely happened before. Notably, the binding of SARS-CoV-2 to human receptor ACE2 via RBD happens in a dynamic airway environment, where mechanical force caused by coughing or sneezing occurs. Thus, we used atomic force microscopy-based single-molecule force spectroscopy (AFM-SMFS) to measure the stability of RBDs and found that the mechanical stability of Omicron RBD increased by ∼20% compared with the wild type. Molecular dynamics (MD) simulations revealed that Omicron RBD showed more hydrogen bonds in the β-core region due to the closing of the α-helical motif caused primarily by the S373P mutation. In addition to a higher unfolding force, we showed a higher dissociation force between Omicron RBD and ACE2. This work reveals the mechanically stabilizing effect of the conserved mutation S373P for Omicron and the possible evolution trend of the β-core region of RBD.

Structural basis for receptor binding and broader interspecies receptor recognition of currently circulating Omicron sub-variants

Multiple SARS-CoV-2 Omicron sub-variants, such as BA.2, BA.2.12.1, BA.4, and BA.5, emerge one after another. BA.5 has become the dominant strain worldwide. Additionally, BA.2.75 is significantly increasing in some countries. Exploring their receptor binding and interspecies transmission risk is urgently needed. Herein, we examine the binding capacities of human and other 28 animal ACE2 orthologs covering nine orders towards S proteins of these sub-variants. The binding affinities between hACE2 and these sub-variants remain in the range as that of previous variants of concerns (VOCs) or interests (VOIs). Notably, R493Q reverse mutation enhances the bindings towards ACE2s from humans and many animals closely related to human life, suggesting an increased risk of cross-species transmission. Structures of S/hACE2 or RBD/hACE2 complexes for these sub-variants and BA.2 S binding to ACE2 of mouse, rat or golden hamster are determined to reveal the molecular basis for receptor binding and broader interspecies recognition.