Abstract:BackgroundVarious computational models have been of interest due to their use in the modelling of gene regulatory networks (GRNs). As a logical model, probabilistic Boolean networks (PBNs) consider molecular and genetic noise, so the study of PBNs provides significant insights into the understanding of the dynamics of GRNs. This will ultimately lead to advances in developing therapeutic methods that intervene in the process of disease development and progression. The applications of PBNs, however, are hindered by the complexities involved in the computation of the state transition matrix and the steady-state distribution of a PBN. For a PBN with n genes and N Boolean networks, the complexity to compute the state transition matrix is O(nN 22n) or O(nN 2n) for a sparse matrix.ResultsThis paper presents a novel implementation of PBNs based on the notions of stochastic logic and stochastic computation. This stochastic implementation of a PBN is referred to as a stochastic Boolean network (SBN). An SBN provides an accurate and efficient simulation of a PBN without and with random gene perturbation. The state transition matrix is computed in an SBN with a complexity of O(nL 2n), where L is a factor related to the stochastic sequence length. Since the minimum sequence length required for obtaining an evaluation accuracy approximately increases in a polynomial order with the number of genes, n, and the number of Boolean networks, N, usually increases exponentially with n, L is typically smaller than N, especially in a network with a large number of genes. Hence, the computational efficiency of an SBN is primarily limited by the number of genes, but not directly by the total possible number of Boolean networks. Furthermore, a time-frame expanded SBN enables an efficient analysis of the steady-state distribution of a PBN. These findings are supported by the simulation results of a simplified p53 network, several randomly generated networks and a network inferred from a T cell immune response dataset. An SBN can also implement the function of an asynchronous PBN and is potentially useful in a hybrid approach in combination with a continuous or single-molecule level stochastic model.ConclusionsStochastic Boolean networks (SBNs) are proposed as an efficient approach to modelling gene regulatory networks (GRNs). The SBN approach is able to recover biologically-proven regulatory behaviours, such as the oscillatory dynamics of the p53-Mdm2 network and the dynamic attractors in a T cell immune response network. The proposed approach can further predict the network dynamics when the genes are under perturbation, thus providing biologically meaningful insights for a better understanding of the dynamics of GRNs. The algorithms and methods described in this paper have been implemented in Matlab packages, which are attached as Additional files.

A Random Boolean Network shifted toward a critical point

Dynamical Interplay Between Random Boolean Networks and Awareness Propagation

Truth-Content and Phase Transitions of Random Boolean Networks with Generic Logics

Dynamics Of Boolean Networks With Small-World Topology

Highly-sensitive measure of complexity captures boolean networks regimes and temporal order more optimally

Fixed-points in Random Boolean Networks: The impact of parallelism in the scale-free topology case

Dynamical regimes in non-ergodic random Boolean networks

Random Pseudofractal Networks with Competition

Dynamical PN Evolution with Magnetic Fields

Intrinsic noise and deviations from criticality in Boolean gene-regulatory networks

Generic Properties of Random Gene Regulatory Networks.

The Stabilization of Random Boolean Networks Through Edge Immunization

Stochastic Boolean networks: An efficient approach to modeling gene regulatory networks

On Pinning Reachability of Probabilistic Boolean Control Networks

Random Boolean Network Models and the Yeast Transcriptional Network

Exploring biological network structure with clustered random networks

Criticality in Random Threshold Networks: Annealed Approximation and Beyond

Dynamical Properties of a Boolean Model of Gene Regulatory Network with Memory

Binary threshold networks as a natural null model for biological networks

Biological Regulatory Networks are Minimally Frustrated

Critical excitation-inhibition balance in dense neural networks