A partition function algorithm for interacting nucleic acid strands
1Lab for Computational Biology, School of Computing Science, Simon Fraser University, Burnaby, BC, Canada and 2Department of Bioinformatics, Institute of Computer Science, Albert-Ludwigs-Universität, Freiburg, Germany
*To whom correspondence should be addressed.
 |
Abstract |
|---|
Recent interests, such as RNA interference and antisense RNA regulation, strongly motivate the problem of predicting whether two nucleic acid strands interact.
Motivation: Regulatory non-coding RNAs (ncRNAs) such as microRNAs play an important role in gene regulation. Studies on both prokaryotic and eukaryotic cells show that such ncRNAs usually bind to their target mRNA to regulate the translation of corresponding genes. The specificity of these interactions depends on the stability of intermolecular and intramolecular base pairing. While methods like deep sequencing allow to discover an ever increasing set of ncRNAs, there are no high-throughput methods available to detect their associated targets. Hence, there is an increasing need for precise computational target prediction. In order to predict base-pairing probability of any two bases in interacting nucleic acids, it is necessary to compute the interaction partition function over the whole ensemble. The partition function is a scalar value from which various thermodynamic quantities can be derived. For example, the equilibrium concentration of each complex nucleic acid species and also the melting temperature of interacting nucleic acids can be calculated based on the partition function of the complex.
Results: We present a model for analyzing the thermodynamics of two interacting nucleic acid strands considering the most general type of interactions studied in the literature. We also present a corresponding dynamic programming algorithm that computes the partition function over (almost) all physically possible joint secondary structures formed by two interacting nucleic acids in O(n6) time. We verify the predictive power of our algorithm by computing (i) the melting temperature for interacting RNA pairs studied in the literature and (ii) the equilibrium concentration for several variants of the OxyS–fhlA complex. In both experiments, our algorithm shows high accuracy and outperforms competitors.
Availability: Software and web server is available at http://compbio.cs.sfu.ca/taverna/pirna/
Contact: cenk{at}cs.sfu.ca; backofen{at}informatik.uni-freiburg.de
Supplementary information: Supplementary data are avaliable at Bioinformatics online.
The authors wish it to be known that in their opinion, the first two authors should be regarded as joint First Authors.
CiteULike 
Connotea 
Del.icio.us What's this?
This article has been cited by other articles:
|
|
 |
|
  F. W. D. Huang, J. Qin, C. M. Reidys, and P. F. Stadler Target prediction and a statistical sampling algorithm for RNA-RNA interaction Bioinformatics, January 15, 2010; 26(2): 175 - 181. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
F. W. D. Huang, J. Qin, C. M. Reidys, and P. F. Stadler Partition function and base pairing probabilities for RNA-RNA interaction prediction Bioinformatics, October 15, 2009; 25(20): 2646 - 2654. [Abstract] [Full Text] [PDF]
|
|