本研究では、計算流体力学(CFD)解析のために、飛翔小管網(LCN)内の骨細胞の3次元(3D)モデルを可視化・開発する方法について概説します。この方法を使用して生成されたモデルは、健康な骨または病気の骨の骨細胞メカニズムを理解するのに役立ちます。
方法論記事
本研究では、計算流体力学(CFD)解析のために、飛翔小管網(LCN)内の骨細胞の3次元(3D)モデルを可視化・開発する方法について概説します。この方法を使用して生成されたモデルは、健康な骨または病気の骨の骨細胞メカニズムを理解するのに役立ちます。
骨細胞は、機械伝達として知られるプロセスでさまざまな生物学的経路を活性化することにより、機械的ひずみと流体流せん断応力(FFSS)に応答すると考えられている骨細胞です。骨細胞ネットワークの共焦点画像から得られるモデルは、直接測定では決定できない骨細胞膜のせん断応力を評価するための計算流体力学(CFD)解析を行うための貴重なツールです。骨の微細構造構造のこれらの高解像度画像を使用した計算モデリングを使用して、骨にかかる機械的負荷を数値的にシミュレートし、骨細胞の負荷誘発刺激を理解しました。
本研究では、Lacunar-Canalicular Network(LCN)の共焦点顕微鏡画像を用いて3次元単一骨細胞モデルを開発し、各種計算モデリングソフトウェアを用いたCFD解析を行う方法について詳しく説明します。共焦点顕微鏡検査の前に、マウスの骨を切片化し、フルオレセインイソチオシアネート(FITC)色素で染色してLCNを標識します。100倍の解像度で、Zスタック画像を共焦点顕微鏡で収集し、MIMICSソフトウェア(3D画像処理ソフトウェア)にインポートして、LCNおよび骨細胞樹状突起プロセスの表面モデルを構築します。
次に、これらの表面を3-Maticソフトウェア(3Dデータ最適化ソフトウェア)のブール演算を使用して減算し、骨細胞体周辺のラクナ流路空間とラクノカナル液を含む樹状突起の周りの小管空間をモデル化します。3D体積流体形状は、CFD解析のためにANSYSソフトウェア(シミュレーションソフトウェア)にインポートされます。ANSYS CFX(CFDソフトウェア)を使用して、流体圧力として骨に生理学的負荷を加え、骨細胞と樹状突起にかかる壁せん断応力を決定します。LCNの形態は、骨細胞膜と細胞プロセスによって感知されるせん断応力値に影響を与えます。したがって、共焦点画像ベースのモデルがどのように開発されるかの詳細は、骨細胞のメカニズムを理解する上で価値があり、この分野での将来の研究の基礎を築くことができます。
骨細胞は、身体運動に応答して骨量を調節すると仮定されています1。骨細胞の膜変形とその樹状突起は、機械的負荷によって骨細胞にFFSSをさらし、これは骨細胞によって検出され、細胞内シグナル伝達を誘発します2,3,4。骨の微細構造は、加齢や骨粗鬆症や糖尿病などの骨疾患、および骨細胞の機械応答性障害を引き起こすペルルカン欠乏症などの状態により、そのラクナ-小管形態の劣化または変化を経て受けます5,6。骨構造のこれらの変化により、骨細胞はさまざまなレベルのFFSSとひずみを経験する7,8。重要なことに、機械的負荷に応答して骨細胞が経験するFFSSは、石灰化した骨基質に埋め込まれているため、in vivoで定量化するのが困難です。
共焦点画像ベースのモデリングは、LCN 9,10のコンピューターモデルを複製することにより、自然環境でアクセスできない骨細胞を研究するという制限を克服する強力な手法です。LCNの相互接続されたネットワークを3Dで処理し、モデル化することは困難でした。透過型電子顕微鏡法(TEM)、走査型電子顕微鏡法(SEM)、シリアルブロック面切片法、およびシリアル集束イオンビーム走査型電子顕微鏡法(FIB / SEM)2,11,12など、いくつかのイメージング技術があります。骨13,14,15を視覚化し、共焦点レーザー走査型顕微鏡(CLSM)を介して3D骨細胞モデルを生成するための貴重な技術が開発されました。CLSMは、すべてのラクナボリュームとほとんどの小管を3Dでイメージングできるため、他のイメージング技術ではなく、計算モデリングのためにここで選択されました16,17。LCN形状は、骨細胞有限要素解析(FEA)のCLSMを使用して生成し、骨のひずみを予測できます。ただし、骨細胞が経験するFFSSを予測するための流体解析は、骨細胞の細胞膜とLCN内の樹状突起のモデリングが必要であり、間質液が18の周りを移動する狭いラクナ-小管腔のモデリングを可能にする必要があるため、より複雑です。
このプロトコルでは、フルオレセインイソチオシアネート(FITC)色素を脱灰していない厚い骨切片に適用してから、共焦点顕微鏡で骨内のLCNを標識し、骨細胞樹状突起膜をLCNからのイメージングデータに基づいてモデル化します。ラクナ-小管腔は計算モデリングを使用してシミュレートされ、身体活動による生理学的負荷はCFDアプローチを使用してモデル化されます。骨細胞は、CFDソフトウェアで流体圧力勾配を受け、LCN内の流体プロファイルを分析し、骨細胞および樹状膜上のFFSSを測定します。さらに、FEAアプローチでは、圧縮力学的負荷を加えることにより、骨細胞のひずみや応力を測定することができます。
また、若くて健康な骨の画像から得られた微細構造を修正し、老齢動物や骨疾患のある動物における裂孔小管の形態の変化をシミュレートするための形状修正技術も開発されました。骨の微細構造の変更には、老化に伴う小管の数の減少、ペルルカン欠損症で何が起こるかをモデル化するためのラクナ-小管腔面積の減少、および老化効果をモデル化するためにそれを増やすこと、および糖尿病骨をモデル化するための小管および樹状突起壁面積の減少が含まれていました5,6。形状修正技術により、若年動物と老人動物、健康な動物と病気の動物の骨など、さまざまな微細構造を持つ骨の骨細胞が経験するFFSSを比較できます。
全体として、共焦点画像ベースのモデリングは、健康な骨の骨細胞の形態だけでなく、骨細胞の形態の老化/疾患関連の変化をシミュレートするための貴重なツールです。さらに、ラクナ-小管腔の表面積や体積などの骨細胞の形態学的パラメータをさまざまな骨で測定および比較して、機械的ひずみに対する細胞応答を予測できます。
動物実験は、ミズーリ大学カンザスシティ校(UMKC)の動物管理・使用委員会の承認を得て実施され、関連する連邦ガイドラインに準拠していました。
1. 骨の準備プロセス
2. 共焦点顕微鏡
3. コンピュータモデリング
4. 3次元画像処理ソフトウェアと3次元データ最適化ソフトウェアによる形状修正技術
注:ジオメトリ修正技術は、加齢や骨疾患による小管密度や直径、ラクナ小管の厚さなど、骨細胞の形態の変化をモデル化するために使用されます。
5. CFD解析
注: 体積骨細胞モデルを生成した後、シミュレーション ソフトウェアの CFX モジュールで、ジオメトリ、メッシュ、セットアップなどのいくつかの手順が実行されます。
6. CFD後処理
このプロトコルでは、機械的負荷により骨細胞とその樹状突起が受ける流体の流れせん断応力の量を調査するために、共焦点由来の骨細胞モデルを開発する方法について説明します。老化したC57BL6マウスと若いC57BL6マウスを選択し、若くて老化した共焦点画像ベースの骨細胞モデルを作製しました。他の6つのシミュレートされた骨細胞モデルは、同じ若い骨細胞モデルからジオメトリ修正技術を使用して生成され、加齢または骨疾患によるLCN形態の変化を研究しました。ジオメトリ修正されたパラメータには、ラクナ-小管の厚さ、小管密度、および直径が含まれていました。改変骨細胞モデルは、一度に1つのパラメータを変更し、他のパラメータを同じにすることで開発されました。したがって、FFSSに対する各形態学的パラメータの影響は、単独で調べられます。その結果、これらの骨細胞モデルでは、骨細胞が経験するFFSSの量が異なることが明らかになり、FFSSがLCNの形態に依存していることが示されました。FFSS値は、骨細胞体や樹状膜など、骨細胞モデルのさまざまな部分で予測されました。すべてのモデルの樹状突起は、骨細胞体膜よりも大きなFFSS値を示しました。
図5は、共焦点FITC画像スタックから直接生成された若年骨細胞モデルと老化骨細胞モデルのFFSS輪郭と、ジオメトリ修正技術を使用してモデル1から開発された他の6つの骨細胞モデルのFFSS輪郭を示しています。若年骨細胞モデル(モデル1)の骨細胞および樹状膜上の平均FFSSは0.42Paであり、これは老化骨細胞モデル(モデル2)の平均FFSS0.13Paよりも有意であった。得られたデータは、ラクナ小管腔が大きい骨細胞(モデル7およびモデル8)で平均FFSSが増加し、小管密度が低い骨細胞(モデル3およびモデル4)でFFSSが減少するのに対し、モデル5およびモデル6では樹状突起径を変更してもFFSSは有意に変化しなかったことを示しています。モデル2とモデル4のFFSS値は、それぞれ0.19 Paと0.13 Paと最も低い値でした。モデル8のFFSSの大きさは、他の骨細胞モデルよりも大きいです。モデル8では、樹状突起のほとんどが赤色であるだけでなく、骨細胞体も他のモデルよりも高いFFSSを経験します。
また、ANSYS CFXでの後処理により得られたラクナ-小管腔の体積と表面積に基づいて、骨細胞が経験するFFSSが提示され、FFSSとLCNの構造変化との相関関係が明らかになりました。また、骨形成を引き起こすFFSSのレベルと考えられる0.8Paを超えるFFSSに曝露された各モデルの表面積を測定することもできた21,22。骨細胞モデリングは、表1に示すように、すべての骨細胞モデルについて、平均FFSSとラクナ-小管腔の体積または表面積との間に相関関係があることを示しました。老齢骨細胞(モデル2)と模擬中年骨細胞(モデル4)で小管密度が最も低いものは、ラクナ小管表面積が1372.7μm2と1122.6μm2と最も低かった。また、モデル2とモデル4は、ラクナ-小管腔の体積がそれぞれ229.4 μm3と203.5 μm3と最も低く、FFSSの量が少ないことと相関していると考えられる。樹状突起の直径が最大のモデル6は、ラクナ-小管表面積が最大でした。モデル8は、最大のラクナ-小管腔体積と最大のFFSS値を持っていました。

図1:5倍、20倍、100倍のオイル対物レンズで撮影されたFITC染色された大腿骨の横断面の単一Z平面共焦点画像。 Lacunar-canalicularネットワークは、100倍対物レンズで撮影した画像に写っています。 この図の拡大版を表示するには、ここをクリックしてください。

図2:生後4ヶ月と22ヶ月齢のC57BL6マウスのFITC染色大腿骨から作製した若年骨細胞(モデル1)と老齢骨細胞(モデル2)の共焦点画像ベースモデル。 拡大画像は、若い骨細胞のラクナ-小管(LC)の厚さと、老化した骨細胞の同じ厚さを示しています。モデル1の画像は、Niroobakhsh et al.23の許可を得て変更されています。 この図の拡大版を表示するには、ここをクリックしてください。

図3:モデル1からジオメトリ修正技術を使用して開発された6つの骨細胞モデル。 小管密度を減少させてモデル3(シミュレートされた中年)とモデル4(シミュレートされた老化)を生成しましたが、樹状突起の直径とLCの厚さはモデル1と同じでした。モデル5とモデル6は、モデル1と同じ小管密度とLC厚さで開発されましたが、樹状突起の直径はモデル1の2倍と4倍大きかった。モデル7とモデル8は、それぞれ2倍と4倍のラクナ-小管の厚さを持っていましたが、小管密度と樹状突起の直径はモデル1と同じままでした。モデル1の画像は、Niroobakhsh et al.23の許可を得て変更されています。 この図の拡大版を表示するには、ここをクリックしてください。

図4:3Dビュー共焦点画像ベースの若い骨細胞モデル。 骨細胞と樹状突起の周囲の流体解析に使用されるラクナ-小管腔のモデリングを描いた共焦点画像ベースの若い骨細胞モデルの断面積の3Dビュー。 この図の拡大版を表示するには、ここをクリックしてください。

図5:すべての骨細胞モデルのFFSS分布は、FFSSの大きさが骨細胞の形態に依存していることを示しています。 赤色の領域は、特に骨細胞体ではなく樹状膜で最も高いFFSS値を示しています。この画像は、Niroobakhsh et al.23の許可を得て変更されています。 この図の拡大版を表示するには、ここをクリックしてください。
| モデル | 腹膜小管表面積 (μm2) | 外堀小管腔容積 (μm3) | アベニューFFSS(PA) | |
| 共焦点画像ベースの骨細胞モデル | モデル1 | 2652.9 | 261 | 0.42 |
| ヤング(生後4ヶ月) | ||||
| モデル2 | 1372.7 | 229.4 | 0.19 | |
| 年齢(22ヶ月) | ||||
| ジオメトリ修正技術を使用して生成されたシミュレートされた骨細胞モデル | モデル3 | 1530.4 | 242.6 | 0.27 |
| (シミュレートされた中年) | ||||
| モデル4 | 1122.6 | 203.5 | 0.13 | |
| (シミュレーションエイジング) | ||||
| モデル5 | 3724.9 | 302.9 | 0.49 | |
| (2 x 樹状突起直径) | ||||
| モデル6 | 4791 | 337 | 0.57 | |
| (4 x樹状突起直径) | ||||
| モデル7 | 3116.3 | 440.6 | 1.17 | |
| (2 x LC厚さ) | ||||
| モデル8 | 3542.2 | 658.2 | 1.89 | |
| (4 x LC厚さ) |
表1:骨細胞モデルの比較。 2つの共焦点画像ベースのモデルと、ジオメトリ修正技術を使用してモデル1から生成された6つのシミュレーションモデルを含む、8つの異なる骨細胞モデルにおける平均FFSS、ラクナ-小管空間体積、および表面積の比較。
このプロトコルは、骨細胞の視覚化と計算モデリングのための共焦点イメージング技術の概要を示しています。共焦点イメージングの前に、骨サンプルを切片化および染色するための骨調製プロセスが行われます。100倍の倍率の共焦点画像は、骨細胞とラクナ-小管腔のコンピューターモデルを開発するために、さまざまなソフトウェアにインポートされます。最後に、骨細胞や樹状膜を囲む身体活動によるFFSSをモデル化するための共焦点画像ベースモデルを用いてCFD解析を行いますが、これは骨基質中の骨細胞にアクセスできないためin vivo では技術的に困難です。提供されたプロトコルは、マルチフィジックスモデリングで有限要素法または流体-構造連成(FSI)アプローチを使用して他の応力やひずみを決定するのにも役立ちます。
共焦点イメージングは、骨のFITC染色を使用して達成できる高解像度画像により、LCN 7,9,24,25の視覚化に使用されてきました。小管の直径は顕微鏡の解像度限界に近いですが、信頼性の高い同焦点の小眼小管画像を取得でき、小管の大部分を捉えることができます。骨細胞の樹状突起の直径が小さいため、樹状突起の画像を正確に撮影することは技術的により困難になります。そのため、ソフトウェアのスムージング操作を使用して骨細胞と樹状突起の体積を減らすことにより、LCNの共焦点画像から骨細胞と樹状突起をモデル化するための有用な技術が開発されました。他のいくつかの発表された研究では、同様のアプローチを使用してLCNをオフセットし、骨細胞樹状膜を生成しました9。スムージング技術のもう1つの利点は、光学イメージングの特徴であるZ軸のイメージング解像度の低下によりZ軸で発生する小管径の伸びなど、CLSMの一部のイメージングエラーが改善されることです。デコンボリューションは、元のデータセットの信号対雑音比を改善するため、Z 方向の光学ストレッチを減少させる代替方法とも考えられています。ラクナ-小管腔は、LCNと骨細胞-樹状膜からモデル化され、LCN内の流体プロファイルを分析します。
特に、既存の方法に関して、別のマウスを犠牲にすることなく幾何学的に改変された骨細胞モデルを開発するための新しい技術が導入されました。ジオメトリ修正技術は、骨粗鬆症、糖尿病、ペルルカン欠乏症などの状態や骨の老化や骨疾患による骨の微細構造の変化をシミュレートするのにも役立ちます。この手法では、若いマウスの共焦点画像ベースの骨細胞モデルを使用して、小管密度、ラクナ-小管腔、樹状突起径などのパラメータを変更することにより、老化した骨細胞または他の骨疾患の骨細胞の形態をモデル化しました。このアプローチは、22か月齢のマウスの実際の骨細胞の共焦点画像ベースのモデリングと、4か月齢のマウスから最初に画像化された骨細胞の小管密度を低下させることによる形状修正モデリングの2つの老化骨細胞モデルの類似性によって検証されました。これら2つの老化骨細胞モデルは、すべての骨細胞モデルの中で最も低い同様のFFSSの大きさを示しました。ジオメトリ修正技術は、動物の実際の骨細胞を使用して生成された骨細胞モデルに匹敵する結果を予測できるため、さまざまな生物学分野で多くの応用が期待できます。この作業の限界には、平滑化技術による表面粗さの詳細の一部が失われることや、小径の小管の分解能の限界にあるという事実、小径の小管の精度が低下する可能性があること、振動流を使用したモデリングではなく流体の流れの理想化などがあります。
この技術の機能を実証するために、若いマウスと老齢のマウスから生成された、または形状修正技術を使用して、骨細胞のFFSS値を異なる骨細胞モデルで測定しました。これらのコンピュータモデルを使用した注目すべき発見の1つは、ラクナ-小管腔の体積が骨細胞が経験する平均FFSSに正比例することであった。この新しい発見は、これまで調査されていなかったものですが、骨細胞ネットワークの形態を骨細胞とそのプロセスによって検出されたFFSSに関連付けることができるため、骨細胞の機械応答性を決定する上でさまざまな用途がある可能性があります。
著者は何も開示していません。
著者らは、全米科学財団(NSF、賞番号NSF-CMMI-1662284 PI:T Ganesh)、国立衛生研究所(NIH - NIA P01 AG039355 PI:LF Bonewald)、(NIH/SIG S10OD021665およびS10RR027668 PI:SL Dallas)、およびミズーリ大学カンザスシティ大学院研究大学院研究助成プログラムに感謝します。
| 名前 | 会社 | カタログ番号 | コメント |
|---|---|---|---|
| 1,200グリットサンドペーパー | Buehler | 30-5170-012-100 | |
| 3-Maticソフトウェア | Materialise | https://www.materialise.com/en/industrial/software/3-matic | 3Dデータ最適化ソフトウェア |
| 600グリットサンドペーパー | Buehler | 30-5118-600-100 | |
| 800グリットサンドペーパー | Buehler | 30-5170-800-100 | |
| ANSYSソフトウェア | ANSYS | https://www.ansys.com/ | シミュレーション ソフトウェア |
| フルオレセイン イソチオシアネート (FITC) | Sigma-Aldrich | F7250 | |
| ImageJ ソフトウェア | https://imagej.net/ij/ | ||
| 微鏡用浸漬油 | ライカ マイクロシステムズ | 195371-10-9 | |
| ライカ TCS Sp5 II 共焦点顕微鏡 | ライカ マイクロシステムズ | TCS Sp5 II | |
| ライツ 1600 インナーホール ダイヤモンド ソー | ライカ | ||
| MIMICS Innovation Suiteソフトウェア | Materialise | https://www.materialise.com/en/healthcare/mimics-innovation-suite | 3D画像ベースの処理ソフトウェア |
| パーマウントマウント培地 | Fisher scientific | SP15-500 | |
| Sampl-Kwick Fast Cure Acrylic Kit | Buehler | 20-3560 | |
| Single Platform Laboratory Shaker | Reliable scientific INC | モデル55S |
このJoVE記事のテキストまたは図の再利用許可をリクエスト
許可をリクエスト