The paradigmatic example of a complex system is the web of biochemical interactions that make up life. We still know very little about the organization of life as a dynamical, interacting network of genes, proteins and biochemical reactions. How do biochemical networks—containing many regulatory, signaling, and metabolic processes—achieve reliability and robustness? Cells function reliably despite noisy dynamic environments, which is all the more impressive given that **control strategies** implemented by intra and inter-cellular processes cannot rely on a centralized, global view of the relevant networks. Are the resulting complex dynamics made up of relatively autonomous **modules**? If so, what is their functional role and how can they be identified? How robust is the collective computation performed by intra-cellular networks to mutations, delays and stochastic noise? To address these questions, we are focused on developing both novel **methodologies** and **informatics tools** to study **control** and **collective computation** in automata networks used to model gene regulation and biochemical signaling.

The study of network structure has uncovered organizational principles in complex systems. However, there is also a need to understand how to control them; for example, to revert a diseased cell to a healthy state, or a mature cell to a pluripotent state. While a recent interest in linear control has lead much recent research to infer control from the structure of network interactions, we have shown that both structural controllability and minimum dominating sets theory fail to properly characterize controllability in systems biology models of biochemical regulation, and even in small network motifs. Indeed, structure-only methods both undershoot and overshoot the number and which sets of variables actually control these models, highlighting the importance of the system dynamics in determining control [Gates and Rocha, 2014][Gates and Rocha, 2016]. We have also shown that the logic of automata transition functions, namely how **canalizing** they are, plays a role in the extent to which structure predicts dynamics (see more on canalization below).

Considering the possible dynamics that can unfold on specific network structure is a focus of research in our group. In addition to control, we study prediction [Kolchinsky and Rocha, 2011], modularity [Kolchinsky, Gates and Rocha, 2015] and multi-scale integration in the dynamics of complex network, such as brain networks [Kolchinsky et al 2014].

*Schema redescription with two symbols* is a method to eliminate redundancy in the transition tables of Boolean automata. One symbol is used to capture redundancy of individual input variables, and another to capture permutability in sets of input variables: fully characterizing the **canalization** present in Boolean functions [Marques-Pita and Rocha, 2011] [Marques-Pita and Rocha, 2013] (see figure). In our formulation, canalization becomes synonymous with redundancy present in the logic of automata. This results in straightforward measures to quantify canalization in an automaton (micro-level), which is in turn integrated into a highly scalable way to characterize the collective dynamics of large-scale automata networks (macro-level). This way, our approach provides a method to link micro- to macro-level dynamics---a crux of complexity. Our methodology is applicable to any complex network that can be modelled using automata, but we focus on biochemical regulation and signalling, towards a better understanding of the (decentralized) control that orchestrates cellular activity---with the ultimate goal of explaining how do cells and tissues "compute".

By removing redundancy (canalization) from discrete models of biochemical regulation, we can extract the effective structure that **controls** their dynamics, revealing their **dynamical modularity** (modules in the dynamics rather than in the structure of networks) and **robustness** [Marques-Pita and Rocha, 2013]. In particular, we can extract the minimal conditions (as schemata or motifs) and critical nodes that control convergence to attractors---associated with phenotypic behavior in these models. The approach is scalable because it only needs to compute the redundancy of the transition functions of each node in the network, rather than the entire dynamical landscape of the multivariate dynamical system. This has lead us, for instance, to obtain a better understanding of a well-known 60-variable model of the intra- and inter cellular genetic regulation of body segmentation in Drosophila Melanogaster. We were able to measure more accurately the size of its wild-type attractor basin (larger than previously thought), to the identify novel minimal conditions and critical nodes that control wild-type behaviour, and estimate its resilience to stochastic interventions.

We have also used our method to show that despite having very different collective behavior, Cellular Automata (CA) rules can be very similar at the local interaction level [Marques-Pita and Rocha, 2011]—leading us to question the tendency in complexity research to pay much more attention to emergent patterns than to local interactions. Additionally, schema redescription allows us to obtain more amenable search spaces of CA rules for the Density Classification Task—obtaining some of the best known rules for this task. [Marques-Pita and Rocha, 2008, Marques-Pita, Mitchell, and Rocha, 2008].

We have been interested on the problem of how information, symbols, representations and the like can arise from a purely dynamical system of many components. This is a topic of particular interest in Cognitive Science, where the notions of *representation and symbol* often divide the field into opposing camps. Often, in the area of Embodied Cognition the idea of self-organization in dynamical systems leads many researchers to reject representational or semiotic elements in their models of cognition. This attitude seems not only excessive, but indeed absurd as it ignores the informational processes so important for biological organisms. Therefore, we have been working both on a re-formulation of the concept of representation for *embodied cognition*, as well as on simulations of dynamical systems (using Celular Automata) where one can study the origin of representations.

The Evolving Cellular Automata experiments of Crutchfield, Mitchell et al, in the late 1990’s were very exciting, as the ability of evolved cellular automata to solve non-trivial computation tasks seemed to provide clues about the origin of representations and information from dynamical systems [Mitchell, 1998] [Rocha ,1998b]. We conducted additional experiments which extended the density classification task with more difficult logical tasks [Rocha, 2000; Rocha, 2004]. Later, we proposed a re-formulation of the concept of representation in cognitive science and artificial life which is based on this work, but argues that the type of emergent computations observed in these experiments do not produce representations quite as rich as those as observed in biology and cognition [Rocha and Hordijk ,2005]. These experiments allow us to think about how to evolve symbols from artificial matter in computational environments. The figure above, depicts a space-time diagram and particle model of a CA rule evolved to solve the AND task. Some additional Figures and experiment details of CA rules for logical tasks in our experiments are also available.

- Fundacao para a Ciencia e Tecnologia, Portugal. PTDC/MEC-AND/30221/2017. Project title: “
*The sperm cell core genetic program: combined clinical and research approach to the diagnosis of male infertility*”. 2018-2021 - National Science Foundation, Research Traineeship Program, NSF1735095: Interdisciplinary Training in Complex Networks and Systems, 2017-2022
- Fulbright U.S. Scholar grant, J. William Fulbright Foreign Scholarship Board (FFSB). 2016-2017
- Fundacao para a Ciencia e Tecnologia, Portugal. PTDC/EIA-CCO/114108/2009. Project title: “
*Collective Computation and Control in Complex Biochemical Systems*” - Fundação Luso-Americana para o Desenvolvimento (Portugal) and National Science Foundation (USA), 2012-2014. Project title: “Network Mining For Gene Regulation And Biochemical Signaling.” (171/11)

Luis Rocha

Melanie Mitchell

Santiago Schnell

Manuel Marques Pita

Santosh Manicka

Alex Gates

Artemy Kolchinsky

Thomas Parmer

- R.B. Correia, A.J. Gates, X. Wang, L.M. Rocha [2018]. CANA: A python package for quantifying control and canalization in Boolean Networks.
*Frontiers in Physiology*.**9**: 1046. DOI: 10.3389/fphys.2018.01046 - Correia, R.B, Gates, A., Manicka, S., M. Marques-Pita, X. Wang, and L.M. Rocha [2017]. "The effective structure of complex networks: Canalization in the dynamics of complex networks drives dynamics, criticality and control".
*Complex Networks 2017. The 6th International Workshop on Complex Networks & Their Applications*. Nov. 29 - Dec. 01, 2017, Lyon, pp. 354-355 - A. Gates and L.M. Rocha. [2016] "Control of complex networks requires both structure and dynamics."
*Scientific Reports***6**, 24456. doi: 10.1038/srep24456. - A. Kolchinsky, A. Gates and L.M. Rocha. [2015] "Modularity and the spread of perturbations in complex dynamical systems."
*Phys. Rev. E Rapid Communications*.**92**, 060801(R). - A. Kolchinsky, M. P. Van Den Heuvel, A. Griffa, P. Hagmann, L.M. Rocha, O. Sporns, J. Goni [2014]. "Multi-scale Integration and Predictability in Resting State Brain Activity".
*Frontiers in Neuroinformatics*, 8:66. doi: 10.3389/fninf.2014.00066. - A. Gates and L.M. Rocha [2014]. “Structure and dynamics affect the controllability of complex systems: a Preliminary Study”.
*Artificial Life 14: Proceedings of the Fourteenth International Conference on the Synthesis and Simulation of Living Systems*: 429-430, MIT Press. - M. Marques-Pita and L.M. Rocha [2013]. “Canalization and control in automata networks: body segmentation in Drosophila Melanogaster”.
*PLoS ONE*,**8**(3): e55946. doi:10.1371/journal.pone.0055946. - A. Kolchinsky, and L.M. Rocha [2011].”Prediction and Modularity in Dynamical Systems“.In:
*Advances in Artificial Life, Proceedings of the Eleventh European Conference on the Synthesis and Simulation of Living Systems (ECAL 2011). August 8 – 12, 2011, Paris, France,*. MIT Press, pp. 423-430. - M. Mourao [2011].
*Reverse engineering the mechanisms and dynamical behavior of complex biochemical pathways*. PhD Dissertation, Indiana University - M. Marques-Pita and L.M. Rocha [2011]. “Schema Redescription in Cellular Automata: Revisiting Emergence in Complex Systems“. In:
*The 2011 IEEE Symposium on Artificial Life, at the IEEE Symposium Series on Computational Intelligence 2011. April 11 – 15, 201, Paris, France,*. IEEE Press, pp: 233-240. - M. Marques-Pita, M. Mitchell, and L.M. Rocha [2008]. “The Role of Conceptual Structure in Learning Cellular Automata to Perform Collective Computation“. In:
*Unconventional Computation: 7th International Conference (UC 2008)*.*Lecture Notes in Computer Science*. Springer-Verlag,**5204**: 146-163. - Rocha, Luis M. and W. Hordijk [2005]. "Material Representations: From the Genetic Code to the Evolution of Cellular Automata".
*Artificial Life*.**11**(1-2), pp. 189 - 214 - Rocha, Luis M. [2000]. "Syntactic autonomy, cellular automata, and RNA editing: or why self-organization needs symbols to evolve and how it might evolve them". In:
*Closure: Emergent Organizations and Their Dynamics*. Chandler J.L.R. and G, Van de Vijver (Eds.)*Annals of the New York Academy of Sciences*. Vol. 901, pp 207-223.