Analysis of SARS-CoV-2 mutations in the main viral protease (NSP5) and its implications on the vaccine designing strategies

Comparative Atlas of SARS-CoV-2 Substitution Mutations: A Focus on Iranian Strains Amidst Global Trends

BACKGROUND

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a new emerging coronavirus that caused coronavirus disease 2019 (COVID-19). Whole-genome tracking of SARS-CoV-2 enhanced our understanding of the mechanism of the disease, control, and prevention of COVID-19.

METHODS

we analyzed 3368 SARS-CoV-2 protein sequences from Iran and compared them with 15.6 million global sequences in the GISAID database, using the Wuhan-Hu-1 strain as a reference.

RESULTS

Our investigation revealed that NSP12-P323L, ORF9c-G50N, NSP14-I42V, membrane-A63T, Q19E, and NSP3-G489S were found to be the most frequent mutations among Iranian SARS-CoV-2 sequences. Furthermore, it was observed that more than 94% of the SARS-CoV-2 genome, including NSP7, NSP8, NSP9, NSP10, NSP11, and ORF8, had no mutations when compared to the Wuhan-Hu-1 strain. Finally, our data indicated that the ORF3a-T24I, NSP3-G489S, NSP5-P132H, NSP14-I42V, envelope-T9I, nucleocapsid-D3L, membrane-Q19E, and membrane-A63T mutations might be responsible factors for the surge in the SARS-CoV-2 Omicron variant wave in Iran.

CONCLUSIONS

real-time genomic surveillance is crucial for detecting new SARS-CoV-2 variants, updating diagnostic tools, designing vaccines, and understanding adaptation to new environments.

Low Prevalence of Nirmatrelvir-Ritonavir Resistance-Associated Mutations in SARS-CoV-2 Lineages From Botswana

Background

We evaluated naturally occurring nirmatrelvir-ritonavir (NTV/r) resistance-associated mutations (RAMs) among severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) strains from Botswana, a country with no NTV/r use to date, in order to recommend the usage of the agent for high-risk patients with coronavirus disease 2019 (COVID-19).

Methods

We conducted a retrospective analysis using 5254 complete SARS-CoV-2 sequences from Botswana (September 2020-September 2023). We evaluated the mutational landscape of SARS-CoV-2 3-Chymotrypsin-like protease (3CLpro) relative to the highlighted list of RAMs granted Food and Drug Administration Emergency Use Authorization in 2023.

Results

The sequenced 5254 samples included Beta variants of concerns (VOCs; n = 323), Delta VOCs (n = 1314), and Omicron VOCs (n = 3354). Overall, 77.8% of the sequences exhibited at least 1 polymorphism within 76/306 amino acid positions in the nsp5 gene. NTV/rRAMs were identified in 34/5254 (0.65%; 95% CI, 0.43%-0.87%) and occurred at 5 distinct positions. Among the NTV/r RAMS detected, A191V was the most prevalent (24/34; 70.6%). Notably, T21I mutation had a prevalence of 20.6% (7/34) and coexisted with either K90R (n = 3) polymorphism in Beta sequences with RAMs or P132H (n = 3) polymorphism for Omicron sequences with RAMs. Other NTV/r RAMs detected included P108S, with a prevalence of 5.88% (2/34), and L50F, with a prevalence of 2.94% (1/34). NTV/r RAMs were significantly higher (P < .001) in Delta (24/35) compared with Beta (4/34) and Omicron (6/34) sequences.

Conclusions

The frequency of NTV/r RAMs in Botswana was low. Higher rates were observed in Delta VOCs compared to Omicron and Beta VOCs. As NTV/r use expands globally, continuous surveillance for drug-resistant variants is essential, given the RAMs identified in our study.

Quinazolines and thiazolidine-2,4-dions as SARS-CoV-2 inhibitors: repurposing, in silico molecular docking and dynamics simulation

This paper presents an extensive analysis of COVID-19 with a specific focus on VEGFR-2 inhibitors as potential treatments. The investigation includes an overview of computational methodologies employed in drug repurposing and highlights in silico research aimed at developing treatments for SARS-CoV-2. The study explores the possible effects of twenty-eight established VEGFR-2 inhibitors, which include amide and urea linkers, against SARS-CoV-2. Among these, nine inhibitors exhibit highly promising in silico outcomes (designated as 3-6, 11, 24, 26, 27, and sorafenib) and are subjected to extensive molecular dynamics (MD) simulations to evaluate the binding modes and affinities of these inhibitors to the SARS-CoV-2 Mpro across a 100 ns timeframe. Additionally, MD simulations are conducted to ascertain the binding free energy of the most compelling ligand-pocket complexes identified through docking studies. The findings provide valuable understanding regarding the dynamic and thermodynamic properties of the interactions between ligands and pockets, reinforcing the outcomes of the docking studies and presenting promising prospects for the creation of therapeutic treatments targeting COVID-19.

Molecular characterization of SARS-CoV-2 Omicron clade and clinical presentation in children

Since its emergence, SARS-CoV-2 Omicron clade has shown a marked degree of variability and different clinical presentation compared with previous clades. Here we demonstrate that at least four Omicron lineages circulated in children since December 2021, and studied until November 2022: BA.1 (33.6%), BA.2 (40.6%), BA.5 (23.7%) and BQ.1 (2.1%). At least 70% of infections concerned children under 1 year, most of them being infected with BA.2 lineages (n = 201, 75.6%). Looking at SARS-CoV-2 genetic variability, 69 SNPs were found to be significantly associated in pairs, (phi <  - 0.3 or > 0.3 and p-value < 0.001). 16 SNPs were involved in 4 distinct clusters (bootstrap > 0.75). One of these clusters (A23040G, A27259C, T23617G, T23620G) was also positively associated with moderate/severe COVID-19 presentation (AOR [95% CI] 2.49 [1.26-4.89] p-value: 0.008) together with comorbidities (AOR [95% CI] 2.67 [1.36-5.24] p-value: 0.004). Overall, these results highlight the extensive SARS-CoV-2 Omicron circulation in children, mostly aged < 1 year, and provide insights on viral diversification even considering low-abundant SNPs, finally suggesting the potential contribution of viral diversification in affecting disease severity.

Molnupiravir maintains antiviral activity against SARS-CoV-2 variants and exhibits a high barrier to the development of resistance

Molnupiravir, an oral prodrug of N-hydroxycytidine (NHC), previously demonstrated broad in vitro antiviral activity against multiple RNA viruses and has shown a high barrier to the development of resistance. Here, we present the antiviral activity of NHC against recent SARS-CoV-2 variants and the results of resistance selection studies to better understand the potential for viral resistance to NHC. NHC activity against SARS-CoV-2 variants omicron (BA.1, BA.1.1, BA.2, BA.4, BA.4.6, BA.5, BQ.1.1, XBB.1, and XBB.1.5), alpha (B.1.1.7), beta (B.1.351), gamma (P.1), delta (B.1.617.2), lambda (C.37), and mu (B.1.621) was evaluated in Vero E6 cells using cytopathic effect assays. Resistance selection studies were performed by passaging SARS-CoV-2 (WA1) in the presence of NHC or a 3C-like protease inhibitor (MRK-A) in Vero E6 cells. Supernatants from cultures exhibiting a cytopathic effect score of ≥2 were re-passaged, and IC50 values were estimated. Whole-genome deep sequencing was performed on viral RNA isolated at each passage. NHC demonstrated similar potency against all SARS-CoV-2 variants evaluated. No evidence of SARS-CoV-2 phenotypic or genotypic resistance to NHC was observed following 30 passages. A random pattern of nucleotide changes was observed in NHC cultures, consistent with the drug's mechanism of action. In contrast, resistance was readily selected in all three MRK-A control cultures with the selection of a T21I substitution in the 3C-like protease. In conclusion, molnupiravir maintains antiviral activity across all major SARS-CoV-2 variants. Furthermore, no evidence of viral resistance to NHC was observed, supporting previous reports that NHC has a high barrier to developing resistance.

SARS-CoV-2: analysis of the effects of mutations in non-structural proteins

A worldwide pandemic that started in China in late 2019 was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a single-stranded RNA virus belonging to the family Coronaviridae. Due to its structural variability and mutability, this virus continues to evolve and pose a major health threat around the world. Its characteristics, such as transmissibility, antigenicity, and resistance to drugs and vaccines, are continually altered through mutations. Examining mutational hotspots and their structural repercussions can thus aid in the development of more-effective vaccinations and treatment plans. In this study, we used full genome sequences of SARS-CoV-2 variants to predict structural changes in viral proteins. These sequences were obtained from the Global Initiative on Sharing Avian Influenza Data (GISAID), and a set of significant mutations were identified in each of the non-structural proteins (NSP1-16) and structural proteins, including the envelope, nucleocapsid, membrane, and spike proteins. The mutations were characterized as stabilizing or destabilizing based on their effect on protein dynamics and stability, and their impact on structure and function was evaluated. Among all of the proteins, NSP6 stands out as especially variable. The results of this study augment our understanding of how mutational events influence virus pathogenicity and evolution.

SARS-CoV-2 Mpro Protease Variants of Concern Display Altered Viral Substrate and Cell Host Target Galectin-8 Processing but Retain Sensitivity toward Antivirals

The main protease of SARS-CoV-2 (M^pro) is the most promising drug target against coronaviruses due to its essential role in virus replication. With newly emerging variants there is a concern that mutations in M^pro may alter the structural and functional properties of protease and subsequently the potency of existing and potential antivirals. We explored the effect of 31 mutations belonging to 5 variants of concern (VOCs) on catalytic parameters and substrate specificity, which revealed changes in substrate binding and the rate of cleavage of a viral peptide. Crystal structures of 11 M^pro mutants provided structural insight into their altered functionality. Additionally, we show M^pro mutations influence proteolysis of an immunomodulatory host protein Galectin-8 (Gal-8) and a subsequent significant decrease in cytokine secretion, providing evidence for alterations in the escape of host-antiviral mechanisms. Accordingly, mutations associated with the Gamma VOC and highly virulent Delta VOC resulted in a significant increase in Gal-8 cleavage. Importantly, IC50s of nirmatrelvir (Pfizer) and our irreversible inhibitor AVI-8053 demonstrated no changes in potency for both drugs for all mutants, suggesting M^pro will remain a high-priority antiviral drug candidate as SARS-CoV-2 evolves.

Comparative Immunoreactivity Analyses of Hantaan Virus Glycoprotein-Derived MHC-I Epitopes in Vaccination

MHC-I antigen processes and presentation trigger host-specific anti-viral cellular responses during infection, in which epitope-recognizing cytotoxic T lymphocytes eliminate infected cells and contribute to viral clearance through a cytolytic killing effect. In this study, Hantaan virus (HTNV) GP-derived 9-mer dominant epitopes were obtained with high affinity to major HLA-I and H-2 superfamilies. Further immunogenicity and conservation analyses selected 11 promising candidates, and molecule docking (MD) was then simulated with the corresponding MHC-I alleles. Two-way hierarchical clustering revealed the interactions between GP peptides and MHC-I haplotypes. Briefly, epitope hotspots sharing good affinity to a wide spectrum of MHC-I molecules highlighted the biomedical practice for vaccination, and haplotype clusters represented the similarities among individuals during T-cell response establishment. Cross-validation proved the patterns observed through both MD simulation and public data integration. Lastly, 148 HTNV variants yielded six types of major amino acid residue replacements involving four in nine hotspots, which minimally influenced the general potential of MHC-I superfamily presentation. Altogether, our work comprehensively evaluates the pan-MHC-I immunoreactivity of HTNV GP through a state-of-the-art workflow in light of comparative immunology, acknowledges present discoveries, and offers guidance for ongoing HTNV vaccine pursuit.