[论文解读] Output-Oblivious Stochastic Chemical Reaction Networks

Ben Chugg, Hooman Hashemi+1 more|arXiv (Cornell University)|Dec 7, 2018

Gene Regulatory Network Analysis被引用 2

一句话总结

本文刻画了在存在领导者的情况下，可通过输出无关的随机化学反应网络（CRNs）稳定计算的函数 f : N² → N 的类别。证明表明，此类函数恰好是那些单调递增且满足以下条件之一的函数：在模某个周期的同余类上为分段仿射（即网格仿射）的函数，或为有限个部分分裂函数的最小值函数（即行为类似于最小值函数的函数）。该结果解决了CRN组合中的关键挑战，确保计算过程中输出物种不会被消耗。

ABSTRACT

We classify the functions $f:\mathbb{N}^2 ightarrow \mathbb{N}$ which are stably computable by output-oblivious Stochastic Chemical Reaction Networks (CRNs), i.e., systems of reactions in which output species are never reactants. While it is known that precisely the semilinear functions are stably computable by CRNs, such CRNs sometimes rely on initially producing too many output species and then consuming the excess in order to reach a correct stable state. These CRNs may be difficult to integrate into larger systems: if the output of a CRN $\mathcal{C}$ becomes the input to a downstream CRN $\mathcal{C}'$, then $\mathcal{C}'$ could inadvertently consume too many outputs before $\mathcal{C}$ stabilizes. If, on the other hand, $\mathcal{C}$ is output-oblivious then $\mathcal{C}'$ may consume $\mathcal{C}$'s output as soon as it is available. In this work we prove that a semilinear function $f:\mathbb{N}^2 ightarrow \mathbb{N}$ is stably computable by an output-oblivious CRN with a leader if and only if it is both increasing and either grid-affine (intuitively, its domains are congruence classes), or the minimum of a finite set of fissure functions (intuitively, functions behaving like the min function).

研究动机与目标

• 确定哪些函数 f : N² → N 可通过输出无关的CRNs稳定计算，其中输出物种不会作为反应物被消耗。
• 解决CRN组合中的挑战，即下游网络可能在稳定前过早消耗输出。
• 刻画在更强的输出无关性约束下可计算的函数类别，超越标准的半线性函数。
• 通过引入输出单调性的实际要求和可组合性，扩展先前关于CRNs中稳定计算的研究。

提出的方法

• 作者在输出物种永不作为反应物的约束下，分析半线性函数，并使用领导者物种启动计算。
• 他们定义并分析了两类关键函数类：网格仿射函数（在模某个周期的同余类上为分段仿射）和部分分裂函数最小值（具有类似最小值函数的分段线性行为）。
• 证明利用了在网格上对半仿射函数的结构分解，借助周期性和将定义域划分为楔形和线段区域。
• 引理建立了在网格上的输出单调性和递增性，表明非单调行为将导致矛盾。
• 分析基于子域上仿射函数的性质，并使用阈值集来定义不同网格偏移下的定义域交集。
• 一个关键技术工具是引理18，证明在单个网格上，单调递增的半仿射函数必须在所有点上非递减；引理19则将其扩展至网格间的转换。

实验结果

研究问题

• RQ1 哪些函数 f : N² → N 可通过输出无关的CRNs稳定计算，其中输出物种不会作为反应物被消耗？
• RQ2 在更强的输出无关性约束下，可稳定计算的函数类是否可以被刻画，超越已知的半线性函数类？
• RQ3 此类函数必须满足哪些结构性质（例如单调性、分段线性）才能在此模型中被计算？
• RQ4 网格仿射函数和分裂函数最小值如何自然地出现在输出单调半线性函数的分解中？
• RQ5 该刻画是否可推广至多于两个输入的函数？

主要发现

• 当且仅当函数 f : N² → N 是单调递增的，且为网格仿射函数或有限个部分分裂函数的最小值时，该函数可通过具有领导者的输出无关CRN稳定计算。
• 输出无关可计算函数的类严格小于所有半线性函数的类，因为它同时要求单调性和特定的分段结构。
• 网格仿射函数是指在模某个周期 p 的每个同余类上为仿射的函数，且在不同网格偏移下保持一致。
• 分裂函数是具有‘拐角’结构的分段线性函数，其最小值可模拟类似最小值函数的行为。
• 该刻画在存在单个领导者分子的假设下成立，这使得在随机模型中能够实现精确控制和稳定。
• 结果可推广至多于两个输出的函数，因为在输出无关约束下，每个输出可独立处理。