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

Abstract

Thermodynamics constrains the flow of matter in a reaction network to occur through routes along which the Gibbs energy decreases, implying that viable steady-state flux patterns should be void of closed reaction cycles. Identifying and removing cycles in large reaction networks can unfortunately be a highly challenging task from a computational viewpoint. We propose here a method that accomplishes it by combining a relaxation algorithm and a Monte Carlo procedure to detect loops, with ad hoc rules (discussed in detail) to eliminate them. As test cases, we tackle (a) the problem of identifying infeasible cycles in the E. coli metabolic network and (b) the problem of correcting thermodynamic infeasibilities in the Flux-Balance-Analysis solutions for 15 human cell-type specific metabolic networks. Results for (a) are compared with previous analyses of the same issue, while results for (b) are weighed against alternative methods to retrieve thermodynamically viable flux patterns based on minimizing specific global quantities. Our method on one hand outperforms previous techniques and, on the other, corrects loopy solutions to Flux Balance Analysis. As a byproduct, it also turns out to be able to reveal possible inconsistencies in model reconstructions.

10 pages; see http://chimera.roma1.infn.it/SYSBIO/ for supporting files

References in corpus (4)

• A scalable algorithm to explore the Gibbs energy landscape of genome-scale metabolic networks
• On the Clausius formulation of the second law in stationary chemical networks through the theorems of the alternative
• The solution space of metabolic networks: producibility, robustness and fluctuations
• Reaction networks as systems for resource allocation: A variational principle for their non-equilibrium steady states

Cited by in corpus (8)

• Irreversible thermodynamics of open chemical networks I: Emergent cycles and broken conservation laws
• Quantitative constraint-based computational model of tumor-to-stroma coupling via lactate shuttle
• Relationship between fitness and heterogeneity in exponentially growing microbial populations
• Inference for dynamics of continuous variables: the Extended Plefka Expansion with hidden nodes
• Microenvironmental cooperation promotes early spread and bistability of a Warburg-like phenotype
• The free lunch of a scale-free metabolism
• The dual of the space of interactions in neural network models
• Phase transitions in integer linear problems