ABSTRACT
Recent advances in next-generation sequencing technology have significantly promoted high-throughput experimental probing of RNA secondary structures.
The resulting enzymatic or chemical probing information is then incorporated into a minimum free energy folding algorithm to predict more accurate RNA secondary structures.
A drawback of this approach is that it does not consider the presence of alternative RNA structures.
In addition, the alternative RNA structures may contaminate experimental probing information of each other and direct the minimum free energy folding to a wrong direction.
In this article, we present a combinatorial solution for this problem, where two alternative structures can be folded simultaneously given the experimental probing information regarding the mixture of these two alternative structures.
We have tested our algorithm with artificially generated mixture probing data on adenine riboswitch and TPP riboswitch.
The experimental results show that our algorithm can successfully recover the ON and OFF structures of these riboswitches.
Full Article
README
Installation
'make'
Usage
./MutualFold -i sequence_file -r reactivity_file [-l min_stack_len -b min_stack_hydro_bonds -f frequency]
sequence_file (-i):
The file containning the sequence to be folded in FASTA format.
reactivity_file (-r):
The file containning the SHAPE reactivity for the RNA transcript.
min_stack_len (-l):
length for a stack to be considered significant. Default: 4
min_stack_hydro_bonds (-b):
Minimum number of hydrogen bonds required for a stack to be considered significant. Default: 8
frequency (-f):
Estimated frequency for the first transcript in percentage (between 1 to 99). Default: 50
DOWNLOADS