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.

The causes of major and rapid transitions observed in biological macroevolution as well as in the evolution of social systems are a subject of much debate. Here we identify the proximate causes of crashes and recoveries that arise dynamically in a model system in which populations of (molecular) species coevolve with their network of chemical interactions. Crashes are events that involve the rapid extinction of many species, and recoveries the assimilation of new ones. These are analyzed and classified in terms of the structural properties of the network. We find that in the absence of large external perturbation, “innovation” is a major cause of large extinctions and the prime cause of recoveries. Another major cause of crashes is the extinction of a “keystone species.” Different classes of causes produce crashes of different characteristic sizes.

Recently there has been a lot of interest in identifying modules at the level of genetic and metabolic networks of organisms, as well as in identifying single genes and reactions that are essential for the organism. A goal of computational and systems biology is to go beyond identification towards an explanation of specific modules and essential genes and reactions in terms of specific structural or evolutionary constraints. In the metabolic networks of Escherichia coli, Saccharomyces cerevisiae and Staphylococcus aureus, we identified metabolites with a low degree of connectivity, particularly those that are produced and/or consumed in just a single reaction. Using flux balance analysis (FBA) we also determined reactions essential for growth in these metabolic networks. We find that most reactions identified as essential in these networks turn out to be those involving the production or consumption of low degree metabolites. Applying graph theoretic methods to these metabolic networks, we identified connected clusters of these low degree metabolites. The genes involved in several operons in E. coli are correctly predicted as those of enzymes catalyzing the reactions of these clusters. Furthermore, we find that larger sized clusters are over-represented in the real network and are analogous to a 'network motif. Using FBA for the above mentioned three organisms we independently identified clusters of reactions whose fluxes are perfectly correlated. We find that the composition of the latter 'functional clusters' is also largely explained in terms of clusters of low degree metabolites in each of these organisms. Our findings mean that most metabolic reactions that are essential can be tagged by one or more low degree metabolites. Those reactions are essential because they are the only ways of producing or consuming their respective tagged metabolites. Furthermore, reactions whose fluxes are strongly correlated can be thought of as 'glued together' by these low degree metabolites. The methods developed here could be used in predicting essential reactions and metabolic modules in other organisms from the list of metabolic reactions.

Experiments on steady-state bacterial cultures have uncovered several quantitative regularities at the system level. These include, first, the exponential growth of cell size with time and the balanced growth of intracellular chemicals between cell birth and division, which are puzzling given the nonlinear and decentralized chemical dynamics in the cell. We model a cell as a set of chemical populations undergoing nonlinear mass action kinetics in a container whose volume is a linear function of the chemical populations. This turns out to be a special class of dynamical systems that generically has attractors in which all populations grow exponentially with time at the same rate. This explains exponential balanced growth of bacterial cells without invoking any regulatory mechanisms and suggests that this could be a robust property of protocells as well. Second, we consider the hypothesis that cells commit themselves to division when a certain internal chemical population reaches a threshold of $N$ molecules. We show that this hypothesis leads to a simple explanation of some of the variability observed across cells in a bacterial culture. In particular, it reproduces the adder property of cell size fluctuations observed recently in E. coli; the observed correlations among interdivision time, birth volume, and added volume in a generation; and the observed scale of the fluctuations (CV $\approx$ 10–30%) when $N$ is between 10 and 100. Third, upon including a suitable regulatory mechanism that optimizes the growth rate of the cell, the model reproduces the observed bacterial growth laws including the dependence of the growth rate and ribosomal protein fraction on the medium. Thus, the models provide a framework for unifying diverse aspects of bacterial growth physiology under one roof. They also suggest new questions for experimental and theoretical enquiry.

Experiments have found that the growth rate and certain other macroscopic properties of bacterial cells in steady-state cultures depend upon the medium in a surprisingly simple manner; these dependencies are referred to as ‘growth laws’. Here we construct a dynamical model of interacting intracellular populations to understand some of the growth laws. The model has only three population variables: an amino acid pool, a pool of enzymes that transport an external nutrient and produce the amino acids, and ribosomes that catalyze their own and the enzymes’ production from the amino acids. We assume that the cell allocates its resources between the enzyme sector and the ribosomal sector to maximize its growth rate. We show that the empirical growth laws follow from this assumption and derive analytic expressions for the phenomenological parameters in terms of the more basic model parameters. Interestingly, the maximization of the growth rate of the cell as a whole implies that the cell allocates resources to the enzyme and ribosomal sectors in inverse proportion to their respective ‘efficiencies’. The work introduces a mathematical scheme in which the cellular growth rate can be explicitly determined and shows that two large parameters, the number of amino acid residues per enzyme and per ribosome, are useful for making approximations.

A model of an evolving network of interacting molecular species is shown to exhibit repeated rounds of crashes in which several species get rapidly depopulated, followed by recoveries. The network inevitably self- organizes into an autocatalytic structure, which consists of an irreducible “core” surrounded by a parasitic “periphery.” Crashes typically occur when the existing autocatalytic set becomes fragile and suffers a “core shift,” defined graph theoretically. The nature of the recovery after a crash, in particular, the time of recovery, depends upon the organizational structure that survives the crash. The largest eigenvalue of the adjacency matrix of the graph is an important signal of network fragility or robustness.