Low-frequency output fluctuations in an open exclusion process with particle pausing

技術要約：粒子の一時停止を伴う開放型排除過程における低周波出力変動

問題提起
相互作用する粒子系（TASEPなどの全非対称単純排除過程）における駆動輸送は、局所的な確率的事象を集団的なゆらぎへと変換する。標準的なTASEPはリボソームやポリメラーゼの交通を記述するが、生物学的システムでは、粒子が活性状態と一時停止（ポーズ）状態の間を可逆的に切り替わる（例：転写のポーズ）ような、遅い内部状態を伴うことが多い。先行研究では、可逆的なポーズが定常電流を減少させ、平均ポーズ粒子数（$N_p$）が重要な制御変数となる有限サイズ領域が存在することを確立しているが、特に $N_p$ がオーダー1である場合である。しかし、この領域における出力（完了した脱出）の時系列統計については未解決のままである。具体的には、平均電流によって特定された「有限サイズ・ポーズ領域」が、完了した脱出の定常統計において明確なシグネチャを持つのか、また、遅い内部状態が単にスループットを低下させるだけでなく、出力の時系列的組織化をどのように再形成するのかが不明である。

手法
本研究では、格子長 $L$ ($50 \le L \le 500$) の有限開放型TASEPを用いる。粒子はレート $\epsilon$ で一方向にホップし、レート $\alpha$ で流入し、レート $\beta$ で脱出する。極めて重要な点は、粒子が活性状態（$A$）とポーズ状態（$P$）の間をレート $k_p$（ポーズ）および $k_u$（解除）で可逆的に切り替わることである。脱出サイトにいるポーズ中の粒子は、脱出する前に活動を再開しなければならない。

動力学は、イベント駆動型キネティック・モンテカルロ・アルゴリズムを用いてシミュレーションされる。主な観測量は以下の通りである：

1. 出力計数統計： 累積脱出数 $N(t)$ を観測窓の長さ $T_w$ にわたって分析する。本研究では、短時間の計数ノイズと遅い出力変調を区別するために、窓依存のファノ因子 $F(T_w) = \text{Var}(N_w)/\langle N_w \rangle$ に焦点を当てる。
2. 構造的観測量： ポーズ粒子数（$N_p = L\rho_{\text{paused}}$）、最大連続粒子クラスターのサイズ（$C_{\text{max}}$）、および特定の状態（ポーズなし vs ポーズあり、または弱混雑 vs 強混雑）の滞在時間。
3. 解析的近似： ポーズ数の動力学を「流入－死亡過程（immigration-death process）」として近似する最小限の「定数出生－線形死亡」モデルを構築し、ノイズ極大値の所在を合理化することを目的とする。

主な貢献と結果

• 非単調な低周波ノイズ： ポーズ集団（$N_p$）が増加すると、平均定常電流（$J$）は滑らかに低下するが、低周波ファノ因子（$F_\infty$）は強い非単調挙動を示す。$F_\infty$ は、平均ポーズ粒子数がオーダー1（$N_p \simeq 1.5\text{--}2$）のときに顕著な極大値に達する。この極大値は、格子長 $L=50$ から $500$ において普遍的に観察される。
• 有限サイズ・スケーリング： このノイズ極大に達するために必要な微視的なポーズ率（$k_p$）は、$k_p^{\text{max}} \propto L^{-1}$ とスケーリングする。これは、長い格子ほどポーズにさらされる粒子が多くなるため、臨界的なオーダー1のポーズ粒子数（$N_p^\star$）を維持するためには、より小さな粒子あたりのレートで十分であることを意味している。
• ゆらぎのメカニズム： 最大ノイズは、ポーズが稀な場合（規則的な出力）でも、持続的な混雑が生じている場合（抑制されているが安定した出力）でもなく、中間的な領域で発生する。それは、生産的な間隔と、弱あるいは消失した出力が続く間隔が交互に現れる「間欠的レジーム」において生じる。この交互作用は、ポーズによって開始される大規模な粒子クラスター（ジャム）の核形成と溶解によって駆動される。
• 解除率（$k_u$）の役割： $N_p$ はクロスオーバーの「位置」を規定するが、$F_\infty$ の「振幅」を完全に決定するわけではない。固定された $N_p$ において、解除が遅い（小さい $k_u$）ほど、相関時間（$\tau_c$）とノイズ振幅（$F_\infty$）の両方が増大する。これは、ノイズの振幅が、ポーズの生存時間とジャムの形成および緩和のタイムスケールの比に依存していることを示唆している。
• 構造的相関：
  • 相関時間： ポーズ粒子の出現と消失に関連する滞在時間のスケール（$N_{\text{paused}}=0$ vs $\ge 1$）は、フィッティングされた出力相関時間（$\tau_c$）を密接に追跡している。対照的に、閾値に基づいた混雑の定義（$C_{\text{max}}$に基づく）は、$N_p$ が増加しても $\tau_c$ を追跡できない。
  • ノイズ振幅： 低周波ノイズの振幅は、最大クラスターの平均サイズではなく、最大クラスターサイズの分散（$\text{Var}(C_{\text{max}})$）と強く相関している。最大ノイズは、支配的なクラスターが広範なサイズを探索するとき、すなわち、繰り返される成長、収縮、および再編成を反映して発生する。
• 最小理論の限界： 独立したポーズ（ポアソン分布）を仮定する最小限の解析モデルは、$N_p^\star \simeq 1.50$ で極大値を予測する。シミュレーションでは、実際の極大値は $N_p \simeq 1.9$ 付近で発生しており、ポーズ数は過分散（分散 $>$ 平均）である。これらの乖離は、空間的相関、境界効果、および排除による混雑に起因する。

意義
本論文は、有限開放型排除過程における可逆的なポーズが、平均スループットが単に減少するだけでなく、分子生産のタイミングが高度に間欠的となる特有のゆらぎレジームを生み出すことを示している。本研究は、最も重大な出力ノイズは、欠陥の平均密度からではなく、少数の欠陥（オーダー1）が交通ジャムと相互作用する際の集団的な動力学から生じることを特定した。

本研究は、以下の2つのレベルの組織化を区別している：

1. ポーズ数ダイナミクスが、出力のデコリレーション（脱相関）のタイムスケールを設定する。
2. ジャムの再編成ダイナミクス（最大クラスターのゆらぎ）が、ノイズの振幅を決定する。

このことは、転写完了時間に敏感なプロセス（例：下流の遺伝子発現）にとって、プロモーターのスイッチングや配列特異的なポーズ部位を必要とせずとも、伸長中の交通量のみによって、生成物の時間的デリバリーが極めて変動し得ることを示唆している。本知見は、有限サイズ効果の理解を拡張するものであり、「有限サイズ・レジーム」が、単なる電流の減少だけでなく、遅い内部状態と集団的なジャミングの相互作用によって駆動される、特有の非単調な低周波出力ノイズの増大によって特徴付けられることを示している。