Simulation study in Probabilistic Boolean Network models for genetic regulatory networks

Probabilistic Boolean Network (PBN) is widely used to model genetic regulatory networks. Evolution of the PBN is according to the transition probability matrix. Steady-state (long-run behaviour) analysis is a key aspect in studying the dynamics of genetic regulatory networks. In this paper, an efficient method to construct the sparse transition probability matrix is proposed, and the power method based on the sparse matrix-vector multiplication is applied to compute the steady-state probability distribution. Such methods provide a tool for us to study the sensitivity of the steady-state distribution to the influence of input genes, gene connections and Boolean networks. Simulation results based on a real network are given to illustrate the method and to demonstrate the steady-state analysis.