NOX family NADPH oxidases in mammals: Evolutionary conservation and isoform-defining sequences

NADPH oxidases are superoxide-producing enzymes that play a role in host defense, biosynthetic pathways, as well as cellular signaling. Humans have 7 NOX isoforms (NOX1-5, DUOX1,2), while mice and rats lack NOX5 and therefore have only 6 NOX isoforms. Whether all human NOX isoforms or their subunits (CYBA, NCF1, 2, 4, NOXO1, NOXA1, DUOXA1, 2) are present and conserved in other mammalian species is unknown. In this study, we have analyzed the conservation of the NOX family during mammalian evolution using an in-silico approach. Complete genomic sequences of 164 mammalian species were available. The possible absence of genes coding for NOX isoforms was investigated using the NCBI orthologs database followed by manual curation. Conservation of a given NOX isoform during mammalian evolution was evaluated by multiple alignment and identification of highly conserved sequences. There was no convincing evidence for the absence of NOX2, 3, 4, and DUOX1, 2 in all the available mammalian genome. However, NOX5 was absent in 27 of 31 rodent, in 2 of 3 lagomorph and in 2 out of 18 bat species. NOX1 was absent in all sequenced Afrotheria and Monotremata species, as well as in 3 of 18 bat species. NOXA1 was absent in all Afrotheria and in 3 out of 4 Eulipotyphla species. We also investigated amino acid sequence conservation among given NOX isoforms. Highly conserved sequences were observed for most isoforms except for NOX5. Interestingly, the highly conserved region of NOX2 sequence was relatively small (11 amino acids), as compared to NOX1, 3, 4. The highly conserved domains are different from one NOX isoform to the other, raising the possibility of distinct evolutionary conserved functional domains. Our results shed a new light on the essentiality of different NOX isoforms. We also identified isoform-defining sequences, i.e., hitherto undescribed conserved domains within specific NOX isoforms.

Mustguseal and Sister Web-Methods: A Practical Guide to Bioinformatic Analysis of Protein Superfamilies

Bioinformatic analysis of functionally diverse superfamilies can help to study the structure-function relationship in proteins, but represents a methodological challenge. The Mustguseal web-server can build large structure-guided sequence alignments of thousands of homologs that cover all currently available sequence variants within a common structural fold. The input to the method is a PDB code of the query protein, which represents the protein superfamily of interest. The collection and subsequent alignment of protein sequences and structures is fully automated and driven by the particular choice of parameters. Four integrated sister web-methods-the Zebra, pocketZebra, visualCMAT, and Yosshi-are available to further analyze the resulting superimposition and identify conserved, subfamily-specific, and co-evolving residues, as well as to classify and study disulfide bonds in protein superfamilies. The integration of these web-based bioinformatic tools provides an out-of-the-box easy-to-use solution, first of its kind, to study protein function and regulation and design improved enzyme variants for practical applications and selective ligands to modulate their functional properties. In this chapter, we provide a step-by-step protocol for a comprehensive bioinformatic analysis of a protein superfamily using a web-browser as the main tool and notes on selecting the appropriate values for the key algorithm parameters depending on your research objective. The web-servers are freely available to all users at https://biokinet.belozersky.msu.ru/m-platform with no login requirement.

Methods for making multiple alignment of genomic sequences for severe acute respiratory syndrome coronavirus 2

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in December 2019 and caused a pandemic. To monitor the global transmission pattern of SARS-CoV-2, it is required to constantly update the phylogenetic tree of genomic sequences with 29.9 kb, which may be time consuming. Phylogenetic analysis of SARS-CoV-2 may be accelerated by making a multiple alignment of nucleotide sequences using the CPA (combining pairwise alignments) method, in which a pairwise alignment is made for a reference and each of other sequences, and the pairwise alignments are combined into a multiple alignment. Here it is shown from the analysis of 3729 genomic sequences for SARS-CoV-2 and outgroup strains that the CPA method can produce a multiple alignment with an elevated or a reduced number of variable sites depending on the reference compared to the OMA (ordinary multiple alignment) method, which was considered to be the most reliable. In particular, the topology of the phylogenetic tree constructed from the multiple alignment made using the CPA method adopting the outgroup sequence as the reference was considerably different from that using the OMA method, suggesting that the outgroup sequence may not be suitable as the reference in the CPA method.

DeepShape: estimating isoform-level ribosome abundance and distribution with Ribo-seq data

Background

Ribosome profiling brings insight to the process of translation. A basic step in profile construction at transcript level is to map Ribo-seq data to transcripts, and then assign a huge number of multiple-mapped reads to similar isoforms. Existing methods either discard the multiple mapped-reads, or allocate them randomly, or assign them proportionally according to transcript abundance estimated from RNA-seq data.

Results

Here we present DeepShape, an RNA-seq free computational method to estimate ribosome abundance of isoforms, and simultaneously compute their ribosome profiles using a deep learning model. Our simulation results demonstrate that DeepShape can provide more accurate estimations on both ribosome abundance and profiles when compared to state-of-the-art methods. We applied DeepShape to a set of Ribo-seq data from PC3 human prostate cancer cells with and without PP242 treatment. In the four cell invasion/metastasis genes that are translationally regulated by PP242 treatment, different isoforms show very different characteristics of translational efficiency and regulation patterns. Transcript level ribosome distributions were analyzed by “Codon Residence Index (CRI)” proposed in this study to investigate the relative speed that a ribosome moves on a codon compared to its synonymous codons. We observe consistent CRI patterns in PC3 cells. We found that the translation of several codons could be regulated by PP242 treatment.

Conclusion

In summary, we demonstrate that DeepShape can serve as a powerful tool for Ribo-seq data analysis.

Advanced Genetic Methodologies in Tracking Evolution and Spread of Chikungunya Virus

Recent advances in genetic methodologies have substantially expanded our ability to track evolution and spatio-temporal distribution of rapidly evolving pathogens. The information gathered from such analyses can be used to decipher host adaptations that shape disease epidemiology. In this chapter, we demonstrate the utilization of freely available resources to track the evolution and spread of Chikungunya virus.

Comparative study of structural models of Leishmania donovani and human GDP-mannose pyrophosphorylases

Leishmania is the parasite responsible for the neglected disease leishmaniasis. Its virulence and survival require biosynthesis of glycoconjugates, whose guanosine diphospho-d-mannose pyrophosphorylase (GDP-MP) is a key player. However, experimentally resolved structures of this enzyme are still lacking. We herein propose structural models of the GDP-MP from human and Leishmania donovani. Based on a multiple sequences alignment, the models were built with MODELLER and then carefully refined with all atom molecular dynamics simulations in explicit solvent. Their quality was evaluated against several standard criteria, including their ability to bind GDP-mannose assessed by redocking calculations. Special attention was given in this study to interactions of the catalytic site residues with the enzyme substrate and competitive inhibitors, opening the perspective of medicinal chemistry developments.

Progressive alignment of genomic signals by multiple dynamic time warping

This paper presents the utilization of progressive alignment principle for positional adjustment of a set of genomic signals with different lengths. The new method of multiple alignment of signals based on dynamic time warping is tested for the purpose of evaluating the similarity of different length genes in phylogenetic studies. Two sets of phylogenetic markers were used to demonstrate the effectiveness of the evaluation of intraspecies and interspecies genetic variability. The part of the proposed method is modification of pairwise alignment of two signals by dynamic time warping with using correlation in a sliding window. The correlation based dynamic time warping allows more accurate alignment dependent on local homologies in sequences without the need of scoring matrix or evolutionary models, because mutual similarities of residues are included in the numerical code of signals.