Fast and accurate structural RNA alignment by progressive lagrangian optimization

  • Authors:
  • Markus Bauer;Gunnar W. Klau;Knut Reinert

  • Affiliations:
  • Algorithmic Bioinformatics Group, Institute of Computer Science, Free University of Berlin, Germany;Mathematics in Life Sciences Group, Institute of Mathematics, Free University of Berlin, Germany;Algorithmic Bioinformatics Group, Institute of Computer Science, Free University of Berlin, Germany

  • Venue:
  • CompLife'05 Proceedings of the First international conference on Computational Life Sciences
  • Year:
  • 2005

Quantified Score

Hi-index 0.00

Visualization

Abstract

During the last few years new functionalities of RNA have been discovered, renewing the need for computational tools for their analysis. To this respect, multiple sequence alignment is an essential step in finding structurally conserved regions in related RNA sequences. In contrast to proteins, many classes of functionally related RNA molecules show a rather weak sequence conservation but instead a fairly well conserved secondary structure. Hence, any method that relates RNA sequences in form of multiple alignments should take structural features into account, which has been verified in recent studies. Progress has been made in developing new structural alignment algorithms, however, current methods are computationally costly or do not have the desired accuracy to make them an everyday tool. In this paper we present a fast, practical, and accurate method for computing multiple, structural RNA alignments. The method is based on combining a new pairwise structural alignment method with the popular program T-Coffee. Our pairwise method is based on an integer linear programming (ILP) formulation resulting from a graph-theoretic reformulation of the structural alignment problem. We find provably optimal or near-optimal solutions of the ILP with a Lagrangian approach. Tests on a recently published benchmark set show that our Lagrangian approach outperforms current programs in quality and in the length of the sequences it can align.