RNA-RNA Interaction Prediction and Antisense RNA Target Search

Recent studies demonstrating the existence of special noncoding "antisense" RNAs used in post transcriptional gene regulation have received considerable attention. These RNAs are synthesized naturally to control gene expression in C. elegans, Drosophila, and other organisms; they are known to regulate plasmid copy numbers in E. coli as well. Small RNAs have also been artificially constructed to knock out genes of interest in humans and other organisms for the purpose of finding out more about their functions. Although there are a number of algorithms for predicting the secondary structure of a single RNA molecule, no such algorithm exists for reliably predicting the joint secondary structure of two interacting RNA molecules or measuring the stability of such a joint structure. In this paper, we describe the RNA-RNA interaction prediction (RIP) problem between an antisense RNA and its target mRNA and develop efficient algorithms to solve it. Our algorithms minimize the joint free energy between the two RNA molecules under a number of energy models with growing complexity. Because the computational resources needed by our most accurate approach is prohibitive for long RNA molecules, we also describe how to speed up our techniques through a number of heuristic approaches while experimentally maintaining the original accuracy. Equipped with this fast approach, we apply our method to discover targets for any given antisense RNA in the associated genome sequence.

[1]  I. Tinoco,et al.  A retroviral RNA kissing complex containing only two G.C base pairs. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

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

[3]  M. Zuker On finding all suboptimal foldings of an RNA molecule. , 1989, Science.

[4]  L. Lim,et al.  An Abundant Class of Tiny RNAs with Probable Regulatory Roles in Caenorhabditis elegans , 2001, Science.

[5]  Kaizhong Zhang,et al.  A new algorithm for computing similarity between RNA structures , 2001, Inf. Sci..

[6]  E. Moss,et al.  MicroRNAs: Hidden in the Genome , 2002, Current Biology.

[7]  Anne Condon,et al.  RNAsoft: a suite of RNA secondary structure prediction and design software tools , 2003, Nucleic Acids Res..

[8]  Michael T. McManus,et al.  Gene silencing in mammals by small interfering RNAs , 2002, Nature Reviews Genetics.

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

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

[11]  E Westhof,et al.  Progression of a loop–loop complex to a four‐way junction is crucial for the activity of a regulatory antisense RNA , 2000, The EMBO journal.

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

[13]  Bin Ma,et al.  Computing similarity between RNA structures , 1999, Theor. Comput. Sci..

[14]  Christian N. S. Pedersen,et al.  Fast evaluation of internal loops in RNA secondary structure prediction , 1999, Bioinform..

[15]  T. Tuschl,et al.  Identification of Novel Genes Coding for Small Expressed RNAs , 2001, Science.

[16]  Anne Condon,et al.  Classifying RNA pseudoknotted structures , 2004, Theor. Comput. Sci..

[17]  D. Higgins,et al.  RAGA: RNA sequence alignment by genetic algorithm. , 1997, Nucleic acids research.

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

[19]  R. Nussinov,et al.  Fast algorithm for predicting the secondary structure of single-stranded RNA. , 1980, Proceedings of the National Academy of Sciences of the United States of America.

[20]  E. Moss,et al.  RNA interference: It's a small RNA world , 2001, Current Biology.

[21]  J. SantaLucia,et al.  Nearest-neighbor thermodynamics and NMR of DNA sequences with internal A.A, C.C, G.G, and T.T mismatches. , 1999, Biochemistry.

[22]  B. Reinhart,et al.  The 21-nucleotide let-7 RNA regulates developmental timing in Caenorhabditis elegans , 2000, Nature.

[23]  Tatsuya Akutsu,et al.  Dynamic programming algorithms for RNA secondary structure prediction with pseudoknots , 2000, Discret. Appl. Math..