[论文解读] Striking a balance with Recon 2.1

Kieran Smallbone|arXiv (Cornell University)|Nov 22, 2013

Microbial Metabolic Engineering and Bioproduction参考文献 18被引用 6

一句话总结

本文提出了Recon 2.1，即人类代谢网络重建Recon 2的碳平衡版本，通过用特定分子形式替换通用代谢物以消除元素不平衡。该方法解决了先前允许虚假碳生成的不平衡反应，从而实现可靠的约束基建模；最终模型在元素上完全平衡，反应数扩展至71,159个，显著提升了系统生物学应用中的预测准确性。

ABSTRACT

Recon 2 is a highly curated reconstruction of the human metabolic network. Whilst the network is state of the art, it has shortcomings, including the presence of unbalanced reactions involving generic metabolites. By replacing these generic molecules with each of their specific instances, we can ensure full elemental balancing, in turn allowing constraint-based analyses to be performed. The resultant model, called Recon 2.1, is an order of magnitude larger than the original.

研究动机与目标

• 解决由通用代谢物引起的Recon 2中元素不平衡问题，特别是在脂肪酸反应中导致的虚假碳生成。
• 通过为通用R基团分配特定碳原子数并调整化学计量系数，确保所有反应中的碳平衡。
• 通过将所有通用代谢物替换为其特定对应物，构建一个完全元素平衡的模型，提升约束基分析的模型可靠性。
• 定义一种最小培养基，以总碳摄取作为生长限制因素，替代仅以葡萄糖为限制条件，以支持在脂肪酸上生长。
• 在模型完整性和计算复杂性之间取得平衡，仅在扩展中包含最常见的脂肪酸。

提出的方法

• 通过用特定脂肪酰辅酶A形式（如十四酰辅酶A、棕榈酰辅酶A）替换通用代谢物（如R-CoA），消除元素不平衡。
• 为通用R基团分配固定碳原子数（16），并通过调整化学计量系数实现反应的碳平衡。
• 将交换反应拆分为输入和输出组分，并引入“碳元反应物”以追踪总碳摄取量。
• 将总碳摄取限制为六个通量单位，以模拟生长限制条件，替代仅以葡萄糖为限制的条件。
• 通过为每种可能的脂肪酸组合创建独立反应，扩展所有涉及通用代谢物的反应，仅排除心磷脂合成酶和生物量反应。
• 对剩余在心磷脂和生物量中的通用物质使用“是”类型反应，以在不完全扩展的情况下维持模型完整性。

实验结果

研究问题

• RQ1 Recon 2中，特别是在脂肪酸代谢中，是否可以解决元素不平衡问题，以实现可靠的约束基建模？
• RQ2 将通用代谢物替换为特定分子形式对模型规模和计算可行性有何影响？
• RQ3 如何重新定义最小培养基，以支持在葡萄糖和脂肪酸上生长，同时仅保持单一生长限制营养素？
• RQ4 在将通用代谢物扩展为其特定形式时，模型完整性和计算复杂性之间的权衡如何？
• RQ5 移除通用代谢物在多大程度上提升了人类代谢网络重建的预测准确性？

主要发现

• Recon 2.1通过用特定脂肪酸形式替换通用R基团并调整化学计量系数，实现了完全的碳平衡，消除了虚假碳生成。
• 模型规模从Recon 2的7,440个反应增加至Recon 2.1的8,089个反应，进一步扩展至Recon 2.1x的71,159个反应，因元素平衡实现约10倍的扩展。
• 所有元素均实现平衡，Recon 2.1x中未发现不平衡或未知反应（0个），而Recon 2中分别为611个不平衡和1,325个未知反应。
• 独特代谢物数量增加了三倍，主要由于脂质分辨率提高，但仍远低于HMDB中定义的40,000种。
• 生物量反应和心磷脂合成酶反应仍部分保持通用形式，使用“是”类型反应以在不完全扩展的情况下保持模型功能。
• 最终模型Recon 2.1x完全实现元素平衡，可通过BioModels获取，访问编号为MODEL1311110001，支持高级系统生物学应用。