[论文解读] Concordant Chemical Reaction Networks

Guy Shinar, Martin Feinberg|arXiv (Cornell University)|Sep 13, 2011

Gene Regulatory Network Analysis参考文献 26被引用 4

一句话总结

本文引入了同调化学反应网络——一种具有结构特性、可在弱单调动力学下确保正平衡态的单射性和唯一性的网络。论文证明，同调性可保证不存在多个正平衡态或化学计量相容的平衡态，从而排除双稳态行为，并建立了持久性和稳定性特性，即使在无降解反应的情况下也成立。

ABSTRACT

We describe a large class of chemical reaction networks, those endowed with a subtle structural property called concordance. We show that the class of concordant networks coincides precisely with the class of networks which, when taken with any weakly monotonic kinetics, invariably give rise to kinetic systems that are injective --- a quality that, among other things, precludes the possibility of switch-like transitions between distinct positive steady states. We also provide persistence characteristics of concordant networks, instability implications of discordance, and consequences of stronger variants of concordance. Some of our results are in the spirit of recent ones by Banaji and Craciun, but here we do not require that every species suffer a degradation reaction. This is especially important in studying biochemical networks, for which it is rare to have all species degrade.

研究动机与目标

• 识别一类具有固有受限动力学行为的化学反应网络，确保在不同正稳态之间不会发生类似开关的跃迁。
• 建立同调性——一种网络的结构属性——可保证在任意弱单调动力学下实现单射性和平衡态的唯一性，前提是动力学满足弱单调性。
• 将结果扩展至不含降解反应的网络，这对真实生化建模至关重要，因为在这些模型中并非所有物种都会降解。
• 表征非同调网络的持久性和不稳定性特性，并与同调网络进行对比。
• 提供一种基于化学计量子空间和物种-反应图分析的框架，用于分析复杂网络中的稳定性和唯一性。

提出的方法

• 将同调性定义为反应网络的一种结构属性，与动力学无关，通过反应化学计量向量和网络矩阵核中非负向量的条件来定义。
• 采用弱单调动力学，其中反应速率随反应物浓度增加而上升，以模拟真实的生化速率定律。
• 应用线性代数技术分析物种生成速率函数的雅可比矩阵，特别是其核和特征值。
• 通过速率常数和反应级数的参数化方法，显式构造具有正平衡态的动力学系统，确保所需的平衡态和稳定性特性。
• 利用网络的完全开放扩展来分析稳定性和持久性，特别是在边界组成情况下的表现。
• 使用物种-反应图及软件工具（例如参考文献[24]中的工具）通过观察或计算验证同调性。

实验结果

研究问题

• RQ1 化学反应网络的何种结构特性可确保在任意弱单调动力学下，不会共存两个不同的正平衡态？
• RQ2 即使在缺乏降解反应的复杂网络中，同调性如何防止双稳态行为？
• RQ3 非同调网络与同调网络相比，其持久性和不稳定性特性有何不同？
• RQ4 在光滑动力学下，弱可逆且保守的同调网络中的唯一平衡态是否可证明为全局稳定？
• RQ5 如何高效检查同调性？物种-反应图和软件工具在此过程中起到何种作用？

主要发现

• 同调网络正是那些所有弱单调动力学均产生单射动力学系统的网络，从而排除了多个正平衡态的存在。
• 若一个同调网络是弱可逆且保守的，则在每个化学计量相容类中恰好存在一个正平衡态，且边界组成是排斥的。
• 对于光滑动力学和同调完全开放扩展，唯一正平衡态的所有特征值均为负实数，意味着局部渐近稳定。
• 非同调网络可能具有退化或不稳定的平衡态；具体而言，若完全开放扩展是非同调的，则正平衡态可能具有正实特征值，意味着不稳定。
• 缺乏降解反应并非理论的障碍，因为同调性定义时并不要求存在此类反应。
• 存在一个非零向量 σ 满足特定核条件和符号条件，可表征非同调性，从而通过显式动力学参数化方法构造不稳定平衡态。