A partition function algorithm for interacting nucleic acid strands

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@cs.sfu.ca; backofen@informatik.uni-freiburg.de Supplementary information: Supplementary data are avaliable at Bioinformatics online.

[1]  Phillip D. Zamore,et al.  RNA Interference , 2000, Science.

[2]  Klas Flärdh,et al.  Antisense RNAs everywhere? , 2002, Trends in genetics : TIG.

[3]  Niles A. Pierce,et al.  A partition function algorithm for nucleic acid secondary structure including pseudoknots , 2003, J. Comput. Chem..

[4]  Erik Winfree,et al.  Thermodynamic Analysis of Interacting Nucleic Acid Strands , 2007, SIAM Rev..

[5]  F. Simmel,et al.  DNA nanodevices. , 2005, Small.

[6]  J. Steitz,et al.  The expanding universe of noncoding RNAs. , 2006, Cold Spring Harbor symposia on quantitative biology.

[7]  E Rivas,et al.  A dynamic programming algorithm for RNA structure prediction including pseudoknots. , 1998, Journal of molecular biology.

[8]  Phillip D. Zamore,et al.  Ribo-gnome: The Big World of Small RNAs , 2005, Science.

[9]  A. Condon,et al.  Secondary structure prediction of interacting RNA molecules. , 2005, Journal of molecular biology.

[10]  D. Turner,et al.  Experimentally derived nearest-neighbor parameters for the stability of RNA three- and four-way multibranch loops. , 2002, Biochemistry.

[11]  G. Storz.  An Expanding Universe of Noncoding RNAs , 2002, Science.

[12]  S. Gottesman.  Micros for microbes: non-coding regulatory RNAs in bacteria. , 2005, Trends in genetics : TIG.

[13]  Philip S Lukeman,et al.  Nucleic acid nanostructures: bottom-up control of geometry on the nanoscale , 2005, Reports on progress in physics. Physical Society.

[14]  N. Seeman.  From genes to machines: DNA nanomechanical devices. , 2005, Trends in biochemical sciences.

[15]  Robert M. Dirks,et al.  An autonomous polymerization motor powered by DNA hybridization , 2007, Nature Nanotechnology.

[16]  L. Argaman,et al.  fhlA repression by OxyS RNA: kissing complex formation at two sites results in a stable antisense-target RNA complex. , 2000, Journal of molecular biology.

[17]  Harry M. T. Choi,et al.  Programming DNA Tube Circumferences , 2008, Science.

[18]  E. Wagner,et al.  Bulged-out nucleotides in an antisense RNA are required for rapid target RNA binding in vitro and inhibition in vivo. , 1995, Nucleic acids research.

[19]  Tatsuya Akutsu,et al.  A grammatical approach to RNA-RNA interaction prediction , 2009, Pattern Recognit..

[20]  D. Turner,et al.  Thermodynamic parameters for an expanded nearest-neighbor model for formation of RNA duplexes with Watson-Crick base pairs. , 1998, Biochemistry.

[21]  J. Sabina,et al.  Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. , 1999, Journal of molecular biology.

[22]  M. Waterman,et al.  RNA secondary structure: a complete mathematical analysis , 1978 .

[23]  I. Tinoco,et al.  Absorbance melting curves of RNA. , 1989, Methods in enzymology.

[24]  J. McCaskill.  The equilibrium partition function and base pair binding probabilities for RNA secondary structure , 1990, Biopolymers.

[25]  Michael Zuker,et al.  Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information , 1981, Nucleic Acids Res..

[26]  Peter F. Stadler,et al.  Partition function and base pairing probabilities of RNA heterodimers , 2006, Algorithms for Molecular Biology.

[27]  Michael Zuker,et al.  UNAFold: software for nucleic acid folding and hybridization. , 2008, Methods in molecular biology.

[28]  Jerrold R. Griggs,et al.  Algorithms for Loop Matchings , 1978 .

[29]  D. Mathews.  Using an RNA secondary structure partition function to determine confidence in base pairs predicted by free energy minimization. , 2004, RNA.

[30]  G. Hannon.  RNA interference : RNA , 2002 .

[31]  Kaizhong Zhang,et al.  RNA-RNA Interaction Prediction and Antisense RNA Target Search , 2006, J. Comput. Biol..

[32]  Peter F. Stadler,et al.  Thermodynamics of RNA-RNA Binding , 2006, German Conference on Bioinformatics.

[33]  M. Zuker,et al.  Prediction of hybridization and melting for double-stranded nucleic acids. , 2004, Biophysical journal.

[34]  D. Haar,et al.  Statistical Physics , 1971, Nature.

[35]  Song Cao,et al.  Predicting RNA pseudoknot folding thermodynamics , 2006, Nucleic acids research.

[36]  R. Giegerich,et al.  Fast and effective prediction of microRNA/target duplexes. , 2004, RNA.

[37]  S. P. Walton,et al.  Thermodynamic and kinetic characterization of antisense oligodeoxynucleotide binding to a structured mRNA. , 2002, Biophysical journal.

[38]  J. M. Diamond,et al.  Thermodynamics of three-way multibranch loops in RNA. , 2001, Biochemistry.

[39]  D. Pervouchine.  IRIS: intermolecular RNA interaction search. , 2004, Genome informatics. International Conference on Genome Informatics.

[40]  D. Bartel.  MicroRNAs Genomics, Biogenesis, Mechanism, and Function , 2004, Cell.

[41]  Rolf Backofen,et al.  IntaRNA: efficient prediction of bacterial sRNA targets incorporating target site accessibility and seed regions , 2008, Bioinform..

[42]  S. Brantl,et al.  Antisense-RNA regulation and RNA interference. , 2002, Biochimica et biophysica acta.