Morphogenesis driven by nematic defects in active biological networks

技術要約：活性生物学的ネットワークにおけるネマチック欠陥による形態形成

問題提起
本論文は、細胞の形態形成、具体的には生物学的組織がいかにして無秩序な状態から特定の形状を持つ秩序構造へと遷移するかという、根本的な課題に取り組んでいる。細胞骨格のフィラメントが部分的な配向秩序を持つネマチック相へと組織化されることは確立されているが、これらの相におけるトポロジカル欠陥が、組織の変形やパターン形成を駆動する上での正確な機械的役割については依然として不明である。既存のモデルは、線形近似や定義済みのモルフォゲン勾配に依存することが多く、これではHydra（ヒドラ）のような再生能力を持つ生物に見られる、ネマチック秩序、弾性、およびストレス駆動型の成長の間の複雑で非線形な相互作用を捉えきれない可能性がある。著者らは、成長パターンをあらかじめ規定することなく、ネマチック秩序、弾性、および適応的緩和の結合から形状変化が創発される「活性生物学的ネットワーク（ABN）」として生体組織を扱う、連続体フレームワークの開発を目指している。

手法
著者らは、横断結合されたフィラメント鎖から構成され、ネマチック相を形成し得る活性生物学的ネットワークとして生体組織をモデル化する、準静的な連続体フレームワークを提案している。核となる仮定と数学的定式化は以下の通りである：

• 運動学と成長： 組織は局所的に非圧縮であるとモデル化される。形態形成は平衡状態の連続として記述される。成長は規定されるものではなく、局所的な機械的ストレスによって駆動される。変形勾配 $\mathbf{F}$ は、ネマチックステップ長テンソル $\mathbf{L}$、弾性歪み $\mathbf{F}_e$、成長テンソル $\mathbf{F}_g$、および緩和テンソル $\mathbf{F}_{rel}$ に乗法的に分解される。
• エネルギー定式化： 全自由エネルギー $\mathcal{F}$ は、主に2つの寄与からなる：
  1. ネマチックエネルギー ($\mathcal{F}_n$)： エリクセンの連続体理論に基づいており、ディレクター場 $\mathbf{n}$ およびスカラー秩序パラメータ $s$ の空間的変化を考慮する。決定的な点として、$s$ は空間的に変化させることができ、トポロジカル欠陥において消失することで特異性を正規化し、局所的な等方相への遷移を可能にする。
  2. 弾性エネルギー ($\mathcal{F}_e$)： 成長と緩和を含むように拡張されたネオ・クラシカル歪みエネルギー関数であり、フィラメントネットワークの歪みに対する応答を捉えるために、弾性的部分に依存する。
• 結合メカニズム： モデルは、(a) トポロジカル欠陥付近の機械的ストレス、(b) 活動誘起型の秩序化（$s$ の増加）、(c) ストレス駆動型の成長（弾性エネルギー密度に比例した新物質の追加）、および (d) 構造的緩和（内部ストレスを軽減するためのストレスフリーな構成の更新）という4つのメカニズムを統合している。
• 数値実装： 結合された非線形問題は、FEniCSx で実装された混合有限要素法を用いて解かれる。ネマチックディレクターの頭尾対称性と、それに伴う半整数欠陥の分岐カットを扱うため、ディレクター場は対称ランク2テンソル場 $\mathbf{M} = \mathbf{n} \otimes \mathbf{n}$ を用いて再定式化される。シミュレーションでは、初期の等方性回転楕円体から成熟した形態への進化を追跡するために、ラインサーチを伴うニュートン法と増分負荷を用いている。

主要な貢献

1. ABNのための連続体フレームワーク： 本論文は、ネマチック秩序、弾性、およびストレス駆動型の成長・緩和を結合させた数学的モデルを導入しており、線形シェル理論を超えて非線形な形態形成の結果を捉えている。
2. 欠陥駆動型の形態形成： ネマチック欠陥のトポロジー（特にその電荷と空間配置）が、形状進化の主要な駆動因子として機能することを実証している。本モデルは、欠陥が単なる特異点ではなく、ストレスを局在化させ成長を方向付ける組織化センターであることを示している。
3. 複雑な形態の再現： 初期欠陥構成を変化させることで、Hydra の再生で見られる3つの異なる形態形成の結果を再現することに成功した：
   • 一軸（頭部ー足部）再生： 対極にある2つの $+1$ 欠陥が伸長を駆動し、口（頂端部の破裂）と足部の形成をもたらす。
   • 二頭（二つの頭部）発生： 上半球にある2つの $+1$ 欠陥が $-1/2$ 欠陥によってバランスを取ることで、2つの異なる頭部を形成する。
   • 触手形成： 触手の先端にある $+1$ 欠陥と、その基部に位置する $-1/2$ 欠陥の構成が、触手の突出を駆動する。
4. ストレス駆動型成長メカニズム： シミュレーションは、成長が（欠陥付近の）高い弾性ストレスが存在する領域に局在すること、一方で緩和によって（幾何学的にストレスを軽減できる）開口部（口）が拡大すること、これは足部や触手先端における持続的なストレスとは対照的であることを明示的に示している。

結果
Hydra 再生に関する有限要素シミュレーションにより、以下のことが明らかになった：

• ストレスの局在化： ネマチックディレクター場の歪みが、トポロジカル欠陥付近で高い機械的ストレスを生成する。これらのストレスが、その近傍での新物質の蓄積（成長）を駆動する。
• 非対称な結果： 一軸構成では、ストレス緩和によって頂端部（口）が破裂して拡大するが、基部（足部）は閉じたまま構造的なストレスを保持する。
• 欠陥の役割： $+1$ 欠陥は突出形成を駆動し、組織の破裂を引き起こす。$-1/2$ 欠陥は、初期段階では主に構造的安定化装置または「トポロジカルなバランサー」として機能するが、触手構成においては、最終的に触手の基部における局所的な成長を誘発するストレス蓄積を引き起こす。
• 形態学的進化： モデルは、ネマチック秩序パラメータとそれによって生じるストレス駆動型の成長の進化のみを通じて、球状の等方性組織片から、特定の構造（口腔、触手など）を持つ複雑で伸長した構造への遷移を捉えている。

意義と主張
著者らは、本研究が、あらかじめ指定された成長パターンや外部のモルフォゲン勾配を必要とせずに、活性ネマチックネットワークにおけるトポロジカル欠陥がいかにして複雑な形態形成プロセスを駆動できるかについての、強固な理論的基礎を提供すると主張している。組織を、ネマチックの歪みによって生成される機械的ストレスに直接反応する活性生物学的ネットワークとして扱うことで、本モデルは、細胞骨格の微視的な組織化と、組織の形態という巨視的な形状との間の溝を埋めている。結果は、Hydra における実験的観察と一致しており、身体軸の極性の「記憶」や特定の構造の創発が、ネマチック欠陥のトポロジカル電荷と配置の中にエンコードされていることを示唆している。本論文は、純粋な生化学的シグナリングモデルとは異なる、生きている組織におけるパターン形成の固有の機械的駆動力を理解するためのステップとして、このフレームワークを位置づけている。