Mechanistic insights into the effects of key mutations on SARS-CoV-2 RBD-ACE2 binding

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.

Most accurate mutations in SARS-CoV-2 genomes identified in Uzbek patients show novel amino acid changes

Purpose

The rapid changes in the coronavirus genomes created new strains after the first variation was found in Wuhan in 2019. SARS-CoV-2 genotypes should periodically undergo whole genome sequencing to control it because it has been extremely helpful in combating the virus. Many diagnoses, treatments, and vaccinations have been developed against it based on genome sequencing. With its practical implications, this study aimed to determine changes in the delta variant of SARS-CoV-2 widespread in Uzbekistan during the pandemic by genome sequencing, thereby providing crucial insights for developing effective control strategies that can be directly applied in the field.

Design

We meticulously generated 17 high-quality whole-genome sequence data from 48 SARS-CoV-2 genotypes of COVID-19 patients who tested positive by PCR in Tashkent, Uzbekistan. Our rigorous approach, which includes stringent quality control measures and multiple rounds of verification, ensures the accuracy and reliability of our findings.

Methods

Our study employed a unique combination of genome sequencing and bioinformatics web tools to analyze amino acid (AA) changes in the virus genomes. This approach allowed us to understand the genetic changes in the delta variant of SARS-CoV-2 widespread in Uzbekistan during the pandemic.

Results

Our study revealed significant nucleotide polymorphisms, including non-synonymous (missense) and synonymous mutations in the coding regions of the sequenced sample genomes. These findings, categorized by phylogenetic analysis into the G clade (or GK sub-clade), contribute to our understanding of the delta variant of SARS-CoV-2 widespread in Uzbekistan during the pandemic. A total of 134 mutations were identified, consisting of 65 shared and 69 unique mutations. These nucleotide changes, including one frameshift mutation, one conservative and disruptive insertion-deletion, four upstream region mutations, four downstream region mutations, 39 synonymous mutations, and 84 missense mutations, are crucial in the ongoing battle against the virus.

Conclusion

The comprehensive whole-genome sequencing data presented in this study aids in tracing the origins and sources of circulating SARS-CoV-2 variants and analyzing emerging variations within Uzbekistan and globally. The genome sequencing of SARS-CoV-2 from samples collected in Uzbekistan in late 2021, during the peak of the pandemic's second wave nationwide, is detailed here. Following acquiring these sequences, research efforts have focused on developing DNA and plant-based edible vaccines utilizing prevalent SARS-CoV-2 strains in Uzbekistan, which are currently undergoing clinical trials.

Evaluating mAbs binding abilities to Omicron subvariant RBDs: implications for selecting effective mAb therapies

The ongoing evolution of the Omicron lineage of SARS-CoV-2 has led to the emergence of subvariants that pose challenges to antibody neutralization. Understanding the binding dynamics between the receptor-binding domains (RBD) of these subvariants spike and monoclonal antibodies (mAbs) is pivotal for elucidating the mechanisms of immune escape and for advancing the development of therapeutic antibodies. This study focused on the RBD regions of Omicron subvariants BA.2, BA.5, BF.7, and XBB.1.5, employing molecular dynamics simulations to unravel their binding mechanisms with a panel of six mAbs, and subsequently analyzing the origins of immune escape from energetic and structural perspectives. Our results indicated that the antibody LY-COV1404 maintained binding affinities across all studied systems, suggesting the resilience of certain antibodies against variant-induced immune escape, as seen with the mAb 1D1-Fab. The newly identified mAb 002-S21F2 showed a similar efficacy profile to LY-COV1404, though with a slightly reduced binding to BF.7. In parallel, mAb REGN-10933 emerged as a potential therapeutic candidate against BF.7 and XBB.1.5, reflecting the importance of identifying variant-specific antibody interactions, akin to the binding optimization observed in BA.4/5 and XBB.1.5. And key residues that facilitate RBD-mAb binding were identified (T345, L441, K444, V445, and T500), alongside residues that hinder protein-protein interactions (D420, L455, K440, and S446). Particularly noteworthy was the inhibited binding of V445 and R509 with mAbs in the presence of mAb 002-S21F2, suggesting a mechanism for immune escape, especially through the reduction of V445 hydrophobicity. These findings enhance our comprehension of the binding interactions between mAbs and RBDs, contributing to the understanding of immune escape mechanisms. They also lay the groundwork for the design and optimization of antiviral drugs and have significant implications for the development of treatments against current and future coronaviruses.

Unraveling the Dynamics of SARS-CoV-2 Mutations: Insights from Surface Plasmon Resonance Biosensor Kinetics

Surface Plasmon Resonance (SPR) technology is known to be a powerful tool for studying biomolecular interactions because it offers real-time and label-free multiparameter analysis with high sensitivity. This article summarizes the results that have been obtained from the use of SPR technology in studying the dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations. This paper will begin by introducing the working principle of SPR and the kinetic parameters of the sensorgram, which include the association rate constant (ka), dissociation rate constant (kd), equilibrium association constant (KA), and equilibrium dissociation constant (KD). At the end of the paper, we will summarize the kinetic data on the interaction between angiotensin-converting enzyme 2 (ACE2) and SARS-CoV-2 obtained from the results of SPR signal analysis. ACE2 is a material that mediates virus entry. Therefore, understanding the kinetic changes between ACE2 and SARS-CoV-2 caused by the mutation will provide beneficial information for drug discovery, vaccine development, and other therapeutic purposes.

Investigation of protein-protein interactions and hotspot region on the NSP7-NSP8 binding site in NSP12 of SARS-CoV-2

Background: The RNA-dependent RNA polymerase (RdRp) complex, essential in viral transcription and replication, is a key target for antiviral therapeutics. The core unit of RdRp comprises the nonstructural protein NSP12, with NSP7 and two copies of NSP8 (NSP81 and NSP82) binding to NSP12 to enhance its affinity for viral RNA and polymerase activity. Notably, the interfaces between these subunits are highly conserved, simplifying the design of molecules that can disrupt their interaction. Methods: We conducted a detailed quantum biochemical analysis to characterize the interactions within the NSP12-NSP7, NSP12-NSP81, and NSP12-NSP82 dimers. Our objective was to ascertain the contribution of individual amino acids to these protein-protein interactions, pinpointing hotspot regions crucial for complex stability. Results: The analysis revealed that the NSP12-NSP81 complex possessed the highest total interaction energy (TIE), with 14 pairs of residues demonstrating significant energetic contributions. In contrast, the NSP12-NSP7 complex exhibited substantial interactions in 8 residue pairs, while the NSP12-NSP82 complex had only one pair showing notable interaction. The study highlighted the importance of hydrogen bonds and π-alkyl interactions in maintaining these complexes. Intriguingly, introducing the RNA sequence with Remdesivir into the complex resulted in negligible alterations in both interaction energy and geometric configuration. Conclusion: Our comprehensive analysis of the RdRp complex at the protein-protein interface provides invaluable insights into interaction dynamics and energetics. These findings can guide the design of small molecules or peptide/peptidomimetic ligands to disrupt these critical interactions, offering a strategic pathway for developing effective antiviral drugs.

Modelling how increased Cathepsin B/L and decreased TMPRSS2 usage for cell entry by the SARS-CoV-2 Omicron variant may affect the efficacy and synergy of TMPRSS2 and Cathepsin B/L inhibitors

The SARS-CoV-2 Omicron variant harbours many mutations in its spike protein compared to the original SARS-CoV-2 strain, which may alter its ability to enter cells, cell tropism, and response to interventions blocking virus entry. To elucidate these effects, we developed a mathematical model of SARS-CoV-2 entry into target cells and applied it to analyse recent in vitro data. SARS-CoV-2 can enter cells via two pathways, one using the host proteases Cathepsin B/L and the other using the host protease TMPRSS2. We found enhanced entry efficiency of the Omicron variant in cells where the original strain preferentially used Cathepsin B/L and reduced efficiency where it used TMPRSS2. The Omicron variant thus appears to have evolved to use the Cathepsin B/L pathway better but at the expense of its ability to use the TMPRSS2 pathway compared to the original strain. We estimated >4-fold enhanced efficiency of the Omicron variant in entry via the Cathepsin B/L pathway and >3-fold reduced efficiency via the TMPRSS2 pathway compared to the original or other strains in a cell type-dependent manner. Our model predicted that Cathepsin B/L inhibitors would be more efficacious and TMPRSS2 inhibitors less efficacious in blocking Omicron variant entry into cells than the original strain. Furthermore, model predictions suggested that drugs simultaneously targeting the two pathways would exhibit synergy. The maximum synergy and drug concentrations yielding it would differ for the Omicron variant compared to the original strain. Our findings provide insights into the cell entry mechanisms of the Omicron variant and have implications for intervention targeting these mechanisms.

Virological characteristics of the SARS-CoV-2 XBB variant derived from recombination of two Omicron subvariants

In late 2022, SARS-CoV-2 Omicron subvariants have become highly diversified, and XBB is spreading rapidly around the world. Our phylogenetic analyses suggested that XBB emerged through the recombination of two cocirculating BA.2 lineages, BJ.1 and BM.1.1.1 (a progeny of BA.2.75), during the summer of 2022. XBB.1 is the variant most profoundly resistant to BA.2/5 breakthrough infection sera to date and is more fusogenic than BA.2.75. The recombination breakpoint is located in the receptor-binding domain of spike, and each region of the recombinant spike confers immune evasion and increases fusogenicity. We further provide the structural basis for the interaction between XBB.1 spike and human ACE2. Finally, the intrinsic pathogenicity of XBB.1 in male hamsters is comparable to or even lower than that of BA.2.75. Our multiscale investigation provides evidence suggesting that XBB is the first observed SARS-CoV-2 variant to increase its fitness through recombination rather than substitutions.

Convergent evolution of SARS-CoV-2 Omicron subvariants leading to the emergence of BQ.1.1 variant

In late 2022, various Omicron subvariants emerged and cocirculated worldwide. These variants convergently acquired amino acid substitutions at critical residues in the spike protein, including residues R346, K444, L452, N460, and F486. Here, we characterize the convergent evolution of Omicron subvariants and the properties of one recent lineage of concern, BQ.1.1. Our phylogenetic analysis suggests that these five substitutions are recurrently acquired, particularly in younger Omicron lineages. Epidemic dynamics modelling suggests that the five substitutions increase viral fitness, and a large proportion of the fitness variation within Omicron lineages can be explained by these substitutions. Compared to BA.5, BQ.1.1 evades breakthrough BA.2 and BA.5 infection sera more efficiently, as demonstrated by neutralization assays. The pathogenicity of BQ.1.1 in hamsters is lower than that of BA.5. Our multiscale investigations illuminate the evolutionary rules governing the convergent evolution for known Omicron lineages as of 2022.

Interactions of angiotensin-converting enzyme-2 (ACE2) and SARS-CoV-2 spike receptor-binding domain (RBD): a structural perspective

BACKGROUND

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has caused millions of infections and deaths worldwide since its discovery in late 2019 in Wuhan, China. The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein binds to the human angiotensin-converting enzyme-2 (ACE2) receptor, a critical component of the renin-angiotensin system (RAS) that initiates the viral transmission. Most of the critical mutations found in SARS-CoV-2 are associated with the RBD of the spike protein. These mutations have the potential to reduce the efficacy of vaccines and neutralizing antibodies.

METHODS

In this review, the structural details of ACE2, RBD and their interactions are discussed. In addition, some critical mutations of RBD and their impact on ACE2-RBD interactions are also discussed.

CONCLUSION

Preventing the interaction between Spike RBD and ACE2 is considered a viable therapeutic strategy since ACE2 binding by RBD is the first step in virus infection. Because the interactions between the two entities are critical for both viral transmission and therapeutic development, it is essential to understand their interactions in detail.

Aloin A inhibits SARS CoV-2 replication by targeting its binding with ACE2 - Evidence from modeling-supported molecular dynamics simulation

The current study aimed to expand on the recently published results and assess the inhibitory efficacy of aloin A against SARS CoV-2. In vitro testing of aloin A against SARS CoV-2 proteases (i.e., MPro and PLPro) showed weak to moderate activity (IC50 = 68.56 ± 1.13 µM and 24.77 ± 1.57 µM, respectively). However, aloin A was able to inhibit the replication of SARS CoV-2 in Vero E6 cells efficiently with an IC50 of 0.095 ± 0.022 µM. Depending on the reported poor permeability of aloin A alongside its insignificant protease inhibitory activities presented in this study, we ran a number of extensive virtual screenings and physics-based simulations to determine the compound's potential mode of action. As a result, RBD-ACE2 was identified as a key target for aloin A. Results from 600 ns-long molecular dynamics (MD) simulation experiments pointed to aloin A's role as an RBD-ACE2 destabilizer. Therefore, the results of this work may pave the way for further development of this scaffold and the eventual production of innovative anti-SARS CoV-2 medicines with several mechanisms of action.Communicated by Ramaswamy H. Sarma.

Omicron BA.1 and BA.2 variants increase the interactions of SARS-CoV-2 spike glycoprotein with ACE2

SARS-CoV-2 infection is initiated by binding of the receptor-binding domain (RBD) of its spike glycoprotein to the peptidase domain (PD) of angiotensin-converting enzyme 2 (ACE2) receptors in host cells. Recently detected Omicron variant of SARS-CoV-2 (B.1.1.529) is heavily mutated on RBD. First the BA.1 and later the BA.2 variant became the most dominant strains of the Omicron variant. To investigate how the mutations of these strains affect RBD-PD interactions, we performed all-atom molecular dynamics simulations of the BA.1 and BA.2 RBD-PD in the presence of full-length glycans, explicit water, and ions. Simulations revealed that RBDs of BA.1 and BA.2 variants exhibit a more dispersed interaction network and make an increased number of salt bridges and hydrophobic interactions with PD compared to wild-type RBD. Although BA.1 and BA.2 differ in two residues at the RBD-ACE2 interface, no major difference in RBD-PD interactions and binding strengths were observed between these variants. Using the conformations sampled in each trajectory, the Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) method estimated ∼34% and ∼51% stronger binding free energies to PD for BA.1 and BA.2 RBD, respectively, than wild-type RBD, which may result in higher binding efficiency of the Omicron variant to infect host cells.

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).

Computational design of stapled peptide inhibitor against SARS‐CoV ‐2 receptor binding domain

Since its first detection in 2019, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS‐CoV‐2) has been the cause of millions of deaths worldwide. Despite the development and administration of different vaccines, the situation is still worrisome as the virus is constantly mutating to produce newer variants some of which are highly infectious. This raises an urgent requirement to understand the infection mechanism and thereby design therapeutic‐based treatment for COVID‐19. The gateway of the virus to the host cell is mediated by the binding of the receptor binding domain (RBD) of the virus spike protein to the angiotensin‐converting enzyme 2 (ACE2) of the human cell. Therefore, the RBD of SARS‐CoV‐2 can be used as a target to design therapeutics. The α1 helix of ACE2, which forms direct contact with the RBD surface, has been used as a template in the current study to design stapled peptide therapeutics. Using computer simulation, the mechanism and thermodynamics of the binding of six stapled peptides with RBD have been estimated. Among these, the one with two lactam stapling agents has shown binding affinity, sufficient to overcome RBD‐ACE2 binding. Analyses of the mechanistic detail reveal that a reorganization of amino acids at the RBD‐ACE2 interface produces favorable enthalpy of binding whereas conformational restriction of the free peptide reduces the loss in entropy to result higher binding affinity. The understanding of the relation of the nature of the stapling agent with their binding affinity opens up the avenue to explore stapled peptides as therapeutic against SARS‐CoV‐2.

Pharmaceutical Prospects of Curcuminoids for the Remedy of COVID-19: Truth or Myth

Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is a positive-strand RNA virus, and has rapidly spread worldwide as a pandemic. The vaccines, repurposed drugs, and specific treatments have led to a surge of novel therapies and guidelines nowadays; however, the epidemic of COVID-19 is not yet fully combated and is still in a vital crisis. In repositioning drugs, natural products are gaining attention because of the large therapeutic window and potent antiviral, immunomodulatory, anti-inflammatory, and antioxidant properties. Of note, the predominant curcumoid extracted from turmeric (Curcuma longa L.) including phenolic curcumin influences multiple signaling pathways and has demonstrated to possess anti-inflammatory, antioxidant, antimicrobial, hypoglycemic, wound healing, chemopreventive, chemosensitizing, and radiosensitizing spectrums. In this review, all pieces of current information related to curcumin-used for the treatment and prevention of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection through in vitro, in vivo, and in silico studies, clinical trials, and new formulation designs are retrieved to re-evaluate the applications based on the pharmaceutical efficacy of clinical therapy and to provide deep insights into knowledge and strategy about the curcumin's role as an immune booster, inflammatory modulator, and therapeutic agent against COVID-19. Moreover, this study will also afford a favorable application or approach with evidence based on the drug discovery and development, pharmacology, functional foods, and nutraceuticals for effectively fighting the COVID-19 pandemic.

Deciphering binding mechanism of inhibitors to SARS-COV-2 main protease through multiple replica accelerated molecular dynamics simulations and free energy landscapes

The pneumonia outbreak caused by the SARS-CoV-2 virus poses a serious threat to human health and the world economy. Development of safe and highly effective antiviral drugs is of great...