Summary
このプロトコルは、心臓のマルチスケール電気機械モデルを適用することにより、患者固有の超音波画像から左心室のパラメトリックモデルを自動的に生成するためのSILICOFCMプラットフォームのワークフローを示しています。このプラットフォームは、実際の臨床試験を減らし、肯定的な治療結果を最大化することを目的とした インシリコ 臨床試験を可能にします。
Abstract
SILICOFCMプロジェクトは、主に家族性心筋症(FCM)の インシリコ 臨床試験のための計算プラットフォームを開発することを目的としています。このプラットフォームのユニークな特徴は、患者固有の生物学的、遺伝的、および臨床画像データの統合です。このプラットフォームにより、治療の検査と最適化が可能になり、肯定的な治療結果を最大化できます。したがって、副作用および薬物相互作用を回避し、心臓突然死を予防することができ、薬物治療の開始から所望の結果までの時間を短縮することができる。本稿では、心臓の電気機械モデルを適用することにより、患者固有の超音波画像から自動的に生成された左心室のパラメトリックモデルを紹介します。薬物効果は、入口と出口の流れ、ECG測定、および心筋特性のカルシウム機能の特定の境界条件を通じて処方されました。患者からの遺伝子データは、心室壁の物質的特性を通じて組み込まれました。頂端ビュー分析では、事前にトレーニングされたU-netフレームワークを使用して左心室をセグメント化し、拡張期および収縮期サイクルの左心室の長さに基づいて境界の長方形を計算します。Mモードビュー分析には、Mモードビューの左心室の特徴的な領域の境界が含まれます。左心室の寸法を抽出した後、メッシュオプションに基づいて有限要素メッシュを生成し、ユーザーが指定した入口速度と出口速度で有限要素解析シミュレーションを実行しました。ユーザーは、圧力-体積、圧力-ひずみ、心筋作業時間図などのさまざまなシミュレーション結果や、変位、圧力、速度、せん断応力などのさまざまな分野のアニメーションをプラットフォーム上で直接視覚化できます。
Introduction
近年、情報技術、シミュレーションソフトウェアパッケージ、医療機器の急速な発展により、大量の臨床情報を収集する機会が提供されています。したがって、包括的で詳細な計算ツールを作成することは、利用可能な豊富なデータから特定の情報を処理するために不可欠になっています。
医師の観点からは、疾患の進行、治療反応、および将来のリスクを推定するために、特定の患者の「正常」表現型と「異常」表現型を区別することが最も重要です。最近の計算モデルは、肥大型心筋症(HCM)および拡張型(DCM)心筋症における心筋の挙動の統合的理解を大幅に改善しました1。心臓全体の電気的活動の高解像度、詳細、解剖学的に正確なモデルを使用することが重要であり、これには膨大な計算時間、専用ソフトウェア、およびスーパーコンピューターが必要です1,2,3。ホルツァプフェル実験に基づく線形弾性および直交異方性材料モデルを使用して、実際の3D心臓モデルの方法論が最近開発され、心臓4内の電気信号輸送および変位場を正確に予測できます。新しい統合モデリングアプローチの開発は、多遺伝子障害患者の症状の種類と重症度を区別し、正常な身体活動の障害の程度を評価するための効果的なツールになる可能性があります。
ただし、患者固有のモデリングには多くの新しい課題があります。人間の心臓の物理的および生物学的特性を完全に決定することは不可能です。非侵襲的な測定には通常、個々の患者の特定のパラメータを推定するのが難しいノイズの多いデータが含まれます。大規模な計算は実行に多くの時間を必要としますが、臨床的な時間枠は限られています。患者の個人データは、生成されたメタデータを患者の機密性を損なうことなく再利用できるように管理する必要があります。これらの課題にもかかわらず、マルチスケール心臓モデルには、観察された過渡反応に厳密に従う予測を達成するのに十分な詳細レベルを含めることができ、それによって将来の臨床応用の可能性を提供します。
しかし、複数の研究室による多大な科学的努力と多額の助成金支援にもかかわらず、現在、SIMULIA Living Heart Model5と呼ばれるマルチスケールおよび全心臓シミュレーション用の市販のソフトウェアパッケージは1つだけです。これには、動的電気機械シミュレーション、洗練された心臓形状、血流モデル、および受動的および能動的特性、線維性、および電気経路を含む完全な心臓組織の特性評価が含まれます。このモデルは個別化医療での使用を目標としていますが、活物質の特性評価はGuccioneらによって導入された現象論的モデルに基づいています6,7。したがって、SIMULIAは、多くの心臓病で観察される収縮性タンパク質の機能特性の変化を直接かつ正確に翻訳することはできません。これらの変化は、分子レベルおよび細胞内レベルでの突然変異やその他の異常によって引き起こされます6。臨床現場での少数のアプリケーションに対するSIMULIAソフトウェアの使用が限られていることは、より高いレベルのマルチスケールヒト心臓モデルの開発における今日の苦労の好例です。一方、分子スケールから臓器スケールへの突然変異の影響を追跡できる新世代のマルチスケールプログラムパッケージの開発を動機付けています。
心臓の電気生理学の主な目的は、胴体内の信号伝播とすべてのコンパートメントの特性を決定することです4,5,6。SILICOFCM8プロジェクトは、患者固有の生物学的、遺伝的、および臨床画像データを使用して心筋症の発症を予測します。これは、現実的なサルコメア系、患者の遺伝的プロファイル、筋繊維の方向、流体-構造相互作用、および電気生理学結合のマルチスケールモデリングによって達成されます。左心室変形、僧帽弁の動き、および複雑な血行動態の影響は、特定の患者の心臓状態の詳細な機能的挙動を与えます。
この記事では、電気機械結合を備えた流体構造心臓モデルを使用して、患者固有の超音波画像から自動的に生成された左心室(LV)のパラメトリックモデルにSILICOFCMプラットフォームを使用する方法を示します。LVの頂端ビューとMモードビュー解析は、深層学習アルゴリズムで生成されました。次に、メッシュジェネレータを使用して有限要素モデルを自動的に構築し、LV収縮9の全サイクルのさまざまな境界条件をシミュレートしました。このプラットフォームでは、圧力-体積図、圧力-ひずみ図、心筋作業時間図などのシミュレーション結果や、変位、圧力、速度、せん断応力などのさまざまな分野のアニメーションを直接視覚化できます。特定の患者からの入力パラメータは、超音波画像からの幾何学的形状、LVの入力および出力境界流条件における速度プロファイル、および特定の薬物療法(例えば、エンレスト、ジゴキシン、マバカムテンなど)である。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
この研究のプロトコルは、2019年2月6日に英国国民保健サービス健康研究局ノースイーストタイン&ウェアサウス研究倫理委員会によって参照番号18/NE/0318で承認され、各参加センターの治験審査委員会によって採択されました。この研究は、グッドクリニカルプラクティスの原則の範囲内で、ヘルシンキ宣言に従って実施されました。インフォームドコンセントは、研究に関与したすべての被験者から得られた。患者情報は匿名に保たれます。
1.超音波MモードまたはアピカルビューDICOM画像分析とLVパラメータ抽出のワークフロー
注:このプロトコルを開始するには、ユーザーはSILICOFCMプラットフォーム8 にログインし、適切なワークフロー(つまり、Mモードまたはアピカルビューを使用した超音波画像分析)を選択する必要があります。たとえば、超音波MモードDICOM画像解析とLVパラメータ抽出のワークフローには、いくつかのステップが含まれます。最初のステップはテンプレートマッチングです。照合するテンプレートには、関連するすべての必要な境界線が含まれます。新しいデータセットごとに、特定の超音波装置に対応するデータセットごとにテンプレートを手動で抽出する必要があることを強調しておく必要があります。ワークフロー図を 図 1 に示します。テンプレートに一致する領域は、Felsenszwalbの効率的なグラフベースの画像セグメンテーションアルゴリズムがセプタムの境界とLVの壁に対応する「強い」境界を見つけた後、分析された画像から抽出されます。 これらの境界と、心臓の直径が最大(拡張期の直径に対応)および最小(収縮期の直径に対応)である場所に基づいて、 さまざまなLV寸法が計算されます。ユーザーは、ビューが M モード ビューかアピカル ビューかを定義する必要があります。
- ユーザー名とパスワードを使用してプラットフォームにログインします。 仮想母集団モジュールで、超音波Mモードまたは頂端ビューワークフローを選択します。
- 使用可能なワークフローのリストで、 超音波マージ ワークフローを選択します。
- [ ファイルのアップロード] セクションで、ユーザーのコンピューターにローカルに保存されている画像とDICOMファイルをアップロードします(図2)。
注意: 特定の患者からのモジュールMモードビューDICOM画像分析とLVパラメーター抽出への入力は、テスト画像(Mモードまたはアピカルビューを表すDICOM形式でユーザーから提供されます)と静的ファイル(図1)です。静的ファイルは、a)モジュールから計算された出力メジャーを上書きするために使用されるテンプレートファイルである 入力.txt ファイル(新しい出力ファイルはこのファイルから順番に作成されます)、およびb)テンプレートマッチングに使用されるテンプレート画像ファイル(つまり、 template_dicom_GEMS.jpg)で構成されます。 - ファイルの 保存先 フォルダとして[プライベートフォルダ]または[パブリックフォルダ]を選択します。ワークフローを開始する前に、[ コメント ] セクションに必要なコメントまたはメモを入力します。以前にアップロードしたアピカルビューとMモードLV超音波画像とDICOMファイルを選択します(ステップ1.3)。(図3)。
注意: 超音波心臓セグメンテーションツールは、頂端ビューとMモードビューの2つのサブモジュールで構成されています。両方のモジュールは、プラットフォームを介してアクセスされます。分析には、対応する断面に表示される特性パラメータの計算が含まれます。どちらの場合も、必要な値は、ユーザー側からパラメトリック心臓モデルへの入力として取得されるファイルに書き込まれます(詳細については、説明を参照してください)。 - [ 実行] ボタンをクリックします。対話型プラットフォームは、実行中のワークフローが終了するとユーザーに通知します。
- 作成されたLVのジオメトリをプラットフォーム上で直接視覚化します(3Dモデルが画面に自動的に表示され、マウスを使用して回転できます)。使用可能なオプションには、視覚化用の シェーディング モデルと ワイヤフレームモデル があります。
注: LV 寸法が抽出された後、メッシュ オプションに基づいて有限要素のメッシュが生成され、ユーザーが指定した入口速度と出口速度で有限要素解析シミュレーションが実行されます(以下を参照)。
2. 流体-構造(FS)シミュレーション用のPAK有限要素ソルバー
注: このツールは、固体流体連成問題の有限要素解析に使用できます。固体と流体の間の強結合と緩連成の両方をサポートします。要素は六面体または四面体で、要素の中心に追加の節点がある場合とない場合があります。このソルバーには、ホルツァプフェルモデル、ハンターマッスルモデルなどの材料モデルが組み込まれています。PAK-FS の情報フロー図を 図 4 に示します。入力ファイルと PAK プリプロセッサから開始されます。PAK プリプロセッサ ツールは、有限要素ソルバーの入力ファイルとして機能する DAT ファイルを出力します。ソルバーからの最終出力は、有限要素シミュレーションの結果(左心室の速度、圧力、ひずみ、応力)を含むVTKファイルです。
- 下部セクションのボタンを使用して、入口と出口の速度の境界条件のテンプレートファイルをダウンロードします(図3)。患者固有の流れ境界条件が利用可能な場合は、これらのファイルをダウンロードして使用してください。下部の対応するボタンをクリックして、メッシュオプションをダウンロードします。これらのファイルをプライベート フォルダまたはパブリック フォルダに保存します。
- これらのファイルは、画像のアップロードと同様の方法でアップロードします (図 3)。規定の入口速度と出口速度は薬物の状態をシミュレートし、メッシュオプションは有限要素メッシュの密度を制御します。患者固有の状態をシミュレートするには、圧力、流量、材料特性、カルシウム機能のデフォルト値を変更します。
- [ 実行] ボタンをクリックします。実行中の新しいワークフローがリストに表示されます。ワークフローのいずれかのセクションが明確でない場合は、[ ヘルプファイル ]ボタン(図3)をクリックして、このワークフローの使用方法の詳細な手順を表示し、結果を解釈します。
注: すべてがエラーなしで実行されると、ワークフローの状態が "実行中" から "完了 OK" に変わります。 - いくつかのオプションで結果を表示します。または、結果をダウンロードします。結果フォルダーには、VTK ファイル、CSV ファイル、およびアニメーションが含まれています。
- 目のボタンをクリックすると、駆出率と全体的な作業効率の値、および圧力対体積、圧力対ひずみ、心筋作業対時間の図が表示されます。
- カメラボタンをクリックして、変位、圧力、せん断応力、速度場のアニメーションをプレビューして再生します。
- 3Dビジュアライゼーションボタンをクリックして、 ParaView Glanceで出力をオンラインで視覚化します。
- 結果として以前にダウンロードした複数のVTKファイルをロードします。関心のあるいくつかのパラメーターを確認し、たとえば、視覚化のためにフィールドを速度に変更します。
- モデルを回転するか、配色を変更します。サーフェスの表現には「 エッジ付きサーフェス 」(Surface) または 「ワイヤフレーム」(Wireframe ) を選択します。ロードされたすべてのVTKファイルに同じ方法を適用します。
3. 心電図計測による心室活性化シーケンスの決定
注:心臓細胞の修正フィッツヒュー-南雲モデルが実装されました。前胸部リードは、標準の6つの電極でモデル化されました。心臓の可能性は逆ECGで最適化されました。洞房結節(時間の関数)の活性化から始まり、心臓と胴体を通る不均一な活動電位で、ユーザーは胴体全体のモデルで電気的活動を得ることができます。PAK-TORSOの情報フロー図を 図5に示します。ユーザーは、胴体モデルの全方向(x、y、z)とECG信号機能のスケーリングを提供します。次に、スケーリングされたモデルが作成され、PAK-FSソルバーを使用してその動作がシミュレートされます。ユーザーは、これらの入力値をテキスト ファイルで指定します。シミュレーションの出力は、胴体埋め込み環境での心臓の電気的活動を含むVTKファイルです。
- ホームページで、[ ワークフローの実行] に移動し、使用可能なワークフローの一覧で torse-cwl を選択します。ワークフローを実行する前に、[ コメント ] セクションにコメントまたはメモを追加します。
- [ 入力テンプレートファイル ]ボタンをクリックし、Webページに表示されているコンテンツを、胴体モデルに使用される 入力.txt ファイルとして保存します。[ 入力ファイル ] フィールドで、ダウンロードした 入力.txt ファイルを選択します。ファイルがインポートされたら、[ 実行 ]ボタンをクリックして計算を開始します。
注:新しいワークフローが左下に表示され、ステータスが「実行中」と表示されます。計算は約0.5時間続き、その後ステータスが「完了OK」に変わります。 - 左下隅にある目のボタンまたはカメラのボタンをクリックして、利用可能なシミュレーションレポートまたはアニメーションをプラットフォーム上で直接視覚化します。
注: 結果には、各時間ステップにおける電界、速度、圧力、せん断応力、および変形の分布が含まれます。これらの各分布のアニメーションは、心周期全体にわたって表示することができます。 - または、3Dビジュアライゼーションボタンをクリックして、 ParaView Glanceで出力をオンラインで視覚化します。
- [ファイルを開く] ボタンを選択し、[GIRDER] タブに移動し、プロンプトが表示されたらユーザーの資格情報を入力して、[プライベート] フォルダーを開きます。
- 次のページで、 ワークフロー出力 フォルダーを選択し、 torso-cwl フォルダーを開きます。リストの最初のフォルダを開きます。
- シミュレーション結果を表すVTKファイルのリストを参照してください。1つ以上のファイルを選択し、[ 選択 ]ボタンをクリックして、 ファイルをParaViewGlanceにロードします。
- マウスを使用してモデルジオメトリを操作(つまり、移動、回転、ズームインまたはズームアウトなど)します。
- 次のように、モデル表現のさまざまなオプションを選択します。
- [ ワイヤフレーム] オプションを選択すると、胴体内にハートが組み込まれた胴体の内部が表示されます。[ ポイント ]オプションを選択すると、完全なハートメッシュを持つ胴体モデルの点線表現が表示されます。
- [ ポイント サイズ ] の値を調整して、表示結果を変更します。 [不透明度](Opacity )の値を調整して、胴体の内部を確認し、ハートメッシュ内に結果を表示します。
- [ 色] ドロップダウンメニューをクリックして、電位などの目的のオプションを選択します。デフォルトのカラースケールをリストされたオプションのいずれかに変更します。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
例として、超音波MモードDICOM画像解析とLVパラメータ抽出のワークフローを 図1に示します。Mモードと頂端ビューは、関心のあるパラメータに応じて、別々にテストすることも、次々にテストすることもできます。次々にテストした場合、結果は1つの共通ファイルに追加されます(収縮期と拡張期では別々に)。1つのビューのみがテストされる場合、不明なパラメータの値はデフォルトのファイル入力から取得されます .txt (図1)。SILICOFCMツールは、表示モードに応じて次のように画像を解析します。
アピカルビュー分析には、以前にトレーニングされたU-netフレームワークを使用してLVをセグメント化し、それに基づいてパラメータLV長、拡張期、2D(LVLd[cm])およびLV長さ、収縮期、2D(LVLs[cm])を計算する境界長方形を計算することが含まれます。ユーザーは、ビューが収縮期または拡張期のどちらを表すかを定義する必要があります。出力ファイルは次のとおりです:a) input_parametric_diastole.stl または input_parametric_systole.stlは、画像が拡張期または収縮期を表すかどうかのユーザー入力に応じてLVLd[cm]またはLVLs[cm]の計算値を含む、作成されたジオメトリを持つPAKソルバーの入力ファイルを表します。b) 画像が拡張期または収縮期を表すかどうかのユーザー入力に応じて、LVLd[cm]またはLVLs[cm]のデータを含むoutput_Apical_view.txt。
Mモードビュー分析には、このビューのLVの特性領域の境界が含まれます。これらの領域に基づいて、さまざまなパラメータが計算されます。出力ファイルは次のとおりです:a) input_parametric_diastole.stl および input_parametric_systole.stlは、心室中隔の厚さ、拡張期、Mモード(IVSd[cm])、心室中隔の厚さ、収縮期、Mモード(IVSs[cm])、LV内部寸法、拡張期、Mモード(LVIDd[cm])、LV内部寸法、収縮期、Mモード(LVIDs[cm])、LV後壁の厚さの計算値を含む、PAKソルバーの入力ジオメトリを表します。 拡張期、Mモード(LVPWd[cm])、およびLV後壁の厚さ、収縮期、Mモード(LVPW [cm]);b) IVSd [cm]、 IVS [cm]、LVIDd [cm]、LVID [cm]、LVPWd [cm]、およびLVPW [cm]のデータを含むoutput_M_mode.txt。
LVの例を生成するために、可能な限り現実的に、特定の部品、すなわちベースパーツ、バルブ(大動脈と僧帽弁)、および接続パーツ(つまり、ベースとバルブ間の接続)を含むパラメトリックタイプのモデルが、LVジオメトリの計算されたパラメータを使用してPAK有限要素シミュレーション(ステップ2)を介して生成されました。すべてのパーツには、長さと半径、およびレイヤーの数(分割)があります。有限要素モデルは、ベースと接続層の一部に接続されたソリッドウォールに囲まれた流体領域で構成されます。LVモデルの現実的な挙動をシミュレートするために、入口(僧帽; 図6A)および出口(大動脈; 図6B)バルブを使用した。繊維方向の自動計算アルゴリズムをこの有限要素モデルに適用しました。1層と3層のソリッド表現の結果を 図7に示します。
図8 は、1秒の時間サイクルにおけるパラメトリックモデル内部の圧力分布を10の時間ステップに分割したものです。最初の5つのステップでは、収縮がスケールの最高値に達し、モデルの体積が同時に増加し始めるまで、大きな圧力変化はありません。 図9 は、LVモデルの流体部分内で速度がどのように分布しているかを示しています。示されているように、分岐内には、ロード/アンロードサイクル中の流体の流れによって引き起こされる顕著な値のピークがあります。 図 10 は、変位がモデルに沿ってどのように分布しているかを示しています。圧力変化と同様に、最初の2つのステップでは、モデルの下部で最大になる収縮まで変位は無視できます。残りの時間を通して、モデルはゆっくりと変形していない状態に戻ります。図 11 は、結果として得られる圧力-体積(PV)図を示しています。 図12 は、心室モデルの流体部分の内部の速度のベクトル表現を示しています。位相速度ベクトルは、流体の流入中に僧帽弁の内部に存在し、心室が0.7秒からの流体で満たされた後、流体は大動脈部分を通って出て行きます。
電気生理学をシミュレートするために、完全結合心臓胴体システムのモノドメイン方程式を適用しました。心筋および胴体の伝導率は、表1に定義される10、11、12、13、14である。ゼロフラックスは、肺、胴体、および心臓と接触するすべての内部境界の境界条件としてVmに使用されました。したがって、−n ·Γ = 0、ここでnは境界上の単位外向き法線ベクトル、Γは細胞内電圧に対するその境界を通る磁束ベクトルであり、Γ = −σ·∂Vm/∂n.次に、心外膜電位定式化を用いてECG逆問題を解くための古典的なアプローチを実施した。使用した方法は、Tikhonov法およびL正則化ベースの方法のファミリー10,11,12,13,14であった。
胴体埋め込み環境における心臓全体の電気的活動を 図13に示します。心臓全体に不均一な活動電位形態を有する特殊な伝導系が、洞房結節の活性化の開始に組み込まれた。 図14 は、測定されたECG信号に対応する9つのシーケンスにおける心室活性化の進行中の健常対象における体表面電位のマップを示す。
図1:超音波MモードDICOM画像解析とLVパラメータ抽出の情報フロー図。 この図の拡大版を表示するには、ここをクリックしてください。
図2:分析用の新しいファイルのアップロード。 この図の拡大版を表示するには、ここをクリックしてください。
図3:超音波アピカルビューDICOM画像分析とLVパラメータ抽出のユーザーインターフェイス。 この図の拡大版を表示するには、ここをクリックしてください。
図4:PAK-FSの情報フロー図。 この図の拡大版を表示するには、ここをクリックしてください。
図5:PAK-TORSO情報フロー図。 この図の拡大版を表示するには、ここをクリックしてください。
図6:速度プロファイル 。 (A)僧帽弁断面における流速の入口関数、(B)大動脈弁断面における出口速度関数。 この図の拡大版を表示するには、ここをクリックしてください。
図7:繊維方向を持つパラメトリックLVモデル 。 (A)単層ソリッド壁、および(B)3層ソリッド壁表現。 この図の拡大版を表示するには、ここをクリックしてください。
図8:パラメトリックLVモデル内の圧力場。 5つの異なる時間ステップが提示されます。 この図の拡大版を表示するには、ここをクリックしてください。
図9:パラメトリックLVモデル内の速度場。 5つの異なる時間ステップが提示されます。 この図の拡大版を表示するには、ここをクリックしてください。
図10:パラメトリックLVモデルの変位。 5つの異なる時間ステップが提示されます。 この図の拡大版を表示するには、ここをクリックしてください。
図11:LV流体-構造相互作用モデルの圧力対体積(PV)図1秒の規定サイクル期間は、10の時間ステップに分割されます。5 つの代表的な時間ステップが示されています。この図の拡大版を表示するには、ここをクリックしてください。
図12:左心室のパラメトリックモデルにおける速度のベクトル表現。 4 つの代表的な時間ステップが示されています。 この図の拡大版を表示するには、ここをクリックしてください。
図13:リードIIのECG信号のさまざまな時点での心臓全体の活性化のシミュレーション 。 (B)の1〜9の活性化シーケンスは、(A)のECG信号に対応します。mV単位の膜貫通電位はカラーバーで示されます。 この図の拡大版を表示するには、ここをクリックしてください。
図14:健康な被験者の体表面電位のマップ。 心室活性化の進行は、ECG信号に対応する9つのシーケンス(下パネル)(上パネル)。mV単位の心臓活動範囲はカラーバーで示されます。 この図の拡大版を表示するには、ここをクリックしてください。
パラメーター | サン | 心房 | ティッカー | 彼の | ティッカー | プルキンエ | 心室 | |
ある | -0.6 | 0.13 | 0.13 | 0.13 | 0.13 | 0.13 | 0.13 | |
B | -0.3 | 0 | 0 | 0 | 0 | 0 | 0 | |
c1(AsV−1 m−3) | 1000 | 2.6 | 2.6 | 2.6 | 2.6 | 2.6 | 2.6 | |
c2 (AsV−1 m−3) | 1 | 1 | 1 | 1 | 1 | 1 | 1 | |
D | 0 | 1 | 1 | 1 | 1 | 1 | 1 | |
e | 0.066 | 0.0132 | 0.0132 | 0.005 | 0.0022 | 0.0047 | 0.006 | |
A (mV) | 33 | 140 | 140 | 140 | 140 | 140 | 140 | |
B (mV) | -22 | -85 | -85 | -85 | -85 | -85 | -85 | |
k | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | 1000 | |
σ (mS·m-1) | 0.5 | 8 | 0.5 | 10 | 15 | 35 | 8 |
表1:修正フィッツヒュー・ナフモ方程式を用いたモノドメインモデルのパラメータ。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
SILICOFCMプロジェクトは、リスク予測のための仮想患者集団を設計し、薬理学的治療の効果をテストし、動物実験とヒト臨床試験を減らすための インシリコ 臨床試験プラットフォームです。薬理学的処理の効果をテストすることは、規定の入口/出口境界流条件、カルシウム機能、および材料壁特性でモデル化されました。このプラットフォームは、心筋症の進行中の患者固有の状態のリスク予測のための、心臓全体のパフォーマンスおよび機能最適化レベルを備えたサルコメアレベルのマルチスケールメソッドを統合します。
SILICOFCM8 プロジェクトにおける心筋症とLVの電気機械的結合の心臓モデリングが発表されました。心臓モデルの幾何学的形状は7つの異なる領域を含んでいた:1)洞房結節;2)心房;3)房室結節;4)彼の束。5)束繊維;6)プルキンエ繊維;7)心室心筋。心室活性化の進行中の健常者における体表面電位マップは、ECG信号に対応する9つの配列で提示されている。
メッシュオプションに基づいて左心室の寸法が抽出された後、有限要素のメッシュが生成され、ユーザーが指定した入口速度と出口速度で有限要素解析シミュレーションが実行されます。ユーザーは、利用可能なアニメーションや図を見ることで、プラットフォーム上で直接ソリューションを視覚化できます。ユーザーは、目の記号をクリックすることで、圧力-体積、圧力-ひずみ、および心筋の作業時間図を視覚化できます。ユーザーがカメラボタンをクリックすると、さまざまなフィールド(変位、圧力、速度、せん断応力)の使用可能なアニメーションのリストが表示されます。
研究にはいくつかの制限があります。パラメトリックLV形状は超音波画像から抽出した。将来のバージョンでは、より詳細な幾何学的再構成に入ります。心筋症の特定の薬との直接的な関係は、この原稿には示されていません。本研究では、流量と圧力の境界条件によって制御します。現在、カルシウム心筋機能と壁の非線形材料特性を通じて遺伝子データが組み込まれています。将来のバージョンでは、患者固有の心筋症(肥大型および拡張した遺伝的変異)の画像再構成と遺伝データの詳細が考慮されます。
SILICOFCM計算プラットフォーム8 は、特に心臓病と患者固有の状態のリスク予測のための インシリコ 臨床試験の新しい道を開きます。リスク予測のための今日の臨床診療のゴールドスタンダードは、心筋症患者の標準的な生存計算機です。ここで使用されるプラットフォームは、モデリングにはバイオマーカーだけでなく、患者固有の形状、流れと圧力の血行動態条件、壁の材料特性(画像の変形から)、および境界条件、カルシウム機能、および材料特性のさまざまな組み合わせによる薬物反応も含まれるため、現在の医療標準と比較してより多くの情報を提供できます。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
著者には利益相反はありません。
Acknowledgments
この研究は、契約番号451-03-68 / 2022-14 / 200107を通じて、助成金契約SILICOFCM 777204に基づく欧州連合のHorizon 2020研究およびイノベーションプログラムおよびセルビア共和国教育科学技術開発省によってサポートされています。この記事は著者の見解のみを反映しています。欧州委員会は、記事に含まれる情報の使用について責任を負いません。
Materials
Name | Company | Catalog Number | Comments |
SILICOFCM project | www.silicofcm.eu | open access for registered users |
References
- Gibbons Kroeker, C. A., Adeeb, S., Tyberg, J. V., Shrive, N. G. A 2D FE model of the heart demonstrates the role of the pericardium in ventricular deformation. American Journal of Physiology. 291 (5), 2229-2236 (2006).
- Pullan, A. J., Buist, M. L., Cheng, L. K. Mathematically Modelling the Electrical Activity of the Heart - From Cell To Body Surface and Back Again. , World Scientific. (2005).
- Trudel, M. -C., Dub´e, B., Potse, M., Gulrajani, R. M., Leon, L. J. Simulation of QRST integral maps with a membrane based computer heart model employing parallel processing. IEEE Transactions on Biomedical Engineering. 51 (8), 1319-1329 (2004).
- Kojic, M., et al. Smeared multiscale finite element models for mass transport and electrophysiology coupled to muscle mechanics. Frontiers in Bioengineering and Biotechnology. 7, 381 (2019).
- Baillargeon, B., Rebelo, N., Fox, D. D., Taylor, R. L., Kuhl, E. The Living Heart Project: A robust and integrative simulator for human heart function. European Journal of Mechanics - A/Solids. 48, 38-47 (2014).
- Guccione, J. M., McCulloch, A. D. Mechanics of active contraction in cardiac muscle: Part I--Constitutive relations for fiber stress that describe deactivation. TheJournal of Biomechanical Engineering. 115, 72-81 (1993).
- Guccione, J. M., Waldman, L. K., McCulloch, A. D. Mechanics of active contraction in cardiac muscle: Part II--Cylindrical models of the systolic left ventricle. The Journal of Biomechanical Engineering. 115, 82-90 (1993).
- H2020 project SILICOFCM: In Silico. trials for drug tracing the effects of sarcomeric protein mutations leading to familial cardiomyopathy. , Available from: www.silicofcm.eu 2018-2022 (2022).
- Kojic, M., et al. Smeared multiscale finite element model for electrophysiology and ionic transport in biological tissue. Computers in Biology and Medicine. 108, 288-304 (2019).
- Wang, Y., Rudy, Y. Application of the method of fundamental solutions to potential-based inverse electrocardiography. Annals of Biomedical Engineering. 34 (8), 1272-1288 (2006).
- Van Oosterom, A. The use of the spatial covariance in computing pericardial potentials. IEEE Transactions on Biomedical Engineering. 46 (7), 778-787 (1999).
- Van Oosterom, A. The spatial covariance used in computing the pericardial potential distribution. Computational Inverse Problems in Electrocardiography. , 1-50 (2001).
- Van Oosterom, A.
Source models in inverse electrocardiography. International Journal of Bioelectromagnetism. 5, 211-214 (2003). - Van Oosterom, A. The equivalent double layer: source models for repolarization. Comprehensive Electrocardiology. , Springer. 227-246 (2010).