Consensus Secondary Structure Research Articles

Non-Coding RNA Analysis Using the Rfam Database.

Rfam is a database of non-coding RNA families in which each family is represented by a multiple sequence alignment, a consensus secondary structure, and a covariance model. Using a combination of manual and literature-based curation and a custom software pipeline, Rfam converts descriptions of RNA families found in the scientific literature into computational models that can be used to annotate RNAs belonging to those families in any DNA or RNA sequence. Valuable research outputs that are often locked up in figures and supplementary information files are encapsulated in Rfam entries and made accessible through the Rfam Web site. The data produced by Rfam have a broad application, from genome annotation to providing training sets for algorithm development. This article gives an overview of how to search and navigate the Rfam Web site, and how to annotate sequences with RNA families. The Rfam database is freely available at http://rfam.org. © 2018 by John Wiley & Sons, Inc.

Rfam 13.0: shifting to a genome-centric resource for non-coding RNA families.

The Rfam database is a collection of RNA families in which each family is represented by a multiple sequence alignment, a consensus secondary structure, and a covariance model. In this paper we introduce Rfam release 13.0, which switches to a new genome-centric approach that annotates a non-redundant set of reference genomes with RNA families. We describe new web interface features including faceted text search and R-scape secondary structure visualizations. We discuss a new literature curation workflow and a pipeline for building families based on RNAcentral. There are 236 new families in release 13.0, bringing the total number of families to 2687. The Rfam website is http://rfam.org.

Estrogen Receptor α- and β-Interacting Proteins Contain Consensus Secondary Structures: An Insilico Study

Background: Estrogen receptor (ER)α and ERβ are ligand-activated transcription factors that regulate gene expression by binding to estrogen-responsive elements and interacting with several coregulators through protein-protein interactions. Usually, these coregulators bind to the various conserved and functional domains of the receptor through a consensus LXXLL sequence, although variations can be found. The interaction of receptor domains and the consensus motif can be a possible target for nuclear receptor (NR) pharmacology, since modifications in these are responsible for possible pathogenesis of various diseases. Purpose: The present study focuses on the secondary structure and conserved domains of the ERα and ERβ interacting proteins, using bioinformatics tools and their relation to the function of the coregulators. Methods: Bioinformatics-based prediction tools like STRING, PSIPRED, PROTPARAM and Conserved Domain Database (CDD) were used. The prediction tools utilized in this study basically determines the characteristics of a possible coregulator by using an already existing protein as a template and determines the presence of any conserved consensus sequence. Coregulators have been enlisted with the help of NCBI, STRING and iHOP. The secondary structures were analyzed using PSIPRED and conserved domains were determined using CDD. Results: The analysis of the structure has shown the presence of conserved domains and homology between the various coregulators. Each interacting protein contains conserved domains like the nuclear coactivators’ domain, the helix-loop-helix domain and the SRC domain. Conclusion: Such studies give the characteristic features of ERα and ERβ interacting proteins and maybe useful to determine their family and uses in NR pharmacology in health and diseases.

Molecular phylogeny of 21 tropical bamboo species reconstructed by integrating non-coding internal transcribed spacer (ITS1 and 2) sequences and their consensus secondary structure.

The unavailability of the reproductive structure and unpredictability of vegetative characters for the identification and phylogenetic study of bamboo prompted the application of molecular techniques for greater resolution and consensus. We first employed internal transcribed spacer (ITS1, 5.8S rRNA and ITS2) sequences to construct the phylogenetic tree of 21 tropical bamboo species. While the sequence alone could grossly reconstruct the traditional phylogeny amongst the 21-tropical species studied, some anomalies were encountered that prompted a further refinement of the phylogenetic analyses. Therefore, we integrated the secondary structure of the ITS sequences to derive individual sequence-structure matrix to gain more resolution on the phylogenetic reconstruction. The results showed that ITS sequence-structure is the reliable alternative to the conventional phenotypic method for the identification of bamboo species. The best-fit topology obtained by the sequence-structure based phylogeny over the sole sequence based one underscores closer clustering of all the studied Bambusa species (Sub-tribe Bambusinae), while Melocannabaccifera, which belongs to Sub-Tribe Melocanneae, disjointedly clustered as an out-group within the consensus phylogenetic tree. In this study, we demonstrated the dependability of the combined (ITS sequence+structure-based) approach over the only sequence-based analysis for phylogenetic relationship assessment of bamboo.

Reconstruction of ancestral RNA sequences under multiple structural constraints.

BackgroundSecondary structures form the scaffold of multiple sequence alignment of non-coding RNA (ncRNA) families. An accurate reconstruction of ancestral ncRNAs must use this structural signal. However, the inference of ancestors of a single ncRNA family with a single consensus structure may bias the results towards sequences with high affinity to this structure, which are far from the true ancestors.MethodsIn this paper, we introduce achARNement, a maximum parsimony approach that, given two alignments of homologous ncRNA families with consensus secondary structures and a phylogenetic tree, simultaneously calculates ancestral RNA sequences for these two families.ResultsWe test our methodology on simulated data sets, and show that achARNement outperforms classical maximum parsimony approaches in terms of accuracy, but also reduces by several orders of magnitude the number of candidate sequences. To conclude this study, we apply our algorithms on the Glm clan and the FinP-traJ clan from the Rfam database.ConclusionsOur results show that our methods reconstruct small sets of high-quality candidate ancestors with better agreement to the two target structures than with classical approaches. Our program is freely available at: http://csb.cs.mcgill.ca/acharnement.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3105-4) contains supplementary material, which is available to authorized users.

Structural and dynamic insights into the C-terminal extension of cysteine proteinase B from Leishmania amazonensis

Cysteine proteinase B (CPB) is a significant virulence factor for Leishmania infections. Upon processing from its zymogen form, it happens a release of the immunomodulatory CPB C-terminal extension (cyspep) into the cytoplasm of the macrophage. Epitopes derived from this fragment were shown to influence the proportion of lymphocytes CD8+ upon infection, favoring the parasite escaping from the host́s immune system. At present, there is no available structural data of cyspep, which impairs a proper understanding of its biological functions. Here, we attempted to build molecular models for this fragment and subsequently evaluate their stabilities in aqueous solution from molecular dynamics simulations analysis. Characterization of our models obtained with distinct techniques (comparative modeling, threading, and ab initio) indicates a prevalence of β−sheets in agreement with consensus secondary structure predictions. Simulation data supported this finding since the formation of new strands, along with a rapid disruption of helical content, were observed. Overall, this study provides a rationalization of epitope mapping data and an improved understanding of cyspep antigenicity.

ConSSert: Consensus SVM Model for Accurate Prediction of Ordered Secondary Structure.

Accurate prediction of protein secondary structure remains a crucial step in most approaches to the protein-folding problem, yet the prediction of ordered secondary structure, specifically beta-strands, remains a challenge. We developed a consensus secondary structure prediction method, conSSert, which is based on support vector machines (SVM) and provides exceptional accuracy for the prediction of beta-strands with QE accuracy of over 0.82 and a Q2-EH of 0.86. conSSert uses as input probabilities for the three types of secondary structure (helix, strand, and coil) that are predicted by four top performing methods: PSSpred, PSIPRED, SPINE-X, and RAPTOR. conSSert was trained/tested using 4261 protein chains from PDBSelect25, and 8632 chains from PISCES. Further validation was performed using targets from CASP9, CASP10, and CASP11. Our data suggest that poor performance in strand prediction is likely a result of training bias and not solely due to the nonlocal nature of beta-sheet contacts. conSSert is freely available for noncommercial use as a webservice: http://ares.tamu.edu/conSSert/.

Structural and Dynamic Basis for Low-Affinity, High-Selectivity Binding of L-Glutamine by the Glutamine Riboswitch.

Naturally occurring L-glutamine riboswitches occur in cyanobacteria and marine metagenomes, where they reside upstream of genes involved in nitrogen metabolism. By combining X-ray, NMR, and MD, we characterized an L-glutamine-dependent conformational transition in the Synechococcus elongatus glutamine riboswitch from tuning fork to L-shaped alignment of stem segments. This transition generates an open ligand-binding pocket with L-glutamine selectivity enforced by Mg(2+)-mediated intermolecular interactions. The transition also stabilizes the P1 helix through a long-range "linchpin" Watson-Crick G-C pair-capping interaction, while melting a short helix below P1 potentially capable of modulating downstream readout. NMR data establish that the ligand-free glutamine riboswitch in Mg(2+) solution exists in a slow equilibrium between flexible tuning fork and a minor conformation, similar, but not identical, to the L-shaped bound conformation. We propose that an open ligand-binding pocket combined with a high conformational penalty for forming the ligand-bound state provide mechanisms for reducing binding affinity while retaining high selectivity.

Molecular organization of the 5S rDNA gene type II in elasmobranchs

ABSTRACTThe 5S rDNA gene is a non-coding RNA that can be found in 2 copies (type I and type II) in bony and cartilaginous fish. Previous studies have pointed out that type II gene is a paralog derived from type I. We analyzed the molecular organization of 5S rDNA type II in elasmobranchs. Although the structure of the 5S rDNA is supposed to be highly conserved, our results show that the secondary structure in this group possesses some variability and is different than the consensus secondary structure. One of these differences in Selachii is an internal loop at nucleotides 7 and 112. These mutations observed in the transcribed region suggest an independent origin of the gene among Batoids and Selachii. All promoters were highly conserved with the exception of BoxA, possibly due to its affinity to polymerase III. This latter enzyme recognizes a dT4 sequence as stop signal, however in Rajiformes this signal was doubled in length to dT8. This could be an adaptation toward a higher efficiency in the termination process. Our results suggest that there is no TATA box in elasmobranchs in the NTS region. We also provide some evidence suggesting that the complexity of the microsatellites present in the NTS region play an important role in the 5S rRNA gene since it is significantly correlated with the length of the NTS.

RNA motif discovery: a computational overview.

Genomic studies have greatly expanded our knowledge of structural non-coding RNAs (ncRNAs). These RNAs fold into characteristic secondary structures and perform specific-structure dependent biological functions. Hence RNA secondary structure prediction is one of the most well studied problems in computational RNA biology. Comparative sequence analysis is one of the more reliable RNA structure prediction approaches as it exploits information of multiple related sequences to infer the consensus secondary structure. This class of methods essentially learns a global secondary structure from the input sequences. In this paper, we consider the more general problem of unearthing common local secondary structure based patterns from a set of related sequences. The input sequences for example could correspond to 3′ or 5′ untranslated regions of a set of orthologous genes and the unearthed local patterns could correspond to regulatory motifs found in these regions. These sequences could also correspond to in vitro selected RNA, genomic segments housing ncRNA genes from the same family and so on. Here, we give a detailed review of the various computational techniques proposed in literature attempting to solve this general motif discovery problem. We also give empirical comparisons of some of the current state of the art methods and point out future directions of research.ReviewersThis article was reviewed by Dr. Erez Levanon, Dr. Sebastian Maurer-Stroh and Dr. Weixiong Zhang.

Biochemical analysis of hatchet self-cleaving ribozymes.

Hatchet RNAs are members of a novel self-cleaving ribozyme class that was recently discovered by using a bioinformatics search strategy. The consensus sequence and secondary structure of this class includes 13 highly conserved and numerous other modestly conserved nucleotides interspersed among bulges linking four base-paired substructures. A representative hatchet ribozyme from a metagenomic source requires divalent ions such as Mg2+ to promote RNA strand scission with a maximum rate constant of ∼4 min−1. As with all other small self-cleaving ribozymes discovered to date, hatchet ribozymes employ a general mechanism for catalysis involving the nucleophilic attack of a ribose 2′-oxygen atom on an adjacent phosphorus center. Kinetic characteristics of the reaction demonstrate that members of this ribozyme class have an essential requirement for divalent metal ions and that they might have a complex active site that employs multiple catalytic strategies to accelerate RNA cleavage by internal phosphoester transfer.

Structure-Based Alignment and Consensus Secondary Structures for Three HIV-Related RNA Genomes.

HIV and related primate lentiviruses possess single-stranded RNA genomes. Multiple regions of these genomes participate in critical steps in the viral replication cycle, and the functions of many RNA elements are dependent on the formation of defined structures. The structures of these elements are still not fully understood, and additional functional elements likely exist that have not been identified. In this work, we compared three full-length HIV-related viral genomes: HIV-1NL4-3, SIVcpz, and SIVmac (the latter two strains are progenitors for all HIV-1 and HIV-2 strains, respectively). Model-free RNA structure comparisons were performed using whole-genome structure information experimentally derived from nucleotide-resolution SHAPE reactivities. Consensus secondary structures were constructed for strongly correlated regions by taking into account both SHAPE probing structural data and nucleotide covariation information from structure-based alignments. In these consensus models, all known functional RNA elements were recapitulated with high accuracy. In addition, we identified multiple previously unannotated structural elements in the HIV-1 genome likely to function in translation, splicing and other replication cycle processes; these are compelling targets for future functional analyses. The structure-informed alignment strategy developed here will be broadly useful for efficient RNA motif discovery.

Model-Free RNA Sequence and Structure Alignment Informed by SHAPE Probing Reveals a Conserved Alternate Secondary Structure for 16S rRNA.

Discovery and characterization of functional RNA structures remains challenging due to deficiencies in de novo secondary structure modeling. Here we describe a dynamic programming approach for model-free sequence comparison that incorporates high-throughput chemical probing data. Based on SHAPE probing data alone, ribosomal RNAs (rRNAs) from three diverse organisms – the eubacteria E. coli and C. difficile and the archeon H. volcanii – could be aligned with accuracies comparable to alignments based on actual sequence identity. When both base sequence identity and chemical probing reactivities were considered together, accuracies improved further. Derived sequence alignments and chemical probing data from protein-free RNAs were then used as pseudo-free energy constraints to model consensus secondary structures for the 16S and 23S rRNAs. There are critical differences between these experimentally-informed models and currently accepted models, including in the functionally important neck and decoding regions of the 16S rRNA. We infer that the 16S rRNA has evolved to undergo large-scale changes in base pairing as part of ribosome function. As high-quality RNA probing data become widely available, structurally-informed sequence alignment will become broadly useful for de novo motif and function discovery.

P147 A trust-wide audit on pelvic inflammatory disease management in a genitourinary medicine setting

<h3>Abstract</h3> SARS-CoV-2 is the causative agent of worldwide pandemic disease COVID-19. SARS-CoV-2 bears positive sense RNA genome, that have organized and complex pattern of replication/transcription process including the generation of subgenomic RNAs. Transcription regulatory sequences (TRS) have important role in the pausing of replication/transcription and generation of subgenomic RNAs. In the present bioinformatics analysis a consensus secondary structure was identified among negative sense subgenomic RNAs at the adjacent of initiation codon. This study proposed that consensus structured domain could involve in mediating the long pausing of replication/transcription complex and responsible for subgenomic RNA production.

Thermodynamic matchers for the construction of the cuckoo RNA family

RNA family models describe classes of functionally related, non-coding RNAs based on sequence and structure conservation. The most important method for modeling RNA families is the use of covariance models, which are stochastic models that serve in the discovery of yet unknown, homologous RNAs. However, the performance of covariance models in finding remote homologs is poor for RNA families with high sequence conservation, while for families with high structure but low sequence conservation, these models are difficult to built in the first place. A complementary approach to RNA family modeling involves the use of thermodynamic matchers. Thermodynamic matchers are RNA folding programs, based on the established thermodynamic model, but tailored to a specific structural motif. As thermodynamic matchers focus on structure and folding energy, they unfold their potential in discovering homologs, when high structure conservation is paired with low sequence conservation. In contrast to covariance models, construction of thermodynamic matchers does not require an input alignment, but requires human design decisions and experimentation, and hence, model construction is more laborious. Here we report a case study on an RNA family that was constructed by means of thermodynamic matchers. It starts from a set of known but structurally different members of the same RNA family. The consensus secondary structure of this family consists of 2 to 4 adjacent hairpins. Each hairpin loop carries the same motif, CCUCCUCCC, while the stems show high variability in their nucleotide content. The present study describes (1) a novel approach for the integration of the structurally varying family into a single RNA family model by means of the thermodynamic matcher methodology, and (2) provides the results of homology searches that were conducted with this model in a wide spectrum of bacterial species.

Cloning and characterization of a novel sigma-like glutathione S-transferase from the giant panda parasitic nematode, Baylisascaris schroederi.

BackgroundBaylisascaris schroederi, an intestinal nematode of the giant panda, is the cause of the often fatal disease, baylisascariasis. Glutathione S-transferases (GSTs) are versatile enzymes that can affect parasite survival and parasite-host interactions and, are therefore, potential targets for the development of diagnostic tests and vaccines.MethodsIn this study, we identified a full-length cDNA that encoded a novel, secretory sigma-like GST (Bsc-GSTσ) from a B. schroederi-omic dataset. Following cloning and sequencing, sequence and structural analyses and comparative modeling were performed using online-bioinformatics and proteomics tools. The recombinant Bsc-GSTσ (rBsc-GSTσ) protein was prokaryotically expressed and then used to detect antigenicity and reactivity using immunoblotting assays. In addition, the native protein in female adult B. schroederi was located via immunofluorescence techniques, while the preliminary ELISA-based serodiagnostic potential of rBsc-GSTσ was assessed in native and infected mouse sera.ResultsBsc-GSTσ contained a 621-bp open reading frame that encoded a polypeptide of 206 amino acids with two typical sigma GST domain profiles, including a GST_N_Sigma_like at the N-terminus and a GST_C_Sigma_like at the C-terminus. The presence of an N-terminal signal sequence indicated that Bsc-GSTσ was a secretory protein. Sequence alignment and phylogenetic analyses showed that Bsc-GSTσ was a nematode-specific member of the Sigma class GSTs and shared the closest genetic distance with its homologue in Ascaris suum. Further comparative structure analyses indicated that Bsc-GSTσ possessed the essential structural motifs (e.g., βαβαββα) and the consensus secondary or tertiary structure that is typical for other characterized GSTσs. Immunolocalization revealed strong distributions of native Bsc-GSTσ in the body hypodermis, lateral chords, gut epithelium, gut microvilli, oviduct epithelium, and ovaries of adult female worms, similar to its homologue in A. suum. Building on good immunogenic properties, rBsc-GSTσ-based ELISA exhibited a sensitivity of 79.1% and a specificity of 82.0% to detect anti-B. schroederi IgG antibodies in the sera of experimentally infected mice.ConclusionThis study presents a comprehensive demonstration of sequence and structural-based analysis of a new, secretory sigma-like GST from a nematode, and its good serodiagnostic performance suggests that rBsc-GSTσ has the potential to detect B. schroederi and, therefore, could be used to develop an ELISA-based serological test to diagnose baylisascariasis in giant pandas.

Structure-Based Alignment and Consensus Secondary Structures for Three HIV-Related RNA Genomes

Structure-Based Alignment and Consensus Secondary Structures for Three HIV-Related RNA Genomes

ITS1, 5.8S and ITS2 secondary structure modelling for intra-specific differentiation among species of the Colletotrichum gloeosporioides sensu lato species complex.

The Colletotrichum gloeosporioides species complex is among the most destructive fungal plant pathogens in the world, however, identification of member species which are of quarantine importance is impacted by a number of factors that negatively affect species identification. Structural information of the rRNA marker may be considered to be a conserved marker which can be used as supplementary information for possible species identification. The difficulty in using ITS rDNA sequences for identification lies in the low level of sequence variation at the intra-specific level and the generation of artificially-induced sequence variation due to errors in polymerization of the ITS array during DNA replication. Type and query ITS sequences were subjected to sequence analyses prior to generation of predicted consensus secondary structures, including the pattern of nucleotide polymorphisms and number of indel haplotypes, GC content, and detection of artificially-induced sequence variation. Data pertaining to structure stability, the presence of conserved motifs in secondary structures and mapping of all sequences onto the consensus C. gloeosporioides sensu stricto secondary structure for ITS1, 5.8S and ITS2 markers was then carried out. Motifs that are evolutionarily conserved among eukaryotes were found for all ITS1, 5.8S and ITS2 sequences. The sequences exhibited conserved features typical of functional rRNAs. Generally, polymorphisms occurred within less conserved regions and were seen as bulges, internal and terminal loops or non-canonical G–U base-pairs within regions of the double stranded helices. Importantly, there were also taxonomic motifs and base changes that were unique to specific taxa and which may be used to support intra-specific identification of members of the C. gloeosporioides sensu lato species complex.Electronic supplementary materialThe online version of this article (doi:10.1186/2193-1801-3-684) contains supplementary material, which is available to authorized users.

Rfam 12.0: updates to the RNA families database.

The Rfam database (available at http://rfam.xfam.org) is a collection of non-coding RNA families represented by manually curated sequence alignments, consensus secondary structures and annotation gathered from corresponding Wikipedia, taxonomy and ontology resources. In this article, we detail updates and improvements to the Rfam data and website for the Rfam 12.0 release. We describe the upgrade of our search pipeline to use Infernal 1.1 and demonstrate its improved homology detection ability by comparison with the previous version. The new pipeline is easier for users to apply to their own data sets, and we illustrate its ability to annotate RNAs in genomic and metagenomic data sets of various sizes. Rfam has been expanded to include 260 new families, including the well-studied large subunit ribosomal RNA family, and for the first time includes information on short sequence- and structure-based RNA motifs present within families.

ITS rDNA sequence comparisons resolve phylogenetic relationships in Orostachys subsection Appendiculatae (Crassulaceae)

Orostachys (Crassulaceae) is a small genus of succulent plants having a predominantly East Asian distribution. Recent DNA sequence comparisons revealed polyphyletic nature of the genus and found distant relationship between its infrageneric taxa. Here we present the first molecular phylogeny of Orostachys subsection Appendiculatae based on a large number of ITS rDNA sequences representing most currently recognized members of the subsection and utilizing secondary structure information. Ribosomal spacer was a highly informative marker and provided a phylogenetic signal sufficient to resolve relationships at different scales, from affinities between species to a fine geographic structure in broadly sampled species. It was also conservative enough to allow unambiguous alignment and construction of consensus secondary structure models for ITS1 and ITS2. These models displayed a number of molecular synapomorphies defining most lineages established in our analyses. We revealed a major split in the subsection placing three species, O. spinosa, O. japonica and O. chanetii, into a strongly supported clade to the exclusion of O. thyrsiflora. Phenotypically distinct monotypic genus Meterostachys was also resolved as a part of the subsection’s clade and showed affinity to O. thyrsiflora. Our data suggested that morphology-based species concept for O. thyrsiflora requires reassessment.