A Toolbox Model of Evolution of Metabolic Pathways on Networks of Arbitrary Topology

In prokaryotic genomes the number of transcriptional regulators is known to be proportional to the square of the total number of protein-coding genes. A toolbox model of evolution was recently proposed to explain this empirical scaling for metabolic enzymes and their regulators. According to its rules, the metabolic network of an organism evolves by horizontal transfer of pathways from other species. These pathways are part of a larger “universal” network formed by the union of all species-specific networks. It remained to be understood, however, how the topological properties of this universal network influence the scaling law of functional content of genomes in the toolbox model. Here we answer this question by first analyzing the scaling properties of the toolbox model on arbitrary tree-like universal networks. We prove that critical branching topology, in which the average number of upstream neighbors of a node is equal to one, is both necessary and sufficient for quadratic scaling. We further generalize the rules of the model to incorporate reactions with multiple substrates/products as well as branched and cyclic metabolic pathways. To achieve its metabolic tasks, the new model employs evolutionary optimized pathways with minimal number of reactions. Numerical simulations of this realistic model on the universal network of all reactions in the KEGG database produced approximately quadratic scaling between the number of regulated pathways and the size of the metabolic network. To quantify the geometrical structure of individual pathways, we investigated the relationship between their number of reactions, byproducts, intermediate, and feedback metabolites. Our results validate and explain the ubiquitous appearance of the quadratic scaling for a broad spectrum of topologies of underlying universal metabolic networks. They also demonstrate why, in spite of “small-world” topology, real-life metabolic networks are characterized by a broad distribution of pathway lengths and sizes of metabolic regulons in regulatory networks.

[1]  Masanori Arita.  The metabolic world of Escherichia coli is not small. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[2]  Oliver Ebenhöh,et al.  A cross species comparison of metabolic network functions. , 2005, Genome informatics. International Conference on Genome Informatics.

[3]  J. Hammersley,et al.  Percolation processes , 1957, Mathematical Proceedings of the Cambridge Philosophical Society.

[4]  Mikko Alava,et al.  Branching Processes , 2009, Encyclopedia of Complexity and Systems Science.

[5]  E. Meléndez-Hevia,et al.  Economy of design in metabolic pathways: further remarks on the game of the pentose phosphate cycle. , 1988, Journal of theoretical biology.

[6]  Oliver Ebenhöh,et al.  Expanding Metabolic Networks: Scopes of Compounds, Robustness, and Evolution , 2005, Journal of Molecular Evolution.

[7]  Paul L. Krapivsky,et al.  Signatures of Arithmetic Simplicity in Metabolic Network Architecture , 2009, PLoS Comput. Biol..

[8]  Arend Hintze,et al.  Evolution of Complex Modular Biological Networks , 2007, PLoS Comput. Biol..

[9]  I. Ispolatov,et al.  Propagation of large concentration changes in reversible protein-binding networks , 2007, Proceedings of the National Academy of Sciences.

[10]  R. Albert,et al.  The large-scale organization of metabolic networks , 2000, Nature.

[11]  A. Wagner.  The yeast protein interaction network evolves rapidly and contains few redundant duplicate genes. , 2001, Molecular biology and evolution.

[12]  Sebastian Bonhoeffer,et al.  The Evolution of Connectivity in Metabolic Networks , 2005, PLoS biology.

[13]  E. M.,et al.  Statistical Mechanics , 2021, Manual for Theoretical Chemistry.

[14]  Oliver Ebenhöh,et al.  An environmental perspective on metabolism. , 2008, Journal of theoretical biology.

[15]  Hiroyuki Ogata,et al.  KEGG: Kyoto Encyclopedia of Genes and Genomes , 1999, Nucleic Acids Res..

[16]  E. Nimwegen.  Scaling Laws in the Functional Content of Genomes , 2003, physics/0307001.

[17]  R. Milo,et al.  Central carbon metabolism as a minimal biochemical walk between precursors for biomass and energy. , 2010, Molecular cell.

[18]  Sergei Maslov,et al.  Toolbox model of evolution of prokaryotic metabolic networks and their regulation , 2009, Proceedings of the National Academy of Sciences.

[19]  D. Fell,et al.  The small world inside large metabolic networks , 2000, Proceedings of the Royal Society of London. Series B: Biological Sciences.