Dataset Information

Frnakenstein: multiple target inverse RNA folding.

ABSTRACT: RNA secondary structure prediction, or folding, is a classic problem in bioinformatics: given a sequence of nucleotides, the aim is to predict the base pairs formed in its three dimensional conformation. The inverse problem of designing a sequence folding into a particular target structure has only more recently received notable interest. With a growing appreciation and understanding of the functional and structural properties of RNA motifs, and a growing interest in utilising biomolecules in nano-scale designs, the interest in the inverse RNA folding problem is bound to increase. However, whereas the RNA folding problem from an algorithmic viewpoint has an elegant and efficient solution, the inverse RNA folding problem appears to be hard.In this paper we present a genetic algorithm approach to solve the inverse folding problem. The main aims of the development was to address the hitherto mostly ignored extension of solving the inverse folding problem, the multi-target inverse folding problem, while simultaneously designing a method with superior performance when measured on the quality of designed sequences. The genetic algorithm has been implemented as a Python program called Frnakenstein. It was benchmarked against four existing methods and several data sets totalling 769 real and predicted single structure targets, and on 292 two structure targets. It performed as well as or better at finding sequences which folded in silico into the target structure than all existing methods, without the heavy bias towards CG base pairs that was observed for all other top performing methods. On the two structure targets it also performed well, generating a perfect design for about 80% of the targets.Our method illustrates that successful designs for the inverse RNA folding problem does not necessarily have to rely on heavy biases in base pair and unpaired base distributions. The design problem seems to become more difficult on larger structures when the target structures are real structures, while no deterioration was observed for predicted structures. Design for two structure targets is considerably more difficult, but far from impossible, demonstrating the feasibility of automated design of artificial riboswitches. The Python implementation is available at http://www.stats.ox.ac.uk/research/genome/software/frnakenstein.

SUBMITTER: Lyngso RB

PROVIDER: S-EPMC3534541 | biostudies-literature | 2012 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Frnakenstein: multiple target inverse RNA folding.

Lyngsø Rune B RB Anderson James W J JW Sizikova Elena E Badugu Amarendra A Hyland Tomas T Hein Jotun J

BMC bioinformatics 20121009

<h4>Background</h4>RNA secondary structure prediction, or folding, is a classic problem in bioinformatics: given a sequence of nucleotides, the aim is to predict the base pairs formed in its three dimensional conformation. The inverse problem of designing a sequence folding into a particular target structure has only more recently received notable interest. With a growing appreciation and understanding of the functional and structural properties of RNA motifs, and a growing interest in utilising ...[more]

PMID: 23043260

Similar Datasets

Project description:Saccharomyces cerevisiae telomerase RNA, TLC1, is an 1157 nt non-coding RNA that functions as both a template for DNA synthesis and a flexible scaffold for telomerase RNP holoenzyme protein subunits. The tractable budding yeast system has provided landmark discoveries about telomere biology in vivo, but yeast telomerase research has been hampered by the fact that the large TLC1 RNA subunit does not support robust telomerase activity in vitro. In contrast, 155-500 nt miniaturized TLC1 alleles comprising the catalytic core domain and lacking the RNA's long arms do reconstitute robust activity. We hypothesized that full-length TLC1 is prone to misfolding in vitro. To create a full-length yeast telomerase RNA, predicted to fold into its biologically relevant structure, we took an inverse RNA-folding approach, changing 59 nucleotides predicted to increase the energetic favorability of folding into the modeled native structure based on the p-num feature of Mfold software. The sequence changes lowered the predicted ∆G of this "determined-arm" allele, DA-TLC1, by 61 kcal/mol (-19%) compared to wild-type. We tested DA-TLC1 for reconstituted activity and found it to be ~5-fold more robust than wild-type TLC1, suggesting that the inverse-folding design indeed improved folding in vitro into a catalytically active conformation. We also tested if DA-TLC1 functions in vivo, discovering that it complements a tlc1∆ strain, allowing cells to avoid senescence and maintain telomeres of nearly wild-type length. However, all inverse-designed RNAs that we tested had reduced abundance in vivo. In particular, inverse-designing nearly all of the Ku arm caused a profound reduction in telomerase RNA abundance in the cell and very short telomeres. Overall, these results show that the inverse design of S. cerevisiae telomerase RNA increases activity in vitro, while reducing abundance in vivo. This study provides a biochemically and biologically tested approach to inverse-design RNAs using Mfold that could be useful for controlling RNA structure in basic research and biomedicine.

Dataset Information

Frnakenstein: multiple target inverse RNA folding.

Publications

Frnakenstein: multiple target inverse RNA folding.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets