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.

Biologically meaningful regulatory logic enhances the convergence rate in Boolean networks and bushiness of their state transition graph

Preponderance of generalized chain functions in reconstructed Boolean models of biological networks

Relative stability of network states in Boolean network models of gene regulation in development

Choices of regulatory logic class modulate the dynamical regime in random Boolean networks

A meta-analysis of Boolean network models reveals design principles of gene regulatory networks

Leveraging developmental landscapes for model selection in Boolean gene regulatory networks

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

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

Basins of Attraction, Commitment Sets and Phenotypes of Boolean Networks

Exploring attractor bifurcations in Boolean networks

Generic Properties of Random Gene Regulatory Networks.

Mapping the attractor landscape of Boolean networks

On Pinning Reachability of Probabilistic Boolean Control Networks

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

Detecting Small Attractors of Large Boolean Networks by Function-Reduction-based Strategy.

Boolean regulatory network reconstruction using literature based knowledge with a genetic algorithm optimization method

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

Dynamics of unperturbed and noisy generalized Boolean networks

Learning Asynchronous Boolean Networks from Single-Cell Data Using Multiobjective Cooperative Genetic Programming

Boolean function metrics can assist modelers to check and choose logical rules