An analytic approximation of the feasible space of metabolic networks
arXiv:1702.05400 · doi:10.1038/ncomms14915
Abstract
Assuming a steady-state condition within a cell, metabolic fluxes satisfy an under-determined linear system of stoichiometric equations. Characterizing the space of fluxes that satisfy such equations along with given bounds (and possibly additional relevant constraints) is considered of utmost importance for the understanding of cellular metabolism. Extreme values for each individual flux can be computed with Linear Programming (as Flux Balance Analysis), and their marginal distributions can be approximately computed with Monte-Carlo sampling. Here we present an approximate analytic method for the latter task based on Expectation Propagation equations that does not involve sampling and can achieve much better predictions than other existing analytic methods. The method is iterative, and its computation time is dominated by one matrix inversion per iteration. With respect to sampling, we show through extensive simulation that it has some advantages including computation time, and the ability to efficiently fix empirically estimated distributions of fluxes.
References in corpus (1)
Cited by in corpus (11)
- Statistical mechanics for metabolic networks during steady-state growth
- Characterizing steady states of genome-scale metabolic networks in continuous cell cultures
- Maximum entropy and population heterogeneity in continuous cell cultures
- Relationship between fitness and heterogeneity in exponentially growing microbial populations
- Compressed sensing reconstruction using Expectation Propagation
- Spin Glass Theory of Interacting Metabolic Networks
- The free lunch of a scale-free metabolism
- A maximum entropy network reconstruction of macroeconomic models
- Constraint-based inverse modeling of metabolic networks: a proof of concept
- Universal scaling relation and criticality in metabolism and growth of Escherichia coli
- Expectation propagation on the diluted Bayesian classifier