Waiting
Login processing...

Trial ends in Request Full Access Tell Your Colleague About Jove
Click here for the English version

Bioengineering

心筋細胞のシステム生物学における細胞構造の役割を研究する心筋細胞の構造的に現実的な有限要素の幾何学的なモデルを作成します。

Published: April 18, 2018 doi: 10.3791/56817

Summary

このプロトコルでは、電子顕微鏡や共焦点顕微鏡画像から心筋細胞の細胞内構造の空間的有限要素モデルを作成する手法について説明します。カルシウム シグナルと生体エネルギー転換の事例を使用してこの空間詳細モデルの力を発揮します。

Abstract

電子線トモグラフィーなど三次元 (3 D) 画像処理技術の出現により、シリアル ブロック面走査電子顕微鏡と共焦点顕微鏡は、科学的なコミュニティは、大規模なデータセットへの前例のないアクセス サブミクロン建築改造を特徴づけるの解像度は、健康と病気心筋細胞機能の変化を伴います。ただし、これらのデータセットは、心筋細胞の機能細胞アーキテクチャの役割を調査するため利用されています。このプロトコルの目的は、高分解能電子顕微鏡と共焦点顕微鏡画像を用いた心筋細胞の高精度有限要素モデルを作成する方法の概要を説明することです。細胞アーキテクチャの詳細かつ正確なモデルには、心筋細胞生物学実験だけを集めることができます以上に新しい洞察力を提供するために重要な可能性があります。このメソッドの力は、計算上、心筋細胞の 1 つの統一された、詳細なモデルを開発する心筋細胞の微細構造の 2 つの異なる画像モダリティからの情報を融合する能力であります。このプロトコルは、電子線トモグラフィーと、心筋細胞の半筋節の有限要素モデルを開発する男性の大人の Wistar (白鼠の特定の品種の名前) ラット心筋細胞の共焦点顕微鏡画像を統合する手順をについて説明します。手順は正確、高解像度の描写を含む三次元有限要素モデルを生成する (〜 35 程度 nm) ミトコンドリア、筋原繊維と心筋細胞の必要なカルシウムを解放するリアノジン受容体クラスターの分布の筋原線維と細胞内コンパートメントに漿網状のネットワーク (SR) から収縮します。生成されたモデルとおりイラストが横行尿細管アーキテクチャまたは漿の網状のネットワークの詳細を取り入れていない、したがって、心筋細胞の最小モデル。それにもかかわらず、モデル既に適用できるカルシウム シグナル伝達とミトコンドリア生体エネルギー、示されていますと、次が表示されます 2 つの事例を用いて説明したセル構造の役割に基づくシミュレーション調査で、詳しいプロトコル。

Introduction

興奮収縮連関 (ECC) 中心部には、重要な複雑な心筋細胞膜の電気的興奮と各ハートビートの中にセルの後続の機械的収縮カップリングを指します。数理モデルは活動電位1,2、生体エネルギー3、ゾル性細胞質カルシウムを調節する連結された生化学的なプロセスの定量的理解の開発に重要な役割を果たしているとその後収縮力発生。このようなモデルは、1 つのハートビートの変更を予測しているも正常にまたはこれらの生化学的なプロセスのいくつかは変更4,5を受けます。心筋細胞の非常に組織の微細構造は、セルと心臓全体の正常な収縮機能に重要な役割を果たすにますます認識されています。確かに、肥大6心不全7、および糖尿病性心筋症8疾患の生化学の変化と並行して形態と心筋の微細構造のコンポーネントの構成の変更が発生します。これらの構造変化が変化しマイナー、適応、または病理学的応答かどうか生化学的な条件はまだ未知9です。生物の形と機能の間の本質的に密結合を意味する実験的研究だけでは構造を改造と心筋細胞の機能間の相関関係よりも深い洞察を提供できません。よく研究の生化学的なプロセスと、細胞内のコンポーネントの構造のアセンブリを組み込むことができる数理モデルの新世代の関係の定量的、総合的な理解を開発する必要が構造, 生物化学と心筋収縮力。このプロトコルでは、このような調査に使用することができます心筋細胞の構造的に高精度有限要素モデルの生成に使用できるメソッドについて説明します。

最後の十年は、3次元電子顕微鏡10、共焦点11、およびナノ ・ マイクロ スケールの組み立てに前例のない、高解像度の洞察力を提供する超解像顕微鏡12の重要な進歩を見ている、心筋細胞の細胞内コンポーネント。最近では、これらのデータセットは、心筋の微細構造13,14,15,16の計算モデルを生成する使用されています。これらのモデルでは、有限要素法17と呼ばれる確立されたエンジニア リング シミュレーション法を使用して、作成有限要素計算メッシュする生化学的なプロセスおよび心筋細胞の収縮をシミュレートできます。ただし、これらのモデルは解像度とディテール画像データセットに提供できる顕微鏡法によって制限されます。たとえば、電子顕微鏡は、セル構造のナノメートル レベルの詳細を生成できますが、モデルを作成する必要がありますイメージ内の特定の蛋白質を識別することは困難です。超解像光学顕微鏡で 50 程度の解像度で高コントラストな画像を提供できる一方でのみ、選択数分子コンポーネントのセルの nm。これら画像モダリティから補完的な情報を統合することによってのみすることができます 1 つ現実的に感度を調べる関数の構造の変化に。相関光学および電子顕微鏡はまだルーチンの手順と制限コンポーネント数が限られている可能性があります蛍光表示で染色、電子顕微鏡ビューと相関も苦しむことになります。

このプロトコルは、統計的手法19を使用して分析し、電子顕微鏡について他の心筋イオン チャネルの空間分布については光顕を融合計算手法18を提示します。筋原線維などミトコンドリアの微細構造コンポーネント。これは、心筋細胞の収縮を調節する生化学的なプロセスの心筋細胞内組織の役割を研究する生化学的なプロセスの生物物理のモデルで使用することができます有限要素モデルを生成します。たとえば、このプロトコルを使用してから健康的なモデルを作成できる、ストレプトゾトシン誘発糖尿病心筋細胞は糖尿病動物で観察される心筋細胞機能構造リフォームの効果を研究するモデル8。提案手法の統計的性質の付加的な利点はプロトコルにも示されている: メソッドは、セル構造の実験的に観察された変化を忠実に再現した有限要素のジオメトリの複数のインスタンスを生成することができます。

概要については、としてプロトコルの手順が含まれます: (i) 十分な解像度とコントラストの 3 D イメージを生成する電子顕微鏡検査のための心筋組織の準備(ii) 復興 3次元電子顕微鏡再構成および画像解析ソフトウェアを使用して電子顕微鏡データから 3 D 画像のスタックのセグメンテーションと呼ばれる IMOD20;(iii) iso2mesh21を使用して入力としてセグメントのデータを使用して有限要素メッシュを生成するには(iv) 新規アルゴリズムとコードを使用して、有限要素メッシュにイオン チャネルの分布をマップします。

各ステップへのアプローチの前提は、プロトコルに記載されている、代表的な結果は付属の数字で提供されます。生成された空間の詳細なモデルを使用するミトコンドリアバイオエナジェティックスと同様に、ECC、間にカルシウムの空間動態を勉強する方法を指定する概要を説明します。それらを克服し、さらに心臓システム生物学にセル構造の役割を定量的に理解を進める進められている新たな展開だけでなく、プロトコルの現在の制限のいくつか、説明します。他の細胞型の有限要素モデルを作成するこれらのメソッドを一般化する可能性がどのようにも解決されます。

このプロトコルのユーザーは、既存の電子顕微鏡画像のスタックへのアクセスを有すればステップ 1 とステップ 2 の復興の一部飛ばしてください。経験豊富な電子と共同で、データを取得しようとすると、ユーザーは、議論し、固定および汚損プロシージャ獲得のための最適なプロトコルを決定する専門家で手順 1 で比較するがあります。

Subscription Required. Please recommend JoVE to your librarian.

Protocol

ここで説明したすべてのメソッドは、オークランド動物倫理委員会の大学、カリフォルニア大学、サンディエゴ動物介護施設ティッシュのプロトコルが開発された当初、利用委員会によって承認されています。

1. 実験準備

  1. 0.15 M と 0.3 M ナトリウムの溶液 ph 7.4の表 1によると cacodylate のバッファーを準備します。
  2. グルタルアルデヒド固定液 (pH 7.4 で 0.15 M ナトリウム cacodylate バッファーの 0.2% タンニン酸 2.5% グルタルアルデヒド + 2% パラホルムアルデヒド (PFA)) を準備コンポーネントを使用して表 1 に示します。
    1. 一定攪拌しながらホット プレート上 60 ° C で蒸留水 20 mL に PFA 粉末 2 g を溶解します。
    2. ソリューションがクリアされるまでに 1 M 水酸化ナトリウム溶液を追加します。
    3. 溶液の温度が 40 ° C を下回るまで待つ、グルタルアルデヒドと 0.3 M ナトリウム cacodylate 20 mL 10 mL を追加します。
    4. タンニン酸 0.2 g を追加し、ソリューションに溶解します。
    5. 0.15 M ナトリウム cacodylate の 50 mL を追加します。
    6. 同じ日に使用する場合、4 ° C で氷または冷蔵庫の中に定着剤の 3 つまたは 4 つの 20 mL シンチレーション バイアルを格納します。
  3. 表 1のレシピによると 20 mM 2, 3 butanedione ジアセチルモノオキシム (BDM) とタイロード液を準備します。
  4. シリカフューム戸棚内部蛍光装置を構築します。
    1. 2 つの異なるレトルト スタンド、スタンド ベースから高さ約 70 cm に 2 つ 100 mL プラスチック注射器チューブをクランプしてください。
    2. プラスチック注射器管と 3 mm 内径のチューブと三方活栓を使用して、図 1に示すように、コネクタ チューブを接続します。
    3. チューブ、心臓を灌流固定接続されますコネクタの下に 3 mm 外径カニューレをフィットします。
  5. タイロード液と 37 ° C で固定液注射器管を埋める
    1. それらを通して注射器ソリューションを渡すことによりチューブ内の空気の泡を削除します。

2. 心筋細胞微細構造の電子顕微鏡データを取得します。

  1. 切除と心の化学固定の動物を準備します。
    注: このプロトコルは、大人男性 Wistar (白鼠の特定の品種の名前) ラット心筋組織を処理する手順を説明します。プロトコルは、心筋組織が他の種から供給されるとき変更を必要があります。
    1. 腹腔内投与で 250 U/mL のヘパリンの 0.5 ml (体重 200-250 g) ラットを扱うことによって動物を麻酔します。ペントバルビ タール (210 mg/kg 体重) の腹腔内注入ラットを扱うまで、10 分間待ちます。
    2. つま先ピンチ反射をテストすることによって麻酔の適切な学位を確認します。
    3. 頚部転位によって、動物を安楽死させます。
    4. 以前公開されたゼウス記事に似たような郭清術を使用して心22,23、それからそれを置き氷冷生理食塩水での収穫。
    5. 上行大動脈を隔離し、cannulate。カニューレの先端は、冠状動脈が分岐を半月弁のすぐ上座っていることを確認します。上行大動脈の周りにしっかりと糸を結ぶ。
    6. カニューレを 〜 70 cm システム (図 1) の基盤上で重力駆動蛍光装置に接続します。
    7. 2 〜 3 分の 20 mM BDM (ミオシン阻害剤) を含むタイロード液で心臓を灌流して蛍光システムの活栓をねじる。
    8. 0.15 M ナトリウム cacodylate 37 ° c 10 分の 0.2% タンニン酸、2.5% グルタルアルデヒド 2% パラホルムアルデヒド固定液による血液灌流を開始するコックをひねます。成功した場合、中心部が淡い茶色を回すし、硬くなりゴムのような。
    9. 左心室 (側の) 自由壁からティッシュのブロックを解剖し、かみそりの刃を使用すると、組織ブロック約 1 mm3の大きさ。
    10. ストア 2 h. のための氷に同じ固定液の予冷 20 mL シンチレーション バイアルのサンプル、バイアルにできるだけ多くのサンプルを保存してソリューションのサンプルが水中に沈むこと。
      注意: サンプルの下の 100% の脱水のステップの後まで寒さを維持する重要です。
  2. さらに固定および電子顕微鏡用試料の汚損します。
    1. ガラス ビーカーに 0.15 M ナトリウム cacodylate バッファーを 100 mL 注ぎ、氷で冷やします。
    2. 4% 四酸化オスミウムと 0.3 M ナトリウム cacodylate、フェロシアン化カリウム 2% 2% 四酸化オスミウムと 0.15 M ナトリウム cacodylate を含む重金属染色溶液にフェロシアン化カリウム 0.08 g の添加に続いての平等なボリュームを準備します。氷の上の解決策を冷やして下さい。
    3. ピペットのうち定着剤、バイアルに試料が水没する十分な冷たい 0.15 M のナトリウム cacodylate バッファーを置き換えます。5 分間氷の上には、バイアルを配置します。
    4. 手順 2.2.3 0.15 M ナトリウム cacodylate 余分な定着剤を削除するとあと 4 回。
      注意: 小さな破片がサンプルに得ることができるティッシュ ブロック中にダイヤモンド ナイフを台無しにすることができますのでガラス ピペットを使用しないでください。事前のサンプルにそれらを追加する前に氷の上のすべてのソリューションを冷却することが重要です。サンプルに (部分的なカバレッジ以上数秒の非常に短い期間容認されるかもしれない簡単にソリューションの交換中に) すべての時間ソリューションによって覆われて残ることが重要です。洗浄またはソリューションを置き換えると、以前のソリューションを削除した後新しいソリューションを迅速に追加します。シンチレーション バイアル洗浄の間氷の上を保ちます。注意この手順は時間に依存しません、ティッシュのブロックが必要な場合は、少し長く洗浄が可能します。
    5. 冷たい重金属染色ソリューションとナトリウム cacodylate バッファーを交換し、一晩氷の上それを格納します。フェロシアン化はバイアルを手で振ることによってよく混合することを確認します。ソリューションが黒く必要があります。
      注意: は、オスミウムは有毒であるために、この手順を実行するとき、発煙のフードで動作します。
    6. 在庫 UA ソリューションの平等なボリュームと二重蒸留水を使用して 2% 4% 水溶液ウラニル酢酸 (UA) 原液 (表 1) を希釈し、氷でクールダウンします。
    7. 氷冷却 0.15 M ナトリウム cacodylate バッファー、重金属汚れに置き換えるし、氷で 5 分間インキュベート サンプルを聞かせてください。
    8. 過剰な重金属染色ソリューションを洗い流す氷の上 4 回 2.2.7 のステップを繰り返します。
    9. リンス サンプル 4 回 2 分の待ち時間の間に精製水 (ダブル蒸留で十分です)。
    10. 精製水を交換して冷えた 2% と UA とサンプル分間インキュベート 60 120 氷の上。
    11. 脱水の手順で使用する氷の上 5 の 20 mL シンチレーション バイアル (十分なサンプルが水没する) エタノールのクールダウンします。六つのバイアルがある蒸留水中のエタノールの割合が高まって: 20%、50%、70%、90%、100%。
    12. エタノール瓶と一緒に氷の上、別の 20 mL シンチレーション バイアルとクールに純粋なアセトンを注ぐ。
    13. 過剰な UA を洗浄する氷の上 4 回 2 分の待ち時間と精製水のサンプルをすすいでください。
  3. エタノールのサンプルを脱水します。
    1. 連続して氷の上、次のようにサンプル バイアル内のソリューションを置き換えることによって冷たいエタノール シリーズの脱水: 20% のエタノールの 10 分。10 分間; 50% エタノール70% のエタノールの 10 分。90% のエタノールの 10 分。その後、10 分間 2 回の 100% エタノール。
    2. 最後の 100% エタノールを冷却の純粋なアセトンに置き換えることにより部屋の温度にサンプルを移行し、10 分室温実験室ベンチにサンプル瓶を配置します。
    3. バイアル内で観察可能である結露を削除する室温で純粋なアセトンの 1 以上 10 分の洗浄を実行します。抱卵中、アセトン/樹脂溶液を作る。
    4. バイアルに純粋なアセトンを純粋なアセトンとエポキシ樹脂 (エポキシ樹脂部品メーカーによって前述の比率) を 50: 50 に置き換えます。同じバッチからすべてのコンポーネントを使用します。ローターで一晩樹脂充填サンプル瓶を残します。
  4. 樹脂に試料を埋め込みます。
    1. 75% 樹脂 (25% アセトン) 50: 50 樹脂を交換し、さらに潜伏/室温で 4 h の 100% 樹脂含浸により 3 h インキュベートします。
    2. 新鮮な 100% 樹脂 100% 樹脂に置き換え、組織サンプルに樹脂のより深い浸透のために夜通しサンプル瓶のまま。
    3. フラット ベースをバイアルからサンプルを取るおよびレンジ可の容器に置いて、アルミや銀の箔をベーキング カップはたとえばこの目的のため使用できます。
    4. (サンプルの変位を避けるため) にそっと注ぐ新鮮な (従って樹脂サンプルの埋め込み) サンプルを樹脂し、樹脂・ 48 h の 60 ° C のオーブンでサンプルを重合。
  5. 断面の画像セルに樹脂製ブロックを向き。
    1. 24セルの向きを評価するためにガラスのナイフで、ウルトラミクロトームを使用して樹脂ブロックから 1 μ m 厚テストのセクションを取得します。
    2. 20 厚テストのセクションを染色 1% トルエン ブルーと 1% ホウ砂 s。
    3. 明視野顕微鏡下で筋肉細胞の向きを確認します。
    4. 染色テスト セクションの画像から推定される方向に基づき、再向き縦または斜め指向 (図 2 aに類似) 細胞で断面 (カットされますガラス ナイフに公開されようにブロックを樹脂同様に図 2 b)。
  6. 電子線トモグラフィーを用いた画像ティッシュ セクション。
    1. 厚切切片を得る (~ 300 nm) 向きを変える組織のダイヤモンドを使用してブロックがナイフし、銅グリッド25にセクションを移動します。
    2. 2 %ua と佐藤リード25厚のセクションを染色します。
    3. セクション25の両側には、金コロイド粒子を適用します。
    4. +70 °25-70 ° から投影画像のシングルまたはデュアル軸傾斜シリーズのセットを取得します。
  7. 3 D 画像のスタック (単一電子断層レントゲン写真撮影のセクションの 3 D 情報を提供する, 二次元画像のシリーズ) を再構築する IMOD20を使用してセルの。この再建されたスタックは、次の手順でセグメント化されます。

3. セグメント筋原線維と EM 3 D 画像データセットからミトコンドリア地域

  1. IMOD プログラム"3dmod"を実行します。
  2. 3dmod グラフィカル ユーザー インターフェイス内「.rec」のアドレスを入力するか、3 D を含む".mrc"ファイル エントリ ボックスに (2.7 の手順で生成された) セルの画像データセットを再構築"画像ファイル:"、"OK"を押します。
  3. 「ファイル」メニューの下で新しいモデルを選択し、「保存としてモデル」を使用して適切な名前を持つファイルを保存...「ファイル」の下のメニュー項目。
    注: IMOD 内のオブジェクトは、「オブジェクト」を構成するセグメント化されたコンポーネントのコレクションを含めることができます。既定では、新しいモデルには id「#1」の新しいオブジェクトが含まれています。
  4. 「特別な」メニューの下には、「描画ツール」を選択します。図 3 aに示すような新しいツール バーが開きます。
    注: 異なるオルガネラの境界の周りの輪郭を作成するのに使用できるいくつかのツールがあります。これらのツールを使用する方法の詳細なヘルプは、IMOD20ドキュメントで参照できます。
  5. 「描画ツール」メニューで「スカルプト」オプションを選択し、イメージ ウィンドウへマウスを上に移動マウス ポインターを中心とした円形の輪郭が表示されます。
  6. ミトコンドリア (分界図 3Aで緑の輪郭線によって示されているように、画像ファイルの暗い領域) をマウスの中央ボタンを押し、その境界の形に円の輪郭の境界をドラッグします。
  7. ミトコンドリア境界の輪郭が完了すると、中ボタンを放すし、データで各ミトコンドリアに対して、手順 3.6。各輪郭は、同じオブジェクト内の新しい輪郭として IMOD によって自動的に認識されます。
  8. 下、「編集 |オブジェクト」メニュー、新しいオブジェクトを作成するには「新しい」を選択します。これは自動的に 1 つずつオブジェクトの合計数と、この番号を新しいオブジェクトに割り当てます。
  9. 3.6、3.7 セグメントし、筋原線維の輪郭を保存する手順を繰り返します。
  10. 同様にセルの境界を分割する手順 3.6、続いて手順 3.8 を繰り返します。
  11. 「ファイル」メニューの下でモデル ファイルを保存します。
  12. 各オブジェクトのバイナリ マスク図 4A ~ Cに示すようにグループ化に変換するコマンド ライン ウィンドウで次のコマンドを実行します。
    1. モデル「imodextract オブジェクト modelfile outputmodelfile"、"オブジェクト"が抽出するオブジェクトの番号から特定のオブジェクトを抽出します。
    2. そのオブジェクトのマスクを作成する: imodmop-マスク 255 outputmodelfile imagefile outputmask.mrc
    3. ファイルを tif スタックに変換: mrc2tif -s outputmask.mrc outputmask.tif

4. セグメント化されたコンポーネントから有限要素メッシュを作成します。

  1. Iso2mesh は、TIFF 画像のスタックを体積の四面体有限要素メッシュに変換する自由に利用できる MATLAB プログラムです。ダウンロードし、iso2mesh.sourceforge.net から iso2mesh を MATLAB のパスに追加します。
  2. ソース コードと github のウェブサイト https://github.com/CellSMB/RyR-simulator からメッシュの RyR クラスターをシミュレートするためにデータをダウンロードします。
  3. CardiacCellMeshGenerator MATLAB アプリケーション (図 5) を開始します。
  4. GUI の左上側の 3 つのプッシュ ボタンを使用して MATLAB に異なるオルガネラ コンポーネント マスクを読み込みます。
  5. 筋原繊維および図 6 aに示すように、ミトコンドリアの間のギャップを設定する別のバイナリ画像のスタックを作成します。
    1. ImageJ を開きます。
    2. ファイルを使用して |ダイアログ ボックスを開く、プログラムに筋原繊維およびミトコンドリアの tiff スタックをロードします。
    3. イメージの追加プラグインを開始を選択して"プロセス |電卓プラス」。
    4. I1 として原画像のスタックを選択し、ミトコンドリアは画像のスタックと i2 は、「追加」演算子を選択します。「OK」をクリックします。
    5. 4.5.4 の結果を表す新しい画像のスタックが表示されたら、選択"編集 |図 6Aに似た画像のスタックを生成する"を反転します。
  6. CardiacCellMeshGenerator プログラムの"RyRGapsFile"ボタンを押して筋原繊維およびミトコンドリアの間のギャップのバイナリ画像のスタックを含むファイルをロードします。
  7. GUI の「メッシュを生成」を押します。これは図 4 eに類似した四面体メッシュを生成する 'cgalmesh' オプションを指定して iso2mesh コマンド v2m がトリガーされます。3 つのファイルがこの手順の終わりで出力されます: .ele ファイル、.face ファイル、および要素、面を構成するノードを構成するノードの一覧、およびノードの座標を含む .node ファイル、それぞれ。

5. 数学的に有限要素メッシュに関心のイオン チャネルの空間的に変化の密度をマップします。

  1. RyR シミュレータの GUI の生成 RyR シミュレータ入力のラベル ボタンを押して必要な入力を生成します。
    注: ボタンは、ステップ 4 から次の情報を使用して、入力を生成する関数 generateRyRsimulatorInputs.m をトリガーします。
    (1) outDir: RyR クラスター シミュレーションに必要な出力ファイルの場所。
    (2) imres: 画像のスタックの 3 つの方向のピクセルの解像度。
    (myofibril_file 3):図 4 bに示すようなバイナリ画像のスタックを含むファイル。
    (sarcolemma_file 4): 4 a を図に示すようなバイナリ画像のスタックを含むファイル。
    (ryrgaps_file 5):図 5 aに示すようなバイナリ画像のスタックを含むファイル。
    1. この関数を実行すると、outDir パスとして指定されたディレクトリ以下のファイルが作成されていることを確認します。
    • d_axial_micron.txt は、z ディスクの位置と画像のスタック内のピクセルの残りの部分の軸方向の距離を表します。
    • d_radial_micron.txt、RyR クラスターの可能な場所のセット内の各ピクセルから z ディスク平面上のピクセルにユークリッド距離 (軸のコンポーネントを除く) を表す。
    • W_micron.txt は、RyR クラスターに存在するすべての利用できる位置の空間座標のリストを表します。
    • フォルダーに残りの 3 つのファイルは、これらのファイル内の値はピクセル座標形式書かれたことを示すために"_micron"ではなく接尾辞"_pixel"を含まれています。
  2. 筋原繊維のバイナリ画像のスタックの RyR クラスター分布をシミュレートします。
    1. 「R でオープンの RyR シミュレータ」のラベル ボタンを押して R プログラムを開始します。
    2. R-gui で選択"ファイル |オープン」と「設定ファイルを見つけるR"RyR シミュレーター パッケージ (RyR-シミュレータ/ソース/設定します。R)。
    3. またファイル ryr シミュレータ並列を開きます。R (RyR-シミュレータ/ソース/ryr-シミュレータ-平行に位置します。R). Rstudio (https://www.rstudio.com) またはプレーンのコマンド ライン インターフェイスを使用できるような環境では、例えばシェルまたはコマンド ウィンドウで R64 コマンドを使用します。
    4. 設定のパラメーターを変更します。R 以下のファイルのとおり:
      1. Path2 を RyR クラスターと筋原線維の実験的に得られた共焦点画像のスタックの 5.1.2 に示すファイルが含まれているフォルダーのアドレスに設定します。
        注: github のリポジトリ フォルダー入力ファイル/マスター-細胞/以前に収集した画像のスタックの生成されたファイルが既に含まれています。
      2. 5.1.2 のステップで生成されたファイルが格納されているフォルダーのアドレスに Path4 を設定します。
      3. Path3 ポイントをユーザーが保存するシミュレートされた RyR クラスターの場所を望んでいるフォルダーに設定します。
      4. (通常は 200 ~ 300 の範囲) モデルでシミュレートする RyR クラスターの数 N のセット。
      5. Etol RyR クラスター分布測定と RyR クラスター分布シミュレーション モデルの違いの設定許容範囲を設定します。
      6. RyR シミュレータ必要があります etol 値を満たす模擬 RyR クラスター パターンを見つけるための試行回数を制限する numIter を設定します。
        注: etol および numIter の値は、設定内の典型的な値に設定されています。R ファイル。
      7. NumPatterns をシミュレートするユーザー別の RyR クラスター パターン数に設定 (通常、統計的信頼性を高めるのための 99 のパターンをシミュレートする有用な練習なのです)。
      8. 高速並列処理点パターンをシミュレートするために R を複数の CPU スレッド (コア) の使用を有効にするのに numCores を設定します。
    5. 以下のパッケージは、パッケージ インストーラーを使用してインストールされていることを確認してください r: 雪、doSNOW、doparallel、foreach、反復子、および保険の gui。
    6. R のコマンド ウィンドウで次のコマンドを入力して、シミュレータを実行: ソース (' ryr シミュレータ並列へのパス。R'、chdir = TRUE)
      メモ: RyR シミュレータ プログラムの出力 (numPatterns ファイルが生成されます) N 行での座標のリストを含む .txt ファイルの一覧では、x、y、および z 座標の n を表す 3 つの列、シミュレート RyR クラスター。
  3. CardiacCellMeshGenerator を用いた計算モデル上に空間密度としてポイントをマップします。
    1. ボタンのラベルを選択して「RyR ポイント ファイルの選択」、RyR シミュレータによる出力をしたものから模擬 RyR クラスター分布のテキスト ファイルを選択します。
    2. MATLAB で GUI で「RyR 密度マッパー」を実行します。これは、4.8 球状カーネルの強度推定法26と呼ばれる方法を使用しての手順で生成された有限要素メッシュ歩 5.3.1 で .txt ファイルでシミュレートされた RyR クラスターの空間位置をマップします。
      注: この手順の出力は各メッシュ ノードで球状体積比イオン チャネルの数値の密度の値の一覧を含む .txt 拡張子を持つファイルです。

Subscription Required. Please recommend JoVE to your librarian.

Representative Results

図 2から図 7は、このプロトコルのいくつかの主要ステップの代表の結果を提供する: (i) の可視化と組織ブロック断面電子顕微鏡ビューの方向(ii) 生成する 3次元電子顕微鏡画像のスタック;(関心のオルガネラの細胞内微細構造 iii) 分割(iv) iso2mesh; を使用して有限要素メッシュの生成(v) 現実的なメッシュの RyR クラスター分布をシミュレートします。(vi) 計算メッシュ空間密度としてこの空間分布をマッピングが続きます。

図 2縦方向のセルの表示方法を示す典型的な明視野画像が表示されます (図 2 aL のラベル)、斜め (図 2 aO をラベル) と重し (図 2 b) に関して、切断面。重し方向づけられたサンプルが縞を示さない間、斜めと縦方向のサンプルは縞を展示します。毛細血管もより円形よりも斜めビューの 2 次元断面ビューで表示されます。

図 3 aは、手順 1 で組織の準備のプロトコルが続く時に獲得できる質の良い断層スタックの代表的なイメージを示しています。顕微鏡画像の傾きシリーズから 3次元復元を実行するときに注意が必要があります。金コロイドの粒子の追跡または十分な金粒子の欠如の品質に影響することができます、画像のコントラストやフォーカス-断層のかなり。これは、大幅に異なる構造を分割する機能を影響します。ケアと経験がティッシュのブロックを十分に汚れるし、組織ボリュームを金コロイド粒子の均一な配分があることを確保するため必要です。疑問がある場合は、専門的な研究や現地機関で施設に相談すると便利があります。モデル生成と低コントラストを試みてデータを正しく再構築、モデルを生成が可能性がありますが、モデリングをシミュレートするのに生物学的システムのプロパティと結果の対応が悪い場合があります。図 3 bは、筋原繊維とミトコンドリアの分割模型のように見えるものの代表的なイメージを示しています。

図 4E iso2mesh によって生成される四面体メッシュの 3 D レンダリングを示しています。良質のオリジナル画像スタックとセグメンテーション タスクは、限り、この手順はかなり堅牢です。図 6B図 6は、3 D の細胞のトポロジ内にある (赤) は典型的な RyR クラスター分布を示します。この段階で、メッシュ自体は RyR シミュレータでの入力ではありません。セルの境界と筋原繊維およびミトコンドリアの間のギャップの画像のスタック-図 4に分割された画像から生成された-RyR シミュレータに主要な入力を形成します。筋原繊維およびミトコンドリアの境界をできるだけ; 正確にセグメントには注意が必要その結果メッシュ内 RyR 配置に影響されます、データのセル トポロジの不正確な表現で不正確な区切りになります。この、カルシウム (アプリケーション 1) 心筋細胞におけるシグナル伝達の時空間的ダイナミクスに影響を与えるでしょう。

図 7は、メッシュ ・ トポロジーに球状カーネル強度推定アルゴリズムを使用した後同じ模擬 RyR クラスターの 3 つのインスタンスを示しています。RyR クラスターの組織の変化に注意してください。実験データ18で発見された変動内にあるこれらの変化を測定しました。カーネルの球の半径とどのカーネル サンプル RyR クラスター密度ステップ サイズを選ぶ際に注意する必要があります。これらの 2 つの要因に影響を与える急激またはメッシュ上びまん性 RyR クラスターの密度の表現が描かれています。これらのパラメーターの値が異なるの効果の分析は合理的なメッシュのパラメーターの選択肢に絞り込むのに役立ちます。

RyR クラスター、筋原繊維とミトコンドリアの結果のメッシュは、心筋細胞構造の最小モデルです。このメッシュはカルシウム動態と生体エネルギーを勉強する他の細胞構造データから派生したバリエーションだけでなく、適用されています。これらの 2 つのアプリケーションの概要構造と機能の関係に洞察力を与えることで細胞構造の有限要素モデリングの力を説明するために以下の通りです。

アプリケーション 118: 心筋細胞にカルシウムの時空間ダイナミクスを解放します。このモデルは、RyR クラスター、筋原繊維とカルシウム シグナルのミトコンドリアの不均質分布の役割を調べるために使用されています。図 8は、このプロトコルから生成されるメッシュ トポロジ上でシミュレートするときの一時的な上昇段階ゾル性細胞質カルシウムの時空間動態を示します。図 8 bは、異機種混在の自由なカルシウムと蛍光体結合カルシウム濃度の経年を示しています。矢印は、RyR クラスターまたはそれらの領域におけるミトコンドリアの密度が高いため高カルシウムの小さな斑点をポイントします。図 8は、通常ライブの共焦点顕微鏡を使用して収集されるカルシウム動態の実験ライン スキャン画像と比較することができますシミュレートされたライン スキャン画像を示しています。モデル シミュレーションは実験的に得られた画像よりもカルシウムの空間的不均一性が多く高い解像度表示されます。さらに、構造顕微鏡データから派生しているモデルのためのままにする必要がありますクラスター数不活化 (4 〜 6) 共焦点ライン スキャン画像で発生する異種の電気的活性化を的確に判断することができる後。似たような不均質、心不全18の動物モデルにおけるカルシウム動態の実験的に得られたライン スキャン画像で観察されます。

アプリケーション 227: 心筋ミトコンドリアのエネルギー代謝の時空間ダイナミクス。この作品の包括的な目的は、カルシウム動態、ミトコンドリアの生体エネルギーと筋収縮との関係を研究することです。最初のステップとして分離、心筋細胞内のミトコンドリアの酸化的リン酸化のシミュレーションが行われています。RyR クラスターの分布が必要でないし 1 つは単に手順 3.4 でメッシュを使用できます。複雑なアプリケーションとして生体エネルギー (応用 2) がカルシウム動態 (アプリケーション 1) メッシュ モデルの構築、すなわちミトコンドリア、筋原繊維とリアノジンと見なされますすべてのコンポーネントを使用するモデルと結合できます。受容体。

図 9は、3.4 の手順を使用して生成されたメッシュのさまざまな地域でこれらの空間的代謝シミュレーションから結果を表します。図 9 aは、図 9 b細胞の筋原線維コンパートメントでのみ ADP の濃度を示していますセル断面全体の酸素濃度を示しています。図 9は、ミトコンドリアのネットワークを介してミトコンドリア膜電位の分布を示しています。これらの数字は、ミトコンドリアと細胞断面における筋原線維の非一様な分布を明らかにします。彼らはまたこの不均一微細構造から生じる不均質代謝風景を明らかにします。これは、顕微鏡データの力に基づいて有限要素シミュレーションだけで実験的手段により取得することは困難ではない構造と機能の関係に洞察力を与えることを示しています。

Figure 1
図 1: 化学固定のため心の蛍光準備します。(A) 蛍光灌流装置の模式図。活栓にはタイロードを切り替えるバルブが含まれますと固定液。基地でビーカーを使用して、心臓を灌流するソリューションをキャッチします。この図の拡大版を表示するのにはここをクリックしてください

Figure 2
図 2: 明視野顕微鏡下で細胞の向きを調べるします。(A) 縦方向 (L) と画像平面に対して斜め方向 (O) のセルを示す心筋のトルエン青いステンド厚いセクションの明るい視野たとえば縦配向セルの痕を注意してください。(B) 画像の平面に揃えセルの横方向の断面組織切片の明るい視野。この図の拡大版を表示するのにはここをクリックしてください

Figure 3
図 3: 電子顕微鏡データのセグメンテーション。(A) マニュアルでのスカルプト ツールを使用してミトコンドリアのセグメンテーションのスナップショット。(B) A の完全な 3 D ビューには、ミトコンドリア、(いろいろな色) の筋原線維筋鞘 (ピンク) で手動で分割されます。この図の拡大版を表示するのにはここをクリックしてください

Figure 4
図 4: 3 D の四面体メッシュ バイナリ画像のスタックに変換する。(A ~ C)ミトコンドリア ・筋原繊維筋鞘のバイナリ マスクそれぞれ。(D) 合併 (赤) は筋原繊維および (グリーン) におけるミトコンドリアのバイナリ マスクの表示します。(E) (D) の画像のスタックが iso2mesh を使用してから生成された代表・四面体・ メッシュ。この図の拡大版を表示するのにはここをクリックしてください

Figure 5
図 5: CardiacCellMeshGenerator.手順 4 と 5 を実行するユーザーのためのこのプロトコル文書に付随するグラフィカル ユーザー インターフェイスのスクリーン ショット。この図の拡大版を表示するのにはここをクリックしてください

Figure 6
図 6: 心臓微細構造のセグメント化された電子顕微鏡画像のスタック上に現実的な RyR クラスター分布をシミュレートします。(A) のバイナリ画像のスタックが筋原繊維および RyR クラスターを見つけることがすべてのピクセルのセットを表すミトコンドリアの間のギャップを描いた.(B) 3 D レンダリング (イメージ スケールを描画されません) 模擬 RyR の 1 つのクラスターの 3 D 画像スタックに示されている 3 D ジオメトリ内の分布 (赤い球体として示されて)グリーンの表面は、画像のスタック内のミトコンドリアを表しています。(C) A 高倍率 (従って小さい部分を切り出してセル断面の境界で)、RyR のオーバーレイのクラスター計算メッシュを 3D 電子断層レントゲン写真撮影の画像のスタックの 1 つのスライスにステップ 3 からステップ 4 から。(B) および (C) Rajagopalから変更されています。18この図の拡大版を表示するのにはここをクリックしてください

Figure 7
図 7: 3 から RyR クラスターの数値密度の空間分布シミュレーション RyR クラスター パターン。すべてのメッシュ ノードが数値密度 0.0 RyR クラスター/μ m3のための白から 0.99 RyR クラスター/μ m3の最大数値密度の赤に色分け。この図の拡大版を表示するのにはここをクリックしてください

Figure 8
図 8: [Ca2 +]i Ca2 +一過性の上昇段階で不(A) 自由に拡散 [Ca2 +](左) と蛍光色 4 ゾル性細胞質の平均がバインドされた Ca2 +[F4Ca](右)。(B) および (C) 横断面; z ディスクで Ca2 +と F4Ca を自由に拡散の時間経過のスナップショットを表示黒い矢印は、高 [Ca2 +]ミトコンドリア クラスタ リングのための領域を示す赤い矢印は、近くに高 [Ca2 +]iいくつかの RyR クラスターのための領域を示します。(C) 模擬ライン スキャン [Ca2 +](中央) [F4Ca](右)。行の位置は、左側のセル断面に表示されます。セル断面の水平方向の寸法は画像の断面の空間スケールの指標として提供されています。この図は、Rajagopalから変更されています。18この図の拡大版を表示するのにはここをクリックしてください

Figure 9
図 9: 心臓エネルギー新陳代謝の空間のシミュレーションから生成される代謝とミトコンドリア膜電位の空間分布。(A) μ M の単位でセル断面全体の酸素濃度の単位 μ m で細胞の筋原線維部分だけの ADP 濃度の (B) (C) ミトコンドリア膜、ミトコンドリアの間で潜在的な分布mV の単位でネットワーク。この図の拡大版を表示するのにはここをクリックしてください

Figure 10
図 10: R 統計パッケージにシミュレートされた RyR クラスター ポイント パターンの可視化します。このウィンドウは、シミュレートされた点パターンの最初を描いた RyR クラスター シミュレーションの結論に表示されます。これは、RyR シミュレータなどを締結しているシミュレートされたパターンの品質を問い合わせる場合があることをインジケーターとして使用できます。この図の拡大版を表示するのにはここをクリックしてください

化学物質 数量
0.2 N 塩酸
1 N の HCl 5 mL
蒸留水 20 mL
0.15 M ナトリウム cacodylate 緩衝在庫
ナトリウム cacodylate 32.1 g [3 H2O]
0.2 N 塩酸 25 mL
蒸留水 1 L を作る
pH 7.4
0.3 M ナトリウム cacodylate 緩衝在庫
ナトリウム cacodylate 64.2 g [3 H2O]
0.2 N 塩酸 25 mL
蒸留水 1 L を作る
pH 7.4
Tyrode 液組成
塩化ナトリウム 139 mM
KCl 3 mM
CaCl2 3 mM
MgCl2 1 mM
グルコース 12 mM
NaHCO3 17 mM
プロベネシド 1 mM
BDM 20 mM
NaCl、KCl、NaHCO3 1 L ソリューション
塩化ナトリウム (分子量 58.44 g/mol) 81.2 g
KCl (分子量 74.55 g/mol) 2.24 g
NaHCO3 (分子量 84.01 g/mol) 14.28 g
蒸留水 1 L の水に溶かす
MgCl2/CaCl2 1 L ソリューション
CaCl2 (2H2O) (分子量 147 g/mol) 1.764 g
MgCl2 (6H2O) (203.31 g/mol) 0.814 g
蒸留水 1 L の水に溶かす
1m グルコース 200 mL ソリューション
ブドウ糖 (分子量 180.16 g/mol) 36.03 g
200 mL 二重蒸留水に溶解します。
Tyrode ソリューション 500 mL ソリューション
NaCl、KCl、NaHCO3 50 mL
グルコース 6 mL
二重蒸留水 400 mL
5 分のためのソリューションでバブル酸素
2 %pfa + 2.5% グルタルアルデヒド 0.2% タニック酸の定着剤 0.15 M ナトリウム cacodylate
粉末パラホルムアルデヒド (PFA) 2 g
酸にタニック 0.2 g
蒸留水 20 mL
25% グルタルアルデヒド 10 mL
0.3 M ナトリウム cacodylate 20 mL
0.15 M ナトリウム cacodylate 50 mL
1 N NaOH 4 グラム NaOH/100 mL の蒸留水
4% ウラニル酢酸原液
酢酸ウラニル 4 g
二重沸騰近く蒸留水 96 mL
化学で、撹拌機を使用して bioling 近く水で UA を溶解します。
光保護のため 200 mL 茶色のボトルに刈り取る #1 フィルターを介してフィルター UA
しっかりとキャップしラベル、4 ° C で保存

表 1: 試薬と電子顕微鏡検査のためのティッシュの準備の量。

Subscription Required. Please recommend JoVE to your librarian.

Discussion

上記のプロトコルでは、心筋細胞微細構造の新しい有限要素の幾何学的モデルを生成する重要な手順について説明します。メソッドにより空間細胞アーキテクチャの詳細を含む心筋細胞動態のより包括的な計算モデルを開発する各種の顕微鏡 (または、原則として、その他のデータ) 計算の融合療法です。現在では、心筋細胞のようなモデルの作成に利用可能な他のプロトコルはありません。

ステップ 1 では、灌流固定のためのプロトコルについて説明します。また、心臓摘出より小さな部分に解剖でき固定液に浸漬します。ただし、このプロセスはすぐに行う必要があります、サンプルのサイズによって結果を変えることが。筆者の経験では、灌流は組織や細胞固定液の徹底した浸透を確認します。著者も以前中心にクレブス Henseliet ソリューションと心臓を灌流固定前にビート。これは、固定する前に作業中心の測定を有効になります。

ステップ 2 の電子顕微鏡検査処理手順は電子線トモグラフィーの典型的であります。シリアル ブロック面走査電子顕微鏡、イオンビーム顕微鏡28を焦点を当てて、他のイメージング技術を使って染色ソリューション、およびタイミングの特に処理ソリューション数量が異なる場合があります。取得の電子顕微鏡画像は、ミトコンドリア、筋原線維、横管などの異なるオルガネラ領域に分割する必要があります。これは、手動または自動化された方法で実施できます。画像のスタックのコントラストは汚損のプロトコルの成功に部分的に依存してが、断層レントゲン写真撮影の際、このコントラストのいくつかが失われる。このプロトコルは、自動分割29を実行する新しいメソッドは、データのスループットを改善するために一般により有利になりますが、異なる構造のセグメント化のための手動ステップを記述します。

電子顕微鏡データセット内を調べる最も重要な資質は、画像のコントラストと構造の保護です。染色レシピこのプロトコルへの混合物の埋め込み樹脂筋原線維、ミトコンドリア心筋細胞の z ディスクとハイコントラストの洗練されたきた。たとえば、著者のステップ 1.2 で筋の網状のネットワークを区別するためタンニン酸の代わりに塩化カルシウム使用している以前。これは細胞小器官セグメンテーション中に原束の外側の境界を識別する助けた。追加染色化学成分および回も筋小胞体や横行軸細管30などの膜構造を可視化するため使用されます。

画像は、構造上の損傷の検査されなければなりません。質の悪い固定または電子顕微鏡による処理は、ミトコンドリアのクリステの深刻な損失、剥がれや腫れのミトコンドリアの割合が高いと細胞膜 (筋鞘) の部分の崩壊で起因できます。近いと慎重な検討はまた t 鋼管の損失を表示ことがあります SR 膜が、素人目にはわかりにくいのと。この問題の解決策があります: (i) 確保徹底した灌流;(ii) 小さなサンプルも定着剤の汚れの浸透を確実に組織を解剖します。(iii) 電子顕微鏡処理こと氷またはコールド」ソリューションで指定場所を確保します。(iv) 染色と脱水 (特に脱水) のタイミングが厳密に続いていることを確認します。

典型的な心筋細胞は、断面直径 ~ 20 μ m、長さ ~ 100 μ m です。これは画像全体の心臓の細胞を重要な課題になります。さらに、心31更に内筋線維の変化の方向は軸の興味のセルの横と縦の軸と一直線に並ぶイメージ ボリュームの買収を複雑にします。IMOD など画像解析プログラムは、所望の方向にデータの合成の再配置を有効にします。シリアル ブロック面走査電子顕微鏡28と関連する画像処理アルゴリズムの29の最近の進歩は、これらの問題を解決すること。それにもかかわらず、光学顕微鏡下で厚いセクションの入念な検査およびブロック (ステップ 2.5) の後続の再配向イメージングと縦・横の軸の合理的な配置と細胞のイメージングの可能性が高くなります軸。これは細胞の容積のほとんどをビューのフィールド内でキャプチャされることが保証されます。

プロトコル手順 3、4 および 5 をインストールする材料の表に記載されているソフトウェアが必要です。IMOD、R 統計パッケージ、iso2mesh、フィジー (顕微鏡データの ImageJ のバージョン) それらを使用する方法で広範なマニュアルとチュートリアルがあります。CardiacCellMeshGenerator は、単純なワークフローのモデルを構築するユーザーを容易に開発されている MATLAB アプリケーションです。ユーザーは、全体のソース コード、およびアプリケーション パッケージを github RyR シミュレータ リポジトリ リンク (https://github.com/CellSMB/RyR-simulator/tree/JoVE) からダウンロードできます。アプリケーション、CardiacCellMeshGenerator.miappinstall は、MATLAB のアプリケーションとしてインストールできます。このプロトコルでは生成されるメッシュの入力見つけることができます RyR-シミュレータ/入力-ファイル/ターゲット-トモ-セル内/追加入力中 RyR-シミュレータ/入力-ファイル/マスター-セル内 RyR クラスター シミュレータを発見することができます/。アプリケーションを開始する前にユーザーがインストールされている iso2mesh フォルダーとそのサブフォルダーを MATLAB パスに追加されていることを確認してください。同様に、r. 内もクラスター シミュレータのための R パッケージ必要 (材料表) がインストールされていることを確認します。

RyR 密度マッパー; を使用してメッシュの RyR クラスター密度をマッピングするとき RyR クラスター シミュレータの入力を生成するときにメッシュを生成するとき、進行状況バーがアプリに含まれています。メッシュ生成のステップが完結したとき、アプリケーションは GUI で 3 D メッシュをレンダリングします。R ユーザー インターフェイスの RyR クラスター シミュレーションを完了すると、グラフィカルなウィンドウが表示されます (図 10)、シミュレートされた 3 D RyR クラスター ポイント パターンの一つを描いたします。

GUI の「メッシュ生成」ボタンは、四面体メッシュを生成する iso2mesh 関数をトリガーします。モデルの解像度は、最大の四面体要素のボリュームとエッジの長さを調整して GUI で「iso2mesh メーター」を設定することによって調整できます。プログラムは現在原地域およびミトコンドリア地域を構成する要素を区別します。この機能は、将来的にセルの他のコンポーネントを含めるに拡張されます。

ユーザーは、次の手順をトラブルシューティングするのにメッシュおよび RyR クラスターの 3 D レンダリングを使用できます。RyR クラスター シミュレーション設定、etol および設定中の numIter。R は、RyR クラスター シミュレーション; の精度に影響します。これらの 2 つのパラメーターは、RyR クラスター パターン シミュレーションを終了するときを決定します。確認する [R] ウィンドウ (図 10) RyR クラスター分布の検査時にこれらのパラメーターを調整ことができます: (i) すべてのクラスターが z ディスクの近くにあります(ii) クラスターに分散して、1 つの領域に集約するのではなく、全体の断面図。Etol と numIter の適切な組み合わせを決定するため、これらの 2 つのパラメーターをシミュレーションのアドホック感度解析を実施できます。

手順 5 で球面カーネル強度推定は、興味 (このコンテキストでは 3 D 空間モデル) のボリュームを渡されるカーネルとして選ばれた直径の球を使用します。これらのパラメーター調整できるこのため図 5の GUI 上でアルゴリズムを平滑化: (1) r_sphere は、球状のカーネルのスムージングの半径を定義します。(2) hval は、球状のカーネルをボリュームの繰り返し渡される必要があります反復空間ステップ サイズを定義します。

空間密度パラメーターが十分であるかどうかを判断する簡単な手法は、図 7に示すように、RyR クラスター密度値の分布を検査することです。高密度値の近所が顕微鏡データ RyR クラスターのサイズと同様のサイズになるよう、密度値を分散する必要があります。

現在のプロトコルは、現実的な筋原線維、ミトコンドリア、および RyR クラスター分布のみを含む有限要素モデルを生成します。このステップへの拡張は、サルコ小網様体カルシウム ポンプ (SERCA2a)、ナトリウム カルシウム交換器 (NCX), カルシウム ポンプなど、ECC の役割を果たすイオン チャンネル他の分布をシミュレートするでしょう。これらの拡張機能にキーは、蛍光抗体法データからイオン チャネルの空間分布の分析です。たとえば、0.67 μ m の平均間隔で t 管状の膜に均一に分散する NCX チャンネルが測定されているが、RyR クラスターとの関係ではないのでクリア32。その一方で、SERCA2a ターゲット抗体は以前 SR33をマークに使用されているが、SERCA2a の空間密度は文献で報告されていません。したがって、この段階で注意して SR ネットワークを介して SERCA 2 a と t-細管上 NCX の均一な分布を想定可能性があります。

有限要素ソルバー34結果有限要素メッシュをカルシウム反応拡散として18と生体エネルギーの27などの生化学的なプロセスをシミュレートする使用ことができます代表結果に示されているように処理します。このため、イオン チャンネルは、メッシュ トポロジ内のノードの反応サイトとして表されます。反応サイトは、ソースまたはセルの別のコンパートメントにこれらのイオン チャネルを通しての物質フローを表すための別の化学種のシンクとして機能します。出力メッシュとこの議定書から RyR 密度ファイルは使用ソフトウェアは、これらのテキスト ファイルを認識する機能を持っている他のインテグレーターまたは他の環境の中でをすることができます。テキスト ファイルは、スクリプトや選択の環境の入力ファイルの形式要件を満たすサード パーティ製ソフトウェアによって操作できます。カルシウム動態に最初のアプリケーションのみカルシウム脱分極相のこのモデルの有用性を示しますが、長い時間スケールでのメッシュを使えないというわけではないです。メッシュとクラスター密度は、恒常性と使用される有限要素ソルバーにのみ依存して安定性の維持と、複数のカルシウム サイクルをシミュレートするために使用できます。

長期的な目的は、このパイプライン モデル建物の心筋細胞と他の細胞型の構造的に高精度有限要素モデルを作成しやすくを改善するためにです。うちは最小限のユーザーと全体の細胞の筋原繊維およびミトコンドリア分布のモデルを作成する心筋細胞のシリアル ブロック面走査電子顕微鏡データに適用できる自動セグメンテーション アルゴリズムは最近、開発されました。入力29。このメソッドは、将来の t-細管と SR 膜ネットワークをセグメント化する拡張されます。分析およびイオン チャネルと細胞小器官の間だけでなく、別のイオン チャネル型、高スループット方法間の共同の場所をシミュレートするユーザーを有効にする RyR シミュレータ アルゴリズムのより一般的なバージョンは、また開発中です。セルの完全なモデルを構築するためのシームレスなプロトコルと科学者、実験的アプローチを現在よりも日常的により多くの細胞の構造と細胞機能の関係に関するモデルに基づく仮説テストを実行するを有効にします。一人で。

Subscription Required. Please recommend JoVE to your librarian.

Disclosures

著者は、彼らは競合する金銭的な利益があることを宣言します。

Acknowledgments

この作品は、ロイヤル会のニュージーランド マースデン高速起動許可 11-UOA-184、人間フロンティア科学プログラム研究助成 2013 年 RGP0027 月とオーストラリアの研究評議会発見プロジェクトの助成金 DP170101358 によって支持されました。

Materials

Name Company Catalog Number Comments
Materials
Sodium chloride Sigma-Aldrich 746398
Calcium chloride Sigma-Aldrich C8106
Magnesium chloride Sigma-Aldrich M2393
Sodium bicarbonate Sigma-Aldrich S5761
Potassium chloride Sigma-Aldrich P5405
Dextrose Sigma-Aldrich D9434
Sodium hydroxide Sigma-Aldrich S8045
Probenecid Sigma-Aldrich P8761
2,3-Butanedione monoxime Sigma-Aldrich B0753
25% Glutaraldehyde EM Grade (500 ml bottle) Merck 354400-500ML
Paraformaldehyde Sigma-Aldrich P6148
Tannic Acid Sigma-Aldrich 403040-500G 100g EM grade
Sodium cacodylate Sigma-Aldrich C0250
Phosphate-buffered saline (PBS) Sigma-Aldrich P4593
Osmium Tetroxide Sigma-Aldrich 75632-10ML 4% in water, 5 ml bottle (or 10 ml bottle also available)
Uranyl Acetate EM Sciences 22400 25g bottle
Potassium Ferrocyanide Merck Millipore 104973
Toluene blue Sigma-Aldrich T3260
Borax Sigma-Aldrich S9640 also termed sodium borate
Ethanol Sigma-Aldrich 792780 Diluted to different percentages with pure water
Acetone EM Sciences RT10017
Resin kit EM Sciences 14040 ACM Durcupan works well
Hydrochloric acid Sigma-Aldrich H9892 1Normal solution
Equipment
Ultramicrotome Leica EM UC7
Transmission electron microscope ThermoFisher Scientific Tecnai F30 http://www.leica-microsystems.com/
Retort stand Proscitech T752
Tubing BioStrategy 75831-346 for langendorff perfusion apparatus, 3 mm diameter is recommended but not essential
Stopcocks SDR QP13813 for langendorff tubing; product is only an example, user can select any
retort stand clamps Proscitech T715
Plastic syringes SDR QPC1108 for solutions on langendorff apparatus
Cannulation silk suture, 7-0 TeleFlex 15B051000 for tieing heart on langedorff apparatus
Cannula Made from 3 mm outer-diameter steel needle
Rubber petri dish mat Proscitech H068 for use as cutting board during fixed-heart dissection
Razor blades Proscitech L056 for cutting fixed-heart into small blocks for EM processing
Glass bottles BioStrategy 89000-236 for storing solutions during tissue fixation and processing for EM
Beakers BioStrategy 213-0477 for storing solutions temporarily and during perfusion
Scintillation vials BioStrategy 548-2170 for tissue samples during EM processing
Dissection kit Proscitech T161 for animal dissection
Syringe Filters Proscitech WS3-02225S for purification of Uranyl Acetate
Aluminium/silver foil baking cups From any baking products store
Dupont Diamond knife BioStrategy 102680-780 35 degree angle version produces best sections.
Colloidal Gold BBI Solutions EM. GC15 15 nm colloidal gold
EM mesh grids Proscitech GCU150 a variety of sizes can be tested: GCU150h, GCU200h for example
Plastic disposal pippettes Proscitech LCH20 best to use plastic disposables especially when working with resin
Software
SerialEM University of Boulder tomography acquisition
MATLAB MathWorks https://www.mathworks.com/products/matlab.html
IMOD University of Boulder image alignment and segmentation
iso2mesh available at http://iso2mesh.sourceforge.net
Fiji or similar image processing software ImageJ Fiji is Just Image J available at https://fiji.sc for manipulation of binary image stacks
RyR-Simulator codes/data CellSMB group available at https://github.com/CellSMB/RyR-simulator
CardiacCellMeshGenerator CellSMB group comes with RyR-Simulator under folder "gui-version"
R-statistics software R-project Download from https://www.r-project.org
spatstat R-project install via R program
rgl R-project install via R program
doparallel R-project install via R program
foreach R-project install via R program
doSNOW R-project install via R program
iterators R-project install via R program

DOWNLOAD MATERIALS LIST

References

  1. Noble, D., Rudy, Y. Models of cardiac ventricular action potentials: iterative interaction between experiment and simulation. Phil. Trans. R. Soc. A: Math., Phys. and Eng. Sci. 359 (1783), 1127-1142 (2001).
  2. Williams, G. S. B., Smith, G. D., Sobie, E. A., Jafri, M. S. Models of cardiac excitation-contraction coupling in ventricular myocytes. Math. Biosci. 226 (1), 1-15 (2010).
  3. Beard, D. A., Vendelin, M. Systems biology of the mitochondrion. Am. J. Phys. - Cell Phys. 291 (6), C1101-C1103 (2006).
  4. Crampin, E. J., Smith, N. P. A Dynamic Model of Excitation-Contraction Coupling during Acidosis in Cardiac Ventricular Myocytes. Biophys. J. 90 (9), 3074-3090 (2006).
  5. Li, L., Louch, W. E., et al. Calcium Dynamics in the Ventricular Myocytes of SERCA2 Knockout Mice: A Modeling Study. Biophys. J. 100 (2), 322-331 (2011).
  6. Shimizu, I., Minamino, T. Physiological and pathological cardiac hypertrophy. J. Mol. Cell. Cardiol. 97, 245-262 (2016).
  7. Wei, S., Guo, A., et al. T-tubule remodeling during transition from hypertrophy to heart failure. Circ. Res. 107 (4), 520-531 (2010).
  8. Jarosz, J., Ghosh, S., et al. Changes in mitochondrial morphology and organization can enhance energy supply from mitochondrial oxidative phosphorylation in diabetic cardiomyopathy. Am. J. Phys. - Cell Phys. 312 (2), C190-C197 (2017).
  9. González, A., Ravassa, S., Beaumont, J., López, B., Díez, J. New Targets to Treat the Structural Remodeling of the Myocardium. J. Am. Coll. Cardiol. 58 (18), 1833-1843 (2011).
  10. Hayashi, T., Martone, M. E., Yu, Z., Thor, A., Doi, M. Three-dimensional electron microscopy reveals new details of membrane systems for Ca2+ signaling in the heart. J. Cell Sci. , (2009).
  11. Soeller, C., Crossman, D., Gilbert, R., Cannell, M. B. Analysis of ryanodine receptor clusters in rat and human cardiac myocytes. Proc. Natl. Acad. Sci. 104 (38), 14958-14963 (2007).
  12. Soeller, C., Baddeley, D. Super-resolution imaging of EC coupling protein distribution in the heart. J. Mol. Cell. Cardiol. 58 (1), 32-40 (2013).
  13. Yu, Z., Holst, M. J., et al. Three-dimensional geometric modeling of membrane-bound organelles in ventricular myocytes: bridging the gap between microscopic imaging and mathematical simulation. J. Struct. Biol. 164 (3), 304-313 (2008).
  14. Hake, J., Edwards, A. G., et al. Modelling cardiac calcium sparks in a three-dimensional reconstruction of a calcium release unit. J. Physiol. 590 (18), 4403-4422 (2012).
  15. Soeller, C., Jayasinghe, I. D., Li, P., Holden, A. V., Cannell, M. B. Three-dimensional high-resolution imaging of cardiac proteins to construct models of intracellular Ca2+ signalling in rat ventricular myocytes. Exp. Physiol. 94 (5), 496-508 (2009).
  16. Kekenes-Huskey, P. M., Cheng, Y., Hake, J. E. Modeling effects of L-type Ca2+ current and Na+-Ca2+ exchanger on Ca2+ trigger flux in rabbit myocytes with realistic t-tubule geometries. Front. in Physiol. 3, 1-14 (2012).
  17. Zienkiewicz, O. C., Taylor, R. L. The finite element method. 1, Butterworth-Heinemann. (2000).
  18. Rajagopal, V., Bass, G., et al. Examination of the effects of heterogeneous organization of RyR clusters, myofibrils and mitochondria on Ca2+ release patterns in cardiomyocytes. PLoS Comp. Biol. 11 (9), e1004417 (2015).
  19. Illian, J., Penttinen, A., Stoyan, H., Stoyan, D. Statistical Analysis and Modelling of Spatial Point Patterns. Statistical Analysis and Modelling of Spatial Point Patterns. , John Wiley & Sons. Chichester, UK. 1-534 (2008).
  20. Kremer, J. R., Mastronarde, D. N., McIntosh, J. R. Computer Visualization of Three-Dimensional Image Data Using IMOD. J. Struct. Biol. 116, 71-76 (1996).
  21. Fang, Q., Boas, D. A. Tetrahedral mesh generation from volumetric binary and grayscale images. Proc. ISBI. , 1142-1145 (2009).
  22. Aune, D. J., Herr, S. E., Menick, D. R. Induction and assessment of ischemia-reperfusion injury in langendorff perfused rat hearts. J. Vis. Exp. (101), e52908 (2015).
  23. Judd, J., Lovas, J., Huang, G. N. Isolation, culture and transduction of adult mouse cardiomyocytes. J. Vis. Exp. (114), (2016).
  24. Hagler, H. K. Ultramicrotomy for biological electron microscopy. Electron Microscopy: Methods and Protocols. 369 (Chapter 5), 67-96 (2007).
  25. He, W., He, Y. Electron tomography for organelles, cells, and tissues. Electron Microscopy: Methods and Protocols. 1117 (20), 445-483 (2014).
  26. Diggle, P., Marron, J. S. Equivalence of smoothing parameter selectors in density and intensity estimation. J. Am. Stat. Assoc. 83 (403), 793-800 (1988).
  27. Ghosh, S., Crampin, E. J., Hanssen, E., Rajagopal, V. A computational study of the role of mitochondrial organization on cardiac bioenergetics. Proc. EMBC. , 2696-2699 (2017).
  28. Pinali, C., Kitmitto, A. Serial block face scanning electron microscopy for the study of cardiac muscle ultrastructure at nanoscale resolutions. J. Mol. Cell. Cardiol. 76, 1-11 (2014).
  29. Hussain, A., Hanssen, E., Rajagopal, V. A Semi-Automated Workflow for Segmenting Contents of Single Cardiac Cells from Serial-Block-Face Scanning Electron Microscopy Data. Microsc Microanal. 23 (S1), 240-241 (2017).
  30. Pinali, C., Bennett, H., Davenport, J. B., Trafford, A. W., Kitmitto, A. Three-dimensional reconstruction of cardiac sarcoplasmic reticulum reveals a continuous network linking transverse-tubules: this organization is perturbed in heart failure. Circ. Res. 113 (11), 1219-1230 (2013).
  31. LeGrice, I. J., Hunter, P. J., Smaill, B. H. Laminar structure of the heart: a mathematical model. Am. J. Physiol. 272 (5 Pt 2), H2466-H2476 (1997).
  32. Jayasinghe, I. D., Cannell, M. B., Soeller, C. Organization of ryanodine receptors, transverse tubules, and sodium-calcium exchanger in rat myocytes. Biophys. J. 97 (10), 2664-2673 (2009).
  33. Jayasinghe, I. D., Crossman, D. J., Soeller, C., Cannell, M. B. Comparison of the organization of t-tubules, sarcoplasmic reticulum and ryanodine receptors in rat and human ventricular myocardium. Clinic. Exp. Pharmacol. P. 39 (5), 469-476 (2012).
  34. Bradley, C., Bowery, A., et al. OpenCMISS: a multi-physics & multi-scale computational infrastructure for the VPH/Physiome project. Prog. Biophys. Mol. Bio. 107 (1), 32-47 (2011).

Tags

バイオ エンジニア リング、問題 134、心筋細胞、心筋の微細構造、電子顕微鏡、共焦点顕微鏡、計算力学モデル、有限要素モデルの作成、システム生物学、カルシウム シグナル、ミトコンドリア、生体エネルギー
心筋細胞のシステム生物学における細胞構造の役割を研究する心筋細胞の構造的に現実的な有限要素の幾何学的なモデルを作成します。
Play Video
PDF DOI DOWNLOAD MATERIALS LIST

Cite this Article

Rajagopal, V., Bass, G., Ghosh, S.,More

Rajagopal, V., Bass, G., Ghosh, S., Hunt, H., Walker, C., Hanssen, E., Crampin, E., Soeller, C. Creating a Structurally Realistic Finite Element Geometric Model of a Cardiomyocyte to Study the Role of Cellular Architecture in Cardiomyocyte Systems Biology. J. Vis. Exp. (134), e56817, doi:10.3791/56817 (2018).

Less
Copy Citation Download Citation Reprints and Permissions
View Video

Get cutting-edge science videos from JoVE sent straight to your inbox every month.

Waiting X
Simple Hit Counter