Home | About | Journals | Submit | Contact Us | Français |

**|**Scientific Reports**|**PMC4834509

Formats

Article sections

- Abstract
- Quantifying Control in Boolean Networks
- Control Portraits of Complex Systems
- Discussion
- Additional Information
- Supplementary Material
- References

Authors

Related links

Sci Rep. 2016; 6: 24456.

Published online 2016 April 18. doi: 10.1038/srep24456

PMCID: PMC4834509

Received 2015 October 27; Accepted 2016 March 21.

Copyright © 2016, Macmillan Publishers Limited

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

This article has been cited by other articles in PMC.

The study of network structure has uncovered signatures of the organization of complex systems. However, there is also a need to understand how to control them; for example, identifying strategies to revert a diseased cell to a healthy state, or a mature cell to a pluripotent state. Two recent methodologies suggest that the controllability of complex systems can be predicted solely from the graph of interactions between variables, without considering their dynamics: structural controllability and minimum dominating sets. We demonstrate that such structure-only methods fail to characterize controllability when dynamics are introduced. We study Boolean network ensembles of network motifs as well as three models of biochemical regulation: the segment polarity network in *Drosophila melanogaster*, the cell cycle of budding yeast *Saccharomyces cerevisiae*, and the floral organ arrangement in *Arabidopsis thaliana*. We demonstrate that structure-only methods both undershoot and overshoot the number and which sets of critical variables best control the dynamics of these models, highlighting the importance of the actual system dynamics in determining control. Our analysis further shows that the logic of automata transition functions, namely how canalizing they are, plays an important role in the extent to which structure predicts dynamics.

Complex systems are typically understood as large nonlinear systems. Their organization and behavior can be modeled by representations such as graphs and collections of automata. Graphs are useful to capture the *structure* of interactions between variables: the static organization of complex systems. However, nodes representing variables in graphs lack intrinsic dynamics. The simplest way to study nonlinear *dynamics* is to allow network nodes to have discrete states and update them with automata; for instance, Boolean Networks (BNs) are canonical models of complex systems which exhibit a wide range of interesting behaviors^{1}.

The study of network structure has uncovered several organizing principles of complex systems — such as scale-free networks and community structure — and how they constrain system behavior, without explicit dynamical rules for node variables^{2}. There is, however, a need to *control* complex systems, in addition to characterizing their organization. This is particularly true in systems biology and medicine, where increasingly accurate models of biochemical regulation have been produced^{3}^{,4}^{,5}^{,6}. More than understanding the organization of biochemical regulation, we need to derive control strategies that allow us, for instance, to revert a mutant cell to a wild-type state^{7}, or a mature cell to a pluripotent state^{8}. While the identification of such control strategies occurs for a given model, not the real system, predictions from control theory can be used for model verification and thus also aid the separate question of the accuracy of that model in predicting the real system.

Network structure has been reported to predict properties of dynamics, such as the synchronization of connected limit-cycle oscillators^{9}, or the likelihood of robust attractors^{10}. On the other hand, there are important system attributes which depend on dynamical characteristics of variables and their interactions; e.g. the critical transition between ordered and chaotic dynamics in BNs depends both on structural (mean connectivity) and dynamical properties of nodes (bias and canalization)^{11}^{,12}^{,13}^{,14}. Indeed, we already know that such dynamical properties strongly impact the stability, robustness, and controllability of existing models of gene regulation and biochemical signaling in a number of organisms^{7}^{,15}^{,16}^{,17}^{,18}. Therefore, a question of central importance remains: *How well does network structure predict the dynamics of the underlying complex system, especially from the viewpoint of control?*

Recently, two related methodologies were used to predict the controllability of complex networks based solely on network structure without consideration of the dynamical properties of variables: *structural controllability* (SC)^{19}^{,20} and *minimum dominating set* (MDS)^{21}^{,22}. Both techniques reduce dynamical systems to graphs where edges denote an interaction between a pair of variables. Using only graph connectivity, the goal is to identify a minimal set of *driver variables* (a.k.a. driver nodes) which can fully control system dynamics^{23}.

SC assumes that, in the absence of cycles, a variable can control at most one of its neighbors in the structural interaction graph^{19}^{,20}. The influence from an intervention on a node then propagates along a backbone of directed paths, where the number of necessary paths to cover the network dictates the minimum set of driver variables (see Supplemental Material, SM). Cycles are considered to be self-regulatory and do not require an external control signal. SC has become an influential method, having been used to suggest that biological systems are harder to control and have appreciably different control profiles than social or technological systems^{24}^{,25}. The methodology has also been used to identify key banks in interbank lending networks^{26}, and to relate circular network motifs to control in transcription regulatory networks^{27}. However, despite its successful characterization of observability (a dual notion to controllability) in several nonlinear dynamical systems^{28}, SC’s application to models of biological and social systems has been heavily critiqued due to its stringent assumptions^{29}^{,30}^{,31}.

MDS starts from the different assumption that each node can influence all of its neighbors simultaneously, but this signal cannot propagate any further. Driver variables are then identified by the minimal set such that every variable is separated by at most one interaction^{21}^{,22}. It has been used to identify control variables in protein interaction networks^{32} and characterize how disease genes perturb the Human regulatory network^{33}.

Because both MDS and SC use only the interaction graph of complex systems, unless otherwise specified, we use *structural control* to refer to both methods. Since these methods are increasingly used in a variety of scientific domains, it is important to study how much network structure predicts the controllability of realistic, nonlinear dynamical systems.

Here, we explore this problem using ensembles of BNs. These canonical models of complex systems are defined by a network of interconnected automata (the structure), and exhibit a wide range of dynamical behaviors^{1}. They have been used to model biochemical regulation in organisms, where dynamical attractors represent cell types, disease and healthy states^{6}^{,34}. It is well known that when the set of system variables is large, enumeration of the state-spaces of BNs becomes difficult, making the control problem for general deterministic BNs computationally intractable (NP-hard)^{35}. However, for small systems we can fully enumerate the state-space and compute the actual controllability (as measured by three proposed measures of controllability) for parameterized ensembles of BNs.

Our analysis is not meant to introduce alternative techniques to uncover control variables in BNs, since methods based on system dynamics already exist^{7}^{,8}^{,36}^{,37}^{,38}^{,39}^{,40}^{,41}^{,42}. The goal is to quantify the discrepancy between control as uncovered by approximate methods that use structure alone, from how actual control unfolds in BNs. Additionally, we characterize critical variables for the control of three models of biochemical regulation: the single-cell segment polarity network in *Drosophila melanogaster*, the eukaryotic cell cycle of budding yeast *Saccharomyces cerevisiae*, and the floral organ arrangement in the flowering plant *Arabidopsis thaliana*. Our results demonstrate that network structure is not sufficient to characterize the controllability of complex systems; predictions based on structural control can both under- and over-estimate the number and set of necessary driver variables. Therefore, previous assertions about the controllability of biochemical systems reached from analyses based on structural control methods do not offer a realistic portrayal of control^{24}^{,25}.

Boolean Networks (BNs) are discrete dynamical systems *X*{*x*_{i}} of *N* Boolean variables *x*_{i}{0,1}. Interactions between variables are represented as a directed adjacency graph, the *structural network*: *G*=(*X*,*E*), where edges *e*_{ji}*E* denote that variable *x*_{j} is an input to variable *x*_{i}. Furthermore, *X*_{i}{*x*_{j}*X*:*e*_{ji}*E*} and |*X*_{i}|=*k*_{i} denote the input set and the *in-degree* of variable *x*_{i}, respectively. Here, variables are updated synchronously according to deterministic logical functions: , such that , where denotes the state of the inputs to x_i at time .

At time *t*, the network is in a *configuration* of states **X**^{t}, which is a vector of all variable states at *t*. The set of all possible network configurations is denoted by , where . The complete dynamical behavior of the system for all initial conditions is captured by the *state-transition graph* (STG): , where each node is a configuration , and an edge denotes that a system in configuration **X**_{α} at time *t* will be in configuration **X**_{β} at time *t*+1. Under deterministic dynamics, only a single transition edge *T*_{α,β} is allowed out of every configuration node **X**_{α}. Because is finite, it contains at least one attractor, as some configuration or cycle of configurations must repeat in time^{43}. An exemplar STG is shown in Fig. 1 (top, left).

We study the control exerted on the dynamics of a BN by a subset of *driver variables D**X*. Here, control *interventions* are instantaneous bit-flip perturbations to the state of the variables in *D*^{44}. To capture all possible trajectories due to controlled interventions on *D*, we introduce the *controlled state transition graph* (CSTG): . The CSTG is an extension of the STG, where a set of additional edges denotes transitions from every configuration to each of its possible 2^{|D|}−1 perturbed counterparts. In Fig. 1, three examples of CSTG are shown with interventions to only one of the three variables: *D*={*x*_{1}},{*x*_{2}},{*x*_{3}}.

From the point of view of control theory^{19}^{,45}, the dynamics of a network of variables *X* is *controllable* by interventions to a subset of driver variables *D**X* when every configuration is reachable from every other configuration in . A configuration **X**_{β} is reachable from **X**_{α} if a directed path from **X**_{α} to **X**_{β} exists^{45}. For BN this is equivalent to requiring that the CSTG be strongly connected. To measure how much control *D* can exert, we tally the fraction of configurations that are reachable by interventions to *D*. Given a configuration **X**_{α}, the *fraction of reachable configurations*
is the number of other configurations **X**_{β} lying on all directed paths from **X**_{α}, normalized by the total number of other configurations 2^{N−1}. The *mean fraction of reachable configurations*:

measures the proportion of configurations which are on average reachable by controlling the set of driver variables *D*. When a network is fully controlled by *D*, , but for partially controlled networks .

Notice that , because the STG of a network (*D*) naturally contains transitions between configurations. Therefore, it is useful to measure the control exerted by a set of driver variables *D* beyond the uncontrolled dynamics. To this end, we introduce the *mean fraction of controlled configurations*:

It measures the fraction of configurations which are on average reachable by controlling the driver variables *D* that were not already reachable via the natural dynamics. By definition, for any system and set of driver variables.

In practice, only certain subsets of configurations are meaningful. These subsets are typically cast as either attractors for the system dynamics or specific trajectories through the state space. Consider the case of BNs as models of biochemical regulation; attractors represent different cell types^{1}^{,29}^{,46}, diseased or normal conditions^{47}, and wild-type or mutant phenotypes^{7}. In this context, the formal sense of controllability is well beyond what is necessary. What is most relevant for some systems is to uncover the driver variables which can steer dynamics from attractor to attractor; transient configurations are irrelevant.

To measure this more realistic sense of control, we introduce the *controlled attractor graph* (CAG): . In this graph, each node represents an attractor. A *basin edge*
, denotes the existence of at least one path from attractor **A**_{κ} to attractor **A**_{γ}. In Fig. 1 (right-side), three examples of CAGs are shown. The *mean fraction of reachable attractors* is then given by

where . It measures the fraction of attractors which are on average reachable by controlling the driver variables in *D*. A network which can be controlled from any of its attractors to any of its attractors must have ; when , all attractors reside in disconnected basins in the original STG so . Naturally, if a network is fully controllable by *D* in the control theory sense , .

Given the structural network *G*=(*X*,*E*) for a BN, many different logical functions *f*_{i} can be assigned to each Boolean variable *x*_{i} (see Background). An *ensemble* of BNs is constructed by considering all possible logical functions constrained by the fixed structure *G*^{48}^{,49} (see SM). However, since non-contingent functions (e.g. tautology and contradiction) are not found in most biological models, we divide the full ensemble into *contingent* and *non-contingent* subsets as follows: those BNs which only contain contingent functions and those BNs which contain at least one non-contingent transition function (NC).

Within the set of contingent functions, there are *canalizing* functions which depend only on a subset of their input variables^{16}^{,50}. These functions are ubiquitous in BN models of gene regulation and contribute with mechanisms of functional redundancy and degeneracy^{7}^{,18}. The redundancy of some logical functions means that the *effective* structure of interactions is reduced^{7}^{,14}: some edges of the structural graph *G* play no role in determining the transitions between configurations.

Since control methodologies based on network structure assume that all interactions (edges) in the structural network are relevant for system dynamics, we further subdivide the contingent subset into two disjoint subsets: BNs which contain fully canalizing functions and thus possess a *reduced effective structure* (RES), and those without canalizing functions retaining a *full effective structure* (FES). Naturally, the FES subset is the scenario most coherent with the idea of using structure to predict controllability, since all interactions in the underlying structural graph *G* are dynamically relevant.

We first consider the entire ensemble of BNs with simple structural graphs known as *network motifs*^{51}. These prototype networks have been useful for exploring the relationship between structure and dynamics of complex networks^{52}^{,53}. The motifs considered in our analysis are depicted in Fig. 2.

Directed network structure motifs used in ensemble study: (**A**) Feed-Forward motif, (**B**) Chain motif, (**C**) Loop motif, (**D**) Loop motif with self-interactions, (**E**) Fan motif, (**F**) Co-regulated motif, (**G**) Co-regulating motif, (**H**) BiParallel motif, (**I**) BiFan motif **...**

Consider the Feed-Forward network motif of *N*=3 variables^{54} shown in Fig. 2A. In this case, the full ensemble consists of 64 distinct BNs of which 36 are NC, 8 have RES, and 20 have FES. Figure 1 depicts the logic of one FES network instance for this motif, along with its STG, CSTGs, and CAGs for various driver sets *D*. The control portrait of the full BN ensemble is shown in Fig. 3; control measures and are shown for all possible driver sets of one or two variables.

Using solely this motif’s interaction network, structural control (both the SC and MDS methods) predicts that variable *x*_{1} is capable of fully controlling the network. However, our analysis reveals that this driver variable can fully control only 8 networks from the ensemble (4RES and 4FES), while the other 56BNs (with the same structure) are not fully controlled (Fig. 3). It is noteworthy that even when considering the FES subset — the scenario most coherent with the idea of using structure to predict the controllability of the dynamics — only 4 out of 20BNs are fully controlled by interventions on *x*_{1}. It is clear that even in the case of such a simple motif, structure does not predict the control of dynamics. An extended analysis of the controlled Feed-Forward BN ensemble is provided in the SM.

Let us now consider the *N*=3 variable loop motif with self-interactions (Fig. 2D). The full ensemble of BNs constrained by this motif is much larger than the previous example (every variable has *k*_{i}=2 inputs); it consists of 4096 networks of which 1352 are NC, 1744 have RES, and 1000 have FES. Figure 4A shows the control portrait of this motif’s BN ensemble for a single (*D*{*x*_{i}}) or pair (*D*{*x*_{i},*x*_{j}}) of driver variables. The control portrait of the STG illustrates the difference between the two measures of controllability. While varies greatly, for all BNs. This means that in some BNs, many configurations can be reached simply because the transient dynamics move through many network configurations. Structural control methodologies ignore this natural propensity for control (self-organization). Thus we use the measure to tally only the proportion of transitions that result from control interventions.

The control portraits in Fig. 4 again demonstrate that structure fails to characterize network control. In this case, SC predicts that any single variable is sufficient for full controllability, while MDS requires any two variables to achieve the same. Yet controllability varies greatly for both cases, depending on the particular transition functions of each BN in the ensemble. For 77% of the BN in the ensemble a single variable is not capable of fully controlling dynamics; even two-variable driver sets fail to control 44% of the BNs.

Similar results hold for the mean fraction of reachable attractors shown in Fig. 4B (middle, right). For 36% of the BNs in the ensemble, a single variable is not capable of fully controlling the system between attractors; even two-variable driver sets fail to control 20% of the BNs, regardless of the dynamical subset. Discounting the 1868 networks (Fig. 4B, left) with only one attractor (hence ) further emphasizes the variation in attractor control, increasing the above proportions to 65% and 36% for one and two driver variables, respectively. Therefore, even if we analyze controllability from the point of view of attractor control rather than the stringent criteria of full controllability, single- and two-variable driver sets fail to achieve controllability of all networks in this ensemble.

The control portraits of the other network motifs analyzed are presented in the Supplemental Material. Their analysis supports the same conclusion: predictions made from structure-only methods are only true for a small number of possible BNs. In general, they fail to predict the actual controllability of all the BN dynamics that can occur for a given motif structure.

To better understand the interplay between structure and dynamics in the context of controlling complex systems, we study three BN models from systems biology which are considerably larger than the network motifs of the previous section.

During the early ontogenesis of the fruit fly, the specification of adult cell types is controlled by a hierarchy of a few genes. The Albert and Othmer segment polarity network (SPN) is a BN model^{55} capable of predicting the steady-state patterns experimentally observed in wild-type and mutant embryonic development with significant accuracy. Here, we analyze the single-cell SPN consisting of 17 gene and protein variables (see SM).

Previous analysis has shown that the SPN model is controlled by the upstream value of the Sloppy Pair Protein (SLP) and the extra-cellular signals of the Hedgehog and Wingless proteins from neighboring cells nhh/nHH and nWG^{55}. The control portrait of this model also demonstrates that these three variables (driver set in Fig. 5) are capable of fully controlling the dynamics from any attractor to any other attractor. This is to be expected in segment polarity regulation since it is a highly orchestrated developmental process. The attractor control ability of individual nodes of the SPN in the inset of Fig. 5 further highlights this behavior, only the 3 chemical species mentioned above have a high when controlled alone, while all internal variables have negligible influence.

The SC analysis of the SPN’s structural graph identifies 4 subsets of |*D*|=4 driver variables, indicated in Fig. 5 by enlarged red circles and labeled , , and (details in SM). is a subset of these 4 variable subsets, so naturally they also achieve , but they all include an additional variable which is redundant for this purpose. However, none of these subsets are sufficient for fully controlling the BN as predicted by SC, these driver sets can control dynamics only to a very small proportion of configurations; is the maximum value attained. These 4 driver sets also show considerable variation in , demonstrating that predictions with equivalent support from the point of view of the SC theory, lead to distinct amounts of real controllability. Interestingly, there are 5 driver variable sets of size |*D*|=4 that lead to greater controllability (with a maximum of ) than predicted by SC. Thus, SC fails to even correctly predict the 4-variable driver sets with greatest controllability.

The MDS analysis of the SPN model predicts that |*D*|=7 variables are required to fully control the system dynamics and uncovers 8 equivalent driver variable sets of this size (see SM). Not surprisingly, all of the MDS driver variable sets achieve full attractor control since they contain ; however, none can fully control the network dynamics achieving only a maximum . Thus, the driver sets predicted by both SC and MDS are not sufficient to control dynamics in the control theory sense, and predict more variables than necessary to achieve attractor control.

The eukaryotic cell cycle process of the budding yeast *Saccharomyces cerevisiae* reflects the cyclical gene expression activity that leads to cell division. Here, we use the 12 variable simplified Boolean model of the yeast Cell-Cycle Network (CCN) derived by Li *et al.*^{17}. The SC analysis of the CCN interaction graph identifies only one driver variable to be sufficient for fully controlling the BN’s dynamics. Yet, as demonstrated in Fig. 6A, it only achieves negligible configuration control () and very weak attractor control (A_y0=0.19). Similarly, MDS analysis identifies 8 driver variable sets of size |*D*|=4 ( to ), none of which achieve full control. It is particularly interesting that the driver sets predicted by MDS lead to values of both and that are essentially random, demonstrating once again that predictions with equivalent support from the point of view of the structure-only theories lead to widely different amounts of real controllability. Our analysis finds 3 driver sets of |*D*|=4 variables that achieve full attractor control (highlighted in yellow in Fig. 6A and detailed in SM). Neither SC nor MDS predict those specific driver sets, which ultimately provide the most useful form of control in such systems. Unlike the SPN, there are no “chief controller” variables in this network, as most variables achieve a similar value of when controlled alone (see inset in Fig. 6A).

(**A**) Control of the eukaryotic cell cycle of budding yeast *Saccharomyces cerevisiae* (CCN) for all driver variable subsets of size |*D*|=1, |*D*|=2, |*D*|=3 and |*D*|=4. (inset) The mean fraction **...**

The CCN was designed such that there is a large attractor basin towards a wild-type attractor which is robust to perturbation^{17}^{,44}. However, our analysis illuminates the tradeoff between robustness and flexibility in relation to system controllability. While a large basin of attraction facilitates controlling the system towards the wild-type behavior (high wild-type robustness), it also reduces the ability to control the system to other smaller basins of attraction (mutant phenotypic behavior), reflecting a tradeoff between wild-type robustness and low flexibility for potential evolvability (a property that was not initially designed into the model to begin with). This tradeoff is further elaborated by the CAGs for all single-variable driver sets, shown in Fig. 6B. Some variables have a propensity to control the system towards the wild-type attractor (green node) or allow the system to remain there (e.g. *Cln*3, *Clb*5,6, *Clb*1,2, *Mcm*1/*SFF*, *Cdc*20/14), while only a few can control the system out of this attractor (e.g. *CellSize*, *SBF*, *Cln*1,2). See SM for more details. A third model of biochemical regulation in the floral organ arrangement in the flowering plant *Arabidopsis thaliana* was analyzed, leading to a similar failure to predict actual control (see SM for details).

When fully canalizing functions are present in a BN, not all of the edges in the structural graph contribute to the collective dynamics; there exists a subgraph that fully captures the dynamically relevant interactions (an *effective* structural graph)^{7}^{,14}. Moreover, most Boolean functions are *partially* canalizing^{7}^{,50} whereby in some input conditions a subset of inputs is redundant, but in other conditions it is not. This means that most edges in the underlying structural graph of a random BN are either entirely or partially redundant.

Since structural controllability methods assume that every edge of the underlying structure fully contributes to the dynamics, it is reasonable to suspect that the larger the *mismatch* between the structural graph and the effective structural graph, the more the predictions from SC and MDS will fail. To study this hypothesis, we constructed several ensembles of BNs where there is a perfect match between the structural graph and the effective structure graph.

First consider the ensemble of BNs with the structural graph of the CCN, but with transition functions chosen from the set of two non-canalizing functions that exist for each variable’s in-degree. This constitutes a *Full Effective Connectivity* (FEC) ensemble of BNs whose effective structure perfectly matches the original structural graph of the CCN — there is no canalization in the dynamics of these networks.

Even though both SC and MDS fail to predict controllability correctly for a sample of 50 networks from the FEC ensemble, our analysis reveals that they are more easily controlled by smaller driver sets than the original CCN model. Specifically, and averaged over all driver variable sets is larger for every FEC sample than for the original CCN model (details in SM). Many networks in the FEC ensemble were fully controllable by 2 driver variables and all networks could be fully controlled by 3 driver variables—whereas the original CCN requires 4 variables for full attractor control.

Interestingly, canalization can also be used to improve controllability if selected appropriately. To see how, we compare BN ensembles with no canalization whatsoever to those with only fully canalizing functions for each motif (see SM for details). This uncovers the cases where canalization actually improves BN controllability, even beyond the controllability attained by networks with no canalization. In all such cases, the resulting effective structure reduces the original structural graph to simpler linear chain motifs (Fig. S36 in SM). This way, canalization of the individual variable transition functions is orchestrated to obtain pathways that channel the collective dynamics towards greater control (macro-level canalization^{7}). Because these linear chain effective structures match the assumptions of structure-only methods more accurately, their predictions are correct in such cases. Thus, canalization can enhance the accuracy of structure-only control methodologies if transition functions are appropriately selected to reduce the effective structure to a linear chain. Naturally, when the size of the network increases from simple motifs to realistic networks, BNs with such precise effective structure become extremely rare in the ensembles.

We studied the interplay between structure and dynamics in the control of complex systems using ensembles of BNs and existing models of biochemical regulation. The analysis of the BN ensembles constrained by network motifs demonstrates that structure-only methods fail to properly characterize control; there is a large variation of possible dynamics that can occur for even the simplest network. The situation only gets worse for structure-only methods when we scale up to real models of biochemical regulation. Our analysis demonstrates that structural control predictions can both underestimate or overestimate the number of driver variables in these systems. These approaches also fail to predict which sets of variables best control dynamics as evaluated by: how much of the total configuration space is accessible , how much of the configuration space is accessible beyond the natural system dynamics , and the ability to transition between attractors . Often, arguments made about how easy it is to control network types (e.g. biological vs. social^{24}) hinge on how many driver variables are predicted by structural control theories. Yet, our analysis reveals that much variation in real control occurs for the same structure and number of driver variables.

Our approach also lays the groundwork for understanding which restrictions must be enforced on the transition functions of BNs such that structure may suffice for predicting controllability or at least improve the accuracy of structure-only methods in predicting control. In our experiments with ensembles of network motifs, canalizing transition functions generally rendered structure-only methods less effective at predicting the control of dynamics. Given the generality of motifs as network building blocks, this suggests our results will generalize to larger systems, as already observed in the three larger gene-regulation models considered here. On the other hand, we showed that it is possible to orchestrate canalization such that the effective structure matches the assumptions of structure-only methods, leading to more accurate predictions about control. This effect was identifiable in small networks, where it is easy to find the necessary effective structures, however, such structures are rare in the space of all possible dynamics for larger networks. Nonetheless, in principle, evolution or human design could select for such networks.

Crucially, without more information about variable dynamics, we certainly cannot assume that a given multi-variate dynamical system meets the assumptions of structure-only methods. For instance, the CCN model uses canalization to make controllability harder than predicted by structure-only methods, while the SPN model uses canalization to control dynamics to the wild-type attractor more easily than suggested by the same methods. All this suggests that canalization plays an important, nontrivial role in determining structure-dynamic relationships. Further research can explore this interplay in greater detail. But our current analysis suggests that, without more information about variable dynamics, structure-only methods cannot be accepted as even an approximation of how control occurs in complex systems.

The control measures we introduced here for BNs provide a complementary viewpoint to those developed to study system robustness^{44}^{,56}. Both concepts are based on the response of the system to perturbations. However, robustness focuses on the quantity of perturbations to which the system’s dynamics is invariant, whereas control tracks the perturbations which alter the system’s dynamics. Future research will also explore other characteristics of the controlled state transition graph and controlled attractor graph so that the relationship between robustness and control can be better studied.

Boolean Networks are ideal, parsimonious systems for our study since they are defined by both a clear interaction structure and rich nonlinear dynamics using only binary variables. However, our conclusions are not necessarily limited to this type of network. The control measures used in our study are formulated with respect to a state transition graph, and are therefore applicable to any discrete, deterministic dynamical system. Our conclusions are thus likely to extend to other classes of complex systems. Indeed, several recent papers have also questioned the validity of structure-only arguments for control of other non-linear systems^{42}. These arguments are grounded in the treatment of finite time constants and self-interactions^{30}, the numerical limitations of nonlocal controlled trajectories^{31}, or the role of symmetry in the non-linear dynamics^{57}. Understanding the discrepancy between network structure and control is also important for specific applications where methods which construct a specific controller (i.e. an algorithm that identifies a specific sequence of controlled interventions given a set of constraints) are desired. Structure-only predictions do not aim to predict controllers, rather they focus on the mere identification of driver variables. The identification of controllers is the subject of much research in systems biology and complex systems; in this case, a greater disparity between structure-only predictions and actual control is expected^{42}.

Ultimately, methodologies that can help us predict control in complex networks while avoiding computational complexity should be developed, but they must combine characteristics of both the structural and dynamical properties of the system. Promising methods are already being developed which include both structure and dynamics, such as monotone control systems^{58}, master stability functions^{59}, schema redescription^{7}, and stabilization subgraphs^{60}. Understanding how such simplifications scale-up while providing a reasonable account of how control operates is very important, especially in real-world systems. This can be accomplished via the type of study we undertook here to analyze the effectiveness of structure-only methods in predicting the controllability of complex systems.

**How to cite this article**: Gates, A. J. and Rocha, L. M. Control of complex networks requires both structure and dynamics. *Sci. Rep.*
**6**, 24456; doi: 10.1038/srep24456 (2016).

We thank Randall Beer, Artemy Kolchinsky, Santosh Manicka, Eran Agmon, Ian Wood, and three anonymous reviewers for helpful conversations and feedback. This work was partially supported by a grant from the National Institutes of Health, National Library of Medicine Program, grant 01LM011945-01 “BLR: Evidence-based Drug-Interaction Discovery: *In-Vivo*, *In-Vitro* and Clinical”, a fellowship from NSF IGERT. *The Dynamics of Brain-Body-Environment Systems in Behavior and Cognition* a grant from the Fundação para a Ciencia e a Tecnologia (Portugal), PTDC/EIA-CCO/114108/2009 “Collective Computation and Control in Complex Biochemical Systems”, as well as a grant from the joint program between the Fundação Luso-Americana para o Desenvolvimento (Portugal) and the National Science Foundation (USA), 2012-2014, “Network Mining For Gene Regulation And Biochemical Signaling.” The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. The source code is available upon request.

**Author Contributions** A.J.G. and L.M.R. designed the study, analyzed the results, and wrote the manuscript.

- Kauffman S. A. The Origins of Order: self-organization and selection in evolution (Oxford University Press, New York, 1993).
- Newman M. The structure and function of complex networks. SIAM Rev. 167–256 (2003).
- Huang S. Gene expression profiling, genetic networks, and cellular states: an integrating concept for tumorigenesis and drug discovery. J. Mol. Med. 77, 469–480 (1999). [PubMed]
- Barabási A.-L. & Oltvai Z. N. Network biology: understanding the cell’s functional organization. Nat. Rev. Genet. 5, 101–113 (2004). [PubMed]
- Zhu X., Gerstein M. & Snyder M. Getting connected: analysis and principles of biological networks. Genes Dev. 21, 1010–1024 (2007). [PubMed]
- Assmann S. M. & Albert R. Discrete dynamic modeling with asynchronous update, or how to model complex systems in the absence of quantitative information. In Belostotsky D. A. (ed.) Plant Systems Biology vol. 553 of Methods in Molecular Biology 207–225 (Humana Press, 2009). [PubMed]
- Marques-Pita M. & Rocha L. M. Canalization and control in automata networks: body segmentation in Drosophila melanogaster. Plos One 8, e55946 (2013). [PMC free article] [PubMed]
- Wang R.-S. & Albert R. Elementary signaling modes predict the essentiality of signal transduction network components. BMC systems biology 5 (2011). [PMC free article] [PubMed]
- Strogatz S. H. Exploring complex networks. Nature 410, 268–276 (2001). [PubMed]
- Klemm K. & Bornholdt S. Topology of biological networks and reliability of information processing. Proc. Natl. Acad. Sci. USA 102, 18414–18419 (2005). [PubMed]
- Shmulevich I., Kauffman S. A. & Aldana M. Eukaryotic cells are dynamically ordered or critical but not chaotic. Proc. Natl. Acad. Sci. USA 102, 13439–13444 (2005). [PubMed]
- Nykter M.
*et al.*Gene expression dynamics in the macrophage exhibit criticality. Proc. Natl. Acad. Sci. USA 105, 1897–1900 (2008). [PubMed] - Hossein S., Reichl M. D. & Bassler K. E. Symmetry in critical random boolean network dynamics. Phys. Rev. E 89, 042808 (2014). [PubMed]
- Marques-Pita M., Manicka S., Teuscher C. & Rocha L. M. Effective Connectivity as an Order Parameter in Random Boolean Networks Submitted (2016).
- 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 100, 10734–10739 (2003). [PubMed]
- Kauffman S., Peterson C., Samuelsson B. & Troein C. Genetic networks with canalyzing Boolean rules are always stable. Proc. Natl. Acad. Sci. USA 101, 17102–17107 (2004). [PubMed]
- Li F., Long T., Lu Y., Ouyang Q. & Tang C. The yeast cell-cycle network is robustly designed. Proc. Natl. Acad. Sci. USA 101, 4781–4786 (2004). [PubMed]
- Gershenson C., Kauffman S. A. & Shmulevich I. The role of redundancy in the robustness of random boolean networks. In Artificial Life X (MIT Press, 2006).
- Lin C. Structural controllability. IEEE Trans. Automat. Contr. 19, 201–208 (1974).
- Liu Y.-Y., Slotine J.-J. & Barabási A.-L. Controllability of complex networks. Nature 473, 167–173 (2011). [PubMed]
- Nacher J. C. & Akutsu T. Dominating scale-free networks with variable scaling exponent: heterogeneous networks are not difficult to control. New J. Phys. 14, 073005 (2012).
- Nacher J. C. & Akutsu T. Structural controllability of unidirectional bipartite networks. Scientific Reports 3 (2013). [PMC free article] [PubMed]
- Valente T. Network Interventions. Science 337, 49–53 (2012). [PubMed]
- Egerstedt M. Complex networks: Degrees of control. Nature 473, 158–159 (2011). [PubMed]
- Ruths J. & Ruths D. Control profiles of complex networks. Science 343, 1373–1376 (2014). [PubMed]
- Delpini D.
*et al.*Evolution of controllability in interbank networks. Scientific Reports 3 (2013). [PMC free article] [PubMed] - Österlund T., Bordel S. & Nielsen J. Controllability analysis of transcriptional regulatory networks reveals circular control patterns among transcription factors. Integrative Biology 7, 560–568 (2015). [PubMed]
- Liu Y.-Y., Slotine J.-J. & Barabási A.-L. Observability of complex systems. Proc. Natl. Acad. Sci. USA 110, 2460–2465 (2013). [PubMed]
- Müller F.-J. & Schuppert A. Few inputs can reprogram biological networks. Nature 478, E4–E5 (2011). [PubMed]
- Cowan N. J., Chastain E. J., Vilhena D. A., Freudenberg J. S. & Bergstrom C. T. Nodal Dynamics, Not Degree Distributions, Determine the Structural Controllability of Complex Networks. Plos One 7, e38398 (2012). [PMC free article] [PubMed]
- Sun J. & Motter A. E. Controllability Transition and Nonlocality in Network Control. Phys. Rev. Lett. 110, 208701 (2013). [PubMed]
- Wuchty S. Controllability in protein interaction networks. Proc. Natl. Acad. Sci. USA 111, 7156–7160 (2014). [PubMed]
- Wang B.
*et al.*Diversified control paths: A significant way disease genes perturb the human regulatory network. PLos One 10 (2015). [PMC free article] [PubMed] - Bornholdt S. Boolean network models of cellular regulation: prospects and limitations. J. R. Soc. Interface 5, S85–S94 (2008). [PubMed]
- Akutsu T., Hayashida M., Ching W.-K. & Ng M. K. Control of Boolean networks: hardness results and algorithms for tree structured networks. J. Theor. Biol. 244, 670–679 (2007). [PubMed]
- Langmead C. J. & Jha S. K. Symbolic approaches for finding control strategies in boolean networks. J. Bioinform. Comput. Biol. 7, 323–338 (2009). [PubMed]
- Cheng D. & Qi H. Controllability and observability of Boolean control networks. Automatica 45, 1659–1667 (2009).
- Srihari S., Raman V., Leong H. W. & Ragan M. A. Evolution and Controllability of Cancer Networks: A Boolean Perspective. IEEE Trans. Control Netw. Syst. 11, 83–94 (2013). [PubMed]
- Jia T. & Barabási A.-L. Control Capacity and A Random Sampling Method in Exploring Controllability of Complex Networks. Scientific Reports 3 (2013). [PMC free article] [PubMed]
- Li R., Yang M. & Chu T. Controllability and observability of boolean networks arising from biology. Chaos 25, 023104 (2015). [PubMed]
- Lu W., Tamura T., Song J. & Akutsu T. Computing smallest intervention strategies for multiple metabolic networks in a boolean model. J. Comp. Biol. 22, 85–110 (2015). [PMC free article] [PubMed]
- Motter A. E. Networkcontrology. Chaos 25, 097621 (2015). [PubMed]
- Wuensche A.
Discrete dynamical networks and their attractor basins. In Standish R.
*et al.*(eds) Complex Systems’98 (University of New South Wales, Sydney, Australia, 1998). - Willadsen K. & Wiles J. Robustness and state-space structure of Boolean gene regulatory models. J. Theor. Biol. 249, 749–765 (2007). [PubMed]
- Sontag E. D. Mathematical control theory: deterministic finite dimensional systems. Springer, New York (1998).
- Kauffman S. A. Metabolic stability and epigenesis in randomly constructed genetic nets. J. Theor. Biol. 22, 437–467 (1969). [PubMed]
- Zhang R.
*et al.*Network model of survival signaling in large granular lymphocyte leukemia. Proc. Natl. Acad. Sci. USA 105, 16308–16313 (2008). [PubMed] - Kauffman S. A., Peterson C., Samuelsson B. & Troein C. Random Boolean network models and the yeast transcriptional network. Proc. Natl. Acad. Sci. USA 100, 14796–14799 (2003). [PubMed]
- Ciliberti S., Martin O. C. & Wagner A. Innovation and robustness in complex regulatory gene networks. Proc. Natl. Acad. Sci. USA 104, 13591–13596 (2007). [PubMed]
- Reichhardt C. J. O. & Bassler K. Canalization and symmetry in boolean models for genetic regulatory networks. Physica A 40, 4339–4350 (2007).
- Shen-Orr S. S., Milo R., Mangan S. & Alon U. Network motifs in the transcriptional regulation network of Escherichia coli. Nature Genetics 31, 64–68 (2002). [PubMed]
- Ingram P. J., Stumpf M. P. & Stark J. Network motifs: structure does not determine function. BMC Genomics 7, 108 (2006). [PMC free article] [PubMed]
- Prill R. J., Iglesias P. A. & Levchenko A. Dynamic Properties of Network Motifs Contribute to Biological Network Organization. PLos Biology 3, e343 (2005). [PubMed]
- Mangan S. & Alon U. Structure and function of the feed-forward loop network motif. Proc. Natl. Acad. Sci. USA 100, 11980–11985 (2003). [PubMed]
- Albert R. & Othmer H. G. The topology of the regulatory interactions predicts the expression pattern of the segment polarity genes in Drosophila melanogaster. J. Theor. Biol. 223, 1–18 (2003). [PubMed]
- Chaves M., Sontag E. D. & Albert R. Methods of robustness analysis for boolean models of gene control networks. IEE P. Syst. Biol. 153, 154–167 (2006). [PubMed]
- Whalen A. J., Brennan S. N., Sauer T. D. & Schiff S. J. Observability and controllability of nonlinear networks: The role of symmetry. Phys. Rev. X 5, 011005 (2015).
- Angeli D. & Sontag E. D. Monotone control systems. IEEE Trans. Automat. Contr. 48, 1684–1698 (2003).
- Gutiérrez R., Sendiña-Nadal I., Zanin M., Papo D. & Boccaletti S. Targeting the dynamics of complex networks. Scientific Reports 2 (2012). [PMC free article] [PubMed]
- Zañudo J. G. & Albert R. Cell fate reprogramming by control of intracellular network dynamics. PLos Comput. Biol. 11, e1004193 (2015). [PMC free article] [PubMed]

Articles from Scientific Reports are provided here courtesy of **Nature Publishing Group**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |