Designing overall stoichiometric conversions and intervening metabolic reactions
- 作者
- Anupam Chowdhury, Costas D. Maranas
- 杂志
- Scientific Reports（2015）
- DOI
- 10.1038/srep16009
- 所属领域
- 代谢工程/合成生物学
- 关键词
- 阅读原文
- 查看 PDF 原文
解决的核心问题
现有从头代谢途径设计工具通常只寻找连接底物与产物的线性路径，不考虑总体转化化学计量、辅因子/共底物/共产物组合，导致碳和能量效率无法在热力学极限附近优化。
研究策略
提出两步计算流程：第一步(optStoic)通过线性规划优化总体转化化学计量，同时满足元素、电荷和热力学可行性约束，并探索共底物/共产物组合；第二步(minRxn/minFlux)从反应数据库中提取符合该化学计量的最小反应集，可附加限制正ΔG反应数、异源反应数或特定宿主来源。
核心内容
该研究提出了一种代谢途径设计的两步计算流程，以突破传统线性路径设计的局限。第一步通过线性规划优化总体反应化学计量，同时满足元素、电荷和热力学可行性约束，并将共底物与共产物纳入设计变量；第二步从反应数据库中筛选符合该化学计量的最小反应集，并支持限制正ΔG反应数或异源反应数。以葡萄糖产乙酸为例，传统途径最大碳效率为67%，而设计的非氧化糖酵解途径可达100%碳效率。在甲醇与CO2共利用场景中，系统识别出496种热力学可行的C2+产物，其中乙酸盐合成网络通过甲醇氧化驱动Wood-Ljungdahl途径固定CO2；甲烷转化案例则显示，连四硫酸盐体系可氧化4.5 mol甲烷生成2.25 mol乙酸，而铁(III)体系仅能氧化0.25 mol，表明电子受体选择对底物利用率影响显著。该框架通过优化化学计量与共代谢物组合，能够在热力学极限附近最大化碳和能量效率，为设计非直观的高效代谢网络提供了系统化工具。
创新点
['将总体化学计量优化作为途径设计的第一步，直接最大化碳/能量/价格效率', '允许共底物和共产物作为设计变量，突破线性路径限制，识别循环网络', '在热力学约束下同时优化反应网络和化学计量系数，而非事后用FBA校正', '提出minFlux（LP）和minRxn（MILP）两种网络筛选策略，可控制正ΔG反应数及宿主多样性']
研究对象（13）
| 底盘细胞 | 代谢通路 | 基因 | 蛋白 | 功能 | 说明 |
|---|---|---|---|---|---|
| Escherichia coli | 非氧化糖酵解（NOG） | — | XPK（木酮糖-5-磷酸磷酸酮酶） | 催化 | 来自B. adolescentis，异源表达，用于无碳损失转化葡萄糖 |
| Escherichia coli | 非氧化糖酵解（NOG） | — | FPK（果糖-6-磷酸磷酸酮酶） | 催化 | 来自B. adolescentis，与XPK可互换，参与乙酰磷酸生成 |
| Escherichia coli | 戊糖磷酸途径 | — | RPI（核酮糖-5-磷酸异构酶） | 催化 | NOG设计中需正ΔG的反应，可通过替代途径避免 |
| Mycobacterium sp. | 乙二醇代谢 | — | DDT（二醇脱水酶） | 催化 | 将乙二醇脱水为乙醛，用于葡萄糖至乙酸的非磷酸酮酶途径 |
| Clostridium ljungdahlii / Moorella thermoacetica | Wood-Ljungdahl途径 + EMP糖酵解 | — | PFOR（丙酮酸：铁氧还蛋白氧化还原酶） | 催化 | 催化丙酮酸氧化脱羧并固定CO2，实现葡萄糖全碳转化 |
| methylotrophs | 甲醇氧化 / RuMP循环 | — | MEDH（甲醇脱氢酶） | 催化 | 氧化甲醇提供还原力，驱动CO2固定 |
| methylotrophs | 核酮糖单磷酸途径（RuMP） | — | HPS（己酮糖-6-磷酸合酶） | 催化 | 固定甲醛到Ru5P骨架 |
| methylotrophs | 核酮糖单磷酸途径（RuMP） | — | PHI（6-磷酸-3-己酮异构酶） | 催化 | 与HPS偶联，参与甲醛固定 |
| methanogenic archaea | 产甲烷途径逆转 / Wood-Ljungdahl | — | MTA（甲基转移酶） | 催化 | 参与甲烷活化生成甲基-CoM |
| methanogenic archaea | Wood-Ljungdahl途径 | — | ACS（乙酰-CoA合酶） | 催化 | 催化CO2与甲基合成乙酰-CoA |
| Clostridium sp. | 甲酸代谢 | — | FDH（甲酸脱氢酶） | 催化 | 固定CO2生成甲酸，用于乙酸盐合成设计 |
| Escherichia coli | 脂肪酸合成 / 丙二酰-CoA途径 | — | ACCOAC（乙酰-CoA羧化酶） | 催化 | 间苯三酚传统途径需额外ATP，在optStoic中被能量平衡排除 |
| Escherichia coli | 间苯三酚合成新途径 | — | 3HXCT（3-羟基-5-氧代己酰-CoA硫解酶） | 催化 | 催化乙酰-CoA与3-羟基丁酰-CoA缩合，生成间苯三酚前体 |
关键发现
- 葡萄糖至乙酸传统途径最大碳效率67%，设计的非氧化糖酵解（NOG）途径可达100%碳效率（图2a-f六种网络）
- 甲醇+CO2共利用时，以甲醇为1计，CO2化学计量系数随产物O:H比从1.25（丙二酸）降至0.07（2-甲基丁醛），共496种热力学可行的C2+产物，其中乙酸盐、3-羟基丁酸、2-酮异戊酸、间苯三酚的合成网络被详细构建（如乙酸盐设计通过甲醇氧化驱动Wood-Ljungdahl固定CO2）
- 甲烷转化为乙酸盐时，连四硫酸盐/硫化氢作为电子受体可氧化4.5 mol甲烷生成2.25 mol乙酸，而铁(III)/铁(II)仅氧化0.25 mol甲烷
- 硝酸盐/亚硝酸盐体系在允许额外共底物/共产物时甲烷利用率远高于亚硫酸氢盐/硫化氢体系（如异丁酸产物甲烷系数9.40 vs 5.62）。
研究结论
两步法（optStoic + minRxn/minFlux）能够设计出比现有工具更高效、非直观的代谢转化网络，通过优化总体化学计量和共代谢物组合，突破线性路径限制，在热力学可行范围内最大化碳/能量效率，并已在三个案例中验证其优越性。