Visualizing RNA base-pairing probabilities with RNAbow diagrams.
Ontology highlight
ABSTRACT:
SUBMITTER: Aalberts DP
PROVIDER: S-EPMC3677257 | biostudies-other | 2013 Apr
REPOSITORIES: biostudies-other
Publications
Visualizing RNA base-pairing probabilities with RNAbow diagrams.
Aalberts Daniel P DP Jannen William K WK
RNA (New York, N.Y.) 20130213 4
There are many effective ways to represent a minimum free energy RNA secondary structure that make it easy to locate its helices and loops. It is a greater challenge to visualize the thermal average probabilities of all folds in a partition function sum; dot plot representations are often puzzling. Therefore, we introduce the RNAbows visualization tool for RNA base pair probabilities. RNAbows represent base pair probabilities with line thickness and shading, yielding intuitive diagrams. RNAbows ...[more]
Similar Datasets
Project description:MotivationRNA secondary structure prediction is widely used to understand RNA function. Recently, there has been a shift away from the classical minimum free energy methods to partition function-based methods that account for folding ensembles and can therefore estimate structure and base pair probabilities. However, the classical partition function algorithm scales cubically with sequence length, and is therefore prohibitively slow for long sequences. This slowness is even more severe than cubic-time free energy minimization due to a substantially larger constant factor in runtime.ResultsInspired by the success of our recent LinearFold algorithm that predicts the approximate minimum free energy structure in linear time, we design a similar linear-time heuristic algorithm, LinearPartition, to approximate the partition function and base-pairing probabilities, which is shown to be orders of magnitude faster than Vienna RNAfold and CONTRAfold (e.g. 2.5 days versus 1.3 min on a sequence with length 32 753 nt). More interestingly, the resulting base-pairing probabilities are even better correlated with the ground-truth structures. LinearPartition also leads to a small accuracy improvement when used for downstream structure prediction on families with the longest length sequences (16S and 23S rRNAs), as well as a substantial improvement on long-distance base pairs (500+ nt apart).Availability and implementationCode: http://github.com/LinearFold/LinearPartition; Server: http://linearfold.org/partition.Supplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC7355276 | biostudies-literature
Project description:An RNA switch triggers biological functions by toggling between two conformations. RNA switches include bacterial riboswitches, where ligand binding can stabilize a bound structure. For RNAs with only one stable structure, structural prediction usually just requires a straightforward free energy minimization, but for an RNA switch, the prediction of a less stable alternative structure is often computationally costly and even problematic. The current sampling-clustering method predicts stable and alternative structures by partitioning structures sampled from the energy landscape into two clusters, but it is very time-consuming. Instead, we predict the alternative structure of an RNA switch from conditional probability calculations within the energy landscape. First, our method excludes base pairs related to the most stable structure in the energy landscape. Then, it detects stable stems ("seeds") in the remaining landscape. Finally, it folds an alternative structure prediction around a seed. While having comparable riboswitch classification performance, the conditional-probability computations had fewer adjustable parameters, offered greater predictive flexibility, and were more than one thousand times faster than the sampling step alone in sampling-clustering predictions, the competing standard. Overall, the described approach helps traverse thermodynamically improbable energy landscapes to find biologically significant substructures and structures rapidly and effectively.
| S-EPMC6561571 | biostudies-literature
Project description:Natural RNAs utilize extensive chemical modifications to diversify their structures and functions. 2-Thiouridine geranylation is a special hydrophobic tRNA modification that has been discovered very recently in several bacteria, such as Escherichia coli, Enterobacter aerogenes, Pseudomonas aeruginosa and Salmonella Typhimurium The geranylated residues are located in the first anticodon position of tRNAs specific for lysine, glutamine and glutamic acid. This big hydrophobic terpene functional group affects the codon recognition patterns and reduces frameshifting errors during translation. We aimed to systematically study the structure, function and biosynthesis mechanism of this geranylation pathway, as well as answer the question of why nature uses such a hydrophobic modification in hydrophilic RNA systems. Recently, we have synthesized the deoxy-analog of S-geranyluridine and showed the geranylated T-G pair is much stronger than the geranylated T-A pair and other mismatched pairs in the B-form DNA duplex context, which is consistent with the observation that the geranylated tRNA(Glu) UUC recognizes GAG more efficiently than GAA. In this manuscript we report the synthesis and base pairing specificity studies of geranylated RNA oligos. We also report extensive molecular simulation studies to explore the structural features of the geranyl group in the context of A-form RNA and its effect on codon-anticodon interaction during ribosome binding.
| S-EPMC5291276 | biostudies-literature
Project description:The ribonuclease Dicer produces microRNAs (miRNAs) and small interfering RNAs that are handed over to Ago proteins to control gene expression by targeting complementary sequences within transcripts. Interestingly, a growing number of reports have demonstrated that the activity of Dicer may extend beyond the biogenesis of small regulatory RNAs. Among them, a report from our latest studies revealed that human Dicer facilitates base pairing of complementary sequences present in two nucleic acids, thus acting as a nucleic acid annealer. Accordingly, in this manuscript, we address how RNA structure influences the annealing activity of human Dicer. We show that Dicer supports hybridization between a small RNA and a complementary sequence of a longer RNA in vitro, even when both complementary sequences are trapped within secondary structures. Moreover, we show that under applied conditions, human Ago2, a core component of RNA-induced silencing complex, displays very limited annealing activity. Based on the available data from new-generation sequencing experiments regarding the RNA pool bound to Dicer in vivo, we show that multiple Dicer-binding sites within mRNAs also contain miRNA targets. Subsequently, we demonstrate in vitro that Dicer but not Ago2 can anneal miRNA to its target present within mRNA. We hypothesize that not all miRNA duplexes are handed over to Ago proteins. Instead, miRNA-Dicer complexes could target specific sequences within transcripts and either compete or cooperate for binding sites with miRNA-Ago complexes. Thus, not only Ago but also Dicer might be directly involved in the posttranscriptional control of gene expression.
| S-EPMC7391396 | biostudies-literature
Project description:RNA-small molecule binding is a key regulatory mechanism which can stabilize 3D structures and activate molecular functions. The discovery of RNA-targeting compounds is thus a current topic of interest for novel therapies. Our work is a first attempt at bringing the scalability and generalization abilities of machine learning methods to the problem of RNA drug discovery, as well as a step towards understanding the interactions which drive binding specificity. Our tool, RNAmigos, builds and encodes a network representation of RNA structures to predict likely ligands for novel binding sites. We subject ligand predictions to virtual screening and show that we are able to place the true ligand in the 71st-73rd percentile in two decoy libraries, showing a significant improvement over several baselines, and a state of the art method. Furthermore, we observe that augmenting structural networks with non-canonical base pairing data is the only representation able to uncover a significant signal, suggesting that such interactions are a necessary source of binding specificity. We also find that pre-training with an auxiliary graph representation learning task significantly boosts performance of ligand prediction. This finding can serve as a general principle for RNA structure-function prediction when data is scarce. RNAmigos shows that RNA binding data contains structural patterns with potential for drug discovery, and provides methodological insights for possible applications to other structure-function learning tasks. The source code, data and a Web server are freely available at http://rnamigos.cs.mcgill.ca.
| S-EPMC7430648 | biostudies-literature
Project description:Elucidation of the folded structures that RNA forms in vivo is vital to understanding its functions. Chemical reagents that modify the Watson-Crick (WC) face of unprotected nucleobases are particularly useful in structure elucidation. Dimethyl sulfate penetrates cell membranes and informs on RNA base-pairing and secondary structure but only modifies the WC face of adenines and cytosines. We present glyoxal, methylglyoxal, and phenylglyoxal as potent in vivo reagents that target the WC face of guanines as well as cytosines and adenines. Tests on rice (Oryza sativa) 5.8S rRNA in vitro read out by reverse transcription and gel electrophoresis demonstrate specific modification of almost all guanines in a time- and pH-dependent manner. Subsequent in vivo tests on rice, a eukaryote, and Bacillus subtilis and Escherichia coli, Gram-positive and Gram-negative bacteria, respectively, showed that all three reagents enter living cells without prior membrane permeabilization or pH adjustment of the surrounding media and specifically modify solvent-exposed guanine, cytosine, and adenine residues.
| S-EPMC5733565 | biostudies-literature
Project description:This paper presents the first empirical investigation that compares Euler and linear diagrams when they are used to represent set cardinality. A common approach is to use area-proportional Euler diagrams but linear diagrams can exploit length-proportional straight-lines for the same purpose. Another common approach is to use numerical annotations. We first conducted two empirical studies, one on Euler diagrams and the other on linear diagrams. These suggest that area-proportional Euler diagrams with numerical annotations and length-proportional linear diagrams without numerical annotations support significantly better task performance. We then conducted a third study to investigate which of these two notations should be used in practice. This suggests that area-proportional Euler diagrams with numerical annotations most effectively supports task performance and so should be used to visualize set cardinalities. However, these studies focused on data that can be visualized reasonably accurately using circles and the results should be taken as valid within that context. Future work needs to determine whether the results generalize both to when circles cannot be used and for other ways of encoding cardinality information.
| S-EPMC6438608 | biostudies-literature
Project description:RNA structure and dynamics are critical to biological function. However, strategies for determining RNA structure in vivo are limited, with established chemical probing and newer duplex detection methods each having deficiencies. Here we convert the common reagent dimethyl sulfate into a useful probe of all 4 RNA nucleotides. Building on this advance, we introduce PAIR-MaP, which uses single-molecule correlated chemical probing to directly detect base-pairing interactions in cells. PAIR-MaP has superior resolution compared to alternative experiments, can resolve multiple sets of pairing interactions for structurally dynamic RNAs, and enables highly accurate structure modeling, including of RNAs containing multiple pseudoknots and extensively bound by proteins. Application of PAIR-MaP to human RNase MRP and 2 bacterial messenger RNA 5' untranslated regions reveals functionally important and complex structures undetected by prior analyses. PAIR-MaP is a powerful, experimentally concise, and broadly applicable strategy for directly visualizing RNA base pairs and dynamics in cells.
| S-EPMC6900531 | biostudies-literature
Project description:Pairs of nucleotides within functional nucleic acid secondary structures often display evidence of coevolution that is consistent with the maintenance of base-pairing. Here, we introduce a sequence evolution model, MESSI (Modeling the Evolution of Secondary Structure Interactions), that infers coevolution associated with base-paired sites in DNA or RNA sequence alignments. MESSI can estimate coevolution while accounting for an unknown secondary structure. MESSI can also use graphics processing unit parallelism to increase computational speed. We used MESSI to infer coevolution associated with GC, AU (AT in DNA), GU (GT in DNA) pairs in noncoding RNA alignments, and in single-stranded RNA and DNA virus alignments. Estimates of GU pair coevolution were found to be higher at base-paired sites in single-stranded RNA viruses and noncoding RNAs than estimates of GT pair coevolution in single-stranded DNA viruses. A potential biophysical explanation is that GT pairs do not stabilize DNA secondary structures to the same extent that GU pairs do in RNA. Additionally, MESSI estimates the degrees of coevolution at individual base-paired sites in an alignment. These estimates were computed for a SHAPE-MaP-determined HIV-1 NL4-3 RNA secondary structure. We found that estimates of coevolution were more strongly correlated with experimentally determined SHAPE-MaP pairing scores than three nonevolutionary measures of base-pairing covariation. To assist researchers in prioritizing substructures with potential functionality, MESSI automatically ranks substructures by degrees of coevolution at base-paired sites within them. Such a ranking was created for an HIV-1 subtype B alignment, revealing an excess of top-ranking substructures that have been previously identified as having structure-related functional importance, among several uncharacterized top-ranking substructures.
| S-EPMC6993869 | biostudies-literature
Project description:The carbon storage regulator system and base-pairing small RNAs (sRNAs) represent two predominant modes of bacterial post-transcriptional regulation, which globally influence gene expression. Binding of CsrA protein to the 5' UTR or initial mRNA coding sequences can affect translation, RNA stability, and/or transcript elongation. Base-pairing sRNAs also regulate these processes, often requiring assistance from the RNA chaperone Hfq. Transcriptomics studies in Escherichia coli have identified many new CsrA targets, including Spot 42 and other base-pairing sRNAs. Spot 42 synthesis is repressed by cAMP-CRP, induced by the presence of glucose, and Spot 42 post-transcriptionally represses operons that facilitate metabolism of nonpreferred carbon sources. CsrA activity is also increased by glucose via effects on CsrA sRNA antagonists, CsrB/C. Here, we elucidate a mechanism wherein CsrA binds to and protects Spot 42 sRNA from RNase E-mediated cleavage. This protection leads to enhanced repression of srlA by Spot 42, a gene required for sorbitol uptake. A second, independent mechanism by which CsrA represses srlA is by binding to and inhibiting translation of srlM mRNA, encoding a transcriptional activator of srlA. Our findings demonstrate a novel means of regulation, by CsrA binding to a sRNA, and indicate that such interactions can help to shape complex bacterial regulatory circuitry.
| S-EPMC10000020 | biostudies-literature