Henry Christopher S, Broadbelt Linda J, Hatzimanikatis Vassily
Department of Chemical and Biological Engineering, McCormick School of Engineering and Applied Sciences, Northwestern University, Evanston, Illinois, USA.
Biophys J. 2007 Mar 1;92(5):1792-805. doi: 10.1529/biophysj.106.093138. Epub 2006 Dec 15.
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.
一种新的代谢通量分析(MFA)形式,即基于热力学的代谢通量分析(TMFA)被引入,它能够在基因组规模上生成热力学可行的通量和代谢物活性谱。TMFA除了使用MFA中通常使用的质量平衡约束外,还涉及一组线性热力学约束的使用。TMFA产生的通量分布不包含任何热力学不可行的反应或途径,并且除了反应通量外,还提供有关反应自由能变化和代谢物活性范围的信息。TMFA被应用于研究由帕尔森及其同事开发的大肠杆菌基因组规模代谢模型中通量的热力学可行范围、反应的吉布斯自由能变化Δ(r)G'以及代谢物的活性。在基因组规模模型的TMFA中,代谢物活性和反应Δ(r)G'在最佳生长时能够达到广泛的值。反应二氢乳清酸酶被确定为大肠杆菌代谢中可能的热力学瓶颈,其Δ(r)G'被限制在接近零的水平,而在整个代谢过程中发现许多反应的Δ(r)G'无论代谢物浓度如何总是高度为负。正如之前所提出的,这些具有唯一负Δ(r)G'的反应可能是细胞调控的候选者,并且我们发现这些反应中有相当数量似乎是众多生物合成途径线性部分的第一步。还确定了ATP/ADP、NAD(P)/NAD(P)H和H(细胞外)(+)/H(细胞内)(+)浓度比的热力学可行范围,并且发现在每种情况下都包含实验观察到的值。此外,我们发现细胞中维持的NAD/NADH和NADP/NADPH比值分别接近最小可行比值和最大可行比值。