Computing expectation values for RNA motifs using discrete convolutions.
BACKGROUND
Computational biologists use Expectation values (E-values) to estimate the number of solutions that can be expected by chance during a database scan. Here we focus on computing Expectation values for RNA motifs defined by single-strand and helix lod-score profiles with variable helix spans. Such E-values cannot be computed assuming a normal score distribution and their estimation previously required lengthy simulations.


RESULTS
We introduce discrete convolutions as an accurate and fast mean to estimate score distributions of lod-score profiles. This method provides excellent score estimations for all single-strand or helical elements tested and also applies to the combination of elements into larger, complex, motifs. Further, the estimated distributions remain accurate even when pseudocounts are introduced into the lod-score profiles. Estimated score distributions are then easily converted into E-values.


CONCLUSIONS
A good agreement was observed between computed E-values and simulations for a number of complete RNA motifs. This method is now implemented into the ERPIN software, but it can be applied as well to any search procedure based on ungapped profiles with statistically independent columns.
Lambert
André
A
CNRS UMR 6207, Université de la Méditerranée, Luminy Case 907, 13288 Marseille cedex 9, France. lambert@cpt.univ-mrs.fr
Legendre
Matthieu
M
Fontaine
Jean-Fred
JF
Gautheret
Daniel
D
D000465:Algorithms; D000818:Animals; D001483:Base Sequence; D019295:Computational Biology / Q000379:methods*; D003627:Data Interpretation, Statistical; D006801:Humans; D015233:Models, Statistical; D008962:Models, Theoretical; D009690:Nucleic Acid Conformation; D012313:RNA / Q000737:chemistry*; D038621:Regulatory Sequences, Ribonucleic Acid; D051140:Selenoproteins / Q000737:chemistry*; D016415:Sequence Alignment; D017423:Sequence Analysis, RNA / Q000379:methods*; D012689:Sequence Homology, Nucleic Acid; D012984:Software
D038621:Regulatory Sequences, Ribonucleic Acid; D051140:Selenoproteins; D012313:RNA