[论文解读] Dynamical Equivalence and Linear Conjugacy of Chemical Reaction Networks: New Results and Methods

Matthew D. Johnston, David Siegel+1 more|arXiv (Cornell University)|Nov 6, 2011

Computational Drug Discovery Methods参考文献 29被引用 13

一句话总结

本文提出混合整数线性规划（MILP）方法，用于在质量作用动力学下计算密集、稀疏、可逆、详细平衡及复合平衡的线性共轭反应网络。证明了密集线性共轭网络由复杂数集唯一确定，且平衡选择不影响复合平衡或详细平衡实现的存在性，从而可在无需预先知道速率常数的情况下，基于结构分析反应机理。

ABSTRACT

In the first part of this paper, we propose new optimization-based methods for the computation of preferred (dense, sparse, reversible, detailed and complex balanced) linearly conjugate reaction network structures with mass action dynamics. The developed methods are extensions of previously published results on dynamically equivalent reaction networks and are based on mixed-integer linear programming. As related theoretical contributions we show that (i) dense linearly conjugate networks define a unique super-structure for any positive diagonal state transformation if the set of chemical complexes is given, and (ii) the existence of linearly conjugate detailed balanced and complex balanced networks do not depend on the selection of equilibrium points. In the second part of the paper it is shown that determining dynamically equivalent realizations to a network that is structurally fixed but parametrically not can also be written and solved as a mixed-integer linear programming problem. Several examples illustrate the presented computation methods.

研究动机与目标

• 开发基于优化的方法，用于计算具有特定动力学性质（如可逆性、详细平衡、复合平衡）的首选线性共轭反应网络结构。
• 扩展现有的MILP框架，以处理具有最小和最大复杂数、反应数及结构约束下的网络。
• 解决当网络结构固定但速率常数未知时的动力学等价性问题，实现仅基于结构的反应机理分析。
• 证明密集线性共轭网络的唯一性以及复合/详细平衡存在性在不同平衡点下的不变性。
• 通过实现无参数的结构分析，拓展反应网络理论的应用范围，以支持动力学等价性和共轭性的研究。

提出的方法

• 将线性共轭网络的搜索问题表述为混合整数线性规划（MILP）问题，利用稀疏性和结构约束。
• 使用化学计量矩阵和复杂数矩阵表示来编码网络结构和动力学，将速率常数作为优化中的变量。
• 通过线性等式和不等式条件施加弱可逆性、可逆性、详细平衡和复合平衡的约束。
• 提出基于变换的方法，通过非平凡的线性共轭在不同网络间映射动力学，保持定性行为不变。
• 利用正对角状态变换关联动力学等价性和共轭性，实现与平衡点无关的分析。
• 将MILP框架扩展至速率常数未固定的情形，实现独立于特定动力学参数的结构分析。

实验结果

研究问题

• RQ1 我们能否通过优化方法计算出具有指定动力学性质（如可逆性、详细平衡）的最稀疏和最密集的线性共轭网络？
• RQ2 最密集的线性共轭网络的结构是否由复杂数集和一个正对角变换唯一确定？
• RQ3 线性共轭复合平衡或详细平衡网络的存在性是否依赖于平衡点的选择？
• RQ4 当仅知道网络结构而未知速率常数时，能否确定并计算动力学等价性？
• RQ5 如何系统性地生成反应网络的替代实现形式，使其在质量作用动力学下保持动力学特性但反应结构不同？

主要发现

• 密集线性共轭网络由复杂数集和一个正对角状态变换唯一确定，与速率常数无关。
• 线性共轭详细平衡和复合平衡网络的存在性独立于平衡点的选择，从而可从任意正平衡点灵活构建。
• 任何线性共轭网络的无权有向反应图均为最密集此类网络的子图，确保了结构层次关系。
• 在弱可逆性、可逆性以及复合/详细平衡约束下，寻找稀疏和密集实现的问题可表述为MILP，从而实现高效计算。
• 当速率常数未固定时，计算具有最小或最大复杂数的动态等价实现仍为MILP问题，支持仅基于结构的分析。
• 通过重新分配速率常数，可构建与动态等价实现具有相同结构的复合平衡网络，表明网络设计中具有参数灵活性。