SARA-Coffee web server, a tool for the computation of RNA sequence and structure multiple alignments

Unraveling RNA dynamical behavior of TPP riboswitches: a comparison between Escherichia coli and Arabidopsis thaliana

Riboswitches are RNA sensors that affect post-transcriptional processes through their ability to bind to small molecules. Thiamine pyrophosphate (TPP) riboswitch class is the most widespread riboswitch occurring in all three domains of life. Even though it controls different genes involved in the synthesis or transport of thiamine and its phosphorylated derivatives in bacteria, archaea, fungi, and plants, the TPP aptamer has a conserved structure. In this study, we aimed at understanding differences in the structural dynamics of TPP riboswitches from Escherichia coli and Arabidopsis thaliana, based on their crystallographic structures (TPPsw^ec and TPPsw^at, respectively) and dynamics in aqueous solution, both in apo and holo states. A combination of Molecular Dynamics Simulations and Network Analysis empowered to find out slight differences in the dynamical behavior of TPP riboswitches, although relevant for their dynamics in bacteria and plants species. Our results suggest that distinct interactions in the microenvironment surrounding nucleotide U36 of TPPsw^ec (and U35 in TPPsw^at) are related to different responses to TPP. The network analysis showed that minor structural differences in the aptamer enable enhanced intramolecular communication in the presence of TPP in TPPsw^ec, but not in TPPsw^at. TPP riboswitches of plants present subtler and slower regulation mechanisms than bacteria do.

Identification and characterization of two selenium-dependent glutathione peroxidase 1 isoforms from Larimichthys crocea

Glutathione peroxidases, a vital family of antioxidant enzymes in oxybiotic organisms, are involved in anti-pathogen immune response. In this study, two complete selenium-dependent glutathione peroxidase 1 cDNAs (designated as LcGPx1a and LcGPx1b) were obtained from the large yellow croaker Larimichthys crocea by rapid amplification of cDNA ends. The full-length sequence of LcGPx1a was 917 bp with a 5'-untranslated region (UTR) of 52 bp, a 3'-UTR of 289 bp, and an open reading frame of 576 bp encoding 191 amino acid (aa) polypeptides. The cDNA of LcGPx1b was composed of 884 bp with a 5'-UTR of 59 bp, a 3'-UTR of 258 bp, and an open reading frame of 567 bp encoding 188 aa polypeptides. The conserved selenocysteine insertion sequence was detected in the 3'-UTR of both isoforms, which can classify types I and II. Protein sequence analysis revealed that both isoforms included a selenocysteine encoded by an opal codon (TGA) and formed the functioning tetrad site with glutamine, tryptophan, and asparagine. Three conservative motifs, including one active site motif ("GKVVLIENVASLUGTT") and two signature site motifs ("LVILGVPCNQFGHQENC" and "V(A/S)WNFEKFLI"), were conserved both in sequence and location. Multiple alignments revealed that they exhibited a high level of identities with GPx1 from other organisms, especially in the abovementioned conserved amino acid sequence motifs. Tissue expression analysis indicated that LcGPx1a and LcGPx1b had a wide distribution in nine tissues with various abundances. The transcript level of LcGPx1a was not significantly different among the nine tissues, whereas that of LcGPx1b was higher in the kidney and head kidney than in the other tissues. After Vibrio parahaemolyticus stimulation, the expression levels of LcGPx1a and LcGPx1b were unanimously altered in the liver, spleen, kidney, and head kidney but with different magnitudes and response time. LcGPx1a and LcGPx1b showed distinct expression trends in the liver, where LcGPx1b was induced and LcGPx1a was depressed in response to pathogen infection. These results indicate that LcGPx1a and LcGPx1b display functional diversities and play crucial roles in mediating the immune response of fish.

SupeRNAlign: a new tool for flexible superposition of homologous RNA structures and inference of accurate structure-based sequence alignments

RNA has been found to play an ever-increasing role in a variety of biological processes. The function of most non-coding RNA molecules depends on their structure. Comparing and classifying macromolecular 3D structures is of crucial importance for structure-based function inference and it is used in the characterization of functional motifs and in structure prediction by comparative modeling. However, compared to the numerous methods for protein structure superposition, there are few tools dedicated to the superimposing of RNA 3D structures. Here, we present SupeRNAlign (v1.3.1), a new method for flexible superposition of RNA 3D structures, and SupeRNAlign-Coffee—a workflow that combines SupeRNAlign with T-Coffee for inferring structure-based sequence alignments. The methods have been benchmarked with eight other methods for RNA structural superposition and alignment. The benchmark included 151 structures from 32 RNA families (with a total of 1734 pairwise superpositions). The accuracy of superpositions was assessed by comparing structure-based sequence alignments to the reference alignments from the Rfam database. SupeRNAlign and SupeRNAlign-Coffee achieved significantly higher scores than most of the benchmarked methods: SupeRNAlign generated the most accurate sequence alignments among the structure superposition methods, and SupeRNAlign-Coffee performed best among the sequence alignment methods.