Analysis of complex metabolic behavior through pathway decomposition

<p>Abstract</p> <p>Background</p> <p>Understanding complex systems through decomposition into simple interacting components is a pervasive paradigm throughout modern science and engineering. For cellular metabolism, complexity can be reduced by decomposition into pathways with particular biochemical functions, and the concept of elementary flux modes provides a systematic way for organizing metabolic networks into such pathways. While decomposition using elementary flux modes has proven to be a powerful tool for understanding and manipulating cellular metabolism, its utility, however, is severely limited since the number of modes in a network increases exponentially with its size.</p> <p>Results</p> <p>Here, we present a new method for decomposition of metabolic flux distributions into elementary flux modes. Our method can easily operate on large, genome-scale networks since it does not require all relevant modes of the metabolic network to be generated. We illustrate the utility of our method for metabolic engineering of <it>Escherichia coli </it>and for understanding the survival of <it>Mycobacterium tuberculosis </it>(MTB) during infection.</p> <p>Conclusions</p> <p>Our method can achieve computational time improvements exceeding 2000-fold and requires only several seconds to generate elementary mode decompositions on genome-scale networks. These improvements arise from not having to generate all relevant elementary modes prior to initiating the decomposition. The decompositions from our method are useful for understanding complex flux distributions and debugging genome-scale models.</p

Organising metabolic networks: cycles in flux distributions

Metabolic networks are among the most widely studied biological systems. The topology and interconnections of metabolic reactions have been well described for many species, but are not sufficient to understand how their activity is regulated in living organisms. The principles directing the dynamic organisation of reaction fluxes remain poorly understood. Cyclic structures are thought to play a central role in the homeostasis of biological systems and in their resilience to a changing environment. In this work, we investigate the role of fluxes of matter cycling in metabolic networks. First, we introduce a methodology for the computation of cyclic and acyclic fluxes in metabolic networks, adapted from an algorithm initially developed to study cyclic fluxes in trophic networks. Subsequently, we apply this methodology to the analysis of three metabolic systems, including the central metabolism of wild type and a deletion mutant of Escherichia coli, erythrocyte metabolism and the central metabolism of the bacterium Methylobacterium extorquens. The role of cycles in driving and maintaining the performance of metabolic functions upon perturbations is unveiled through these examples. This methodology may be used to further investigate the role of cycles in living organisms, their pro-activity and organisational invariance, leading to a better understanding of biological entailment and information processing

Enzyme economy in metabolic networks

Metabolic systems are governed by a compromise between metabolic benefit and enzyme cost. This hypothesis and its consequences can be studied by kinetic models in which enzyme profiles are chosen by optimality principles. In enzyme-optimal states, active enzymes must provide benefits: a higher enzyme level must provide a metabolic benefit to justify the additional enzyme cost. This entails general relations between metabolic fluxes, reaction elasticities, and enzyme costs, the laws of metabolic economics. The laws can be formulated using economic potentials and loads, state variables that quantify how metabolites, reactions, and enzymes affect the metabolic performance in a steady state. Economic balance equations link them to fluxes, reaction elasticities, and enzyme levels locally in the network. Economically feasible fluxes must be free of futile cycles and must lead from lower to higher economic potentials, just like thermodynamics makes them lead from higher to lower chemical potentials. Metabolic economics provides algebraic conditions for economical fluxes, which are independent of the underlying kinetic models. It justifies and extends the principle of minimal fluxes and shows how to construct kinetic models in enzyme-optimal states, where all enzymes have a positive influence on the metabolic performance

Global organization of metabolic fluxes in the bacterium, Escherichia coli

Cellular metabolism, the integrated interconversion of thousands of metabolic substrates through enzyme-catalyzed biochemical reactions, is the most investigated complex intercellular web of molecular interactions. While the topological organization of individual reactions into metabolic networks is increasingly well understood, the principles governing their global functional utilization under different growth conditions pose many open questions. We implement a flux balance analysis of the E. coli MG1655 metabolism, finding that the network utilization is highly uneven: while most metabolic reactions have small fluxes, the metabolism's activity is dominated by several reactions with very high fluxes. E. coli responds to changes in growth conditions by reorganizing the rates of selected fluxes predominantly within this high flux backbone. The identified behavior likely represents a universal feature of metabolic activity in all cells, with potential implications to metabolic engineering.Comment: 15 pages 4 figure

How enzyme economy shapes metabolic fluxes

Metabolic fluxes are governed by physical and economic principles. Stationarity constrains them to a subspace in flux space and thermodynamics makes them lead from higher to lower chemical potentials. At the same time, fluxes in cells represent a compromise between metabolic performance and enzyme cost. To capture this, some flux prediction methods penalise larger fluxes by heuristic cost terms. Economic flux analysis, in contrast, postulates a balance between enzyme costs and metabolic benefits as a necessary condition for fluxes to be realised by kinetic models with optimal enzyme levels. The constraints are formulated using economic potentials, state variables that capture the enzyme labour embodied in metabolites. Generally, fluxes must lead from lower to higher economic potentials. This principle, which resembles thermodynamic constraints, can complement stationarity and thermodynamic constraints in flux analysis. Futile modes, which would be incompatible with economic potentials, are defined algebraically and can be systematically removed from flux distributions. Enzymes that participate in potential futile modes are likely targets of regulation. Economic flux analysis can predict high-yield and low-yield strategies, and captures preemptive expression, multi-objective optimisation, and flux distributions across several cells living in symbiosis. Inspired by labour value theories in economics, it justifies and extends the principle of minimal fluxes and provides an intuitive framework to model the complex interplay of fluxes, metabolic control, and enzyme costs in cells

Rigidity and flexibility of biological networks

The network approach became a widely used tool to understand the behaviour of complex systems in the last decade. We start from a short description of structural rigidity theory. A detailed account on the combinatorial rigidity analysis of protein structures, as well as local flexibility measures of proteins and their applications in explaining allostery and thermostability is given. We also briefly discuss the network aspects of cytoskeletal tensegrity. Finally, we show the importance of the balance between functional flexibility and rigidity in protein-protein interaction, metabolic, gene regulatory and neuronal networks. Our summary raises the possibility that the concepts of flexibility and rigidity can be generalized to all networks.Comment: 21 pages, 4 figures, 1 tabl

Using in silico models to simulate dual perturbation experiments: procedure development and interpretation of outcomes.

BackgroundA growing number of realistic in silico models of metabolic functions are being formulated and can serve as 'dry lab' platforms to prototype and simulate experiments before they are performed. For example, dual perturbation experiments that vary both genetic and environmental parameters can readily be simulated in silico. Genetic and environmental perturbations were applied to a cell-scale model of the human erythrocyte and subsequently investigated.ResultsThe resulting steady state fluxes and concentrations, as well as dynamic responses to the perturbations were analyzed, yielding two important conclusions: 1) that transporters are informative about the internal states (fluxes and concentrations) of a cell and, 2) that genetic variations can disrupt the natural sequence of dynamic interactions between network components. The former arises from adjustments in energy and redox states, while the latter is a result of shifting time scales in aggregate pool formation of metabolites. These two concepts are illustrated for glucose-6 phosphate dehydrogenase (G6PD) and pyruvate kinase (PK) in the human red blood cell.ConclusionDual perturbation experiments in silico are much more informative for the characterization of functional states than single perturbations. Predictions from an experimentally validated cellular model of metabolism indicate that the measurement of cofactor precursor transport rates can inform the internal state of the cell when the external demands are altered or a causal genetic variation is introduced. Finally, genetic mutations that alter the clinical phenotype may also disrupt the 'natural' time scale hierarchy of interactions in the network