Bioinformatics Advance Access originally published online on August 11, 2009
Bioinformatics 2009 25(20):2646-2654; doi:10.1093/bioinformatics/btp481
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Partition function and base pairing probabilities for RNA–RNA interaction prediction
1 Center for Combinatorics, LPMC-TJKLC, 2 College of Life Science, Nankai University Tianjin 300071, P.R. China, 3 Bioinformatics Group, Department of Computer Science, and Interdisciplinary Center for Bioinformatics, University of Leipzig, Härtelstrasse 16-18, D-04107 Leipzig, 4 Max Planck Institute for Mathematics in the Sciences, Inselstrasse 22, 5 RNomics Group, Fraunhofer Institut for Cell Therapy and Immunology, Perlickstraße 1,D-04103 Leipzig, Germany, 6 Institute for Theoretical Chemistry, University of Vienna, Währingerstrasse 17, A-1090 Vienna, Austria and 7 The Santa Fe Institute, 1399 Hyde Park Rd., Santa Fe, NM, USA
* To whom correspondence should be addressed.
| Abstract |
|---|
Motivation: The RNA–RNA interaction problem (RIP) consists in finding the energetically optimal structure of two RNA molecules that bind to each other. The standard model allows secondary structures in both partners as well as additional base pairs between the two RNAs subject to certain restrictions that ensure that RIP is solvabale by a polynomial time dynamic programming algorithm. RNA–RNA binding, like RNA folding, is typically not dominated by the ground state structure. Instead, a large ensemble of alternative structures contributes to the interaction thermodynamics.
Results: We present here an O(N6) time and O(N4) dynamics programming algorithm for computing the full partition function for RIP which is based on the combinatorial notion of tight structures. Albeit equivalent to recent work by H. Chitsaz and collaborators, our approach in addition provides a full-fledged computation of the base pairing probabilities, which relies on the notion of a decomposition tree for joint structures. In practise, our implementation is efficient enough to investigate, for instance, the interactions of small bacterial RNAs and their target mRNAs.
Availability: The program rip is implemented in C. The source code is available for download from http://www.combinatorics.cn/cbpc/rip.html and http://www.bioinf.uni-leipzig.de/Software/rip.html.
Contact: duck{at}santafe.edu
Supplementary information: Supplementary data are available at Bioinformatics online.
Associate Editor: Ivo Hofacker
Received on April 16, 2009; revised on June 21, 2009; accepted on August 2, 2009