Predicting a set of minimal free energy RNA secondary structures common to two sequences

Abstract

Motivation: Function derives from structure, therefore, there is need for methods to predict functional RNA structures. Results: The Dynalign algorithm, which predicts the lowest free energy secondary structure common to two unaligned RNA sequences, is extended to the prediction of a set of low-energy structures. Dot plots can be drawn to show all base pairs in structures within an energy increment. Dynalign predicts more well-defined structures than structure prediction using a single sequence; in 5S rRNA sequences, the average number of base pairs in structures with energy within 20% of the lowest energy structure is 317 using Dynalign, but 569 using a single sequence. Structure prediction with Dynalign can also be constrained according to experiment or comparative analysis. The accuracy, measured as sensitivity and positive predictive value, of Dynalign is greater than predictions with a single sequence. Availability: Dynalign can be downloaded at http://rna.urmc.rochester.edu Contact:David_Mathews@urmc.rochester.edu Supplementary information:http://rna.urmc.rochester.edu/supplementary.html