Insights into the structure and dynamics of SARS-CoV-2 spike glycoprotein double mutant L452R-E484Q

Impact of mutations in SARS-CoV-2 recombinant sub-variant XBB.1.16 on the binding affinity with human ACE2 receptor

Despite the waning threat of the COVID-19 pandemic, its detrimental impact on global health persists. Regardless of natural immunity or immunity obtained through vaccination, emerging variants of the virus continue to undergo mutations and propagate globally. The persistent mutations in SARS-CoV-2, along with the subsequent formation of recombinant sub-variants has become a challenge for researchers and health professionals, raising concerns about the efficacy of current vaccines. Gaining a better understanding of the biochemical interactions between the Spike Protein (RBD) of SARS-CoV-2 variants and the human ACE2 receptor can prove to be beneficial in designing and developing antiviral therapeutics that are equally effective against all strains and emerging variants. Our objective in this study was to investigate the interfacial binding pattern of the SARS-CoV-2 RBD-ACE2 complex of the Wild Type (WT), Omicron, and the Omicron recombinant sub-variant XBB.1.16. We aimed to examine the atomic level factors and observe how mutations influence the interaction between the virus and its host using Molecular Dynamics simulation, MM/GBSA energy calculations, and Principal Component Analysis. Our findings reveal a higher degree of structural deviation and flexibility in XBB.1.16 compared to WT and Omicron. PCA indicated a wider cluster and significant flexibility in the movements of XBB.1.16 which can also be observed in free energy landscapes, while the normal mode analysis revealed converging motions within the RBD-ACE2 complexes which can facilitate the interaction between them. A pattern of decreased binding affinity was observed in case of XBB.1.16 when compared to the WT and Omicron. These observed deviations in XBB.1.16 when compared to its parent lineage Omicron, and WT can be attributed to the mutations specific to it. Collectively, these results enhance our understanding of the impact of mutations on the interaction between this strain and the host, taking us one step closer to designing effective antiviral therapeutics against the continually mutating strains.

Computational insights into dynamics and conformational stability of N-acetylmannosamine kinase mutations

The activity of UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE) is essential for the biosynthesis of sialic acid, which is involved in cellular processes in health and diseases. GNE contains an N-terminal epimerase domain and a C-terminal kinase domain (N-acetylmannosamine kinase, MNK). Mutations of the GNE protein led to hypoactivity of the enzyme and cause sialurea or autosomal recessive inclusion body myopathy/Nonaka myopathy. Here, we used all-atom molecular dynamics (MD) simulations to comprehend the folding, dynamics and conformational stability of MNK variants, including the wild type (WT) and three mutants (H677R, V696M and H677R/V696M). The deleterious and destabilizing nature of MNK mutants were predicted using different prediction tools. Results predicted that mutations modulate the stability, flexibility and function of MNK. The effect of mutations on the conformational stability and dynamics of MNK was next studied through the free-energy landscape (FEL), hydrogen-bonds and secondary structure changes. The FEL results show that the mutations interfere with various conformational transitions in both WT and mutants, exposing the structural underpinnings of protein destabilization and unfolding brought on by mutation. We discover that, when compared to the other two mutations, V696M and H677R/V696M, H677R has the most harmful effects. These findings have a strong correlation with published experimental studies that demonstrate how these mutations disrupt MNK activity. Hence, this computational study describes the structural details to unravel the mutant effects at the atomistic resolution and has implications for understanding the GNE's physiological and pathological role.Communicated by Ramaswamy H. Sarma.

Comprehensive deep mutational scanning reveals the pH induced stability and binding differences between SARS-CoV-2 spike RBD and human ACE2

The SARS-CoV-2 spike (S) glycoprotein with its mobile receptor-binding domain (RBD), binds to the human ACE2 receptor and thus facilitates virus entry through low-pH-endosomal pathways. The high degree of SARS-CoV-2 mutability has raised concern among scientists and medical professionals because it created doubt about the effectiveness of drugs and vaccinations designed specifically for COVID-19. In this study, we used computational saturation mutagenesis approach, including structure-based free energy calculations to analyse the effects of the missense mutations on the SARS-CoV-2 S-RBD stability and the S-RBD binding affinity with ACE2 at three different pH (pH 4.5, pH 6.5, and pH 7.4). A total of 3705 mutations in the S-RBD protein were analyzed, and we discovered that most of these mutations destabilize the RBD protein. Specifically, residues G404, G431, G447, A475, and G526 were important for RBD protein stability. In addition, RBD residues Y449, Y489, Y495, Q498, and N487 were critical for the RBD-ACE2 interaction. Next, we found that the distribution of the mean stability changes and mean binding energy changes of RBD due to mutations at both serological and endosomal pH correlated well, indicating the similar effects of mutations. Overall, this computational analysis is useful for understanding the effects of missense mutations in SARS-CoV-2 pathogenesis at different pH.Communicated by Ramaswamy H. Sarma.

In silico identification of microRNAs targeting the PPARα/γ: promising therapeutics for SARS-CoV‑2 infection

The SARS-CoV-2 lifecycle is dependent on the host metabolism machinery. It upregulates the PPARα and PPARγ genes in lipid metabolism, which supports the essential viral replication complex including lipid rafts and palmitoylation of viral protein. The use of PPAR ligands in SARS-CoV-2 infection may have positive effects by preventing cytokine storm and the ensuing inflammatory cascade. The inhibition of PPARα and PPARγ genes may alter the metabolism and may disrupt the lifecycle of SARS-CoV-2 and COVID-19 progression. In the present work, we have identified possible miRNAs targeting PPARα and PPARγ in search of modulators of PPARα and PPARγ genes expression. The identified miRNAs could possibly be viewed as new therapeutic targets against COVID-19 infection.

The rapid emergence of multiple sublineages of Omicron (B.1.1.529) variant: Dynamic profiling via molecular phylogenetics and mutational landscape studies

PURPOSE

The recent Omicron (B.1.1.529) variant poses a significant threat to global health. This variant has spread worldwide, and several sublineages have rapidly emerged. Study tried to analyze the microevolution of this variant.

METHODS

We studied the molecular phylogenetics, divergence, geographical distributions, frequencies, risk mutations for antibody affinity, and mutational landscape for Omicron sublineages using in silico analysis and statistical models. The risk mutation of spike for nAb affinity was analyzed and illustrated by statistical plots. Finally, the mutational properties of the spike mutations and their stability were predicted and demonstrated.

RESULTS

First, we studied the microevolutionary Omicron sublineages using molecular phylogenetics. Simultaneously, we revealed divergence events of the Omicron sublineages and observed the lowest minimum divergence of 51 in clade 21K and the highest maximum divergence of 90 in clade 21L. We have demonstrated cluster analyses, geographical distributions, frequencies of Omicron and its sublineages. Finally, we evaluated the mutational landscape of the Omicron sublineages. In this mutational study, we performed a genome-wide analysis of general mutations, mutations in the non-spike genome, and spike mutations of Omicron sublineages. The risk mutation of S-glycoprotein for nAb affinity has been analyzed for Omicron sublineages. Here, we found that some sublineages have all four significant highly destabilizing mutations. Such sublineages are BA.1 (G446S, E484A, T95I, and D614G), BA.2 (H655Y, Q493R, G493S, and D614G), BA.4 (N501Y, Y505H, N969K, and D614G), and BA.2.75 (Q454H, T547K, N764K, D614G and G446S). Finally, from the mutation stability prediction through ΔΔG, we observed that BA.1 and BA.4 had two destabilizing and two stabilizing mutations. Similarly, BA.2, BA.5, and BA.2.12.1 have one destabilizing and three stabilizing mutations. However, all four mutations in BA.2.75 are stabilizing mutations.

CONCLUSIONS

Our molecular phylogenetic studies provided a deeper understanding of the microevolution of sublineages and the creation of Omicron. Similarly, this study might help scientists develop pan-coronavirus vaccines that consider their mutational properties.

Isolation and In Silico Inhibitory Potential against SARS-CoV-2 RNA Polymerase of the Rare Kaempferol 3-O-(6″-O-acetyl)-Glucoside from Calligonum tetrapterum

The phytochemical constituents of Calligonum tetrapterum Jaub. & Spach (Family Polygonaceae) were studied for the first time. The study resulted in the isolation of the rare flavonol glycoside, kaempferol 3-O-(6″-O-acetyl)-glucoside,(K3G-A). The potential inhibitive activity of K3G-A toward SARS-CoV-2 was investigated utilizing several in silico approaches. First, molecular fingerprints and structural similarity experiments were carried out for K3G-A against nine co-crystallized ligands of nine proteins of SARS-CoV-2 to reveal if there is a structural similarity with any of them. The conducted studies showed the high similarity of K3G-A and remdesivir, the co-crystallized ligand of SARS-CoV-2 RNA-dependent RNA polymerase (PDB ID: 7BV2), RdRp. To validate these findings, a DFT study was conducted and confirmed the proposed similarity on the electronic and orbital levels. The binding of K3G-A against RdRp was confirmed through molecular docking studies exhibiting a binding energy of −27.43 kcal/mol, which was higher than that of remdesivir. Moreover, the RdRp-K3G-A complex was subjected to several MD studies at 100 ns that authenticated the accurate mode of binding and the correct dynamic behavior. Finally, in silico ADMET and toxicity evaluation of K3G-A was conducted and denoted the safety and the drug-likeness of K3G-A. In addition to K3G-A, two other metabolites were isolated and identified to be kaempferol (K) and β-sitosterol (β-S).