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.

Inferring probabilistic Boolean networks from steady-state gene data samples

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

A new method and application for controlling the steady-state probability distributions of probabilistic Boolean networks

Relationships between probabilistic Boolean networks and dynamic Bayesian networks as models of gene regulatory networks

A novel constrained genetic algorithm-based Boolean network inference method from steady-state gene expression data

Generating Probabilistic Boolean Networks from a Prescribed Stationary Distribution

An Approximation Method for Solving the Steady-State Probability Distribution of Probabilistic Boolean Networks

Probabilistic Boolean networks: a rule-based uncertainty model for gene regulatory networks

Reviving the Two-state Markov Chain Approach (Technical Report)

Simulation Study in Probabilistic Boolean Network Models for Genetic Regulatory Networks.

Comparison of probabilistic Boolean network and dynamic Bayesian network approaches for inferring gene regulatory networks

Synthesising Executable Gene Regulatory Networks from Single-cell Gene Expression Data

A matrix perturbation method for computing the steady-state probability distributions of probabilistic Boolean networks with gene perturbations

From Boolean to probabilistic Boolean networks as models of genetic regulatory networks

LogicGep: Boolean networks inference using symbolic regression from time-series transcriptomic profiling data

On Pinning Reachability of Probabilistic Boolean Control Networks

Optimal Reconstruction of Probabilistic Boolean Networks

On Computation of the Steady-State Probability Distribution of Probabilistic Boolean Networks with Gene Perturbation.

Gene network inference from incomplete expression data: transcriptional control of hematopoietic commitment

Inferring signaling pathways with probabilistic programming

Toward Intracellular Delivery And Drug Discovery: Stochastic Logic Networks As Efficient Computational Models For Gene Regulatory Networks