title: Computing expectation values for RNA motifs using discrete convolutions.
abstract: 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.
authors.affiliation: CNRS UMR 6207, Université de la Méditerranée, Luminy Case 907, 13288 Marseille cedex 9, France. lambert@cpt.univ-mrs.fr
authors.forename: André
authors.initials: A
authors.lastname: Lambert
authors.forename: Matthieu
authors.initials: M
authors.lastname: Legendre
authors.forename: Jean-Fred
authors.initials: JF
authors.lastname: Fontaine
authors.forename: Daniel
authors.initials: D
authors.lastname: Gautheret
mesh_terms: 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
chemical_list: D038621:Regulatory Sequences, Ribonucleic Acid; D051140:Selenoproteins; D012313:RNA