Title: Searching for repeats, as an example of using the generalised Ruzzo-Tompa algorithm to find optimal subsequences with gaps
Authors: John L. Spouge; Leonardo Mariño-Ramírez; Sergey L. Sheetlin
Addresses: Computational Biology Branch, National Center for Biotechnology Information, Bethesda, MD 20894, USA ' Computational Biology Branch, National Center for Biotechnology Information, Bethesda, MD 20894, USA ' Computational Biology Branch, National Center for Biotechnology Information, Bethesda, MD 20894, USA
Abstract: Some biological sequences contain subsequences of unusual composition; e.g. some proteins contain DNA binding domains, transmembrane regions and charged regions, and some DNA sequences contain repeats. The linear-time Ruzzo-Tompa (RT) algorithm finds subsequences of unusual composition, using a sequence of scores as input and the corresponding 'maximal segments' as output. In principle, permitting gaps in the output subsequences could improve sensitivity. Here, the input of the RT algorithm is generalised to a finite, totally ordered, weighted graph, so the algorithm locates paths of maximal weight through increasing but not necessarily adjacent vertices. By permitting the penalised deletion of unfavourable letters, the generalisation therefore includes gaps. The program RepWords, which finds inexact simple repeats in DNA, exemplifies the general concepts by out-performing a similar extant, ad hoc tool. With minimal programming effort, the generalised Ruzzo-Tompa algorithm could improve the performance of many programs for finding biological subsequences of unusual composition.
Keywords: generalised Ruzzo-Tompa algorithm; optimal subsequences; unusual composition; gaps; repeats; biological sequences; DNA sequences.
DOI: 10.1504/IJBRA.2014.062991
International Journal of Bioinformatics Research and Applications, 2014 Vol.10 No.4/5, pp.384 - 408
Published online: 24 Oct 2014 *
Full-text access for editors Full-text access for subscribers Purchase this article Comment on this article