Impact of new variants on SARS-CoV-2 infectivity and neutralization: A molecular assessment of the alterations in the spike-host protein interactions

Antibody drugs targeting SARS-CoV-2: Time for a rethink?

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) heavily burdens human health. Multiple neutralizing antibodies (nAbs) have been issued for emergency use or tested for treating infected patients in the clinic. However, SARS-CoV-2 variants of concern (VOC) carrying mutations reduce the effectiveness of nAbs by preventing neutralization. Uncoding the mutation profile and immune evasion mechanism of SARS-CoV-2 can improve the outcome of Ab-mediated therapies. In this review, we first outline the development status of anti-SARS-CoV-2 Ab drugs and provide an overview of SARS-CoV-2 variants and their prevalence. We next focus on the failure causes of anti-SARS-CoV-2 Ab drugs and rethink the design strategy for developing new Ab drugs against COVID-19. This review provides updated information for the development of therapeutic Ab drugs against SARS-CoV-2 variants.

Molecular evolutionary characteristics of severe acute respiratory syndrome coronavirus 2 and the relatedness of epidemiological and socio-environmental factors

After the first outbreak, SARS-CoV-2 infection continues to occur due to the emergence of new variants. There is limited information available on the comparative evaluation of evolutionary characteristics of SARS-CoV-2 among different countries over time, and its relatedness to epidemiological and socio-environmental factors within those countries. We assessed comparative Bayesian evolutionary characteristics for SARS-CoV-2 in eight countries from 2020 to 2022 using BEAST version 2.6.7. Additionally, the relatedness between virus evolution factors and both epidemiological and socio-environmental factors was analyzed using Pearson's correlation coefficient. The estimated substitution rates in the gene encoding S protein of SARS-CoV-2 exhibited a continuous increase from 2020 to 2022 and were divided into two distinct groups in 2022 (p value < 0.05). Effective population size (Ne) generally showed decreased patterns by time. Notably, the change rates of the substitution rates were negatively correlated with the cumulative vaccination rates in 2021. A strict and rapid vaccination policy in the United Arab Emirates dramatically reduced the evolution of the virus, compared to other countries. Also, the average yearly temperature in countries were negatively correlated with the substitution rates. The changes of six epitopes in SARS-CoV-2 were related to various socio-environmental factors. We figured out comparative virus evolutionary traits and the association of epidemiological and socio-environmental factors especially cumulative vaccination rates and average temperature.

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.

Targeting the receptor binding domain and heparan sulfate binding for antiviral drug development against SARS-CoV-2 variants

The emergence of SARS-CoV-2 variants diminished the efficacy of current antiviral drugs and vaccines. Hence, identifying highly conserved sequences and potentially druggable pockets for drug development was a promising strategy against SARS-CoV-2 variants. In viral infection, the receptor-binding domain (RBD) proteins are essential in binding to the host receptor. Others, Heparan sulfate (HS), widely distributed on the surface of host cells, is thought to play a central role in the viral infection cycle of SARS-CoV-2. Therefore, it might be a reasonable strategy for antiviral drug design to interfere with the RBD in the HS binding site. In this study, we used computational approaches to analyze multiple sequences of coronaviruses and reveal important information about the binding of HS to RBD in the SARS-CoV-2 spike protein. Our results showed that the potential hot-spots, including R454 and E471, in RBD, exhibited strong interactions in the HS-RBD binding region. Therefore, we screened different compounds in the natural product database towards these hot-spots to find potential antiviral candidates using LibDock, Autodock vina and furthermore applying the MD simulation in AMBER20. The results showed three potential natural compounds, including Acetoside (ACE), Hyperoside (HYP), and Isoquercitrin (ISO), had a strong affinity to the RBD. Our results demonstrate a feasible approach to identify potential antiviral agents by evaluating the binding interaction between viral glycoproteins and host receptors. The present study provided the applications of the structure-based computational approach for designing and developing of new antiviral drugs against SARS-CoV-2 variants.

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 new paradigm for molecular dynamics databases: the COVID-19 database, the legacy of a titanic community effort

Molecular dynamics (MD) simulations are keeping computers busy around the world, generating a huge amount of data that is typically not open to the scientific community. Pioneering efforts to ensure the safety and reusability of MD data have been based on the use of simple databases providing a limited set of standard analyses on single-short trajectories. Despite their value, these databases do not offer a true solution for the current community of MD users, who want a flexible analysis pipeline and the possibility to address huge non-Markovian ensembles of large systems. Here we present a new paradigm for MD databases, resilient to large systems and long trajectories, and designed to be compatible with modern MD simulations. The data are offered to the community through a web-based graphical user interface (GUI), implemented with state-of-the-art technology, which incorporates system-specific analysis designed by the trajectory providers. A REST API and associated Jupyter Notebooks are integrated into the platform, allowing fully customized meta-analysis by final users. The new technology is illustrated using a collection of trajectories obtained by the community in the context of the effort to fight the COVID-19 pandemic. The server is accessible at https://bioexcel-cv19.bsc.es/#/. It is free and open to all users and there are no login requirements. It is also integrated into the simulations section of the BioExcel-MolSSI COVID-19 Molecular Structure and Therapeutics Hub: https://covid.molssi.org/simulations/ and is part of the MDDB effort (https://mddbr.eu).

Immunology of SARS-CoV-2 Infection

According to the data from the World Health Organization, about 800 million of the world population had contracted coronavirus infection caused by SARS-CoV-2 by mid-2023. Properties of this virus have allowed it to circulate in the human population for a long time, evolving defense mechanisms against the host immune system. Severity of the disease depends largely on the degree of activation of the systemic immune response, including overstimulation of macrophages and monocytes, cytokine production, and triggering of adaptive T- and B-cell responses, while SARS-CoV-2 evades the immune system actions. In this review, we discuss immune responses triggered in response to the SARS-CoV-2 virus entry into the cell and malfunctions of the immune system that lead to the development of severe disease.

SARS-CoV-2S-Protein-Ace2 Binding Analysis Using Surface Plasmon Resonance

Surface plasmon resonance (SPR) allows for the label-free determination of the binding affinity and rate constants of bimolecular interactions. Here, we describe the method used for the analysis of the Ace2-SARS-CoV2 S-protein interaction using indirect capture of the S-protein onto the SPR surface, and flowing monomeric Ace2. This method will allow for the determination of the rate constants for affinity, with additional analysis that is achievable using S-protein capture levels in conjunction with the sensorgram response for relative activity benchmarking.

Structural basis of G protein-Coupled receptor CMKLR1 activation and signaling induced by a chemerin-derived agonist

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

Dynamics of water-mediated interaction effects on the stability and transmission of Omicron

SARS-Cov-2 Omicron variant and its highly transmissible sublineages amidst news of emerging hybrid variants strengthen the evidence of its ability to rapidly spread and evolve giving rise to unprecedented future waves. Owing to the presence of isolated RBD, monomeric and trimeric Cryo-EM structures of spike protein in complex with ACE2 receptor, comparative analysis of Alpha, Beta, Gamma, Delta, and Omicron assist in a rational assessment of their probability to evolve as new or hybrid variants in future. This study proposes the role of hydration forces in mediating Omicron function and dynamics based on a stronger interplay between protein and solvent with each Covid wave. Mutations of multiple hydrophobic residues into hydrophilic residues underwent concerted interactions with water leading to variations in charge distribution in Delta and Omicron during molecular dynamics simulations. Moreover, comparative analysis of interacting moieties characterized a large number of mutations lying at RBD into constrained, homologous and low-affinity groups referred to as mutational drivers inferring that the probability of future mutations relies on their function. Furthermore, the computational findings reveal a significant difference in angular distances among variants of concern due 3 amino acid insertion (EPE) in Omicron variant that not only facilitates tight domain organization but also seems requisite for characterization of mutational processes. The outcome of this work signifies the possible relation between hydration forces, their impact on conformation and binding affinities, and viral fitness that will significantly aid in understanding dynamics of drug targets for Covid-19 countermeasures. The emerging scenario is that hydration forces and hydrophobic interactions are crucial variables to probe in mutational analysis to explore conformational landscape of macromolecules and reveal the molecular origins of protein behaviors.

Chronicling the 3-year evolution of the COVID-19 pandemic: analysis of disease management, characteristics of major variants, and impacts on pathogenicity

Announced on December 31, 2019, the novel coronavirus arising in Wuhan City, Hubei Province resulted in millions of cases and lives lost. Following intense tracking, coronavirus disease 2019 (COVID-19) was declared a pandemic by the World Health Organization (WHO) in 2020. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified as the cause of COVID-19 and the continuous evolution of the virus has given rise to several variants. In this review, a comprehensive analysis of the response to the pandemic over the first three-year period is provided, focusing on disease management, development of vaccines and therapeutics, and identification of variants. The transmissibility and pathogenicity of SARS-CoV-2 variants including Alpha, Beta, Gamma, Delta, and Omicron are compared. The binding characteristics of the SARS-CoV-2 spike protein to the angiotensin-converting enzyme 2 (ACE2) receptor and reproduction numbers are evaluated. The effects of major variants on disease severity, hospitalisation, and case-fatality rates are outlined. In addition to the spike protein, open reading frames mutations are investigated. We also compare the pathogenicity of SARS-CoV-2 with SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV). Overall, this study highlights the strengths and weaknesses of the global response to the pandemic, as well as the importance of prevention and preparedness. Monitoring the evolution of SARS-CoV-2 is critical in identifying and potentially predicting the health outcomes of concerning variants as they emerge. The ultimate goal would be a position in which existing vaccines and therapeutics could be adapted to suit new variants in as close to real-time as possible.

Comparative Study of the Mutations Observed in the SARS-CoV-2 RBD Variants of Concern and Their Impact on the Interaction with the ACE2 Protein

SARS-CoV-2 strains have made an appearance across the globe, causing over 757 million cases and over 6.85 million deaths at the time of writing. The emergence of these variants shows the amplitude of genetic variation to which the wild-type strains have been subjected. The rise of the different SARS-CoV-2 variants resulting from such genetic modification has significantly affected COVD-19's major impact on proliferation, virulence, and clinics. With the emergence of the variants of concern, the spike protein has been identified as a possible therapeutic target due to its critical role in binding to human cells and pathogenesis. These mutations could be linked to functional heterogeneity and use a different infection strategy. For example, the Omicron variant's multiple mutations should be carefully examined, as they represent one of the most widely spread strains and hint to us that there may be more genetic changes in the virus. As a result, we applied a common protocol where we reconstructed SARS-CoV-2 variants of concern and performed molecular dynamics simulations to study the stability of the ACE2-RBD complex in each variant. We also carried out free energy calculations to compare the binding and biophysical properties of the different SARS-CoV-2 variants when they interact with ACE2. Therefore, we were able to obtain consistent results and uncover new crucial residues that were essential for preserving a balance between maintaining a high affinity for ACE2 and the capacity to evade RBD-targeted antibodies. Our detailed structural analysis showed that SARS-CoV-2 variants of concern show a higher affinity for ACE2 compared to the Wuhan strain. Additionally, residues K417N and E484K/A might play a crucial role in antibody evasion, whereas Q498R and N501Y are specifically mutated to strengthen RBD affinity to ACE2 and, thereby, increase the viral effect of the COVID-19 virus.

Surfaces: a software to quantify and visualize interactions within and between proteins and ligands

SUMMARY

Computational methods for the quantification and visualization of the relative contribution of molecular interactions to the stability of biomolecular structures and complexes are fundamental to understand, modulate and engineer biological processes. Here, we present Surfaces, an easy to use, fast and customizable software for quantification and visualization of molecular interactions based on the calculation of surface areas in contact. Surfaces calculations shows equivalent or better correlations with experimental data as computationally expensive methods based on molecular dynamics.

AVAILABILITY AND IMPLEMENTATION

All scripts are available at https://github.com/NRGLab/Surfaces. Surface's documentation is available at https://surfaces-tutorial.readthedocs.io/en/latest/index.html.

Insights into Omicron's Low Fusogenicity through In Silico Molecular Studies on Spike-Furin Interactions

The Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant and its subvariants have a unique set of mutations. Two of those mutations (N679 K and P681 H) reside close to the S1 /S2 furin cleavage site (FCS; 685-686). When these mutations reside together, they exert less-efficient membrane fusion than wild type and most other variants of concern such as the Delta variant. Here, we in silico targeted these mutations to find out which of the amino acids and interactions change plays the key role in fusion. To comprehend the epistatic effect of N679 K and P681 H mutations on the spike protein, we in silico constructed three types of spike protein sequences by changing the respective amino acids on 679 and 681 positions (P681 H, N679 K, K679 N-H681 P variants). We then analyzed the binding affinity of furin and spike (Furin-Wild, Furin-Omicron, Furin-P681 H, Furin-N679 K, and Furin-K679 N/H681 P) complexes. Omicron and P681 H variants showed a similar higher binding energy trend compared to the wild type and N679 K. The variation in hydrogen, hydrophobic, and salt bridge bonds between spike protein and furin provided an explanation for the observed low fusogenicity of Omicron. The fate of the epistasis in furin binding and possible cleavage depends on the efficient interaction between FCS in spike and furin catalytic triad, and in addition, the loss of the hydrogen bond between Arg 681 (spike) and Asn 295 (furin) along with inhibitor-like ineffective higher affinity plays an important role in the enzymatic activity.

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

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

Structural basis of CMKLR1 signaling induced by chemerin9

Chemokine-like receptor 1 (CMKLR1), also known as chemerin receptor 23 (ChemR23) or chemerin receptor 1, is a chemoattractant G protein-coupled receptor (GPCR) that responds to the adipokine chemerin and is highly expressed in innate immune cells, including macrophages and neutrophils. The signaling pathways of CMKLR1 can lead to both pro- and anti-inflammatory effects depending on the ligands and physiological contexts. To understand the molecular mechanisms of CMKLR1 signaling, we determined a high-resolution cryo-electron microscopy (cryo-EM) structure of the CMKLR1-Gi signaling complex with chemerin9, a nanopeptide agonist derived from chemerin, which induced complex phenotypic changes of macrophages in our assays. The cryo-EM structure, together with molecular dynamics simulations and mutagenesis studies, revealed the molecular basis of CMKLR1 signaling by elucidating the interactions at the ligand-binding pocket and the agonist-induced conformational changes. Our results are expected to facilitate the development of small molecule CMKLR1 agonists that mimic the action of chemerin9 to promote the resolution of inflammation.

Understanding Mutations in Human SARS-CoV-2 Spike Glycoprotein: A Systematic Review & Meta-Analysis

Genetic variant(s) of concern (VoC) of SARS-CoV-2 have been emerging worldwide due to mutations in the gene encoding spike glycoprotein. We performed comprehensive analyses of spike protein mutations in the significant variant clade of SARS-CoV-2, using the data available on the Nextstrain server. We selected various mutations, namely, A222V, N439K, N501Y, L452R, Y453F, E484K, K417N, T478K, L981F, L212I, N856K, T547K, G496S, and Y369C for this study. These mutations were chosen based on their global entropic score, emergence, spread, transmission, and their location in the spike receptor binding domain (RBD). The relative abundance of these mutations was mapped with global mutation D614G as a reference. Our analyses suggest the rapid emergence of newer global mutations alongside D614G, as reported during the recent waves of COVID-19 in various parts of the world. These mutations could be instrumentally imperative for the transmission, infectivity, virulence, and host immune system's evasion of SARS-CoV-2. The probable impact of these mutations on vaccine effectiveness, antigenic diversity, antibody interactions, protein stability, RBD flexibility, and accessibility to human cell receptor ACE2 was studied in silico. Overall, the present study can help researchers to design the next generation of vaccines and biotherapeutics to combat COVID-19 infection.

Structural evolution of Delta lineage of SARS-CoV-2

One of the main obstacles in prevention and treatment of COVID-19 is the rapid evolution of the SARS-CoV-2 Spike protein. Given that Spike is the main target of common treatments of COVID-19, mutations occurring at this virulent factor can affect the effectiveness of treatments. The B.1.617.2 lineage of SARS-CoV-2, being characterized by many Spike mutations inside and outside of its receptor-binding domain (RBD), shows high infectivity and relative resistance to existing cures. Here, utilizing a wide range of computational biology approaches, such as immunoinformatics, molecular dynamics (MD), analysis of intrinsically disordered regions (IDRs), protein-protein interaction analyses, residue scanning, and free energy calculations, we examine the structural and biological attributes of the B.1.617.2 Spike protein. Furthermore, the antibody design protocol of Rosetta was implemented for evaluation the stability and affinity improvement of the Bamlanivimab (LY-CoV55) antibody, which is not capable of interactions with the B.1.617.2 Spike. We observed that the detected mutations in the Spike of the B1.617.2 variant of concern can cause extensive structural changes compatible with the described variation in immunogenicity, secondary and tertiary structure, oligomerization potency, Furin cleavability, and drug targetability. Compared to the Spike of Wuhan lineage, the B.1.617.2 Spike is more stable and binds to the Angiotensin-converting enzyme 2 (ACE2) with higher affinity.

CoVEffect: interactive system for mining the effects of SARS-CoV-2 mutations and variants based on deep learning

BACKGROUND

Literature about SARS-CoV-2 widely discusses the effects of variations that have spread in the past 3 years. Such information is dispersed in the texts of several research articles, hindering the possibility of practically integrating it with related datasets (e.g., millions of SARS-CoV-2 sequences available to the community). We aim to fill this gap, by mining literature abstracts to extract-for each variant/mutation-its related effects (in epidemiological, immunological, clinical, or viral kinetics terms) with labeled higher/lower levels in relation to the nonmutated virus.

RESULTS

The proposed framework comprises (i) the provisioning of abstracts from a COVID-19-related big data corpus (CORD-19) and (ii) the identification of mutation/variant effects in abstracts using a GPT2-based prediction model. The above techniques enable the prediction of mutations/variants with their effects and levels in 2 distinct scenarios: (i) the batch annotation of the most relevant CORD-19 abstracts and (ii) the on-demand annotation of any user-selected CORD-19 abstract through the CoVEffect web application (http://gmql.eu/coveffect), which assists expert users with semiautomated data labeling. On the interface, users can inspect the predictions and correct them; user inputs can then extend the training dataset used by the prediction model. Our prototype model was trained through a carefully designed process, using a minimal and highly diversified pool of samples.

CONCLUSIONS

The CoVEffect interface serves for the assisted annotation of abstracts, allowing the download of curated datasets for further use in data integration or analysis pipelines. The overall framework can be adapted to resolve similar unstructured-to-structured text translation tasks, which are typical of biomedical domains.

How an ACE2 mimicking epitope-MIP nanofilm recognizes template-related peptides and the receptor binding domain of SARS-CoV-2

Here we aim to gain a mechanistic understanding of the formation of epitope-imprinted polymer nanofilms using a non-terminal peptide sequence, i.e. the peptide GFNCYFP (G485 to P491) of the SARS-CoV-2 receptor binding domain (RBD). This epitope is chemisorbed on the gold surface through the central cysteine 488 followed by the electrosynthesis of a ∼5 nm thick polyscopoletin film around the surface confined templates. The interaction of peptides and the parent RBD and spike protein with the imprinted polyscopoletin nanofilm was followed by electrochemical redox marker gating, surface enhanced infrared absorption spectroscopy and conductive AFM. Because the use of non-terminal epitopes is especially intricate, here we characterize the binding pockets through their interaction with 5 peptides rationally derived from the template sequence, i.e. implementing central single amino acid mismatch as well as elongations and truncations at its C- and N- termini. Already a single amino acid mismatch, i.e. the central Cys488 substituted by a serine, results in ca. 15-fold lower affinity. Further truncation of the peptides to tetrapeptide (EGFN) and hexapeptide (YFPLQS) results also in a significantly lower affinity. We concluded that the affinity towards the different peptides is mainly determined by the four amino acid motif CYFP present in the sequence of the template peptide. A higher affinity than that for the peptides is found for the parent proteins RBD and spike protein, which seems to be due to out of cavity effects caused by their larger footprint on the nanofilm surface.

SARS-CoV-2 Delta Variant: Interplay between Individual Mutations and Their Allosteric Synergy

Since its first appearance in April 2021, B.1.617.2, also termed variant Delta, catalyzed one major worldwide wave dominating the second year of coronavirus disease 2019 (COVID-19) pandemic. Despite its quick disappearance worldwide, the strong virulence caused by a few point mutations remains an unsolved problem largely. Along with the other two sublineages, the Delta variant harbors an accumulation of Spike protein mutations, including the previously identified L452R, E484Q, and the newly emerged T478K on its receptor binding domain (RBD). We used molecular dynamics (MD) simulations, in combination with free energy perturbation (FEP) calculations, to examine the effects of two combinative mutation sets, L452R + E484Q and L452R + T478K. Our dynamic trajectories reveal an enhancement in binding affinity between mutated RBD and the common receptor protein angiotensin converting enzyme 2 (ACE2) through a net increase in the buried molecular surface area of the binary complex. This enhanced binding, mediated through Gln493, sets the same stage for all three sublineages due to the presence of L452R mutation. The other mutation component, E484Q or T478K, was found to impact the RBD-ACE2 binding and help the variant to evade several monoclonal antibodies (mAbs) in a distinct manner. Especially for L452R + T478K, synergies between mutations are mediated through a complex residual and water interaction network and further enhance its binding to ACE2. Taking together, this study demonstrates that new variants of SARS-CoV-2 accomplish both "attack" (infection) and "defense" (antibody neutralization escape) with the same "polished sword" (mutated Spike RBD).

Immunological study of COVID-19 vaccine candidate based on recombinant spike trimer protein from different SARS-CoV-2 variants of concern

The emergency of new SARS-CoV-2 variants that feature increased immune escape marks an urgent demand for better vaccines that will provide broader immunogenicity. Here, we evaluated the immunogenic capacity of vaccine candidates based on the recombinant trimeric spike protein (S) of different SARS-CoV-2 variants of concern (VOC), including the ancestral Wuhan, Beta and Delta viruses. In particular, we assessed formulations containing either single or combined S protein variants. Our study shows that the formulation containing the single S protein from the ancestral Wuhan virus at a concentration of 2µg (SW2-Vac 2µg) displayed in the mouse model the highest IgG antibody levels against all the three (Wuhan, Beta, and Delta) SARS-CoV-2 S protein variants tested. In addition, this formulation induced significantly higher neutralizing antibody titers against the three viral variants when compared with authorized Gam-COVID-Vac-rAd26/rAd5 (Sputnik V) or ChAdOx1 (AstraZeneca) vaccines. SW2-Vac 2µg was also able to induce IFN-gamma and IL-17, memory CD4 populations and follicular T cells. Used as a booster dose for schedules performed with different authorized vaccines, SW2-Vac 2µg vaccine candidate also induced higher levels of total IgG and IgG isotypes against S protein from different SARS-CoV-2 variants in comparison with those observed with homologous 3-dose schedule of Sputnik V or AstraZeneca. Moreover, SW2-Vac 2µg booster induced broadly strong neutralizing antibody levels against the three tested SARS-CoV-2 variants. SW2-Vac 2µg booster also induced CD4+ central memory, CD4+ effector and CD8+ populations. Overall, the results demonstrate that SW2-Vac 2 µg is a promising formulation for the development of a next generation COVID-19 vaccine.

Surface charge changes in spike RBD mutations of SARS-CoV-2 and its variant strains alter the virus evasiveness via HSPGs: A review and mechanistic hypothesis

With the COVID-19 pandemic continuing, more contagious SARS-CoV-2 variants, including Omicron, have been emerging. The mutations, especially those that occurred on the spike (S) protein receptor-binding domain (RBD), are of significant concern due to their potential capacity to increase viral infectivity, virulence, and breakthrough antibodies' protection. However, the molecular mechanism involved in the pathophysiological change of SARS-CoV-2 mutations remains poorly understood. Here, we summarized 21 RBD mutations and their human angiotensin-converting enzyme 2 (hACE2) and/or neutralizing antibodies' binding characteristics. We found that most RBD mutations, which could increase surface positive charge or polarity, enhanced their hACE2 binding affinity and immune evasion. Based on the dependence of electrostatic interaction of the epitope residue of virus and docking protein (like virus receptors or antibodies) for its invasion, we postulated that the charge and/or polarity changes of novel mutations on the RBD domain of S protein could affect its affinity for the hACE2 and antibodies. Thus, we modeled mutant S trimers and RBD-hACE2 complexes and calculated their electrotactic distribution to study surface charge changes. Meanwhile, we emphasized that heparan sulfate proteoglycans (HSPGs) might play an important role in the hACE2-mediated entry of SARS-CoV-2 into cells. Those hypotheses provide some hints on how SARS-CoV-2 mutations enhance viral fitness and immune evasion, which may indicate potential ways for drug design, next-generation vaccine development, and antibody therapies.

Role of the Pangolin in Origin of SARS-CoV-2: An Evolutionary Perspective

After the recent emergence of SARS-CoV-2 infection, unanswered questions remain related to its evolutionary history, path of transmission or divergence and role of recombination. There is emerging evidence on amino acid substitutions occurring in key residues of the receptor-binding domain of the spike glycoprotein in coronavirus isolates from bat and pangolins. In this article, we summarize our current knowledge on the origin of SARS-CoV-2. We also analyze the host ACE2-interacting residues of the receptor-binding domain of spike glycoprotein in SARS-CoV-2 isolates from bats, and compare it to pangolin SARS-CoV-2 isolates collected from Guangdong province (GD Pangolin-CoV) and Guangxi autonomous regions (GX Pangolin-CoV) of South China. Based on our comparative analysis, we support the view that the Guangdong Pangolins are the intermediate hosts that adapted the SARS-CoV-2 and represented a significant evolutionary link in the path of transmission of SARS-CoV-2 virus. We also discuss the role of intermediate hosts in the origin of Omicron.

Impact of the temperature on the interactions between common variants of the SARS-CoV-2 receptor binding domain and the human ACE2

Several key mutations in the Spike protein receptor binding domain (RBD) have been identified to influence its affinity for the human Angiotensin-Converting Enzyme 2 (ACE2). Here, we perform a comparative study of the ACE2 binding to the wild type (Wuhan) RBD and some of its variants: Alpha B.1.1.7, Beta B.1.351, Delta B.1.617.2, Kappa B.1.617.1, B.1.1.7 + L452R and Omicron B.1.1.529. Using a coiled-coil mediated tethering approach of ACE2 in a novel surface plasmon resonance (SPR)-based assay, we measured interactions at different temperatures. Binding experiments at 10 °C enhanced the kinetic dissimilarities between the RBD variants and allowed a proper fit to a Langmuir 1:1 model with high accuracy and reproducibility, thus unraveling subtle differences within RBD mutants and ACE2 glycovariants. Our study emphasizes the importance of SPR-based assay parameters in the acquisition of biologically relevant data and offers a powerful tool to deepen our understanding of the role of the various RBD mutations in ACE2 interaction binding parameters.

Molecular characteristics, immune evasion, and impact of SARS-CoV-2 variants

The persistent COVID-19 pandemic since 2020 has brought an enormous public health burden to the global society and is accompanied by various evolution of the virus genome. The consistently emerging SARS-CoV-2 variants harboring critical mutations impact the molecular characteristics of viral proteins and display heterogeneous behaviors in immune evasion, transmissibility, and the clinical manifestation during infection, which differ each strain and endow them with distinguished features during populational spread. Several SARS-CoV-2 variants, identified as Variants of Concern (VOC) by the World Health Organization, challenged global efforts on COVID-19 control due to the rapid worldwide spread and enhanced immune evasion from current antibodies and vaccines. Moreover, the recent Omicron variant even exacerbated the global anxiety in the continuous pandemic. Its significant evasion from current medical treatment and disease control even highlights the necessity of combinatory investigation of the mutational pattern and influence of the mutations on viral dynamics against populational immunity, which would greatly facilitate drug and vaccine development and benefit the global public health policymaking. Hence in this review, we summarized the molecular characteristics, immune evasion, and impacts of the SARS-CoV-2 variants and focused on the parallel comparison of different variants in mutational profile, transmissibility and tropism alteration, treatment effectiveness, and clinical manifestations, in order to provide a comprehensive landscape for SARS-CoV-2 variant research.

Point mutations in SARS-CoV-2 variants induce long-range dynamical perturbations in neutralizing antibodies

Monoclonal antibodies are emerging as a viable treatment for the coronavirus disease 19 (COVID-19). However, newly evolved variants of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) can reduce the efficacy of currently available antibodies and can diminish vaccine-induced immunity. Here, we demonstrate that the microscopic dynamics of neutralizing monoclonal antibodies can be profoundly modified by the mutations present in the spike proteins of the SARS-COV-2 variants currently circulating in the world population. The dynamical perturbations within the antibody structure, which alter the thermodynamics of antigen recognition, are diverse and can depend both on the nature of the antibody and on the spatial location of the spike mutation. The correlation between the motion of the antibody and that of the spike receptor binding domain (RBD) can also be changed, modulating binding affinity. Using protein-graph-connectivity networks, we delineated the mutant-induced modifications in the information-flow along allosteric pathway throughout the antibody. Changes in the collective dynamics were spatially distributed both locally and across long-range distances within the antibody. On the receptor side, we identified an anchor-like structural element that prevents the detachment of the antibodies; individual mutations there can significantly affect the antibody binding propensity. Our study provides insight into how virus neutralization by monoclonal antibodies can be impacted by local mutations in the epitope via a change in dynamics. This realization adds a new layer of sophistication to the efforts for rational design of monoclonal antibodies against new variants of SARS-CoV2, taking the allostery in the antibody into consideration.

Mutations in the new variants of SARS-CoV-2 spike protein modulates the dynamics of the neutralizing antibodies. Capturing such modulations from MD simulations and graph network model identifies the role of mutations in facilitating immune evasion.