[論文レビュー] Linking Cellular and Mechanical Processes in Articular Cartilage Lesion Formation: A Mathematical Model

G Kapitanov, Xiayi Wang+3 more|arXiv (Cornell University)|Jan 6, 2016

Osteoarthritis Treatment and Mechanisms参考文献 48被引用数 19

ひとこと要約

本研究では、有限要素法を用いた機械的シミュレーションと年齢構造付きの細胞および反応拡散化学ダイナミクスを統合した新規のマルチスケール数学的モデルを提示し、関節軟骨における外傷後変形性関節症病変形成を模擬する。モデルは、初期の機械的負荷が壊死性細胞死を誘発し、ROS、DAMPs、PIC、EPOといった炎症性および抗炎症性シグナルが組織内を拡散する連鎖反応を引き起こすことを示している。14日間の進行において、PICおよびDAMPの分解速度（λF、λM）といった重要なパラメータが病変進行および細胞状態遷移に顕著に影響することが判明した。

ABSTRACT

A severe application of stress on articular cartilage can initiate a cascade of biochemical reactions that can lead to the development of osteoarthritis. We constructed a multiscale mathematical model of the process with three components: cellular, chemical, and mechanical. The cellular component describes the different chondrocyte states according to the chemicals these cells release. The chemical component models the change in concentrations of those chemicals. The mechanical component contains a simulation of pressure application onto a cartilage explant and the resulting strains that initiate the biochemical processes. The model creates a framework for incorporating explicit mechanics, simulated by finite element analysis, into a theoretical biology framework.

研究の動機と目的

• 関節軟骨における機械的損傷と生化学的・細胞応答を結びつけるメカニズム的枠組みを構築すること。
• 既存のモデルにおける欠落を補うために、明示的な有限要素法を用いた機械的応力シミュレーションと、細胞および化学ダイナミクスのバイオマス的モデルを統合すること。
• 損傷後、軟骨細胞死と炎症の拡散を駆動する主要な閾値およびパラメータ感受性を特定すること。
• さまざまな機械的および生化学的条件下での治療効果がPTOA発症に与える影響を評価する予測ツールを提供すること。

提案手法

• 機械的部品では、Abaqusを用いた有限要素解析を用い、金属インデンタがシリンダ型軟骨エクスプラントに圧力を加えるのをシミュレートし、応力分布および初期の壊死性細胞密度を計算する。
• 細胞部品では、年齢構造付きモデルを用い、各状態（CU、CT、CE、ST、SA、DA、DN）に滞在した時間を基に軟骨細胞状態を追跡する。遷移は化学的シグナルによって駆動される。
• 化学部品では、ROS、DAMPs、プロインフレマトリー接着性サイトコイニン（PIC）、およびEPOの濃度を、状態依存的生産および分解を伴う反応拡散偏微分方程式でモデル化する。
• モデルは、機械的応力と初期の壊死性細胞密度を結びつけ、DAMP放出を引き起こし、炎症および抗炎症応答を含むフィードバックループを通じて生化学的連鎖を開始する。
• 主要パラメータ（例：λF、λM、β11）の感度分析を実施し、細胞および化学ダイナミクスに与える影響を評価する。
• 数値解法として、陰解法を用いた差分法を用い、2次元径方向軸方向空間で損傷後14日間の進行をシミュレートする。

実験結果

リサーチクエスチョン

• RQ1 初期の機械的負荷が関節軟骨における壊死性細胞死をどのように誘発し、その後の生化学的シグナル伝達をどのように開始するか？
• RQ2 ROS、DAMPs、PIC、EPOといった主要なシグナル分子が損傷後、軟骨細胞死および炎症の拡散に果たす役割は何か？
• RQ3 分解速度（λF、λM）および生成速度（β11、β13）のパラメータ変動が病変形成の拡散に与える影響はどのようなものか？
• RQ4 有限要素法を用いた応力データの統合が、マルチスケール軟骨損傷モデルの生物学的妥当性および予測能力をどの程度向上させるか？
• RQ5 炎症連鎖における臨界閾値または転換点を同定できるか？それらは治療的標的として有効であるか？

主な発見

• モデルは定性的に予想される挙動をうまく再現した：炎症メディエーター（PIC、ROS、DAMPs）および分解代謝細胞状態（ST、SA）のピーク濃度は、圧力負荷領域の中心部で最も高くなることが確認され、損傷部位の局在と整合的であった。
• PIC分解速度λFを0.1に低下させた場合、ST、SA、およびROSレベルに顕著な上昇が観察された。これは、炎症性サイトコイニンのクリアランス能力の低下が病変進行を加速することを示唆している。
• DAMP分解速度λMを0.1に低下させた場合、14日目までに健全な軟骨細胞（CU）がほぼ完全に消失し、炎症性化学物質の顕著な拡散が観察された。これは、DAMPの持続的保持が持続的炎症の主要な駆動要因であることを示している。
• 感度分析により、生物学的に妥当な反応が確認された：β11（PIC生成）およびβ13（EPO抑制）の摂動は、細胞状態遷移およびサイトコイニン濃度のダイナミクスに顕著な影響を及ぼした。
• EPO濃度が細胞状態の逆転（EPOR活性状態のSAから健全状態のCUへの移行）に濃度依存的効果を示した。これは、EPOが内因性の抗炎症調節因子としての役割を果たす可能性を支持する。
• 有限要素法を用いた応力データのバイオマス的フレームワークへの統合により、生化学的連鎖の現実的な開始が可能となり、機械的入力からの損傷誘因病態をシミュレートするモデルの有効性が検証された。