The importance of accessory protein variants in the pathogenicity of SARS-CoV-2.

The coronavirus disease 2019 (COVID-19) is caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS- CoV-2) with an estimated fatality rate of less than 1%. The SARS-CoV-2 accessory proteins ORF3a, ORF6, ORF7a, ORF7b, ORF8, and ORF10 possess putative functions to manipulate host immune mechanisms. These involve interferons, which appear as a consensus function, immune signaling receptor NLRP3 (NLR family pyrin domain-containing 3) inflammasome, and inflammatory cytokines such as interleukin 1ß (IL-1ß) and are critical in COVID-19 pathology. Outspread variations of each of the six accessory proteins were observed across six continents of all complete SARS-CoV-2 proteomes based on the data reported before November 2020. A decreasing order of percentage of unique variations in the accessory proteins was determined as ORF3a > ORF8 > ORF7a > ORF6 > ORF10 > ORF7b across all continents. The highest and lowest unique variations of ORF3a were observed in South America and Oceania, respectively. These findings suggest that the wide variations in accessory proteins seem to affect the pathogenicity of SARS-CoV-2.

Periodically aperiodic pattern of SARS-CoV-2 mutations underpins the uncertainty of its origin and evolution.

Various lineages of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) have contributed to prolongation of the Coronavirus Disease 2019 (COVID-19) pandemic. Several non-synonymous mutations in SARS-CoV-2 proteins have generated multiple SARS-CoV-2 variants. In our previous report, we have shown that an evenly uneven distribution of unique protein variants of SARS-CoV-2 is geo-location or demography-specific. However, the correlation between the demographic transmutability of the SARS-CoV-2 infection and mutations in various proteins remains unknown due to hidden symmetry/asymmetry in the occurrence of mutations. This study tracked how these mutations are emerging in SARS-CoV-2 proteins in six model countries and globally. In a geo-location, considering the mutations having a frequency of detection of at least 500 in each SARS-CoV-2 protein, we studied the country-wise percentage of invariant residues. Our data revealed that since October 2020, highly frequent mutations in SARS-CoV-2 have been observed mostly in the Open Reading Frame (ORF) 7b and ORF8, worldwide. No such highly frequent mutations in any of the SARS-CoV-2 proteins were found in the UK, India, and Brazil, which does not correlate with the degree of transmissibility of the virus in India and Brazil. However, we have found a signature that SARS-CoV-2 proteins were evolving at a higher rate, and considering global data, mutations are detected in the majority of the available amino acid locations. Fractal analysis of each protein's normalized factor time series showed a periodically aperiodic emergence of dominant variants for SARS-CoV-2 protein mutations across different countries. It was noticed that certain high-frequency variants have emerged in the last couple of months, and thus the emerging SARS-CoV-2 strains are expected to contain prevalent mutations in the ORF3a, membrane, and ORF8 proteins. In contrast to other beta-coronaviruses, SARS-CoV-2 variants have rapidly emerged based on demographically dependent mutations. Characterization of the periodically aperiodic nature of the demographic spread of SARS-CoV-2 variants in various countries can contribute to the identification of the origin of SARS-CoV-2.

Emerging genetic diversity of SARS-CoV-2 RNA dependent RNA polymerase (RdRp) alters its B-cell epitopes.

The RNA dependent RNA polymerase (RdRp) plays crucial role in virus life cycle by replicating the viral genome. The SARS-CoV-2 is an RNA virus that rapidly spread worldwide and acquired mutations. This study was carried out to identify mutations in RdRp as the SARS-CoV-2 spread in India. We compared 50217 RdRp sequences reported from India with the first reported RdRp sequence from Wuhan, China to identify 223 mutations acquired among Indian isolates. Our protein modelling study revealed that several mutants can potentially alter stability and flexibility of RdRp. We predicted the potential B cell epitopes contributed by RdRp and identified thirty-six linear continuous and twenty-five discontinuous epitopes. Among 223 RdRp mutants, 44% of them localises in the B cell epitopes region. Altogether, this study highlights the need to identify and characterize the variations in RdRp to understand the impact of these mutations on SARS-CoV-2.

Possible Transmission Flow of SARS-CoV-2 Based on ACE2 Features.

Angiotensin-converting enzyme 2 (ACE2) is the cellular receptor for the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) that is engendering the severe coronavirus disease 2019 (COVID-19) pandemic. The spike (S) protein receptor-binding domain (RBD) of SARS-CoV-2 binds to the three sub-domains viz. amino acids (aa) 22-42, aa 79-84, and aa 330-393 of ACE2 on human cells to initiate entry. It was reported earlier that the receptor utilization capacity of ACE2 proteins from different species, such as cats, chimpanzees, dogs, and cattle, are different. A comprehensive analysis of ACE2 receptors of nineteen species was carried out in this study, and the findings propose a possible SARS-CoV-2 transmission flow across these nineteen species.

Partial purification of histone H3 proteolytic activity from the budding yeast Saccharomyces cerevisiae.

The proteolytic clipping of histone tails has recently emerged as a novel form of irreversible post-translational modification (PTM) of histones. Histone clipping has been implicated as a regulatory process leading to the permanent removal of PTMs from histone proteins. However, there is scarcity of literature that describes the identification and characterization of histone-specific proteases. Here, we employed various biochemical methods to report histone H3-specific proteolytic activity from budding yeast. Our results demonstrate that H3 proteolytic activity was associated with sepharose bead matrices and activity was not affected by a variety of stress conditions. We have also identified the existence of an unknown protein that acts as a physiological inhibitor of the H3-clipping activity of yeast H3 protease. Moreover, through protease inhibition assays, we have also characterized yeast H3 protease as a serine protease. Interestingly, unlike glutamate dehydrogenase (GDH), yeast H3 proteolytic activity was not inhibited by Stefin B. Together, our findings suggest the existence of a novel H3 protease in yeast that is different from other reported histone H3 proteases. The presence of histone H3 proteolytic activity, along with the physiological inhibitor in yeast, suggests an interesting molecular mechanism that regulates the activity of histone proteases. Copyright © 2016 John Wiley & Sons, Ltd.

Signaling of chloroquine-induced stress in the yeast Saccharomyces cerevisiae requires the Hog1 and Slt2 mitogen-activated protein kinase pathways.

Chloroquine (CQ) has been under clinical use for several decades, and yet little is known about CQ sensing and signaling mechanisms or about their impact on various biological pathways. We employed the budding yeast Saccharomyces cerevisiae as a model organism to study the pathways targeted by CQ. Our screening with yeast mutants revealed that it targets histone proteins and histone deacetylases (HDACs). Here, we also describe the novel role of mitogen-activated protein kinases Hog1 and Slt2, which aid in survival in the presence of CQ. Cells deficient in Hog1 or Slt2 are found to be CQ hypersensitive, and both proteins were phosphorylated in response to CQ exposure. CQ-activated Hog1p is translocated to the nucleus and facilitates the expression of GPD1 (glycerol-3-phosphate dehydrogenase), which is required for the synthesis of glycerol (one of the major osmolytes). Moreover, cells treated with CQ exhibited an increase in intracellular reactive oxygen species (ROS) levels and the effects were rescued by addition of reduced glutathione to the medium. The deletion of SOD1, the superoxide dismutase in yeast, resulted in hypersensitivity to CQ. We have also observed P38 as well as P42/44 phosphorylation in HEK293T human cells upon exposure to CQ, indicating that the kinds of responses generated in yeast and human cells are similar. In summary, our findings define the multiple biological pathways targeted by CQ that might be useful for understanding the toxicity modulated by this pharmacologically important molecule.

Mitogen-activated protein kinase Hog1 is activated in response to curcumin exposure in the budding yeast Saccharomyces cerevisiae.

BACKGROUND: Curcumin (CUR), an active polyphenol derived from the spice turmeric, has been traditionally used for centuries in ancient Indian medicine to treat a number of diseases. The physiological effects of CUR have been shown to be diverse; however, the target molecules and pathways that CUR affects have yet to be fully described. RESULTS: Here, we demonstrate for the first time that the budding yeast mitogen-activated protein kinase (MAPK) Hog1 is essential for the response to CUR. Moreover, CUR-induced Hog1 phosphorylation was rescued by supplementation of iron to the growth medium. Hog1 was rapidly phosphorylated upon CUR treatment, but unlike the response to hyperosmotic shock (0.8 M NaCl), it remains activated for an extended period of time. A detailed analysis of HOG pathway mutants revealed that Pbs2p, Ptc2p, and Ssk2p are required for optimal CUR-induced Hog1 phosphorylation. We also observed a Hog1 dependent transcriptional response to CUR treatment that involved the up-regulation of glycerol-3-phosphate dehydrogenase 1 (GPD1), a factor that is essential for the hyperosmotic stress response. CONCLUSIONS: Our present finding revealed the role of Hog1 MAPK in regulation of CUR-induced transcriptional response. We anticipate that our finding will enhance the understanding on the molecular mode of action of CUR on S. cerevisiae.

An ebselen like catalyst with enhanced GPx activity via a selenol intermediate.

The reaction of KSeO(t)Bu with 2-iodo-arylbenzamides gave benzamide ring-substituted, quinine-derived isoselenazolones 1b­1d. The reaction of PhSH with ortho-methyl-substituted isoselenazolone 1b gave selenol 3b, which is oxidized by H2O2 to regenerate 1b. Isoselenazolone 1b shows a high rate (0.33 × 103 µM min(−1)) of oxidation of PhSH with H2O2, which is ∼103-fold more active than ebselen (1a) and ≥30-fold more active than the other isoselenazolones of this study. Compound 1b shows less inhibition of the growth of yeast cells than 1a.

Ebselen, a promising antioxidant drug: mechanisms of action and targets of biological pathways.

Ebselen, an organoselenium compound, mimics glutathione peroxidase activity. It is a multifunctional compound, which catalyzes several essential reactions for the protection of cellular components from oxidative and free radical damage. Based on a number of in vitro and in vivo studies, various mechanisms are proposed to understand the biomedical actions of ebselen in health and diseases. It modulates metallo-proteins, enzymatic cofactors, gene expression, epigenetics, antioxidant defenses and immune systems. Owing to these properties, ebselen is currently under clinical trials for the prevention and treatment of various disorders such as cardiovascular diseases, arthritis, stroke, atherosclerosis, and cancer. A few ebselen-based pharmaceutical agents are under extensive investigation. As ebselen has been shown to have significant cellular toxicity, appropriate studies are needed to redesign the ebselen-based therapy for clinical trials. This review summarizes current understanding of the biochemical and molecular properties, and pharmacological applications of ebselen and future directions in this area of research.

Proteolytic clipping of histone tails: the emerging role of histone proteases in regulation of various biological processes.

Chromatin is a dynamic DNA scaffold structure that responds to a variety of external and internal stimuli to regulate the fundamental biological processes. Majority of the cases chromatin dynamicity is exhibited through chemical modifications and physical changes between DNA and histones. These modifications are reversible and complex signaling pathways involving chromatin-modifying enzymes regulate the fluidity of chromatin. Fluidity of chromatin can also be impacted through irreversible change, proteolytic processing of histones which is a poorly understood phenomenon. In recent studies, histone proteolysis has been implicated as a regulatory process involved in the permanent removal of epigenetic marks from histones. Activities responsible for clipping of histone tails and their significance in various biological processes have been observed in several organisms. Here, we have reviewed the properties of some of the known histone proteases, analyzed their significance in biological processes and have provided future directions.

Correction: Assessment of the Biological Pathways Targeted by Isocyanate Using N-Succinimidyl N-Methylcarbamate in Budding Yeast Saccharomyces cerevisiae.

[This corrects the article DOI: 10.1371/journal.pone.0092993.].

Alleviation of arsenic-induced neurobehavioral defects with selenium in the larvae of Zaprionus indianus.

Selenium is an essential antioxidative micronutrient. This study was conducted to characterize the arsenic toxicity induced on the African fig fly, Zaprionus indianus, and its possible amelioration by selenium. We used computational tools and in vivo experiments to elucidate the mechanism of action of arsenic and selenium on Z. indianus larvae. We conducted experiments to study neurobehavioral parameters including learning and memory ability test and crawling and contraction assays. Our in silico study revealed twelve primary targets of arsenic trioxide. The gene ontology annotation of primary and secondary targets of arsenic trioxide revealed selenocysteine metabolic processes as one of the most reliable targets. To validate our in silico data, we analyzed the effect of arsenic trioxide on larvae of Z. indianus and tested the possible amelioration by sodium selenite supplementation. Our data demonstrated that the arsenic trioxide deteriorated the learning and memory ability of 2nd instar larvae of Z. indianus and such effect was reversed by sodium selenite supplementation. Furthermore, crawling and contraction assay done on 3rd instar larvae showed that there was reduction in both parameters upon arsenic trioxide exposure, which was restored with sodium selenite supplementation. Altogether, our computational and in vivo results strongly indicated that the neurobehavioral defects induced by arsenic trioxide on the larvae of Z. indianus can be successfully alleviated in the presence of sodium selenite.

A computational study on mitogenome-encoded proteins of Pavo cristatus and Pavo muticus identifies key genetic variations with functional implications.

BACKGROUND: The Pavo cristatus population, native to the Indian subcontinent, is thriving well in India. However, the Pavo muticus population, native to the tropical forests of Southeast Asia, has reduced drastically and has been categorised as an endangered group. To understand the probable genetic factors associated with the decline of P. muticus, we compared the mitogenome-encoded proteins (13 proteins) between these two species. RESULTS: Our data revealed that the most frequent variant between these two species was mtND1, which had an alteration in 9.57% residues, followed by mtND5 and mtATP6. We extended our study on the rest of the proteins and observed that cytochrome c oxidase subunits 1, 2, and 3 do not have any change. The 3-dimensional structure of all 13 proteins was modeled using the Phyre2 programme. Our data show that most of the proteins are alpha helical, and the variations observed in P. muticus reside on the surface of the respective proteins. The effect of variation on protein function was also predicted, and our results show that amino acid substitution in mtND1 at 14 sites could be deleterious. Similarly, destabilising changes were observed in mtND1, 2, 3, 4, 5, and 6 and mtATP6-8 due to amino acid substitution in P. muticus. Furthermore, protein disorder scores were considerably altered in mtND1, 2, and 5 of P. muticus. CONCLUSIONS: The results presented here strongly suggest that variations in mitogenome-encoded proteins of P. cristatus and P. muticus may alter their structure and functions. Subsequently, these variations could alter energy production and may correlate with the decline in the population of P. muticus.

An immunoinformatics approach to study the epitopes of SARS-CoV-2 helicase, Nsp13.

Introduction and objective: Vaccines are administered worldwide to control on-going coronavirus disease-19 (COVID-19) pandemic caused by SARS-CoV-2. Vaccine efficacy is largely contributed by the epitopes present on the viral proteins and their alteration might help emerging variants to escape host immune surveillance. Therefore, this study was designed to study SARS-CoV-2 Nsp13 protein, its epitopes and evolution. Methods: Clustal Omega was used to identify mutations in Nsp13 protein. Secondary structure and disorder score was predicted by CFSSP and PONDR-VSL2 webservers. Protein stability was predicted by DynaMut webserver. B cell epitopes were predicted by IEDB DiscoTope 2.0 tools and their 3D structures were represented by discovery studio. Antigenicity and allergenicity of epitopes were predicted by Vaxijen2.0 and AllergenFPv.1.0. Physiochemical properties of epitopes were predicted by Toxinpred, HLP webserver tool. Results: Our data revealed 182 mutations in Nsp13 among Indian SARS-CoV-2 isolates, which were characterised by secondary structure and per-residue disorderness, stability and dynamicity predictions. To correlate the functional impact of these mutations, we characterised the most prominent B cell and T cell epitopes contributed by Nsp13. Our data revealed twenty-one epitopes, which exhibited antigenicity, stability and interactions with MHC class-I and class-II molecules. Subsequently, the physiochemical properties of these epitopes were analysed. Furthermore, eighteen mutations reside in these Nsp13 epitopes. Conclusions: We report appearance of eighteen mutations in the predicted twenty-one epitopes of Nsp13. Among these, at least seven epitopes closely matches with the functionally validated epitopes. Altogether, our study shows the pattern of evolution of Nsp13 epitopes and their probable implications.

Introducción y objetivo: Las vacunas se administran a nivel mundial para controlar la pandemia en curso de la enfermedad por coronavirus de 2019 (COVID-19) causada por SARS-CoV-2. A la eficacia de la vacuna contribuyen ampliamente los epítopes presentes en las proteínas virales, y su alteración puede contribuir a que las variantes emergentes se escapen de la vigilancia inmunológica del huésped. Por tanto, este estudio fue diseñado para estudiar la proteína Nsp13 de SARS-CoV-2, sus epítopes y su evolución. Métodos: Se utilizó Clustal Omega para identificar las mutaciones de la proteína Nsp13. La estructura secundaria y la tasa de desorden se predijeron mediante los servidores web CFSSP y PONDR-VSL2. La estabilidad de la proteína fue predicha mediante el servidor web DynaMut. Los epítopes de las células B fueron predichos mediante las herramientas DiscoTope 2.0 de IEDB, y sus estructuras en 3D fueron representadas mediante Discovery Studio.La antigenicidad y alergenicidad de los epítopes fueron predichas mediante Vaxijen2.0 y AlergenFPv.1.0. Las propiedades fisioquímicas de los epítopes fueron predichas mediante Toxinpred, la herramienta del servidor web HLP. Resultados: Nuestros datos revelaron 182 mutaciones en Nsp13 entre los aislados indios de SARS-CoV-2, que fueron caracterizadas mediante las predicciones de la estructura secundaria y la capacidad de desorden por residuo, la estabilidad y la dinamicidad. Para correlacionar el impacto funcional de estas mutaciones, caracterizamos los epítopes más prominentes de las células B y las células T a los que contribuyó Nsp13. Nuestros datos revelaron veintiún epítopes, que exhibieron antigenicidad, estabilidad e interacciones con las moléculas MHC de clase I y clase II. Seguidamente se analizaron las propiedades fisioquímicas de estos epítopes. Además, en estos epítopes de Nsp13 residen ocho mutaciones. Conclusiones: Reportamos el aspecto de ocho mutaciones en los veintiún epítopes de Nsp13 predichos. Entre estos, al menos siete epítopes concuerdan estrechamente con los epítopes funcionalmente validados. En su conjunto, nuestro estudio refleja el patrón evolutivo de los epítopes de Nsp13 y sus implicaciones probables.

RNA degradation eliminates developmental transcripts during murine embryonic stem cell differentiation via CAPRIN1-XRN2.

Embryonic stem cells (ESCs) are self-renewing and pluripotent. In recent years, factors that control pluripotency, mostly nuclear, have been identified. To identify non-nuclear regulators of ESCs, we screened an endogenously labeled fluorescent fusion-protein library in mouse ESCs. One of the more compelling hits was the cell-cycle-associated protein 1 (CAPRIN1). CAPRIN1 knockout had little effect in ESCs, but it significantly altered differentiation and gene expression programs. Using RIP-seq and SLAM-seq, we found that CAPRIN1 associates with, and promotes the degradation of, thousands of RNA transcripts. CAPRIN1 interactome identified XRN2 as the likely ribonuclease. Upon early ESC differentiation, XRN2 is located in the nucleus and colocalizes with CAPRIN1 in small RNA granules in a CAPRIN1-dependent manner. We propose that CAPRIN1 regulates an RNA degradation pathway operating during early ESC differentiation, thus eliminating undesired spuriously transcribed transcripts in ESCs.

Single-cell RNA-seq reveals early heterogeneity during aging in yeast.

The budding yeast Saccharomyces cerevisiae (S. cerevisiae) has relatively short lifespan and is genetically tractable, making it a widely used model organism in aging research. Here, we carried out a systematic and quantitative investigation of yeast aging with single-cell resolution through transcriptomic sequencing. We optimized a single-cell RNA sequencing (scRNA-seq) protocol to quantitatively study the whole transcriptome profiles of single yeast cells at different ages, finding increased cell-to-cell transcriptional variability during aging. The single-cell transcriptome analysis also highlighted key biological processes or cellular components, including oxidation-reduction process, oxidative stress response (OSR), translation, ribosome biogenesis and mitochondrion that underlie aging in yeast. We uncovered a molecular marker of FIT3, indicating the early heterogeneity during aging in yeast. We also analyzed the regulation of transcription factors and further characterized the distinctive temporal regulation of the OSR by YAP1 and proteasome activity by RPN4 during aging in yeast. Overall, our data profoundly reveal early heterogeneity during aging in yeast and shed light on the aging dynamics at the single cell level.

Emergence of unique SARS-CoV-2 ORF10 variants and their impact on protein structure and function.

The devastating impact of the ongoing coronavirus disease 2019 (COVID-19) on public health, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has made targeting the COVID-19 pandemic a top priority in medical research and pharmaceutical development. Surveillance of SARS-CoV-2 mutations is essential for the comprehension of SARS-CoV-2 variant diversity and their impact on virulence and pathogenicity. The SARS-CoV-2 open reading frame 10 (ORF10) protein interacts with multiple human proteins CUL2, ELOB, ELOC, MAP7D1, PPT1, RBX1, THTPA, TIMM8B, and ZYG11B expressed in lung tissue. Mutations and co-occurring mutations in the emerging SARS-CoV-2 ORF10 variants are expected to impact the severity of the virus and its associated consequences. In this article, we highlight 128 single mutations and 35 co-occurring mutations in the unique SARS-CoV-2 ORF10 variants. The possible predicted effects of these mutations and co-occurring mutations on the secondary structure of ORF10 variants and host protein interactomes are presented. The findings highlight the possible effects of mutations and co-occurring mutations on the emerging 140 ORF10 unique variants from secondary structure and intrinsic protein disorder perspectives.

Non-uniform aspects of the SARS-CoV-2 intraspecies evolution reopen question of its origin.

Several hypotheses have been presented on the origin of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) from its identification as the agent causing the current coronavirus disease 19 (COVID-19) pandemic. So far, no solid evidence has been found to support any hypothesis on the origin of this virus, and the issue continue to resurface over and over again. Here we have unfolded a pattern of distribution of several mutations in the SARS-CoV-2 proteins in 24 geo-locations across different continents. The results showed an evenly uneven distribution of the unique protein variants, distinct mutations, unique frequency of common conserved residues, and mutational residues across these 24 geo-locations. Furthermore, ample mutations were identified in the evolutionarily conserved invariant regions in the SARS-CoV-2 proteins across almost all geo-locations studied. This pattern of mutations potentially breaches the law of evolutionary conserved functional units of the beta-coronavirus genus. These mutations may lead to several novel SARS-CoV-2 variants with a high degree of transmissibility and virulence. A thorough investigation on the origin and characteristics of SARS-CoV-2 needs to be conducted in the interest of science and for the preparation of meeting the challenges of potential future pandemics.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2.

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

Identification and molecular characterization of mutations in nucleocapsid phosphoprotein of SARS-CoV-2.

SARS-CoV-2 genome encodes four structural proteins that include the spike glycoprotein, membrane protein, envelope protein and nucleocapsid phosphoprotein (N-protein). The N-protein interacts with viral genomic RNA and helps in packaging. As SARS-CoV-2 spread to almost all countries worldwide within 2-3 months, it also acquired mutations in its RNA genome. Therefore, this study was conducted with an aim to identify the variations present in N-protein of SARS-CoV-2. Here, we analysed 4,163 reported sequence of N-protein from United States of America (USA) and compared them with the first reported sequence from Wuhan, China. Our study identified 107 mutations that reside all over the N-protein. Further, we show the high rate of mutations in intrinsically disordered regions (IDRs) of N-protein. Our study show 45% residues of IDR2 harbour mutations. The RNA-binding domain (RBD) and dimerization domain of N-protein also have mutations at key residues. We further measured the effect of these mutations on N-protein stability and dynamicity and our data reveals that multiple mutations can cause considerable alterations. Altogether, our data strongly suggests that N-protein is one of the mutational hotspot proteins of SARS-CoV-2 that is changing rapidly and these mutations can potentially interferes with various aspects of N-protein functions including its interaction with RNA, oligomerization and signalling events.