[论文解读] Existence of Positive Steady States for Mass Conserving and Mass-Action Chemical Reaction Networks with a Single Terminal-Linkage Class

Santiago Akle, Onkar Dalal+4 more|arXiv (Cornell University)|May 12, 2011

Microbial Metabolic Engineering and Bioproduction参考文献 11被引用 4

一句话总结

该论文证明了在质量作用动力学下，具有单一终态链接类的质量守恒化学反应网络，无论网络缺陷度或速率常数如何，始终存在正稳态。作者提出了一种凸优化公式和不动点算法，该算法在弱可逆同质系统中收敛，可通过迭代精炼通量和浓度向量可靠计算这些稳态。

ABSTRACT

We establish that mass conserving single terminal-linkage networks of chemical reactions admit positive steady states regardless of network deficiency and the choice of reaction rate constants. This result holds for closed systems without material exchange across the boundary, as well as for open systems with material exchange at rates that satisfy a simple sufficient and necessary condition. Our proof uses a fixed point of a novel convex optimization formulation to find the steady state behavior of chemical reaction networks that satisfy the law of mass-action kinetics. A fixed point iteration can be used to compute these steady states, and we show that it converges for weakly reversible homogeneous systems. We report the results of our algorithm on numerical experiments.

研究动机与目标

• 建立质量守恒化学反应网络中正稳态存在性的必要与充分条件。
• 开发一种计算上可行的方法，以独立于动力学速率常数的方式识别正稳态。
• 将闭系统的结果扩展至具有外部物质交换的开放系统，仅需对交换速率施加一个简单条件。
• 设计并验证一种用于计算弱可逆、同质网络中稳态的不动点算法。
• 研究该算法在不同网络拓扑和规模下的收敛行为与可扩展性。

提出的方法

• 利用新颖的凸优化框架，将稳态问题表述为对通量和浓度向量的不动点迭代。
• 以对数形式定义质量作用条件：$ Y^T \log(c^*\,\!)=\log(v^*) $，将物种浓度与反应通量联系起来。
• 利用通量平衡条件 $ YA_k v^* = b $ 和质量作用条件 $ \psi(c^*) = v^* $ 来表征稳态。
• 通过埃拉德-雷尼图模型生成随机网络，同时施加约束以确保质量守恒和所需的终态链接类结构。
• 对化学计量矩阵 $ Y $ 进行归一化，使得 $ Y^T \mathbf{1} = \mathbf{1} $，通过 $ e^T Y A_k = 0 $ 确保质量守恒。
• 实现算法1，迭代更新 $ v_k $ 和 $ c_k $，最小化 $ \|YA_k v_k\|_\infty $ 中的不可行性，直至收敛。

实验结果

研究问题

• RQ1 质量守恒且具有单一终态链接类的化学反应网络在何种条件下存在正稳态？
• RQ2 不动点算法能否可靠计算弱可逆、同质网络中的正稳态？
• RQ3 该算法的收敛行为如何随网络规模和拓扑结构变化？
• RQ4 正稳态的存在性是否可推广至具有多个终态链接类的网络？
• RQ5 外部物质交换速率 $ b $ 在决定正稳态的存在性与可计算性方面起什么作用？

主要发现

• 所有质量守恒且仅具有单一终态链接类的网络，无论缺陷度或速率常数如何，均存在正稳态。
• 对于弱可逆、同质系统，该不动点算法单调收敛，且每次迭代中通量平衡条件的不可行性均减小。
• 对于包含50个物种和500个复杂物的网络，当网络规模扩大五十倍时，迭代次数增加不足10%。
• 平均迭代次数在不同网络规模下保持稳定，表明该算法对大规模系统具有稳健的可扩展性。
• 数值实验表明，即使对于具有多个终态链接类的网络，正稳态很可能也存在，尽管理论证明尚待完成。
• 凸优化公式可直接计算稳态，无需事先知晓速率常数或网络缺陷度。