Next Article in Journal
Minirobots Moving at Different Partial Speeds
Next Article in Special Issue
On the Periodic Structure of Parallel Dynamical Systems on Generalized Independent Boolean Functions
Previous Article in Journal
Implicit-Explicit Methods for a Convection-Diffusion-Reaction Model of the Propagation of Forest Fires
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Lyapunov Exponent of Monotone Boolean Networks

Institute for Systems Biology, Seattle, WA 98103, USA
Mathematics 2020, 8(6), 1035; https://doi.org/10.3390/math8061035
Submission received: 8 June 2020 / Revised: 20 June 2020 / Accepted: 23 June 2020 / Published: 24 June 2020
(This article belongs to the Special Issue Boolean Networks Models in Science and Engineering)

Abstract

:
Boolean networks are discrete dynamical systems comprised of coupled Boolean functions. An important parameter that characterizes such systems is the Lyapunov exponent, which measures the state stability of the system to small perturbations. We consider networks comprised of monotone Boolean functions and derive asymptotic formulas for the Lyapunov exponent of almost all monotone Boolean networks. The formulas are different depending on whether the number of variables of the constituent Boolean functions, or equivalently, the connectivity of the Boolean network, is even or odd.

1. Introduction

Boolean networks are complex dynamical systems that were proposed as models of genetic regulatory networks [1,2] and have since then been used to model a range of complex phenomena. Random Boolean networks (RBNs) are ensembles of randomly generated Boolean networks with random topology and random updating functions. It is known that RBNs undergo a phase transition. In the thermodynamic limit, meaning that the number of nodes goes to infinity, the phase transition curve for RBNs is given by λ = log ( 2 K p 1 p ) = 0 , where p is the so-called bias of the random Boolean functions, which is the probability that the function takes on the value 1, and K is the connectivity of the network, or equivalently, the number of inputs to each Boolean function [3].
Under a synchronous updating scheme, whereby all Boolean functions get updated simultaneously at each time step, this phase transition curve separates two qualitatively distinct dynamical regimes. Below the curve, when λ < 0 , infinitesimally small perturbations to the state of the network decay, while above the curve, when λ > 0 , small perturbations spread throughout the network. Thus, the two regimes are often referred to as ordered and chaotic regimes, respectively. The quantity λ coincides with the Lyapunov exponent [4,5]. Networks that are known to operate in the dynamically critical regime, meaning that their Lyapunov exponent is 0, are widely known to have many optimal properties and are thought to be a hallmark of living systems [2,6]. In the context of Boolean networks, these properties include: the ability to balance resilience to random mutations with emergence of new phenotypes for adapting to new environmental challenges [7], maximization of the associative memory of the network [8], maximization of the diversity in structure-dynamics relationships [9], maximization of communication among nodes [10], maximization of so-called set complexity of the dynamics of the network [11], emergence of diversity in a spatial arrangement of Boolean networks representing cells in a tissue [12], and optimum network learning and generalization when the Boolean network has input and output signals [13,14].
We consider monotone Boolean networks in which all updating rules belong to the class of monotone Boolean functions. This class of functions is one of the most widely studied classes of Boolean functions [15]. Monotone Boolean networks have been primarily studied in the asynchronous updating scheme setting, whereby only one node is updated at a time. Some work has focused on long-term dynamics, such as fixed points and limit cycles [16,17,18]. It is also known that certain classes of fully asynchronous Boolean networks can be simulated by monotone Boolean networks [19]. However, in the context of the Lyapunov exponent, asynchronous updating schemes appear to be less relevant than synchronous updating schemes [20]. Our main results are asymptotic formulas, depending on whether n is even or odd, for the so-called expected average sensitivity of a monotone Boolean function, first given in [21]. In light of the results in [5], the logarithm of the expected average sensitivities, s ^ f , in Theorems 4 and 5, can be directly interpreted as the Lyapunov exponent λ .
In the remainder, we will use n, rather than the conventional K in the Boolean network literature, to denote the number of variables of the Boolean functions. Our results concerning almost all monotone Boolean networks can be understood in a probabilistic manner, meaning that the asymptotic formulas are valid with probability almost 1 if a monotone Boolean network is chosen at random from the set of all such networks. As the formulas are asymptotic, they should not be interpreted for small n. At the same time, they are quite accurate even for n = 7 , which follows from the results originally published in [22]. Absent additional constraints on the Boolean functions, such as the classes of canalizing functions [5,23,24,25] or Post classes [26], random monotone Boolean networks quickly enter the disordered regime relative to n, but slower than RBNs. For example, for n = 4 , monotone Boolean networks have the expected average sensitivity s ^ f = 1.125 , which is already slightly in the chaotic regime ( λ > 0 ). For n = 20 , the expected average sensitivity is s f = 3.2 , whereas for random Boolean networks, it is s f = 10 (we can also see this if we plug n = 20 and p = 1 / 2 into 2 n p 1 p ).

2. Definitions and Preliminaries

Let f : 0 , 1 n 0 , 1 be a Boolean function of n variables x 1 , , x n . Let
f ( x ˜ ) / x j = f ( x ˜ ( j , 0 ) ) f ( x ˜ ( j , 1 ) )
be the partial derivative of f with respect to x j , where ⊕ is addition modulo 2 (exclusive OR) and x ˜ j , k = x 1 , , x j 1 , k , x j + 1 , x n , k = 0 , 1 . Clearly, the partial derivative is a Boolean function itself that specifies whether a change in the jth input causes a change in the original function f. The activity of variable x j in function f can be defined as
α j f = 1 2 n x ˜ 0 , 1 n f x ˜ / x j .
Note that although the vector x ˜ consists of n components (variables), the jth variable is fictitious in f ( x ˜ ) / x j . A variable x j is fictitious in f if f ( x ˜ ( j , 0 ) ) = f ( x ˜ ( j , 1 ) ) for all x ˜ ( j , 0 ) and x ˜ ( j , 1 ) . For a n-variable Boolean function f, we can form its activity vector α ˜ f = [ α 1 f , , α n f ] . It is easy to see that 0 α j f 1 , for any j = 1 , , n . In fact, we can consider α j f to be a probability that toggling the jth input bit changes the function value, when the input vectors x ˜ are distributed uniformly over { 0 , 1 } n . Since we’re in the binary setting, the activity is also the expectation of the partial derivative with respect to the uniform distribution: α j f = E [ f ( x ˜ ) / x j ] . The activity of a fictitious variable x j is α j f = 0 . Under an arbitrary distribution, α j f is referred to as the influence of variable x j on the function f [27]. The influence of variables was used in the context of genetic regulatory network modeling in [28].
Another important quantity is the sensitivity of a Boolean function f, which measures how sensitive the output of the function is to changes in the inputs. The sensitivity s f x ˜ of f on vector x ˜ is defined as the number of Hamming neighbors of x ˜ on which the function value is different than on x ˜ (two vectors are Hamming neighbors if they differ in only one component). That is,
s f x ˜ = i 1 , , n : f x ˜ e i f x ˜ = i = 1 n χ f x ˜ e i f x ˜ ,
where e i is the unit vector with 1 in the ith position and 0s everywhere else and χ A is an indicator function that is equal to 1 if and only if A is true. The average sensitivity s f is defined by taking the expectation of s f x ˜ with respect to the distribution of x ˜ . It is easy to see that under the uniform distribution, the average sensitivity is equal to the sum of the activities:
s f = E s f x ˜ = i = 1 n E χ f x ˜ e i f x ˜ = i = 1 n α i f .
Therefore, s f is a number between 0 and n.
The average sensitivity has been studied intensively by a number of authors [5,29,30,31,32,33,34,35,36,37]. For example, it was shown by Friedgut [33] that if the average sensitivity of f is k then f can be approximated by a function depending on only c k variables where c is a constant depending only on the accuracy of the approximation, but not on n. Shi [34] showed that the average sensitivity can serve as a lower bound of quantum query complexity. Average sensitivity was used to characterize the noise sensitivity of monotone Boolean functions by Mossel and O’Donnell [35]. Zhang [36] gives lower and upper bounds of the average sensitivity of a monotone Boolean function. The upper bound is asymptotic to n , which has been shown by Bshouty and Tamon [37]. Shmulevich and Kauffman [5] have shown that the average sensitivity determines the critical phase transition curve in random Boolean networks and thus coincides with the Lyapunov exponent as λ = log s f . We now turn to monotone Boolean functions.

2.1. Monotone Boolean Functions

Let α ˜ = α 1 , , α n and β ˜ = β 1 , , β n be two different n-element binary vectors. We say that α ˜ precedes β ˜ , denoted as α ˜ β ˜ , if α i β i for every i, 1 i n . If α ˜ β ˜ and β ˜ α ˜ , then α ˜ and β ˜ are said to be incomparable. Relative to the predicate ≺, the set of all binary vectors of a given length is a partially ordered set. A Boolean function f x 1 , , x n is called monotone if for any two vectors α ˜ and β ˜ such that α ˜ β ˜ , we have f α ˜ f ( β ˜ ) .
We denote by M n the set of all monotone Boolean functions of n variables. Let E n denote the Boolean n-cube, that is, a graph with 2 n vertices each of which is labeled by an n-element binary vector. Two vertices α ˜ = α 1 , , α n and β ˜ = β 1 , , β n are connected by an edge if and only if the Hamming distance ρ ( α ˜ , β ˜ ) = i = 1 n α i β i = 1 . The set of those vectors from E n in which there are exactly k units, 0 k n , is called the kth layer of E n and is denoted by E n , k .
A vector α ˜ E n is called a minimal one or minimal unit of monotone Boolean function f ( x 1 , , x n ) if f ( α ˜ ) = 1 and f ( β ˜ ) = 0 for any β ˜ α ˜ . A vector α ˜ E n is called an maximal zero of monotone Boolean function f ( x 1 , , x n ) if f ( α ˜ ) = 0 and f ( β ˜ ) = 1 for any β ˜ α ˜ . The minimal ones correspond directly to the terms in the minimal disjunctive normal form (DNF) representation of the monotone Boolean function. In [38], asymptotic formulae for the number of monotone Boolean functions of n variables with a most probable number of minimal ones were derived, confirming the conjecture in [39] that the number of monotone Boolean functions relative to the number of minimal ones asymptotically follows a normal distribution, with the assumption of all monotone Boolean functions being equiprobable.

2.2. The Structure of Special Monotone Boolean Functions

We now briefly review some known results concerning the structure of so-called special monotone Boolean functions. Let M 0 n denote the set of functions in M n possessing the following properties. If n is even, then M 0 n contains only functions f M n such that all minimal ones of f are situated in E n , n / 2 1 , E n , n / 2 , and E n , n / 2 + 1 while function f is equal to 1 on all vectors in E n , n / 2 + 2 , , E n , n . For odd n, M 0 n contains only functions f M n such that all minimal ones of f are situated in either E n , n 3 / 2 , E n , n 1 / 2 , and E n , n + 1 / 2 or E n , n 1 / 2 , E n , n + 1 / 2 , and E n , n + 3 / 2 . In the first case, f α ˜ = 1 for all α ˜ in E n , n + 3 / 2 , , E n , n while in the second case, f α ˜ = 1 for all α ˜ in E n , n + 5 / 2 , , E n , n .
Then, as shown in [22],
lim n M 0 n M n = 1 ,
which we denote by M 0 n M n . In [40], asymptotic formulae for the number of special functions from M 0 n were established and subsequently used to characterize statistical properties of a popular class of nonlinear digital filters called stack filters [41]. The set of these special functions is denoted by M 0 1 n and, depending on whether n is even or odd, is defined differently. While we shall omit the rather lengthy definitions of special functions, the result from [40] that will be important to us is that M 0 1 n M n . In other words, almost all monotone Boolean functions are special. We shall also need the following results.
Let us start with the case of even n. Let
r 0 = r 0 n = v 0 = v 0 n = n n / 2 1 2 n / 2 1 , z 0 = 1 2 n n / 2 .
Let M 0 1 n , r , z , v denote the set of functions f M 0 1 n such that f has r minimal ones in E n , n / 2 1 , v maximal zeros in E n , n / 2 + 1 , and f is equal to 1 on z vertices in E n , n / 2 . In [40], the following result was proved.
Theorem 1.
Let n be even,
r = r 0 + k , z = z 0 + u , v = v 0 + t ,
where r 0 , z 0 , v 0 are defined in (6). Then, for any k, t, and u such that k n 2 n / 4 , t n 2 n / 4 , u n 2 n / 2 ,
M 0 1 n , r , z , v 2 n + 1 π 3 n n / 2 3 M n × exp 2 n / 2 n n / 2 1 k 2 + t 2 2 u 2 n n / 2 .
For any odd n, we use the parameters r 1 , z 1 , v 1 which are given by
r 1 = r 1 n = n n 3 / 2 2 n + 3 / 2 , v 1 = v 1 n = n n + 1 / 2 2 n + 1 / 2 ,
z 1 = 1 2 n n 1 / 2 + r 1 n + 3 / 2 v 1 n + 1 / 2
and parameters r 2 , z 2 , v 2 , which are given by
r 2 = r 2 n = n n 1 / 2 2 n + 1 / 2 , v 2 = v 2 n = n n + 3 / 2 2 n + 3 / 2 ,
z 2 = 1 2 n n + 1 / 2 + r 2 n 1 / 2 v 2 n + 3 / 2 .
Let M 0 , 1 1 n , r , z , v denote the set of functions f M 0 1 n such that f has r minimal ones in E n , n 3 / 2 , v maximal zeros in E n , n + 1 / 2 , and f is equal to 1 on z vertices in E n , n 1 / 2 . Similarly, let M 0 , 2 1 n , r , z , v denote the set of functions f M 0 1 n such that f has r minimal ones in E n , n 1 / 2 , v maximal zeros in E n , n + 3 / 2 , and f is equal to 1 on z vertices in E n , n + 1 / 2 . Then, in [40], the following two Theorems were proved.
Theorem 2.
Let n be odd,
r = r 1 + k , z = z 1 + u , v = v 1 + t ,
where r 1 , z 1 , v 1 are defined in (8) and (9). Then, for any k, t, and u such that k n 2 n / 4 , t n 2 n / 4 , u n 2 n / 2 ,
M 0 , 1 1 n , r , z , v 1 2 2 n + 1 π 3 n n 1 / 2 3 M n × exp 2 n + 1 / 2 n n 3 / 2 k 2 2 n 1 / 2 n n + 1 / 2 t 2 2 u 2 n n 1 / 2 .
Theorem 3.
Let n be odd,
r = r 2 + k , z = z 2 + u , v = v 2 + t ,
where r 2 , z 2 , v 2 are defined in (10) and (11). Then, for any k, t, and u such that k n 2 n / 4 , t n 2 n / 4 , u n 2 n / 2 ,
M 0 , 2 1 n , r , z , v 1 2 2 n + 1 π 3 n n 1 / 2 3 M n × exp 2 n 1 / 2 n n 1 / 2 k 2 2 n + 1 / 2 n n + 3 / 2 t 2 2 u 2 n n + 1 / 2 .
Prior to presenting the main results, we summarize some of the notation in Table 1.

3. Main Results

Since M 0 n M n , we can focus our attention on functions in M 0 n and derive the average sensitivity of a typical function from M 0 n . By ‘typical’ we mean the most probable Boolean function relative to the parameters k, t, and u in Theorems 1–3. It can easily be seen that the most probable special Boolean functions will have k = t = u = 0 . This will imply, to take the n-even case as an example, that the most probable function f has r 0 minimal ones in E n , n / 2 1 , v 0 maximal zeros in E n , n / 2 + 1 , and f is equal to 1 on z 0 vertices in E n , n / 2 , where r 0 , v 0 , and z 0 are given in (6). Our proofs are thus based on the derivation of the average sensitivity of such a function. Whenever we make probabilistic assertions using words such as ‘most probable’ or ‘typical’ or talk about expectations, we are implicitly endowing the set M 0 1 n , r , z , v with a uniform probability distribution for fixed parameters n , r , z , v . This should not be confused with the Gaussian-like distribution of M 0 1 n , r , z , v relative to its parameters n , r , z , v [38]. We will also omit the floor notation · as the results are asymptotic.
Theorem 4.
Let n be even and let f M 0 1 n be a typical monotone Boolean function. Then, the expected average sensitivity s ^ f = E [ s f ] of f is
s ^ f n 2 n n n / 2 1 2 n / 2 1 + 1 .
Proof. 
We will proceed by first focusing on determining the activity of an arbitrary variable x j of a typical function f . By simple symmetry arguments, if we were to sample randomly from the set M n of monotone Boolean functions, the expected activities would be equal for all the variables. It will follow by (4) that the expected average sensitivity will be equal to n multiplied by the expected activity. Since the function f is such that its minimal ones are situated in E n , n / 2 1 , E n , n / 2 , and E n , n / 2 + 1 while it is equal to 1 on all vectors in E n , n / 2 + 2 , , E n , n , the only non-trivial behavior occurs between the layers E n , n / 2 2 and E n , n / 2 + 2 .
Let us consider the minimal ones, and hence all of the ones, on E n , n / 2 1 . Since we are considering variable x j , half of these minimal ones will have x j = 0 (i.e., x ˜ j , 0 ) and the other half will have x j = 1 (i.e., x ˜ j , 1 ). It is easy to see that if x ˜ j , 0 E n , n / 2 1 is a minimal one, then by monotonicity, f x ˜ j , 1 = 1 . Consequently, the Hamming neighbors x ˜ j , 0 and x ˜ j , 1 contribute nothing to the sum in (2). On the other hand, if x ˜ j , 1 E n , n / 2 1 is a minimal one, then f x ˜ / x j = 1 , since f x ˜ = 0 for all x ˜ E n , n / 2 2 . The most probable number of such minimal ones on E n , n / 2 1 contributing to the sum in (2) is thus equal to
1 2 n n / 2 1 2 n / 2 1 .
The number of zeros on E n , n / 2 1 is equal to
n n / 2 1 r 0 = n n / 2 1 1 2 n / 2 1 .
As above, half of these will have x j = 0 and half will have x j = 1 . We need not consider vectors x ˜ j , 1 E n , n / 2 1 , since f x ˜ = 0 for all x ˜ E n , n / 2 2 . However, we should consider the number of ones situated on the middle layer E n , n / 2 . The middle layer contains
1 2 n n / 2
ones and an equal number of zeros. Thus, half of the vectors x ˜ j , 1 E n , n / 2 will be ones and the other half will be zeros. In total, the number of vectors x ˜ E n , n / 2 1 such that f x ˜ = 0 , x j = 0 , and f x ˜ j , 1 = 1 , is equal to
1 4 n n / 2 1 1 2 n / 2 1 .
We have now examined all partial derivatives above and below the layer E n , n / 2 1 .
Let us now jump to layer E n , n / 2 + 1 , as it will be similar by duality considerations. The number of maximal zeros on that layer is equal to v 0 = r 0 . Half of these will have x j = 1 and thus f x ˜ / x j = 0 due to monotonicity. The other half will have x j = 0 and since f x ˜ = 1 for all x ˜ E n , n / 2 + 2 , the same total as in (14) will result. Similarly, the number of ones on E n , n / 2 + 1 is the same as in (15). We are only concerned with x ˜ j , 1 E n , n / 2 + 1 , such that f x ˜ j , 1 = 1 and f x ˜ j , 0 = 0 . As above, because the middle layer contains the same number of ones and zeros, the total number of such pairs of vectors is the same as in (17). Thus, having accounted for all partial derivatives above and below E n , n / 2 + 1 and having convinced ourselves that their contribution to the overall activity of variable x j is the same as in (14) and (17), we can multiply (14) and (17) by 2 and add them together to obtain
n n / 2 1 2 n / 2 1 + 1 2 n n / 2 1 1 2 n / 2 1
= 1 2 n n / 2 1 2 n / 2 1 + 1
Since there are 2 n + 1 Hamming neighbors x ˜ j , 0 x ˜ j , 1 and since the average sensitivity is n times the activity, we must multiply (18) by n 2 n + 1 , resulting in the statement of the theorem. □
The case of odd n is somewhat more involved because there is no “middle” layer E n , n / 2 . Instead, typical functions break up into two sets: M 0 , 1 1 n , r , z , v and M 0 , 2 1 n , r , z , v . In the first case, all minimal ones are on layers E n , n 3 / 2 , E n , n 1 / 2 , and E n , n + 1 / 2 while in the second case, all minimal ones are situated on E n , n 1 / 2 , E n , n + 1 / 2 , and E n , n + 3 / 2 . Under random sampling, these two cases will occur with equal probabilities. As we shall see, the results will be different for each case. Thus, the expected average sensitivity will be the average of the the expected average sensitivities corresponding to these two cases.
Theorem 5.
Let n be odd and let f M 0 1 n be a typical monotone Boolean function. Then, the expected average sensitivity s ^ f = E [ s f ] of f is
s ^ f 1 2 s ^ 1 f + s ^ 2 f ,
where s ^ 1 f and s ^ 2 f are given in (28) and (37), respectively.
Proof. 
Let us first address the set M 0 , 1 1 n , r , z , v . As in Theorem 4, we only need to concern ourselves with four cases. These are:
  • x ˜ E n , n 3 / 2 , f x ˜ = 1 , x j = 1
  • x ˜ E n , n 3 / 2 , f x ˜ = 0 , x j = 0
  • x ˜ E n , n + 1 / 2 , f x ˜ = 1 , x j = 1
  • x ˜ E n , n + 1 / 2 , f x ˜ = 0 , x j = 0
All other situations, such as x ˜ E n , n 3 / 2 , f x ˜ = 1 , x j = 0 , will result in the partial derivatives being zero due to monotonicity of f, hence will make no contribution to the activity of variable x j . There are
1 2 n n 3 / 2 2 n + 3 / 2
minimal ones on E n , n 3 / 2 such that x j = 1 . At the same time, there are
1 2 n n 3 / 2 1 2 n + 3 / 2
zeros on E n , n 3 / 2 such that x j = 0 . Unlike in the n-even case, where the middle layer contains an equal number of ones and zeros, the number of ones and zeros on E n , n 1 / 2 is not equal. The number of ones on E n , n 1 / 2 is given in (9). Thus, if x ˜ j , 0 is a zero on E n , n 3 / 2 , then the probability that f x ˜ j , 1 = 1 is
n n 1 / 2 1 1 2 n n 1 / 2 + n n 3 / 2 2 n + 3 / 2 n + 3 / 2                         n n + 1 / 2 2 n + 1 / 2 n + 1 / 2
where we have simply divided the number of ones on E n , n 1 / 2 by the total number of vectors on that layer. Thus, multiplying (21) by (22) and adding to (20) gives us the total contribution to the activity of variable x j from cases 1 and 2 above, which is equal to
1 2 n n 3 / 2 2 n + 3 / 2 + 1 2 n n 3 / 2 1 2 n + 3 / 2 ×      n n 1 / 2 1 1 2 n n 1 / 2 + n n 3 / 2 2 n + 3 / 2 n + 3 / 2                        n n + 1 / 2 2 n + 1 / 2 n + 1 / 2
Cases 3 and 4 are similar. The number of (maximal) zeros on E n , n + 1 / 2 for which x j = 0 is equal to
1 2 n n + 1 / 2 2 n + 1 / 2
and the number of ones on E n , n + 1 / 2 for which x j = 1 is equal to
1 2 n n + 1 / 2 1 2 n + 1 / 2 .
The proportion of zeros on E n , n 1 / 2 is
1 n n 1 / 2 1 1 2 n n 1 / 2 +        n n 3 / 2 2 n + 3 / 2 n + 3 / 2 n n + 1 / 2 2 n + 1 / 2 n + 1 / 2
Thus, as before, multiplying (25) by (26) and adding to (24), we get the total contribution to the activity from above and below the layer E n , n + 1 / 2 (cases 3 and 4), resulting in
1 2 n n + 1 / 2 2 n + 1 / 2 + 1 2 n n + 1 / 2 1 2 n + 1 / 2 ×              1 n n 1 / 2 1 1 2 n n 1 / 2 +         n n 3 / 2 2 n + 3 / 2 n + 3 / 2 n n + 1 / 2 2 n + 1 / 2 n + 1 / 2
Finally, adding (23) and (27) and then multiplying by n 2 n + 1 as in Theorem 4, we obtain the expected average sensitivity s ^ 1 f of a typical function from M 0 , 1 1 n , r , z , v :
s ^ 1 f n 2 n + 1 1 2 n n 3 / 2 2 n + 3 / 2 + 1 2 n n 3 / 2 1 2 n + 3 / 2 ×       n n 1 / 2 1 1 2 n n 1 / 2 + n n 3 / 2 2 n + 3 / 2 n + 3 / 2             n n + 1 / 2 2 n + 1 / 2 n + 1 / 2 +        1 2 n n + 1 / 2 2 n + 1 / 2 + 1 2 n n + 1 / 2 1 2 n + 1 / 2 ×            1 n n 1 / 2 1 1 2 n n 1 / 2 +      n n 3 / 2 2 n + 3 / 2 n + 3 / 2 n n + 1 / 2 2 n + 1 / 2 n + 1 / 2 .
We now proceed to derive the expected average sensitivity s ^ 2 f of a typical function from M 0 , 2 1 n , r , z , v , following the same steps.
Again, we only need to concern ourselves with four cases:
  • x ˜ E n , n 1 / 2 , f x ˜ = 1 , x j = 1
  • x ˜ E n , n 1 / 2 , f x ˜ = 0 , x j = 0
  • x ˜ E n , n + 3 / 2 , f x ˜ = 1 , x j = 1
  • x ˜ E n , n + 3 / 2 , f x ˜ = 0 , x j = 0
There are
1 2 n n 1 / 2 2 n + 1 / 2
minimal ones on E n , n 1 / 2 such that x j = 1 . At the same time, there are
1 2 n n 1 / 2 1 2 n + 1 / 2
zeros on E n , n 1 / 2 such that x j = 0 . As for M 0 , 1 1 n , r , z , v , the number of ones and zeros on E n , n + 1 / 2 is not equal. The number of ones on E n , n + 1 / 2 is given in (11). Thus, if x ˜ j , 0 is a zero on E n , n 1 / 2 , then the probability that f x ˜ j , 1 = 1 is
n n + 1 / 2 1 1 2 n n + 1 / 2 + n n 1 / 2 2 n + 1 / 2 n 1 / 2                        n n + 3 / 2 2 n + 3 / 2 n + 3 / 2
where we have simply divided the number of ones on E n , n + 1 / 2 by the total number of vectors on that layer. Thus, multiplying (30) by (31) and adding to (29) gives us the total contribution to the activity of variable x j from cases 1 and 2 above, which is equal to
1 2 n n 1 / 2 2 n + 1 / 2 + 1 2 n n 1 / 2 1 2 n + 1 / 2 ×       n n + 1 / 2 1 1 2 n n + 1 / 2 + n n 1 / 2 2 n + 1 / 2 n 1 / 2                        n n + 3 / 2 2 n + 3 / 2 n + 3 / 2
Cases 3 and 4 are similar. The number of (maximal) zeros on E n , n + 3 / 2 for which x j = 0 is equal to
1 2 n n + 3 / 2 2 n + 3 / 2
and the number of ones on E n , n + 3 / 2 for which x j = 1 is equal to
1 2 n n + 3 / 2 1 2 n + 3 / 2 .
The proportion of zeros on E n , n + 1 / 2 is
1 n n + 1 / 2 1 1 2 n n + 1 / 2 +        n n 1 / 2 2 n + 1 / 2 n 1 / 2 n n + 3 / 2 2 n + 3 / 2 n + 3 / 2
Thus, as before, multiplying (34) by (35) and adding to (33), we get the total contribution to the activity from above and below the layer E n , n + 3 / 2 (cases 3 and 4), resulting in
1 2 n n + 3 / 2 2 n + 3 / 2 + 1 2 n n + 3 / 2 1 2 n + 3 / 2 ×             1 n n + 1 / 2 1 1 2 n n + 1 / 2 +        n n 1 / 2 2 n + 1 / 2 n 1 / 2 n n + 3 / 2 2 n + 3 / 2 n + 3 / 2
Finally, adding (32) and (36) and then multiplying by n 2 n + 1 as in Theorem 4, we obtain the expected average sensitivity s ^ 2 f of a typical function from M 0 , 2 1 n , r , z , v :
s ^ 2 f n 2 n + 1 1 2 n n 1 / 2 2 n + 1 / 2 + 1 2 n n 1 / 2 1 2 n + 1 / 2 ×       n n + 1 / 2 1 1 2 n n + 1 / 2 + n n 1 / 2 2 n + 1 / 2 n 1 / 2             n n + 3 / 2 2 n + 3 / 2 n + 3 / 2 +         1 2 n n + 3 / 2 2 n + 3 / 2 + 1 2 n n + 3 / 2 1 2 n + 3 / 2 ×             1 n n + 1 / 2 1 1 2 n n + 1 / 2 +       n n 1 / 2 2 n + 1 / 2 n 1 / 2 n n + 3 / 2 2 n + 3 / 2 n + 3 / 2 .
Given that a function is equally likely to be picked from M 0 , 1 1 n , r , z , v as from M 0 , 2 1 n , r , z , v , the expected average sensitivity is the average of Equations (28) and (37). □

Funding

This research received no external funding.

Acknowledgments

I am grateful to the Institute for Systems Biology (ISB) for their support.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Kauffman, S.A. Metabolic stability and epigenesis in randomly constructed genetic nets. J. Theor. Biol. 1969, 22, 437–467. [Google Scholar] [CrossRef]
  2. Kauffman, S.A. The Origins of Order: Self-Organization and Selection in Evolution; Oxford University Press: Oxford, UK, 1993. [Google Scholar]
  3. Derrida, B.; Pomeau, Y. Random networks of automata: A simple annealed approximation. EPL (Europhys. Lett.) 1986, 1, 45. [Google Scholar] [CrossRef]
  4. Luque, B.; Solé, R.V. Lyapunov exponents in random Boolean networks. Physica A 2000, 284, 33–45. [Google Scholar] [CrossRef] [Green Version]
  5. Shmulevich, I.; Kauffman, S.A. Activities and sensitivities in Boolean network models. Phys. Rev. Lett. 2004, 93, 048701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Munoz, M.A. Colloquium: Criticality and dynamical scaling in living systems. Rev. Mod. Phys. 2018, 90, 031001. [Google Scholar] [CrossRef] [Green Version]
  7. Torres-Sosa, C.; Huang, S.; Aldana, M. Criticality Is an Emergent Property of Genetic Networks that Exhibit Evolvability. PLoS Comput. Boil. 2012, 8, e1002669. [Google Scholar] [CrossRef] [Green Version]
  8. Krawitz, P.M.; Shmulevich, I. Basin Entropy in Boolean Network Ensembles. Phys. Rev. Lett. 2007, 98. [Google Scholar] [CrossRef] [Green Version]
  9. Nykter, M.; Larjo, A.; Aho, T.; Kauffman, S.A.; Yli-Harja, O.; Shmulevich, I.; Price, N.D. Critical Networks Exhibit Maximal Information Diversity in Structure-Dynamics Relationships. Phys. Rev. Lett. 2008, 100. [Google Scholar] [CrossRef] [Green Version]
  10. Ribeiro, A.S.; Kauffman, S.A.; Lloyd-Price, J.; Samuelsson, B.; Socolar, J.E.S. Mutual information in random Boolean models of regulatory networks. Phys. Rev. E 2008, 77. [Google Scholar] [CrossRef] [Green Version]
  11. Galas, D.J.; Nykter, M.; Carter, G.W.; Price, N.D.; Shmulevich, I. Biological Information as Set-Based Complexity. IEEE Trans. Inf. Theory 2010, 56, 667–677. [Google Scholar] [CrossRef] [Green Version]
  12. Kang, C.; Aguilar, B.; Shmulevich, I. Emergence of diversity in homogeneous coupled Boolean networks. Phys. Rev. E 2018, 97, 052415. [Google Scholar] [CrossRef] [PubMed]
  13. Goudarzi, A.; Teuscher, C.; Gulbahce, N.; Rohlf, T. Emergent Criticality through Adaptive Information Processing in Boolean Networks. Phys. Rev. Lett. 2012, 108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Echlin, M.; Aguilar, B.; Notarangelo, M.; Gibbs, D.L.; Shmulevich, I. Flexibility of Boolean Network Reservoir Computers in Approximating Arbitrary Recursive and Non-Recursive Binary Filters. Entropy 2018, 20, 954. [Google Scholar] [CrossRef] [Green Version]
  15. Korshunov, A.D. Monotone boolean functions. Russ. Math. Surv. 2003, 58, 929–1001. [Google Scholar] [CrossRef]
  16. Aracena, J.; Gadouleau, M.; Richard, A.; Salinas, L. Fixing monotone Boolean networks asynchronously. Inf. Comput. 2020, 104540. [Google Scholar] [CrossRef] [Green Version]
  17. Aracena, J.; Richard, A.; Salinas, L. Number of fixed points and disjoint cycles in monotone Boolean networks. SIAM J. Discret. Math. 2017, 31, 1702–1725. [Google Scholar] [CrossRef] [Green Version]
  18. Aracena, J.; Demongeot, J.; Goles, E. On limit cycles of monotone functions with symmetric connection graph. Theor. Comput. Sci. 2004, 322, 237–244. [Google Scholar] [CrossRef] [Green Version]
  19. Melliti, T.; Regnault, D.; Richard, A.; Sené, S. Asynchronous simulation of Boolean networks. In International Conference on Cellular Automata; Springer: Cham, Switzerland, 2016; pp. 182–191. [Google Scholar]
  20. Mesot, B.; Teuscher, C. Critical values in asynchronous random boolean networks. In European Conference on Artificial Life; Springer: Berlin/Heidelberg, Germany, 2003; pp. 367–376. [Google Scholar]
  21. Shmulevich, I. Average sensitivity of typical monotone Boolean functions. arXiv 2005, arXiv:math/0507030. [Google Scholar]
  22. Korshunov, A. On the number of monotone Boolean functions. Probl. Kibern. 1981, 38, 5–108. [Google Scholar]
  23. Harris, S.E.; Sawhill, B.K.; Wuensche, A.; Kauffman, S. A model of transcriptional regulatory networks based on biases in the observed regulation rules. Complexity 2002, 7, 23–40. [Google Scholar] [CrossRef]
  24. Serra, R.; Villani, M.; Semeria, A. Genetic network models and statistical properties of gene expression data in knock-out experiments. J. Theor. Biol. 2004, 227, 149–157. [Google Scholar] [CrossRef] [PubMed]
  25. Li, Y.; Adeyeye, J.O.; Murrugarra, D.; Aguilar, B.; Laubenbacher, R. Boolean nested canalizing functions: A comprehensive analysis. Theor. Comput. Sci. 2013, 481, 24–36. [Google Scholar] [CrossRef]
  26. Shmulevich, I.; Lähdesmäki, H.; Dougherty, E.R.; Astola, J.; Zhang, W. The role of certain Post classes in Boolean network models of genetic networks. Proc. Natl. Acad. Sci. USA 2003, 100, 10734–10739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Kahn, J.; Kalai, G.; Linial, N. The Influence of Variables on Boolean Functions. In Proceedings of the 29th Annual Symposium on Foundations of Computer Science, SFCS ’88, White Plains, NY, USA, 24–26 October 1988; pp. 68–80. [Google Scholar] [CrossRef]
  28. Shmulevich, I.; Dougherty, E.R.; Kim, S.; Zhang, W. Probabilistic Boolean networks: A rule-based uncertainty model for gene regulatory networks. Bioinformatics 2002, 18, 261–274. [Google Scholar] [CrossRef] [PubMed]
  29. Friedgut, E.; Kalai, G. Every monotone graph property has a sharp threshold. Proc. Am. Math. Soc. 1996, 124, 2993–3002. [Google Scholar] [CrossRef]
  30. Bernasconi, A.; Damm, C.; Shparlinski, I. The average sensitivity of square-freeness. Comput. Complex. 2000, 9, 39–51. [Google Scholar] [CrossRef]
  31. Bernasconi, A. Sensitivity vs. block sensitivity (an average-case study). Inf. Process. Lett. 1996, 59, 151–157. [Google Scholar] [CrossRef]
  32. Boppana, R.B. The average sensitivity of bounded-depth circuits. Inf. Process. Lett. 1997, 63, 257–261. [Google Scholar] [CrossRef]
  33. Friedgut, E. Boolean functions with low average sensitivity depend on few coordinates. Combinatorica 1998, 18, 27–35. [Google Scholar] [CrossRef]
  34. Shi, Y. Lower bounds of quantum black-box complexity and degree of approximating polynomials by influence of Boolean variables. Inf. Process. Lett. 2000, 75, 79–83. [Google Scholar] [CrossRef] [Green Version]
  35. Mossel, E.; O’Donnell, R. On the noise sensitivity of monotone functions. Random Struct. Algorithms 2003, 23, 333–350. [Google Scholar] [CrossRef] [Green Version]
  36. Zhang, S. Note on the average sensitivity of monotone Boolean functions. Available online: http://www.cse.cuhk.edu.hk/~syzhang/papers/Sensitivity.pdf (accessed on 1 June 2020).
  37. Bshouty, N.H.; Tamon, C. On the Fourier spectrum of monotone functions. J. ACM (JACM) 1996, 43, 747–770. [Google Scholar] [CrossRef] [Green Version]
  38. Korshunov, A.D.; Shmulevich, I. On the distribution of the number of monotone Boolean functions relative to the number of lower units. Discret. Math. 2002, 257, 463–479. [Google Scholar] [CrossRef] [Green Version]
  39. Shmulevich, I.; Sellke, T.; Gabbouj, M.; Coyle, E. Stack filters and free distributive lattices. In Proceedings of the 1995 IEEE Workshop on Nonlinear Signal and Image Processing, Neos Marmaras, Greece, 20–22 June 1995; pp. 829–930. [Google Scholar]
  40. Korshunov, A.D.; Shmulevich, I. The number of special monotone Boolean functions and statistical properties of stack filters. Diskretn. Anal. Issled. Oper. 2000, 7, 17–44. [Google Scholar]
  41. Shmulevich, I.; Yli-Harja, O.; Astola, J.; Korshunov, A. On the robustness of the class of stack filters. IEEE Trans. Signal Process. 2002, 50, 1640–1649. [Google Scholar] [CrossRef]
Table 1. Table of notation.
Table 1. Table of notation.
fBoolean function
nnumber of variables of f
f ( x ˜ ) / x j partial derivative of f with respect to x j
x ˜ = x 1 , , x n vector of n binary values
α j f activity of variable x j
s f average sensitivity of a Boolean function f
ρ ( α ˜ , β ˜ ) Hamming distance between α ˜ and β ˜
M n set of monotone Boolean functions of n variables
M 0 n set of special monotone Boolean functions of n variables
E n Boolean n-cube or 0 , 1 n with vertices x ˜
E n , k { x ˜ E n | ρ ( x ˜ , 0 ˜ ) = k } or the kth layer of E n

Share and Cite

MDPI and ACS Style

Shmulevich, I. On the Lyapunov Exponent of Monotone Boolean Networks . Mathematics 2020, 8, 1035. https://doi.org/10.3390/math8061035

AMA Style

Shmulevich I. On the Lyapunov Exponent of Monotone Boolean Networks . Mathematics. 2020; 8(6):1035. https://doi.org/10.3390/math8061035

Chicago/Turabian Style

Shmulevich, Ilya. 2020. "On the Lyapunov Exponent of Monotone Boolean Networks " Mathematics 8, no. 6: 1035. https://doi.org/10.3390/math8061035

APA Style

Shmulevich, I. (2020). On the Lyapunov Exponent of Monotone Boolean Networks . Mathematics, 8(6), 1035. https://doi.org/10.3390/math8061035

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop