A Faster Algorithm for RNA Co-folding

  • Authors:
  • Michal Ziv-Ukelson;Irit Gat-Viks;Ydo Wexler;Ron Shamir

  • Affiliations:
  • Computer Science Department, Ben Gurion University of the Negev, Beer-Sheva,;Computational Molecular Biology Department, Max Planck Institute for Molecular Genetics, Berlin, Germany;Microsoft Research, Microsoft Corporation, Redmond,;School of Computer Science, Tel Aviv University,

  • Venue:
  • WABI '08 Proceedings of the 8th international workshop on Algorithms in Bioinformatics
  • Year:
  • 2008

Quantified Score

Hi-index 0.00

Visualization

Abstract

The current pairwise RNA (secondary) structural alignment algorithms are based on Sankoff's dynamic programming algorithm from 1985. Sankoff's algorithm requires O(N6) time and O(N4) space, where Ndenotes the length of the compared sequences, and thus its applicability is very limited. The current literature offers many heuristics for speeding up Sankoff's alignment process, some making restrictive assumptions on the length or the shape of the RNA substructures. We show how to speed up Sankoff's algorithm in practice via non-heuristic methods, without compromising optimality. Our analysis shows that the expected time complexity of the new algorithm is O(N4驴(N)), where 驴(N) converges to O(N), assuming a standard polymer folding model which was supported by experimental analysis. Hence our algorithm speeds up Sankoff's algorithm by a linear factor on average. In simulations, our algorithm speeds up computation by a factor of 3-12 for sequences of length 25-250.Availability:Code and data sets are available, upon request.