Abstract

Abstract: High-level, mathematically precise descriptions of the global organisation of complex metabolic networks are necessary for understanding the global structure of metabolic networks, the interpretation and integration of large amounts of biologic data (sequences, various -omics) and ultimately for rational design of therapies for disease processes. Metabolic networks are highly organised to execute their function efficiently while tolerating wide variation in their environment. These networks are constrained by physical requirements (e.g. conservation of energy, redox and small moieties) but are also remarkably robust and evolvable. The authors use well-known features of the stoichiometry of bacterial metabolic networks to demonstrate how network architecture facilitates such capabilities, and to develop a minimal abstract metabolism which incorporates the known features of the stoichiometry and respects the constraints on enzymes and reactions. This model shows that the essential functionality and constraints drive the tradeoffs between robustness and fragility, as well as the large-scale structure and organisation of the whole network, particularly high variability. The authors emphasise how domainspecific constraints and tradeoffs imposed by the environment are important factors in shaping stoichiometry. Importantly, the consequence of these highly organised tradeoffs and tolerances is an architecture that has a highly structured modularity that is self-dissimilar and scale-rich. Introduction Metabolic networks, which have been extensively studied for decades, are emblematic of how evolution has sculpted biologic systems for optimal function. In addition to unambiguous functional descriptions of core metabolism, this conserved network has been recently described in detail in terms of its stoichiometry (mass and energy balance). A higher level, mathematically defined description of the global organisation of complex metabolic networks is critical for a deep understanding of metabolism, from the interpretation of huge amounts of biologic data (sequences, various -omics) to design of therapies for disease processes. The stakes are high for obtaining the big picture right: biologic data plugged into a distorted model or interpreted in the context of a flawed universal law propagates misinterpretations. In flux analyses [1], stoichiometry is considered as a constraint, and fluxes are optimised to satisfy a global objective, typically growth. Previous studies, however, have not directly addressed whether the stoichiometry itself is highly optimal or organised in any sense and contributes to the origins and purpose of complexity in biological networks. Yet biochemistry textbooks describe metabolism as having evolved to be 'highly integrated' with the appearance of a 'coherent design' [2]. Here we explore both important 'design' (with no implication of a 'designer') features of metabolism and the sense in which stoichiometry itself has highly organised and optimised tolerances and tradeoffs (HOT) Basic features of metabolic networks Metabolism is essentially a linked series of chemical reactions, which function to synthesise building blocks for usable cellular components and to extract energy and reducing power from the cellular environment, in the context of total organism homeostasis. Constraints on the network are imposed by highly unpredictable intracellular and extracellular environments as well as the details of enzyme molecular structure, the cost of making enzymes and the conservation of atoms, energy and small moieties. The simplest model of metabolic networks is a stoichiometry matrix (s-matrix for short) of chemical reactions with the metabolites in rows and reactions in columns and is defined unambiguously except for permutations of rows # IEE, 2005