Emergence and Growth of Complex Networks in Adaptive Systems

We consider the population dynamics of a set of species whose network of catalytic interactions is described by a directed graph. The relationship between the attractors of this dynamics and the underlying graph theoretic structures like cycles and autocatalytic sets is discussed. It is shown that when the population dynamics is suitably coupled to a slow dynamics of the graph itself, the network evolves towards increasing complexity driven by autocatalytic sets. Some quantitative measures of network complexity are described.Comment: 10 pages (including figures), 3 Postscript figure

Autocatalytic Sets and the Growth of Complexity in an Evolutionary Model

A model of $s$ interacting species is considered with two types of dynamical variables. The fast variables are the populations of the species and slow variables the links of a directed graph that defines the catalytic interactions among them. The graph evolves via mutations of the least fit species. Starting from a sparse random graph, we find that an autocatalytic set (ACS) inevitably appears and triggers a cascade of exponentially increasing connectivity until it spans the whole graph. The connectivity subsequently saturates in a statistical steady state. The time scales for the appearance of an ACS in the graph and its growth have a power law dependence on $s$ and the catalytic probability. At the end of the growth period the network is highly non-random, being localized on an exponentially small region of graph space for large $s$.Comment: 13 pages REVTEX (including figures), 4 Postscript figure

Network Models of Phage-Bacteria Coevolution

Bacteria and their bacteriophages are the most abundant, widespread and diverse groups of biological entities on the planet. In an attempt to understand how the interactions between bacteria, virulent phages and temperate phages might affect the diversity of these groups, we developed a novel stochastic network model for examining the co-evolution of these ecologies. In our approach, nodes represent whole species or strains of bacteria or phages, rather than individuals, with "speciation" and extinction modelled by duplication and removal of nodes. Phage-bacteria links represent host-parasite relationships and temperate-virulent phage links denote prophage-encoded resistance. The effect of horizontal transfer of genetic information between strains was also included in the dynamical rules. The observed networks evolved in a highly dynamic fashion but the ecosystems were prone to collapse (one or more entire groups going extinct). Diversity could be stably maintained in the model only if the probability of speciation was independent of the diversity. Such an effect could be achieved in real ecosystems if the speciation rate is primarily set by the availability of ecological niches.Comment: 8 pages, 6 figure

Symbolic dynamics of biological feedback networks

We formulate general rules for a coarse-graining of the dynamics, which we term `symbolic dynamics', of feedback networks with monotone interactions, such as most biological modules. Networks which are more complex than simple cyclic structures can exhibit multiple different symbolic dynamics. Nevertheless, we show several examples where the symbolic dynamics is dominated by a single pattern that is very robust to changes in parameters and is consistent with the dynamics being dictated by a single feedback loop. Our analysis provides a method for extracting these dominant loops from short time series, even if they only show transient trajectories.Comment: 4 pages, 4 figure

Structure and function of negative feedback loops at the interface of genetic and metabolic networks

The molecular network in an organism consists of transcription/translation regulation, protein-protein interactions/modifications and a metabolic network, together forming a system that allows the cell to respond sensibly to the multiple signal molecules that exist in its environment. A key part of this overall system of molecular regulation is therefore the interface between the genetic and the metabolic network. A motif that occurs very often at this interface is a negative feedback loop used to regulate the level of the signal molecules. In this work we use mathematical models to investigate the steady state and dynamical behaviour of different negative feedback loops. We show, in particular, that feedback loops where the signal molecule does not cause the dissociation of the transcription factor from the DNA respond faster than loops where the molecule acts by sequestering transcription factors off the DNA. We use three examples, the bet, mer and lac systems in E. coli, to illustrate the behaviour of such feedback loops.Comment: 8 pages, 4 figure

Flux-based classification of reactions reveals a functional bow-tie organization of complex metabolic networks

Unraveling the structure of complex biological networks and relating it to their functional role is an important task in systems biology. Here we attempt to characterize the functional organization of the large-scale metabolic networks of three microorganisms. We apply flux balance analysis to study the optimal growth states of these organisms in different environments. By investigating the differential usage of reactions across flux patterns for different environments, we observe a striking bimodal distribution in the activity of reactions. Motivated by this, we propose a simple algorithm to decompose the metabolic network into three sub-networks. It turns out that our reaction classifier which is blind to the biochemical role of pathways leads to three functionally relevant sub-networks that correspond to input, output and intermediate parts of the metabolic network with distinct structural characteristics. Our decomposition method unveils a functional bow-tie organization of metabolic networks that is different from the bow-tie structure determined by graph-theoretic methods that do not incorporate functionality.Comment: 11 pages, 6 figures, 1 tabl

Low Degree Metabolites Explain Essential Reactions and Enhance Modularity in Biological Networks

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 E. coli, S. cerevisiae and S. aureus, we identified metabolites with a low degree of connectivity, particularly those that are produced and/or consumed in just a single reaction. Using 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. 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.Comment: 12 pages main text with 2 figures and 2 tables. 16 pages of Supplementary material. Revised version has title changed and contains study of 3 organisms instead of 1 earlie