Enzyme allocation problems in kinetic metabolic networks: optimal solutions are elementary flux modes

Including glutamine in a resource allocation model of energy metabolism in cancer and yeast cells

Energy metabolism is crucial for all living cells, especially during fast growth or stress scenarios. Many cancer and activated immune cells (Warburg effect) or yeasts (Crabtree effect) mostly rely on aerobic glucose fermentation leading to lactate or ethanol, respectively, to generate ATP. In recent years, several mathematical models have been proposed to explain the Warburg effect on theoretical grounds. Besides glucose, glutamine is a very important substrate for eukaryotic cells-not only for biosynthesis, but also for energy metabolism. Here, we present a minimal constraint-based stoichiometric model for explaining both the classical Warburg effect and the experimentally observed respirofermentation of glutamine (WarburQ effect). We consider glucose and glutamine respiration as well as the respective fermentation pathways. Our resource allocation model calculates the ATP production rate, taking into account enzyme masses and, therefore, pathway costs. While our calculation predicts glucose fermentation to be a superior energy-generating pathway in human cells, different enzyme characteristics in yeasts reduce this advantage, in some cases to such an extent that glucose respiration is preferred. The latter is observed for the fungal pathogen Candida albicans, which is a known Crabtree-negative yeast. Further, optimization results show that glutamine is a valuable energy source and important substrate under glucose limitation, in addition to its role as a carbon and nitrogen source of biomass in eukaryotic cells. In conclusion, our model provides insights that glutamine is an underestimated fuel for eukaryotic cells during fast growth and infection scenarios and explains well the observed parallel respirofermentation of glucose and glutamine in several cell types.

Understanding and computational design of genetic circuits of metabolic networks

The expression of metabolic proteins is controlled by genetic circuits, matching metabolic demands and changing environmental conditions. Ideally, this regulation brings about a competitive level of metabolic fitness. Understanding how cells can achieve a robust (close-to-optimal) functioning of metabolism by appropriate control of gene expression aids synthetic biology by providing design criteria of synthetic circuits for biotechnological purposes. It also extends our understanding of the designs of genetic circuitry found in nature such as metabolite control of transcription factor activity, promoter architectures and transcription factor dependencies, and operon composition (in bacteria). Here, we review, explain and illustrate an approach that allows for the inference and design of genetic circuitry that steers metabolic networks to achieve a maximal flux per unit invested protein across dynamic conditions. We discuss how this approach and its understanding can be used to rationalize Escherichia coli's strategy to regulate the expression of its ribosomes and infer the design of circuitry controlling gene expression of amino-acid biosynthesis enzymes. The inferred regulation indeed resembles E. coli's circuits, suggesting that these have evolved to maximize amino-acid production fluxes per unit invested protein. We end by an outlook of the use of this approach in metabolic engineering applications.

Whole-cell metabolic control analysis

Since its conception some fifty years ago, metabolic control analysis (MCA) aims to understand how cells control their metabolism by adjusting the activity of their enzymes. Here we extend its scope to a whole-cell context. We consider metabolism in the evolutionary context of growth-rate maximisation by optimisation of protein concentrations. This framework allows for the prediction of flux control coefficients from proteomics data or stoichiometric modelling. Since genes compete for finite biosynthetic resources, we treat all protein concentrations as interdependent. We show that elementary flux modes (EFMs) emerge naturally as the optimal metabolic networks in the whole-cell context and we derive their control properties. In the evolutionary optimum, the number of expressed EFMs is determined by the number of protein-concentration constraints that limit growth rate. We use published glucose-limited chemostat data of S. cerevisiae to illustrate that it uses only two EFMs prior to the onset of fermentation and that it uses four EFMs during fermentation. We discuss published enzyme-titration data to show that S. cerevisiae and E. coli indeed can express proteins at growth-rate maximising concentrations. Accordingly, we extend MCA to elementary flux modes operating at an optimal state. We find that the expression of growth-unassociated proteins changes results from classical metabolic control analysis. Finally, we show how flux control coefficients can be estimated from proteomics and ribosome-profiling data. We analyse published proteomics data of E. coli to provide a whole-cell perspective of the control of metabolic enzymes on growth rate. We hope that this paper stimulates a renewed interest in metabolic control analysis, so that it can serve again the purpose it once had: to identify general principles that emerge from the biochemistry of the cell and are conserved across biological species.

Mathematical properties of optimal fluxes in cellular reaction networks at balanced growth

The physiology of biological cells evolved under physical and chemical constraints, such as mass conservation across the network of biochemical reactions, nonlinear reaction kinetics, and limits on cell density. For unicellular organisms, the fitness that governs this evolution is mainly determined by the balanced cellular growth rate. We previously introduced growth balance analysis (GBA) as a general framework to model and analyze such nonlinear systems, revealing important analytical properties of optimal balanced growth states. It has been shown that at optimality, only a minimal subset of reactions can have nonzero flux. However, no general principles have been established to determine if a specific reaction is active at optimality. Here, we extend the GBA framework to study the optimality of each biochemical reaction, and we identify the mathematical conditions determining whether a reaction is active or not at optimal growth in a given environment. We reformulate the mathematical problem in terms of a minimal number of dimensionless variables and use the Karush-Kuhn-Tucker (KKT) conditions to identify fundamental principles of optimal resource allocation in GBA models of any size and complexity. Our approach helps to identify from first principles the economic values of biochemical reactions, expressed as marginal changes in cellular growth rate; these economic values can be related to the costs and benefits of proteome allocation into the reactions' catalysts. Our formulation also generalizes the concepts of Metabolic Control Analysis to models of growing cells. We show how the extended GBA framework unifies and extends previous approaches of cellular modeling and analysis, putting forward a program to analyze cellular growth through the stationarity conditions of a Lagrangian function. GBA thereby provides a general theoretical toolbox for the study of fundamental mathematical properties of balanced cellular growth.

Shaping bacterial gene expression by physiological and proteome allocation constraints

Networks of molecular regulators are often the primary objects of focus in the study of gene regulation, with the machinery of protein synthesis tacitly relegated to the background. Shifting focus to the constraints imposed by the allocation of protein synthesis flux reveals surprising ways in which the actions of molecular regulators are shaped by physiological demands. Using carbon catabolite repression as a case study, we describe how physiological constraints are sensed through metabolic fluxes and how flux-controlled regulation gives rise to simple empirical relations between protein levels and the rate of cell growth.

Interrogating the effect of enzyme kinetics on metabolism using differentiable constraint-based models

Metabolic models are typically characterized by a large number of parameters. Traditionally, metabolic control analysis is applied to differential equation-based models to investigate the sensitivity of predictions to parameters. A corresponding theory for constraint-based models is lacking, due to their formulation as optimization problems. Here, we show that optimal solutions of optimization problems can be efficiently differentiated using constrained optimization duality and implicit differentiation. We use this to calculate the sensitivities of predicted reaction fluxes and enzyme concentrations to turnover numbers in an enzyme-constrained metabolic model of Escherichia coli. The sensitivities quantitatively identify rate limiting enzymes and are mathematically precise, unlike current finite difference based approaches used for sensitivity analysis. Further, efficient differentiation of constraint-based models unlocks the ability to use gradient information for parameter estimation. We demonstrate this by improving, genome-wide, the state-of-the-art turnover number estimates for E. coli. Finally, we show that this technique can be generalized to arbitrarily complex models. By differentiating the optimal solution of a model incorporating both thermodynamic and kinetic rate equations, the effect of metabolite concentrations on biomass growth can be elucidated. We benchmark these metabolite sensitivities against a large experimental gene knockdown study, and find good alignment between the predicted sensitivities and in vivo metabolome changes. In sum, we demonstrate several applications of differentiating optimal solutions of constraint-based metabolic models, and show how it connects to classic metabolic control analysis.

Thermodynamic limitations of PHB production from formate and fructose in Cupriavidus necator

The chemolithotroph Cupriavidus necator H16 is known as a natural producer of the bioplastic-polymer PHB, as well as for its metabolic versatility to utilize different substrates, including formate as the sole carbon and energy source. Depending on the entry point of the substrate, this versatility requires adjustment of the thermodynamic landscape to maintain sufficiently high driving forces for biological processes. Here we employed a model of the core metabolism of C. necator H16 to analyze the thermodynamic driving forces and PHB yields from formate for different metabolic engineering strategies. For this, we enumerated elementary flux modes (EFMs) of the network and evaluated their PHB yields as well as thermodynamics via Max-min driving force (MDF) analysis and random sampling of driving forces. A heterologous ATP:citrate lyase reaction was predicted to increase driving force for producing acetyl-CoA. A heterologous phosphoketolase reaction was predicted to increase maximal PHB yields as well as driving forces. These enzymes were then verified experimentally to enhance PHB titers between 60 and 300% in select conditions. The EFM analysis also revealed that PHB production from formate may be limited by low driving forces through citrate lyase and aconitase, as well as cofactor balancing, and identified additional reactions associated with low and high PHB yield. Proteomics analysis of the engineered strains confirmed an increased abundance of aconitase and cofactor balancing. The findings of this study aid in understanding metabolic adaptation. Furthermore, the outlined approach will be useful in designing metabolic engineering strategies in other non-model bacteria.

Nonlinear multi-objective flux balance analysis of the Warburg Effect

Due to its implication in cancer treatment, the Warburg Effect has received extensive in silico investigation. Flux Balance Analysis (FBA), based on constrained optimization, was successfully applied in the Warburg Effect modelling. Yet, the assumption that cell types have one invariant cellular objective severely limits the applicability of the previous FBA models. Meanwhile, we note that cell types with different objectives show different extents of the Warburg Effect. To extend the applicability of the previous model and model the disparate cellular pathway preferences in different cell types, we built a Nonlinear Multi-Objective FBA (NLMOFBA) model by including three key objective terms (ATP production rate, lactate generation rate and ATP yield) into one objective function through linear scalarization. By constructing a cellular objective map and iteratively varying the objective weights, we showed disparate cellular pathway preferences manifested by different cell types driven by their unique cellular objectives, and we gained insights about the causal relationship between cellular objectives and the Warburg Effect. In addition, we obtained other biologically consistent results by using our NLMOFBA model. For example, augmented with the constraint associated with inefficient mitochondria function, low oxygen availability, or limited substrate, NLMOFBA predicts cellular pathways supported by the biology literature. Collectively, our NLMOFBA model can help build a complete understanding towards the Warburg Effect in different cell types. Finally, we investigated the impact of glutaminolysis, an important pathway related to glycolysis, on the occurrence of the Warburg Effect by using linear programming.

Soil enzymes as indicators of soil function: A step toward greater realism in microbial ecological modeling

Soil carbon (C) and nitrogen (N) cycles and their complex responses to environmental changes have received increasing attention. However, large uncertainties in model predictions remain, partially due to the lack of explicit representation and parameterization of microbial processes. One great challenge is to effectively integrate rich microbial functional traits into ecosystem modeling for better predictions. Here, using soil enzymes as indicators of soil function, we developed a competitive dynamic enzyme allocation scheme and detailed enzyme-mediated soil inorganic N processes in the Microbial-ENzyme Decomposition (MEND) model. We conducted a rigorous calibration and validation of MEND with diverse soil C-N fluxes, microbial C:N ratios, and functional gene abundances from a 12-year CO₂  × N grassland experiment (BioCON) in Minnesota, USA. In addition to accurately simulating soil CO₂ fluxes and multiple N variables, the model correctly predicted microbial C:N ratios and their negative response to enriched N supply. Model validation further showed that, compared to the changes in simulated enzyme concentrations and decomposition rates, the changes in simulated activities of eight C-N-associated enzymes were better explained by the measured gene abundances in responses to elevated atmospheric CO₂ concentration. Our results demonstrated that using enzymes as indicators of soil function and validating model predictions with functional gene abundances in ecosystem modeling can provide a basis for testing hypotheses about microbially mediated biogeochemical processes in response to environmental changes. Further development and applications of the modeling framework presented here will enable microbial ecologists to address ecosystem-level questions beyond empirical observations, toward more predictive understanding, an ultimate goal of microbial ecology.

Elementary vectors and autocatalytic sets for resource allocation in next-generation models of cellular growth

Traditional (genome-scale) metabolic models of cellular growth involve an approximate biomass “reaction”, which specifies biomass composition in terms of precursor metabolites (such as amino acids and nucleotides). On the one hand, biomass composition is often not known exactly and may vary drastically between conditions and strains. On the other hand, the predictions of computational models crucially depend on biomass. Also elementary flux modes (EFMs), which generate the flux cone, depend on the biomass reaction. To better understand cellular phenotypes across growth conditions, we introduce and analyze new classes of elementary vectors for comprehensive (next-generation) metabolic models, involving explicit synthesis reactions for all macromolecules. Elementary growth modes (EGMs) are given by stoichiometry and generate the growth cone. Unlike EFMs, they are not support-minimal, in general, but cannot be decomposed “without cancellations”. In models with additional (capacity) constraints, elementary growth vectors (EGVs) generate a growth polyhedron and depend also on growth rate. However, EGMs/EGVs do not depend on the biomass composition. In fact, they cover all possible biomass compositions and can be seen as unbiased versions of elementary flux modes/vectors (EFMs/EFVs) used in traditional models. To relate the new concepts to other branches of theory, we consider autocatalytic sets of reactions. Further, we illustrate our results in a small model of a self-fabricating cell, involving glucose and ammonium uptake, amino acid and lipid synthesis, and the expression of all enzymes and the ribosome itself. In particular, we study the variation of biomass composition as a function of growth rate. In agreement with experimental data, low nitrogen uptake correlates with high carbon (lipid) storage.

Next-generation, genome-scale metabolic models allow to study the reallocation of cellular resources upon changing environmental conditions, by not only modeling flux distributions, but also expression profiles of the catalyzing proteome. In particular, they do no longer assume a fixed biomass composition. Methods to identify optimal solutions in such comprehensive models exist, however, an unbiased understanding of all feasible allocations is missing so far. Here we develop new concepts, called elementary growth modes and vectors, that provide a generalized definition of minimal pathways, thereby extending classical elementary flux modes (used in traditional models with a fixed biomass composition). The new concepts provide an understanding of all possible flux distributions and of all possible biomass compositions. In other words, elementary growth modes and vectors are the unique functional units in any comprehensive model of cellular growth. As an example, we show that lipid accumulation upon nitrogen starvation is a consequence of resource allocation and does not require active regulation. Our work puts current approaches on a theoretical basis and allows to seamlessly transfer existing workflows (e.g. for the design of cell factories) to next-generation metabolic models.

Selection for Cell Yield Does Not Reduce Overflow Metabolism in Escherichia coli

Overflow metabolism is ubiquitous in nature, and it is often considered inefficient because it leads to a relatively low biomass yield per consumed carbon. This metabolic strategy has been described as advantageous because it supports high growth rates during nutrient competition. Here, we experimentally evolved bacteria without nutrient competition by repeatedly growing and mixing millions of parallel batch cultures of Escherichia coli. Each culture originated from a water-in-oil emulsion droplet seeded with a single cell. Unexpectedly we found that overflow metabolism (acetate production) did not change. Instead, the numerical cell yield during the consumption of the accumulated acetate increased as a consequence of a reduction in cell size. Our experiments and a mathematical model show that fast growth and overflow metabolism, followed by the consumption of the overflow metabolite, can lead to a higher numerical cell yield and therefore a higher fitness compared with full respiration of the substrate. This provides an evolutionary scenario where overflow metabolism can be favorable even in the absence of nutrient competition.

Searching for principles of microbial physiology

Why do evolutionarily distinct microorganisms display similar physiological behaviours? Why are transitions from high-ATP yield to low(er)-ATP yield metabolisms so widespread across species? Why is fast growth generally accompanied with low stress tolerance? Do these regularities occur because most microbial species are subject to the same selective pressures and physicochemical constraints? If so, a broadly-applicable theory might be developed that predicts common microbiological behaviours. Microbial systems biologists have been working out the contours of this theory for the last two decades, guided by experimental data. At its foundations lie basic principles from evolutionary biology, enzyme biochemistry, metabolism, cell composition and steady-state growth. The theory makes predictions about fitness costs and benefits of protein expression, physicochemical constraints on cell growth and characteristics of optimal metabolisms that maximise growth rate. Comparisons of the theory with experimental data indicates that microorganisms often aim for maximisation of growth rate, also in the presence of stresses; they often express optimal metabolisms and metabolic proteins at optimal concentrations. This review explains the current status of the theory for microbiologists; its roots, predictions, experimental evidence and future directions.

Answer Set Programming for Computing Constraints-Based Elementary Flux Modes: Application to Escherichia coli Core Metabolism

Elementary Flux Modes (EFMs) provide a rigorous basis to systematically characterize the steady state, cellular phenotypes, as well as metabolic network robustness and fragility. However, the number of EFMs typically grows exponentially with the size of the metabolic network, leading to excessive computational demands, and unfortunately, a large fraction of these EFMs are not biologically feasible due to system constraints. This combinatorial explosion often prevents the complete analysis of genome-scale metabolic models. Traditionally, EFMs are computed by the double description method, an efficient algorithm based on matrix calculation; however, only a few constraints can be integrated into this computation. They must be monotonic with regard to the set inclusion of the supports; otherwise, they must be treated in post-processing and thus do not save computational time. We present aspefm, a hybrid computational tool based on Answer Set Programming (ASP) and Linear Programming (LP) that permits the computation of EFMs while implementing many different types of constraints. We apply our methodology to the Escherichia coli core model, which contains 226×106 EFMs. In considering transcriptional and environmental regulation, thermodynamic constraints, and resource usage considerations, the solution space is reduced to 1118 EFMs that can be computed directly with aspefm. The solution set, for E. coli growth on O2 gradients spanning fully aerobic to anaerobic, can be further reduced to four optimal EFMs using post-processing and Pareto front analysis.

Dynamic metabolic resource allocation based on the maximum entropy principle

Organisms have evolved a variety of mechanisms to cope with the unpredictability of environmental conditions, and yet mainstream models of metabolic regulation are typically based on strict optimality principles that do not account for uncertainty. This paper introduces a dynamic metabolic modelling framework that is a synthesis of recent ideas on resource allocation and the powerful optimal control formulation of Ramkrishna and colleagues. In particular, their work is extended based on the hypothesis that cellular resources are allocated among elementary flux modes according to the principle of maximum entropy. These concepts both generalise and unify prior approaches to dynamic metabolic modelling by establishing a smooth interpolation between dynamic flux balance analysis and dynamic metabolic models without regulation. The resulting theory is successful in describing 'bet-hedging' strategies employed by cell populations dealing with uncertainty in a fluctuating environment, including heterogenous resource investment, accumulation of reserves in growth-limiting conditions, and the observed behaviour of yeast growing in batch and continuous cultures. The maximum entropy principle is also shown to yield an optimal control law consistent with partitioning resources between elementary flux mode families, which has important practical implications for model reduction, selection, and simulation.

An analytical theory of balanced cellular growth

The biological fitness of microbes is largely determined by the rate with which they replicate their biomass composition. Mathematical models that maximize this balanced growth rate while accounting for mass conservation, reaction kinetics, and limits on dry mass per volume are inevitably non-linear. Here, we develop a general theory for such models, termed Growth Balance Analysis (GBA), which provides explicit expressions for protein concentrations, fluxes, and growth rates. These variables are functions of the concentrations of cellular components, for which we calculate marginal fitness costs and benefits that are related to metabolic control coefficients. At maximal growth rate, the net benefits of all concentrations are equal. Based solely on physicochemical constraints, GBA unveils fundamental quantitative principles of cellular resource allocation and growth; it accurately predicts the relationship between growth rates and ribosome concentrations in E. coli and yeast and between growth rate and dry mass density in E. coli.

Thermodynamic Approaches in Flux Analysis

Networks of reactions inside the cell are constrained by the laws of mass and energy balance. Constrained-based modelling (CBM) is the most used method to describe the mass balance of metabolic network. The main key concepts in CBM are stoichiometric analysis such as elementary flux mode analysis or flux balance analysis. Some of these methods have focused on adding thermodynamics constraints to eliminate non-physical fluxes or inconsistencies in the metabolic system. Here, we review the main different approaches and how they tackle the different class of problems.

Elementary Growth Modes provide a molecular description of cellular self-fabrication

In this paper we try to describe all possible molecular states (phenotypes) for a cell that fabricates itself at a constant rate, given its enzyme kinetics and the stoichiometry of all reactions. For this, we must understand the process of cellular growth: steady-state self-fabrication requires a cell to synthesize all of its components, including metabolites, enzymes and ribosomes, in proportions that match its own composition. Simultaneously, the concentrations of these components affect the rates of metabolism and biosynthesis, and hence the growth rate. We here derive a theory that describes all phenotypes that solve this circular problem. All phenotypes can be described as a combination of minimal building blocks, which we call Elementary Growth Modes (EGMs). EGMs can be used as the theoretical basis for all models that explicitly model self-fabrication, such as the currently popular Metabolism and Expression models. We then use our theory to make concrete biological predictions. We find that natural selection for maximal growth rate drives microorganisms to states of minimal phenotypic complexity: only one EGM will be active when growth rate is maximised. The phenotype of a cell is only extended with one more EGM whenever growth becomes limited by an additional biophysical constraint, such as a limited solvent capacity of a cellular compartment. The theory presented here extends recent results on Elementary Flux Modes: the minimal building blocks of cellular growth models that lack the self-fabrication aspect. Our theory starts from basic biochemical and evolutionary considerations, and describes unicellular life, both in growth-promoting and in stress-inducing environments, in terms of EGMs.

Regulation underlying hierarchical and simultaneous utilization of carbon substrates by flux sensors in Escherichia coli

Many microbes exhibit nutrient preferences, exemplified by the “hierarchical” consumption of certain carbon substrates. Here we systematically investigate under which physiological conditions hierarchical substrate utilization occurs and its mechanisms of implementation. We show utilization hierarchy of Escherichia coli to be ordered by the carbon-uptake flux rather than the identity of the substrates. A detailed study of glycerol uptake finds that it is fully suppressed if the uptake flux of another glycolytic substrate exceeds a threshold, set to the influx obtained when grown on glycerol alone. Below this threshold, limited glycerol uptake is “supplemented” such that the total carbon uptake is maintained at the threshold. This behavior results from total-flux feedback mediated by cAMP-Crp signaling, but also requires inhibition by regulator fructose-1,6-bisphosphate, which senses the upper glycolytic flux and ensures that glycerol uptake defers to other glycolytic substrates but not to gluconeogenic ones. A quantitative model reproduces all observed utilization patterns including those of key mutants. The proposed mechanism relies on differential regulation of uptake enzymes and requires a specific operon organization. This organization is found conserved across related species for several uptake systems, suggesting the deployment of similar mechanisms for hierarchical substrate utilization by a spectrum of microbes.

Modelling microbial communities using biochemical resource allocation analysis

To understand the functioning and dynamics of microbial communities is a fundamental challenge in current biology. To tackle this challenge, the construction of computational models of interacting microbes is an indispensable tool. There is, however, a large chasm between ecologically motivated descriptions of microbial growth used in many current ecosystems simulations, and detailed metabolic pathway and genome-based descriptions developed in the context of systems and synthetic biology. Here, we seek to demonstrate how resource allocation models of microbial growth offer the potential to advance ecosystem simulations and their parametrization. In particular, recent work on quantitative resource allocation allow us to formulate mechanistic models of microbial growth that are physiologically meaningful while remaining computationally tractable. These models go beyond Michaelis-Menten and Monod-type growth models, and are capable of accounting for emergent properties that underlie the remarkable plasticity of microbial growth. We outline the utility and advantages of using biochemical resource allocation models by considering a coarse-grained model of cyanobacterial growth and demonstrate how the model allows us to address specific questions of relevance for the simulation of marine microbial ecosystems, including the physiological acclimation of protein expression to different environments, the description of co-limitation by several nutrients and the differential use of alternative nutrient sources, as well as the description of metabolic diversity based on our increasing knowledge about quantitative cell physiology.

Growth strategy of microbes on mixed carbon sources

A classic problem in microbiology is that bacteria display two types of growth behavior when cultured on a mixture of two carbon sources: the two sources are sequentially consumed one after another (diauxie) or they are simultaneously consumed (co-utilization). The search for the molecular mechanism of diauxie led to the discovery of the lac operon. However, questions remain as why microbes would bother to have different strategies of taking up nutrients. Here we show that diauxie versus co-utilization can be understood from the topological features of the metabolic network. A model of optimal allocation of protein resources quantitatively explains why and how the cell makes the choice. In case of co-utilization, the model predicts the percentage of each carbon source in supplying the amino acid pools, which is quantitatively verified by experiments. Our work solves a long-standing puzzle and provides a quantitative framework for the carbon source utilization of microbes.

Bacteria grown on two carbon sources either consume both sources simultaneously or consume them sequentially. Here the authors use a metabolic network model of E. coli to show that optimal protein resource allocation and topological features of the network can explain the choice of carbon acquisition.

The number of active metabolic pathways is bounded by the number of cellular constraints at maximal metabolic rates

Growth rate is a near-universal selective pressure across microbial species. High growth rates require hundreds of metabolic enzymes, each with different nonlinear kinetics, to be precisely tuned within the bounds set by physicochemical constraints. Yet, the metabolic behaviour of many species is characterized by simple relations between growth rate, enzyme expression levels and metabolic rates. We asked if this simplicity could be the outcome of optimisation by evolution. Indeed, when the growth rate is maximized-in a static environment under mass-conservation and enzyme expression constraints-we prove mathematically that the resulting optimal metabolic flux distribution is described by a limited number of subnetworks, known as Elementary Flux Modes (EFMs). We show that, because EFMs are the minimal subnetworks leading to growth, a small active number automatically leads to the simple relations that are measured. We find that the maximal number of flux-carrying EFMs is determined only by the number of imposed constraints on enzyme expression, not by the size, kinetics or topology of the network. This minimal-EFM extremum principle is illustrated in a graphical framework, which explains qualitative changes in microbial behaviours, such as overflow metabolism and co-consumption, and provides a method for identification of the enzyme expression constraints that limit growth under the prevalent conditions. The extremum principle applies to all microorganisms that are selected for maximal growth rates under protein concentration constraints, for example the solvent capacities of cytosol, membrane or periplasmic space.

Maximum entropy and population heterogeneity in continuous cell cultures

Continuous cultures of mammalian cells are complex systems displaying hallmark phenomena of nonlinear dynamics, such as multi-stability, hysteresis, as well as sharp transitions between different metabolic states. In this context mathematical models may suggest control strategies to steer the system towards desired states. Although even clonal populations are known to exhibit cell-to-cell variability, most of the currently studied models assume that the population is homogeneous. To overcome this limitation, we use the maximum entropy principle to model the phenotypic distribution of cells in a chemostat as a function of the dilution rate. We consider the coupling between cell metabolism and extracellular variables describing the state of the bioreactor and take into account the impact of toxic byproduct accumulation on cell viability. We present a formal solution for the stationary state of the chemostat and show how to apply it in two examples. First, a simplified model of cell metabolism where the exact solution is tractable, and then a genome-scale metabolic network of the Chinese hamster ovary (CHO) cell line. Along the way we discuss several consequences of heterogeneity, such as: qualitative changes in the dynamical landscape of the system, increasing concentrations of byproducts that vanish in the homogeneous case, and larger population sizes.

OptMDFpathway: Identification of metabolic pathways with maximal thermodynamic driving force and its application for analyzing the endogenous CO2 fixation potential of Escherichia coli

Constraint-based modeling techniques have become a standard tool for the in silico analysis of metabolic networks. To further improve their accuracy, recent methodological developments focused on integration of thermodynamic information in metabolic models to assess the feasibility of flux distributions by thermodynamic driving forces. Here we present OptMDFpathway, a method that extends the recently proposed framework of Max-min Driving Force (MDF) for thermodynamic pathway analysis. Given a metabolic network model, OptMDFpathway identifies both the optimal MDF for a desired phenotypic behavior as well as the respective pathway itself that supports the optimal driving force. OptMDFpathway is formulated as a mixed-integer linear program and is applicable to genome-scale metabolic networks. As an important theoretical result, we also show that there exists always at least one elementary mode in the network that reaches the maximal MDF. We employed our new approach to systematically identify all substrate-product combinations in Escherichia coli where product synthesis allows for concomitant net CO₂ assimilation via thermodynamically feasible pathways. Although biomass synthesis cannot be coupled to net CO₂ fixation in E. coli we found that as many as 145 of the 949 cytosolic carbon metabolites contained in the genome-scale model iJO1366 enable net CO₂ incorporation along thermodynamically feasible pathways with glycerol as substrate and 34 with glucose. The most promising products in terms of carbon assimilation yield and thermodynamic driving forces are orotate, aspartate and the C4-metabolites of the tricarboxylic acid cycle. We also identified thermodynamic bottlenecks frequently limiting the maximal driving force of the CO₂-fixing pathways. Our results indicate that heterotrophic organisms like E. coli hold a possibly underestimated potential for CO₂ assimilation which may complement existing biotechnological approaches for capturing CO₂. Furthermore, we envision that the developed OptMDFpathway approach can be used for many other applications within the framework of constrained-based modeling and for rational design of metabolic networks.

When analyzing metabolic networks, one often searches for metabolic pathways with certain (desired) properties, for example, conversion routes that maximize the yield of a product from a given substrate. While those problems can be solved with established methods of constraint-based modeling, no algorithm is currently available for genome-scale models to identify the pathway that has the highest possible thermodynamic driving force among all solutions with predefined stoichiometric properties. This gap is closed with our new approach OptMDFpathway which is based on the recently introduced concept of Max-min Driving Force (MDF). OptMDFpathway offers various applications, especially in the context of metabolic design of cell factories. To demonstrate the power and usefulness of OptMDFpathway, we employed it to analyze the endogenous CO₂ fixation potential of Escherichia coli. While E. coli cannot assimilate CO₂ into biomass, net CO₂ fixation can take place along linear pathways from substrate to product and we show that thermodynamically feasible pathways with net CO₂ assimilation exist for 145 (34) products when choosing glycerol (glucose) as substrate. Our results indicate that heterotrophic organisms like E. coli hold a possibly underestimated potential for CO₂ assimilation which may complement existing biotechnological approaches for capturing CO₂.

Maintaining maximal metabolic flux by gene expression control

One of the marvels of biology is the phenotypic plasticity of microorganisms. It allows them to maintain high growth rates across conditions. Studies suggest that cells can express metabolic enzymes at tuned concentrations through adjustment of gene expression. The associated transcription factors are often regulated by intracellular metabolites. Here we study metabolite-mediated regulation of metabolic-gene expression that maximises metabolic fluxes across conditions. We developed an adaptive control theory, qORAC (for ‘Specific Flux (q) Optimization by Robust Adaptive Control’), and illustrate it with several examples of metabolic pathways. The key feature of the theory is that it does not require knowledge of the regulatory network, only of the metabolic part. We derive that maximal metabolic flux can be maintained in the face of varying N environmental parameters only if the number of transcription-factor binding metabolites is at least equal to N. The controlling circuits appear to require simple biochemical kinetics. We conclude that microorganisms likely can achieve maximal rates in metabolic pathways, in the face of environmental changes.

To attain high growth rates, microorganisms need to sustain high activities of metabolic reactions. Since the catalysing enzymes are in finite supply, cells need to carefully tune their concentrations. When conditions change, cells need to adjust those concentrations. How cells maintain high metabolism rates across conditions by way of gene regulatory mechanisms and whether they can maximise metabolic activity is far from clear. Here we present a general theory that solves this metabolic control problem, which we have called qORAC for specific flux (q) Optimisation by Robust Adaptive Control. It considers that external changes are sensed by internal “sensor” metabolites that bind to transcription factors in order to regulate enzyme-synthesis rates. We show that such a combined system of metabolism and its gene network can self-optimise its metabolic activity across conditions. We present the mathematical conditions for the required adaptive control for robust system-steering to optimal states across conditions. We provide explicit examples of such self-optimising coupled metabolism and gene network systems. We prove that a cell can be robust to changes in K parameters, e.g. external conditions, if at least K internal metabolite concentrations act transcription-factor binding sensors. We find that the optimal relation of the enzyme synthesis rates of self-optimising systems and the concentration of the sensor metabolites can generally be implemented by basic biochemistry. Our results indicate how cells are able to maintain maximal reaction rates, even in changing conditions.

Linear programming model can explain respiration of fermentation products

Many differentiated cells rely primarily on mitochondrial oxidative phosphorylation for generating energy in the form of ATP needed for cellular metabolism. In contrast most tumor cells instead rely on aerobic glycolysis leading to lactate to about the same extent as on respiration. Warburg found that cancer cells to support oxidative phosphorylation, tend to ferment glucose or other energy source into lactate even in the presence of sufficient oxygen, which is an inefficient way to generate ATP. This effect also occurs in striated muscle cells, activated lymphocytes and microglia, endothelial cells and several mammalian cell types, a phenomenon termed the "Warburg effect". The effect is paradoxical at first glance because the ATP production rate of aerobic glycolysis is much slower than that of respiration and the energy demands are better to be met by pure oxidative phosphorylation. We tackle this question by building a minimal model including three combined reactions. The new aspect in extension to earlier models is that we take into account the possible uptake and oxidation of the fermentation products. We examine the case where the cell can allocate protein on several enzymes in a varying distribution and model this by a linear programming problem in which the objective is to maximize the ATP production rate under different combinations of constraints on enzymes. Depending on the cost of reactions and limitation of the substrates, this leads to pure respiration, pure fermentation, and a mixture of respiration and fermentation. The model predicts that fermentation products are only oxidized when glucose is scarce or its uptake is severely limited.

Metabolic enzyme cost explains variable trade-offs between microbial growth rate and yield

Microbes may maximize the number of daughter cells per time or per amount of nutrients consumed. These two strategies correspond, respectively, to the use of enzyme-efficient or substrate-efficient metabolic pathways. In reality, fast growth is often associated with wasteful, yield-inefficient metabolism, and a general thermodynamic trade-off between growth rate and biomass yield has been proposed to explain this. We studied growth rate/yield trade-offs by using a novel modeling framework, Enzyme-Flux Cost Minimization (EFCM) and by assuming that the growth rate depends directly on the enzyme investment per rate of biomass production. In a comprehensive mathematical model of core metabolism in E. coli, we screened all elementary flux modes leading to cell synthesis, characterized them by the growth rates and yields they provide, and studied the shape of the resulting rate/yield Pareto front. By varying the model parameters, we found that the rate/yield trade-off is not universal, but depends on metabolic kinetics and environmental conditions. A prominent trade-off emerges under oxygen-limited growth, where yield-inefficient pathways support a 2-to-3 times higher growth rate than yield-efficient pathways. EFCM can be widely used to predict optimal metabolic states and growth rates under varying nutrient levels, perturbations of enzyme parameters, and single or multiple gene knockouts.

Resource allocation in living organisms

Quantitative prediction of resource allocation for living systems has been an intensive area of research in the field of biology. Resource allocation was initially investigated in higher organisms by using empirical mathematical models based on mass distribution. A challenge is now to go a step further by reconciling the cellular scale to the individual scale. In the present paper, we review the foundations of modelling of resource allocation, particularly at the cellular scale: from small macro-molecular models to genome-scale cellular models. We enlighten how the combination of omic measurements and computational advances together with systems biology has contributed to dramatic progresses in the current understanding and prediction of cellular resource allocation. Accurate genome-wide predictive methods of resource allocation based on the resource balance analysis (RBA) framework have been developed and ensure a good trade-off between the complexity/tractability and the prediction capability of the model. The RBA framework shows promise for a wide range of applications in metabolic engineering and synthetic biology, and for pursuing investigations of the design principles of cellular and multi-cellular organisms.

From elementary flux modes to elementary flux vectors: Metabolic pathway analysis with arbitrary linear flux constraints

Elementary flux modes (EFMs) emerged as a formal concept to describe metabolic pathways and have become an established tool for constraint-based modeling and metabolic network analysis. EFMs are characteristic (support-minimal) vectors of the flux cone that contains all feasible steady-state flux vectors of a given metabolic network. EFMs account for (homogeneous) linear constraints arising from reaction irreversibilities and the assumption of steady state; however, other (inhomogeneous) linear constraints, such as minimal and maximal reaction rates frequently used by other constraint-based techniques (such as flux balance analysis [FBA]), cannot be directly integrated. These additional constraints further restrict the space of feasible flux vectors and turn the flux cone into a general flux polyhedron in which the concept of EFMs is not directly applicable anymore. For this reason, there has been a conceptual gap between EFM-based (pathway) analysis methods and linear optimization (FBA) techniques, as they operate on different geometric objects. One approach to overcome these limitations was proposed ten years ago and is based on the concept of elementary flux vectors (EFVs). Only recently has the community started to recognize the potential of EFVs for metabolic network analysis. In fact, EFVs exactly represent the conceptual development required to generalize the idea of EFMs from flux cones to flux polyhedra. This work aims to present a concise theoretical and practical introduction to EFVs that is accessible to a broad audience. We highlight the close relationship between EFMs and EFVs and demonstrate that almost all applications of EFMs (in flux cones) are possible for EFVs (in flux polyhedra) as well. In fact, certain properties can only be studied with EFVs. Thus, we conclude that EFVs provide a powerful and unifying framework for constraint-based modeling of metabolic networks.

Mathematical models for explaining the Warburg effect: a review focussed on ATP and biomass production

For producing ATP, tumour cells rely on glycolysis leading to lactate to about the same extent as on respiration. Thus, the ATP synthesis flux from glycolysis is considerably higher than in the corresponding healthy cells. This is known as the Warburg effect (named after German biochemist Otto H. Warburg) and also applies to striated muscle cells, activated lymphocytes, microglia, endothelial cells and several other cell types. For similar phenomena in several yeasts and many bacteria, the terms Crabtree effect and overflow metabolism respectively, are used. The Warburg effect is paradoxical at first sight because the molar ATP yield of glycolysis is much lower than that of respiration. Although a straightforward explanation is that glycolysis allows a higher ATP production rate, the question arises why cells do not re-allocate protein to the high-yield pathway of respiration. Mathematical modelling can help explain this phenomenon. Here, we review several models at various scales proposed in the literature for explaining the Warburg effect. These models support the hypothesis that glycolysis allows for a higher proliferation rate due to increased ATP production and precursor supply rates.

Resource allocation in metabolic networks: kinetic optimization and approximations by FBA

Based on recent theoretical results on optimal flux distributions in kinetic metabolic networks, we explore the congruences and differences between solutions of kinetic optimization problems and results obtained by constraint-based methods. We demonstrate that, for a certain resource allocation problem, kinetic optimization and standard flux balance analysis (FBA) give rise to qualitatively different results. Furthermore, we introduce a variant of FBA, called satFBA, whose predictions are in qualitative agreement with kinetic optimization.

The steady-state assumption in oscillating and growing systems

The steady-state assumption, which states that the production and consumption of metabolites inside the cell are balanced, is one of the key aspects that makes an efficient analysis of genome-scale metabolic networks possible. It can be motivated from two different perspectives. In the time-scales perspective, we use the fact that metabolism is much faster than other cellular processes such as gene expression. Hence, the steady-state assumption is derived as a quasi-steady-state approximation of the metabolism that adapts to the changing cellular conditions. In this article we focus on the second perspective, stating that on the long run no metabolite can accumulate or deplete. In contrast to the first perspective it is not immediately clear how this perspective can be captured mathematically and what assumptions are required to obtain the steady-state condition. By presenting a mathematical framework based on the second perspective we demonstrate that the assumption of steady-state also applies to oscillating and growing systems without requiring quasi-steady-state at any time point. However, we also show that the average concentrations may not be compatible with the average fluxes. In summary, we establish a mathematical foundation for the steady-state assumption for long time periods that justifies its successful use in many applications. Furthermore, this mathematical foundation also pinpoints unintuitive effects in the integration of metabolite concentrations using nonlinear constraints into steady-state models for long time periods.

Elementary Vectors and Conformal Sums in Polyhedral Geometry and their Relevance for Metabolic Pathway Analysis

A fundamental result in metabolic pathway analysis states that every flux mode can be decomposed into a sum of elementary modes. However, only a decomposition without cancelations is biochemically meaningful, since a reversible reaction cannot have different directions in the contributing elementary modes. This essential requirement has been largely overlooked by the metabolic pathway community. Indeed, every flux mode can be decomposed into elementary modes without cancelations. The result is an immediate consequence of a theorem by Rockafellar which states that every element of a linear subspace is a conformal sum (a sum without cancelations) of elementary vectors (support-minimal vectors). In this work, we extend the theorem, first to “subspace cones” and then to general polyhedral cones and polyhedra. Thereby, we refine Minkowski's and Carathéodory's theorems, two fundamental results in polyhedral geometry. We note that, in general, elementary vectors need not be support-minimal; in fact, they are conformally non-decomposable and form a unique minimal set of conformal generators. Our treatment is mathematically rigorous, but suitable for systems biologists, since we give self-contained proofs for our results and use concepts motivated by metabolic pathway analysis. In particular, we study cones defined by linear subspaces and nonnegativity conditions — like the flux cone — and use them to analyze general polyhedral cones and polyhedra. Finally, we review applications of elementary vectors and conformal sums in metabolic pathway analysis.

Which sets of elementary flux modes form thermodynamically feasible flux distributions?

Elementary flux modes (EFMs) are non-decomposable steady-state fluxes through metabolic networks. Every possible flux through a network can be described as a superposition of EFMs. The definition of EFMs is based on the stoichiometry of the network, and it has been shown previously that not all EFMs are thermodynamically feasible. These infeasible EFMs cannot contribute to a biologically meaningful flux distribution. In this work, we show that a set of thermodynamically feasible EFMs need not be thermodynamically consistent. We use first principles of thermodynamics to define the feasibility of a flux distribution and present a method to compute the largest thermodynamically consistent sets (LTCSs) of EFMs. An LTCS contains the maximum number of EFMs that can be combined to form a thermodynamically feasible flux distribution. As a case study we analyze all LTCSs found in Escherichia coli when grown on glucose and show that only one LTCS shows the required phenotypical properties. Using our method, we find that in our E. coli model < 10% of all EFMs are thermodynamically relevant.

Mathematical models of cancer metabolism

Metabolism is essential for life, and its alteration is implicated in multiple human diseases. The transformation from a normal to a cancerous cell requires metabolic changes to fuel the high metabolic demands of cancer cells, including but not limited to cell proliferation and cell migration. In recent years, there have been a number of new discoveries connecting known aberrations in oncogenic and tumour suppressor pathways with metabolic alterations required to sustain cell proliferation and migration. However, an understanding of the selective advantage of these metabolic alterations is still lacking. Here, we review the literature on mathematical models of metabolism, with an emphasis on their contribution to the identification of the selective advantage of metabolic phenotypes that seem otherwise wasteful or accidental. We will show how the molecular hallmarks of cancer can be related to cell proliferation and tissue remodelling, the two major physiological requirements for the development of a multicellular structure. We will cover different areas such as genome-wide gene expression analysis, flux balance models, kinetic models, reaction diffusion models and models of the tumour microenvironment. We will also highlight current challenges and how their resolution will help to achieve a better understanding of cancer metabolism and the metabolic vulnerabilities of cancers.

Causes of upregulation of glycolysis in lymphocytes upon stimulation. A comparison with other cell types

In this review, we revisit the metabolic shift from respiration to glycolysis in lymphocytes upon activation, which is known as the Warburg effect in tumour cells. We compare the situation in lymphocytes with those in several other cell types, such as muscle cells, Kupffer cells, microglia cells, astrocytes, stem cells, tumour cells and various unicellular organisms (e.g. yeasts). We critically discuss and compare several explanations put forward in the literature for the observation that proliferating cells adopt this apparently less efficient pathway: hypoxia, poisoning of competitors by end products, higher ATP production rate, higher precursor supply, regulatory effects, and avoiding harmful effects (e.g. by reactive oxygen species). We conclude that in the case of lymphocytes, increased ATP production rate and precursor supply are the main advantages of upregulating glycolysis.

Metabolism at evolutionary optimal States

Metabolism is generally required for cellular maintenance and for the generation of offspring under conditions that support growth. The rates, yields (efficiencies), adaptation time and robustness of metabolism are therefore key determinants of cellular fitness. For biotechnological applications and our understanding of the evolution of metabolism, it is necessary to figure out how the functional system properties of metabolism can be optimized, via adjustments of the kinetics and expression of enzymes, and by rewiring metabolism. The trade-offs that can occur during such optimizations then indicate fundamental limits to evolutionary innovations and bioengineering. In this paper, we review several theoretical and experimental findings about mechanisms for metabolic optimization.

Elucidating the adaptation and temporal coordination of metabolic pathways using in-silico evolution

Cellular metabolism, the interconversion of small molecules by chemical reactions, is a tightly coordinated process that requires integration of diverse environmental and intracellular cues. While for many organisms the topology of the network of metabolic reactions is increasingly known, the regulatory principles that shape the network's adaptation to diverse and changing environments remain largely elusive. To investigate the principles of metabolic adaptation and regulation in metabolic pathways, we propose a computational approach based on in-silico evolution. Rather than analyzing existing regulatory schemes, we let a population of minimal, prototypical metabolic cells evolve rate constants and appropriate regulatory schemes that allow for optimal growth in static and fluctuating environments. Applying our approach to a small, but already sufficiently complex, minimal system reveals intricate transitions between metabolic modes. These results have implications for trade-offs in resource allocation. Going from static to varying environments, we show that for fluctuating nutrient availability, active metabolic regulation results in a significantly increased overall rate of metabolism.