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 the 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.
References
[1]
Price, N.; Famili, I.; Beard, D.; Palsson, B. Extreme pathways and kirchhoff's second law. Biophys. J. 2002, 83, 2879–2882, doi:10.1016/S0006-3495(02)75297-1. 12425318
[2]
Soh, K.; Hatzimanikatis, V. Network thermodynamics in the post-genomic era. Curr. Opin. Microbiol. 2010, 13, 350–357, doi:10.1016/j.mib.2010.03.001. 20378394
[3]
Beard, D.; Babson, E.; Curtis, E.; Qian, H. Thermodynamic constraints for biochemical networks. J. Theor. Biol. 2004, 228, 327–333, doi:10.1016/j.jtbi.2004.01.008. 15135031
[4]
Hoppe, A.; Hoffmann, S.; Holzhutter, H. Including metabolite concentrations into flux balance analysis: Thermodynamic realizability as a constraint on flux distributions in metabolic networks. BMC Syst. Biol. 2007, 1, 23:1–23:12.
[5]
Qian, H.; Beard, D. Thermodynamics of stoichiometric biochemical networks in living systems far from equilibrium. Biophys. Chem. 2005, 114, 213–220, doi:10.1016/j.bpc.2004.12.001. 15829355
[6]
Beard, D.; Liang, S.; Qian, H. Energy balance for analysis of complex metabolic networks. Biophys. J. 2002, 83, 79–86, doi:10.1016/S0006-3495(02)75150-3. 12080101
[7]
Palsson, B.O. Systems Biology: Properties of Reconstructed Networks; Cambridge University Press: Cambridge, NY, USA, 2006.
[8]
Bowden, A.C. Fundamentals of Enzyme Kinetics; Wiley-Blackwell: Weinheim, Germany, 2013.
[9]
Ge, H.; Qian, M.; Qian, H. Stochastic theory of nonequilibrium steady states. Part II: Applications in chemical biophysics. Phys. Rep. 2012, 510, 87–118, doi:10.1016/j.physrep.2011.09.001.
[10]
Frey, E.; Kroy, K. Brownian motion: A paradigm of soft matter and biological physics. Ann. Phys. 2005, 14, 20–50, doi:10.1002/andp.200410132.
[11]
Beg, Q.; Vazquez, A.; Ernst, J.; de Menezes, M.; Bar-Joseph, Z.; Barabási, A.L.; Oltvai, Z.N. Intracellular crowding defines the mode and sequence of substrate uptake by Escherichia coli and constrains its metabolic activity. Proc. Natl. Acad. Sci. USA 2007, 104, 12663–12668, doi:10.1073/pnas.0609845104. 17652176
[12]
De Martino, A.; Marinari, E. The solution space of metabolic networks: Producibility, robustness and fluctuations. J. Phys. Conf. Ser. 2010, 233, 012019, doi:10.1088/1742-6596/233/1/012019.
[13]
Schrijver, A. Theory of Linear and Integer Programming; Wiley: New York, NY, USA, 1986.
[14]
Orth, J.; Thiele, I.; Palsson, B.O. What is flux balance analysis? Nat. Biotechnol. 2010, 28, 245–248, doi:10.1038/nbt.1614. 20212490
[15]
Segrè, D.; Vitkup, D.; Church, G. Analysis of optimality in natural and perturbed metabolic networks. Proc. Natl. Acad. Sci. USA 2002, 99, 15112–15117, doi:10.1073/pnas.232349399. 12415116
[16]
Jankowski, M.; Henry, C.; Broadbelt, L.; Hatzimanikatis, V. Group contribution method for thermodynamic analysis of complex metabolic networks. Biophys. J. 2008, 95, 1487–1499, doi:10.1529/biophysj.107.124784. 18645197
[17]
Fleming, R.; Thiele, I.; Nasheuer, H. Quantitative assignment of reaction directionality in constraint-based models of metabolism: Application to Escherichia coli. Biophys. Chem. 2009, 145, 47–56, doi:10.1016/j.bpc.2009.08.007. 19783351
[18]
Kummel, A.; Panke, S.; Heinemann, M. Systematic assignment of thermodynamic constraints in metabolic network models. BMC Bioinforma. 2006, 7, 512:1–512:12.
[19]
Alberty, R.A. Thermodynamics of Biochemical Reactions; Wiley: Hoboken, NJ, USA, 2003.
[20]
Schellenberger, J.; Lewis, N.; Palsson, B.O. Elimination of thermodynamically infeasible loops in steady-state metabolic models. Biophys. J. 2011, 100, 544–553, doi:10.1016/j.bpj.2010.12.3707. 21281568
[21]
Henry, C.; Broadbelt, L.; Hatzimanikatis, V. Thermodynamics-based metabolic flux analysis. Biophys. J. 2007, 92, 1792–1805, doi:10.1529/biophysj.106.093138. 17172310
[22]
Müller, A.; Brockmayr, A. Fast thermodynamically constrained flux variability analysis. Bioinformatics 2013, 29, 903–909, doi:10.1093/bioinformatics/btt059. 23390138
[23]
Beard, D.A.; Qian, H. Chemical Biophysics; Cambridge University Press: Cambridge, NY, USA, 2008.
[24]
De Martino, D.; Figliuzzi, M.; de Martino, A.; Marinari, E. A scalable algorithm to explore the gibbs energy landscape of genome-scale metabolic networks. PLoS Comp. Biol. 2012, 8, e1002562, doi:10.1371/journal.pcbi.1002562.
[25]
De Martino, D. Thermodynamics of biochemical networks and duality theorems. Phys. Rev. E 2013, 87, 053108, doi:10.1103/PhysRevE.87.053108.
[26]
Johnson, D.B. Finding all the elemtary circuits of a directed graph. SIAM J. Comput. 1975, 4, 77–84, doi:10.1137/0204007.
[27]
Wright, J.; Wagner, A. Exhaustive identification of steady state cycles in large stoichiometric networks. BMC Syst. Biol. 2008, 2, 61:1–61:9.
[28]
Wiback, S.; Famili, I.; Greenberg, H.; Palsson, B.O. Monte Carlo sampling can be used to determine the size and shape of the steady-state flux space. J. Theor. Biol. 2004, 228, 437–447, doi:10.1016/j.jtbi.2004.02.006. 15178193
[29]
Price, N.; Schellenberger, J.; Palsson, B.O. Uniform sampling of steady-state flux spaces: Means to design experiments and to interpret enzymopathies. Biophys. J. 2004, 87, 2172–2186, doi:10.1529/biophysj.104.043000. 15454420
[30]
Mezard, M.; Montanari, A. Information, Physics, and Computation; Oxford University Press: Oxford, NY, USA, 2009.
[31]
Feist, A.; Henry, C.; Reed, J.; Krummenacker, M.; Joyce, A.; Karp, P.; Broadbelt, L.; Hatzimanikatis, V.; Palsson, B.O. A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol. Syst. Biol. 2007, 3, 121:1–121:18.
[32]
Thiele, I.; Swainston, N.; Fleming, R.M.T.; Hoppe, A.; Sahoo, S.; Aurich, M.K.; Haraldsdottir, H.; Mo, M.L.; Rolfsson, O.; Stobbe, M.D.; et al. A community-driven global reconstruction of human metabolism. Nat. Biotechnol. 2013, 31, 419–425, doi:10.1038/nbt.2488. 23455439
[33]
Schellenberger, J.; Que, R.; Fleming, R.M.T.; Thiele, I.; Orth, J.D.; Feist, A.M.; Zielinski, D.C.; Bordbar, A.; Lewis, N.E.; Rahmanian, S.; et al. Quantitative prediction of cellular metabolism with constraint-based models: The COBRA Toolbox v2.0. Nat. Protoc. 2011, 6, 1290–1307, doi:10.1038/nprot.2011.308. 21886097
[34]
Shlomi, T.; Cabili, M.N.; Herrg？rd, M.J.; Palsson, B.O.; Ruppin, E. Network-based prediction of human tissue-specific metabolism. Nat. Biotechnol. 2008, 26, 1003–1010, doi:10.1038/nbt.1487. 18711341
[35]
Krauth, W.; Mezard, M. Learning algorithms with optimal stability in neural networks. J. Phys. A 1987, 20, L745–L752, doi:10.1088/0305-4470/20/11/013.
[36]
Binder, K.; Heermann, D.W. Monte Carlo Simulation in Statistical Physics; Springer: Heidelberg, Germany, 2002.
[37]
De Martino, A.; de Martino, D.; Mulet, R.; Uguzzoni, G. Reaction networks as systems for resource allocation: A variational principle for their non-equilibrium steady states. PLoS One 2012, 7, e39849, doi:10.1371/journal.pone.0039849. 22815715
[38]
Schilling, C.H.; Letscher, D.; Palsson, B.O. Theory for the systemic definition of metabolic pathways and their use in interpreting metabolic function from a pathway-oriented perspective. J. Theor. Biol. 2000, 203, 229–248, doi:10.1006/jtbi.2000.1073. 10716907
[39]
Mahadevan, R.; Schilling, C. The effects of alternate optimal solutions in constraint-based genome-scale metabolic models. Metab. Eng. 2003, 5, 264–276, doi:10.1016/j.ymben.2003.09.002. 14642354
