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.

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.

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.

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.

Implications derived from S-protein variants of SARS-CoV-2 from six continents

The spike (S) protein is a critical determinant of the infectivity and antigenicity of SARS-CoV-2. Several mutations in the S protein of SARS-CoV-2 have already been detected, and their effect in immune system evasion and enhanced transmission as a cause of increased morbidity and mortality are being investigated. From pathogenic and epidemiological perspectives, S proteins are of prime interest to researchers. This study focused on the unique variants of S proteins from six continents: Asia, Africa, Europe, Oceania, South America, and North America. In comparison to the other five continents, Africa had the highest percentage of unique S proteins (29.1%). The phylogenetic relationship implies that unique S proteins from North America are significantly different from those of the other five continents. They are most likely to spread to the other geographic locations through international travel or naturally by emerging mutations. It is suggested that restriction of international travel should be considered, and massive vaccination as an utmost measure to combat the spread of the COVID-19 pandemic. It is also further suggested that the efficacy of existing vaccines and future vaccine development must be reviewed with careful scrutiny, and if needed, further re-engineered based on requirements dictated by new emerging S protein variants.

Autoimmunity roots of the thrombotic events after COVID-19 vaccination

Although vaccination represents the most promising way to stop or contain the coronavirus disease 2019 (COVID-19) pandemic and safety and effectiveness of available vaccines were proven, a small number of individuals who received anti-SARS-CoV-2 vaccines developed a prothrombotic syndrome. Vaccine-induced immune thrombotic thrombocytopenia (VITT) can be triggered by the adenoviral vector-based vaccine, whereas lipid nanoparticle-mRNA-based vaccines can induce rare cases of deep vein thrombosis (DVT). Although the main pathogenic mechanisms behind this rare phenomenon have not yet been identified, both host and vaccine factors might be involved, with pathology at least in part being related to the vaccine-triggered autoimmune reaction. In this review, we are considering some aspects related to pathogenesis, major risk factors, as well as peculiarities of diagnosis and treatment of this rare condition.

COVID-19 Vaccines and Thrombosis-Roadblock or Dead-End Street?

Two adenovirus-based vaccines, ChAdOx1 nCoV-19 and Ad26.COV2.S, and two mRNA-based vaccines, BNT162b2 and mRNA.1273, have been approved by the European Medicines Agency (EMA), and are invaluable in preventing and reducing the incidence of coronavirus disease-2019 (COVID-19). Recent reports have pointed to thrombosis with associated thrombocytopenia as an adverse effect occurring at a low frequency in some individuals after vaccination. The causes of such events may be related to SARS-CoV-2 spike protein interactions with different C-type lectin receptors, heparan sulfate proteoglycans (HSPGs) and the CD147 receptor, or to different soluble splice variants of the spike protein, adenovirus vector interactions with the CD46 receptor or platelet factor 4 antibodies. Similar findings have been reported for several viral diseases after vaccine administration. In addition, immunological mechanisms elicited by viral vectors related to cellular delivery could play a relevant role in individuals with certain genetic backgrounds. Although rare, the potential COVID-19 vaccine-induced immune thrombotic thrombocytopenia (VITT) requires immediate validation, especially in risk groups, such as the elderly, chronic smokers, and individuals with pre-existing incidences of thrombocytopenia; and if necessary, a reformulation of existing vaccines.

Carbon-Based Nanomaterials: Promising Antiviral Agents to Combat COVID-19 in the Microbial-Resistant Era

Therapeutic options for the highly pathogenic human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causing the current pandemic coronavirus disease (COVID-19) are urgently needed. COVID-19 is associated with viral pneumonia and acute respiratory distress syndrome causing significant morbidity and mortality. The proposed treatments for COVID-19 have shown little or no effect in the clinic so far. Additionally, bacterial and fungal pathogens contribute to the SARS-CoV-2-mediated pneumonia disease complex. The antibiotic resistance in pneumonia treatment is increasing at an alarming rate. Therefore, carbon-based nanomaterials (CBNs), such as fullerene, carbon dots, graphene, and their derivatives constitute a promising alternative due to their wide-spectrum antimicrobial activity, biocompatibility, biodegradability, and capacity to induce tissue regeneration. Furthermore, the antimicrobial mode of action is mainly physical (e.g., membrane distortion), characterized by a low risk of antimicrobial resistance. In this Review, we evaluated the literature on the antiviral activity and broad-spectrum antimicrobial properties of CBNs. CBNs had antiviral activity against 13 enveloped positive-sense single-stranded RNA viruses, including SARS-CoV-2. CBNs with low or no toxicity to humans are promising therapeutics against the COVID-19 pneumonia complex with other viruses, bacteria, and fungi, including those that are multidrug-resistant.

Pathogenic perspective of missense mutations of ORF3a protein of SARS-CoV-2

One of the most important proteins for COVID-19 pathogenesis in SARS-CoV-2 is the ORF3a which is the largest accessory protein among others coded by the SARS-CoV-2 genome. The major roles of the protein include virulence, infectivity, ion channel activity, morphogenesis, and virus release. The coronavirus, SARS-CoV-2 is mutating rapidly, therefore, critical study of mutations in ORF3a is certainly important from the pathogenic perspective. Here, a sum of 175 non-synonymous mutations in the ORF3a of SARS-CoV-2 were identified from 7194 complete genomes of SARS-CoV-2 available from NCBI database. Effects of these mutations on structural stability, and functions of ORF3a were also studied. Broadly, three different classes of mutations, such as neutral, disease, and mixed (neutral and disease) types of mutations were observed. Consecutive phenomena of mutations in ORF3a protein were studied based on the timeline of detection of the mutations. Considering the amino acid compositions of the ORF3a protein, twenty clusters were detected using the K-means clustering method. The present findings on 175 novel mutations of ORF3a proteins will extend our knowledge on ORF3a, a vital accessory protein in SARS-CoV-2, to enlighten the pathogenicity of this life-threatening virus.

Notable sequence homology of the ORF10 protein introspects the architecture of SARS-CoV-2

The current Coronavirus Disease 19 (COVID-19) pandemic, caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) shows similar pathology to MERS and SARS-CoV, with a current estimated fatality rate of 1.4%. Open reading frame 10 (ORF10) is a unique SARS-CoV-2 accessory protein, which contains eleven cytotoxic T lymphocyte (CTL) epitopes each of nine amino acids in length. Twenty-two unique SARS-CoV-2 ORF10 variants have been identified based on missense mutations found in sequence databases. Some of these mutations are predicted to decrease the stability of ORF10 in silico physicochemical and structural comparative analyses were carried out on SARS-CoV-2 and Pangolin-CoV ORF10 proteins, which share 97.37% amino acid (aa) homology. Though there is a high degree of ORF10 protein similarity of SARS-CoV-2 and Pangolin-CoV, there are differences of these two ORF10 proteins related to their sub-structure (loop/coil region), solubility, antigenicity and shift from strand to coil at aa position 26 (tyrosine). SARS-CoV-2 ORF10, which is apparently expressed in vivo since reactive T cell clones are found in convalescent patients should be monitored for changes which could correlate with the pathogenesis of COVID-19.

A unique view of SARS-CoV-2 through the lens of ORF8 protein

Immune evasion is one of the unique characteristics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) attributed to its ORF8 protein. This protein modulates the adaptive host immunity through down-regulation of MHC-1 (Major Histocompatibility Complex) molecules and innate immune responses by surpassing the host's interferon-mediated antiviral response. To understand the host's immune perspective in reference to the ORF8 protein, a comprehensive study of the ORF8 protein and mutations possessed by it have been performed. Chemical and structural properties of ORF8 proteins from different hosts, such as human, bat, and pangolin, suggest that the ORF8 of SARS-CoV-2 is much closer to ORF8 of Bat RaTG13-CoV than to that of Pangolin-CoV. Eighty-seven mutations across unique variants of ORF8 in SARS-CoV-2 can be grouped into four classes based on their predicted effects (Hussain et al., 2021) [1]. Based on the geo-locations and timescale of sample collection, a possible flow of mutations was built. Furthermore, conclusive flows of amalgamation of mutations were found upon sequence similarity analyses and consideration of the amino acid conservation phylogenies. Therefore, this study seeks to highlight the uniqueness of the rapidly evolving SARS-CoV-2 through the ORF8.

Urgent Need for Field Surveys of Coronaviruses in Southeast Asia to Understand the SARS-CoV-2 Phylogeny and Risk Assessment for Future Outbreaks

Phylogenetic analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is focused on a single isolate of bat coronaviruses (bat CoVs) which does not adequately represent genetically related coronaviruses (CoVs) [...].

Rare mutations in the accessory proteins ORF6, ORF7b, and ORF10 of the SARS-CoV-2 genomes

A total number of 3080 SARS-CoV-2 genomes from all continents are considered from the NCBI database. Every accessory protein ORF6, ORF7b, and ORF10 of SARS-CoV-2 possess a single missense mutation in less than 1.5% of the 3080 genomes. It has now been observed that different non-synonymous mutations occurred in these three accessory proteins. Most of these rare mutations are changing the amino acids such as hydrophilic to hydrophobic, acidic or basic to hydrophobic, and vice versa etc. So these highly conserved proteins might play an essential role in virus pathogenicity. This study opens a question whether it carries some messages about the virus rapid replications, and virulence.

Clade GR and clade GH isolates of SARS-CoV-2 in Asia show highest amount of SNPs

Clades are monophyletic groups composed of a common ancestor and all its lineal descendants. As the propensity of virulence of a disease depends upon the type of clade the virus belongs to and it causes different fatality rates of disease in different countries, so the clade-wise analysis of SARS-CoV-2 isolates collected from different countries can illuminate the actual evolutionary relationships between them. In this study, 1566 SARS-CoV-2 genome sequences across ten Asian countries are collected, clustered, and characterized based on the clade they belong to. The isolates are compared to the Wuhan reference sequence” hCoV-19/Wuhan/WIV04/19″ to identify the mutations that occurred at different protein regions. Structural changes in amino acids due to mutations lead to functional instability of the proteins. Detailed clade-wise functional assessments are carried out to quantify the stability and vulnerability of the mutations occurring in SARS-CoV-2 genomes which can shade light on personalized prevention and treatment of the disease and encourage towards the invention of clade-specific vaccines.

The structural basis of accelerated host cell entry by SARS-CoV-2†

Severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) is the causative agent of the pandemic coronavirus disease 2019 (COVID‐19) that exhibits an overwhelming contagious capacity over other human coronaviruses (HCoVs). This structural snapshot describes the structural bases underlying the pandemic capacity of SARS‐CoV‐2 and explains its fast motion over respiratory epithelia that allow its rapid cellular entry. Based on notable viral spike (S) protein features, we propose that the flat sialic acid‐binding domain at the N‐terminal domain (NTD) of the S1 subunit leads to more effective first contact and interaction with the sialic acid layer over the epithelium, and this, in turn, allows faster viral ‘surfing’ of the epithelium and receptor scanning by SARS‐CoV‐2. Angiotensin‐converting enzyme 2 (ACE‐2) protein on the epithelial surface is the primary entry receptor for SARS‐CoV‐2, and protein–protein interaction assays demonstrate high‐affinity binding of the spike protein (S protein) to ACE‐2. To date, no high‐frequency mutations were detected at the C‐terminal domain of the S1 subunit in the S protein, where the receptor‐binding domain (RBD) is located. Tight binding to ACE‐2 by a conserved viral RBD suggests the ACE2‐RBD interaction is likely optimal. Moreover, the viral S subunit contains a cleavage site for furin and other proteases, which accelerates cell entry by SARS‐CoV‐2. The model proposed here describes a structural basis for the accelerated host cell entry by SARS‐CoV‐2 relative to other HCoVs and also discusses emerging hypotheses that are likely to contribute to the development of antiviral strategies to combat the pandemic capacity of 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.

Variability of Accessory Proteins Rules the SARS-CoV-2 Pathogenicity.. [DOI: 10.1101/2020.11.06.372227]

The coronavirus disease 2019 (COVID-19) is caused by the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) which is pandemic with an estimated fatality rate less than 1% is ongoing. SARS-CoV-2 accessory proteins ORF3a, ORF6, ORF7a, ORF7b, ORF8, and ORF10 with putative functions to manipulate host immune mechanisms such as interferons, immune signaling receptor NLRP3 (NOD-, LRR-, and pyrin domain-containing 3) inflammasome, inflammatory cytokines such as interleukin 1β(IL-1β) are critical in COVID-19 pathology. Outspread variations of each of the six accessory proteins of all complete proteomes (available as of October 26, 2020, in the National Center for Biotechnology Information depository) of SARS-CoV-2, were observed across six continents. Across all continents, the decreasing order of percentage of unique variations in the accessory proteins was found to be ORF3a>ORF8>ORF7a>ORF6>ORF10>ORF7b. The highest and lowest unique variations of ORF3a were observed in South America and Oceania, respectively. This finding suggests that the wide variations of accessory proteins seem to govern the pathogenicity of SARS-CoV-2, and consequently, certain propositions and recommendations can be made in the public interest.

SARS-CoV2 envelope protein: non-synonymous mutations and its consequences

In the NCBI database, as on June 6, 2020, total number of available complete genome sequences of SARS-CoV2 across the world is 3617. The envelope (E) protein of SARS-CoV2 possesses several non-synonymous mutations over the transmembrane and C-terminus domains in 15 (0.414%) genomes among 3617 SARS-CoV2 genomes, analyzed. More precisely, 10(0.386%) out of 2588 genomes from the USA, 3(0.806%) from Asia, 1 (0.348%) from Europe and 1 (0.274%) from Oceania contained the missense mutations over the E-protein of SARS-CoV2 genomes. The C-terminus motif DLLV has been to DFLV and YLLV in the proteins from QJR88103 (Australia: Victoria) and QKI36831 (China: Guangzhou) respectively, which might affect the binding of this motif with the host protein PALS1.

Missense mutations in SARS-CoV2 genomes from Indian patients

The genetic diversity of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV2) in several countries sums up to worldwide genetic diversity. In this present study, variations in terms of missense mutations among the SARS-CoV2 genomes from 128 Indian patients, as of May 2020, are accounted and thereby some key findings with some hypotheses were made. These mutations across various genes of these genomes show wide genetic variations in sequence and rapid evolution of SARS-CoV2 virus. The presence of unique mutations in the studied SARS-CoV2 genomes may led to their attenuation. Few Genes such as ORF6, ORF10 are free from any mutations in the Indian context of 339 genomes observed as of 14th July 2020. Further, E protein contains only one mutation. This may suggest that designing a therapeutic approach against ORF6, ORF10 and E genes may have a beneficial effect in controlling COVID-19 pandemic especially in India.

Molecular phylogeny and missense mutations at envelope proteins across coronaviruses

Envelope (E) protein is one of the structural viroporins (76-109 amino acids long) present in the coronavirus. Sixteen sequentially different E-proteins were observed from a total of 4917 available complete SARS-CoV-2 genomes as on 18th June 2020 in the NCBI database. The missense mutations over the envelope protein across various coronaviruses of the β-genus were analyzed to know the immediate parental origin of the envelope protein of SARS-CoV-2. The evolutionary origin is also endorsed by the phylogenetic analysis of the envelope proteins comparing sequence homology as well as amino acid conservations.

Viewpoint: Origin of SARS-CoV-2

The origin of the severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) virus causing the COVID-19 pandemic has not yet been fully determined. Despite the consensus about the SARS-CoV-2 origin from bat CoV RaTG13, discrepancy to host tropism to other human Coronaviruses exist. SARS-CoV-2 also possesses some differences in its S protein receptor-binding domain, glycan-binding N-terminal domain and the surface of the sialic acid-binding domain. Despite similarities based on cryo-EM and biochemical studies, the SARS-CoV-2 shows higher stability and binding affinity to the ACE2 receptor. The SARS-CoV-2 does not appear to present a mutational “hot spot” as only the D614G mutation has been identified from clinical isolates. As laboratory manipulation is highly unlikely for the origin of SARS-CoV-2, the current possibilities comprise either natural selection in animal host before zoonotic transfer or natural selection in humans following zoonotic transfer. In the former case, despite SARS-CoV-2 and bat RaTG13 showing 96% identity some pangolin Coronaviruses exhibit very high similarity to particularly the receptor-binding domain of SARS-CoV-2. In the latter case, it can be hypothesized that the SARS-CoV-2 genome has adapted during human-to-human transmission and based on available data, the isolated SARS-CoV-2 genomes derive from a common origin. Before the origin of SARS-CoV-2 can be confirmed additional research is required

Notable sequence homology of the ORF10 protein introspects the architecture of SARS-COV-2.. [DOI: 10.1101/2020.09.06.284976]

The global public health is endangered due to COVID-19 pandemic, which is caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2). Despite having similar pathology to MERS and SARS-CoV, the infection fatality rate of SARS-CoV-2 is likely lower than 1%. SARS-CoV-2 has been reported to be uniquely characterized by the accessory protein ORF10, which contains eleven cytotoxic T lymphocyte (CTL) epitopes of nine amino acids length each, across various human leukocyte antigen (HLA) subtypes. In this study, all missense mutations found in sequence databases were examined across twnety-two unique SARS-CoV-2 ORF10 variants that could possibly alter viral pathogenicity. Some of these mutations decrease the stability of ORF10, e.g. I4L and V6I were found in the MoRF region of ORF10 which may also possibly contribute to Intrinsic protein disorder. Furthermore, a physicochemical and structural comparative analysis was carried out on SARS-CoV-2 and Pangolin-CoV ORF10 proteins, which share 97.37% amino acid homology. The high degree of physicochemical and structural similarity of ORF10 proteins of SARS-CoV-2 and Pangolin-CoV open questions about the architecture of SARS-CoV-2 due to the disagreement of these two ORF10 proteins over their sub-structure (loop/coil region), solubility, antigenicity and change from the strand to coil at amino acid position 26, where tyrosine is present. Altogether, SARS-CoV-2 ORF10 is a promising pharmaceutical target and a protein which should be monitored for changes which correlate to change pathogenesis and clinical course of COVID-19 infection.

Molecular conservation and differential mutation on ORF3a gene in Indian SARS-CoV2 genomes

A global emergency due to the COVID-19 pandemic demands various studies related to genes and genomes of the SARS-CoV2. Among other important proteins, the role of accessory proteins are of immense importance in replication, regulation of infections of the coronavirus in the hosts. The largest accessory protein in the SARS-CoV2 genome is ORF3a which modulates the host response to the virus infection and consequently it plays an important role in pathogenesis. In this study, an attempt is made to decipher the conservation of nucleotides, dimers, codons and amino acids in the ORF3a genes across thirty-two genomes of Indian patients. ORF3a gene possesses single and double point mutations in Indian SARS-CoV2 genomes suggesting the change of SARS-CoV2's virulence property in Indian patients. We find that the parental origin of the ORF3a gene over the genomes of SARS-CoV2 and Pangolin-CoV is same from the phylogenetic analysis based on conservation of nucleotides and so on. This study highlights the accumulation of mutation on ORF3a in Indian SARS-CoV2 genomes which may provide the designing therapeutic approach against SARS-CoV2.

Analysis of Purines and Pyrimidines distribution over miRNAs of Human, Gorilla, Chimpanzee, Mouse and Rat

Meaningful words in English need vowels to break up the sounds that consonants make. The Nature has encoded her messages in RNA molecules using only four alphabets A, U, C and G in which the nine member double-ring bases (adenine (A) and Guanine (G)) are purines, while the six member single-ring bases (cytosine (C) and uracil (U)) are pyrimidines. Four bases A, U, C and G of RNA sequences are divided into three kinds of classifications according to their chemical properties. One of the three classifications, the purine-pyrimidine class is important. In understanding the distribution (organization) of purines and pyrimidines over some of the non-coding RNAs, all miRNAs of the three species (human, gorilla and chimpanzee) of Hominidae family and two species (mouse and rat) from of Muridae family are considered. The distribution of purines and pyrimidines over miRNA shows deviation from randomness. Based on the quantitative metrics (fractal dimension, Hurst exponent, Hamming distance, distance pattern of purine-pyrimidine, density distribution of purine-pyrimidine and Shannon entropy) five different clusters have been made for the five species. We have observed some couple of results including the closeness of different clusters among the five species.

An evolutionary analysis identifies a conserved pentapeptide stretch containing the two essential lysine residues for rice L-myo-inositol 1-phosphate synthase catalytic activity

A molecular evolutionary analysis of a well conserved protein helps to determine the essential amino acids in the core catalytic region. Based on the chemical properties of amino acid residues, phylogenetic analysis of a total of 172 homologous sequences of a highly conserved enzyme, L-myo-inositol 1-phosphate synthase or MIPS from evolutionarily diverse organisms was performed. This study revealed the presence of six phylogenetically conserved blocks, out of which four embrace the catalytic core of the functional protein. Further, specific amino acid modifications targeting the lysine residues, known to be important for MIPS catalysis, were performed at the catalytic site of a MIPS from monocotyledonous model plant, Oryza sativa (OsMIPS1). Following this study, OsMIPS mutants with deletion or replacement of lysine residues in the conserved blocks were made. Based on the enzyme kinetics performed on the deletion/replacement mutants, phylogenetic and structural comparison with the already established crystal structures from non-plant sources, an evolutionarily conserved peptide stretch was identified at the active pocket which contains the two most important lysine residues essential for catalytic activity.

Chemical property based sequence characterization of PpcA and its homolog proteins PpcB-E: A mathematical approach

Periplasmic c7 type cytochrome A (PpcA) protein is determined in Geobacter sulfurreducens along with its other four homologs (PpcB-E). From the crystal structure viewpoint the observation emerges that PpcA protein can bind with Deoxycholate (DXCA), while its other homologs do not. But it is yet to be established with certainty the reason behind this from primary protein sequence information. This study is primarily based on primary protein sequence analysis through the chemical basis of embedded amino acids. Firstly, we look for the chemical group specific score of amino acids. Along with this, we have developed a new methodology for the phylogenetic analysis based on chemical group dissimilarities of amino acids. This new methodology is applied to the cytochrome c7 family members and pinpoint how a particular sequence is differing with others. Secondly, we build a graph theoretic model on using amino acid sequences which is also applied to the cytochrome c7 family members and some unique characteristics and their domains are highlighted. Thirdly, we search for unique patterns as subsequences which are common among the group or specific individual member. In all the cases, we are able to show some distinct features of PpcA that emerges PpcA as an outstanding protein compared to its other homologs, resulting towards its binding with deoxycholate. Similarly, some notable features for the structurally dissimilar protein PpcD compared to the other homologs are also brought out. Further, the five members of cytochrome family being homolog proteins, they must have some common significant features which are also enumerated in this study.

Mathematical Characterization of Protein Sequences Using Patterns as Chemical Group Combinations of Amino Acids

Comparison of amino acid sequence similarity is the fundamental concept behind the protein phylogenetic tree formation. By virtue of this method, we can explain the evolutionary relationships, but further explanations are not possible unless sequences are studied through the chemical nature of individual amino acids. Here we develop a new methodology to characterize the protein sequences on the basis of the chemical nature of the amino acids. We design various algorithms for studying the variation of chemical group transitions and various chemical group combinations as patterns in the protein sequences. The amino acid sequence of conventional myosin II head domain of 14 family members are taken to illustrate this new approach. We find two blocks of maximum length 6 aa as 'FPKATD' and 'Y/FTNEKL' without repeating the same chemical nature and one block of maximum length 20 aa with the repetition of chemical nature which are common among all 14 members. We also check commonality with another motor protein sub-family kinesin, KIF1A. Based on our analysis we find a common block of length 8 aa both in myosin II and KIF1A. This motif is located in the neck linker region which could be responsible for the generation of mechanical force, enabling us to find the unique blocks which remain chemically conserved across the family. We also validate our methodology with different protein families such as MYOI, Myosin light chain kinase (MLCK) and Rho-associated protein kinase (ROCK), Na+/K+-ATPase and Ca2+-ATPase. Altogether, our studies provide a new methodology for investigating the conserved amino acids' pattern in different proteins.

Biodegradation of sulfosulphuron in agricultural soil by Trichoderma sp

Sulfosulphuron-degrading fungus was isolated by enrichment technique from the sulfosulphuron-contaminated soil of wheat rhizosphere. To assess the biodegradation potential of isolated Trichoderma sp., minimal potato dextrose agar broth with different levels of sulfosulphuron (up to 2 g l(-1) ) was evaluated in the growth and biotransformation experiments. ESI LC-MS/MS analysis revealed the presence of degradation products 2-amino-4,6-dimethoxypyrimidine (I) and 2-ethylsulfonyl imidazo{1,2-a} pyridine-3-sulfonamide-2-ethylsulfonyl imidazo{1,2-a} pyridine-3-sulfonamide (II) indicating the cleavage of the urea bridge and the presence of the by-product N-(4,6-dimethoxypyrimidin-2-yl)urea (III) indicating the degradation of sulfonylamide linkage. Two other metabolites, N-(4,6-dimethoxypyrimidin-2-yl)-N'-hydroxyurea (IV) and N, N'-bis(4,6-dimethoxypyrimidin-2-yl)urea (V), were also identified. From the previous reports, it was found that the degradation of sulfonyl urea herbicides took place through the chemical degradation of the sulfonylurea bridge followed by microbial degradation. During this investigation, Trichoderma sp. grew well with and degraded sulfosulphuron via both the decarboxylation on the sulphonyl urea bridge and the hydrolytic cleavage of the sulfonylamide linkage as demonstrated by the formation of metabolites. Trichoderma is nonphytopathogenic in nature, and some species of it restrict the growth of soil-dwelling phytopathogens. Therefore, it is a promising candidate for the decontamination of soil from sulfosulphuron residues.

SIGNIFICANCE AND IMPACT OF THE STUDY

The degradation of sulfosulphuron by any individual fungus is being reported for the first time. Trichoderma sp. isolated from wheat-rhizospheric soil could survive in minimal broth rich in sulfosulphuron. Previous reports have described the complete degradation of any sulfonyl urea herbicides by micro-organisms only after the pH-dependent chemical hydrolysis of the sulfonyl urea bridge of the herbicide. This study demonstrates the novel result that the Trichoderma sp. utilized the sulfosulphuron as a sole carbon source and degraded it by cleaving sulfonyl urea bridge and sulfonylamide linkage. Thus, the application of Trichoderma sp., which is nonphytopathogenic, has the potential to decontaminate agricultural soil from sulfosulphuron load.

Underlying mathematics in diversification of human olfactory receptors in different loci

As per conservative estimate, approximately 51-105 Olfactory Receptors (ORs) loci are present in human genome occurring in clusters. These clusters are apparently unevenly spread as mosaics over 21 pairs of human chromosomes. Olfactory Receptor (OR) gene families which are thought to have expanded for the need to provide recognition capability for a huge number of pure and complex odorants, form the largest known multigene family in the human genome. Recent studies have shown that 388 full length and 414 OR pseudo-genes are present in these OR genomic clusters. In this paper, the authors report a classification method for all human ORs based on their sequential quantitative information like presence of poly strings of nucleotides bases, long range correlation and so on. An L-System generated sequence has been taken as an input into a star-model of specific subfamily members and resultant sequence has been mapped to a specific OR based on the classification scheme using fractal parameters like Hurst exponent and fractal dimensions.

Understanding Genomic Evolution of Olfactory Receptors through Fractal and Mathematical Morphology

Fractals and Mathematical Morphology are immensely used to study many problems in different branches of science and technology including the domain of Biology. There are many more unrevealed facts and figures of genes and genome in Computational Biology. In this paper, our objective is to explore how the evolutionary network is associated among Human, Chimpanzee and Mouse with regards to their genomic information. We are about to explore their genomic evolution through the quantitative measures of fractals and morphology. We have considered olfactory receptors for our case study. These olfactory receptors do function in different species with subtle differences in the structures of DNA sequences. Those subtle differences can be exposed through intricate details of Fractals and Mathematical Morphology.

Designing exons for human olfactory receptor gene subfamilies using a mathematical paradigm

Ligands for only two human olfactory receptors are known. One of them, OR1D2, binds to Bourgeonal, a volatile chemical constituent of the fragrance of the mythical flower, Lily of the valley or Our Lady's tears, Convallaria majalis (also the national flower of Finland). OR1D2, OR1D4 and OR1D5 are three full-length olfactory receptors present in an olfactory locus in the human genome. These receptors are more than 80% identical in DNA sequences and have 108 base pair mismatches among them. Apparently, these mismatch positions show no striking pattern using computer pattern recognition tools. In an attempt to find a mathematical rule in those mismatches, we find that an L-system generated sequence can be inserted into the OR1D2 subfamily-specific star model and novel full-length olfactory receptors can be generated. This remarkable mathematical principle could be utilized for making new subfamily olfactory receptor members from any olfactory receptor subfamily. The aroma and electronic nose industry might utilize this rule in future.

Underlying Mathematics in Diversification of Human Olfactory Receptors in Different Loci

As per conservative estimate, approximately (51-105) Olfactory Receptors (ORs) loci are present in human genome occurring in clusters. These clusters are apparently unevenly spread as mosaics over 21 pair of human chromosomes. Olfactory Receptor (OR) gene families which are thought to have expanded for the need to provide recognition capability for huge number of pure and complex odorants. ORs form the largest known multi-gene family in the human genome. Recent studies have shown that 388 full length and 414 OR pseudo-genes are present in these OR genomic clusters. In this paper, the authors report a classification method for all human ORs based on their sequential quantitative information like presence of poly strings of nucleotides bases, long range correlation and so on. An L-System generated sequence has been taken as an input into a star-model of specific subfamily members and resultant sequence has been mapped to a specific OR based on the classification scheme using fractal parameters like Hurst exponent and fractal dimensions.

Simultaneous determination of chlorimuron-ethyl and its major metabolites in soil and water by reversed phase liquid chromatography

A rapid reversed-phase HPLC method was developed for simultaneously determination of chlorimuron-ethyl and its metabolites such as ethyl-2-aminosulfonyl benzoate, 4-methoxy-6-chloro-2-amino pyrimidine and saccharin in soil and water systems by using RP18 column and methanol-water (70:30 v/v) solvent system with detection at 230 nm. The recoveries of chlorimuron-ethyl and its metabolites ranged from 80 to 95%.