Thermodynamics-based Metabolic Flux Analysis

Overview

Journal Biophys J

Publisher Cell Press

Specialty Biophysics

Date 2006 Dec 19

PMID 17172310

Citations 218

Authors

Christopher S Henry

Linda J Broadbelt

Vassily Hatzimanikatis

Affiliations

Soon will be listed here.

Abstract

A new form of metabolic flux analysis (MFA) called thermodynamics-based metabolic flux analysis (TMFA) is introduced with the capability of generating thermodynamically feasible flux and metabolite activity profiles on a genome scale. TMFA involves the use of a set of linear thermodynamic constraints in addition to the mass balance constraints typically used in MFA. TMFA produces flux distributions that do not contain any thermodynamically infeasible reactions or pathways, and it provides information about the free energy change of reactions and the range of metabolite activities in addition to reaction fluxes. TMFA is applied to study the thermodynamically feasible ranges for the fluxes and the Gibbs free energy change, Delta(r)G', of the reactions and the activities of the metabolites in the genome-scale metabolic model of Escherichia coli developed by Palsson and co-workers. In the TMFA of the genome scale model, the metabolite activities and reaction Delta(r)G' are able to achieve a wide range of values at optimal growth. The reaction dihydroorotase is identified as a possible thermodynamic bottleneck in E. coli metabolism with a Delta(r)G' constrained close to zero while numerous reactions are identified throughout metabolism for which Delta(r)G' is always highly negative regardless of metabolite concentrations. As it has been proposed previously, these reactions with exclusively negative Delta(r)G' might be candidates for cell regulation, and we find that a significant number of these reactions appear to be the first steps in the linear portion of numerous biosynthesis pathways. The thermodynamically feasible ranges for the concentration ratios ATP/ADP, NAD(P)/NAD(P)H, and H(extracellular)(+)/H(intracellular)(+) are also determined and found to encompass the values observed experimentally in every case. Further, we find that the NAD/NADH and NADP/NADPH ratios maintained in the cell are close to the minimum feasible ratio and maximum feasible ratio, respectively.

Citing Articles

Implementation of a genome-scale model for upgrading syngas fermentations.

Scott Jr W, Rockx S, Marien Q, Regueira A, Candry P, Ganigue R Comput Struct Biotechnol J. 2025; 27:649-660.

PMID: 40046987 PMC: 11880719. DOI: 10.1016/j.csbj.2025.01.013.

Generative machine learning produces kinetic models that accurately characterize intracellular metabolic states.

Choudhury S, Narayanan B, Moret M, Hatzimanikatis V, Miskovic L Nat Catal. 2024; 7(10):1086-1098.

PMID: 39463726 PMC: 11499278. DOI: 10.1038/s41929-024-01220-6.

Free-energy transduction mechanisms shape the flux space of metabolic networks.

Pfeuty B Biophys J. 2024; 123(20):3600-3611.

PMID: 39277793 PMC: 11494513. DOI: 10.1016/j.bpj.2024.09.012.

Quantitative principles of microbial metabolism shared across scales.

Sher D, Segre D, Follows M Nat Microbiol. 2024; 9(8):1940-1953.

PMID: 39107418 DOI: 10.1038/s41564-024-01764-0.

Current State, Challenges, and Opportunities in Genome-Scale Resource Allocation Models: A Mathematical Perspective.

Schroeder W, Suthers P, Willis T, Mooney E, Maranas C Metabolites. 2024; 14(7).

PMID: 39057688 PMC: 11278519. DOI: 10.3390/metabo14070365.

References

Albe K, Butler M, Wright B . Cellular concentrations of enzymes and their substrates. J Theor Biol. 1990; 143(2):163-95. DOI: 10.1016/s0022-5193(05)80266-8. View

Heinrich R, Rapoport T . A linear steady-state treatment of enzymatic chains. General properties, control and effector strength. Eur J Biochem. 1974; 42(1):89-95. DOI: 10.1111/j.1432-1033.1974.tb03318.x. View

Beard D, Liang S, Qian H . Energy balance for analysis of complex metabolic networks. Biophys J. 2002; 83(1):79-86. PMC: 1302128. DOI: 10.1016/S0006-3495(02)75150-3. View

Edwards J, Palsson B . Robustness analysis of the Escherichia coli metabolic network. Biotechnol Prog. 2000; 16(6):927-39. DOI: 10.1021/bp0000712. View

Edwards J, Ibarra R, Palsson B . In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data. Nat Biotechnol. 2001; 19(2):125-30. DOI: 10.1038/84379. View

Reed J, Vo T, Schilling C, Palsson B . An expanded genome-scale model of Escherichia coli K-12 (iJR904 GSM/GPR). Genome Biol. 2003; 4(9):R54. PMC: 193654. DOI: 10.1186/gb-2003-4-9-r54. View

Fell D, Sauro H . Metabolic control and its analysis. Additional relationships between elasticities and control coefficients. Eur J Biochem. 1985; 148(3):555-61. DOI: 10.1111/j.1432-1033.1985.tb08876.x. View

Varma A, Palsson B . Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110. Appl Environ Microbiol. 1994; 60(10):3724-31. PMC: 201879. DOI: 10.1128/aem.60.10.3724-3731.1994. View

Jin Y, Jeffries T . Stoichiometric network constraints on xylose metabolism by recombinant Saccharomyces cerevisiae. Metab Eng. 2004; 6(3):229-38. DOI: 10.1016/j.ymben.2003.11.006. View

10.

Andersen K, von Meyenburg K . Charges of nicotinamide adenine nucleotides and adenylate energy charge as regulatory parameters of the metabolism in Escherichia coli. J Biol Chem. 1977; 252(12):4151-6. View

11.

Mavrovouniotis M . Estimation of standard Gibbs energy changes of biotransformations. J Biol Chem. 1991; 266(22):14440-5. View

12.

Yang C, Hua Q, Baba T, Mori H, Shimizu K . Analysis of Escherichia coli anaplerotic metabolism and its regulation mechanisms from the metabolic responses to altered dilution rates and phosphoenolpyruvate carboxykinase knockout. Biotechnol Bioeng. 2003; 84(2):129-44. DOI: 10.1002/bit.10692. View

13.

Beard D, Qian H . Thermodynamic-based computational profiling of cellular regulatory control in hepatocyte metabolism. Am J Physiol Endocrinol Metab. 2004; 288(3):E633-44. DOI: 10.1152/ajpendo.00239.2004. View

14.

Nielsen J . Metabolic control analysis of biochemical pathways based on a thermokinetic description of reaction rates. Biochem J. 1997; 321 ( Pt 1):133-8. PMC: 1218046. DOI: 10.1042/bj3210133. View

15.

Mori M, TATIBANA M . Multi-enzyme complex of glutamine-dependent carbamoyl-phosphate synthetase with aspartate carbamoyltransferase and dihydroorotase from rat ascites-hepatoma cells. Purification, molecular properties and limited proteolysis. Eur J Biochem. 1978; 86(2):381-8. DOI: 10.1111/j.1432-1033.1978.tb12320.x. View

16.

Beard D, Babson E, Curtis E, Qian H . Thermodynamic constraints for biochemical networks. J Theor Biol. 2004; 228(3):327-33. DOI: 10.1016/j.jtbi.2004.01.008. View

17.

Peng L, Arauzo-Bravo M, Shimizu K . Metabolic flux analysis for a ppc mutant Escherichia coli based on 13C-labelling experiments together with enzyme activity assays and intracellular metabolite measurements. FEMS Microbiol Lett. 2004; 235(1):17-23. DOI: 10.1016/j.femsle.2004.04.003. View

18.

Mavrovouniotis M . Group contributions for estimating standard gibbs energies of formation of biochemical compounds in aqueous solution. Biotechnol Bioeng. 1990; 36(10):1070-82. DOI: 10.1002/bit.260361013. View

19.

ALBERTY R . Calculation of standard transformed formation properties of biochemical reactants and standard apparent reduction potentials of half reactions. Arch Biochem Biophys. 1998; 358(1):25-39. DOI: 10.1006/abbi.1998.0831. View

20.

Goldberg R, Tewari Y, Bhat T . Thermodynamics of enzyme-catalyzed reactions--a database for quantitative biochemistry. Bioinformatics. 2004; 20(16):2874-7. DOI: 10.1093/bioinformatics/bth314. View