An approximation method for solving the steady-state probability distribution of probabilistic Boolean networks.
BACKGROUND
Probabilistic Boolean networks (PBNs) have been proposed to model genetic regulatory interactions. The steady-state probability distribution of a PBN gives important information about the captured genetic network. The computation of the steady-state probability distribution usually includes construction of the transition probability matrix and computation of the steady-state probability distribution. The size of the transition probability matrix is 2(n)-by-2(n) where n is the number of genes in the genetic network. Therefore, the computational costs of these two steps are very expensive and it is essential to develop a fast approximation method.


RESULTS
In this article, we propose an approximation method for computing the steady-state probability distribution of a PBN based on neglecting some Boolean networks (BNs) with very small probabilities during the construction of the transition probability matrix. An error analysis of this approximation method is given and theoretical result on the distribution of BNs in a PBN with at most two Boolean functions for one gene is also presented. These give a foundation and support for the approximation method. Numerical experiments based on a genetic network are given to demonstrate the efficiency of the proposed method.
Ching
Wai-Ki
WK
Advanced Modeling and Applied Computing Laboratory, Department of Mathematics, The University of Hong Kong, Hong Kong.
Zhang
Shuqin
S
Ng
Michael K
MK
Akutsu
Tatsuya
T
D000465:Algorithms; D019295:Computational Biology / Q000379:methods*; D003198:Computer Simulation; D005786:Gene Expression Regulation; D008390:Markov Chains; D008957:Models, Genetic; D015233:Models, Statistical; D009010:Monte Carlo Method; D011336:Probability; D013269:Stochastic Processes