Abstract:Abstract Background Cellular metabolism is one of the most investigated system of biological interactions. While the topological nature of individual reactions and pathways in the network is quite well understood there is still a lack of comprehension regarding the global functional behavior of the system. In the last few years flux-balance analysis (FBA) has been the most successful and widely used technique for studying metabolism at system level. This method strongly relies on the hypothesis that the organism maximizes an objective function. However only under very specific biological conditions ( e.g . maximization of biomass for E. coli in reach nutrient medium) the cell seems to obey such optimization law. A more refined analysis not assuming extremization remains an elusive task for large metabolic systems due to algorithmic limitations. Results In this work we propose a novel algorithmic strategy that provides an efficient characterization of the whole set of stable fluxes compatible with the metabolic constraints. Using a technique derived from the fields of statistical physics and information theory we designed a message-passing algorithm to estimate the size of the affine space containing all possible steady-state flux distributions of metabolic networks. The algorithm, based on the well known Bethe approximation, can be used to approximately compute the volume of a non full-dimensional convex polytope in high dimensions. We first compare the accuracy of the predictions with an exact algorithm on small random metabolic networks. We also verify that the predictions of the algorithm match closely those of Monte Carlo based methods in the case of the Red Blood Cell metabolic network. Then we test the effect of gene knock-outs on the size of the solution space in the case of E. coli central metabolism. Finally we analyze the statistical properties of the average fluxes of the reactions in the E. coli metabolic network. Conclusion We propose a novel efficient distributed algorithmic strategy to estimate the size and shape of the affine space of a non full-dimensional convex polytope in high dimensions. The method is shown to obtain, quantitatively and qualitatively compatible results with the ones of standard algorithms (where this comparison is possible) being still efficient on the analysis of large biological systems, where exact deterministic methods experience an explosion in algorithmic time. The algorithm we propose can be considered as an alternative to Monte Carlo sampling methods.

LooplessFluxSampler: an efficient toolbox for sampling the loopless flux solution space of metabolic models

Avoid Internal Loops in Steady State Flux Space Sampling

Imposition of Different Optimizing Object with Non-Linear Constraints on Flux Sampling and Elimination of Free Futile Pathways

Geometric algorithms for sampling the flux space of metabolic networks

Flux sampling in genome-scale metabolic modeling of microbial communities

Estimating the size of the solution space of metabolic networks

FluxAnalyzer: exploring structure, pathways, and flux distributions in metabolic networks on interactive flux maps

A Novel Methodology to Estimate Metabolic Flux Distributions in Constraint-Based Models

Accelerating flux balance calculations in genome-scale metabolic models by localizing the application of loopless constraints

Elimination of thermodynamically infeasible loops in steady-state metabolic models

An Improved Multi-parametric Programming Algorithm for Flux Balance Analysis of Metabolic Networks

Improving Metabolic Flux Estimation Via Evolutionary Optimization For Convex Solution Space

Constraint-based inverse modeling of metabolic networks: a proof of concept

A robust and efficient method for estimating enzyme complex abundance and metabolic flux from expression data

Fluxer: a web application to compute, analyze and visualize genome-scale metabolic flux networks

Fluxomers: a new approach for 13C metabolic flux analysis

E-Flux2 and SPOT: Validated Methods for Inferring Intracellular Metabolic Flux Distributions from Transcriptomic Data

FluxExplorer： A General Platform for Modeling and Analyses of Metabolic Networks Based on Stoichiometry

Conserved moiety fluxomics

Counting and correcting thermodynamically infeasible flux cycles in genome-scale metabolic networks

Metabolite Biomarker Discovery for Metabolic Diseases by Flux Analysis