DNA 複製から見える、情報を時間へ展開する構造

情報という言葉は、しばしば「何が記録されているか」という内容と結び付けて考えられる。文章なら文字列、計算機ならビット列、DNA なら A、T、G、C の塩基配列である。この見方では、情報を保存するとは内容を失わずに保持することであり、同じ内容が残っていれば同じ情報が残っていると考えられる。

しかし、保存されている内容と、その内容が実際にどう使われるかは同じ問題ではない。文章なら、同じ文字列でもどの規則で解釈するかによって処理結果が変わる。計算機では、同じビット列でも命令として実行するのか、数値として読むのか、単なるデータとして保存するのかによって働きが変わる。DNA でも、塩基配列が存在することだけから、その領域が細胞周期のどの時点で複製されるかまでは決まらない。

この違いを扱うには、「内容」と「実行条件」を分ける必要がある。内容は何が保存されているかを規定する。一方、実行条件は、その内容のどの部分が、どの時点で、どの順序で処理されるかを制約する。情報がどこに置かれ、周囲とどのような関係を作り、現在どの状態にあるかが実行条件へ影響するなら、構造は情報を収める容器ではない。内容を実際の処理へ接続する状態変数の一部になる。

既稿では、差異が存在するだけでは意味や機能にはならず、それを読み取り、更新につなげる構造が必要だと考えた[1]。本稿では、その議論を時間方向へ拡張する。内容が変わっていなくても、構造状態が異なれば実行の可否だけでなく、実行される時期や順序まで変わるのかを考える。さらに、構造そのものが時間とともに変化するなら、一時点の配置だけを記録しても実行条件を十分には表せない。

この問いを具体化する題材として、DNA 複製と染色体の 3 次元構造を扱う。細胞では、ゲノムという複製対象が存在するだけでなく、各領域をいつ複製するかという時間的な秩序が成立している。その秩序と核内構造の時間変化を同じ状態記述へ入れることで、情報の実行条件を静的な構造ではなく、構造が時間の中でたどる軌道として表せるかを検討する。最終的には、この関係を数理モデルとして記述し、どの観測結果が仮説を支持し、どの結果なら見直す必要があるかまで明示する。

1. 同じ情報でも、同じようには使われない

細胞が分裂するためには、その前に DNA を複製し、二つの娘細胞へ同じゲノムを渡せる状態にする必要がある。ただし、哺乳類の DNA は複製期間の開始と同時に全領域が一斉に複製されるわけではない。S 期の早い時点で複製される領域もあれば、後半まで複製されない領域もある。この領域ごとの時間的な秩序は複製タイミングと呼ばれ、ゲノム上で一定のまとまりを持つだけでなく、細胞の種類や発生状態とも関係している[2]。

この現象は、DNA 配列だけでは細胞の現在の実行状態を表せないことを示している。塩基配列を読めば、ある領域にどの塩基が並んでいるかは分かる。しかし、その配列自体に「S 期開始直後に複製する」「後半まで待つ」という時刻がそのまま書き込まれているわけではない。同じゲノムを持つ細胞でも、細胞の種類や状態が変われば複製タイミングのパターンは変わり得る[2]。複製対象の内容を知ることと、その対象が現在どの順序で処理されるかを知ることは別である。

理由は二段階に分けられる。第一に、DNA 配列は複製される対象を規定するが、複製装置が現在どの領域へ作用できるかという細胞状態までは記述しない。第二に、実際の複製開始には、細胞周期の進行、複製因子の状態、DNA やクロマチンに付随する化学的状態など、配列以外の条件が関与する。そのため、配列を固定しても、それらの条件が変われば実行時刻まで同じになるとは限らない。

複製タイミングが毎回完全な偶然で決まるのであれば、領域ごとの早い、遅いという秩序を別の状態として考える必要はない。しかし、実際には複製タイミングにはゲノム上の秩序があり、細胞状態との対応も観測される[2]。一定の条件で似た時間順序が繰り返されるなら、その順序を再現するための状態が細胞内に存在すると考える必要がある。

その状態を一つの要因へ還元することはできない。DNA 配列そのものに加えて、化学修飾、複製因子の量や活性、細胞周期上の位置などが候補になる。さらに、DNA は核内で直線状に並んでいるのではなく、折り畳まれた染色体として 3 次元的な配置を作る。あるゲノム領域がどの領域と近接し、核内のどの構造的な区画へ属しているかも、配列とは別に記述できる状態である。

観点	説明できること	それだけでは不足すること
DNA 配列	どの塩基がどの順序で並び、何が複製対象になっているかを記述できる。	その領域が現在の細胞周期でいつ複製されるかという実行時刻までは表せない。
分子状態と細胞周期	複製因子の状態や、細胞が複製可能な段階へ進んでいるかを記述できる。	ゲノムの各領域が核内でどのような位置関係を作っているかまでは表せない。
構造状態	各領域が核内でどの領域と近接し、どの構造的な区画に属しているかを記述できる。	複製される DNA の内容そのものや、構造以外の複製条件までは決まらない。
実行結果	DNA 合成が始まったか、どの領域が先に複製されたかを観測できる。	観測された順序だけから、それを生じさせた原因を一意には特定できない。

この区別を置くと、染色体構造について検証すべき関係も明確になる。構造と複製タイミングが対応しているだけなら、両者が同じ細胞状態を反映しているだけかもしれない。構造が複製を制御する状態の一部だと主張するには、構造に異常を生じさせる条件で、DNA 配列を大きく変えなくても複製の開始や時間的な順序まで変化するかを見る必要がある。

さらに、その変化が構造そのものに由来するのか、構造を変えるために操作した分子が複製装置へ直接作用した結果なのかも分けなければならない。構造の変化と複製の変化が同時に観測されたというだけでは、因果方向は確定しない。必要なのは、配列、分子状態、構造状態、実行結果を分けて観測し、どの状態変化がどの結果に先行するかを追うことである。この条件を満たして初めて、染色体の 3 次元構造を、複製結果に付随する形ではなく、情報を時間的な実行へ接続する状態として評価できる。

2. 配置は、ただの置き場所ではない

DNA は細胞核の中に無秩序に詰め込まれているわけではない。長い DNA 分子は染色体として折り畳まれ、ある領域同士は頻繁に接触し、別の領域同士は接触しにくいという偏りを持つ。この接触関係をゲノム全体で見ると、染色体は大きく A コンパートメントと B コンパートメントに分けられる。A コンパートメントには遺伝子活動が比較的活発な領域が多く、B コンパートメントにはより凝縮した領域が多い。そして、A コンパートメントは早く複製される領域、B コンパートメントは遅く複製される領域と強く対応している。

この対応が意味するのは、ゲノム上の位置と核内の位置が異なる状態を記述するということである。塩基配列を読めば、ある領域がゲノムのどこにあり、何が書かれているかは分かる。しかし、その領域が核内でどの領域と近接し、活動性の高い区画と低い区画のどちらへ組み込まれているかは、配列だけから直接読み取ることはできない。染色体には一次元の配列とは別に、領域同士の関係から成る 3 次元の状態が存在する。

既稿「境界が生命の形を作る」では、組織を構成する細胞の種類だけを列挙しても、その形を十分には説明できないことを扱った[3]。同じ種類の細胞でも、どの細胞と接し、どちら側に境界を持ち、どの方向を向くかによって組織内での役割が変わる。部品そのものの性質に加えて、部品同士の関係が機能条件を作るからである。染色体でも同じ観点を置ける。塩基配列という部品の内容だけでなく、ゲノム領域同士が核内でどのような関係を作っているかが、複製を含む細胞内処理の条件になっている可能性がある。

その可能性を支持する早い証拠の一つが、複製タイミングと染色体の長距離相互作用との対応である。2010 年の研究では、複製タイミングのパターンと Hi-C で測定した長距離の染色体相互作用との間に強い対応が報告された[4]。早く複製される領域は早く複製される領域と、遅く複製される領域は遅く複製される領域と空間的にまとまりやすかった。複製される時刻という時間的な分類と、核内で誰と接触するかという空間的な分類が、独立ではなかったことになる。

ただし、この結果だけでは因果関係は決まらない。第一に、空間構造が複製タイミングを制御している可能性がある。特定の領域同士が同じ核内環境へ集まることで、複製因子へ接近しやすくなったり、逆に複製開始が抑えられたりするなら、配置が実行時刻へ影響する。第二に、空間構造と複製タイミングの双方を別の分子機構が制御している可能性もある。その場合、両者に強い相関があっても、一方が他方の直接原因とは限らない。

2021 年に RIF1 と PP1 の相互作用を操作した研究は、この区別が必要であることを示している。RIF1 は複製タイミングの制御と核内構造の双方に関与するが、実験条件によっては特定の核内構造の変化と複製タイミングの変化を部分的に分離できた[5]。つまり、空間構造と時間制御は同じ細胞状態の中で密接に結び付いていても、常に一対一で対応する変数ではない。A/B コンパートメントを観測しただけで、その細胞の複製タイミングが機械的に決まると考えるのは強すぎる。

この点は、本稿で構造を「実行条件」と呼ぶときにも制約になる。構造が実行条件の一部であることと、構造だけで実行結果が決まることは異なる。複製には、細胞周期、複製因子、RIF1 や PP1 を含む分子状態など複数の条件が関与している。その中で染色体構造が独自の寄与を持つのかを調べるには、単なる相関を越えて、構造状態が変化する条件で複製の可否や時間順序がどう変わるかを追う必要がある。

そこで次に必要になるのが、A/B コンパートメントの形成に異常を生じさせる介入と、その後の DNA 複製を同じ時間軸で観測する実験である。構造形成が正常な軌道から外れたときに DNA 合成そのものが進まなくなるのか、それとも合成は進むが複製される領域の順序だけが変わるのかを分けて測れば、構造が関係している範囲を絞り込める。

ただし、介入に使う分子が DNA 複製にも直接関与する場合、構造変化と複製異常が同時に起きただけでは、構造を原因と断定できない。介入、構造変化、複製開始、複製タイミングという複数の段階を時間的に分け、どの変化がどの結果に先行するかを見る必要がある。次に扱う GINS4 の実験は、この制約を残しながらも、染色体の 3 次元構造と DNA 複製を単なる対応関係より一段強い形で結び付ける材料を与えている。

3. 染色体の 3 次元構造形成が乱れると、複製も乱れる

2026 年に理化学研究所と京都大学などの研究チームが Nature Communications に発表した研究は、A/B コンパートメントの形成と DNA 複製を同じ細胞周期の中で追跡し、両者の関係を実験的に調べた[6]。研究対象はマウス胚性幹細胞であり、細胞分裂後の G1 期から DNA 複製が始まる S 期へ移る過程で、染色体の 3 次元構造がどのように再形成されるか、その変化が DNA 合成や複製タイミングとどう対応するかを解析した。理化学研究所の公式発表でも、DNA 配列上の情報だけでなく、核内で形成される 3 次元構造が複製制御に関与することが研究の主眼として説明されている[7]。

研究の鍵になったのは GINS4 である。GINS4 は DNA 複製ヘリカーゼ複合体の構成要素であり、本来は DNA 複製との関係が強いタンパク質である。ところが、研究チームがゲノム規模の CRISPR スクリーニングによって A/B コンパートメント形成に必要な因子を調べたところ、GINS4 が染色体の空間的な区画化にも関与することが見つかった[6]。一つの分子が、複製装置の構成要素であると同時に、複製開始前の染色体構造形成にも関係していたことになる。

GINS4 を減少させても、染色体構造は最初から全面的に崩れたわけではない。細胞分裂後の初期には、対照群と似た構造形成が進んだ。しかし時間が経過すると、通常とは異なる変化が現れた。A コンパートメントと B コンパートメントの間の A-B 接触が増え、B コンパートメント内部の B-B 接触は減少した。その結果、対照細胞がたどる構造変化とは異なる軌跡へ分岐した[6]。

この時間差は、構造異常を一枚の静止画像として扱えないことを示している。GINS4 が減少した直後に正常構造が消失するのではなく、初期には正常群と似た状態を通り、その後の構造更新が正常な経路から外れていく。異常は「構造があるか、ないか」ではなく、細胞周期が進むにつれて本来到達すべき構造へ到達できなくなる形で現れている。

同じ過程で DNA 複製にも異常が生じた。A/B コンパートメントが十分に組織化されていない細胞では、S 期へ進む分子的な条件が整っていても DNA 合成が効率よく進まなかった[6]。ここでは二つの段階を分けて考える必要がある。細胞が S 期へ進むことと、ゲノム全体で効率的な DNA 合成を実行できることは同じではない。細胞周期上は複製可能な段階へ到達しても、染色体側の構造状態が十分に整っていなければ、複製装置が通常どおり働くとは限らない。

さらに、構造形成が完全には戻っていなくても DNA 合成自体は起こる状態が観測された。この場合には、複製できるかどうかだけでなく、どの領域をいつ複製するかという複製タイミングにも異常が残った[6]。大きな構造異常では DNA 合成そのものが阻害され、構造が部分的に回復した状態では DNA 合成が進んでも時間的な順序が正常化しない。この差は、染色体構造が単なる複製開始の許可条件としてだけ働いているわけではない可能性を示す。

状態	染色体構造	DNA 複製	読み取れること
通常状態	A/B コンパートメントが細胞周期に沿って再形成される。	S 期で DNA 合成が進み、領域ごとの複製タイミングも維持される。	構造形成と複製制御が同じ時間過程の中で成立している。
GINS4 減少	初期には対照群と似ているが、その後 A-B 接触が増え、B-B 接触が減少して正常な構造軌跡から外れる。	S 期へ進む条件があっても、効率的な DNA 合成が阻害される。	正常な構造形成が失われる条件と、複製実行の障害が同じ時間過程で現れる。
部分回復	A/B 分離は回復方向へ進むが、正常状態には十分到達しない。	DNA 合成が起きても複製タイミングに異常が残る。	構造状態は複製の可否だけでなく、領域ごとの時間的な秩序とも対応する。

この結果は、「A コンパートメントは早く、B コンパートメントは遅く複製される」という既知の相関より一段先へ進んでいる。GINS4 を減少させると、まず染色体が正常とは異なる構造変化の軌跡へ入り、その条件で DNA 合成の効率や複製タイミングにも異常が現れる。空間構造と複製は、単に別々の測定値が対応しているだけではなく、同じ細胞周期の進行過程の中で機能的に結び付いている。

ただし、この実験だけから「3 次元構造の変化が複製異常を直接引き起こした」と断定することはできない。GINS4 自体が DNA 複製ヘリカーゼ複合体の構成要素だからである。GINS4 の減少によって構造形成が乱れ、その構造異常を介して複製が乱れた可能性と、GINS4 の減少が複製装置へ直接影響した可能性は、完全には分離されていない。この制約を残した上でも、構造形成の程度と DNA 合成、さらに複製タイミングの異常が同じ実験系で段階的に対応したことには意味がある。

そこで本稿では、この結果を「構造だけが複製を決める」とは読まない。より限定して、正常な染色体構造の形成が、DNA 複製を成立させる複数の条件の一つとして働いている可能性が高まったと捉える。しかも、その条件は単純な正常、異常の二値ではない。構造が大きく外れれば DNA 合成自体が進みにくくなり、部分的に戻れば合成できても時間的な順序が残る。この連続性が、次に構造を静的な配置ではなく、時間の中で変化する状態として扱う理由になる。

4. 構造が戻りきらないと、実行順序がずれる

この研究で注目すべきなのは、染色体構造を「正常」と「異常」の二値で分けるだけでは、DNA 複製の結果を説明しきれない点である。A/B コンパートメントの分離が完全ではないものの、DNA 合成そのものは開始できた細胞では、通常なら S 期の後半に複製される領域の一部が早い時期に複製されるなど、複製タイミングに異常が生じた[6]。構造が一定水準まで回復すれば複製開始は可能になるが、それだけでは正常な時間順序まで保証されなかった。

この結果から、構造状態と複製結果の関係には少なくとも二つの段階があると考えられる。構造の乱れが大きい段階では、DNA 合成そのものが効率よく進まない。一方、A/B コンパートメントが部分的に形成される段階まで戻ると、DNA 合成は可能になるが、どの領域をいつ複製するかという時間的な秩序には異常が残る。複製を開始できる条件と、複製順序を正常に保つ条件は完全には同じではない。

構造状態	DNA 合成	複製タイミング	読み取れること
A/B コンパートメントが正常に形成される	効率よく進む。	領域ごとの通常の時間秩序が維持される。	複製の実行と時間順序の双方が成立する。
構造の乱れが大きい	効率よく進みにくい。	正常な順序を評価できる段階まで複製が進みにくい。	構造形成の不全と複製実行の障害が対応する。
A/B 分離が部分的に回復する	開始できる。	後期複製領域の一部が早期に複製されるなど、異常が残る。	複製開始に十分な状態と、時間秩序の正常化に十分な状態は一致しない。

この差があるため、染色体構造の役割を単純な「開始スイッチ」として表すのは不十分である。もし構造が複製開始の可否だけを決めるのであれば、DNA 合成を開始できた時点で、その後の複製順序まで正常に戻ることが期待される。しかし実際には、DNA 合成が起きても後期複製領域の一部が早期へずれる。構造状態は、実行を許可するかどうかだけでなく、ゲノムの各領域をどの時点で処理可能にするかという時間的な制約にも関係している。

ここで因果関係を一段分ける必要がある。第一に、A/B コンパートメントの形成状態が変わると、各ゲノム領域が置かれる核内の関係状態も変わる。第二に、その関係状態が十分に成立していないと、DNA 合成そのもの、あるいは領域ごとの複製時刻が正常な状態から外れる。観測されているのは、単に構造と複製結果が同時に違うという対応ではなく、構造の成立度合いに応じて複製異常の現れ方も変わるという段階差である。

ただし、この段階差を構造だけの効果として確定することはできない。GINS4 自体が DNA 複製ヘリカーゼ複合体の構成要素であり、その量や機能の回復が DNA 合成へ直接影響した可能性が残るからである。構造の回復度と GINS4 の機能回復は、今回の実験だけでは完全には分離されていない。ここで確実に言えるのは、A/B コンパートメントの形成状態と複製タイミングの正常化が同時には回復せず、構造が不完全な状態では時間秩序にも異常が残ったということである。

この制約を残しても、本稿の中心命題は前節より具体的になる。DNA 配列が複製対象の内容を規定していても、その内容がいつ処理されるかは配列だけでは決まらない。染色体構造と細胞状態が異なれば、複製できるかどうかだけでなく、どの領域が先に処理され、どの領域が後になるかまで変わり得る。構造は情報内容そのものを変更するのではなく、内容を時間的な実行順序へ接続する条件の一部として働く。

しかし、ここまでの説明にはまだ不足がある。「構造が崩れた」「部分的に戻った」という表現では、異なる時点の染色体を静止画のように比較しているからである。実際の染色体は、細胞周期を通じて一つの形を維持しているわけではない。分裂後には構造を組み直し、G1 期から S 期へ進むにつれて状態を変化させる。正常な細胞でも、構造は時間とともに変わり続ける。

そのため、正常性を一つの完成形との一致として定義すると、細胞周期上の途中状態を正しく評価できない。ある構造が G1 初期では正常でも、S 期直前まで同じ状態に留まっていれば、今度は構造形成の遅れとして扱う必要がある。評価対象は「どの形になっているか」だけではなく、「どの時点でその状態に到達し、その後どこへ移るか」である。染色体構造を複製の実行条件として扱うなら、一時点の配置ではなく、時間の中でたどる構造変化まで記述する必要がある。

5. 染色体構造は固定物ではなく、毎周期作り直される

染色体構造を時間変化として扱う必要があるのは、正常な細胞でもその形が一定ではないからである。細胞分裂に入ると、染色体は高度に凝縮した分裂期構造を取る。分裂が終わって核が再形成されると、そこから間期の 3 次元構造が組み直される。A/B コンパートメント、局所的な接触領域、長距離ループなどは、分裂前の状態をそのまま保持しているのではなく、細胞周期の進行に沿って再形成される。

単一細胞 Hi-C を用いた研究では、この再形成が連続的な時間変化として観測されている[8]。染色体コンパートメント、局所的な接触領域、長距離ループは同じ速度で一斉に戻るのではなく、それぞれ異なる時間変化を示す。それでも、各細胞の染色体接触パターンには細胞周期上の位置に応じた規則性があり、接触構造の特徴だけから個々の細胞を細胞周期上に並べられるほど、構造変化には連続した秩序が存在した。

この結果は、正常な染色体構造を一つの完成形として定義できないことを意味する。分裂期の凝縮した構造は、G1 期の基準で見れば大きく異なるが、分裂期に存在する限り異常ではない。逆に、G1 初期としては正常な構造が S 期直前まで変化せずに残っていれば、それは正常な形であっても、正常な時刻には存在していない。構造の評価には、形そのものと細胞周期上の位置を同時に持つ必要がある。

細胞周期の段階	染色体構造	複製との関係	正常性の判断
分裂期	染色体が高度に凝縮し、間期の 3 次元構造は大きく再編される。	DNA 複製は進行しない。	間期構造から大きく外れていても、この時期には正常である。
G1 期	A/B コンパートメントや長距離接触が再形成される。	次の S 期で使われる複製タイミングのプログラムも成立していく。	構造が完成しているかではなく、再形成が適切に進んでいるかが重要になる。
S 期	形成された間期構造の中で各ゲノム領域の関係状態が維持される。	領域ごとの複製タイミングに沿って DNA 合成が進む。	G1 期で正常だった途中状態に留まっていれば、S 期では遅れになる。
次の分裂へ	間期構造から再び凝縮した分裂期構造へ移行する。	一周期の複製を終え、次の構造再編へ移る。	一つの形を保存するのではなく、周期全体を再現することが正常になる。

複製タイミングとの関係では、分裂後の G1 期が特に重要になる。2015 年の研究では、早い G1 期に長距離の染色体接触が再形成される時期と、複製タイミングのプログラムが成立する時期が重なることが示された[9]。S 期に入って DNA 合成が始まってから複製順序が初めて決まるのではない。その前段階で核内構造の再構築と、どの領域を早く複製し、どの領域を遅く複製するかという時間的な秩序の形成が並行して進んでいる。

ここでは二つの時間過程が接続している。第一に、分裂によって大きく変形した染色体が、G1 期を通じて間期の構造へ再構築される。第二に、その再構築と重なる時期に、次の S 期で使われる複製タイミングのプログラムが成立する。染色体構造と複製順序がこの段階でともに組み直されるなら、S 期の複製を説明するときに、S 期開始時点の構造だけを見るのでは不十分である。そこへ至るまでに、構造がどのように変化してきたかも状態記述へ含める必要がある。

この観点を入れると、「構造が複製を制御する」という表現も限定して捉え直せる。細胞内に一つの完成済み構造が存在し、それが固定された制御装置として S 期を支配しているわけではない。分裂後に染色体構造が再構築され、その時間発展の先で複製可能な状態へ到達する。構造が複製条件に関与するとすれば、重要なのは一時点の形だけではなく、細胞周期のどの位置でどの構造状態へ到達したかである。

この違いは、GINS4 を減少させた細胞の結果ともつながる。GINS4 減少細胞は最初から正常細胞と全く異なる構造を示したわけではなく、初期には似た状態を通過した後で、時間とともに別の構造軌跡へ外れた。異常は初期状態の違いだけではなく、その後の更新過程の違いとして現れた。正常性を一時点の構造差だけで測れば、この分岐を十分には捉えられない。

このため、本稿でいう「安定」は、形が変化しないことを意味しない。分裂期には染色体を凝縮し、分裂後には G1 期の構造を再形成し、S 期には複製を進め、その後再び分裂へ向かう。この順序だった変化を周期ごとに再現できることが、細胞周期における構造の安定性である。静止点として同じ形へ留まり続けるのではなく、時間とともに状態を変えながら同じ種類の周期的な経路へ戻ることが安定になる。

したがって、染色体構造を評価する対象も一時点の形から時間的な経路へ移る。一つの「正常構造」からの距離だけを測るのではなく、正常な細胞が細胞周期を通じてたどる一連の構造変化と比較し、現在の細胞がそのどこにあり、正常軌道からの偏差がどう変化しているかを見る必要がある。次節では、この時間的な正常性を、固定した構造ではなく周期軌道として数理的に記述する。

6. 構造振動モデルでは、細胞周期を周期軌道として書ける

前節までの議論では、染色体構造の正常性を一時点の形ではなく、細胞周期に沿った時間変化として評価する必要があることを確認した。この時間変化は、既稿で定義した構造振動モデルへ接続できる。構造振動モデルでは、時刻 \(t\) の構造を \(S_t\) とし、その構造を計算可能な特徴量へ写像したものを \(x_t=\phi(S_t)\) と置く。染色体構造であれば、\(\phi\) の成分には A-A 接触、B-B 接触、A-B 接触、長距離接触、コンパートメント分離の強さなどを含められる。

隣接する二時点の構造差は、構造振幅として次のように定義できる[10]。

\[
A_t=d(S_t,S_{t-1})
=\left\|\phi(S_t)-\phi(S_{t-1})\right\|_2
\]

\(S_t\) は時刻 \(t\) の染色体構造、\(\phi\) は構造を数値ベクトルへ変換する写像、\(A_t\) は一つ前の時点から構造がどれだけ変化したかを表す。この量が大きいこと自体は異常を意味しない。分裂期から G1 期へ移るときには、染色体の凝縮状態から間期の 3 次元構造へ大きく再編されるため、正常な細胞でも構造差は大きくなり得る。

量	表すもの	大きい場合の意味	それだけでは判断できないこと
\(S_t\)	時刻 \(t\) の染色体構造そのもの。	大きさという概念は持たない。	細胞周期上のどの段階に対応する構造か。
\(x_t=\phi(S_t)\)	染色体構造を比較可能な特徴量へ変換した状態。	各成分の値によって構造的特徴の強弱が表れる。	その状態が現在の細胞周期に対して正常かどうか。
\(A_t\)	一つ前の時点からどれだけ構造が変化したか。	構造更新が大きかったことを示す。	正常方向へ進んだのか、異常方向へ外れたのか。

この区別は重要である。構造振幅 \(A_t\) は変化量を測るが、その変化が適切かどうかまでは判定しない。分裂期から G1 期へ大きく構造が変化することは正常である一方、G1 期の途中で本来形成されるべき A/B コンパートメントが崩れる方向へ大きく動くなら異常である。同じ「大きな変化」でも、細胞周期上の位置と変化の方向によって意味が異なる。

構造振動モデルでは、このような系の安定状態を静止点だけに限定していない。状態が時間とともに変化していても、一定周期で近い経路を繰り返すなら、その周期軌道を安定状態として扱える[10]。細胞周期はこの考え方と整合する。染色体構造は分裂、G1、S、G2 を通じて連続的に変化するが、正常な細胞では各周期で同じ種類の構造変化が再現される。

この周期性を表すため、細胞周期上の位置を位相 \(\theta_t\) と書く。最小モデルでは、位相が一定速度 \(\omega\) で進むとして、次のように置ける。

\[
\theta_{t+1}
=
(\theta_t+\omega)\bmod 2\pi
\]

\(\theta_t\) は細胞周期のどの位置にいるかを表し、\(\omega\) は一刻みでどれだけ周期が進むかを表す。位相が \(2\pi\) に達すると 0 に戻るため、細胞周期を一周する運動として記述できる。分裂期、G1 期、S 期、G2 期は、この円周上の異なる領域に対応する。

位相上の位置	対応する細胞状態	期待される染色体構造
分裂期付近	染色体分配の段階。	高度に凝縮した分裂期構造。
G1 期	分裂後に間期構造を再形成する段階。	A/B コンパートメントや長距離接触が再構築されていく。
S 期	DNA 複製が進行する段階。	複製タイミングと対応した間期構造が形成されている。
G2 期	複製後から次の分裂へ向かう段階。	間期構造を保ちながら次の分裂期構造へ移る準備をする。

正常な細胞について、位相 \(\theta\) ごとに期待される染色体構造を \(S^\ast(\theta)\) と置く。その計算表現は次のようになる。

\[
x^\ast(\theta)=\phi(S^\ast(\theta))
\]

ここで \(S^\ast(\theta)\) は一つの固定した「正常構造」ではない。分裂期には分裂期に対応する正常構造があり、G1 期には再形成途中の正常構造があり、S 期には複製進行中の正常構造がある。それらを位相順に並べたものが一周する正常軌道になる。3 次元ゲノムが細胞周期を通じて解体と再形成を繰り返すことは、既存研究でも体系的に整理されている[11]。

この表現に変えると、正常性の基準も変わる。一つの完成した構造へ近いかどうかではなく、現在の位相 \(\theta_t\) に対して期待される \(S^\ast(\theta_t)\) に近いかどうかを評価することになる。同じ染色体構造でも、G1 初期であれば正常な途中状態として扱える一方、S 期直前に同じ状態へ留まっていれば、構造形成の遅れと解釈できる。

ここで、位相 \(\theta_t\) と構造 \(x_t\) を区別しておく必要がある。実データへ適用する場合、位相を同じ染色体構造だけから推定して、その構造を正常軌道と比較すると評価が循環する可能性がある。そのため、細胞周期マーカーや実験上の時間情報など、構造とは独立した情報から \(\theta_t\) を定め、その位相に対応する正常構造と比較する方がよい。

構造振幅 \(A_t\) と周期軌道 \(S^\ast(\theta)\) を分けて定義すると、次に必要な量も明確になる。\(A_t\) は「どれだけ動いたか」を表すが、現在の構造が正常な周期軌道からどれだけ外れているかは表さない。GINS4 減少細胞のように、初期には正常群と似た状態を通りながら、その後別の軌跡へ分岐する現象を記述するには、正常軌道そのものからの距離を別に測る必要がある。

7. 正常構造からの距離ではなく、正常軌道からの距離を測る

正常な染色体構造を位相ごとの周期軌道として定義すると、構造異常も一時点の正常、異常という二値ではなく、正常軌道からの偏差として連続的に扱える。実際の細胞で観測された構造を \(x_t=\phi(S_t)\)、同じ細胞周期位相で正常細胞に期待される構造を \(x^\ast(\theta_t)\) とし、両者の距離を \(D_t\) と置く。

\[
D_t
=
\left\|x_t-x^\ast(\theta_t)\right\|_Q
=
\sqrt{
\left(x_t-x^\ast(\theta_t)\right)^\top
Q
\left(x_t-x^\ast(\theta_t)\right)
}
\]

\(D_t\) は、時刻 \(t\) の実測構造が、その細胞周期位相で期待される正常構造からどれだけ外れているかを表す。\(Q\) は構造特徴の重みや特徴間の関係を表す正定値行列であり、構造振動モデルで定義した重み付き距離をそのまま利用できる[10]。特徴量ごとの重みだけを扱う場合は \(Q\) を正の対角行列とし、たとえば A/B コンパートメントの分離を重視するなら、それに対応する対角要素を大きくできる。

この定義で必要なのは、実測構造を単一の「正常構造」と比較しないことである。分裂期の構造を G1 期の正常構造と比較すれば、正常な細胞でも大きな距離が生じる。比較対象は常に、その細胞が現在いる位相 \(\theta_t\) に対応する \(x^\ast(\theta_t)\) でなければならない。正常性は形そのものではなく、現在の時刻に対して適切な形にあるかどうかとして評価される。

そのため、実データへ適用するときには位相 \(\theta_t\) の決め方にも制約がある。染色体構造そのものから位相を推定し、その同じ構造を正常軌道と比較すると、評価対象と基準が循環する可能性がある。細胞周期マーカー、実験開始からの経過時間、DNA 量など、構造とは独立した情報から位相を推定し、その位相に対応する正常構造と比較する方がよい。

今回の研究へ当てはめるなら、\(x_t\) の成分には A-A 接触、B-B 接触、A-B 接触、長距離接触、コンパートメント分離の強さなどを含められる。GINS4 枯渇細胞では、初期には対照群と似た構造変化を示した後、A-B 接触が増加し、B-B 接触が減少して、主成分分析上でも対照群とは異なる軌跡へ分岐した[6]。構造振動モデルでは、この分岐を、同じ位相の正常軌道に対する \(D_t\) が小さい状態から、より大きな状態へ移る現象として表現できる。

ただし、主成分分析上で軌跡が分かれたという観測だけから、\(D_t\) が時間とともに単調に増加したとまでは言えない。\(D_t\) の時間変化を実際に評価するには、正常軌道 \(x^\ast(\theta)\)、特徴量 \(\phi\)、重み \(Q\) を定義した上で各時点の距離を計算する必要がある。本稿で \(D_t\) を導入する目的は、論文の結果をそのまま数値へ読み替えることではなく、観測された軌道分岐を検証可能な量へ変換することである。

量	定義	測っているもの	単独では判断できないこと
\(A_t\)	\(\left\\|x_t-x_{t-1}\right\\|_2\)	一つ前の時点から構造がどれだけ変化したか。	その変化が正常方向か異常方向か。
\(D_t\)	\(\left\\|x_t-x^\ast(\theta_t)\right\\|_Q\)	現在の構造が、同じ細胞周期位相の正常軌道からどれだけ外れているか。	正常軌道へ近づいているのか、遠ざかっているのか。
\(\Delta D_t\)	\(D_t-D_{t-1}\)	正常軌道からの偏差が前時点より増えたか減ったか。	構造そのものがどれだけ大きく変化したか。

\(D_t\approx0\) であれば、実測構造はその位相で期待される正常構造に近い。\(D_t\) が大きければ、細胞周期は進んでいても染色体構造の再形成が追い付いていないか、正常とは異なる構造状態へ外れている可能性がある。これにより、「構造が戻った」「まだ弱い」「正常群とは異なる」といった定性的な表現を、細胞周期上の位置を考慮した偏差として扱える。

一方、\(D_t\) だけでは構造がどちらへ動いているかは分からない。そのため、隣接時点との差として

\[
\Delta D_t=D_t-D_{t-1}
\]

を導入する。\(\Delta D_t<0\) なら位相ごとの正常基準からの偏差が縮まり、\(\Delta D_t>0\) なら偏差が広がっている。これは構造ベクトルそのものの移動方向ではなく、正常軌道に対する偏差の増減を表す。

たとえば、分裂後の染色体が急速に再構築されている細胞では \(A_t\) が大きくても、同じ位相の正常軌道へ近づいているなら \(\Delta D_t<0\) となる。この場合、大きな構造変化は異常ではなく、正常な再形成である。逆に、構造がほとんど変化せず \(A_t\) が小さくても、正常軌道から離れた位置に留まり続けて \(D_t\) が大きいままなら、構造形成が停止した状態として扱う必要がある。

観測状態	\(A_t\)	\(D_t\)	\(\Delta D_t\)	解釈
大きく構造が変わりながら正常軌道へ戻る	大きい	減少する	負	正常な再構築または回復。
異常位置でほとんど変化しない	小さい	大きいまま	0 付近	構造形成の停滞。
正常軌道から別の方向へ急速に外れる	大きい	増加する	正	構造軌道からの逸脱。

この三つの量を分けることで、一時点の形だけでは区別できなかった状態を記述できる。\(A_t\) は変化の大きさ、\(D_t\) は現在位置の異常度、\(\Delta D_t\) は正常軌道へ近づいているか遠ざかっているかを表す。GINS4 枯渇細胞についても、「構造が正常か異常か」ではなく、正常細胞と似た軌道からいつ分岐し、その後どの程度離れていくかという時間過程として扱える。

この定式化によって、構造異常は固定した基準からの距離ではなく、時間依存の正常軌道からの偏差になる。次に必要なのは、この偏差が単なる構造上の違いにとどまらず、DNA 複製という実行過程へどのように影響するかをモデルへ組み込むことである。

8. 複製開始を構造状態に結合する

正常軌道からの構造偏差 \(D_t\) を定義すると、次に問えるのは、その偏差が DNA 複製へどう接続されるかである。構造が複製の実行条件の一部として働くなら、構造状態は単に観測対象として存在するだけではなく、各ゲノム領域が次の時点で複製へ進む確率へ影響するはずである。そこで、構造の時間発展と複製の時間発展を同じ状態モデルへ入れる。

ゲノムを複数の領域 \(i=1,\ldots,n\) に分け、領域 \(i\) が時刻 \(t\) までに複製済みである割合、またはその期待値を \(r_i(t)\) とする。多数の細胞をまとめて見る場合、\(r_i(t)\) はその領域を既に複製した細胞の割合として読める。個々の細胞に限れば、未複製なら 0、複製済みなら 1 に対応する。

\[
0\le r_i(t)\le1
\]

次に、まだ複製されていない領域 \(i\) が、時刻 \(t\) から次の時間刻みまでに複製へ進む確率を \(p_i(t)\) とする。本稿では、複製開始を決める最小モデルとして次の形を置く。

\[
p_i(t)
=
\sigma\left(
\alpha_i(\theta_t)
+\beta q_i(S_t)
-\lambda D_t
+\gamma^\top C_t
\right)
\]

\(\sigma\) は入力を 0 から 1 の範囲へ変換する関数であり、式全体を複製確率として読めるようにする。この式では、複製開始を一つの原因へ還元せず、領域固有の時間条件、局所構造、細胞全体の構造偏差、構造以外の分子状態を分離している。

項	表す状態	複製モデルでの役割	具体的な意味
\(\alpha_i(\theta_t)\)	領域 \(i\) 固有の時間プログラム	その位相で領域 \(i\) が複製されやすい基礎条件を与える。	早期複製領域と後期複製領域の違いを表す。
\(q_i(S_t)\)	領域 \(i\) の局所的な構造状態	その領域が現在どのような構造環境に置かれているかを反映する。	A/B コンパートメントへの所属や周辺領域との接触状態などを特徴量にできる。
\(D_t\)	細胞全体の正常軌道からの構造偏差	構造全体が正常な細胞周期軌道から外れた影響を表す。	GINS4 枯渇時のような、局所だけではない構造形成全体の異常を扱う。
\(C_t\)	構造以外の分子状態	構造だけでは説明できない複製条件をモデルへ残す。	GINS4、RIF1、PP1、CDK 活性などを含められる。

\(\alpha_i(\theta_t)\) と \(q_i(S_t)\) は似ているように見えるが、役割は異なる。\(\alpha_i(\theta_t)\) は、正常な細胞周期を進んだときに、その領域がいつ複製されやすいかという基準を表す。一方、\(q_i(S_t)\) は、その時点で実際に置かれている局所的な構造環境を表す。正常なら両者は整合する可能性が高いが、構造形成が乱れれば、本来は後期に複製される領域が通常とは異なる構造環境へ置かれ、時間プログラムと局所状態が食い違う場合を記述できる。

\(D_t\) はさらに異なる尺度を持つ。\(q_i(S_t)\) が領域ごとの局所状態を表すのに対し、\(D_t\) は細胞全体の染色体構造が正常軌道からどれだけ外れているかを表す。GINS4 枯渇細胞では A-B 接触の増加や B-B 接触の減少など、複数の構造特徴が同時に変化した[6]。そのような全体的な構造異常が複製効率を低下させるという仮説を、式では \(-\lambda D_t\) として表している。

\(\lambda>0\) とすれば、正常軌道から遠ざかるほど複製へ進む確率は低下する。ただし、これは研究から直接測定された法則ではなく、本稿が置く検証対象である。すべてのゲノム領域が構造偏差から同じ影響を受ける保証もない。実データでは、領域ごとに感受性が異なる可能性を考え、\(\lambda\) を \(\lambda_i\) へ拡張できる。

\(C_t\) を別に残すことも必要である。今回の実験で操作された GINS4 は、染色体構造の形成だけでなく DNA 複製ヘリカーゼ複合体にも直接関与する。そのため、GINS4 枯渇による複製低下をすべて \(D_t\) の効果として吸収すると、構造を介した影響と複製装置への直接作用を区別できなくなる。構造偏差 \(D_t\) と GINS4 などの分子状態 \(C_t\) を別の変数として置くことで、少なくともモデル上では二つの経路を分離できる。

この複製確率を使うと、領域 \(i\) の複製状態は、まだ複製されていない部分だけが次の刻みに進むとして次のように書ける。

\[
r_i(t+1)
=
r_i(t)
+
\left(1-r_i(t)\right)p_i(t)
\]

\(r_i(t)=1\) なら、その領域は既に複製済みなのでそれ以上増えない。\(r_i(t)<1\) なら、未複製部分 \(\left(1-r_i(t)\right)\) のうち \(p_i(t)\) に対応する割合が次の時点までに複製へ進む。この式は個々の細胞で DNA が連続的に少しずつ複製済みになることを意味するのではない。単一細胞では複製済みか未複製かという状態遷移を考え、その多数の細胞について期待値または割合を取った表現である。

この更新式には履歴も入る。一度複製された領域は再び未複製状態へ戻らないため、時刻 \(t+1\) の状態は現在の複製確率 \(p_i(t)\) だけではなく、それ以前にどこまで複製が進んだかを表す \(r_i(t)\) に依存する。同じ構造状態と分子状態に置かれていても、既に複製を終えた細胞と、まだ複製していない細胞では次の状態が異なる。

このモデルが表している因果関係は二段階に分けられる。まず、細胞周期の位相、局所的な染色体構造、正常軌道からの偏差、複製因子の状態が、その時点で領域 \(i\) を複製可能にする条件を作り、\(p_i(t)\) を決める。次に、その確率が既存の複製状態 \(r_i(t)\) に作用し、次の時点の \(r_i(t+1)\) を変える。構造が DNA の内容を書き換えるのではなく、現在の状態から次の実行状態へ移る確率に作用するという位置付けである。

なお、この式を書いただけでは、\(\beta\)、\(\lambda\)、\(\gamma\) を因果効果として識別できるわけではない。局所構造 \(q_i(S_t)\)、全体偏差 \(D_t\)、GINS4 などの分子状態 \(C_t\) は実際の細胞では互いに関連している可能性がある。どの項が独立に複製へ寄与しているかを決めるには、前節までに述べたように、介入条件や時間順序を利用して変数間の影響を分離する必要がある。

この最小モデルの目的は、DNA 複製の分子機構を一つの式で再現することではない。目的は、「同じ DNA 配列が存在していれば次の複製状態も決まる」という記述では不足することを明示することである。複製対象となる内容は DNA 配列によって規定されるが、その内容がいつ実行されるかは、領域固有の時間プログラム、局所構造、構造軌道全体の状態、複製因子、そしてそれまでの複製履歴に依存する。

この形まで書けば、染色体構造を「情報そのもの」と呼ぶ必要はない。構造が担うのは、保存された DNA 配列を変更することではなく、その配列のどの領域を現在処理可能にするかを制約することである。内容と構造の役割を分けたまま、構造状態を次の実行状態へ結び付けられる。

9. GINS4 枯渇は、構造軌道からの逸脱として読める

前節では、染色体構造と DNA 複製を同じ状態モデルへ入れ、構造偏差 \(D_t\) が複製確率へ影響し得る形を置いた。ここでは逆に、GINS4 枯渇という実験介入が染色体構造の時間発展をどう変えるかを、構造振動モデルの更新則として読む。

構造振動モデルでは、観測結果が次の更新へフィードバックする一般的な系を含め、時刻 \(t\) の構造 \(S_t\) が次の構造 \(S_{t+1}\) へ移る過程を次のように書く[10]。

\[
S_{t+1}
=
M(S_t,C_t,O_t,U_t)
\]

ただし、今回の GINS4 実験で Hi-C や DNA 合成量として得られる \(O_t\) は研究者側の観測であり、その観測値自体が細胞の次状態を更新するわけではない。そこで本稿への適用では観測を更新則から分離し、構造の時間発展を

\[
S_{t+1}
=
M(S_t,C_t,U_t)
\]

と書く。\(S_t\) は現在の構造、\(C_t\) は細胞周期や分子状態などの条件、\(U_t\) は実験的な介入、\(M\) はそれらを受けて次の構造状態を決める更新規則である。\(O_t\) は Hi-C や DNA 合成量など、更新された状態を外部から測る観測量として扱う。この式では、染色体構造を外から与えられた固定値として扱わず、現在の状態から次の状態へ逐次更新される変数として扱う。

記号	本稿での意味	GINS4 実験での対応
\(S_t\)	時刻 \(t\) の染色体構造	A-A、B-B、A-B 接触やコンパートメント分離を含む構造状態。
\(C_t\)	構造更新に影響する細胞内条件	細胞周期上の位置や、介入後に実際に成立した GINS4 の量・活性を含む分子状態。
\(O_t\)	実験によって得られる観測	Hi-C、主成分分析、DNA 合成量、複製タイミングなど。今回の適用では更新則への入力ではない。
\(U_t\)	系へ加えられる介入	GINS4 を急速に減少させる実験操作。
\(M\)	現在状態から次状態への更新規則	正常群と枯渇群で異なる軌跡を生じさせる時間発展。

今回の実験では、GINS4 を急速に減少させる操作を \(U_t=U_{\mathrm{KD}}\) という介入として表し、その結果として細胞内に成立した GINS4 の量や活性は \(C_t\) の一部として扱う。対照細胞では、分裂後の G1 期から S 期へ進むにつれて A/B コンパートメントが再形成され、構造は正常細胞に共通する時間的な経路を進む。一方、GINS4 枯渇細胞では、初期には対照群と似た構造を示した後、時間経過とともに A-B 接触が増え、B-B 接触が減少し、主成分分析上でも対照群とは異なる軌跡へ分岐した[6]。

この結果で区別すべきなのは、初期状態と更新規則である。GINS4 枯渇細胞が最初から全く別の構造を持っていたのであれば、異常は単なる初期条件の違いとして説明できる。しかし実際には、初期には対照群と似た状態を通過し、その後の時間発展で差が拡大した。少なくとも観測上は、「どこから始まったか」だけではなく、「その状態から次にどこへ進んだか」が群間差を作っている。

比較時点	対照細胞	GINS4 枯渇細胞	モデル上の解釈
分裂後の初期	G1 期の構造再形成を開始する。	対照群と似た構造変化を示す。	初期の \(D_t\) は必ずしも大きくない。
時間経過後	A/B コンパートメントの分離が正常な経路で進む。	A-B 接触が増え、B-B 接触が減少する。	正常軌道との構造差が拡大する。
9 時間以降	対照群に共通する軌跡を進む。	主成分分析上で異なる軌跡が明瞭になる。	正常軌道からの時間依存的な逸脱として記述できる。

前節で定義した正常軌道からの偏差 \(D_t\) を使えば、この分岐を量として記述できる。ただし、論文で直接測定されたのは主成分分析上の軌跡や接触状態の違いであり、本稿で定義した \(D_t\) そのものではない。そのため、単純に

\[
D_t\uparrow
\]

と実測結果として断定するのではなく、正常軌道 \(x^\ast(\theta)\)、特徴量 \(\phi\)、重み \(Q\) を定義した場合には、枯渇群で \(D_t\) が増加するかを検証できる、と読むのが正確である。

この区別は、構造振動モデルを論文結果の言い換えで終わらせないために必要になる。主成分分析上で二群が離れたことと、位相ごとの正常軌道からの距離 \(D_t\) が増加したことは同じではない。前者は論文で観測された事実であり、後者は本稿のモデルを実データへ適用したときに検証すべき量である。両者を分けることで、観測結果とモデル上の解釈を混同せずに済む。

それでも、GINS4 枯渇実験は軌道という見方に適している。論文では、対照群と枯渇群を一時点の Hi-C マップだけで比較するのではなく、複数時点を通じた構造変化として比較しており、枯渇群では 9 時間以降に対照群とは異なる軌跡が現れると報告されている[6]。構造異常が「ある時点で違っていた」という形ではなく、「途中までは似ていたが、その後の時間発展が分岐した」という形で観測されている。

このとき、異常は少なくとも三つの量へ分解できる。現在どれだけ構造が変化しているかは \(A_t\)、正常軌道からどれだけ離れているかは \(D_t\)、その距離が増えているか減っているかは \(\Delta D_t\) で表せる。GINS4 枯渇細胞で確認したいのは、単に構造変化が大きいかではない。対照群と同じ位相に置いたときに \(D_t\) が増え、その増加がどの時点から始まるかである。

この表現にすると、GINS4 枯渇を「構造を壊した結果」とだけ呼ぶよりも精密になる。介入によって最初の構造がただちに消失したのではなく、細胞周期に伴う構造更新の過程が途中から正常群と異なる方向へ進んだと捉えられるからである。作用点は、現在の構造そのものだけではなく、現在状態から次状態への遷移にもある。

ただし、ここでも GINS4 の作用を構造更新だけへ限定してはいけない。GINS4 は DNA 複製そのものにも関与するため、\(U_{\mathrm{KD}}\) は構造更新則 \(M\) だけでなく、前節で置いた分子状態 \(C_t\) や複製確率 \(p_i(t)\) にも直接影響し得る。本稿のモデルで GINS4 枯渇を扱うなら、「介入によって構造軌道が変わり、その結果として複製が変わる」という一本の経路ではなく、構造を介する経路と複製装置へ直接作用する経路を分けて残す必要がある。

その上で構造振動モデルが与える追加的な問いは明確である。GINS4 枯渇後に正常軌道からの偏差 \(D_t\) がいつ増え始めるのか、その変化が DNA 合成量や複製タイミングの異常より先に現れるのかを測ればよい。構造軌道の逸脱が複製異常に先行するなら、構造が実行条件として働くという仮説を強く支持する。逆に、複製異常が先に現れて構造変化が後から続くなら、構造を主因とするモデルは修正が必要になる。

GINS4 枯渇実験を軌道として読む価値は、この時間順序まで検証対象にできる点にある。正常と異常を一枚の構造図で比較するだけでは、どちらが先に変化したかは分からない。時間発展をモデルへ入れることで、構造差の有無から一歩進み、構造軌道の逸脱が複製異常を予測できるかという因果に近い問いへ変えられる。

10. 構造のずれが大きいほど、時間秩序は乱れるのか

数理モデル化の価値は、既に観測された現象を記号へ置き換えることではない。モデルから、まだ直接測定されていない関係を導き、その関係が成立するかどうかを実験で判定できることにある。ここまでのモデルから最も直接的に導かれる予測は、正常な構造軌道からの偏差が大きい細胞ほど、DNA 複製の時間的な制御も正常状態から外れるというものである。

まず、複製側のずれを量として定義する。ゲノム領域 \(i\) の正常な複製時刻を \(T_i^\ast\)、観測された複製時刻を \(T_i\) とし、全領域について平均した複製時刻の誤差を \(E\) と置く。

\[
E
=
\frac{1}{n}
\sum_{i=1}^{n}
\left|T_i-T_i^\ast\right|
\]

\(E\) が小さければ、各領域は正常時に近い時刻で複製されている。\(E\) が大きければ、本来より早く複製される領域や、逆に遅れて複製される領域が増えていることを表す。ただし、この量が直接測っているのは複製時刻の誤差であり、領域間の順位だけではない。すべての領域が同じ時間だけ遅れた場合にも \(E\) は増えるため、「複製順序の入れ替わり」と「複製時刻全体のずれ」を厳密に分離したい場合には、順位相関など別の指標を追加する必要がある。

構造側では、前節までに定義した \(D_t\) を使う。ただし、複製時刻の誤差 \(E\) と比較するには、どの時点の \(D_t\) を使うかを固定しなければならない。たとえば G1/S 移行付近の構造状態がその後の複製タイミングを制約するという仮説を検証するなら、その位相で測った構造偏差を \(D_{\mathrm{G1/S}}\) とする。

このとき、最も単純な中心仮説は次のように書ける。

\[
E
=
f(D_{\mathrm{G1/S}}),
\qquad
\frac{\partial E}{\partial D_{\mathrm{G1/S}}}
>0
\]

意味するところは単純である。G1/S 移行時に正常な構造軌道から大きく外れているほど、その後に観測される複製時刻の誤差も大きくなる、という予測である。構造が複製の実行条件の一部であるなら、構造偏差は単なる形態上の違いではなく、後続する複製過程の異常量と対応するはずである。

今回の GINS4 研究は、この予測と整合する観察を含んでいる。A/B コンパートメントの形成が不十分な状態で DNA 合成が起きた細胞では、通常なら後期に複製される領域の一部が早期へ移るなど、複製タイミングに異常が残った[6]。構造形成の不完全さと複製時刻の異常が同じ実験条件で対応しているため、本稿のモデルが予測する方向とは一致する。

しかし、この結果だけから \(D_{\mathrm{G1/S}}\) と \(E\) の因果関係が成立したとは言えない。第一に、論文では本稿で定義した \(D_{\mathrm{G1/S}}\) を直接計算していない。第二に、GINS4 は構造形成だけでなく DNA 複製そのものにも関与するため、GINS4 の量や機能状態が構造偏差と複製誤差の双方を同時に変化させた可能性がある。構造と複製の相関だけでは、この共通原因を除外できない。

この代替説明を分けるには、構造以外の条件をできるだけそろえた上で、\(D_{\mathrm{G1/S}}\) と \(E\) の関係を測る必要がある。たとえば、細胞周期位相と GINS4 の量が近い細胞群を比較し、その中で A/B コンパートメントの形成度だけが異なる条件を作る。その条件でも構造偏差が大きい細胞ほど \(E\) が増えるなら、GINS4 の直接作用だけでは説明しにくくなり、構造軌道そのものの寄与をより強く支持できる。

逆に、GINS4 の量などをそろえると \(D_{\mathrm{G1/S}}\) と \(E\) の関係が消えるのであれば、構造偏差は複製異常の原因ではなく、同じ分子異常から生じる付随的な指標だった可能性が高くなる。その場合には、前節で置いた \(-\lambda D_t\) の寄与を小さくするか削除し、GINS4 などを含む \(C_t\) の直接効果を大きくしたモデルへ修正する必要がある。

もう一つの検証軸が時間順序である。構造が後続する複製の実行条件として働くなら、正常軌道からの逸脱は、少なくとも対応する複製タイミング異常より先に観測される必要がある。G1 期で \(D_t\) が増加し、その後の S 期で複製時刻のずれが現れるなら、「構造変化が先、複製異常が後」というモデルと整合する。

逆に、DNA 複製の異常が先に現れ、その後に A/B コンパートメントの分離が崩れるのであれば、因果方向は逆である可能性がある。その場合、構造は複製を制御する状態というより、複製異常の結果として変化する状態として扱う必要がある。一時点の相関では区別できないこの差を、時系列モデルでは変化の順序として検証できる。

検証する仮説	支持する結果	反証または修正につながる結果	必要な統制
構造偏差が大きいほど複製時刻の誤差も大きい	\(D_{\mathrm{G1/S}}\) が大きい細胞ほど \(E\) が一貫して大きい。	\(D_{\mathrm{G1/S}}\) と \(E\) に再現性のある関係がない。	細胞周期位相、GINS4 量、その他の複製因子をできるだけそろえる。
構造軌道の逸脱が複製異常に先行する	G1 期の \(D_t\) 増加が、その後の複製タイミング異常より先に現れる。	複製異常が先行し、その後に \(D_t\) が増加する。	構造と複製を十分な時間分解能で同じ条件から追跡する。
構造の回復に伴って複製時刻も正常化する	\(D_t\) の低下に続いて \(E\) も減少する。	構造が正常軌道へ戻っても \(E\) が改善しない。	構造回復と GINS4 自体の回復を可能な限り分離する。
構造には GINS4 の直接作用とは別の寄与がある	GINS4 などの分子状態を統制しても \(D\) と \(E\) の関係が残る。	分子状態を統制すると \(D\) の説明力が消える。	\(C_t\) と構造指標を別々に測定する。

さらに、モデルは「構造が悪いほど複製も悪い」という単純な単調関係だけを要求するわけではない。第4章で見たように、構造異常が大きい場合には DNA 合成そのものが進みにくく、部分的に回復した場合には DNA 合成は可能でも複製タイミングに異常が残る。この二段階が実際に存在するなら、\(D\) と複製結果の関係には、複製不能となる領域と、複製可能だが時間秩序が乱れる領域が存在する可能性がある。

その場合、\(E=f(D)\) は最初の近似にすぎない。大きな \(D\) では複製自体が十分に進まず \(T_i\) を定義しにくくなるため、複製量と複製時刻を別々の目的変数として扱う必要がある。構造偏差がある閾値を超えると DNA 合成効率が低下し、それより小さい偏差の範囲では複製可能であっても時刻誤差 \(E\) が増える、という二段階の関係も検証対象になる。

ここまで進めると、本稿の中心命題は観察の説明から反証可能な仮説へ変わる。「染色体構造と複製タイミングは関係している」という相関だけではなく、どの位相で構造偏差を測り、その偏差がどの複製指標へ先行し、どの条件を統制しても関係が残るかを指定できる。逆の結果が得られた場合には、構造の寄与を弱め、GINS4 など別の分子状態を中心とするモデルへ組み直すこともできる。

数理モデルの役割は、構造の重要性を証明済みの前提として置くことではない。構造が実行条件であるという主張を、測定可能な量と時間順序へ分解し、どの結果ならその主張を捨てるべきかまで明示することである。

11. 構造は、情報を時間へ展開する状態である

ここまでの数理化を情報という言葉へ戻すと、染色体構造の役割をより限定して記述できる。構造に DNA 配列とは別の文章が書かれていると考える必要はない。A/B コンパートメントに、塩基配列と同じ形式で読み出せる「第二の暗号」が保存されているわけでもない。DNA 配列と染色体構造は、同じ種類の情報を二重に持っているのではなく、細胞内で異なる役割を担っている。

DNA 配列が規定するのは、複製される対象の内容である。一方、前節までで扱った染色体構造は、その内容が現在どのような関係状態に置かれているかを表す。さらに細胞周期や複製因子の状態によって、同じゲノム領域でも複製可能になる時刻が変わる。内容が存在することと、その内容を現在実行できることは別の状態として記述する必要がある。

DNA 配列を \(I\)、時刻 \(t\) の構造を \(S_t\)、その他の細胞状態を \(C_t\)、実際の処理状態を \(R_t\) とすれば、最小限の関係は次のように書ける。

\[
R_{t+1}
=
F(I,S_t,C_t,R_t)
\]

\(I\) は複製される対象の内容、\(S_t\) はその内容が置かれた構造的な関係状態、\(C_t\) は細胞周期、GINS4、RIF1、PP1、CDK 活性などの条件、\(R_t\) は既にどこまで処理が進んだかを表す。次の実行状態 \(R_{t+1}\) は、DNA 配列だけから決まるのではなく、現在の構造、分子状態、それまでの処理履歴を含めて決まる。

状態	保持しているもの	この状態を捨てると失われる説明
\(I\)	DNA 配列という複製対象の内容	何を複製するのかを特定できなくなる。
\(S_t\)	ゲノム領域同士の接触やコンパートメントなど、その時点の構造状態	同じ配列が現在どのような核内環境に置かれているかを区別できなくなる。
\(C_t\)	細胞周期や複製因子など、構造以外の実行条件	構造が同じでも複製結果が異なる条件を説明できなくなる。
\(R_t\)	それまでにどこまで処理が進んだかという履歴	既に複製済みの領域と未複製の領域を区別できなくなる。

この分解には、状態記述の不足を明示できるという意味がある。DNA 配列 \(I\) だけを残して核内構造 \(S_t\) を捨てれば、複製される対象は保存されるが、その領域が現在どのような構造環境に置かれているかは失われる。前節のモデルで言えば、局所構造 \(q_i(S_t)\) や正常軌道からの偏差 \(D_t\) を計算できなくなり、複製時刻の違いを説明する変数の一部を落とすことになる。

構造 \(S_t\) を残しても、細胞周期上の位相を捨てれば別の問題が生じる。同じ構造が G1 初期に現れているなら正常な再形成過程かもしれないが、S 期直前まで同じ位置に留まっていれば構造形成の遅れかもしれない。一時点の構造だけでは、その状態が正常軌道上の途中なのか、正常な進行から取り残された状態なのかを区別できない。

さらに、現在の構造と細胞状態が同じでも、それまでの複製履歴 \(R_t\) が違えば次の状態は同じにならない。既に複製された領域は、同じ条件が続いてももう一度未複製状態から開始しない。第8章の更新式が示すように、現在の入力だけではなく、既にどこまで処理されたかが次状態を制約する。

既稿「単純な法則の中には、過去の影響が残っている」では、観測対象を少数の変数へ縮約すると、捨てた高次自由度の影響が後の時間発展に記憶として現れる構造を扱った[12]。本稿で扱う染色体構造と同じ現象ではない。しかし、状態を少なすぎる変数で表したときに、本来は次状態の決定に必要だった条件が見えなくなるという点では共通している。

その意味で、本稿の主張は「構造そのものが新しい情報内容を保存している」というものではない。DNA 配列だけでは不足する状態変数を、構造、細胞周期、分子状態、履歴へ分解している。情報内容と実行条件を同じものとして扱わないことが、このモデルの要点である。

この観点では、構造は情報の単なる「置き場所」ではないが、固定した「展開規則」そのものでもない。より正確には、現在どの展開が可能かを制約する状態である。同じ DNA 配列が保存されていても、構造軌道が異なれば、どの領域が複製可能になるか、どの領域が後へ回るか、DNA 合成自体が進めるかまで変わり得る。

静的な内容 \(I\) は、それだけでは時間順序を持たない。\(S_t\)、\(C_t\)、\(R_t\) が時間とともに更新されることで、保存された内容のどの部分を次に処理するかが逐次決まる。染色体構造を「情報を時間へ展開する状態」と呼ぶのは、この内容と実行の間を結ぶ役割を指している。

12. 情報の働き方は、構造の軌道にも依存する

本稿の出発点は、同じゲノムを持つ細胞でも、ゲノム領域ごとに複製される時刻が異なるという事実だった。DNA 配列は複製対象を規定するが、その配列だけでは、現在どの領域を先に処理し、どの領域を後へ回すかという細胞周期上の実行状態を記述できない。そこで、配列とは別に存在する状態として、核内の 3 次元構造を取り上げた。

A/B コンパートメントと複製タイミングには強い対応があり、GINS4 を減少させた実験では、A/B コンパートメント形成の異常とともに DNA 合成効率や複製タイミングにも異常が現れた[6]。ただし、GINS4 自体が DNA 複製へ直接関与するため、構造だけを原因とすることはできない。この制約をモデルにも残し、構造状態 \(S_t\) と GINS4 などの分子状態 \(C_t\) を分離した。

さらに、染色体構造は固定された制御装置ではない。分裂期には高度に凝縮し、分裂後の G1 期には A/B コンパートメントや長距離接触を再形成し、その後 S 期へ進む。正常な細胞でも構造は大きく変化しており、細胞周期に沿った連続的な変化として観測できる[8]。G1 期に核内構造が再形成される時期と、複製タイミングのプログラムが成立する時期も重なっている[9]。

このため、正常性を一つの完成した染色体構造として定義することはできない。G1 初期として正常な構造でも、S 期直前まで同じ状態に留まっていれば異常になり得る。正常とは特定の形を維持することではなく、細胞周期の位相に応じた構造へ移り続け、周期ごとに近い時間的な経路を再現することである。

構造振動モデルでは、この違いを周期軌道として表した。隣接時点でどれだけ構造が変わったかを \(A_t\)、同じ細胞周期位相の正常軌道からどれだけ外れているかを \(D_t\)、その距離が増減しているかを \(\Delta D_t\) とした。これにより、大きく構造が変化している正常な再形成と、ほとんど変化せず異常位置へ留まっている停滞を区別できる。

問い	必要な量	本稿での判断
何が保存されているか	DNA 配列 \(I\)	複製対象となる内容を規定する。
現在どの構造にあるか	構造状態 \(S_t\)	ゲノム領域が置かれた関係状態を規定する。
その構造が現在の時刻に適切か	位相 \(\theta_t\) と偏差 \(D_t\)	一時点の形ではなく、正常軌道との対応で判断する。
どの領域が次に複製されるか	\(p_i(t)\)	時間プログラム、局所構造、全体偏差、分子状態から決まると仮定する。
構造が本当に複製へ寄与するか	\(D\) と複製誤差 \(E\) の関係	分子状態を統制しても関係が残り、構造偏差が時間的に先行するかで検証する。

このモデルからは、単なる説明ではなく検証可能な予測が出る。G1/S 移行時の正常軌道からの構造偏差 \(D_{\mathrm{G1/S}}\) が大きいほど、その後の複製時刻の誤差 \(E\) も大きくなるという予測である。さらに、構造が後続する複製の実行条件なら、\(D_t\) の増加は対応する複製異常より先に現れる必要がある。

この予測が、GINS4 の量や細胞周期位相などをそろえても成立するなら、染色体の 3 次元構造は複製結果に付随して変化するだけではなく、時間的な実行秩序を作る状態の一部として位置づけられる。逆に、分子状態を統制すると \(D\) と \(E\) の関係が消える、あるいは複製異常が常に構造偏差より先行するなら、構造の寄与を小さく見積もり、GINS4 などの直接作用を中心とするモデルへ修正しなければならない。

この反証可能性が、本稿で構造を情報へ結び付ける際の境界になる。「構造も情報である」と広く言い換えるだけでは、何を観測すればその主張が誤りになるのか分からない。本稿では、DNA 配列という内容、構造状態、細胞周期位相、分子条件、複製履歴を分け、それらが次の実行状態へどう作用するかを変数として置いた。その結果、構造の役割を、内容そのものではなく実行条件の一部として限定できる。

中心命題は、この限定を含めて表す必要がある。情報の働き方は、保存されている内容だけでは決まらず、その内容を含む構造が時間の中でどの軌道をたどっているかにも依存する。DNA 配列は複製される内容を保存する。染色体構造は、その内容が現在どの関係状態に置かれているかを規定する。細胞周期に沿った構造の更新は、その状態を時間とともに変化させ、各領域を異なる時刻の実行へ接続する。

情報は内容として保存され、構造の時間変化を通じて実行順序へ展開される。

参考文献

id774, 意味は差異の読み取りから生まれる（2026-05-09）. https://blog.id774.net/entry/2026/05/09/4740/
Marchal, C., Sima, J., Gilbert, D. M., Control of DNA replication timing in the 3D genome. Nature Reviews Molecular Cell Biology 20, 721–737 (2019). https://doi.org/10.1038/s41580-019-0162-y
id774, 境界が生命の形を作る（2026-08-12）. https://blog.id774.net/entry/2026/08/12/4961/
Ryba, T. et al., Evolutionarily conserved replication timing profiles predict long-range chromatin interactions and distinguish closely related cell types. Genome Research 20, 761–770 (2010). https://doi.org/10.1101/gr.099655.109
Gnan, S. et al., Nuclear organisation and replication timing are coupled through RIF1–PP1 interaction. Nature Communications 12, 2910 (2021). https://doi.org/10.1038/s41467-021-22899-2
Oji, A., Yusa, K., Noda, I., Ichinose, T., Kondo, Y., Hiratani, I., Nuclear compartmentalization at the G1/S transition plays a key role in DNA replication control. Nature Communications 17, 6961 (2026). https://doi.org/10.1038/s41467-026-75264-6
理化学研究所, 細胞分裂後の染色体の形が DNA 複製を左右する（2026-08-10）. https://www.riken.jp/press/2026/20260810_1/index.html
Nagano, T. et al., Cell-cycle dynamics of chromosomal organization at single-cell resolution. Nature 547, 61–67 (2017). https://doi.org/10.1038/nature23001
Dileep, V., Ay, F., Sima, J., Vera, D. L., Noble, W. S., Gilbert, D. M., Topologically associating domains and their long-range contacts are established during early G1 coincident with the establishment of the replication-timing program. Genome Research 25, 1104–1113 (2015). https://doi.org/10.1101/gr.183699.114
id774, 構造振動モデルを数理モデルとして定義する（2026-04-05）. https://blog.id774.net/entry/2026/04/05/4318/
Miura, H., Hiratani, I., Cell cycle dynamics and developmental dynamics of the 3D genome: toward linking the two timescales. Current Opinion in Genetics & Development 73, 101898 (2022). https://doi.org/10.1016/j.gde.2021.101898
id774, 単純な法則の中には、過去の影響が残っている（2026-08-13）. https://blog.id774.net/entry/2026/08/13/5493/