Efficient and accurate P-value computation for Position Weight Matrices

Helene Touzet 1, 2 Jean-Stéphane Varré 1, 2
2 SEQUOIA - Sequential Learning
LIFL - Laboratoire d'Informatique Fondamentale de Lille, Inria Lille - Nord Europe
Abstract : Background Position Weight Matrices (PWMs) are probabilistic representations of signals in sequences. They are widely used to model approximate patterns in DNA or in protein sequences. The usage of PWMs needs as a prerequisite to knowing the statistical significance of a word according to its score. This is done by defining the P-value of a score, which is the probability that the background model can achieve a score larger than or equal to the observed value. This gives rise to the following problem: Given a P-value, find the corresponding score threshold. Existing methods rely on dynamic programming or probability generating functions. For many examples of PWMs, they fail to give accurate results in a reasonable amount of time. Results The contribution of this paper is two fold. First, we study the theoretical complexity of the problem, and we prove that it is NP-hard. Then, we describe a novel algorithm that solves the P-value problem efficiently. The main idea is to use a series of discretized score distributions that improves the final result step by step until some convergence criterion is met. Moreover, the algorithm is capable of calculating the exact P-value without any error, even for matrices with non-integer coefficient values. The same approach is also used to devise an accurate algorithm for the reverse problem: finding the P-value for a given score. Both methods are implemented in a software called TFM-PVALUE, that is freely available. Conclusion We have tested TFM-PVALUE on a large set of PWMs representing transcription factor binding sites. Experimental results show that it achieves better performance in terms of computational time and precision than existing tools.
Complete list of metadatas

Cited literature [15 references]  Display  Hide  Download

https://hal.inria.fr/inria-00270263
Contributor : Jean-Stéphane Varré <>
Submitted on : Friday, April 4, 2008 - 11:28:34 AM
Last modification on : Monday, September 30, 2019 - 3:36:03 PM
Long-term archiving on : Friday, September 28, 2012 - 12:15:59 PM

File

Algorithms_for_Molecular_Biolo...
Publisher files allowed on an open archive

Identifiers

  • HAL Id : inria-00270263, version 1

Collections

Citation

Helene Touzet, Jean-Stéphane Varré. Efficient and accurate P-value computation for Position Weight Matrices. Algorithms for Molecular Biology, BioMed Central, 2007, 2 (15), pp.1-12. ⟨inria-00270263⟩

Share

Metrics

Record views

235

Files downloads

253