2,145 research outputs found
Flux imbalance analysis and the sensitivity of cellular growth to changes in metabolite pools
Stoichiometric models of metabolism, such as flux balance analysis (FBA), are classically applied to predicting steady state rates - or fluxes - of metabolic reactions in genome-scale metabolic networks. Here we revisit the central assumption of FBA, i.e. that intracellular metabolites are at steady state, and show that deviations from flux balance (i.e. flux imbalances) are informative of some features of in vivo metabolite concentrations. Mathematically, the sensitivity of FBA to these flux imbalances is captured by a native feature of linear optimization, the dual problem, and its corresponding variables, known as shadow prices. First, using recently published data on chemostat growth of Saccharomyces cerevisae under different nutrient limitations, we show that shadow prices anticorrelate with experimentally measured degrees of growth limitation of intracellular metabolites. We next hypothesize that metabolites which are limiting for growth (and thus have very negative shadow price) cannot vary dramatically in an uncontrolled way, and must respond rapidly to perturbations. Using a collection of published datasets monitoring the time-dependent metabolomic response of Escherichia coli to carbon and nitrogen perturbations, we test this hypothesis and find that metabolites with negative shadow price indeed show lower temporal variation following a perturbation than metabolites with zero shadow price. Finally, we illustrate the broader applicability of flux imbalance analysis to other constraint-based methods. In particular, we explore the biological significance of shadow prices in a constraint-based method for integrating gene expression data with a stoichiometric model. In this case, shadow prices point to metabolites that should rise or drop in concentration in order to increase consistency between flux predictions and gene expression data. In general, these results suggest that the sensitivity of metabolic optima to violations of the steady state constraints carries biologically significant information on the processes that control intracellular metabolites in the cell.Published versio
Constrained Allocation Flux Balance Analysis
New experimental results on bacterial growth inspire a novel top-down
approach to study cell metabolism, combining mass balance and proteomic
constraints to extend and complement Flux Balance Analysis. We introduce here
Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic
costs associated to growth are accounted for in an effective way through a
single additional genome-wide constraint. Its roots lie in the experimentally
observed pattern of proteome allocation for metabolic functions, allowing to
bridge regulation and metabolism in a transparent way under the principle of
growth-rate maximization. We provide a simple method to solve CAFBA efficiently
and propose an "ensemble averaging" procedure to account for unknown protein
costs. Applying this approach to modeling E. coli metabolism, we find that, as
the growth rate increases, CAFBA solutions cross over from respiratory,
growth-yield maximizing states (preferred at slow growth) to fermentative
states with carbon overflow (preferred at fast growth). In addition, CAFBA
allows for quantitatively accurate predictions on the rate of acetate excretion
and growth yield based on only 3 parameters determined by empirical growth
laws.Comment: 21 pages, 6 figures (main) + 33 pages, various figures and tables
(supporting); for the supplementary MatLab code, see
http://tinyurl.com/h763es
Enumeration of minimal stoichiometric precursor sets in metabolic networks
Background: What an organism needs at least from its environment to produce a set of metabolites, e.g. target(s) of interest and/or biomass, has been called a minimal precursor set. Early approaches to enumerate all minimal precursor sets took into account only the topology of the metabolic network (topological precursor sets). Due to cycles and the stoichiometric values of the reactions, it is often not possible to produce the target(s) from a topological precursor set in the sense that there is no feasible flux. Although considering the stoichiometry makes the problem harder, it enables to obtain biologically reasonable precursor sets that we call stoichiometric. Recently a method to enumerate all minimal stoichiometric precursor sets was proposed in the literature. The relationship between topological and stoichiometric precursor sets had however not yet been studied.
Results: Such relationship between topological and stoichiometric precursor sets is highlighted. We also present two algorithms that enumerate all minimal stoichiometric precursor sets. The first one is of theoretical interest only and is based on the above mentioned relationship. The second approach solves a series of mixed integer linear programming problems. We compared the computed minimal precursor sets to experimentally obtained growth media of several Escherichia coli strains using genome-scale metabolic networks.
Conclusions: The results show that the second approach efficiently enumerates minimal precursor sets taking stoichiometry into account, and allows for broad in silico studies of strains or species interactions that may help to understand e.g. pathotype and niche-specific metabolic capabilities. sasita is written in Java, uses cplex as LP solver and can be downloaded together with all networks and input files used in this paper at http://www.sasita.gforge.inria.fr
Construction and completion of flux balance models from pathway databases
Motivation: Flux balance analysis (FBA) is a well-known technique for genome-scale modeling of metabolic flux. Typically, an FBA formulation requires the accurate specification of four sets: biochemical reactions, biomass metabolites, nutrients and secreted metabolites. The development of FBA models can be time consuming and tedious because of the difficulty in assembling completely accurate descriptions of these sets, and in identifying errors in the composition of these sets. For example, the presence of a single non-producible metabolite in the biomass will make the entire model infeasible. Other difficulties in FBA modeling are that model distributions, and predicted fluxes, can be cryptic and difficult to understand
Characterizing steady states of genome-scale metabolic networks in continuous cell cultures
We present a model for continuous cell culture coupling intra-cellular
metabolism to extracellular variables describing the state of the bioreactor,
taking into account the growth capacity of the cell and the impact of toxic
byproduct accumulation. We provide a method to determine the steady states of
this system that is tractable for metabolic networks of arbitrary complexity.
We demonstrate our approach in a toy model first, and then in a genome-scale
metabolic network of the Chinese hamster ovary cell line, obtaining results
that are in qualitative agreement with experimental observations. More
importantly, we derive a number of consequences from the model that are
independent of parameter values. First, that the ratio between cell density and
dilution rate is an ideal control parameter to fix a steady state with desired
metabolic properties invariant across perfusion systems. This conclusion is
robust even in the presence of multi-stability, which is explained in our model
by the negative feedback loop on cell growth due to toxic byproduct
accumulation. Moreover, a complex landscape of steady states in continuous cell
culture emerges from our simulations, including multiple metabolic switches,
which also explain why cell-line and media benchmarks carried out in batch
culture cannot be extrapolated to perfusion. On the other hand, we predict
invariance laws between continuous cell cultures with different parameters. A
practical consequence is that the chemostat is an ideal experimental model for
large-scale high-density perfusion cultures, where the complex landscape of
metabolic transitions is faithfully reproduced. Thus, in order to actually
reflect the expected behavior in perfusion, performance benchmarks of
cell-lines and culture media should be carried out in a chemostat
Method for finding metabolic properties based on the general growth law. Liver examples. A General framework for biological modeling
We propose a method for finding metabolic parameters of cells, organs and
whole organisms, which is based on the earlier discovered general growth law.
Based on the obtained results and analysis of available biological models, we
propose a general framework for modeling biological phenomena and discuss how
it can be used in Virtual Liver Network project. The foundational idea of the
study is that growth of cells, organs, systems and whole organisms, besides
biomolecular machinery, is influenced by biophysical mechanisms acting at
different scale levels. In particular, the general growth law uniquely defines
distribution of nutritional resources between maintenance needs and biomass
synthesis at each phase of growth and at each scale level. We exemplify the
approach considering metabolic properties of growing human and dog livers and
liver transplants. A procedure for verification of obtained results has been
introduced too. We found that two examined dogs have high metabolic rates
consuming about 0.62 and 1 gram of nutrients per cubic centimeter of liver per
day, and verified this using the proposed verification procedure. We also
evaluated consumption rate of nutrients in human livers, determining it to be
about 0.088 gram of nutrients per cubic centimeter of liver per day for males,
and about 0.098 for females. This noticeable difference can be explained by
evolutionary development, which required females to have greater liver
processing capacity to support pregnancy. We also found how much nutrients go
to biomass synthesis and maintenance at each phase of liver and liver
transplant growth. Obtained results demonstrate that the proposed approach can
be used for finding metabolic characteristics of cells, organs, and whole
organisms, which can further serve as important inputs for many applications in
biology (protein expression), biotechnology (synthesis of substances), and
medicine.Comment: 20 pages, 6 figures, 4 table
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
Computing with bacterial constituents, cells and populations: from bioputing to bactoputing
The relevance of biological materials and processes to computing—aliasbioputing—has been explored for decades. These materials include DNA, RNA and proteins, while the processes include transcription, translation, signal transduction and regulation. Recently, the use of bacteria themselves as living computers has been explored but this use generally falls within the classical paradigm of computing. Computer scientists, however, have a variety of problems to which they seek solutions, while microbiologists are having new insights into the problems bacteria are solving and how they are solving them. Here, we envisage that bacteria might be used for new sorts of computing. These could be based on the capacity of bacteria to grow, move and adapt to a myriad different fickle environments both as individuals and as populations of bacteria plus bacteriophage. New principles might be based on the way that bacteria explore phenotype space via hyperstructure dynamics and the fundamental nature of the cell cycle. This computing might even extend to developing a high level language appropriate to using populations of bacteria and bacteriophage. Here, we offer a speculative tour of what we term bactoputing, namely the use of the natural behaviour of bacteria for calculating
- …