von Bertalanffy 1.0: a COBRA toolbox extension to thermodynamically constrain metabolic models

MOTIVATION In flux balance analysis of genome scale stoichiometric models of metabolism, the principal constraints are uptake or secretion rates, the steady state mass conservation assumption and reaction directionality. Here, we introduce an algorithmic pipeline for quantitative assignment of reaction directionality in multi-compartmental genome scale models based on an application of the second law of thermodynamics to each reaction. Given experimental or computationally estimated standard metabolite species Gibbs energy and metabolite concentrations, the algorithms bounds reaction Gibbs energy, which is transformed to in vivo pH, temperature, ionic strength and electrical potential. RESULTS This cross-platform MATLAB extension to the COnstraint-Based Reconstruction and Analysis (COBRA) toolbox is computationally efficient, extensively documented and open source. AVAILABILITY http://opencobra.sourceforge.net.

[1]  Matthew D. Jankowski,et al.  Group contribution method for thermodynamic analysis of complex metabolic networks. , 2008, Biophysical journal.

[2]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.

[3]  R. Alberty.  Biochemical thermodynamics: applications of Mathematica. , 2006, Methods of biochemical analysis.

[4]  M. Meehan.  General System Theory: Foundations, Development, Applications , 1969 .

[5]  J. Lopreato,et al.  General system theory : foundations, development, applications , 1970 .

[6]  Michael Hucka,et al.  LibSBML: an API Library for SBML , 2008, Bioinform..

[7]  Ronan M. T. Fleming,et al.  Quantitative assignment of reaction directionality in constraint-based models of metabolism: application to Escherichia coli. , 2009, Biophysical chemistry.

[8]  C. Mathews.  Thermodynamics of biochemical reactions , 2004 .

[9]  Adam M. Feist,et al.  A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information , 2007, Molecular systems biology.

[10]  B. Palsson,et al.  A protocol for generating a high-quality genome-scale metabolic reconstruction , 2010 .

[11]  A. Goldberg.  General System Theory: Foundations, Development, Applications. , 1969 .

[12]  Ronan M. T. Fleming,et al.  Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0 , 2007, Nature Protocols.