[論文レビュー] Rate-Independent Computation in Continuous Chemical Reaction Networks

Ho-Lin Chen, David Doty+2 more|arXiv (Cornell University)|Jul 29, 2021

Advanced biosensing and bioanalysis techniques参考文献 53被引用数 44

ひとこと要約

この論文は、連続化学反応系（CRN）に対するレート独立の計算概念を定義し、任意の時間変動する速度律関数の下で計算可能な実数値関数を正確に特徴づける。セグメント到達性フレームワークと負値のデュアルレール表現を用いて。

ABSTRACT

Coupled chemical interactions in a well-mixed solution are commonly formalized as chemical reaction networks (CRNs). However, despite the widespread use of CRNs in the natural sciences, the range of computational behaviors exhibited by CRNs is not well understood. Here we study the following problem: what functions $f:\mathbb{R}^k o \mathbb{R}$ can be computed by a CRN, in which the CRN eventually produces the correct amount of the "output" molecule, no matter the rate at which reactions proceed? This captures a previously unexplored, but very natural class of computations: for example, the reaction $X_1 + X_2 o Y$ can be thought to compute the function $y = \min(x_1, x_2)$. Such a CRN is robust in the sense that it is correct no matter the kinetic model of chemistry, so long as it respects the stoichiometric constraints. We develop a reachability relation based on "what could happen" if reaction rates can vary arbitrarily over time. We define *stable computation* analogously to probability 1 computation in distributed computing, and connect it with a seemingly stronger notion of rate-independent computation based on convergence under a wide class of generalized rate laws. We also consider the "dual-rail representation" that can represent negative values as the difference of two concentrations and allows the composition of CRN modules. We prove that a function is rate-independently computable if and only if it is piecewise linear (with rational coefficients) and continuous (dual-rail representation), or non-negative with discontinuities occurring only when some inputs switch from zero to positive (direct representation). The many contexts where continuous piecewise linear functions are powerful targets for implementation, combined with the systematic construction we develop for computing these functions, demonstrate the potential of rate-independent chemical computation.

研究の動機と目的

• 反応速度に関係なくCRNが信頼性高く実行できる計算がどれであるかを理解する動機づけ。
• CRNの反応式量論（ストイキオの計算能力）を、速度法則と切り離して特定する。
• 任意の速度変動下で“何が起こりうるか”をとらえる到達性フレームワークを開発する。
• 連続性と分節線形性の観点からレート独立計算可能な関数を特徴づける。
• 負値を扱いモジュールの構成を可能にするデュアルレール符号化を示す。

提案手法

• 反応式（ストイキオメトリ）と因果制約を満たすセグメント（直線）軌道に基づく広範な到達性関係を定義する。
• 分散システムの確率的計算に類似した安定計算を形式化する。
• 一般化された速度律下の公正計算を導入し、前向きCRN（フォワード型CRN）の安定計算との関係を証明する。
• レート独立計算可能性は、正の連続性を持つ分節線形関数（直接表現）または連続分節線形関数（デュアルレール表現）に対応することを証明する。
• デュアルレール表現を用いて負の値とCRNモジュールの組み合わせを可能にする。

実験結果

リサーチクエスチョン

• RQ1 反応速度が時間的に任意に変化する場合、CRNはレート独立の意味でどの実数値関数を計算できるか？
• RQ2 広範な速度法則下でCRNにおける“何が起こりうるか”をとらえる到達性をどう formalizeするか？
• RQ3 直接表現とデュアルレール表現の下でのレート独立計算可能な関数の正確な特徴づけは何か？
• RQ4 レート独立モジュールを計算可能性を損なうことなく組み合わせられるか？
• RQ5 連続CRNにおける安定計算と公正計算の概念は、特にフォワード型ネットワークにおいてどのように関係するか？

主な発見

• 広く許容される速度律に基づく到達性関係は、非負性、反応物の利用可能性、および種の因果的生成という直感的制約と同値である。
• フォワード型CRNでは安定計算と公正計算が一致し、堅牢な収束とレート独立性を結びつける。
• 正確には、レート独立計算可能な関数は正の連続的な分節線形関数（直接表現）または連続的な分節線形関数（デュアルレール表現）である。
• デュアルレール符号化は干渉なしにCRNモジュールの組み合わせを可能にし、負値出力をサポートする。
• 非入力種の初期濃度が非ゼロの場合、レート独立CRNsは表現力においてReLUニューラルネットワークに似る。
• 質量作用法および他の速度法則（Michaelis–MentenおよびHill関数動力学を含む）を、統一された到達性モデルの下で受け入れる。