Foldalign: RNA Structure and Sequence Alignment
IntroductionFOLDALIGN is a tool for aligning RNA structures. It folds and aligns (make a foldalignment) based on a lightweight energy model and sequence similarity. The current version 2.5.0 makes pairwise foldalignments. The foldalignments can be local and global. Local foldalignments can be extended into scans of the two sequences by limiting the local alignment length. When scanning FOLDALIGN program outputs the best local alignment score for each pair of positions in the two sequences. This list can be further processed to produce a list of non-overlapping local alignments between the two sequences.
AvailabilityFOLDALIGN is available as a software package, or as a web-server.
Abstract from Bioinformatics, 21(9): 1815-1824, 2005Motivation: Searching for non-coding RNA (ncRNA) genes and structural RNA elements (eleRNA) are major challenges in gene finding today as these often are conserved in structure rather than in sequence. Even though the number of available methods is growing, it is still of interest to pairwise detect two genes with low sequence similarity, where the genes are part of a larger genomic region.
Results: Here we present such an approach for pairwise local alignment which is based on FOLDALIGN and the Sankoff algorithm for simultaneous structural alignment of multiple sequences. We include the ability to conduct mutual scans of two sequences of arbitrary length while searching for common local structural motifs of some maximum length. This drastically reduces the complexity of the algorithm. The scoring scheme includes structural parameters corresponding to those available for free energy as well as for substitution matrices similar to RIBOSUM. The new FOLDALIGN implementation is tested on a dataset where the ncRNAs and eleRNAs have sequence similarity <40% and where the ncRNAs and eleRNAs are energetically indistinguishable from the surrounding genomic sequence context. The method is tested in two ways: (1) its ability to find the common structure between the genes only and (2) its ability to locate ncRNAs and eleRNAs in a genomic context. In case (1), it makes sense to compare with methods like Dynalign, and the performances are very similar, but FOLDALIGN is substantially faster. The structure prediction performance for a family is typically around 0.7 using Matthews correlation coefficient. In case (2), the algorithm is successful at locating RNA families with an average sensitivity of 0.8 and a positive predictive value of 0.9 using a BLAST-like hit selection scheme.
Availability: The program is available online at http://foldalign.kvl.dk/
Fast Pairwise Structural RNA Alignments by Pruning of the Dynamical Programming Matrix
Havgaard JH, Torarinsson E, Gorodkin J PLoS Computational Biology 2007, 3(10):e193
[ PubMed | Paper | Software | Sup ]
The FOLDALIGN web server for pairwise structural RNA alignment and mutual motif search
Havgaard JH, Lyngso RB, Gorodkin J Nucleic Acids Research 2005, 33(Web Server):W650-W653
[ PubMed | Paper | Webserver | Software | Sup ]