The regulatory network of E. coli metabolism as a boolean dynamical system exhibits both homeostasis and flexibility of response

Abstract

Elucidating the architecture and dynamics of large scale genetic regulatory networks of cells is an important goal in systems biology. We study the system level dynamical properties of the genetic network of Escherichia coli that regulates its metabolism, and show how its design leads to biologically useful cellular properties. Our study uses the database (Covert et al., Nature 2004) containing 583 genes and 96 external metabolites which describes not only the network connections but also the boolean rule at each gene node that controls the switching on or off of the gene as a function of its inputs. We have studied how the attractors of the boolean dynamical system constructed from this database depend on the initial condition of the genes and on various environmental conditions corresponding to buffered minimal media. We find that the system exhibits homeostasis in that its attractors, that turn out to be fixed points or low period cycles, are highly insensitive to initial conditions or perturbations of gene configurations for any given fixed environment. At the same time the attractors show a wide variation when external media are varied implying that the system mounts a highly flexible response to changed environmental conditions. The regulatory dynamics acts to enhance the cellular growth rate under changed media. Our study shows that the reconstructed genetic network regulating metabolism in {\it E. coli} is hierarchical, modular, and largely acyclic, with environmental variables controlling the root of the hierarchy. This architecture makes the cell highly robust to perturbations of gene configurations as well as highly responsive to environmental changes. The twin properties of homeostasis and response flexibility are achieved by this dynamical system even though it is not close to the edge of chaos.