A multiscale cytoskeletal network model for shear rheological property and its evolutionary mechanism

技術要約：せん断レオロジー特性とその進化メカニズムに関するマルチスケール細胞骨格ネットワークモデル

問題提起
細胞骨格は、細胞の生理機能に不可欠な動的なバイオポリマーネットワークであるが、その機械的挙動は複数の時空間スケールにわたる。従来のモデルでは、クロスリンカーのキネティクス（動力学）に起因する動的な再構築と、フィラメントの特性に起因する粘弾性散逸メカニズムを統合的に捉えることが困難であるため、複雑な細胞応答を包括的に把握することには大きな課題が存在する。既存の連続体モデルは局所的な不均一性を捉えきれず、一方で離散構造モデルは、構成要素を純粋な弾性要素として簡略化しがちであり、時間依存的な特性や動的なクロスリンカーの進化を無視していることが多い。さらに、多くの離散モデルは計算上の制約から短時間のタイムスケール（≤ 60秒）に限定されており、定常状態のクリープのような長期的なレオロジー挙動を予測する能力に欠けている。本研究は、局所的なタンパク質の生化学的・力学的特性、アクチンフィラメントの粘弾性、およびそれらの変形状態を統合し、長期間（最大1,000秒）にわたるせん断レオロジー特性を予測するモデルの必要性に対処するものである。

手法
著者らは、電子顕微鏡の観察から得られた有限四角形細胞要素の代表単位に基づく、マルチスケール有限要素（FE）モデルを開発した。このモデルは以下の3つの異なるスケールを統合している。

1. フィラメントスケール： アクチンフィラメントを、前応力がかかった3パラメータ粘弾性固体ロッドとしてモデル化している。その構成挙動は、曲げ支配から張力支配への変形遷移を捉える時間依存的な緩和弾性率によって記述される。
2. ノードスケール： ノードは実験統計に基づき、以下の2種類に分類される。
   • 物理的絡み合いノード： 故障しないヒンジノードとして簡略化される。
   • 化学的クロスリンクノード： Bellモデルに従う、力によって制御された解離および結合速度を持つ回転バネとしてモデル化される。
3. ネットワークスケール： 一様な格子に対して初期角度とランダムな変位を導入することで、ランダムな統計的ネットワークを生成する。モデルは、ミリ秒から数千秒までの準静的平衡を効率的に解くために、入れ子状のニュートン・ラフソン反復スキームと、セグメント化された適応型タイムステップ戦略（対数スケール）を組み合わせた、境界修正型FE法を採用している。

数値解法においては、クロスリンカーの結合密度とフィラメントの粘弾性の変化を反映するように、剛性行列を動的に更新する。モデルは、ネットワークの再構築を特徴付けるために、バルクせん断歪み、見かけのせん断弾性率、および構造的シャノンエントロピーを算出する。

主な貢献

• 拡張された予測タイムスケール： 本研究は、離散ネットワークモデルの予測能力を約1,000秒まで成功裏に拡張し、従来のBrownian/FE-Langevinダイナミクスモデルの計算上の限界を克服した。
• 統合されたマルチスケール結合： 本モデルは、個々のフィラメントの粘弾性変形遷移（曲げから張力へ）と、クロスリンカーの生化学的・力学的キネティクスを、単一の数値フレームワーク内で独自に結合させている。
• 境界修正： 境界フィラメントにかかる偏心的な表面力を扱うために、特定の境界修正型FE法が開発され、せん断応力予測の精度が向上した。
• 病理および薬理学的シミュレーション： 本フレームワークは、癌化による軟化および薬剤介入のシミュレーションに適用され、微視的なパラメータの変化と巨視的なレオロジー結果との関連付けを行った。

結果

• レオロジー挙動： モデルは、短期的および中期的パワーロー（べき乗則）クリープ挙動、および長期的な歪みプラトーを含む、実験的に観察されたせん断レオロジー応答を再現した。パワーロー指数は、フィラメントの体積分率が増加するにつれて減少することが判明した。
• 進化メカニズム： シミュレーションにより、以下の4段階の経時的進化が明らかになった。
  1. 急速変化段階： フィラメントの曲げと高いネットワークエントロピーによって駆動される急速な歪みの成長。
  2. 遷移段階： クロスリンカーが解離し、変形が張力支配モードへと移行してエントロピーが減少することによる、成長の減速。
  3. 中期段階： ネットワークが安定するにつれ、フィラメントの粘弾性が支配的となる緩やかなクリープ。
  4. 長期プラトー： フィラメント固有の粘弾性特性によって決定される定常状態の挙動。
• 癌化現象： （フィラメントの前応力、粘性、密度が減少し、クロスリンカーの解離が増加するという特徴を持つ）癌細胞のシミュレーションは、加速されたネットワーク再構築と、大幅に低下した見かけのせん断弾性率を示した。
• 薬物の効果： 単一パラメータ分析により、以下が示された：
  • フィラメントの前応力を高めることは、定常状態のせん断弾性を効果的に回復させる。
  • フィラメントの粘性を高めることは、プラトーに達するまでの時間を延長させる。
  • クロスリンカーの解離速度を低下させることは、短期的弾性を高め、構造エントロピーのピークを遅らせる。これは、再構築プロセスが長期化することを示唆している。

意義
著者らは、本研究が細胞の機械的応答の時空間的進化と病理学的メカニズムの理解を前進させると主張している。クロスリンカーの生化学的・力学的結合とフィラメントの変形遷移の相乗効果が、ネットワークの再構築とパワーロー応答を支配していることを特定することで、本研究は観察されたレオロジー挙動に対するメカニズム的な説明を提供している。本モデルは、ポリマーネットワークの性能を制御するための理論的基礎を提供しており、具体的には、癌による軟化を効果的に逆転させるためには、フィラメントの前応力とクロスリンカーのキネティクスの両方を標的とした相乗的な介入が必要であることを示唆している。また、著者らは、マイオシンユニット、物理的絡み合いの移動性、および三次元ネットワーク効果の除外に関する現在の限界を認めつつも、このモデリング手法をハイドロゲルやコロイドゲルなどの他のポリマーシステムへ拡張できる可能性についても言及している。