Abstract

It is worth noting that both nodes’ coupling connections and logical updating functions play a vital role in state evolutions of Boolean networks (BNs). In this paper, a new concept named structural controllability (SC) about Boolean control networks (BCNs) with known partial information on nodes’ connections is studied. Then, by referring to semi-tensor product (STP) techniques, several types of SC are presented according to different issues of Boolean functions. Thereafter, several necessary and sufficient conditions are derived for SC of BCNs. Finally, a biological model of the lactose operon in Escherichia coli is simulated to show the effectiveness of the main theoretical results.

1. Introduction

Boolean networks (BNs) are known as discrete-time logical dynamics, firstly proposed in [1], where nodes evolute the corresponding states according to several Boolean functions. In a BN, every gene can be chosen from a set with logical variables 1 and 0. Therefore, a BN is illustrated according to certain Boolean functions . A digraph (named as network structure digraph) with nodes can well show the connections among every node (named as network structure), where there exists a directed edge from node to node if depends on , i.e., if there is a sequence of variables satisfying , where . Recently, many important properties of BNs have been issued [218], and so on.

Recently, great attention has been paid to BN owing to a new concept called semi-tensor product (STP) of matrices [19]. STP is firstly proposed by Cheng and his collaborators [19], which is a generalization of conventional matrix product. STP can be used as an efficient tool to convert BNs into corresponding algebraic representations. Then, many fundamental results about BNs have been issued [2022]. Please refer to [23, 24] for more applications of STP.

It is worth noting that stability and controllability are of importance in different research fields, such as impulsive differential systems [25, 26], functional differential equations [27], switched systems, and Boolean networks [2830]. If we take external control inputs into consideration, BNs will be Boolean control networks (BCNs) established by Akutsu [28] for the first time. Note that it has been increasingly challenging for researches concerning BNs and BCNs [29, 30]. Recently, the problem of controllability has been studied rapidly; some fundamental results have been established. To list a few, one of the main results on controllability has been issued [31], based on an input-state incidence matrix, where controllability is illustrated by the positiveness of a controllability matrix. Moreover, by Perron–Frobenius theory, controllability has been studied in [20]. In the past few years, another important issue (called pinning control method) gets great attention with less control cost and lower computational complexity [3234].

It is noted that both nodes’ coupling connections and logical updating functions play a vital role in state evolutions of BNs. During the past few decades, a relationship between nodes’ connection and some concepts of BNs, such as stability and oscillatority, has been issued in [3539]. Note that in real-world, not all the messages about nodes coupling and the updating Boolean functions are available, that is, only part of the nodes’ connections are available. For instance, a new research concept named structural controllability (SC) concerning BCNs is well addressed in [40], whose full messages of nodes’ connections are available while Boolean functions are unavailable. Therefore, it is important to analyze dynamical properties with only part of the messages on nodes’ connections, i.e., all the messages of BNs are unavailable.

Therefore, based on the above analysis, some contributions of this paper are given: (i) SC of BCNs with partial information of nodes’ connections available has been firstly issued, where partial nodes’ connections are available; (ii) four cases of SC are defined owing to different types of network structure and logical functions. Then, several types of SC matrices are constructed using STP and Hadamard product, and some criteria are issued; (iii) SC of BCNs with partial information of nodes’ connections available is a simple extension of general controllability of BCNs with full information of both network structure and logical functions.

The rest of the paper is organized as follows: in Section 2, some preliminaries are given. In Section 3, the main results of SC are given, and some criteria are obtained. In Section 4, to explain the main results, one biological model is provided. In Section 5, a brief conclusion is established.

Compared with the conference version of this paper [41], in this paper, much more detailed descriptions of the theoretical results about SC have been established. In [41], only two cases of SC are considered. However, in the revision, totally four cases of SC are discussed, and several necessary and sufficient conditions are derived for SC of BCNs, which is a generalization of the theoretical results in [41]. In addition, in this paper, some comments and detailed comparisons with the existing references have been presented, as well as a new figure has been added to better illustrate the implication relationships between four different types of SC.

Here, some notations are given: denotes the set consisting of positive integers. , where and . is the -th column of the identity matrix . . denotes the -th column (columns) of matrix . denotes the set of real matrices (logical matrices). is the set consisting of Boolean matrices. is a logical matrix. Boolean addition is defined as , where . , where , . (A) represents the -th block of matrix A. A Boolean matrix satisfies all the entries belonging to . .” denotes the known Kronecker product, “” denotes the known Khatri–Rao product, “” denotes the known Hadamard product, and “” denotes STP of two given matrices.

2. Preliminaries

Firstly, here, some basic illustrations on STP and BCNs are introduced.

Definition 1. The STP of two given matrices and is defined as , where is the least common multiple of and .

Define the bijective relationship “” between and as , which naturally extends to between and . Equivalently, Boolean functions can be transferred into an algebraic representation.

Lemma 1. Given a logical function , its multilinear form is given as , where is the structure matrix, uniquely determined by .

Some well-known structure matrices of binary operators are presented based on Lemma 1, such as negation “, conjunction “, disjunction “, conditional “, and biconditional “.

Normally, a BCN consisting of nodes and external inputs can be described:where the negation of is denoted by , i.e., . Here, express the states, expresses the control input, and and are the index sets of the in-neighbors of -th node. For each , we call sets as the activating coupling node sets of the -th node. In addition, for each , and are the index sets of coupling nodes of -th node, and we call sets as the inhibiting couplings node sets of the -th node. , , are Boolean functions, which only consist of “” and “.” Therefore, the logical dynamics of (1) is only connected by logical operators “,” “,” and “.” Thus, given a BCN (1), one can always obtain its corresponding algebraic representation using Lemma 1.

2.1. Problem Formulation

It is noted that dynamical evolutions of BN (1) are determined by the Boolean functions , as well as its coupling sets , , , and , . In practical biological applications, due to experimental impacts and other reasons, we might only know the in-neighbors of some nodes, but cannot acquire the coupling relationship for each node. Therefore, assume that only partial information of (1) is known in this paper. Here, we assume that the sets , and , , are given beforehand, while binary operators among different nodes in , , , and , , are unavailable, which implies that for nodes , the corresponding activating coupling nodes, and inhibiting coupling nodes are given, whereas binary operators among nodes are unavailable.

For simple illustrations, suppose only activating coupling nodes and inhibiting coupling nodes for nodes , are given, which is as follows:where , are known and , and , are unknown. The binary operators among all the nodes are unknown, which implies that the logical operators between nodes in , and , are unknown.

Based on Lemma 1, one has equivalent algebraic form of logical dynamics of (2). Given sets , and , , we make the following ordering assumptions to simplify calculation:

Then, under assumption (3), one has thatwhere binary operators , , , , , and . Here, operators , , , and are binary operators connecting the in-neighbor sets , and .

Given binary operators , , and , , one has thatwhere matrix is uniquely calculated according to binary operators , , and , . Thus, there totally exist different kinds of structure matrices for matrices , . Given each , illustrate all the kinds of matrices in (5) by set .

Denoting and , one can further calculate a unique matrix , , based on Lemma 1, satisfying

Thus, for , the following equation is obtained:

Let ; it leads to a simplification illustration of (4) for nodes , using Khatri–Rao product,where . Then, one has . Denoting , assume that , where matrices .

Then, let ; one can obtain algebraic form for the last nodes:where expresses the structure matrix determined by in-neighbor sets , and , , and its binary operators. Therefore, it leads to totally cases for matrix . Denote the set of all the different kinds of by , that is, , where . Therefore, considering (2), one has thatwhere and . Multiplying (10) yields an algebraic form of (2),where matrix is called the state transition matrix of system (2) choosing from a set . In addition, one can calculate that and represent all of the state transition matrices by , where . This implies that each matrix , , is one possible state transition matrix for system (2), where only the corresponding in-neighbor sets , and , , for those first nodes, are known.

3. Main Results

In this section, several types of controllability of (2) are presented, where only the corresponding in-neighbors of nodes , are known.

3.1. First Class of Structural Controllability

At the beginning, only the corresponding in-neighbor sets for , are supposed to be known, the corresponding Boolean functions , , and both the in-neighbor sets and the Boolean functions , , are unknown. Under these assumptions, reachability and controllability for (2) are firstly defined as follows.

Firstly, consider the first class structural reachability (FCSR) between two given states and under any Boolean functions , .

Definition 2. Consider system (2) and the neighbor sets , and for , are available. The destination state achieves FCSR from the initial state , if can be steered to for arbitrary Boolean functions . In addition, (2) achieves first class structural controllability (FCSC) if for any initial state and destination state , achieves FCSR from .

By considering , divide , as follows: , and denotewhere , and matrix , , is the controllability matrix when state transition matrix of system (2) is . Then, introduce the following Boolean matrix in order to issue FCSC,where “” represents the well-known Hadamard matrix product. Here, name as FCSC matrix of (2) with partial in-neighbor sets , and for nodes available, .

By referring to matrix established in (13), the following criteria for FCSC are established.

Theorem 1. Consider system (2) with available sets , and for nodes , then state achieves FCSR from state if and only if the following condition holds.

Proof. According to the definition of Hadamard product “,” one can easily obtain that if and only if for each , . For system , , according to [31], one has that is reachable from if and only if . Thus, by Definition 2, achieves FCSC from if and only if condition holds.
Then, FCSC for (2) with available partial in-neighbor sets , and , , is established.

Theorem 2. Given system (2) with available neighbor sets , and for nodes , (2) achieves FCSC if and only if condition holds.

Proof. From Theorem 1, for any initial state and destination state , achieves FCSR from if and only if . Thus, according to Definition 2, system (2) achieves FCSC.

3.2. Second Class of Structural Controllability

Based on Definition 2, FCSR is introduced on reachability between the initial state and destination state for any Boolean functions . Hence, assume that state can be steered to state under certain possible Boolean mappings , then the following second class of structural controllability is established.

Definition 3. Given system (2) and the neighbor sets , and , for nodes , are available, the destination state achieves second class structural reachability (SCSR) from the initial state , if there is certain Boolean functions satisfying that can be reachable from . In addition, (2) achieves SCSC if for any , achieves second class structural controllability (SCSR) from .

By referring to established matrices in (12), , define

We name as the SCSC matrix of (2) with available sets , and for nodes , . According to defined in (14), the following necessary and sufficient criteria are derived.

Theorem 3. Consider system (2) with available neighbor sets , and for nodes , then state achieves SCSR from state if and only if condition holds.

Proof. if and only if there is a sequence of Boolean matrices , , such that . For system , according to [31], one has that can be reachable from if and only if . Thus, by Definition 3, achieves SCSR from if and only if condition holds.

Theorem 4. Considering system (2), the neighbor sets , and for nodes , , are available, then (2) achieves SCSC if and only if condition holds.

Proof. From Theorem 3, for any two states , state achieves SCSR from state if and only if condition holds. Thus, according to Definition 3, system (2) achieves SCSC.
In the above, controllability of (2) with known in-neighbor sets , and for nodes , under the situations of the arbitrariness and existence of logical functions , has been established, which are shown in Theorems 1 and 3. In the following, controllability of (2) under the situation of arbitrary logical functions , , and certain possible logical functions , , will be studied.

3.3. Third Class of Structural Controllability

In the following sequel, we consider the reachability between two given states and under any Boolean functions , , and certain possible logical functions , .

Definition 4. Considering system (2) and the neighbor sets , and , for nodes , are available, the destination state achieves third class structural reachability (TCSR) from the initial state , if for any Boolean functions , there are certain Boolean mappings satisfying that can be reachable from . In addition, (2) achieves third class structural controllability (TCSC), if for any , achieves TCSR from .

For convenience, denote the set of in (11) by , where , , . Divide into equal blocks as , and denotewhere , and matrix is the corresponding controllability matrix when state transition matrix of system (2) is . In addition, introduce the following Boolean matrix in order to issue TCSC,where , . We name as TCSC matrix of (2) with available sets , and for .

By referring to matrix defined in (16), the following criteria for TCSC are established.

Theorem 5. Given system (2), for , the neighbor sets , and for nodes are available, then state achieves LCSR from state if and only if condition holds.

Proof. if and only if , and there is a sequence of Boolean matrices , such that for every . For system , according to [31], one has that is reachable from if and only if . Thus, by Definition 4, state achieves TCSR from state if and only if condition holds.

Theorem 6. Given system (2), the neighbor sets , and for nodes , , are available, then (2) achieves TCSC if and only if condition holds.

Proof. From Theorem 5, for any two states , state achieves TCSR from state if and only if condition holds. Thus, according to Definition 4, system (2) achieves TCSC.

Remark 1. More strictly, if for , , logical functions are the same, then for any functions , we define the TCSC matrix as , where .

3.4. Last Class of Structural Controllability

Based on Definition 4, TCSR is introduced on reachability between the initial state and destination state for any Boolean functions and certain Boolean functions , . Thus, if destination state is reachable from the initial state under certain Boolean functions , and any Boolean functions , the following last class of structural controllability is established.

Definition 5. Considering system (2), for , the neighbor sets , and for nodes are available. The destination state achieves last class structural reachability (LCSR) from initial state , if there is certain Boolean functions satisfying that can be reachable from for any Boolean functions . Moreover, (2) achieves last class structural controllability (LCSC) if for any , state achieves LCSR from state .

By referring to matrix established in (15), , definewhere . We call as LCSC matrix of (2) with available sets , and for nodes .

By referring to matrix defined in (17), the following criteria for LCSC are obtained.

Theorem 7. Given system (2), for , the neighbor sets , and for nodes are available, then the destination state achieves LCSR from the initial state if and only if condition holds.

Proof. if and only if there exists a sequence of matrices , such that , and for every matrix , . For system , according to [31], one has that is reachable from if and only if . Thus, by Definition 5, achieves LCSR from if and only if condition holds.

Theorem 8. Given system (2), for , the neighbor sets , and for nodes are available, then (2) achieves LCSC if and only if condition holds.

Proof. From Theorem 7, for any initial state and destination state , state achieves LCSR from state if and only if condition holds. Thus, according to Definition 5, system (2) achieves LCSC.

Remark 2. In this paper, four types of SC of BCNs with partial information available have been studied, where one only knows part of messages on nodes’ connections. Among the abovementioned four different kinds of structural controllability, the relationships are shown as follows: FCSC implies TCSC and LCSC; the third and the last class of SC both imply SCSC. In addition, FCSC is the most conservative definition and all other types of SC can be implied by FCSC, and the diagram describing the detailed relationships among different kinds of SC is shown as Figure 1.

Remark 3. By comparing with the recent theoretical results [39, 40], the established theoretical results in this paper are more general. If , that is, one knows all the in-neighbors of (2), then the results in [40] can be directly obtained from FCSC in this paper. Moreover, if , FCSC in this paper implies common controllability [39], which means that system (2) is always controllable no matter which functions are chosen. Thus, the study of these four cases of SC of BCNs with partial information available is a simple extension of common controllability, which provides a new insight for further studying and reducing the high computational complexity of BCNs.

4. Simulations

Here, simulations on a Boolean model of the lactose operon in Escherichia coli are given to illustrate the obtained results.

Example 1. Firstly, a lactose operon model in Escherichia coli modeled by a BCN is considered [42]:where denotes the biological structure of mRNA, denotes the biological structure of lacZ polypeptide, denotes the intracellular lactose, denotes the external glucose, and denotes the external lactose.
Using STP tool, the corresponding algebraic form of system (18) is given as follows: , where . Then, split matrix into four blocks and obtain its controllability matrix , where . After a straightforward calculation, one has that . By referring to [31], one can conclude that system (18) is not controllable.
Suppose that only in-neighbors for nodes and are given, but the detailed logical functions are unknown. Thus, the following system is considered:where and express the unavailable sets of the activating indegrees (inhibiting indegrees) of node .
Based on equations (5)–(9), a simple equivalent algebraic illustration of (19) can be derived:where , , , , matrix , and the cardinality of is . Here, , , , and . Thus, one can obtain its equivalent augmented system of (19): , where matrix . Due to the high computational complexity of the set , the detailed information for matrices is omitted.
A straightforward computation shows that . By referring to Theorems 1, 3, 5, and 7, system (19) achieves SCSC and TCSC, while it cannot achieve FCSC or LCSC.
When , if , then , ; if , then , ; if , then , ; if , then , . Thus, . By referring to Theorems 3 and 5, system (19) achieves SCSC and TCSC. This implies that, for any , and any logical functions , , there is a logical functions satisfying that state can be reached from initial state . According to Definitions 3 and 4, one can conclude that system (19) achieves SCSC and TCSC, respectively.
When , that is, , if , then , ; if , then , ; if , then , ; if , then , . Thus, . According to Theorem 7, system (19) cannot achieve LCSC. In this case, there exists a logical function , such that system (19) is not controllable for any logical functions , . By Definition 5, system (19) cannot achieve LCSC.

5. Conclusion

In this paper, SC of BCNs has been systematically addressed, where only partial information of nodes’ connections is available. By referring to its equivalent algebraic representation of BCNs, four kinds of SC have been established, based on different issues of nodes’ connections. Then, certain criteria for different types of SC have been proposed, by defining its corresponding SC matrix. Finally, a Boolean model of the lactose operon in Escherichia coli has been simulated to show the effectiveness of the main theoretical results. One interesting research issue for future works includes SC of probabilistic BCNs with time-variant information of nodes’ connections.

Data Availability

No data were used to support this study.

Disclosure

This paper was presented in part at the 11th International Conference on Information Science and Technology (ICIST), Chengdu, China, May 21–23, 2021.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was jointly supported by the National Natural Science Foundation of China (61903339), the National Training Programs of Innovation and Entrepreneurship (202010345009), the Guangdong Basic and Applied Basic Research Foundation (2019A1515110234), and the Shenzhen Science and Technology Program (RCBS20200714114921371).