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

Systems approaches to modelling pathways and networks.

Peer reviewedPreprin

Elementary approaches to microbial growth rate maximisation

This thesis, called Elementary approaches to microbial growth rate maximisation, reports on a theoretical search for principles underlying single cell growth, in particular for microbial species that are selected for fast growth rates. First, the optimally growing cell is characterised in terms of its elementary modes. We prove an extremum principle: a cell that maximises a metabolic rate uses few Elementary Flux Modes (EFMs, the minimal pathways that support steady-state metabolism). The number of active EFMs is bounded by the number of growth-limiting constraints. Later, this extremum principle is extended in a theory that explicitly accounts for self-fabrication. For this, we had to define the elementary modes that underlie balanced self-fabrication: minimal self-supporting sets of expressed enzymes that we call Elementary Growth Modes (EGMs). It turns out that many of the results for EFMs can be extended to their more general self-fabrication analogue. Where the above extremum principles tell us that few elementary modes are used by a rate-maximising cell, it does not tell us how the cell can find them. Therefore, we also search for an elementary adaptation method. It turns out that stochastic phenotype switching with growth rate dependent switching rates provides an adaptation mechanism that is often competitive with more conventional regulatory-circuitry based mechanisms. The derived theory is applied in two ways. First, the extremum principles are used to review the mathematical fundaments of all optimisation-based explanations of overflow metabolism. Second, a computational tool is presented that enumerates Elementary Conversion Modes. These elementary modes can be computed for larger networks than EFMs and EGMs, and still provide an overview of the metabolic capabilities of an organism

Exploiting the pathway structure of metabolism to reveal high-order epistasis

<p>Abstract</p> <p>Background</p> <p>Biological robustness results from redundant pathways that achieve an essential objective, e.g. the production of biomass. As a consequence, the biological roles of many genes can only be revealed through multiple knockouts that identify a <it>set </it>of genes as essential for a given function. The identification of such "epistatic" essential relationships between network components is critical for the understanding and eventual manipulation of robust systems-level phenotypes.</p> <p>Results</p> <p>We introduce and apply a network-based approach for genome-scale metabolic knockout design. We apply this method to uncover over 11,000 minimal knockouts for biomass production in an <it>in silico </it>genome-scale model of <it>E. coli</it>. A large majority of these "essential sets" contain 5 or more reactions, and thus represent complex epistatic relationships between components of the <it>E. coli </it>metabolic network.</p> <p>Conclusion</p> <p>The complex minimal biomass knockouts discovered with our approach illuminate robust essential systems-level roles for reactions in the <it>E. coli </it>metabolic network. Unlike previous approaches, our method yields results regarding high-order epistatic relationships and is applicable at the genome-scale.</p

Metabolic Pathway Analysis: from small to genome-scale networks

The need for mathematical modelling of biological processes has grown alongside with the achievements in the experimental field leading to the appearance and development of new fields like systems biology. Systems biology aims at generating new knowledge through modelling and integration of experimental data in order to develop a holistic understanding of organisms. In the first part of my PhD thesis, I compare two different levels of abstraction used for computing metabolic pathways, constraint-based and graph theoretical methods. I show that the current representations of metabolism as a simple graph correspond to wrong mathematical descriptions of metabolic pathways. On the other hand, the use of stoichiometric information and convex analysis as modelling framework like in elementary flux mode analysis, allows to correctly predict metabolic pathways. In the second part of the thesis, I present two of the first methods, based on elementary flux mode analysis, that can compute metabolic pathways in such large metabolic networks: the K-shortest EFMs method and the EFMEvolver method. These methods contribute to an enrichment of the mathematical tools available to model cell biology and more precisely, metabolism. The application of these new methods to biotechnological problems is also explored in this part. In the last part of my thesis, I give an overview of recent achievements in metabolic network reconstruction and constraint-based modelling as well as open issues. Moreover, I discuss possible strategies for integrating experimental data with elementary flux mode analysis. Further improvements in elementary flux mode computation on that direction are put forward

Investigations on the application of complex cell models in the simulation of bioprocesses

[no abstract

Flux cost functions and the choice of metabolic fluxes

Metabolic fluxes in cells are governed by physical, biochemical, physiological, and economic principles. Cells may show "economical" behaviour, trading metabolic performance against the costly side-effects of high enzyme or metabolite concentrations. Some constraint-based flux prediction methods score fluxes by heuristic flux costs as proxies of enzyme investments. However, linear cost functions ignore enzyme kinetics and the tight coupling between fluxes, metabolite levels and enzyme levels. To derive more realistic cost functions, I define an apparent "enzymatic flux cost" as the minimal enzyme cost at which the fluxes can be realised in a given kinetic model, and a "kinetic flux cost", which includes metabolite cost. I discuss the mathematical properties of such flux cost functions, their usage for flux prediction, and their importance for cells' metabolic strategies. The enzymatic flux cost scales linearly with the fluxes and is a concave function on the flux polytope. The costs of two flows are usually not additive, due to an additional "compromise cost". Between flux polytopes, where fluxes change their directions, the enzymatic cost shows a jump. With strictly concave flux cost functions, cells can reduce their enzymatic cost by running different fluxes in different cell compartments or at different moments in time. The enzymactic flux cost can be translated into an approximated cell growth rate, a convex function on the flux polytope. Growth-maximising metabolic states can be predicted by Flux Cost Minimisation (FCM), a variant of FBA based on general flux cost functions. The solutions are flux distributions in corners of the flux polytope, i.e. typically elementary flux modes. Enzymatic flux costs can be linearly or nonlinearly approximated, providing model parameters for linear FBA based on kinetic parameters and extracellular concentrations, and justified by a kinetic model

Elasticity sampling links thermodynamics to metabolic control

Metabolic networks can be turned into kinetic models in a predefined steady state by sampling the reaction elasticities in this state. Elasticities for many reversible rate laws can be computed from the reaction Gibbs free energies, which are determined by the state, and from physically unconstrained saturation values. Starting from a network structure with allosteric regulation and consistent metabolic fluxes and concentrations, one can sample the elasticities, compute the control coefficients, and reconstruct a kinetic model with consistent reversible rate laws. Some of the model variables are manually chosen, fitted to data, or optimised, while the others are computed from them. The resulting model ensemble allows for probabilistic predictions, for instance, about possible dynamic behaviour. By adding more data or tighter constraints, the predictions can be made more precise. Model variants differing in network structure, flux distributions, thermodynamic forces, regulation, or rate laws can be realised by different model ensembles and compared by significance tests. The thermodynamic forces have specific effects on flux control, on the synergisms between enzymes, and on the emergence and propagation of metabolite fluctuations. Large kinetic models could help to simulate global metabolic dynamics and to predict the effects of enzyme inhibition, differential expression, genetic modifications, and their combinations on metabolic fluxes. MATLAB code for elasticity sampling is freely available

Interval and Possibilistic Methods for Constraint-Based Metabolic Models

This thesis is devoted to the study and application of constraint-based metabolic models. The objective was to find simple ways to handle the difficulties that arise in practice due to uncertainty (knowledge is incomplete, there is a lack of measurable variables, and those available are imprecise). With this purpose, tools have been developed to model, analyse, estimate and predict the metabolic behaviour of cells. The document is structured in three parts. First, related literature is revised and summarised. This results in a unified perspective of several methodologies that use constraint-based representations of the cell metabolism. Three outstanding methods are discussed in detail, network-based pathways analysis (NPA), metabolic flux analysis (MFA), and flux balance analysis (FBA). Four types of metabolic pathways are also compared to clarify the subtle differences among them. The second part is devoted to interval methods for constraint-based models. The first contribution is an interval approach to traditional MFA, particularly useful to estimate the metabolic fluxes under data scarcity (FS-MFA). These estimates provide insight on the internal state of cells, which determines the behaviour they exhibit at given conditions. The second contribution is a procedure for monitoring the metabolic fluxes during a cultivation process that uses FS-MFA to handle uncertainty. The third part of the document addresses the use of possibility theory. The main contribution is a possibilistic framework to (a) evaluate model and measurements consistency, and (b) perform flux estimations (Poss-MFA). It combines flexibility on the assumptions and computational efficiency. Poss-MFA is also applied to monitoring fluxes and metabolite concentrations during a cultivation, information of great use for fault-detection and control of industrial processes. Afterwards, the FBA problem is addressed.Llaneras Estrada, F. (2011). Interval and Possibilistic Methods for Constraint-Based Metabolic Models [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/10528Palanci