[论文解读] Finding weakly reversible realizations of chemical reaction networks using optimization

Gábor Szederkényi, Katalin M. Hangos+1 more|arXiv (Cornell University)|Mar 24, 2011

Gene Regulatory Network Analysis参考文献 23被引用 10

一句话总结

本文提出一种基于优化的算法，通过利用密集实现和混合整数线性规划（MILP），计算具有质量作用动力学的化学反应网络（CRNs）的动态等价弱可逆实现。该方法识别与给定CRN相容的最大弱可逆结构，或确定不存在此类实现，利用图论性质和线性规划技术。

ABSTRACT

An algorithm is given in this paper for the computation of dynamically equivalent weakly reversible realizations with the maximal number of reactions, for chemical reaction networks (CRNs) with mass action kinetics. The original problem statement can be traced back at least 30 years ago. The algorithm uses standard linear and mixed integer linear programming, and it is based on elementary graph theory and important former results on the dense realizations of CRNs. The proposed method is also capable of determining if no dynamically equivalent weakly reversible structure exists for a given reaction network with a previously fixed complex set.

研究动机与目标

• 开发一种系统化方法，用于计算具有质量作用动力学的CRNs的动态等价弱可逆实现。
• 确定给定CRN在固定复合物集合下是否可实现任何弱可逆实现。
• 计算与给定CRN动态等价的、反应数最多的弱可逆结构。
• 提供一种计算框架，整合图论与优化技术，用于CRNs的结构分析。
• 解决CRN理论中关于弱可逆性的开放挑战，如20世纪80年代和2000年代基础文献中所强调的。

提出的方法

• 该算法首先使用混合整数线性规划（MILP）计算给定CRN的密集实现，其中包含在固定复合物集合下所有可能的动态等价反应。
• 然后通过分析反应图中的强连通分量，识别并移除违反弱可逆性的有向边。
• 该方法利用图论分解技术隔离非平凡的强连通分量，并迭代删除阻碍弱可逆性的边。
• 关键步骤包括移除连接不同强连通分量的边，随后进一步修剪剩余子图中破坏弱可逆性的边。
• 该算法采用线性规划（LP）和MILP以确保可行性与最优性，最终输出为反应数最多的弱可逆CRN。
• 该过程以分步方式进行，中间可视化图示（如图6–9）展示了删除过程及最终结构。

实验结果

研究问题

• RQ1 对于任意给定的具有质量作用动力学的化学反应网络，能否计算出其动态等价的弱可逆实现？
• RQ2 与给定CRN动态等价的、反应数最多的弱可逆结构是什么？
• RQ3 能否通过算法确定在固定复合物集合下，给定CRN不存在任何弱可逆实现？
• RQ4 优化技术如何系统性地强制实现弱可逆性，同时保持动态等价性？
• RQ5 使用标准优化求解器构建此类实现的计算复杂度和可行性如何？

主要发现

• 该算法成功计算出一个包含80个反应的弱可逆实现，这是在给定复合物集合下可能的最大值。
• 该网络的密集实现包含80个反应，最终的弱可逆结构通过迭代修剪移除53条边和12个额外反应获得。
• 所得弱可逆网络的缺陷数为零，表明其具有有利的动力学性质，如持久性和稳定性。
• 当必要结构条件被违反时，该算法能正确识别出不存在弱可逆实现，从而提供可行性检查。
• 该算法在示例网络上的总运行时间为80.5秒，证明了其实际计算可行性。
• 该方法具有通用性，可扩展至并行实现，尤其适用于计算密集的MILP阶段。