Conditions for Darwinian evolution in compartmentalized autocatalytic reaction networks

技術要約：区画化された自己触媒反応ネットワークにおけるダーウィン的進化の条件

問題提起
生命の起源には、遺伝、変異、および差異的な生殖が可能な「ダーウィン的集団」への移行が含まれていた可能性が高い。テンプレートに基づく複製体（例：RNA）に焦点を当てる「複製優先（replication-first）」シナリオに対し、「代謝優先（metabolism-first）」シナリオは、初期のプロトセルが、化学組成そのものが表現型として機能する自己触媒化学系（ACS）に依存していたことを提案している。一つの未解決の重要な問いは、そのようなシステムにおける組成の遺伝が、初期の区画化サイクルが不完全であった可能性を考慮しても、自然選択を維持できるほど十分に堅牢であるかどうかである。先行研究では、ACSの組成は、分裂プロセスを通じて存続するには不安定すぎるか、あるいは必要な変異を生み出すには安定すぎるかのいずれかである可能性が示唆されてきた。本論文は、マルチステーブル（多安定）な自己触媒ネットワークが、成長・分裂サイクルを通じて継承可能な組成状態を維持し、テンプレート複製体なしにダーウィン的集団を出現させるための具体的な条件を扱う。

方法論
著者らは、成長および分裂する区画内において、共有基質（$S$）を消費する競合する二つの自己触媒種（$X_1$ および $X_2$）のダイナミクスを分析するために、決定論的および確率論的な数学的モデリングを組み合わせて用いている。

1. 反応ネットワークモデル： 本研究は、二つの対称な自己触媒種が競合する最小限の二安定システムに焦理している。繁殖率は関数 $r(x)$ によって制御され、成長率は $r(x)x$ となる。著者らは、線形（$r(x)x = \epsilon + \kappa x$）および非線形（$r(x)x = \epsilon + \kappa x^2$）の繁殖則の両方を分析する。
2. 成長および分裂（GD）プロトコル： ダイナミクスは、以下の3つのレジームの下でモデル化される：
   • 連続希釈（SD）： 区画が時間間隔 $\Delta t$ の間成長し、その後因子 $m$ で希釈される離散的なサイクル。
   • 連続槽型反応器（CSTR）： 一定の速度での流入と希釈を伴う連続的なプロセス。
   • 一般化されたGD： 分裂前に体積成長が冪乗則（$dV/dt \propto V^\alpha$）に従う中間的なプロトコル。
3. 安定性解析： 著者らは、対称な軌跡（$x_1 = x_2$）の安定性を分析することにより、二安定性（二つの異なる安定した組成状態の存在）の条件を導出する。ポアンカレ写像を用いて不動点を決定し、希釈間隔（$\Delta t$）および希釈率（$\bar{\phi}$）の臨界閾値を特定するための分岐図を利用する。
4. 変異と選択の組み込み：
   • 変異： 有限体積における分子ノイズをシミュレートするために、ギレスピー・アルゴリズムを介して確率的ゆらぎを導入し、状態間の遷移時間を評価する。
   • 差異的な生殖： モデルを拡張し、希釈率（したがって実効的な成長率）が内部の化学組成に依存するようにすることで、選択圧をシミュレートする。
5. 集団ダイナミクス： $N$ 個の並行する系統において、様々な選択圧の下でのダーウィン的集団の出現を観察するために、ライト・フィッシャー型のプロセスをシミュレートする。
6. 実験への適用： 理論的枠組みを、既知の実験的ACSであるAzoarcusリボザイム・システムに適用し、組成の遺伝が観察可能なパラメータ領域を予測する。

主な貢献および結果

• 非線形性の必要性： 著者らは、ACSが連続希釈下で遺伝（二安定性）を示すためには、繁殖率関数 $r(x)x$ が非線形（具体的には凸関数）でなければならないことを証明した。線形の成長則は単一の安定な対称状態をもたらし、組成情報を消去してしまう。
• 希釈の臨界閾値： 二安定性は保証されているわけではなく、特定のパラメータ領域内にのみ存在する。本研究は、臨界希釈間隔（$\Delta t_c$）および臨界希釈率（$\bar{\phi}_c$）を特定した。間隔が長すぎるか、あるいは希釈率が高すぎる場合、システムは二安定性を失い、単一の状態へと収束する。
• 一般化されたプロトコルの境界設定： 中心的な発見は、一般の成長・分裂プロトコルにおける二安定性を支持するパラメータ領域が、連続希釈（SD）および連続槽型反応器（CSTR）プロトコルの限界によって厳密に制限されていることである。具体的には、あるシステムがSDの下で二安定であれば、反応キネティクスが凸である限り、一般のプロトコル下でも二安定となる。これにより、SDが遺伝性をテストするための保守的な十分条件であることが確立された。
• ノイズと選択に対する堅牢性：
  • 変異： 確率的ノイズは状態間の遷移を引き起こすが、著者らは、遷移のタイムスケールが成長・分裂サイクルよりも十分に長い限り、遺伝性は堅牢に維持されることを示した。変異と遺伝が共存できるパラメータ空間が存在する。
  • 差異的な生殖： 本研究は、逆説的な結果を明らかにしている。過度に強い選択圧（成長率の大きな差）は、二安定性を破壊し、システムを単一安定な状態へと崩壊させ、進化の可能性を排除してしまう。
• ダーウィン的進化の実証： プロトセルの集団におけるシミュレーションにおいて、遺伝、変異、および差異的な生殖が特定されたパラメータ空間内でバランスが取れているとき、集団は進化することを示した。選択が特定の状態を好む場合、その状態が支配的となり、圧力が逆転すると集団は適応する。逆に、ノイズが高すぎるか選択が強すぎる場合、集団は進化に失敗する。
• 実験的実現可能性： Azoarcusリボザイム・システムにモデルを適用することで、必要な非線形性を導入するための代謝結合が含まれている限り、連続希釈間隔が50–125分、希釈因子が1サイクルあたり2.5–11倍であれば、組成の遺伝が観察可能であると予測している。

意義および主張
本論文は、初期の地球において、テンプレート複製体の存在なしに、自己触媒化学システムを介したダーウィン的集団の出現が実現可能であることを確立している。著者らは、反応の非線形性、希釈プロトコル、およびノイズと選択のバランスに関する具体的なパラメータ領域（特に、組成の遺伝が維持される領域）を特定したと主張している。

本研究は、区画化された自己触媒システムにおける自然選択のための厳格な条件を提供している。著者らは、初期の区画化サイクルは不完全であったかもしれないが、マルチステーブルなネットワークは、幅広い成長・分裂レジームの下で継承可能な状態を維持できると論じている。これらの知見は、オープンエンドな進化（新しい自己触媒コアの出現）にはさらなる複雑さが必要であることを示唆しつつも、生命の起源のための最小要件を理解するための基礎を築いており、「代謝優先」のシナリオが、化学キネティクスと環境的制約が導出された安定性の境界内に収まる限り、実行可能な経路であることを示唆している。結論として、オープンエンドな進化にはさらなる複雑さが必要であるが、組成状態の継承と選択という根本的なメカニズムは、これらのモデルの中で堅牢に確立されている。