Structural basis for SARS-CoV-2 Delta variant recognition of ACE2 receptor and broadly neutralizing antibodies

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

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

Tracking vaccine effectiveness in an evolving pandemic, countering misleading hot takes and epidemiologic fallacies

With the emergence of Omicron during the pandemic and the establishment of antibody waning over time, vaccine effectiveness, especially against infection, declined sharply from the original levels seen after the initial rollout. However, studies have demonstrated that they still provided substantial protection vs severe/fatal disease even with Omicron and after waning. Social media has been rife with reports claiming vaccines provided no benefit and some even claiming they made things worse, often driven by simple presentations of raw observational data using erroneous arguments involving epidemiologic fallacies including the base rate fallacy, Simpson's paradox, and the ecological fallacy and ignoring the extensive bias especially from confounding that is an inherent feature of these data. Similar fallacious arguments have been made by some in promoting vaccination policies, as well. Generally, vaccine effectiveness cannot be accurately estimated from raw population summaries but instead require rigorous, careful studies using epidemiologic designs and statistical analysis tools attempting to adjust for key confounders and sources of bias. This article summarizes what aggregated evidence across studies reveals about effectiveness of the mRNA vaccines as the pandemic has evolved, chronologically summarized with emerging variants and highlighting some of the fallacies and flawed arguments feeding social media-based claims that have obscured society's collective understanding.

Subtle Changes at the RBD/hACE2 Interface During SARS-CoV-2 Variant Evolution: A Molecular Dynamics Study

The SARS-CoV-2 Omicron variants show different behavior compared to the previous variants, especially with respect to the Delta variant, which promotes a lower morbidity despite being much more contagious. In this perspective, we performed molecular dynamics (MD) simulations of the different spike RBD/hACE2 complexes corresponding to the WT, Delta and four Omicron variants. Carrying out a comprehensive analysis of residue interactions within and between the two partners allowed us to draw the profile of each variant by using complementary methods (PairInt, hydrophobic potential, contact PCA). PairInt calculations highlighted the residues most involved in electrostatic interactions, which make a strong contribution to the binding with highly stable interactions between spike RBD and hACE2. Apolar contacts made a substantial and complementary contribution in Omicron with the detection of two hydrophobic patches. Contact networks and cross-correlation matrices were able to detect subtle changes at point mutations as the S375F mutation occurring in all Omicron variants, which is likely to confer an advantage in binding stability. This study brings new highlights on the dynamic binding of spike RBD to hACE2, which may explain the final persistence of Omicron over Delta.

An ACE2-Fc decoy produced in glycoengineered plants neutralizes ancestral and newly emerging SARS-CoV-2 variants and demonstrates therapeutic efficacy in hamsters

Newly emerging SARS-CoV-2 variants of concern (VOCs) continue to drive COVID-19 waves and are typically associated with immune escape and increased resistance to current therapeutics including monoclonal antibodies. By contrast, VOCs still display strong binding to the host cell receptor ACE2. Consistent with these properties, we have now found that a soluble ACE2-Fc decoy produced in glycoengineered plants effectively neutralizes different SARS-CoV-2 isolates and exhibits even increased potency against VOCs as compared to an ancestral virus strain. In a golden Syrian hamster model, therapeutic intranasal delivery of ACE2-Fc effectively reduced weight loss and SARS-CoV-2 replication in the lungs when administered 24 h post-inoculation. This protective effect was not observed upon treatment of the infected animals with a non-binding ACE2-Fc mutant, demonstrating that the plant-derived ACE2-Fc decoy interferes specifically with the attachment of the virus to host cells. The results obtained provide support for further development of decoy-based antiviral approaches by plant molecular pharming.

Structure and dynamics of Alpha B.1.1.7 SARS-CoV-2 S-protein in complex with Fab of neutralizing antibody REGN10987

One of the approaches for treatment of COVID-19 is a use of neutralizing antibodies (nAbs). The study of the mechanisms by which nAbs recognize different strains of SARS-CoV-2 may facilitate the development of new drugs and vaccines against the coronavirus infection. In this work, we present the 3.1 Å resolution cryo-electron microscopy structure of a full-length trimeric spike-protein (S-protein) of the SARS-CoV-2 Alpha (B.1.1.7) variant in complex with the Fab of the REGN10987 nAb. In the complex, two receptor-binding domains (RBDs) of the S-protein were observed in the 'up' state, whereas third RBD was in the 'down' state. This distinguishes the obtained structure from the complex of Delta (B.1.617.2) S-protein with REGN10987-Fab, where only one RBD was in the 'up' state. Probably some of the substituted residues (K478T, A570D, and S982A) located at the interprotomer interfaces are responsible for the greater Alpha S-protein opening upon the REGN10987-Fab binding. The Fab identically binds to the RBD in the both 'up' and 'down' conformations. The RBD-Fab interaction interface was refined to a resolution of 3.6 Å. The antibody binds to the receptor-binding motif (RBM), which prevents the S-protein from the binding to its receptor, angiotensin-converting enzyme 2 (ACE-2). Comparison with the structures of the Wuhan (wild type) and Delta RBD variants in complex with REGN10987-Fab revealed that the N501Y and T478K/L452R mutations presented in the RBM of the Alpha and Delta variants, respectively, do not affect the mode of the RBD-Fab interaction.

SARS-CoV-2: The Interplay Between Evolution and Host Immunity

The persistence of SARS-CoV-2 infections at a global level reflects the repeated emergence of variant strains encoding unique constellations of mutations. These variants have been generated principally because of a dynamic host immune landscape, the countermeasures deployed to combat disease, and selection for enhanced infection of the upper airway and respiratory transmission. The resulting viral diversity creates a challenge for vaccination efforts to maintain efficacy, especially regarding humoral aspects of protection. Here, we review our understanding of how SARS-CoV-2 has evolved during the pandemic, the immune mechanisms that confer protection, and the impact viral evolution has had on transmissibility and adaptive immunity elicited by natural infection and/or vaccination. Evidence suggests that SARS-CoV-2 evolution initially selected variants with increased transmissibility but currently is driven by immune escape. The virus likely will continue to drift to maintain fitness until countermeasures capable of disrupting transmission cycles become widely available.

Allosteric Control and Glycan Shielding Adaptations in the SARS-CoV-2 Spike from Early to Peak Virulence

The SARS-CoV-2 Spike glycoprotein is central to viral infectivity and immune evasion, making it a key target for vaccine and therapeutic design. This trimeric peplomer undergoes dynamic conformational changes, particularly in its Receptor Binding Domain (RBD), which transitions between closed (down) and ACE2-accessible (up) states relative to the rest of the protein, to facilitate host cell entry. Structural understanding of such critical inter-domain motions, as well as epitope exposure quantification, is essential for obtaining an effective molecular handle over this protein and, in turn, exploiting it towards improved immunogen development. Focusing on the early circulating D614G form and the later emerging Delta (B.1.617.2) variant with higher virulence, we performed large-scale molecular dynamics simulations of the soluble form of the Spike in both 'down' and 'up' conformations of the RBD. Guided by differences in overall fluctuations, we described reaction coordinates based on domain rotations and tilting to extract features that distinguish D614G versus Delta structural behavior of the N-terminal Domain (NTD) and RBD. Using reaction coordinate analysis and Principal Component Analysis (PCA), we identify allosteric coupling between the N-terminal Domain (NTD) and RBD, where NTD tilting influences RBD gating. While some of these motions are conserved across variants, Delta exhibits an optimized RBD-gating mechanism that enhances ACE2 accessibility. Additionally, glycan remodeling in Delta enhances shielding at the NTD supersite, contributing to reduced sensitivity to neutralizing antibodies. Finally, we uncover the impact of the D950N mutation in the HR1 region, which modulates downstream Spike dynamics and immune evasion. Together, our findings reveal variant-specific and conserved structural determinants of SARS-CoV-2 Spike function, providing a mechanistic basis for allosteric modulation, glycan-mediated immune evasion, and viral adaptation. These insights offer valuable guidance for rational vaccine and therapeutic design against SARS-CoV-2 and emerging variants.

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

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

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

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

Structural insights into ACE2 interactions and immune activation of SARS-CoV-2 and its variants: an in-silico study

The initial interaction between COVID-19 and the human body involves the receptor-binding domain (RBD) of the viral spike protein with the angiotensin-converting enzyme 2 (ACE2) receptor. Likewise, the spike protein can engage with immune-related proteins, such as toll-like receptors (TLRs) and pulmonary surfactant proteins A (SP-A) and D (SP-D), thereby triggering immune responses. In this study, we utilize computational methods to investigate the interactions between the spike protein and TLRs (specifically TLR2 and TLR4), as well as (SP-A) and (SP-D). The study is conducted on four variants of concern (VOC) to differentiate and identify common virus behaviours. An assessment of the structural stability of various variants indicates slight changes attributed to mutations, yet overall structural integrity remains preserved. Our findings reveal the spike protein's ability to bind with TLR4 and TLR2, prompting immune activation. In addition, our in-silico results reveal almost similar docking scores and therefore affinity for both ACE2-spike and TLR4-spike complexes. We demonstrate that even minor changes due to mutations in all variants, surfactant A and D proteins can function as inhibitors against the spike in all variants, hindering the ACE2-RBD interaction.Communicated by Ramaswamy H. Sarma.

Ultrasensitive detection of intact SARS-CoV-2 particles in complex biofluids using microfluidic affinity capture

Measuring virus in biofluids is complicated by confounding biomolecules coisolated with viral nucleic acids. To address this, we developed an affinity-based microfluidic device for specific capture of intact severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our approach used an engineered angiotensin-converting enzyme 2 to capture intact virus from plasma and other complex biofluids. Our device leverages a staggered herringbone pattern, nanoparticle surface coating, and processing conditions to achieve detection of as few as 3 viral copies per milliliter. We further validated our microfluidic assay on 103 plasma, 36 saliva, and 29 stool samples collected from unique patients with COVID-19, showing SARS-CoV-2 detection in 72% of plasma samples. Longitudinal monitoring in the plasma revealed our device's capacity for ultrasensitive detection of active viral infections over time. Our technology can be adapted to target other viruses using relevant cell entry molecules for affinity capture. This versatility underscores the potential for widespread application in viral load monitoring and disease management.

Differential Impact of Spike Protein Mutations on SARS-CoV-2 Infectivity and Immune Evasion: Insights from Delta and Kappa Variants

SARS-CoV-2 continues to pose a global health challenge due to its high transmissibility and mutability, with new variants emerging that potentially undermine vaccination and therapeutic efforts. Mutations in the spike protein, particularly in the receptor-binding domain (RBD), significantly influence viral transmissibility and immune escape. However, the complex interplay of these mutations and their combined effects on viral fitness remain to be analyzed. In this study, we investigated the functional impact of key mutations found in the Delta and Kappa variants of SARS-CoV-2. Using pseudovirus assays, we demonstrated that the T478K and L452R mutations characteristic of the Delta variant primarily enhance viral infectivity, with minimal effect on antibody-mediated neutralization. Conversely, the E484Q mutation of the Kappa variant, alone or in combination with L452R, significantly improved evasion of antibody-mediated neutralization but appeared to compromise viral fitness and infectivity. Notably, contrary to previous reports, we found that the P681R mutation contributed neither to increased infectivity nor immune evasion at least in the assay system employed in this study. Our findings suggest that the Delta variant's global dominance over the Kappa variant may be attributed to its superior infectivity and transmissibility rather than enhanced immune evasion capabilities. These results provide valuable insights into the functional consequences of spike protein mutations and may aid in predicting the emergence and spread of future SARS-CoV-2 variants. Such understanding is crucial for enhancing public health preparedness and informing the development of next-generation vaccines and therapeutics.

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

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

Berbamine prevents SARS-CoV-2 entry and transmission

Effective antiviral drugs are essential to combat COVID-19 and future pandemics. Although many compounds show antiviral in vitro activity, only a few retain effectiveness in vivo against SARS-CoV-2. Here, we show that berbamine (Berb) is effective against SARS-CoV, MER-CoV, SARS-CoV-2 and its variants, including the XBB.1.16 variant. In hACE2.Tg mice, Berb suppresses SARS-CoV-2 replication through two distinct mechanisms: inhibiting spike-mediated viral entry and enhancing antiviral gene expression during infection. The administration of Berb, in combination with remdesivir (RDV), clofazimine (Clof) and fangchinoline (Fcn), nearly eliminated viral load and promoted recovery from acute SARS-CoV-2 infection and its variants. Co-housed mice in direct contact with either pre-treated or untreated infected mice exhibited negligible viral loads, reduced lung pathology, and decreased viral shedding, suggesting that Berb may effectively hinder virus transmission. This broad-spectrum activity positions Berb as a promising preventive or therapeutic option against betacoronaviruses.

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

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

Ongoing evolution of SARS-CoV-2 drives escape from mRNA vaccine-induced humoral immunity

With the onset of the COVID-19 pandemic 4 years ago, viral sequencing continues to document numerous individual mutations in the viral spike protein across many variants. To determine the ability of vaccine-mediated humoral immunity to combat continued SARS-CoV-2 evolution, we construct a comprehensive panel of pseudoviruses harboring each individual mutation spanning 4 years of the pandemic to understand the fitness cost and resistance benefits of each. These efforts identify numerous mutations that escape from vaccine-induced humoral immunity. Across 50 variants and 131 mutants we construct, we observe progressive loss of neutralization across variants, irrespective of vaccine doses, as well as increasing infectivity and ACE2 binding. Importantly, the recent XBB.1.5 booster significantly increases titers against most variants but not JN.1, KP.2, or KP.3. These findings demonstrate that variants continue to evade updated mRNA vaccines, highlighting the need for different approaches to control SARS-CoV-2 transmission.

Key β1-4 galactosylated glycan receptors of SARS-CoV-2 and its inhibitor from the galactosylated glycoproteins of bovine milk

INTRODUCTION

The binding of the spike (S) protein of SARS-CoV-2 to angiotensin-converting enzyme 2 (ACE2) is a critical stage in the process of infection. While previous studies indicated that the S protein and ACE2 are extensively glycosylated, the functions of glycans in their interactions remain uncertain.

OBJECTIVES

This study aimed to investigate the glycan receptors of SARS-CoV-2 and evaluate the inhibitory effects of galactosylated glycoproteins derived from bovine milk on the attachment of SARS-CoV-2 pseudovirus.

METHODS

An antibody-overlay lectin microarray was used to profile the glycopatterns of the S protein-S1 of SARS-CoV-2 and ACE2. Molecular dynamics simulation was used to mimic the interaction between the S protein and ACE2. The effects of N-glycans and β1-4 galactosylation on the interactions between SARS-CoV-2, its variations (B1.617.2 (Delta) and B1.1.529 (Omicron)), and ACE2 was assessed using molecular docking simulation and protein microarrays. The impact of glycoproteins (specifically sialylated glycoproteins or de-sialylated glycoproteins) derived from bovine milk on the interaction between S1 and ACE2, as well as on pseudoviral attachment and entry, was assessed using protein microarrays and pseudovirus-based microneutralization assays.

RESULTS

Our findings indicated that the galactosylated glycoforms were the most prevalent for both S1 and ACE2. Importantly, we demonstrated that the β1-4 galactosylated N-glycans of ACE2 played a crucial role in the binding of S1 of SARS-CoV-2 and its variations to ACE2. The glycoproteins derived from bovine milk had a large amount of galactosylated glycans, which are comparable to the glycoforms of ACE2. The glycoproteins effectively blocked the attachment and entry of the SARS-CoV-2 pseudovirus by competitively blocking the binding of S1 to ACE2.

CONCLUSIONS

Our findings demonstrated that the β1-4 galactosylated N-glycans of ACE2 play a crucial role as glycan receptors for the binding of S1 of SARS-CoV-2 and its variations. Moreover, the glycoproteins with 'receptor-like' glycoforms could be an effective inhibitor to prevent SARS-CoV-2 infection.

Dynamic expedition of leading mutations in SARS-CoV-2 spike glycoproteins

The continuous evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which caused the recent pandemic, has generated countless new variants with varying fitness. Mutations of the spike glycoprotein play a particularly vital role in shaping its evolutionary trajectory, as they have the capability to alter its infectivity and antigenicity. We present a time-resolved statistical method, Dynamic Expedition of Leading Mutations (deLemus), to analyze the evolutionary dynamics of the SARS-CoV-2 spike glycoprotein. The proposed L -index of the deLemus method is effective in quantifying the mutation strength of each amino acid site and outlining evolutionarily significant sites, allowing the comprehensive characterization of the evolutionary mutation pattern of the spike glycoprotein.

Comparison of the Neutralization Power of Sotrovimab Against SARS-CoV-2 Variants: Development of a Rapid Computational Method

BACKGROUND

The rapid evolution of SARS-CoV-2 imposed a huge challenge on disease control. Immune evasion caused by genetic variations of the SARS-CoV-2 spike protein's immunogenic epitopes affects the efficiency of monoclonal antibody-based therapy of COVID-19. Therefore, a rapid method is needed to evaluate the efficacy of the available monoclonal antibodies against the new emerging variants or potential novel variants.

OBJECTIVE

The aim of this study is to develop a rapid computational method to evaluate the neutralization power of anti-SARS-CoV-2 monoclonal antibodies against new SARS-CoV-2 variants and other potential new mutations.

METHODS

The amino acid sequence of the extracellular domain of the spike proteins of the severe acute respiratory syndrome coronavirus (GenBank accession number YP_009825051.1) and SARS-CoV-2 (GenBank accession number YP_009724390.1) were used to create computational 3D models for the native spike proteins. Specific mutations were introduced to the curated sequence to generate the different variant spike models. The neutralization potential of sotrovimab (S309) against these variants was evaluated based on its molecular interactions and Gibbs free energy in comparison to a reference model after molecular replacement of the reference receptor-binding domain with the variant's receptor-binding domain.

RESULTS

Our results show a loss in the binding affinity of the neutralizing antibody S309 with both SARS-CoV and SARS-CoV-2. The binding affinity of S309 was greater to the Alpha, Beta, Gamma, and Kappa variants than to the original Wuhan strain of SARS-CoV-2. However, S309 showed a substantially decreased binding affinity to the Delta and Omicron variants. Based on the mutational profile of Omicron subvariants, our data describe the effect of the G339H and G339D mutations and their role in escaping antibody neutralization, which is in line with published clinical reports.

CONCLUSIONS

This method is rapid, applicable, and of interest to adapt the use of therapeutic antibodies to the treatment of emerging variants. It could be applied to antibody-based treatment of other viral infections.

In Silico Design of miniACE2 Decoys with In Vitro Enhanced Neutralization Activity against SARS-CoV-2, Encompassing Omicron Subvariants

The COVID-19 pandemic has overwhelmed healthcare systems and triggered global economic downturns. While vaccines have reduced the lethality rate of SARS-CoV-2 to 0.9% as of October 2024, the continuous evolution of variants remains a significant public health challenge. Next-generation medical therapies offer hope in addressing this threat, especially for immunocompromised individuals who experience prolonged infections and severe illnesses, contributing to viral evolution. These cases increase the risk of new variants emerging. This study explores miniACE2 decoys as a novel strategy to counteract SARS-CoV-2 variants. Using in silico design and molecular dynamics, blocking proteins (BPs) were developed with stronger binding affinity for the receptor-binding domain of multiple variants than naturally soluble human ACE2. The BPs were expressed in E. coli and tested in vitro, showing promising neutralizing effects. Notably, miniACE2 BP9 exhibited an average IC50 of 4.9 µg/mL across several variants, including the Wuhan strain, Mu, Omicron BA.1, and BA.2 This low IC50 demonstrates the potent neutralizing ability of BP9, indicating its efficacy at low concentrations.Based on these findings, BP9 has emerged as a promising therapeutic candidate for combating SARS-CoV-2 and its evolving variants, thereby positioning it as a potential emergency biopharmaceutical.

Complex vaccination strategies prevent the emergence of vaccine resistance

Vaccination is the most effective tool to control infectious diseases. However, the evolution of vaccine resistance, exemplified by vaccine resistance in SARS-CoV-2, remains a concern. Here, we model complex vaccination strategies against a pathogen with multiple epitopes-molecules targeted by the vaccine. We found that a vaccine targeting one epitope was ineffective in preventing vaccine escape. Vaccine resistance in highly infectious pathogens was prevented by the full-epitope vaccine, that is, one targeting all available epitopes, but only when the rate of pathogen evolution was low. Strikingly, a bet-hedging strategy of random administration of vaccines targeting different epitopes was the most effective in preventing vaccine resistance in pathogens with the low rate of infection and high rate of evolution. Thus, complex vaccination strategies, when biologically feasible, may be preferable to the currently used single-vaccine approaches for long-term control of disease outbreaks, especially when applied to livestock with near 100% vaccination rates.

The comparison of pathogenicity among SARS-CoV-2 variants in domestic cats

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been detected or isolated from domestic cats. It is unclear whether cats play an important role in the SARS-CoV-2 transmission cycle. In this study, we examined the susceptibility of cats to SARS-CoV-2, including wild type and variants, by animal experiments. Cats inoculated with wild type, gamma, and delta variants secreted a large amount of SARS-CoV-2 for 1 week after the inoculation from nasal, oropharyngeal, and rectal routes. Only 100 TCID50 of virus could infect cats and replicate well without severe clinical symptoms. In addition, one cat inoculated with wild type showed persistent virus secretion in feces for over 28 days post-inoculation (dpi). The titer of virus-neutralizing (VN) antibodies against SARS-CoV-2 increased from 11 dpi, reaching a peak at 14 dpi. However, the omicron variant could not replicate well in cat tissues and induced a lower titer of VN antibodies. It is concluded that cats were highly susceptible to SARS-CoV-2 infection, but not to the Omicron Variant, which caused the attenuated pathogenicity.

Structural review of SARS-CoV-2 antiviral targets

The coronavirus disease 2019 (COVID-19), the disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), represents the most disastrous infectious disease pandemic of the past century. As a member of the Betacoronavirus genus, the SARS-CoV-2 genome encodes a total of 29 proteins. The spike protein, RNA-dependent RNA polymerase, and proteases play crucial roles in the virus replication process and are promising targets for drug development. In recent years, structural studies of these viral proteins and of their complexes with antibodies and inhibitors have provided valuable insights into their functions and laid a solid foundation for drug development. In this review, we summarize the structural features of these proteins and discuss recent progress in research regarding therapeutic development, highlighting mechanistically representative molecules and those that have already been approved or are under clinical investigation.

Visceral Adiposity and Neutralizing Antibody Expression: An Adult-Based Cross-Sectional Study

Purpose

Visceral adiposity is a significant risk factor for severe COVID-19. However, the impact of the Chinese visceral adiposity index (CVAI) on the efficacy of SARS-CoV-2 vaccines remains poorly understood. This study aims to explore the impact of CVAI on the production of neutralizing antibodies (NAb) in inactivated SARS-CoV-2 vaccines and the potential mechanism, thereby optimizing vaccination guidance.

Methods

In this cross-sectional study, 206 health workers (completed two SARS-CoV-2 vaccination on February 8th and March 10th, 2021, respectively) were recruited. All baseline anthropometric parameters of the participants were collected, and venous blood samples were obtained 6 weeks later to measure peripheral innate immune cells, inflammatory cytokines, and NAb titers against SARS-CoV-2. CVAI were calculated according to the formula and divided participants into two groups depending on CVAI median.

Results

The median NAb titer among healthcare workers was 12.94 AU/mL, with an efficacy of 87.86% for the SARS-CoV-2 vaccine. NAb titers were lower in the CVAI dysfunction group than in the CVAI reference group (median: 11.40 AU/mL vs 15.57 AU/mL), the hsCRP levels (median: 0.50 mg/L vs 0.30 mg/L) and peripheral monocyte count (mean: 0.47 × 109/L vs 0.42 × 109/L) in the CVAI dysfunction group were higher than in the CVAI reference group. Additionally, CVAI showed positive correlations with hsCRP, monocytes, lymphocytes, and B-lymphocytes, and a negative correlation with NAb titers.

Conclusion

CVAI may inhibit SARS-CoV-2 neutralizing antibody expression through inducing immune dysfunction and chronic inflammation. Thus, more attention should be paid to the vaccination for high CVAI population to improve the effectiveness of vaccination, which could provide more robust support for COVID-19 epidemic prevention and control.

Adsorption-Driven Deformation and Footprints of the RBD Proteins in SARS-CoV-2 Variants on Biological and Inanimate Surfaces

Respiratory viruses, carried through airborne microdroplets, frequently adhere to surfaces, including plastics and metals. However, our understanding of the interactions between viruses and materials remains limited, particularly in scenarios involving polarizable surfaces. Here, we investigate the role of the receptor-binding domain (RBD) of the spike protein mutations on the adsorption of SARS-CoV-2 to hydrophobic and hydrophilic surfaces employing molecular simulations. To contextualize our findings, we contrast the interactions on inanimate surfaces with those on native biological interfaces, specifically the angiotensin-converting enzyme 2. Notably, we identify a 2-fold increase in structural deformations for the protein's receptor binding motif (RBM) onto inanimate surfaces, indicative of enhanced shock-absorbing mechanisms. Furthermore, the distribution of adsorbed amino acids (landing footprints) on the inanimate surface reveals a distinct regional asymmetry relative to the biological interface, with roughly half of the adsorbed amino acids arranged in opposite sites. In spite of the H-bonds formed at the hydrophilic substrate, the simulations consistently show a higher number of contacts and interfacial area with the hydrophobic surface, where the wild-type RBD adsorbs more strongly than the Delta or Omicron RBDs. In contrast, the adsorption of Delta and Omicron to hydrophilic surfaces was characterized by a distinctive hopping-pattern. The novel shock-absorbing mechanisms identified in the virus adsorption on inanimate surfaces show the embedded high-deformation capacity of the RBD without losing its secondary structure, which could lead to current experimental strategies in the design of virucidal surfaces.

Structural insights into immune escape at killer T cell epitope by SARS-CoV-2 Spike Y453F variants

CD8+ T cell immunity, mediated by human leukocyte antigen (HLA) and T cell receptor (TCR), plays a critical role in conferring immune memory and protection against viral pathogens. The emergence of SARS-CoV-2 variants poses a serious challenge to the efficacy of current vaccines. Whereas numerous SARS-CoV-2 mutations associated with immune escape from CD8+ T cells have been documented, the molecular effects of most mutations on epitope-specific TCR recognition remain largely unexplored. Here, we studied an HLA-A24-restricted NYN epitope (Spike448-456) that elicits broad CD8+ T cell responses in COVID-19 patients characterized by a common TCR repertoire. Four natural mutations, N450K, L452Q, L452R, and Y453F, arose within the NYN epitope and have been transmitted in certain viral lineages. Our findings indicate that these mutations have minimal impact on the epitope's presentation by cell surface HLA, yet they diminish the affinities of their respective peptide-HLA complexes (pHLAs) for NYN peptide-specific TCRs, particularly L452R and Y453F. Furthermore, we determined the crystal structure of HLA-A24 loaded with the Y453F peptide (NYNYLFRLF), and subsequently a ternary structure of the public TCRNYN-I complexed to the original NYN-HLA-A24 (NYNYLYRLF). Our structural analysis unveiled that despite competent presentation by HLA, the mutant Y453F peptide failed to establish a stable TCR-pHLA ternary complex due to reduced peptide: TCR contacts. This study supports the idea that cellular immunity restriction is an important driving force behind viral evolution.

Single-molecule imaging reveals allosteric stimulation of SARS-CoV-2 spike receptor binding domain by host sialic acid

Conformational dynamics of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S) mediate exposure of the binding site for the cellular receptor, angiotensin-converting enzyme 2 (ACE2). The N-terminal domain (NTD) of S binds terminal sialic acid (SA) moieties on the cell surface, but the functional role of this interaction in virus entry is unknown. Here, we report that NTD-SA interaction enhances both S-mediated virus attachment and ACE2 binding. Through single-molecule Förster resonance energy transfer imaging of individual S trimers, we demonstrate that SA binding to the NTD allosterically shifts the S conformational equilibrium, favoring enhanced exposure of the ACE2-binding site. Antibodies that target the NTD block SA binding, which contributes to their mechanism of neutralization. These findings inform on mechanisms of S activation at the cell surface.

Exploring conformational landscapes and binding mechanisms of convergent evolution for the SARS-CoV-2 spike Omicron variant complexes with the ACE2 receptor using AlphaFold2-based structural ensembles and molecular dynamics simulations

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

Identifications of novel host cell factors that interact with the receptor-binding domain of the SARS-CoV-2 spike protein

SARS-CoV-2 entry into host cells is facilitated by the interaction between the receptor-binding domain of its spike protein (CoV2-RBD) and host cell receptor, ACE2, promoting viral membrane fusion. The virus also uses endocytic pathways for entry, but the mediating host factors remain largely unknown. It is also unknown whether mutations in the RBD of SARS-CoV-2 variants promote interactions with additional host factors to promote viral entry. Here, we used the GST pull-down approach to identify novel surface-located host factors that bind to CoV2-RBD. One of these factors, SH3BP4, regulates internalization of CoV2-RBD in an ACE2-independent but integrin- and clathrin-dependent manner and mediates SARS-CoV-2 pseudovirus entry, suggesting that SH3BP4 promotes viral entry via the endocytic route. Many of the identified factors, including SH3BP4, ADAM9, and TMEM2, show stronger affinity to CoV2-RBD than to RBD of the less infective SARS-CoV, suggesting SARS-CoV-2-specific utilization. We also found factors preferentially binding to the RBD of the SARS-CoV-2 Delta variant, potentially enhancing its entry. These data identify the repertoire of host cell surface factors that function in the events leading to the entry of SARS-CoV-2.

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.

Nanobody repertoire generated against the spike protein of ancestral SARS-CoV-2 remains efficacious against the rapidly evolving virus

To date, all major modes of monoclonal antibody therapy targeting SARS-CoV-2 have lost significant efficacy against the latest circulating variants. As SARS-CoV-2 omicron sublineages account for over 90% of COVID-19 infections, evasion of immune responses generated by vaccination or exposure to previous variants poses a significant challenge. A compelling new therapeutic strategy against SARS-CoV-2 is that of single-domain antibodies, termed nanobodies, which address certain limitations of monoclonal antibodies. Here, we demonstrate that our high-affinity nanobody repertoire, generated against wild-type SARS-CoV-2 spike protein (Mast et al., 2021), remains effective against variants of concern, including omicron BA.4/BA.5; a subset is predicted to counter resistance in emerging XBB and BQ.1.1 sublineages. Furthermore, we reveal the synergistic potential of nanobody cocktails in neutralizing emerging variants. Our study highlights the power of nanobody technology as a versatile therapeutic and diagnostic tool to combat rapidly evolving infectious diseases such as SARS-CoV-2.

The Recognition Pathway of the SARS-CoV-2 Spike Receptor-Binding Domain to Human Angiotensin-Converting Enzyme 2

COVID-19 caused by SARS-CoV-2 has spread around the world. The receptor-binding domain (RBD) of the spike protein of SARS-CoV-2 is a critical component that directly interacts with host ACE2. Here, we simulate the ACE2 recognition processes of RBD of the WT, Delta, and OmicronBA.2 variants using our recently developed supervised Gaussian accelerated molecular dynamics (Su-GaMD) approach. We show that RBD recognizes ACE2 through three contact regions (regions I, II, and III), which aligns well with the anchor-locker mechanism. The higher binding free energy in State d of the RBDOmicronBA.2-ACE2 system correlates well with the increased infectivity of OmicronBA.2 in comparison with other variants. For RBDDelta, the T478K mutation affects the first step of recognition, while the L452R mutation, through its nearby Y449, affects the RBDDelta-ACE2 binding in the last step of recognition. For RBDOmicronBA.2, the E484A mutation affects the first step of recognition, the Q493R, N501Y, and Y505H mutations affect the binding free energy in the last step of recognition, mutations in the contact regions affect the recognition directly, and other mutations indirectly affect recognition through dynamic correlations with the contact regions. These results provide theoretical insights for RBD-ACE2 recognition and may facilitate drug design against SARS-CoV-2.

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

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

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

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

SARS-CoV-2 variant with the spike protein mutation F306L in the southern border provinces of Thailand

The southernmost part of Thailand is a unique and culturally diverse region that has been greatly affected by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) outbreak during the coronavirus disease-2019 pandemic. To gain insights into this situation, we analyzed 1942 whole-genome sequences of SARS-CoV-2 obtained from the five southernmost provinces of Thailand between April 2021 and March 2022, together with those publicly available in the Global Initiative on Sharing All Influenza Data database. Our analysis revealed evidence for transboundary transmissions of the virus in and out of the five southernmost provinces during the study period, from both domestic and international sources. The most prevalent viral variant in our sequence dataset was the Delta B.1.617.2.85 variant, also known as the Delta AY.85 variant, with many samples carrying a non-synonymous mutation F306L in their spike protein. Protein-protein docking and binding interface analyses suggested that the mutation may enhance the binding between the spike protein and host cell receptor protein angiotensin-converting enzyme 2, and we found that the mutation was significantly associated with an increased fatality rate. This mutation has also been observed in other SARS-CoV-2 variants, suggesting that it is of particular interest and should be monitored.

Unraveling the binding mechanisms of SARS-CoV-2 variants through molecular simulations

The emergence of SARS-CoV-2 variants like Delta (AY.29) and Omicron (EG.5) poses continued challenges for vaccines and therapeutics. Mutations in the viral spike protein are key in altering infectivity and immune evasion. This study uses computational modeling to investigate the molecular binding mechanisms between spike protein variants and the ACE2 host receptor. Using the MARTNI force field, coarse-grained molecular dynamics (CGMD) simulations and nudged elastic band (NEB) calculations explore spike-ACE2 interactions for the wild type, Delta variant, and Omicron variant. The simulations reveal Omicron has the strongest binding affinity at -128.35 ± 10.91 kcal/mol, followed by Delta and wild type. Key mutations in Delta and Omicron, like Q493R and Q498R, optimize electrostatic contacts, enhancing ACE2 interactions. The wild-type spike has the highest transition state energy barrier at 17.87 kcal/mol, while Delta has the lowest barrier at 9.21 kcal/mol. Despite slightly higher dual barriers, Omicron's increased binding energy lowers its overall barrier to rapidly bind ACE2. These findings provide residue-level insights into mutation effects on SARS-CoV-2 infectivity. The computational modeling elucidates mechanisms underlying spike-ACE2 binding kinetics, aiding the development of vaccines and therapies targeting emerging viral strains.

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.

Synchronous Diagnosis of Respiratory Viruses Variants via Receptonics Based on Modeling Receptor-Ligand Dynamics

The transmission and pathogenesis of highly contagious fatal respiratory viruses are increasing, and the need for an on-site diagnostic platform has arisen as an issue worldwide. Furthermore, as the spread of respiratory viruses continues, different variants have become the dominant circulating strains. To prevent virus transmission, the development of highly sensitive and accurate on-site diagnostic assays is urgently needed. Herein, a facile diagnostic device is presented for multi-detection based on the results of detailed receptor-ligand dynamics simulations for the screening of various viral strains. The novel bioreceptor-treated electronics (receptonics) device consists of a multichannel graphene transistor and cell-entry receptors conjugated to N-heterocyclic carbene (NHC). An ultrasensitive multi-detection performance is achieved without the need for sample pretreatment, which will enable rapid diagnosis and prevent the spread of pathogens. This platform can be applied for the diagnosis of variants of concern in clinical respiratory virus samples and primate models. This multi-screening platform can be used to enhance surveillance and discriminate emerging virus variants before they become a severe threat to public health.

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.

Design of human ACE2 mimic miniprotein binders that interact with RBD of SARS-CoV-2 variants of concerns

The world of medicine demands from the research community solutions to the emerging problem of SARS-CoV-2 variants and other such potential global pandemics. With advantages of specificity over small molecule drugs and designability over antibodies, miniprotein therapeutics offers a unique solution to the threats of rapidly emerging SARS-CoV-2 variants. Unfortunately, most of the promising miniprotein binders are de novo designed and it is not viable to generate molecules for each new variant. Therefore in this study, we demonstrate a method for design of miniprotein mimics from the interaction interphase of human angiotensin converting enzyme 2 (ACE2). ACE2 is the natural interacting partner for the SARS-CoV-2 spike receptor binding domain (RBD) and acts as a recognition molecule for viral entry into the host cells. Starting with ACE2 N-terminal triple helix interaction interphase, we generated more than 70 miniprotein sequences. Employing Rosetta folding and docking scores we selected 10 promising miniprotein candidates amongst which 3 were found to be soluble in lab studies. Further, using molecular mechanics (MM) calculations on molecular dynamics (MD) trajectories we test interaction of miniproteins with RBD from various variants of concern (VOC). Presently, we report two key findings; miniproteins in this study are generated using less than 10 lab testing experiments, yet when tested through in-vitro experiments, they show submicro to nanomolar affinities towards SARS-CoV-2 RBD. Also in simulation studies, when compared with previously developed therapeutics, our miniproteins display remarkable ability to mimic ACE2 interphase; making them an ideal solution to the ever evolving problem of VOCs.

Nanobody repertoire generated against the spike protein of ancestral SARS-CoV-2 remains efficacious against the rapidly evolving virus

To date, all major modes of monoclonal antibody therapy targeting SARS-CoV-2 have lost significant efficacy against the latest circulating variants. As SARS-CoV-2 omicron sublineages account for over 90% of COVID-19 infections, evasion of immune responses generated by vaccination or exposure to previous variants poses a significant challenge. A compelling new therapeutic strategy against SARS-CoV-2 is that of single domain antibodies, termed nanobodies, which address certain limitations of monoclonal antibodies. Here we demonstrate that our high-affinity nanobody repertoire, generated against wild-type SARS-CoV-2 spike protein (Mast, Fridy et al. 2021), remains effective against variants of concern, including omicron BA.4/BA.5; a subset is predicted to counter resistance in emerging XBB and BQ.1.1 sublineages. Furthermore, we reveal the synergistic potential of nanobody cocktails in neutralizing emerging variants. Our study highlights the power of nanobody technology as a versatile therapeutic and diagnostic tool to combat rapidly evolving infectious diseases such as SARS-CoV-2.

A cardiotoxicity-eliminated ACE2 variant as a pan-inhibitor against coronavirus cell invasion

Human recombinant ACE2 (hrACE2) has been highly anticipated as a successful COVID-19 treatment; however, its potential to cause cardiac side effects has given rise to many concerns. Here, we developed a cardiotoxicity-eliminated hrACE2 variant, which had four mutation sites within hrACE2 (H345L, H374L, H378L, H505L) and was named as hrACE2-4mu. hrACE2-4mu has a consistent binding affinity with the variant SARS-CoV-2 spike proteins (SPs) and an efficient ability to block SP-induced SARS-CoV-2 entry into cells. In golden hamsters, injection of purified wild-type (WT) hrACE2 rescues the early stages of pneumonia caused by the SPs of the WT, delta, and omicron variants with reduced inflammatory cell infiltration. However, long-term injection of WT hrACE2 induces undesired cardiac fibrosis, as demonstrated by upregulated fibronectin and collagen expression. Our newly developed hrACE2-4mu showed similar protective abilities against a series of coronavirus cell invasions as WT hrACE2, meanwhile it did not cause apparent cardiac side effects. Thus, we generated a cardiotoxicity-eliminated variant of hrACE2 as a pan-inhibitor against coronavirus cell invasion, providing a potential novel strategy for the treatment of COVID-19 and other coronaviruses.

Dauricine interferes with SARS-CoV-2 variants infection by blocking the interface between RBD and ACE2

The continued viral evolution results in the emergence of various SARS-CoV-2 variants, such as delta or omicron, that are partially resistant to current vaccines and antiviral medicines, posing an increased risk to global public health and raising the importance of continuous development of antiviral medicines. Inhibitor screening targeting the interactions between the viral spike proteins and their human receptor ACE2 represents a promising approach for drug discovery. Here, we demonstrate that the evolutionary trend of the SARS-CoV-2 variants is associated with increased electrostatic interactions between S proteins and ACE2. Virtual screening based on the ACE2-RBD binding interface identified nine monomers of Traditional Chinese medicine (TCM). Furthermore, live-virus neutralization assays revealed that Dauricine, one of the identified monomers, exhibited an antiviral activity with an IC50 range of 18.2 to 33.3 μM for original strain, Delta, and Omicron strains, respectively. The computational study showed that the polycyclic and methoxy groups of Dauricine adhere to the RBD surface through π-π and electrostatic interactions. The discovery of Dauricine is a successful attempt to target viral entry, which will not only help society to respond quickly to viral variants, but also accelerate variant drug development thereby reducing the pressure on health authorities to respond to outbreaks.

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.

Genomic surveillance of SARS-CoV-2 reveals highest severity and mortality of delta over other variants: evidence from Cameroon

While the SARS-CoV-2 dynamic has been described globally, there is a lack of data from Sub-Saharan Africa. We herein report the dynamics of SARS-CoV-2 lineages from March 2020 to March 2022 in Cameroon. Of the 760 whole-genome sequences successfully generated by the national genomic surveillance network, 74% were viral sub-lineages of origin and non-variants of concern, 15% Delta, 6% Omicron, 3% Alpha and 2% Beta variants. The pandemic was driven by SARS-CoV-2 lineages of origin in wave 1 (16 weeks, 2.3% CFR), the Alpha and Beta variants in wave 2 (21 weeks, 1.6% CFR), Delta variants in wave 3 (11 weeks, 2.0% CFR), and omicron variants in wave 4 (8 weeks, 0.73% CFR), with a declining trend over time (p = 0.01208). Even though SARS-CoV-2 heterogeneity did not seemingly contribute to the breadth of transmission, the viral lineages of origin and especially the Delta variants appeared as drivers of COVID-19 severity in Cameroon.

Structural insights into broadly neutralizing antibodies elicited by hybrid immunity against SARS-CoV-2

ABSTRACTIncreasing spread by SARS-CoV-2 Omicron variants challenges existing vaccines and broadly reactive neutralizing antibodies (bNAbs) against COVID-19. Here we determine the diversity, potency, breadth and structural insights of bNAbs derived from memory B cells of BNT162b2-vaccinee after homogeneous Omicron BA.1 breakthrough infection. The infection activates diverse memory B cell clonotypes for generating potent class I/II and III bNAbs with new epitopes mapped to the receptor-binding domain (RBD). The top eight bNAbs neutralize wildtype and BA.1 potently but display divergent IgH/IgL sequences and neuralization profiles against other variants of concern (VOCs). Two of them (P2D9 and P3E6) belonging to class III NAbs display comparable potency against BA.4/BA.5, although structural analysis reveals distinct modes of action. P3E6 neutralizes all variants tested through a unique bivalent interaction with two RBDs. Our findings provide new insights into hybrid immunity on BNT162b2-induced diverse memory B cells in response to Omicron breakthrough infection for generating diverse bNAbs with distinct structural basis.

Fangchinoline inhibits SARS-CoV-2 and MERS-CoV entry

The COVID-19 pandemic caused by SARS-CoV-2, lead to mild to severe respiratory illness and resulted in 6.9 million deaths worldwide. Although vaccines are effective in preventing COVID-19, they may not be sufficient to protect immunocompromised individuals from this respiratory illness. Moreover, novel emerging variants of SARS-CoV-2 pose a risk of new COVID-19 waves. Therefore, identification of effective antivirals is critical in controlling SARS and other coronaviruses, such as MERS-CoV. We show that Fangchinoline (Fcn), a bisbenzylisoquinoline alkaloid, inhibits replication of SARS-CoV, SARS-CoV-2, and MERS-CoV in a range of in vitro assays, by blocking entry. Therapeutic use of Fcn inhibited viral loads in the lungs, and suppressed associated airway inflammation in hACE2. Tg mice and Syrian hamster infected with SARS-CoV-2. Combination of Fcn with remdesivir (RDV) or an anti-leprosy drug, Clofazimine, exhibited synergistic antiviral activity. Compared to Fcn, its synthetic derivative, MK-04-003, more effectively inhibited SARS-CoV-2 and its variants B.1.617.2 and BA.5 in mice. Taken together these data demonstrate that Fcn is a pan beta coronavirus inhibitor, which possibly can be used to combat novel emerging coronavirus diseases.

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.

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.