Scientific Reports. Dec/31/2016; 7

Published online Apr/27/2017

PMID: 28452370

PMC: 5408230

Science Direct: srep46495

doi: 10.1038/srep46495

Observability of Boolean multiplex control networks

Abstract

Human Genome Project, which is an international scientific research project with the goal of determining the sequence of nucleotide base pairs1, inspired a new view of biology called the systems biology. Instead of investigating individual genes, proteins or cells, systems biology studies the behavior and relationships of all cells, proteins, DNAs and RNAs in the same biological system, called a cellular network2. The Boolean Networks (BNs) as a powerful tool in describing, analyzing, and simulating the cellular networks, has been most widely used3 4 5 6 7 8 9 10 11 12 13 14 15 16 17.

From decades ago, Kauffman put forward the theory which can describe the net of cell and gene using BNs4. And he made expatiation about the relationship between BNs and gene as well as life5 6. Because the construction and evolutionary process of cell and gene can be revealed very well by BNs, BNs turn into a hot topic concerned by biologists, physicists and neuroscientists. Huang, S. et al. talked about the Boolean modeling and analyzing of biological system10 11. Aldana, M. et al. studied the topological structure of BNs7. Akutsu, T. et al. and Albert, R. et al. considered the dynamic features of BNs12 13. More detailedly and recently, Lu, J. et al. analyzed the synchronization problem of master-slave probabilistic BNs18.

In recent years accompany with the development of biology, control of biological system becomes into a hot topic19 20 21 22 23 24 25 26 27 28 29. As to the research of genetic regulatory networks (GRNs), one of the major goals is to carry out the therapeutic intervention strategies for diseased targets30 31. Correspondingly, Boolean control networks (BCNs) as a theoretical branch of the above studies provide an efficient way to investigate the control of GRNs based on abstract models. So the interests to the BCNs are increasingly going up. The application of BCNs includes not only GRNs32, but also other various fields, such as man-machine dynamic game33 and internal combustion engines34 35. Recently, based on semi-tensor product (STP) proposed by Cheng, D. et al.36, many basic problems for BCNs have been investigated, for example, realization23, controllability24 26, optimal control15 33, etc.

Observability of a system is a structural property. It is also a fundamental concept in control theory and systematic science and, not surprisingly, it has found many applications in systems biology. As early as 1976, Cobelli et al. studied controllability, observability and structural identifiability of biological compartmental systems of any structure37. In evolutionary dynamics, observability is the key to study whether the genetic process itself can be recovered from measurements of phenotypic characteristics38. Observability analysis is a necessary preliminary step to the design of observers, that is, systems that provide an estimate of the complete internal state based on measurements of the inputs and outputs39. There have been many studies on the observability of BCNs in recent years. Cheng, D. et al. have investigated the controllability and observability of BCNs24. Li, F. et al. have studied the observability of time-delayed BCNs25. Laschov, D. et al. have considered the observability of BNs through a graph-theoretic approach39. Zhang, K. et al. have proposed a unified approach to determine all the four types of observability of BCNs in the literature40.

From the view of systems biology, the analysis in system-level of biological regulation needs to consider the interactions of genes on a holistic level, rather than the independent characteristics of isolated parts of an organism41. To understand the intricate variability of biological systems, where many hierarchical levels and interactions coexist, a new level of description is required. Thereupon, multiplex networks as an extension of complex networks were firstly proposed by Mucha in 201042, which is composed of several layered networks interrelated with each other shown in Fig. 1. The previous description implicitly assumes that all biochemical signals are equivalent and then collapses information from different pathways. Actually, in cellular biochemical networks, many different signaling channels do work in parallel43. Not only in cellular biochemical, multiplex networks have been applied to the natural, social, and information sciences42. As an old concept, multi-layer social networks have been studied from 197544. Transportation systems are natural candidates for a multi-layer network representation. In a recent paper, a two-layer structure has been created by merging the world wide port and airport networks45. In multiplex networks, each layer could have particular features and dynamical processes. Between layers, interconnections are represented by some special nodes on behalf of different roles participating in multiple layers of interactions. The final states of those common nodes at each time are determined by all involved layers, which is different from the traditional sense of coupling.

Recent years more and more researchers studied the BMNs. For example, Xu, M. et al. investigated the synchronizability of two-layer networks46. Cozzo talked about the stability of BMNs43. Luo et al. studied the controllability of BMCNs27. Zhong, J. et al. studied controllability problem for multi-level Boolean control networks47. But when it comes to the observability problem of BMCNs, to our best knowledge, there have been even no results, because there are many differences between BMNs and single-layer BNs. Even for the degenerated BMNs, their observability are different from the single-layer BNs’, for example the BCNs studied by Cheng, D. et al.24 and Li, F. et al.25. Because even when the number of layer is one, our system still has holistic states, which have logical relationship with the states in basic layers. From above discussions, we can know that a study of the observability of the BMNs is meaningful and challenging.

In this article, by following this stream of research, we first address and characterize observability of BCNs defined on multiple topological layers. Based on the model of Boolean multiplex networks presented by Coozo et al.43, we introduce the input controls and the outputs. The model of BMCNs are changed into algebraic representation using STP tools. We consider the observability of BMCNs, following the standard formulation of the observability problem in systems theory, namely, we assume that the BN structure is known and that the goal is to infer the initial condition based on an output sequence. To clearly show the results of our research, we gave observable and unobservable examples in the final part of our essay.

The rest of this article is organized as follows. In Section II, we introduce the dynamic structure of BMCNs. In Section III, some concepts and properties of the STP are introduced, and we change our network into discrete form using STP tool. In Section IV, necessary and sufficient conditions for the observability of the BMCNs are obtained. Examples are given to show the effectiveness of the obtained results in Section V. Finally, a brief summary is given in Section VI.

Model of BMCNs

In this section, we introduce the model of BMCNs. For multiplex networks, different from the single-layer model, some nodes exist in multiple layers, the states which on different layers evolve independently of each other. The multiplex network we defined has N nodes per layer and K layers, and the number of total different nodes is n (N ≤ n ≤ NK). For example in Fig. 1, we have that N = 4, K = 2, n = 5. For statement ease, we define some related notions.

is the set of {0, 1}.
, and a_i,l = 1 if node i in the l layer and 0 otherwise. The layers set of node i is which refers the set of l which has a_i,l = 1.
, and γ_i,j,l = 1 if node j is the incoming neighbors of node i in the l layer. The incoming neighbors set of node i in the l layer is which refers the set of j which has γ_i,j,l = 1. And we set .

In Fig. 1, we have that the layers set of node 1 is , and a_1,1 = 1 and a_1,2 = 1, the layers set of node 2 is , and a_2,1 = 1 and a_2,2 = 0. The incoming neighbors set of node 1 in layer 1 is Γ_in(1)(1) = {1,4}, and .

In each layer, for the specific , if a_i,l = 1, we assume represents the state of node i on layer l at time t, then the update dynamics of state can be described as

where , is the update rule of node i on layer l.

Furthermore, assume represents the holistic states of node i at time t, which means the global states of , see Fig. 2. is influenced by , so we describe it as

where.. is the canalizing function.

When considering control-related problems for BMNs, based on above system structure, we introduce the m-dimensional control as the inputs of the system, correspondingly, then we have the outputs , then the BMCNs can be described as

and

where is the canalizing function of node i with the controls , see Fig. 3.

In finally, the output dynamics of the BMCNs are given by the following equation

where h_j is the output function.

Remark 1. The BMCNs are not the simple superposition of K single-layer BCNs. Because we have the holistic states which are affected by corresponding states in each layers. Between the holistic states with the states in basic layers, we have the canalizing functions, which determine the value of the holistic states. Even when the number of layer is one, our system still has holistic states, which have logical relationship with the states in basic layers through the canalizing functions. So it is still different from the one layer BNs.

Algebraic representation of BMCNs

In this section, we introduce some concepts and properties, changing our BMCNs into algebraic representation.

Concepts and properties of the semi-tensor product

In this subsection, some concepts and properties of the STP will be briefly introduced36.

Definition 1.36

Let X be a row vector of dimension np, and Y = [y₁, y₂,…, y_p]^T be a column vector of p dimension. Then we split X into p equal-size blocks as X¹, X², …, X^p, which are 1 × n rows. Define the STP, denoted by , as
Let . 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 A and B, denoted by , as the following: C consists of m × q blocks as C = (C^ij) and each block is

where Aⁱ is the i–th row of A and B_j is the j–th column of B.

And here we give some fundamental properties of the STP in the following36:

Lemma 1.36 Assume , then (where ⊗ is the Kronecker product, I_t is the identity matrix).

Assume , then

Lemma 2.36 Assume A ∈ M_m×n is given.

Let Z ∈ R^t be a row vector. Then
Let Z ∈ R^t be a column vector. Then

It is easy to find out that STP of matrix can be seen as a generalization of conventional matrix product. All the fundamental properties of conventional matrix product, such as distributive rule, associative rule, remain true. So we can omit the symbol .

Here we defined some notions for statement ease.

, where denotes the i–th column of the identify matrix I_n.
Assume a matrix , where are positive integer constants. We can briefly denoted it as .
A matrix is called a logical matrix if the columns of A, denoted by Col (A), satisfy . And the set of m × n logic matrices is denoted by .

Then we define a swap matrix , which is constructed in the following way: label its columns by (11, 12, …, 1n, …, m1, m2, …, mn) and its rows by (11, 21, …, m1, …, 1n, 2n, …, mn). And its element in the position ((I, J), (i, j)) is assigned as

And we denote W_[n] = W_{[n, n]} when m = n.

Lemma 3.36 Let and be two columns. Then

For the logical function with n arguments , we can convert it into an algebraic function using the STP of matrices. A logical domain, denoted by , is defined as . We identify each element in with a vector as , and . Based on this, we have

Lemma 4.36 Any logical function f(x₁, x₂, …, x_n) with logical arguments x₁, x₂, …, x_n ∈ Δ, can be expressed in a multi-linear form as

where is unique, which is called the structure matrix of logical function f.

And here we give another lemma:

Lemma 5.36Assume, then

where

Here M_r = δ₄[1, 4], which is power-reducing matrix and it can be verified that P² = M_rP, ∀P ∈ Δ.

Based on the above properties of STP, we put forward an obvious proposition.

Proposition 1.For each, if, we can find a matrix R_isuch that

Algebraic structure of the BMCNs

In this subsection, we change our BMCNs into discrete version using STP tool. To express it more clearly, here we give some description of variables.

means the state of layer l.
means the control.
means the state of all layers.

In the following step we will change the given BMCNs (3)-(4) into algebraic representation, as we will find out the algebraic relation between x(t + 1) and x(t) as well as the algebraic relation between and x(t).

At the first place, we will find out the algebraic relation between x(t + 1) and x(t). Using lemma 4 and proposition 1, for each logical rule , we can find its structure matrix , so we obtain that

where , and . Hence

with . And we have defined that , so we obtain that

where .

Subsequently, we will find out the algebraic relation between and x(t). Using the similar steps above, the algebraic representation of (4) can be obtained as

where is the structure matrix of logical function . And we have that . So we obtain that

where .

Means that

Similarly, by letting , we obtain the algebraic expression of the output dynamics (5) as follows:

where , here H_js are the structure matrixes of h_j, j = 1, 2, …, p.

Here we give an example to illustrate this process.

Example 1. Consider following two-layer BMCN, with N = 2, K = 2, n = 4 and m = 4

and we have that

where represent the logical functions of negation, disjunction, conjunction, implication, and equivalence, respectively. Based on Lemma 4, we obtain the corresponding structure matrices of those logical operators, as given in Table 1.

Define . Then we calculate the control-depending network transition matrix of system.

and we have that

Then

Here we have found out the algebraic relation between x(t + 1) and x(t) as well as the algebraic relation between and x(t). Furthermore, we assume that

Then, according to properties of STP, we obtain the matrix expression of output, as follows

Observability of BMCNs

In this section, we will analyze and characterize the observability of the BMCNs, with two different types of controls. We first provide some definitions as follows.

Consider the BMCN (3)-(4) with output dynamics (5). For any initial state and control input sequence , the holistic trajectory at time t is denoted by . Output trajectory at time t denote by y(; u, x(0)).

Definition 2. The BMCN (3)-(5) is observable if there exists a finite control sequence , with s > 0, such that for any with i ≠ j, we have for some .

In other words, there exists a control input sequence for which the initial state can be uniquely determined from the knowledge of the output sequence.

Remark 2. Our definition is motivated by the definition of observability for BCNs proposed in Laschov, et al.39, which is different from the one proposed by Cheng, D. et al.24. In Cheng, D. et al.24, a BCN is said to be observable if the initial state can be uniquely determined from the knowledge of the control inputs (which may depend on the initial state) and the outputs.

We consider two kinds of controls. The first is that the controls are determined by certain logical rules, which called the input networks.

where , are logical function.

According Lemma 4, we know that the input network (11) can be expressed as

where are the structure matrix of logical function g_i, respectively. Then,

with .

Theorem 1. Consider (3)-(5) (or equivalently (8)-(10)) with input network control (11). The system is observable if and only if there exists finite time s, s > 0, such that , for some , where

Proof. By considering the input network, put together (8)-(9) with (12), we can obtain the system

A straightforward computation shows, we calculate the output as follows

Hence, in the matrix form, we obtain

From the solution structure of the system of linear algebraic equations, we know that for some initial control input , the system of linear equations (14) with 2^NK-dimension unknown vector x(0) has a unique solution if and only if the system matrix has rank 2^NK. That is, for some initial control input , the initial state x(0) is uniquely determined by the knowledge of the output sequence if and only if

This completes the proof.

Remark 3. From the proof of above theorem, we obtain that for some if the matrix has full column rank, means that , then the initial state x(0) can be reconstructed by the left inverse of operation on output sequence ,

In the following, we consider the case when the controls are free Boolean sequence. Precisely, m controls are described as and freely designed.

Theorem 2. Consider (3)-(4) and (5) (or equivalently (8), (9) and (10)), with a free Boolean sequence control. The system is observable if and only if there exists a finite control sequence with such that , where

Proof. If the controls come from a free Boolean sequence, the system is that

If free control inputs are given, then a straightforward computation shows the following:

Hence, in the matrix form, we obtain

As a similar analysis discussed in the proof of Theorem 1, we know that for a given free control inputs , the system of linear equations (17) with 2^NK-dimension unknown vector x(0) has a unique solution if and only if the system matrix has rank 2^NK. That is, for a given free control inputs , the initial state x(0) is uniquely determined by the knowledge of the output sequence if and only if . Furthermore, as mentioned in Remark 3, the initial state x(0) can exactly calculated as This completes the proof.

Remark 4. The observability in our paper is the observability of which is the all initial states of all layers in the initial time. Boolean control network (3)-(4) is observable if for the initial state , there exists finite time , such that the initial state can be uniquely determined from the knowledge of the controls and the outputs . Based on the initial state , we can easily obtain the holistic states through the canalizing function . So the holistic states are also observable.

Examples

In this section, we will give some examples to illustrate our results. Example 2 is a observable case and Example 3 is an unobservable case.

Example 2. (Continue to Example 1) Consider the two-layer BMCN given in Example 1. Assume that the control inputs are determined by the following input network

Then we obtain that

If we take By calculation, we have

Then we have that

And we can obtain that . Then from Theorem 1, we know that the system is observable under the input network (18).

Example 3. Consider following two-layer BMCN, with N = 2, K = 2, n = 3 and m = 1

and we have that

Then we calculate the control-depending network transition matrix of system.

and we have that

Then

Furthermore, we assume that

Then, according to properties of STP, we obtain the matrix expression of output, as follows

Now, we analyze the observability of this system, based on Theorem 2. We can calculate that while ,

And while , we have that

Then, by induction, we easy obtain that, for any s > 0, and free control input , , and furthermore, . That implies, for any s > 0, and free control input

So the linear homogeneous equation

has the non-zero solution. Then we obtain that for arbitrary s > 0, we still have that , by Theorem 2, the system is unobservable.

Conclusions

In this paper, input controls were introduced into BMNs. By means of STP approach, the above logical dynamics has been converted into an algebraic form and the observability of dynamics is discussed. Firstly, we gave the theorem about the observability of whole dynamic system. Subsequently, the observability of each node in the special layer has been proved. Finally, we put forward some examples to illustrate our results.

Additional Information

How to cite this article: Wu, Y. et al. Observability of Boolean multiplex control networks. Sci. Rep.7, 46495; doi: 10.1038/srep46495 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figure 1

Schematic of multiple networks with two layers.

Here K = 2 means that the system is a two layers network. N = 4 means that there are four nodes in each layer. And n = 5 means that there are five total different nodes in the system. are the states in layer 1, and are the states in layer 2.

Figure 2

Schematic illustration of the relationship of node states in the fixed layers with the holistic states in BMNs.

Where, represent the holistic states of BMNs. For example is the holistic state of and . It is affected by and through canalizing function . The second node is only existed in layer one. So holistic state is only affected by .

Figure 3

Schematic illustration of BNs with control and output.

The inputs m dimension control have been introduced. denote outputs. From the figure, we see that inputs affect the node states in each layers as well as the abstract holistic states. And outputs are affected by holistic states .

Table 1

Structure matrices of some basic logical functions.

f(x₁, x₂)	¬x₁	x₁ ∨ x₂	x₁∧x₂	x₁ → x₂	x₁ ↔ x₂
M_f	M_n = δ₂[2, 1]	M_d = δ₂[1, 1, 1, 2]	M_c = δ₂[1, 2, 2, 2]	M_i = δ₂[1, 2, 1, 1]	M_e = δ₂[1, 2, 2, 1]

Footnotes

The authors declare no competing financial interests.

Author Contributions Wei Wang put forward the idea, Yuhu Wu and Jingxue Xu demonstrated the results, Xi-Ming Sun pointed examples. All authors reviewed the manuscript.

Acknowledgments

This research is supported by the National Natural Science Foundation of China under Grant 61325014 and 11401142.

References

1. HoodL. & RowenL.The human genome project: big science transforms biology and medicine. Genome Medicine5, 1–8 (2012).[Google Scholar]
2. KitanoH.Systems biology: A brief overview. Science295, 1662–1664 (2002).[PubMed][Google Scholar]
3. ZhongJ., LuJ., LiuY. & CaoJ.Synchronization in an array of output-coupled boolean networks with time delay. IEEE Transactions on Neural Networks & Learning Systems25, 2288–2294 (2014).[PubMed][Google Scholar]
4. KauffmanS. A.Metabolic stability and epigenesis in randomly constructed genetic nets. Journal of Theoretical Biology22, 437–467 (1969).[PubMed][Google Scholar]
5. KauffmanS. A. & (OUP), O. U. P. Origins of order: self-organization and selection in evolution. Journal of Evolutionary Biology13, 133–144 (1993).[Google Scholar]
6. KauffmanS. A.At home in the universe. Mathematical Social Sciences33, 94–95 (1995).[Google Scholar]
7. AldanaM.Boolean dynamics of networks with scale-free topology. Physica D Nonlinear Phenomena185, 45–66 (2003).[Google Scholar]
8. HeidelJ., MaloneyJ., FarrowC. & RogersJ. A.Finding cycles in synchronous boolean networks with applications to biochemical systems. International Journal of Bifurcation & Chaos13, 535–552 (2011).[Google Scholar]
9. FarrowC., HeidelJ., MaloneyJ. & RogersJ.Scalar equations for synchronous boolean networks with biological applications. IEEE Transactions on Neural Networks15, 348–354 (2004).[PubMed][Google Scholar]
10. HuangS. & IngberD. E.Shape-dependent control of cell growth, differentiation, and apoptosis: Switching between attractors in cell regulatory networks. Experimental Cell Research261, 91–103 (2000).[PubMed][Google Scholar]
11. HuangS.Regulation of cellular states in mammalian cells from a genomewide view. In Gene Regulations and Metabolism - Postgenomic Computational Approaches181–220 (2002).[Google Scholar]
12. AkutsuT., MiyanoS. & KuharaS.Inferring qualitative relations in genetic networks and metabolic pathways. Bioinformatics16, 727–734 (2000).[PubMed][Google Scholar]
13. AlbertR. & BarabasiA. L.Dynamics of complex systems: scaling laws for the period of boolean networks. Physical Review Letters84, 5660–5663 (2000).[PubMed][Google Scholar]
14. MengM. & FengJ. E.Synchronization of interconnected multi-valued logical networks. In Chinese Control Conference1659–1669 (2013).[Google Scholar]
15. WuY. & ShenT.An algebraic expression of finite horizon optimal control algorithm for stochastic logical dynamical systems. Systems & Control Letters82, 108–114 (2015).[Google Scholar]
16. VillegasP., Ruiz-FrancoJ., HidalgoJ. & MuñozM. A.Intrinsic noise and deviations from criticality in boolean gene-regulatory networks. Scientific Reports6, 34743 (2016).[PubMed][Google Scholar]
17. ChenH., WangG., SimhaR., DuC. & ChenZ.Boolean models of biological processes explain cascade-like behavior. Scientific Reports7, 20067 (2016).[PubMed][Google Scholar]
18. LuJ., ZhongJ., LiL., HoD. W. & CaoJ.Synchronization analysis of master-slave probabilistic boolean networks. Scientific reports5, 13437 (2015).[PubMed][Google Scholar]
19. ChenH. & SunJ.A new approach for global controllability of higher order boolean control network. Neural Networks the Official Journal of the International Neural Network Society39, 12–17 (2013).[PubMed][Google Scholar]
20. LiF. & SunJ.Controllability of higher order boolean control networks. Applied Mathematics & Computation219, 158–169 (2012).[Google Scholar]
21. LuJ., ZhongJ., HoD. W. C., TangY. & CaoJ.On controllability of delayed boolean control networks. SIAM Journal on Control & Optimization54, 475–494 (2016).[Google Scholar]
22. LuJ., ZhongJ., HuangC. & CaoJ.On pinning controllability of boolean control networks. IEEE Transactions on Automatic Control61, 1658–1663 (2016).[Google Scholar]
23. ChengD., LiZ. & QiH.Realization of boolean control networks. Automatica (Journal of IFAC)46, 62–69 (2010).[Google Scholar]
24. ChengD., QiH. & LiZ.Controllability and observability of boolean control networks. Automatica45, 1659–1667 (2009).[Google Scholar]
25. LiF., SunJ. & WuQ. D.Observability of boolean control networks with state time delays. IEEE Transactions on Neural Networks22, 948–954 (2011).[PubMed][Google Scholar]
26. ZhangL. & ZhangK.Controllability and observability of boolean control networks with time-variant delays in states. IEEE Transactions on Neural Networks & Learning Systems24, 1478–1484 (2013).[PubMed][Google Scholar]
27. LuoC., WangX. & LiuH.Controllability of time-delayed boolean multiplex control networks under asynchronous stochastic update. Scientific Reports4, 7522 (2014).[PubMed][Google Scholar]
28. ZhongJ., LuJ., HuangT. & HoD. W. C.Controllability and synchronization analysis of identical-hierarchy mixed-valued logical control networks. IEEE Transactions on Cybernetics, doi: (2016).[PubMed][Google Scholar]
29. IdekerT., GalitskiT. & HoodL.A new approach to decoding life: systems biology. Annual Review of Genomics & Human Genetics2, 343–372 (2003).[Google Scholar]
30. KadelkaC., MurrugarraD. & LaubenbacherR.Stabilizing gene regulatory networks through feedforward loops. Chaos: An Interdisciplinary Journal of Nonlinear Science23, 025107 (2013).[Google Scholar]
31. ChavesM., SontagE. D. & AlbertR.Methods of robustness analysis for boolean models of gene control networks. Systems Biology153, 154–167 (2006).[PubMed][Google Scholar]
32. FaryabiB., DattaA. & DoughertyE. R.On approximate stochastic control in genetic regulatory networks. IET Systems Biology1, 361–368 (2007).[PubMed][Google Scholar]
33. ChengD., ZhaoY. & XuT.Receding horizon based feedback optimization for mix-valued logical networks. IEEE Transactions on Automatic Control60, 3362–3366 (2015).[Google Scholar]
34. WuY., KumarM. & ShenT.A stochastic logical system approach to model and optimal control of cyclic variation of residual gas fraction in combustion engines. Applied Thermal Engineering93, 251–259 (2016).[Google Scholar]
35. WuY. & ShenT.Policy iteration approach to control residual gas fraction in ic engines under the framework of stochastic logical dynamics. IEEE Transactions on Control Systems Technology25, 1100–1107 (2017).[Google Scholar]
36. ChengD., QiH. & ZhaoY.An Introduction to Semi-Tensor Product of Matrices and Its Applications (World Scientific Publishing Co. Pte. Ltd., 2012).
37. CobelliC. & Romanin-JacurG.Controllability, observability and structural identifiability of multi input and multi output biological compartmental systems. IEEE Transactions on Biomedical Engineering23, 93–100 (1976).[PubMed][Google Scholar]
38. LopezI., GamezM. & CarrenoR.Observability in dynamic evolutionary models. Biosystems73, 99–109 (2004).[PubMed][Google Scholar]
39. LaschovD., MargaliotM. & EvenG.Observability of boolean networks: A graph-theoretic approach. Automatica49, 2351–2362 (2013).[Google Scholar]
40. ZhangK. & ZhangL.Observability of boolean control networks: A unified approach based on finite automata. IEEE Transactions on Automatic Control61, 2733–2738 (2016).[Google Scholar]
41. KitanoH.Systems biology: a brief overview. Science295, 1662–1664 (2002).[PubMed][Google Scholar]
42. MuchaP. J. & OnnelaJ. P.Community structure in time-dependent, multiscale, and multiplex networks. Science328, 876–878 (2010).[PubMed][Google Scholar]
43. CozzoE., ArenasA. & MorenoY.Stability of boolean multilevel networks. Physical Review E Statistical Nonlinear & Soft Matter Physics86, 2569–2575 (2012).[Google Scholar]
44. GoffmanE.Frame Analysis. An Essay on the Organization of Experience (PENGUIN, 1975).
45. ParshaniR., RozenblatC., IetriD., DucruetC. & HavlinS.Inter-similarity between coupled networks. EPL (Europhysics Letters)92, 68002 (2010).[Google Scholar]
46. XuM., ZhouJ., LuJ.-a. & WuX.Synchronizability of two-layer networks. The European Physical Journal B88, 1–6 (2015).[Google Scholar]
47. ZhongJ., HoD. W. C., LuJ. & XuW.Controllability for a special case of multi-level boolean control networks. In 2016 IEEE International Conference on Industrial Technology (ICIT)1378–1383 (2016).[Google Scholar]