化学反应网络中的Cramér-Rao界与绝对灵敏度

Cramer-Rao bound and absolute sensitivity in chemical reaction networks

摘要 Abstract

化学反应网络(CRN)是一类重要的模型,用于理解生物功能,如细胞信息处理、代谢通路的稳健性和控制、昼夜节律等。然而,任何描述特定功能的CRN都不是孤立存在的,而是更大网络的一部分,因此始终受到外部变化的影响。在[Shinar, Alon, 和 Feinberg. "化学反应网络的敏感性和鲁棒性."《SIAM应用数学杂志》(2009): 977-998.]中,研究了CRN对线性守恒量变化的响应,即灵敏度,并提出了如何构建绝对的、即与基无关的灵敏度的问题。本文通过应用信息几何方法,提供了这样的构造。其想法是追踪特定化学物质浓度的变化如何传播到稳态下所有浓度的变化,这被编码在绝对灵敏度矩阵中。作为主要的技术工具,证明了一个针对CRN的多元Cramér-Rao界,该界基于准热力学稳态与指数族概率分布之间的类比。这导致了准热力学CRN的绝对灵敏度矩阵的一个线性代数特征刻画。作为一个例子,IDHKP-IDH乙醛酸旁路调控系统的中心模块通过浓度空间的大规模随机抽样进行了分析,并从实验上确认了IDH酶的稳健性发现,揭示了分布层面上隐藏的对称性,为CRN的稳健性属性分析提供了蓝图。

Chemical reaction networks (CRN) comprise an important class of models to understand biological functions such as cellular information processing, the robustness and control of metabolic pathways, circadian rhythms, and many more. However, any CRN describing a certain function does not act in isolation but is a part of a much larger network and as such is constantly subject to external changes. In [Shinar, Alon, and Feinberg. "Sensitivity and robustness in chemical reaction networks." SIAM J App Math (2009): 977-998.], the responses of CRN to changes in the linear conserved quantities, called sensitivities, were studied in and the question of how to construct absolute, i.e., basis-independent, sensitivities was raised. In this article, by applying information geometric methods, such a construction is provided. The idea is to track how concentration changes in a particular chemical propagate to changes of all concentrations within a steady state. This is encoded in the matrix of absolute sensitivities. As the main technical tool, a multivariate Cramer-Rao bound for CRN is proven, which is based on the the analogy between quasi-thermostatic steady states and the exponential family of probability distributions. This leads to a linear algebraic characterization of the matrix of absolute sensitivities for quasi-thermostatic CRN. As an example, the core module of the IDHKP-IDH glyoxylate bypass regulation system is analyzed analytically and numerically by extensive random sampling of the concentration space. The experimentally known findings for the robustness of the IDH enzyme are confirmed and a hidden symmetry at the level of distributions is revealed, providing a blueprint for the analysis of the robustness properties of CRNs.