RSEARCH: finding homologs of single structured RNA sequences.
BACKGROUND
For many RNA molecules, secondary structure rather than primary sequence is the evolutionarily conserved feature. No programs have yet been published that allow searching a sequence database for homologs of a single RNA molecule on the basis of secondary structure.


RESULTS
We have developed a program, RSEARCH, that takes a single RNA sequence with its secondary structure and utilizes a local alignment algorithm to search a database for homologous RNAs. For this purpose, we have developed a series of base pair and single nucleotide substitution matrices for RNA sequences called RIBOSUM matrices. RSEARCH reports the statistical confidence for each hit as well as the structural alignment of the hit. We show several examples in which RSEARCH outperforms the primary sequence search programs BLAST and SSEARCH. The primary drawback of the program is that it is slow. The C code for RSEARCH is freely available from our lab's website.


CONCLUSIONS
RSEARCH outperforms primary sequence programs in finding homologs of structured RNA sequences.
D000818:Animals; D017360:Arabidopsis / Q000235:genetics; D019611:Archaeoglobus fulgidus / Q000235:genetics; D001482:Base Composition / Q000235:genetics; D019295:Computational Biology / Q000379:methods / Q000592:standards / Q000706:statistics & numerical data; D030541:Databases, Genetic; D035683:MicroRNAs / Q000737:chemistry / Q000235:genetics; D008957:Models, Genetic; D009690:Nucleic Acid Conformation; D043902:Pyrococcus horikoshii / Q000235:genetics; D012313:RNA / Q000737:chemistry* / Q000235:genetics*; D019642:RNA, Archaeal / Q000737:chemistry / Q000235:genetics; D012331:RNA, Fungal / Q000737:chemistry / Q000235:genetics; D018510:RNA, Helminth / Q000737:chemistry / Q000145:classification / Q000235:genetics; D018749:RNA, Plant / Q000737:chemistry / Q000235:genetics; D012344:RNA, Transfer, Ala / Q000737:chemistry / Q000235:genetics; D012348:RNA, Transfer, Asn / Q000737:chemistry / Q000235:genetics; D043262:Ribonuclease P / Q000737:chemistry / Q000235:genetics; D012441:Saccharomyces cerevisiae / Q000235:genetics; D012689:Sequence Homology, Nucleic Acid; D018271:Signal Recognition Particle / Q000235:genetics; D012984:Software / Q000592:standards / Q000706:statistics & numerical data; D012986:Software Validation
D035683:MicroRNAs; D019642:RNA, Archaeal; D012331:RNA, Fungal; D018510:RNA, Helminth; D018749:RNA, Plant; D012344:RNA, Transfer, Ala; D012348:RNA, Transfer, Asn; D018271:Signal Recognition Particle; D012313:RNA; D043262:Ribonuclease P