Boolean network

A Boolean network is a particular kind of sequential dynamical system, where time and states are discrete, i.e. both the set of variables and the set of states in the time series each have a bijection onto an integer series.

A random Boolean network (RBN) is one that is randomly selected from the set of all possible boolean networks of a particular size, N. One then can study statistically, how the expected properties of such networks depend on various statistical properties of the ensemble of all possible networks. For example, one may study how the RBN behavior changes as the average connectivity is changed.

The first Boolean networks were proposed by Stuart A. Kauffman in 1969, as random models of genetic regulatory networks[5] but their mathematical understanding only started in the 2000s.[6][7]

Attractors

Since a Boolean network has only 2^N possible states, a trajectory will sooner or later reach a previously visited state, and thus, since the dynamics are deterministic, the trajectory will fall into a steady state or cycle called an attractor (though in the broader field of dynamical systems a cycle is only an attractor if perturbations from it lead back to it). If the attractor has only a single state it is called a point attractor, and if the attractor consists of more than one state it is called a cycle attractor. The set of states that lead to an attractor is called the basin of the attractor. States which occur only at the beginning of trajectories (no trajectories lead to them), are called garden-of-Eden states[8] and the dynamics of the network flow from these states towards attractors. The time it takes to reach an attractor is called transient time.[4]

With growing computer power and increasing understanding of the seemingly simple model, different authors gave different estimates for the mean number and length of the attractors, here a brief summary of key publications.[9]

Author	Year	Mean attractor length	Mean attractor number	comment
Kauffman [5]	1969	${\displaystyle \langle A\rangle \sim {\sqrt {N}}}$	${\displaystyle \langle \nu \rangle \sim {\sqrt {N}}}$
Bastolla/ Parisi[10]	1998	faster than a power law, ${\displaystyle \langle A\rangle >N^{x}\forall x}$	faster than a power law, ${\displaystyle \langle \nu \rangle >N^{x}\forall x}$	first numerical evidences
Bilke/ Sjunnesson[11]	2002		linear with system size, ${\displaystyle \langle \nu \rangle \sim N}$
Socolar/Kauffman[12]	2003		faster than linear, ${\displaystyle \langle \nu \rangle >N^{x}}$ with ${\displaystyle x>1}$
Samuelsson/Troein[13]	2003		superpolynomial growth, ${\displaystyle \langle \nu \rangle >N^{x}\forall x}$	mathematical proof
Mihaljev/Drossel[14]	2005	faster than a power law, ${\displaystyle \langle A\rangle >N^{x}\forall x}$	faster than a power law, ${\displaystyle \langle \nu \rangle >N^{x}\forall x}$

In dynamical systems theory, the structure and length of the attractors of a network corresponds to the dynamic phase of the network. The stability of Boolean networks depends on the connections of their nodes. A Boolean network can exhibit stable, critical or chaotic behavior. This phenomenon is governed by a critical value of the average number of connections of nodes (${\displaystyle K_{c}}$), and can be characterized by the Hamming distance as distance measure. In the unstable regime, the distance between two initially close states on average grows exponentially in time, while in the stable regime it decreases exponentially. In this, with "initially close states" one means that the Hamming distance is small compared with the number of nodes (${\displaystyle N}$) in the network.

For N-K-model[15] the network is stable if ${\displaystyle K<K_{c}}$, critical if ${\displaystyle K=K_{c}}$, and unstable if ${\displaystyle K>K_{c}}$.

The state of a given node ${\displaystyle n_{i}}$ is updated according to its truth table, whose outputs are randomly populated. ${\displaystyle p_{i}}$ denotes the probability of assigning an off output to a given series of input signals.

If ${\displaystyle p_{i}=p=const.}$ for every node, the transition between the stable and chaotic range depends on ${\displaystyle p}$. According to Bernard Derrida and Yves Pomeau[16] , the critical value of the average number of connections is ${\displaystyle K_{c}=1/[2p(1-p)]}$.

If ${\displaystyle K}$ is not constant, and there is no correlation between the in-degrees and out-degrees, the conditions of stability is determined by ${\displaystyle \langle K^{in}\rangle }$[17][18][19] The network is stable if ${\displaystyle \langle K^{in}\rangle <K_{c}}$, critical if ${\displaystyle \langle K^{in}\rangle =K_{c}}$, and unstable if ${\displaystyle \langle K^{in}\rangle >K_{c}}$.

The conditions of stability are the same in the case of networks with scale-free topology where the in-and out-degree distribution is a power-law distribution: ${\displaystyle P(K)\propto K^{-\gamma }}$, and ${\displaystyle \langle K^{in}\rangle =\langle K^{out}\rangle }$, since every out-link from a node is an in-link to another.[20]

Sensitivity shows the probability that the output of the Boolean function of a given node changes if its input changes. For random Boolean networks, ${\displaystyle q_{i}=2p_{i}(1-p_{i})}$. In the general case, stability of the network is governed by the largest eigenvalue ${\displaystyle \lambda _{Q}}$ of matrix ${\displaystyle Q}$, where ${\displaystyle Q_{ij}=q_{i}A_{ij}}$, and ${\displaystyle A}$ is the adjacency matrix of the network.[21] The network is stable if ${\displaystyle \lambda _{Q}<1}$, critical if ${\displaystyle \lambda _{Q}=1}$, unstable if ${\displaystyle \lambda _{Q}>1}$.

• NK model

• Dubrova, E., Teslenko, M., Martinelli, A., (2005). *Kauffman Networks: Analysis and Applications, in "Proceedings of International Conference on Computer-Aided Design", pages 479-484.

• Analysis of Dynamic Algebraic Models (ADAM) v1.1
• bioasp/bonesis: Synthesis of Most Permissive Boolean Networks from network architecture and dynamical properties
• CoLoMoTo (Consortium for Logical Models and Tools)
• DDLab
• NetBuilder Boolean Networks Simulator
• Open Source Boolean Network Simulator
• JavaScript Kauffman Network
• Probabilistic Boolean Networks (PBN)
• RBNLab
• A SAT-based tool for computing attractors in Boolean Networks