DNA Replication under Thermal, Chemical, and Genotoxic Stress

技術要約：熱、化学、および遺伝毒性ストレス下におけるDNA複製

問題提起
忠実なゲノム複製はゲノムの完全性に不可欠であるが、真核生物のDNA複製は、熱、化学、および遺伝毒性ストレスによる変動にもかかわらず、堅牢性を維持しなければならない。Saccharomyces cerevisiae（出芽酵母）は、レプリケーションを研究するためのモデルシステムであるが、従来の視点では、その複製プログラムは主に決定論的であり、オリジンの発火時刻や一定のフォーク速度を持つものとして扱われることが多い。しかし、実験的証拠は、オリジンの発火効率における実質的な確率性と、レプリケーションフォーク速度の著しい不均一性を示している。既存の確率モデルは、多くの場合、粗い解像度（例：約1 kbp）で動作するか、あるいは一定の平均フォーク速度を仮定しており、微細なダイナミクス、稀に発生する長期の複製イベント、または細胞が多様なストレス条件下で適応する具体的なメカニズムを捉えることができていない。微細なスケールのダイナミクス、稀な長期複製イベント、およびストレス条件下での細胞の適応メカニズムを理解するために、確率的なオリジン発火、フォーク速度分布、およびリソース制限を単一塩基対解像度で統合する定量的な枠組みが必要とされている。

方法論
著者らは、S. cerevisiae の全ゲノム複製を単一塩基対（1 bp）解像度で扱う、格子ベースの確率的モンテカルロ・フレームワークを開発した。このモデルは、ゲノムを未複製状態（0）から複製状態（1）へと遷移する1次元のストランドとして表現する。フレームワークは、実験に基づいた2つの主要なパラメータによって制御される：

1. 確率的オリジン発火： モデルは、適合性（$P$）、中央値発火時刻（$t_{1/2}$）、および発火時間の幅（$t_w$）によって特徴付けられる、459個の潜在的なオリジンに関する実験的に決定されたデータを利用する。発火は、時間依存的なガウス確率分布としてモデル化される。オリジン間の大きな隙間を防ぐため、ゲノムは分割され、オリジンは一定の割合で確率的にサンプリングされる。
2. 制限因子 ($L(t)$)： 時間依存的な制限因子は、オリジン発火を開始するために必要な有限の細胞内リソースを表す。オリジンが発火するためには、2つの制限因子が利用可能である必要がある。これらの因子は、活性なレプリケーションフォークを形成するために動員され、フォークの収束（合流）時に放出され、その後のオリジン発火のために再利用（リサイクル）される。$L(t)$ の時間的ダイナミクスは、実験的な活性フォーク統計量によって制約されたシグモイド関数に従う。
3. フォーク速度分布： 一定の速度を仮定するモデルとは異なり、本フレームワークは、位置（$\xi$）、幅（$\omega$）、および歪度（$\alpha$）によってパラメータ化された歪んだガウス分布（$f(v)$）から抽出されたフォーク速度を割り当てる。これらのパラメータは、温度やストレス条件に基づく実験データまたは線形補外から導出される。

モデルは、熱、化学（ヒドロキシ尿素）、および遺伝毒性ストレスへの適用前に、実験的なレプリケーションプロファイル（活性フォーク数、複製分率、オリジン間距離、および時空間的レプリケーションプロファイル）に対してベンチマーク検証が行われる。

主な貢献

• 単一塩基対解像度： 本フレームワークは、最も物理的に意味のある解像度（1 bp）を達成し、タイミングとフォークの伝播を明示的に追跡することで、微細なダイナミクスや局所的な変異パターンの解析を可能にする。
• 最小限のパラメータ化： モデルは、条件固有の調節メカニズムを導入することなく、時間依存的な制限因子とフォーク速度分布という2つの有効なパラメータのみを用いて、多様なストレス応答を再現する。
• 不均一性の統合： フォーク速度の不均一性と確率的なオリジン発火を明示的に組み込み、決定論的または一定速度の仮定を超越している。

結果

• ベンチマーク： モデルは、23°Cにおける染色体14の活性フォーク数、ゲノム複製分率、および時空間的レプリケーションプロファイル（STRP）の実験データを定量的に再現する。
• Erland分布に従うS期持続時間： フォーク速度の不均一性を導入することで、総複製時間（$T_{rep}$）の分布が、顕著な右側に歪んだ裾を持つアーラン分布に従うようになる。これは、E. coli やヒト細胞株で観察される、稀に発生する異常に長い複製イベントの出現を説明するものであり、S. cerevisiae においても同様の挙動を予測する。一定のフォーク速度を用いたシミュレーションでは、これらの長時間イベントを再現できない。
• 熱ストレス： モデルは、総レプリケーションフォーク数の温度に対する非単調な依存性を予測する。フォーク数は23°C付近でピークに達する。低温では、フォーク速度が遅いために合流と制限因子の隔離が遅れ、新しいオリジン発火が抑制される。高温では、高速なフォーク速度により、多くのオリジンが発火前に受動的に複製される。モデルは、フォーク速度分布のパラメータの線形補外を用いることで、様々な温度における実験的な $T_{rep}$ 値を正確に再現する。
• 化学ストレス（ヒドロキシ尿素）： ヒドロキシ尿素（HU）ストレス下では、フォーク速度の低下だけでは実験的なレプリケーションプロファイルを一致させるには不十分であることをモデルは示す。データを正確に再現するには、制限因子の利用可能性を（元の値の約30%まで）同時に減少させる必要がある。モデルは、フォーク速度とHU濃度の間のべき乗則スケーリング関係を予測する。
• 遺伝毒性ストレス： 本フレームワークは、様々な S. cerevisiae 変異株（ヒストンアセチル化、dNTPプール、およびフォーク保護複合体に影響を与えるもの）におけるレプリケーションダイナミクスを正常にモデル化する。それは、特定の変異がフォーク速度分布をどのように変化させるか（例：rtt109$\Delta$）を捉え、場合によっては制限因子関数を調節し、総複製時間に変化をもたらす様子を捉えている。

意義
本論文は、多様な条件下での堅牢なゲノム複製は、制限的なレプリケーションリソースへの競争とフォーク速度の不均一性によって駆動される最小限のフレームワークから生じるものであると主張している。本研究は、S期の進行の急激な加速やS期の持続時間の分布を含む、グローバルなレプリケーションダイナミクスが、これらの基本的なメカニズムからどのように生じるかについての物理的な解釈を提供している。単一塩基対解像度で動作することにより、本フレームワークは、点変異や疾患に関連する局所的なゲノム領域のレプリケーションダイナミクスを調査するための汎用性の高いプラットフォームを提供する。著者らは、本研究が、熱、化学、および遺伝毒性の条件下でのレプリケーション制御を、共有された一連の有効なパラメータを通じて統一する方法を理解するための、定量的な基礎を確立するものであると結論付けている。