Controllability of time-delayed Boolean multiplex control networks under asynchronous stochastic update.

In this article, the controllability of asynchronous Boolean multiplex control networks (ABMCNs) with time delay is studied. Firstly, dynamical model of Boolean multiplex control networks is constructed, which is assumed to be under Harvey' asynchronous update and time delay is introduced both in states and controls. By using of semi-tensor product (STP) approach, the logical dynamics is converted into an equivalent algebraic form by obtaining the control-depending network transition matrices of delayed system. Secondly, a necessary and sufficient condition is proved that only control-depending fixed points of the studied dynamics can be controlled with probability one. Thirdly, respectively for two types of controls, the controllability of dynamical control system is investigated. When initial states and time delay are given, formulae are obtained to show a) the reachable set at time s under specified controls; b) the reachable set at time s under arbitrary controls; c) the reachable probabilities to different destination states. Furthermore, an approach is discussed to find a precise control sequence which can steer dynamical system into a specified target with the maximum reachable probability. Examples are shown to illustrate the feasibility of the proposed scheme.

Boolean networks (BNs), as a class of simplified discrete models, are widely applied to reveal the generic properties of biological systems in an integrative and holistic manner. In 1969, random Boolean networks (RBNs), also known as N-K models, were originally proposed by Stuart Kauffman1. This classic model consists of N nodes representing genes, each of which receives inputs from K randomly selected neighbors. The nodes on networks are characterized by two qualitative values, usually referred to logical 0 and 1, to present the active and repressing states of genes, respectively. Boolean rules assigned to nodes are employed to indicate the mutual regulations among genes. Based on synchronous update, at each time t, the state of each node on network is determined by the Boolean rule and K inputs at the previous time t-1. During the past few decades, Boolean networks have been used to unveil characteristics of complex systems and abundant of results have been achieved, such as dynamical behaviors of BNs23, efficient attractor-seeking algorithms456, biological control78 and some applications in biological research91011.

Biological information operates on multiple hierarchical levels of living organization12. From the viewpoint of systems biology, system-level analysis of biological regulation requires the interactions of genes on a holistic level, rather than the characteristics of isolated parts of an organism13. As mentioned in Ref. 14, “the same gene or biochemical species can be involved in a regulatory interaction, in a metabolic reaction, or in another signaling pathway”. Therefore, to understand the intricate variability of biological systems, where many hierarchical levels and interactions coexist, a new level of description is required. Meanwhile, multiplex networks as an extension of complex networks were firstly proposed by Mucha in 201015, which is composed of several layered networks interrelated with each other shown in Fig. 1. Each layer in multiplex networks could have particular features and dynamical processes. Interconnections between layers are represented by some special nodes on behave of different roles participating in multiple layers of interactions. Different from the traditional sense of coupling, the final states of those common nodes at each time step are determined by all of involved layers. During the past four years, a variety of studies based on multiplex networks have been achieved, including network topology and dynamic properties16, diffusion dynamics1718 and game theory1920, etc. It's noteworthy that multiplex networks provide a novel way to construct the multilevel models of biochemical systems and be better depict a richer structure of interactions.

Figure 1

Illustration of multiplex networks with two layers.

Nodes a, b, c, d are identical in both layer 1 and layer 2.

To incorporate mutual regulations of genes observed in the real biological system, Boolean networks as abstract models are employed to investigate dynamical properties of systems. For a certain degree of simplification, synchronous update scheme is adopted in the previous studies of RBNs, which is based on an assumption that update scheme isn't an essential factor on the consideration of dynamical behaviors. However, regulated entities can't implement the interactions and renew their states simultaneously at each time step by following a synchronized clock. As the studies in Ref. 21, “factors such as mRNA and protein synthesis, degradation and transport times mean that the system is replete with delays of varying amounts, and genes are activated or inhibited in a fundamentally asynchronous manner”, that means there are multiple timescales should be considered in the biological systems. Asynchronous Boolean networks (ABNs) were firstly proposed by Harvey et al.22, followed by a series of related studies. In23, Greil et al. illustrated that the growth of the mean number and size of attractors in asynchronous critical BNs are in strong contrast to the synchronous version. Furthermore, dynamics of critical BNs under deterministic asynchronous update was also studied24; in25, Saadatpour et al. carried out a comparative study on the attractors of a signal transduction network modelled by BNs under synchronous and asynchronous updating schemes; in26, Tournier and Chaves investigated dynamics of the interconnection of two ABNs by directly analyzing the properties of two individual modules, that can be applied to analyze the multicellular modeling and high dimensional model; in27, asynchronous stochastic Boolean networks were proposed to investigate dynamical behaviors of a T-helper network; in28, Jack et al. simulated quantitative cellular responses of signal transduction in a single cell by means of asynchronous threshold BNs. The previous results were presented to verify the asynchronous update is more plausible for many cases of biological systems. Therefore, studies of BNs under asynchronous stochastic update are meaningful and applicable.

From system-level understanding of biological systems, to find a mechanisms that systematically control the states of regulatory entities can be implemented to minimize malfunctions and provide potential therapeutic targets for treatment of disease2930. Boolean control networks (BCNs) as dynamical control systems provide an efficient approach to carry out theoretical and numerical analysis31. Controllability, roughly speaking, which is to steer a control system from an arbitrary initial state to an arbitrary final state by using the set of admissible controls, is one of the fundamental concepts at the onset of control theory infiltrating into the research of gene regulatory networks (GRNs)3233. Over the past few years, controllability of BCNs have been receiving considerable attention, such as controllability and observability of BCNs34; controllability of BCNs with time-invariant delay35 and time-variant delay36 as well as time delay involved both in states and controls37; controllability of μ-th order BCNs38 and the approach to transform μ-th order BCNs to equivalent time-variant BCNs39; studies on controllability of BCNs via the Perron-Frobenius theory40, etc. The previous works were based on an independent network, which represents the assembly of genes or other entities to fulfil a specified function. With the in-depth of research, it is necessary to further study how the interplay among multiple interdependent networks affects dynamical behaviors of system. Compared with the traditional models, Boolean multiplex networks have more complex topology structure and higher holistic level which can provide a more generalized model to be better in conformity with the development of biology. Moreover, most of the previous studies on controllability of BCNs were assumed to be updated under synchronous scheme, i.e. the studied models are deterministic systems. As the above discussion, asynchronous update is closer to the real situation, based on which studies on BNs can be more likely to obtain the essential properties of biological systems. In41, the reachable sets of Boolean multiplex networks under asynchronous update scheme at time s were revealed, where the asynchronous scheme was based on randomly chosen update nodes at each time step. However, Due to signal propagation delays in the environment, a propagation delay τ can be seen as a particular form of asynchronous phenomenon existing in the processes of transcription and translation in biological systems42 or information propagation in society systems43. Hence, we think it's valuable to extend the related research into the field of asynchronous Boolean multiplex networks with time delay.

In this article, controllability of ABMCNs with time delay is discussed. The dynamical model of Boolean multiplex control networks is constructed by introducing inputs as controls into the model proposed by Cozzo et al.14. For obtaining the more general results, time delay is involved both in states and inputs44. Harvey's update scheme, i.e. only one node could be randomly chosen to renew its state at each time step, is implemented. In45, as a kind of non-deterministic system, the controllability of probabilistic Boolean networks was discussed, in which the concept of controllable probability was firstly proposed. But, authors just showed the sufficiency of the controllability with probability but not verify the necessary. In our work, a necessary and sufficient condition is proved that only control-depending fixed points of asynchronous delayed system can be controlled with probability one, which provide the theoretical basis to discuss the controllability of non-deterministic system from the perspective of probability. Based on the algebraic representation of the studied model, controllability of delayed system is to be analytically discussed respectively for two types of controls, i.e. free Boolean control sequences and the controls satisfying certain logical rule. When initial state sequence and time delay are given, we discuss the formulae to calculate reachable sets at time s under specified or free controls, as well as the reachable probabilities to different destination states. Furthermore, we are to illustrate the method to determine specific controls which can drive dynamical system to a given target with the maximum reachable probability.

This article is organized as follows. In Preliminaries, semi-tensor product as mathematic tools applied in this article is briefly introduced. In Main Results, the studied model of ABMCNs with time delay is firstly proposed and converted into linear form. Based on two types of controls, the controllability of dynamical control system is discussed. Some examples are shown to illustrate the main results. Finally, a concluding remark is given.

Preliminaries

In this section, STP of matrix is briefly introduced, by means of which logical dynamics can be converted into an equivalent algebraic form.

Definition 1 (31):

Let be a row vector of dimension np, and be a column vector of dimension p. Then we split into p equal-size blocks as 1, 2, ..., , which are 1 × n rows. Define the STP, denoted by , as

Let ∈ M and ∈ M. If either n is a factor of p, say nt = p and denote it as , or p is a factor of n, say n = pt and denote it as , then we define the STP of and , denoted by , as the following: consists of m × q blocks as = () and each block is where is the i-th row of and is the j-th column of .

Example 1: Let and. Then, one can obtain

Remark 1: It is noted that when n = p, STP of and turns into the conventional matrix produce. So, STP can be seen as a generalization of the conventional matrix product and all the fundamental properties of matrix product, such as distributive rule, associative rule, etc, still hold.

And, it can be verified that for two column vectors and , .

Some related properties of STP are collected as follows:

Proposition 1: Assume , then (where refers to the Kronecker product, is the identity matrix)

Assume , then

Proposition 2: Assume ∈ M is given,

Let be a row vector. Then,

Let be a column vector. Then,

For statement ease, some notations used in this article are defined as follows.

denotes the r-th column of the n×n identity matrix and , which is the set of all columns of .

A matrix ∈ M can be called a logical matrix if , which is briefly denoted by . And the set of n×m logical matrices is denoted by .

Next, we define the swap matrix , let and be two column vectors where [ is a mn × mn matrix labeled columns by and rows by , the elements in position ((I, J), (i, j)) is [ is briefly denoted by [.

Assume and (t) ∈ Δ2, we can get , where . Here, M = δ4[1,4], which is power-reducing matrix and it can be verified that 2 = M, .

In order to get the matrix expression of logical dynamics, the Boolean values should be denoted as vectors and . And the following lemma is fundamental for the matrix expression of logical functions.

Lemma 131: Any logical function with logical arguments , can be expressed in a multi-linear form as where M ∈ 2×2 is unique, which is called the structure matrix of logical function f.

More details on STP can be found in Ref. 31. In the following, the matrix products are assumed to be STP and the symbol is omitted if no confusion arises.

Main Results

Algebraic expression of asynchronous Boolean multiplex control networks with time delay

Regulatory entities in multiplex take part in several layers of networks, the states of which on different layers evolve independently. However, a final deterministic state of each entity should be obtained at the end of each time step determined by all of values on involved layers.

For a Boolean multiplex network with nodes and layers, assume accounts for the state of node i on layer l at time t. When time delay τ is considered in states, one can obtain where is the update function of node i on layer l. Furthermore, assume represents the overall state of node i at time t. Refer to14, we can get where is the canalizing function. Boolean functions are canalizing if whenever the canalizing variable takes a given value, the function always yields the same output, irrespective of the values of other variables14. Note that, strictly speaking, there exists an interval between the renewal of the value of node i on layer l, say , and the overall state in the whole multiplex. In the following discussion, based on an assumption that the interval between the above two states is instantaneous, the same time step t is used for both of them.

Next, we introduce controls with time delay τ into system (1), the corresponding dynamical control system can be described as where u(t), are controls and is the update rule of node i on layer l with controls.

By means of Lemma 1, a structure matrix can be calculated for each logical rule , based on which one can obtain the algebraic form of Eq.(3) as follows. where , . Subsequently, the algebraic representation of Eq.(2) can be obtained as where is the structure matrix of logical function and .

Under Harvey's asynchronous update, at each time step t, only one node is at random chosen for update. Hence, one can obtain Multiplying all the equations of system (6), one can get wher is called as the control-depending network transition matrix, which involves all of the state transfer information of a dynamical control system.

In the following, respectively for two kinds of controls, the controllability of ABMCNs with time delay is to be discussed:

Controls come from a free Boolean sequence. Precisely, at time t, controls are freely designed and described as .

The controls are determined by certain logical rules, which can be called input control networks: where are logical rules.

Deterministic controllability of asynchronous Boolean multiplex control networks with time delay

Synchronous BNs are deterministic dynamical systems, however, under Harvey's asynchronous update scheme, different update choices can be randomly chosen with the same probability at each time step. Correspondingly, when logical system is converted into linear form, there are different control-depending network transition matrices , . Say, the average probability for each transition matrix is . Then, one can obtain

Definition 2: Consider system (6) with time delay τ, given an initial state sequence , , the destination state is said to be controllable with probability one at time s > 0, if a group of controls (t), can be found such that . Noted that is the smallest integer larger than or equal to a, for instance, .

Remark 2: When time delay τ and time s are given, according to the discussed model, i.e. Eq. (1), the previous location should be s−(1+τ), continue to induce, one can obtain . Since an initial state sequence , is given, one can verify that location should be in the scope of , i.e. the initial state would be , .

Definition 3: As to system (6) with time delay τ, when a control exists such that state holds , is said to be a control-depending fixed point.

Theorem 1: Consider system (6) with time delay τ, when an initial state sequence () is given, the destination state is said to be controllable at time s> 0 with probability one, only and if only state is a control-depending fixed point.

Proof:

(Sufficiency)Assume the destination state is a control-depending fixed point of system (6) with time delay τ. According to Definition 3, a control can be found that . Consequently, we can find a group of controls (t) = , . Then, one can obtain from the initial state . So, can be controllable from itself with probability one at time s.

(Necessity) When the destination state is said to be controllable with probability one at time s>0 from initial state , should be proved to be a control-depending fixed point. Firstly, we assume . According to Definition 2, a control (s−τ−1) can be found that , which means , where and . When , considering the rule of Harvey's update scheme, there should be only one node and the rest elements . And when , there should be only one node and the rest elements . The two results are contradictory. So the above assumption can't be held. Say, should be equal to . Deduce the rest from this, one can obtain , . According to Definition 3, it can be proved that should be a control-depending fixed point of system (6) with time delay τ.

This completes the proof.

From the above results, we can conclude that, as to system (6) with time delay τ, when two states are given and , can't be controllable from with probability one. Hence, it is necessary and reasonable to discuss the controllability of the studied model from the perspective of probability.

Controllability of asynchronous Boolean multiplex control networks with time delay via free Boolean sequence

In this section, controls are assumed to be free Boolean sequences, based on which the controllability of studied model is discussed.

Definition 4: Given the initial state sequence , and the destination state , system (6) with time delay τ is said to be controllable to with probability at time s > 0, if a group of controls (t), can be found such that .

When the initial states and a control sequence are specified, the following approach can be used to calculate the reachable set with probability at time s.

Before the next discussions, we define two operations:

, furthermore, when , . Correspondingly, and . For instance, and .

Let column vector ∈ R, all of row indies of in which row elements aren't equal to zero compose a set denoted by Ω (). For example, = [1,0,2,1]T and Ω() = {1,3,4}.

Theorem 2: For system (6) with time delay τ, given the initial state sequence , and controls (t), , the destination state is reachable with probability at time s, iff where , , , and represent the φ-th column of matrix.

Proof: By means of STP, system (6) with time delay τ can be rewritten as At time t, since each node in multiplex has the same probability to be chosen for update, one can obtain the overall expected value of as Since time delay τ is involved both in states and controls, the location of the initial state which evolves into the destination state is . To expand the above formula and yields

where .

This completes the proof.

When the initial states are given and controls are freely chosen, we provide the following approach to calculate the reachable set with probability at time s. Assume is the set of initial states, we denote by R(X0) the reachable set from set X0 with time delay τ at time s under arbitrary controls.

Lemma 2: For system (6) with time delay τ, is the set of initial states. Controls (t), can be freely chosen, one can obtain where , , and Row(·) represents the i-th row of matrix.

Proof:

1) Assume state , should be proofed.

Since the destination state is reachable with probability at time s, one can find a sequence of controls to steer the system from initial state to the destination states . Correspondingly, based on Theorem 2, it's easy to get the element in the position (r, φ) of matrix should be non-zero, which means the φ-th element of row vector is non-zero.

2) When , we can assume the φ-th element of row vector is non-zero. According to Theorem 2, by means of , the destination states is reachable with probability at time s. Furthermore, can be decomposed into a sequence of controls as .

This completes the proof.

And, when a control sequence is given, we also can obtain the specific reachable probability from certain initial states to a given destination state at time s.

Lemma 3: For system (6) with time delay τ, assume the initial state sequence as , and controls as . The reachable probability from the initial states to the destination state at time s is where , , and (·) is the element at position (i, j) of matrix.

Remark 3: Entry (β, φ) of matrix indicate the state transfer information of dynamics from initial state under control sequence after time steps to destination state .

Controllability of asynchronous Boolean multiplex control networks with time delay via input control networks

Based on STP of matrix, the linear representation of system (7) can be obtained as where is the network transient matrix of input control network.

Definition 5: Consider system (6) with input control network (7) and time delay τ, when initial state and destination state are given, is said to be controllable with probability from at time s, if an initial control can be found such that where .

Theorem 3: For system (6) with input control network (7) and time delay τ, the destination state is controllable with probability from initial state under initial control 0 at time s iff where , , .

Proof:

One can obtain

This completes the proof.

Lemma 4: For system (6) with input control network (7) and time delay τ, when the initial control 0 can be freely chosen, the set of states which are reachable with probability from initial states at time s is where .

Lemma 5: For system (6) with input control network (7) and time delay τ, the probability from the initial states to the destination state under initial control 0 at time s is where .

Examples

Example 1 Consider Boolean multiplex control network (16) with layers, nodes and control shown in Fig 2. Assume system (16) is under Harvey's asynchronous update and time delay τ both in states and controls. where , , , and represent the logical functions of negation, disjunction, conjunction, implication and equivalence, respectively. Correspondingly, one can obtain the algebraic representation of logical functions as , , , and .

Figure 2

An asynchronous Boolean multiplex control network with time delay (16).

Based on the above discussion, we define and . And, as to the canalizing function , without loss of the generality, we choose disjunction function, i.e. . The control 1(t) in system (16) is free Boolean variable. In the following, the controllability of ABMCNs (16) with time delay τ is to be discussed. Firstly, we calculate the control-depending network transition matrix of system. Note that, at time t, .

Case 1: at time t, when node 1 is selected for update,

Case 2: at time t, when node 2 is selected for update,

Case 3: at time t, when node 3 is selected for update,

Case 4: at time t, when node 4 is selected for update,

Therefore, all of the control-depending network transition matrices can be calculated as follows.

Assume time step s = 7 and time delay τ = 2, randomly choose the initial states , and . One can obtain and the initial state . Respectively for the free control sequence is given or arbitrary, the reachable set of system (16) is discussed as follows.

One can obtain Firstly, we assume the control sequence is specified. According to Theorem 2, controls . Hence, when a control sequence is given as , , , one can obtain . By means of Theorem 2, the reachable set from initial state under the given controls at time 7 can be calculated as i.e., as to the nondeterministic system (16), there are totally 8 states which have the possibility to be reached under the specified controls from initial state . Moreover, one can obtain the destination state has the maximum probability 21/64.

When the control sequence is arbitrary, based on Lemma 2, we can calculate the corresponding reachable set R(X0). For matrix (17), the row vectors in 3th, 4th, 7th, 8th and 12th rows are zero, that mean states are unreachable from the initial state with time step s = 7 and time delay τ = 2. Hence, we can obtain the reachable set In Fig 3, the reachable states of system (16) from under free controls with time delay τ = 2 in 3 steps are depicted.

Figure 3

The state transfer graph of system (16) from initial state (1011) in 3 steps.

When a destination state is given, different controls can steer system from the initial states into the target with different probabilities. Since a control sequence can't be found to make system turn into the target with possibility one, hence, controls which can get the maximum probabilities become the focus of attention. In virtue of the matrix (17), we can conveniently obtain the expected control sequence. Assume , the maximum probability 19/64 at (6, 1) of matrix (17), that means control can steer the initial state to destination state (1, 0, 1, 0) at time step 7 with probability 19/64. Subsequently, we can calculate , i.e. 1(−2) = 1, 1(1) = 1, 1(4) = 1.

Example 2 In the process of the cell cycle, the onset of M (mitosis) and S (DNA replication) phases are directed by the periodic activation of cyclin-dependent kinases (cdk's). Romond et al.46 constructed the differential equations model to reflect the above dynamics and Heidel et al.4 proposed the corresponding Boolean model. Based on the previous studies, considering time delay in the process, Boolean control model was further extended into multiplex architecture as follows. where +~M = δ2[2,1,1,2] and ·~M = δ2[1,2,2,2].

The controls in system (18) are produced by input control network as follows By means of STP, one can obtain linear representation of system (19) as (t + 1) = 1(t − τ)2(t − τ) = M2(t−τ)1(t − τ) = M[2](t − τ) and = M[2] = δ4[3,1,4,2].

Case 1: at time t, when node is selected for update,

Case 2: at time t, when node is selected for update, Case 3: at time t, when node is selected for update,

Case 4: at time t, when node is selected for update,

Assume time step s = 9 and time delay τ = 3, randomly choose the initial states as , , and . One can obtain and the initial state . According to Theorem 3, we can calculate that

Using Lemma 4, when the initial control 0 is free, except three unreachable states , and , all of the rest states can be reachable from initial state at time s = 9 with time delay τ = 3. Furthermore, when initial control is assumed to be specified, for instance, , one can obtain the corresponding reachable set as follows.

Note that states , and can be reached with the same probability 0.1406 from initial state under initial control = (1,0) at time s = 9 with time delay τ = 3. Similarly, destination states and can be reached with the same probability 0.1094, which is also the minimum reachable probability compared with the rest reachable states. Correspondingly, we can obtain the maximum reachable probability belonging to state is 0.3594.

In some applications, such as the therapeutic intervention, normally a final target is clear, i.e. an expected state for biological system is given. Hence, we should find a specific control sequence to steer system from the initial state to target with the maximum probability. Based on the above discussion, an approach can be obtained. Assume a required target at time s = 9 with time delay τ = 3 is . By using of matrix (20), we can get the maximum reachable probability is 0.2656 at the Row 2 and Column 3. According to Theorem 3 and Lemma 5, we can calculate the initial control , i.e. u1(0) = 0 and u2(0) = 1.

Conclusions

In this article, inputs as controls are introduced into Boolean multiplex networks under asynchronous stochastic update, meanwhile, time delay as additional factor is considered both in states and inputs of system. By means of STP approach, the above logical dynamics is converted into algebraic form and the controllability of dynamics is discussed. Firstly, it is proved that only control-depending fixed points can be controlled with probability one, which means the discussion of controllability of asynchronous Boolean control networks should be in terms of probabilities. Subsequently, respectively for two kinds of controls, formulae to calculate the reachable set from an initial state to a destination state under specified controls or arbitrary controls are provided, as well as the approach to obtain the specific reachable probabilities from an initial state to different destination states. Moreover, we also present how to find a precise control sequence which can steer dynamics into a given target with the maximum reachable probability.

Author Contributions

L.C., W.X.Y. devised the model and carried out theoretical analysis. L.C. and L.H. implemented numerical simulations. L.C., W.X.Y. and L.H. wrote the main text of the manuscript.