Emerging Mutations in Nsp1 of SARS-CoV-2 and Their Effect on the Structural Stability

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.

SARS-CoV-2 population dynamics in immunocompetent individuals in a closed transmission chain shows genomic diversity over the course of infection

BACKGROUND

SARS-CoV-2 remains rapidly evolving, and many biologically important genomic substitutions/indels have characterised novel SARS-CoV-2 lineages, which have emerged during successive global waves of the pandemic. Worldwide genomic sequencing has been able to monitor these waves, track transmission clusters, and examine viral evolution in real time to help inform healthcare policy. One school of thought is that an apparent greater than average divergence in an emerging lineage from contemporary variants may require persistent infection, for example in an immunocompromised host. Due to the nature of the COVID-19 pandemic and sampling, there were few studies that examined the evolutionary trajectory of SARS-CoV-2 in healthy individuals.

METHODS

We investigated viral evolutionary trends and participant symptomatology within a cluster of 16 SARS-CoV-2 infected, immunocompetent individuals with no co-morbidities in a closed transmission chain. Longitudinal nasopharyngeal swab sampling allowed characterisation of SARS-CoV-2 intra-host variation over time at both the dominant and minor genomic variant levels through Nimagen-Illumina sequencing.

RESULTS

A change in viral lineage assignment was observed in individual infections; however, there was only one indel and no evidence of recombination over the period of an acute infection. Minor and dominant genomic modifications varied between participants, with some minor genomic modifications increasing in abundance to become the dominant viral sequence during infection.

CONCLUSIONS

Data from this cohort of SARS-CoV-2-infected participants demonstrated that long-term persistent infection in an immunocompromised host was not necessarily a prerequisite for generating a greater than average frequency of amino acid substitutions. Amino acid substitutions at both the dominant and minor genomic sequence level were observed in immunocompetent individuals during infection showing that viral lineage changes can occur generating viral diversity.

Progressive Evolutionary Dynamics of Gene-Specific ω Led to the Emergence of Novel SARS-CoV-2 Strains Having Super-Infectivity and Virulence with Vaccine Neutralization

An estimation of the proportion of nonsynonymous to synonymous mutation (dn/ds, ω) of the SARS-CoV-2 genome would indicate the evolutionary dynamics necessary to evolve into novel strains with increased infection, virulence, and vaccine neutralization. A temporal estimation of ω of the whole genome, and all twenty-nine SARS-CoV-2 genes of major virulent strains of alpha, delta and omicron demonstrates that the SARS-CoV-2 genome originally emerged (ω ~ 0.04) with a strong purifying selection (ω < 1) and reached (ω ~ 0.85) in omicron towards diversifying selection (ω > 1). A marked increase in the ω occurred in the spike gene from alpha (ω = 0.2) to omicron (ω = 1.97). The ω of the replication machinery genes including RDRP, NSP3, NSP4, NSP7, NSP8, NSP10, NSP13, NSP14, and ORF9 are markedly increased, indicating that these genes/proteins are yet to be evolutionary stabilized and are contributing to the evolution of novel virulent strains. The delta-specific maximum increase in ω in the immunomodulatory genes of NSP8, NSP10, NSP16, ORF4, ORF5, ORF6, ORF7A, and ORF8 compared to alpha or omicron indicates delta-specific vulnerabilities for severe COVID-19 related hospitalization and death. The maximum values of ω are observed for spike (S), NSP4, ORF8 and NSP15, which indicates that the gene-specific temporal estimation of ω identifies specific genes for its super-infectivity and virulency that could be targeted for drug development.

SARS-CoV-2 targets ribosomal RNA biogenesis

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) hinders host gene expression, curbing defenses and licensing viral protein synthesis and virulence. During SARS-CoV-2 infection, the virulence factor non-structural protein 1 (Nsp1) targets the mRNA entry channel of mature cytoplasmic ribosomes, limiting translation. We show that Nsp1 also restrains translation by targeting nucleolar ribosome biogenesis. SARS-CoV-2 infection disrupts 18S and 28S ribosomal RNA (rRNA) processing. Expression of Nsp1 recapitulates the processing defects. Nsp1 abrogates rRNA production without altering the expression of critical processing factors or nucleolar organization. Instead, Nsp1 localizes to the nucleolus, interacting with precursor-rRNA and hindering its maturation separately from the viral protein's role in restricting mature ribosomes. Thus, SARS-CoV-2 Nsp1 limits translation by targeting ribosome biogenesis and mature ribosomes. These findings revise our understanding of how SARS-CoV-2 Nsp1 controls human protein synthesis, suggesting that efforts to counter Nsp1's effect on translation should consider the protein's impact from ribosome manufacturing to mature ribosomes.

Elucidating the binding mechanism of SARS-CoV-2 NSP6-TBK1 and structure-based designing of phytocompounds inhibitors for instigating the host immune response

SARS-CoV-2, also referred to as severe acute respiratory syndrome coronavirus 2, is the virus responsible for causing COVID-19, an infectious disease that emerged in Wuhan, China, in December 2019. Among its crucial functions, NSP6 plays a vital role in evading the human immune system by directly interacting with a receptor called TANK-binding kinase (TBK1), leading to the suppression of IFNβ production. Consequently, in the present study we used the structural and biophysical approaches to analyze the effect of newly emerged mutations on the binding of NSP6 and TBK1. Among the identified mutations, four (F35G, L37F, L125F, and I162T) were found to significantly destabilize the structure of NSP6. Furthermore, the molecular docking analysis highlighted that the mutant NSP6 displayed its highest binding affinity with TBK1, exhibiting docking scores of -1436.2 for the wildtype and -1723.2, -1788.6, -1510.2, and -1551.7 for the F35G, L37F, L125F, and I162T mutants, respectively. This suggests the potential for an enhanced immune system evasion capability of NSP6. Particularly, the F35G mutation exhibited the strongest binding affinity, supported by a calculated binding free energy of -172.19 kcal/mol. To disrupt the binding between NSP6 and TBK1, we conducted virtual drug screening to develop a novel inhibitor derived from natural products. From this screening, we identified the top 5 hit compounds as the most promising candidates with a docking score of -6.59 kcal/mol, -6.52 kcal/mol, -6.32 kcal/mol, -6.22 kcal/mol, and -6.21 kcal/mol. The molecular dynamic simulation of top 3 hits further verified the dynamic stability of drugs-NSP6 complexes. In conclusion, this study provides valuable insight into the higher infectivity of the SARS-CoV-2 new variants and a strong rationale for the development of novel drugs against NSP6.

Emerging variants of SARS-CoV-2 NSP10 highlight strong functional conservation of its binding to two non-structural proteins, NSP14 and NSP16

The coronavirus SARS-CoV-2 protects its RNA from being recognized by host immune responses by methylation of its 5' end, also known as capping. This process is carried out by two enzymes, non-structural protein 16 (NSP16) containing 2'-O-methyltransferase and NSP14 through its N7 methyltransferase activity, which are essential for the replication of the viral genome as well as evading the host's innate immunity. NSP10 acts as a crucial cofactor and stimulator of NSP14 and NSP16. To further understand the role of NSP10, we carried out a comprehensive analysis of >13 million globally collected whole-genome sequences (WGS) of SARS-CoV-2 obtained from the Global Initiative Sharing All Influenza Data (GISAID) and compared it with the reference genome Wuhan/WIV04/2019 to identify all currently known variants in NSP10. T12I, T102I, and A104V in NSP10 have been identified as the three most frequent variants and characterized using X-ray crystallography, biophysical assays, and enhanced sampling simulations. In contrast to other proteins such as spike and NSP6, NSP10 is significantly less prone to mutation due to its crucial role in replication. The functional effects of the variants were examined for their impact on the binding affinity and stability of both NSP14-NSP10 and NSP16-NSP10 complexes. These results highlight the limited changes induced by variant evolution in NSP10 and reflect on the critical roles NSP10 plays during the SARS-CoV-2 life cycle. These results also indicate that there is limited capacity for the virus to overcome inhibitors targeting NSP10 via the generation of variants in inhibitor binding pockets.

Insights into the evolution of mutations in SARS-CoV-2 non-spike proteins

The COVID-19 pandemic has been driven by the emergence of SARS-CoV-2 variants with mutations across all the viral proteins. Although mutations in the spike protein have received significant attention, understanding the prevalence and potential impact of mutations in other viral proteins is essential for comprehending the evolution of SARS-CoV-2. Here, we conducted a comprehensive analysis of approximately 14 million sequences of SARS-CoV-2 deposited in the GISAID database until December 2022 to identify prevalent mutations in the non-spike proteins at the global and country levels. Additionally, we evaluated the energetics of each mutation to better understand their impact on protein stability. While the consequences of many mutations remain unclear, we discuss potential structural and functional significance of some mutations. Our study highlights the ongoing evolutionary process of SARS-CoV-2 and underscores the importance of understanding changes in non-spike proteins.

Biochemical and HDX Mass Spectral Characterization of the SARS-CoV-2 Nsp1 Protein

A major challenge in defining the pathophysiology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is to better understand virally encoded multifunctional proteins and their interactions with host factors. Among the many proteins encoded by the positive-sense, single-stranded RNA genome, nonstructural protein 1 (Nsp1) stands out due to its impact on several stages of the viral replication cycle. Nsp1 is the major virulence factor that inhibits mRNA translation. Nsp1 also promotes host mRNA cleavage to modulate host and viral protein expression and to suppress host immune functions. To better define how this multifunctional protein can facilitate distinct functions, we characterize SARS-CoV-2 Nsp1 by using a combination of biophysical techniques, including light scattering, circular dichroism, hydrogen/deuterium exchange mass spectrometry (HDX-MS), and temperature-dependent HDX-MS. Our results reveal that the SARS-CoV-2 Nsp1 N- and C-terminus are unstructured in solution, and in the absence of other proteins, the C-terminus has an increased propensity to adopt a helical conformation. In addition, our data indicate that a short helix exists near the C-terminus and adjoins the region that binds the ribosome. Together, these findings provide insights into the dynamic nature of Nsp1 that impacts its functions during infection. Furthermore, our results will inform efforts to understand SARS-CoV-2 infection and antiviral development.

A novel mutation-proof, next-generation vaccine to fight against upcoming SARS-CoV-2 variants and subvariants, designed through AI enabled approaches and tools, along with the machine learning based immune simulation: A vaccine breakthrough

Emerging SARS-CoV-2 variants and subvariants are great concerns for their significant mutations, which are also responsible for vaccine escape. Therefore, the study was undertaken to develop a mutation-proof, next-generation vaccine to protect against all upcoming SARS-CoV-2 variants. We used advanced computational and bioinformatics approaches to develop a multi-epitopic vaccine, especially the AI model for mutation selection and machine learning (ML) strategies for immune simulation. AI-enabled and the top-ranked antigenic selection approaches were used to select nine mutations from 835 RBD mutations. We selected twelve common antigenic B cell and T cell epitopes (CTL and HTL) containing the nine RBD mutations and joined them with the adjuvants, PADRE sequence, and suitable linkers. The constructs' binding affinity was confirmed through docking with TLR4/MD2 complex and showed significant binding free energy (-96.67 kcal mol^-1) with positive binding affinity. Similarly, the calculated eigenvalue (2.428517e-05) from the NMA of the complex reveals proper molecular motion and superior residues' flexibility. Immune simulation shows that the candidate can induce a robust immune response. The designed mutation-proof, multi-epitopic vaccine could be a remarkable candidate for upcoming SARS-CoV-2 variants and subvariants. The study method might guide researchers in developing AI-ML and immunoinformatics-based vaccines for infectious disease.

SARS-CoV-2 Non-structural protein 1(NSP1) mutation virulence and natural selection: Evolutionary trends in the six continents

OBJECTIVE

Rapid transmission and reproduction of RNA viruses prepare conducive conditions to have a high rate of mutations in their genetic sequence. The viral mutations make adapt the severe acute respiratory syndrome coronavirus 2 in the host environment and help the evolution of the virus then also caused a high mortality rate by the virus that threatens worldwide health. Mutations and adaptation help the virus to escape confrontations that are done against it.

METHODS

In the present study, we analyzed 6,510,947 sequences of non-structural protein 1 as one of the conserved regions of the virus to find out frequent mutations and substitute amino acids in comparison with the wild type. NSP1 mutations rate divided into continents were different.

RESULTS

Based on this continental categorization, E87D in global vision and also in Europe notably increased. The E87D mutation has signed up to January 2022 as the first frequent mutation observed. The remarkable mutations, H110Y and R24C have the second and third frequencies, respectively.

CONCLUSION

According to the important role of non-structural protein 1 on the host mRNA translation, developing drug design against the protein could be so hopeful to find more effective ways the control and then treatment of the global pandemic coronavirus disease 2019.

ViralVar: A Web Tool for Multilevel Visualization of SARS-CoV-2 Genomes

The unprecedented growth of publicly available SARS-CoV-2 genome sequence data has increased the demand for effective and accessible SARS-CoV-2 data analysis and visualization tools. The majority of the currently available tools either require computational expertise to deploy them or limit user input to preselected subsets of SARS-CoV-2 genomes. To address these limitations, we developed ViralVar, a publicly available, point-and-click webtool that gives users the freedom to investigate and visualize user-selected subsets of SARS-CoV-2 genomes obtained from the GISAID public database. ViralVar has two primary features that enable: (1) the visualization of the spatiotemporal dynamics of SARS-CoV-2 lineages and (2) a structural/functional analysis of genomic mutations. As proof-of-principle, ViralVar was used to explore the evolution of the SARS-CoV-2 pandemic in the USA in pediatric, adult, and elderly populations (n > 1.7 million genomes). Whereas the spatiotemporal dynamics of the variants did not differ between these age groups, several USA-specific sublineages arose relative to the rest of the world. Our development and utilization of ViralVar to provide insights on the evolution of SARS-CoV-2 in the USA demonstrates the importance of developing accessible tools to facilitate and accelerate the large-scale surveillance of circulating pathogens.

Two Ligand-Binding Sites on SARS-CoV-2 Non-Structural Protein 1 Revealed by Fragment-Based X-ray Screening

The regular reappearance of coronavirus (CoV) outbreaks over the past 20 years has caused significant health consequences and financial burdens worldwide. The most recent and still ongoing novel CoV pandemic, caused by Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) has brought a range of devastating consequences. Due to the exceptionally fast development of vaccines, the mortality rate of the virus has been curbed to a significant extent. However, the limitations of vaccination efficiency and applicability, coupled with the still high infection rate, emphasise the urgent need for discovering safe and effective antivirals against SARS-CoV-2 by suppressing its replication or attenuating its virulence. Non-structural protein 1 (nsp1), a unique viral and conserved leader protein, is a crucial virulence factor for causing host mRNA degradation, suppressing interferon (IFN) expression and host antiviral signalling pathways. In view of the essential role of nsp1 in the CoV life cycle, it is regarded as an exploitable target for antiviral drug discovery. Here, we report a variety of fragment hits against the N-terminal domain of SARS-CoV-2 nsp1 identified by fragment-based screening via X-ray crystallography. We also determined the structure of nsp1 at atomic resolution (0.99 Å). Binding affinities of hits against nsp1 and potential stabilisation were determined by orthogonal biophysical assays such as microscale thermophoresis and thermal shift assays. We identified two ligand-binding sites on nsp1, one deep and one shallow pocket, which are not conserved between the three medically relevant SARS, SARS-CoV-2 and MERS coronaviruses. Our study provides an excellent starting point for the development of more potent nsp1-targeting inhibitors and functional studies on SARS-CoV-2 nsp1.

Evolution of ACE2-independent SARS-CoV-2 infection and mouse adaption after passage in cells expressing human and mouse ACE2

Human ACE2 Human angiotensin converting enzyme 2 (hACE2) is the key cell attachment and entry receptor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with the original SARS-CoV-2 isolates unable to use mouse ACE2 (mACE2). Herein we describe the emergence of a SARS-CoV-2 strain capable of ACE2-independent infection and the evolution of mouse-adapted (MA) SARS-CoV-2 by in vitro serial passaging of virus in co-cultures of cell lines expressing hACE2 and mACE2. MA viruses evolved with up to five amino acid changes in the spike protein, all of which have been seen in human isolates. MA viruses replicated to high titers in C57BL/6J mouse lungs and nasal turbinates and caused characteristic lung histopathology. One MA virus also evolved to replicate efficiently in several ACE2-negative cell lines across several species, including clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9) ACE2 knockout cells. An E484D substitution is likely involved in ACE2-independent entry and has appeared in only ≈0.003 per cent of human isolates globally, suggesting that it provided no significant selection advantage in humans. ACE2-independent entry reveals a SARS-CoV-2 infection mechanism that has potential implications for disease pathogenesis, evolution, tropism, and perhaps also intervention development.

SARS-CoV-2 Non-Structural Protein 1(NSP1) Mutation Virulence and Natural Selection: Evolutionary Trends in the Six Continents

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an unsegmented positivesense single-stranded RNA virus that belongs to the β-coronavirus . This virus was the cause of a novel severe acute respiratory syndrome in 2019 (COVID-19) that emerged in Wuhan, China at the early stage of the pandemic and rapidly spread around the world. Rapid transmission and reproduction of SARS-CoV-2 threaten worldwide health with a high mortality rate from the virus. According to the significant role of non-structural protein 1 (NSP1) in inhibiting host mRNA translation, this study focuses on the link between amino acid sequences of NSP1 and alterations of them spreading around the world. The SARS-CoV-2 NSP1 protein sequences were analyzed and FASTA files were processed by Python language programming libraries. Reference sequences compared with each NSP1 sample to identify every mutation and categorize them were based on continents and frequencies. NSP1 mutations rate divided into continents were different. Based on continental studies, E87D in global vision and also in Europe notably increased. The E87D mutation has significantly risen especially in the last months of the study as the first frequent mutation observed. The remarkable mutations, H110Y and R24C, have the second and third frequencies, respectively. Based on this mutational information, despite NSP1 being a conserved sequence occurrence, these mutations change the rate of flexibility and stability of the NSP1 protein, which can eventually affect inhibiting the host translation.

IMPORTANCE

In this study, we analyzed 6,510,947 sequences of non-structural protein 1 as a conserved region of SARS-CoV-2. According to the obtained results, 93.4819% of samples had no mutant regions on their amino acid sequences. Heat map data of mutational samples demonstrated high percentages of mutations that occurred in the region of 72 to 126 amino acids indicating a hot spot region of the protein. Increased rates of E87D, H110Y, and R24C mutations in the timeline of our study were reported as significant compared to available mutant samples. Analyzing the details of replacing amino acids in the most frequent E87D mutation reveals the role of this alteration in increasing molecule flexibility and destabilizing the structure of the protein.

Identifying Markers of Emerging SARS-CoV-2 Variants in Patients With Secondary Immunodeficiency

Since the end of 2019, the world has been challenged by the coronavirus disease 2019 (COVID-19) pandemic. With COVID-19 cases rising globally, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to evolve, resulting in the emergence of variants of interest (VOI) and of concern (VOC). Of the hundreds of millions infected, immunodeficient patients are one of the vulnerable cohorts that are most susceptible to this virus. These individuals include those with preexisting health conditions and/or those undergoing immunosuppressive treatment (secondary immunodeficiency). In these cases, several researchers have reported chronic infections in the presence of anti-COVID-19 treatments that may potentially lead to the evolution of the virus within the host. Such variations occurred in a variety of viral proteins, including key structural ones involved in pathogenesis such as spike proteins. Tracking and comparing such mutations with those arisen in the general population may provide information about functional sites within the SARS-CoV-2 genome. In this study, we reviewed the current literature regarding the specific features of SARS-CoV-2 evolution in immunocompromised patients and identified recurrent de novo amino acid changes in virus isolates of these patients that can potentially play an important role in SARS-CoV-2 pathogenesis and evolution.

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.

Drug targeting Nsp1-ribosomal complex shows antiviral activity against SARS-CoV-2

The SARS-CoV-2 non-structural protein 1 (Nsp1) contains an N-terminal domain and C-terminal helices connected by a short linker region. The C-terminal helices of Nsp1 (Nsp1-C-ter) from SARS-CoV-2 bind in the mRNA entry channel of the 40S ribosomal subunit and blocks mRNA entry, thereby shutting down host protein synthesis. Nsp1 suppresses host immune function and is vital for viral replication. Hence, Nsp1 appears to be an attractive target for therapeutics. In this study, we have in silico screened Food and Drug Administration (FDA)-approved drugs against Nsp1-C-ter. Among the top hits obtained, montelukast sodium hydrate binds to Nsp1 with a binding affinity (K_D) of 10.8 ± 0.2 µM in vitro. It forms a stable complex with Nsp1-C-ter in simulation runs with –95.8 ± 13.3 kJ/mol binding energy. Montelukast sodium hydrate also rescues the inhibitory effect of Nsp1 in host protein synthesis, as demonstrated by the expression of firefly luciferase reporter gene in cells. Importantly, it shows antiviral activity against SARS-CoV-2 with reduced viral replication in HEK cells expressing ACE2 and Vero-E6 cells. We, therefore, propose montelukast sodium hydrate can be used as a lead molecule to design potent inhibitors to help combat SARS-CoV-2 infection.