Summary
この研究では、一括パラメータアプローチと有限要素分析に基づいて、保存された駆出率を持つ心不全の2つの計算モデルを紹介する。これらのモデルは、圧力過負荷および心室コンプライアンスの低下によって誘発される左心室および関連血管系のヘモダイナミクスの変化を評価するために使用される。
Abstract
心血管疾患の計算モデリング分野における科学的取り組みは、主に駆出率(HFrEF)の減少した心不全に焦点を当てており、最近では世界的に主要な心不全となっている保存された駆出率(HFpEF)を持つ心不全を広く見落としています。この論文では、 インシリコ表現における HFpEFの貧弱性によって動機づけられ、左心室圧過多に起因するHFpEFのヘモダイナミクスをシミュレートする2つの異なる計算モデルが提示される。まず、数値ソルバーを使用してオブジェクト指向の一括パラメータモデルを開発しました。このモデルは、構成要素の幾何学的および機械的特性に依存し、低い計算コストの利点を提供するゼロ次元(0D)Windkesselのようなネットワークに基づいています。第2に、有限要素解析(FEA)ソフトウェアパッケージを多次元シミュレーションの実施に利用した。FEA モデルは、電気機械心応答、構造変形、流体キャビティベースのヘモダイナミクスの 3 次元(3D)マルチフィジックス モデルを組み合わせて、簡略化された一括パラメータ モデルを利用して、異なる流体キャビティ間の流動交換プロファイルを定義します。各アプローチを通じて、圧力過負荷に起因する左心室および近位血管系の急性および慢性血力学的変化の両方をうまくシミュレートした。具体的には、圧力過負荷を大動脈弁のオリフィス面積を減少させることでモデル化し、一方、慢性リモデリングは左心室壁のコンプライアンスを低下させることによりシミュレートした。HFpEFの科学的および臨床的文献と一致して、(i)左心室と大動脈の間の経大動脈圧勾配の急性上昇および脳卒中容積の減少および(ii)拡張期の機能不全を示す末方拡張期左心室容積の慢性的な減少を示す。最後に、FEAモデルは、HFpEF心筋のストレスが心臓周期を通じて健康な心臓組織よりも著しく高いことを示している。
Introduction
心不全は、心臓が身体の代謝要求に追いつくために十分にポンプまたは充填することができないときに起こる世界的な主要な死因である。この放出画分、すなわち各収縮で排出される左心室に記憶される血液の相対量は、(i)心不全を減少した駆出率(HFrEF)および(ii)心不全に分類するために臨床的に使用され、保存された駆出画分(HFpEF)を有する心不全、噴出分画はそれぞれ45%未満、または45%未満、2である。HFpEFの症状は、しばしば左室圧過負荷に反応して発症し、これは大動脈狭窄、高血圧、および左心室流出管閉塞3、4、5、6、7を含むいくつかの状態によって引き起こされる可能性がある。圧力過負荷は、分子および細胞収差のカスケードを駆動し、左心室壁の肥厚(同心リモデリング)を導き、最終的には、壁硬化またはコンプライアンス8、9、10の損失につながる。これらの生体力学的変化は、心血管の断熱圧と気圧の関係が高まり、エンド・拡張期体積11の減少をもたらすため、心血管のヘモダイナミクスに大きな影響を及ぼす。
心血管系の計算モデリングは、生理学と疾患の両方における血圧および流れの理解を進め、診断および治療戦略の開発を促進した12。シリコモデルでは、低次元または高次元モデルに分類され、前者は低い計算需要を有する世界的な血行力学特性を評価する分析法を利用し、後者は2Dまたは3Dドメイン13における心血管力学および血行力学のより広範な多尺度および多物理記述を提供する。一括パラメータの Windkessel 表現は、低次元の記述の中で最も一般的です。電気回路のたとえ(オームの法則)に基づいて、これは抵抗、容量性、および誘導要素14の組み合わせを通して心血管系の全体的な血行力学的挙動を模倣する。このグループの最近の研究では、従来の電気アナログモデルよりも直感的な方法で、大きな血管-心臓室およびバルブの形状と力学の変化をモデル化することを可能にする油圧ドメインの代替Windkesselモデルを提案しています。このシミュレーションは、オブジェクト指向の数値ソルバー( 材料表を参照)で開発され、正常な血行力学、心呼吸結合の生理学的影響、単心生理学における呼吸駆動血流、および大動脈収縮による血行力変化を捉えることができます。この記述は、心不全15を含む病的状態のスペクトルをモデル化するための物理的に直感的なアプローチを提供することにより、一括パラメータモデルの能力を拡張する。
高次元モデルはFEAに基づいて、時空間的な力学と流体構造の相互作用を計算します。これらの表現は、局所血流分野の詳細かつ正確な記述を提供することができます。しかし、その低い計算効率のために、彼らは全体の心血管樹16、17の研究には適していません。ソフトウェアパッケージ(材料表を参照)は、電気機械応答、構造変形、流体キャビティベースのヘモダイナミクスを統合した4室の成人人心臓の解剖学的に正確なFEAプラットフォームとして採用されました。適応されたヒト心臓モデルはまた、異なる流体腔間の流れ交換を定義する単純な一括パラメータモデルと、心臓組織18、19の完全な機械的特徴付けを含む。
心不全のいくつかの一括パラメータおよびFEAモデルは、血行力学的異常を捕捉し、治療戦略を評価するために、特にHFrEF20、21、22、23、24のための機械的循環補助装置の文脈で策定されている。したがって、様々な複雑さの0D一括パラメータモデルの広範な配列は、20、21、23、24の2または3要素電気アナログウィンドケセルシステムの最適化を介して生理学的およびHFrEF条件における人間の心臓のヘモダイナミクスを正常に捕捉している。これらの表現の大部分は、時間変動エラストアンス製剤に基づく単一または二心体モデルであり、心臓の収縮作用を再現し、非線形のエンド拡張期圧力容積関係を使用して心室充填25、26、27を記述する。複雑な心血管ネットワークを捉え、心房および心室ポンプ作用の両方を模倣する包括的なモデルは、デバイステストのプラットフォームとして使用されてきました。それにもかかわらず、HFrEFの分野の周りに重要な文献体が存在するが、HFpEFのインリコモデルでは非常に少数が提案されている20、22、28、29、30、31。
最近、バーコフら32とグラネガーら28によって開発されたHFpEFヘモダイナミクスの低次元モデルは、4室の心臓の圧力体積(PV)ループを捕捉することができ、HFpEFの様々な表現型のヘモダイナミクスを完全に再現する。さらに、インシリコプラットフォームを利用して、HFpEF用の機械循環器の実現可能性を評価し、HFpEFの生理学研究やデバイス開発のための計算研究の先駆的な研究を行っています。しかし、これらのモデルは、病気の進行中に観察された血液の流れと圧力の動的変化を捉えることができないままです。Kadryらら30の最近の研究では、心筋の活発な弛緩と左心室の受動的な剛性を低次元モデルで調整することにより、様々な表現型の拡張期機能不全を捉えています。彼らの研究は、心筋の能動的および受動的特性の両方に基づいて拡張期機能不全の包括的な血行力学分析を提供する。同様に、高次元モデルの文献は、主にHFrEF19、33、34、35、36、37に焦点を当てています。Bakirら33は、HFrEF血動力学的プロファイルと左心室補助装置(LVAD)の有効性を予測する完全結合心液電気力学FEAモデルを提案した。このバイベンチュラ(または2室)モデルは、健康な心臓のヘモダイナミクスをシミュレートするために結合されたウィンドケッセル回路を利用し、LvADサポート33、37を有するHFrEFおよびHFrEF。
同様に、Sackら35は右心室機能障害を調査する二心筋モデルを開発した。彼らの二心学的幾何学は患者の磁気共鳴画像法(MRI)データから得られ、モデルの有限要素メッシュは、VAD支持の失敗した右心室35の血行力学を解析するために画像セグメンテーションを用いて構築された。4室FEA心臓アプローチは、心臓19,34の電気機械挙動のモデルの精度を高めるために開発された。二心筋の記述とは対照的に、人間の心臓のMRI由来4室モデルは、心血管解剖学18のより良い表現を提供する。この研究で採用された心臓モデルは、4室のFEAモデルの確立された例である。一括パラメータおよびバイベンチュラFEAモデルとは異なり、この表現は、疾患進行時に発生する血行力学的変化を34,37にキャプチャします。Genetら34は、例えば、HFrEFおよびHFpEFで観察された改造の数値成長モデルを実施するために同じプラットフォームを使用した。しかし、これらのモデルは、構造力学に対する心臓肥大の影響のみを評価し、関連する血行力学の包括的な説明を提供しない。
この研究におけるインシリコモデルにおけるHFpEFの不足に対処するために、このグループ15 およびFEAモデルによって以前に開発された一括パラメータモデルは、HFpEFの血行力学的プロファイルをシミュレートするために再適応された。このため、ベースラインで心血管のヘモダイナミクスをシミュレートする各モデルの能力が最初に実証されます。狭窄によって引き起こされた左心室圧過負荷の影響と、心臓リモデリングによる左心室コンプライアンスの低下-HFpEF-の典型的な特徴が評価される。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. 0D 集中パラメータモデル
- シミュレーションの設定
注: 数値ソルバ環境 (図1に示すように、表を参照) で、ドメインを構築します。これは、4室の心臓、上半身、腹部、下半身、胸部の区画、大動脈、肺動脈、および上側および下腹部のカバエを含む近位血管系で構成される。このシミュレーションで使用される標準要素は、デフォルトの油圧ライブラリの一部です。詳細については、「補足ファイル」を参照してください。- 油圧ライブラリをナビゲートして、必要な要素を見つけるために:油圧パイプライン、定量油圧チャンバー、線形抵抗、遠心ポンプ、チェックバルブ、可変面積オリフィス、およびカスタム油圧流体。
- 油圧パイプライン要素をワークスペースにドロップします。
注:これらは、摩擦損失だけでなく、血管や心臓室の壁のコンプライアンスと流体圧縮率を説明します。このブロックを通じて、圧力損失はダルシー・ヴァイスバッハ法を使用して計算され、壁のコンプライアンスによる直径の変化は、コンプライアンス比例定数、光圧、および時間定数に依存します。最後に、流体の圧縮性は、媒体のバルク弾性率によって定義される。 - 固定体積油圧チャンバ要素を挿入して、壁のコンプライアンスと流体の圧縮性を定義します。
注: このブロックは摩擦による圧力損失を考慮していません。 - 線形抵抗要素を追加して、流れに対する抵抗を定義します。
注: これは、電気アナログ Windkessel モデルで使用される抵抗要素に似た、脈管構造の幾何学的特性とは無関係です。遠心ポンプ、チェックバルブ、可変面積オリフィス、およびカスタム油圧流体要素などの他のブロックは、システムに所望の圧力入力を生成し、血流に対する心臓弁の影響をモデル化し、血液の機械的特性を定義するために挿入する必要があります。これらの要素を通じて、生理学と疾患の両方における心血管系の挙動を完全に捉えることができる。遠心ポンプの入力信号は 図S1Aにあります。 - カスタムの可変コンプライアンスコンプライアンスチャンバー要素を通じて、各心臓室の収縮性をモデル化します。
注: これは、時間変動のユーザ定義入力信号としてコンプライアンスを受け入れ、時間変動エラストアンスモデル(図S1B-D)に基づいています。
- 油圧パイプライン要素をワークスペースにドロップします。
- 表 S1に示すように、各要素に対するパラメータを指定します。
- 時間変動するユーザー定義入力信号(LVポンプ、可変コンプライアンスコンプライアンス要素、可変領域オリフィスブロック)を必要とする各ブロックに物理信号(PS)繰り返しシーケンス要素を挿入します。
注:このシミュレーションに使用される入力信号は 、図S1にあります。 - デフォルトの ODE 23t 暗黙的ソルバーを 選択し、100 s のシミュレーションを実行して、定常状態に達します。
- 油圧ライブラリをナビゲートして、必要な要素を見つけるために:油圧パイプライン、定量油圧チャンバー、線形抵抗、遠心ポンプ、チェックバルブ、可変面積オリフィス、およびカスタム油圧流体。
2. FEA モデル
- シミュレーションの設定
注: FEA モデルは、連続した電気機械解析を使用します。このモデルでは、電気解析が最初に行われます。次に、得られた電位を、次の機械解析で励起源として使用します。したがって、シミュレーション設定には、FEA シミュレーション ソフトウェア (材料表)18で事前定義されている電気(ELEC) ドメインとメカニカル ドメインの 2 つの作業領域が含まれています。したがって、次のセクションでは、分析ワークフローのみを説明します。FEA モデルでは、電気および機械材料モデリング18に対して、次のユーザ サブルーチンHETVAL、VUANISOHYPER、および UAMPを使用します。- 標準モジュールで事前定義された温度手順を使用して電気解析を実行するには、ELECドメインをナビゲートします。
- BEATという名前の単一分析ステップを使用します。心周期の持続時間を500 msに設定し、sinatrial(SA)ノードを表すノードセットに電位パルスを適用する(ノードセット: R_Atrium-1.SA_NODE)。
- モデルガイド18に記載されているように、200ms以上の-80 mVから20mVまでのデフォルトの電気波形を、滑らかなステップ振幅定義で確認します。電気解析で材料定数の既定値を使用して、AV 遅延を調整します。
- Jobモジュールを起動し、 heart-elecという名前のジョブを作成します。
- 電気解析の設定が完了したら 、MECH 領域をナビゲートして流体キャビティベースの機械的解析を実行します。
注: 機械的なシミュレーションは電気解析の後に実行され、結果として得られる電位は機械解析の励起源として使用されます。機械解析には複数のステップが含まれています。- PRE-LOAD、BEAT1、およびRECOVERY1という名前の 3 つの主要なステップを使用します。PRE-LOADステップで、心臓のプレストレス状態の境界条件を確認します。ステップ時間として0.3 sを使用して、流体チャンバ内の圧力を直線的に上昇させます。
注: 定義済みの流体キャビティ圧力値を 表 S3に示します。心臓のプレストレス状態は、通常の心臓シミュレーション設定ですでに定義されており、初期ノードの状態は 、表S5に示すように、外部シミュレーションファイルに提供されています。ステップ 3.2.2~3.2.4 で説明されているように、逆の機械的シミュレーションを使用したゼロ応力状態の再計算は、境界条件が変更されるたびに必要です。 - BEAT1ステップでは、収縮をシミュレートするステップ時間として 0.5 s を使用します。
- RECOVERY1ステップでは、心臓の弛緩と心室充填の 0.5 s を選択し、心拍数は 60 bpm です。
- 後続のステップである BEATX と RECOVERYXを有効にして、複数の心周期をシミュレートして定常状態に達します。
注:3つの心周期は定常状態に達するのに十分です。シミュレーションの1サイクルは、24コアプロセッサ(3.2 GHz x 24)で約8時間で完了します。 - Jobモジュールを起動し、ハートメカという名前のジョブを作成し、倍精度オプションを有効にします。
- PRE-LOAD、BEAT1、およびRECOVERY1という名前の 3 つの主要なステップを使用します。PRE-LOADステップで、心臓のプレストレス状態の境界条件を確認します。ステップ時間として0.3 sを使用して、流体チャンバ内の圧力を直線的に上昇させます。
- 標準モジュールで事前定義された温度手順を使用して電気解析を実行するには、ELECドメインをナビゲートします。
- 単純化された一括パラメータの Windkessel モデルを確認します。
注: FEA モデルの機械的領域は、単純化された一括パラメータ回路に基づく血流モデルを持ち、 図 218に示すように、表面ベースの流体キャビティと流体交換の組み合わせとして作成されます。- 上記の注記で説明した Windkessel 表現を使用してシミュレーションを実行します。血流モデルの表現を確認して、流れ抵抗と構造コンプライアンスに対する抵抗要素と静電容量性要素の値をそれぞれ調整します。
- 4つの心室の3D有限要素表現を確認し、その幾何学的位置が正確であることを確認します。
- 心臓アセンブリを確認し、 インタラクション モジュールに切り替えて、4つの心臓室のそれぞれのコンプライアンスと収縮性の値を調整します。
注: インタラクション モジュールのデフォルト値は、理想化された健康な人間の心臓の鼓動サイクル18をシミュレートするように設定されています。 - 相互作用モジュール、CAV-AORTA、CAV-LA、CAV-LV、CAV-PULMONARY_TRUNK、CAV-RA、CAV-RV、CAV-SVC、CAV-動脈コンプ、CAV-肺コンプ、およびCAV-VENOUS-COMP(表S3)における次の静水性流体キャビティを確認する。
- コンプライアンスチャンバ(CAV-動脈コンプ、CAV-肺コンプ、およびCAV-VENOUS-COMP)を立方体容積として使用し、動脈、静脈、肺循環のコンプライアンスを表します。
- 3 つのコンプライアンス立方体容積を接地スプリングに取り付け、剛性値を確認して、動脈、静脈、肺循環における圧力容積応答をモデル化します。
- 静水性流体腔間の流体交換定義を確認してください: 動脈静脈、静脈右心房、右心室右心室、右心室肺系、肺系左心房、左心室左心室、および左心室-大動脈(表S4)。
- 粘性抵抗係数を調整して、各流体交換リンクの血流モデルを変更します(粘性抵抗効果の詳細については 、補足ファイル を参照)。
- マルチフィジックスシミュレーション
- 作業ディレクトリで CAE データベース ファイルを見つけます。
注: このプロトコルの FEA モデルはデータベースに配信され 、LH-ヒューマン モデル-β-V2_1.caeと名付けられています。 - 入力ファイル、オブジェクト ファイル、ライブラリ ファイルを作業ディレクトリに挿入してシミュレーションを実行します。入力ファイルとライブラリファイルの完全なリストについては、 表 S5 を参照してください。
- FEA モデル シミュレーション ソフトウェアを起動します( 材料表を参照)。
注 : ソフトウェア プロバイダに問い合わせ、新しいバージョン18との互換性を確認してください。 - セクション 2.2 および 2.3 で説明されているように 、ELEC ドメインと MECH ドメインの両方の部品、アセンブリ、および境界条件を確認します。
- まず、セクション 2.1.1.3 で説明されているように 、heart-elecという名前の電気シミュレーション ジョブを実行します。電位の結果を視覚的に調べて、 心電位 シミュレーションが期待どおりに実行されたことを確認します。次に、結果ファイル heart-elec.odb が作業ディレクトリにあることを確認します。
- MECHドメインに切り替えて、2番目のシミュレーションフェーズに移行します。機械的シミュレーションで使用される材料定数の値を確認して、希望する受動および活動的な心応答をモデル化します。
- 機械解析用の材料ライブラリ ファイルで、HYBRID-文字列名が使用されていることを確認します。心臓室の材料応答を変更するには、適切なハイブリッド材料ファイルを調整するか、CAEモジュールの材料セクションで新しい材料の挙動を定義して材料全体の応答を置き換えます。
注:組み込みの構成法に関する詳細な情報は、ユーザーガイド18に記載されています。 - PRE-LOADステップで、静水キャビティの圧力を設定して、所望の生理学的挙動を得ます。ステップ 2.1.2.1 で説明したように、ゼロから希望の圧力レベルに上昇するには、組み込みのスムーズ振幅オプションを使用します。
- 循環系内の一定の全体的な血液量で血流モデルを実行するには、2.1.2.1 で定義された圧力境界条件を無効にします。セクション 2.1.2.5 で説明されているように 、heart-mechという名前のシミュレーション ジョブを実行します。
- 作業ディレクトリで CAE データベース ファイルを見つけます。
3. 大動脈弁狭窄症
注:大動脈狭窄は、圧力過負荷につながり、最終的には左心室壁の同心円状のリモデリングとコンプライアンス損失につながるため、HFpEFのドライバーであることが多いです。大動脈狭窄症で観察される血行力学は、しばしばHFpEFに見られるものに進行する。
- 一括パラメータ モデル
- PS繰り返しシーケンス要素の入力信号を、左心室コンパートメントにある大動脈弁に対して変更します。ベースラインと比較して70%に等しいオリフィス領域の減少をシミュレートする(表S6)。
注: 入力値は、各ハートビートの間に、ステノミックバルブのオリフィス領域を表します。オリフィス面積値は、大動脈弁PS要素の 開始出力値 ベクトルに、元の値に対する最終的なオリフィス面積に対応する10進値を掛けることによって容易に調整することができる。この研究では、70%の収縮を達成するために0.3の係数を使用した。
- PS繰り返しシーケンス要素の入力信号を、左心室コンパートメントにある大動脈弁に対して変更します。ベースラインと比較して70%に等しいオリフィス領域の減少をシミュレートする(表S6)。
- FEA モデル
- LINK-LV-ARTERIALパラメータの流体交換定義を変更します。
注:このパラメータは、左心室と大間の間の血流に合わせた粘性抵抗係数を有します。有効交換領域を変更して血流を調整し、適切な大動脈狭窄モデルを作成することができる(表S7)。 - ツールボックス フォルダーを見つけて、そのフォルダー内のファイルをメイン作業ディレクトリにコピーします。
- ツールボックスファイル18を実行して逆加工シミュレーションを実行します。この目的のために、左心室と左心房の吸引圧力を流体腔内の6mmHgに変更し、大動脈狭窄モデルの初期容積状態を調整する。 inversePreliminary.py 関数を実行します。
注:境界条件を変更する場合は、必ず逆加工シミュレーションを使用してゼロ応力状態を再計算する必要があります。 - 逆の機械的シミュレーションが完了したら、後処理機能を実行します: calcNodeCoords.py と straight_mv_chordae.py。他の流れパラメータのデフォルト値を使用し、セクション2.1.2.5で説明されているように新しい機械シミュレーションを実行します。
- LINK-LV-ARTERIALパラメータの流体交換定義を変更します。
4. HFpEFヘモダイナミクス
メモ: 慢性リモデリングの効果をシミュレートするために、左心の機械的特性が修正されました。
- 一括パラメータ モデル
- LVコンプライアンス要素の左心室拡張期コンプライアンスを変更して、圧力過負荷による壁硬化を模倣し、 表S8の末方拡張期コンプライアンスの値を使用する。
注:コンプライアンスがエンドシトールからエンドダイストールに直線的に落ちると仮定します。 - HFpEFで観察された高い左心室圧を捕捉するために、LVポンプのリーク抵抗を106 Pa sm-3(表S8)18×に上げます。
- LVコンプライアンス要素の左心室拡張期コンプライアンスを変更して、圧力過負荷による壁硬化を模倣し、 表S8の末方拡張期コンプライアンスの値を使用する。
- FEA モデル
- 左心室ジオメトリの有効な材料特性を編集します。剛性成分を増やして、構成モデルの繊維およびシート方向の応力成分に影響を与える活性組織応答を調整します。
- メカマットLV_ACTIVEファイルの左心室の材料応答を修正します。
注: 左心室の剛性の大きさは適切な拡張期の承諾効果を提供するために調整することができる。 - 異方性超弾性製剤の剛性パラメータ a および b を増やし、HFpEF 生理学に対する剛性応答の増加を捉えます。
- PRE-LOADステップでは、左心室と左心房の流体空洞圧力を20mmHgに設定します。
- 逆の機械的シミュレーションを実行して、左心室と心房の容積状態を取得します。 ハートメック-inverse.odb ファイル18から節点座標をエクスポートします。
- ステップ 3.2.4 の説明に従って、後処理機能 calcNodeCoords.pyとstraight_mv_chordae.pyを実行します。新しい節点入力ファイルを作業ディレクトリ内で見つけ、セクション 2.1.2.5 で説明されているように、新しい機械シミュレーションを実行します。
- メカマットLV_ACTIVEファイルの左心室の材料応答を修正します。
- 左心室ジオメトリの有効な材料特性を編集します。剛性成分を増やして、構成モデルの繊維およびシート方向の応力成分に影響を与える活性組織応答を調整します。
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
ベースライン シミュレーションの結果を 図 3に示します。左心室と大オータ(図3A)の圧力と体積波形、左心室PVループ(図3B)を示しています。インシリコモデルの2つは、生理学的範囲内にある同様の大動脈および左心室血行性を示す。2つのプラットフォームによって予測される応答のわずかな違いは、非線形性が一括パラメータプラットフォームと比較してFEAモデルによってより良く捕捉される心室の空と充填段階の間に気づくことができます。生理学では、このような非線形効果は主に心臓組織の超弾性応答の結果として生じ、したがって、マルチドメインおよび高次計算モデル18によってより正確に再現される。
大動脈狭窄症のために心室および大動脈血行態が得られたが、これはしばしば左心室圧過負荷につながり、最終的にはHFpEFの発達につながる。大動脈弁オリフィス面積の70%減少での圧力および体積波形は、両方のモデルに対して図4に示されている。狭窄は大動脈弁を横切る圧力勾配の上昇をもたらした。この研究で検討した70%の狭窄について、41mmHgおよび54mmHgのピークトランスアオルト圧力勾配を、それぞれ塊状パラメータ(図4A)およびFEA(図4B)モデルで得た。この中程度の変動は、コンプライアンスが単に数値の配列によって定義される、一括パラメータモデルにおける心臓組織の物質特性を定義する構成方程式の欠如の別の結果として生じる可能性が高い。したがって、このモデルは、FEA モデルによって正確に表現される流体と構造の相互作用をキャプチャしません。それにもかかわらず、両方のモデルの結果は、米国心エコー学会(ASE)およびヨーロッパ心エコー学会(EAE)中等大動脈弁狭窄の分類と一致しており、これは約60〜75%の大動脈収縮に対する40〜65mmHgのピーク経大動脈勾配を示す。
左心室PVループはベースラインで、70%大動脈狭窄、および心室壁の硬化後のHFpEFのを図5にまとめた。同様のパターンは、一括パラメータモデルの結果を示す図5Aと、FEAを介して得られたヘモダイナミクスを示す図5Bで観察することができます。これらのPVループは、HFpEF1、11、28、32の科学および臨床文献のものと一致しています。特に、両方のモデルは大動脈狭窄によって誘発された後負荷の上昇による収縮期左心室圧の増加を捕捉することができる。さらに、狭窄PVループにおいてエンドシストリック体積が増加し、ストローク体積の低下を招く。左心室コンプライアンスの改造と喪失時に、エンド拡張期圧力容積関係(EDPVR)が上昇し、末拡張期圧力が高くなり、末端拡張期の体積が下がってしまう。これらの現象は、左心室が十分にリラックスして充填できないことに起因し、低次元モデルと高次元モデルの両方でHFpEF PVループによって正常に捉えられる。
拡張期機能の低下を示すもう1つの指標として、僧帽弁を通る流れが、初期弛緩(E)および心房収縮(A)相の両方を強調する 図S2に示されている。正常および狭窄プロファイルと比較して、HFpEF流動は、わずかに高いピークE相僧帽弁流と有意に減少したピークA相流によって特徴付けられるが、左心室の受動的な硬化は、科学文献30と一致する上昇したE/A比を生じることを強調する。最後に、 図6 は、収縮期および拡張期の両方における正常およびHFpEF心における心筋ストレスマップの変化を示す。左心室の長軸図は、体積平均応力分布を示し、心室コンプライアンスの特徴的な喪失によるHFpEF心臓の高い応力を示す。ピーク収縮期(t = 0.2 s)および末端拡張(t = 1.0 s)の間の健康な心臓のための(61.1 ±49.8)kPaおよび(0.51 ±7.35)kPaのベースライン値から、 HFpEFでは(97.2±205.7)kPaおよび(2.69±16.34)kPaに対応して平均応力が増加し、HFpEFで観察された血行力変化が、障害のある心臓に影響を及ぼす深い構造変化に根ざしていることを示唆している。
図1:オブジェクト指向数値ソルバーにおける解剖学的に導出された一括パラメータモデルの領域(材料表参照)は、4室の心臓、大動脈、および上半身、腹部、下半身、肺循環を示す。略語: LV = 左心室;RV = 右心室;LA = 左心房;RA = 右心房;R1 = 動脈抵抗;R2 = 静脈抵抗;C = コンプライアンス;IVC:下の大静脈。SVC:優れた大静脈。この図の大きなバージョンを表示するには、ここをクリックしてください。
図2:人の心臓の有限要素解析モデル(A)ヒト心臓の有限要素解析モデルの3D表現(B) 構造流体交換モデル18と結合されたモデルにおける血流モデルの単純化された一括パラメータ表現。略語: LV = 左心室;RV = 右心室;LA = 左心房;RA = 右心房;R大動脈=大動脈弁抵抗;R僧帽弁=僧帽弁抵抗;R肺=肺弁抵抗;R三尖=三尖弁抵抗;C動脈=全身動脈コンプライアンス;R系=全身動脈抵抗;C静脈=全身静脈コンプライアンス、R静脈=全身静脈抵抗;C肺=肺コンプライアンス;R肺系=肺抵抗性。この図の大きなバージョンを表示するには、ここをクリックしてください。
図3:人間の心臓の塊状パラメータと有限要素解析モデルのベースラインシミュレーションと圧力体積波形(A)左心室圧と体積波形と、ベースラインでの一括パラメータとFEAモデルによって計算された大動脈圧。(B)ベースラインで両プラットフォームを通して得られる左心室PVループ。略語: FEA = 有限要素解析;LV =左心室;PV = 圧力量。この図の大きなバージョンを表示するには、ここをクリックしてください。
図4:左心室圧と体積波形と大動脈弁の大動脈の減少で計算された大動脈圧は、大動脈弁オリフィス領域の70%減少で計算した。略語: FEA = 有限要素解析;LV = 左心室。この図の大きなバージョンを表示するには、ここをクリックしてください。
図5:健康な心臓の左心室PVループ、急性狭窄誘発圧過負荷下、および慢性リモデリングおよび硬化後のHFpEF心臓の(A)塊状パラメータ(B)FEAモデル。略語: EDPVRH = シミュレートされた健康な心臓におけるエンド拡張期の圧力量関係;EDPVRHFpEF: シミュレーションされた HFpEF 生理学におけるエンド拡張期の圧力容積関係;PV - 圧力量。FEA = 有限要素解析。この図の大きなバージョンを表示するには、ここをクリックしてください。
図6: フォンミーゼス応力(平均: 75%)FEAモデルによって予測されるように、生理学的状態およびピーク収縮期および拡張期のHFpEF心臓の下で。 カラーマップは、MPa の応力レベルを示します。HFpEF(92.7-2.7 kPa)では、ピーク収縮期(t = 0.2 s)およびエンドダイストール(t = 1.0 s)の間の健康な心臓(61.1-0.5 kPa)と比較して、より高いストレスが見られます。 この図の大きなバージョンを表示するには、ここをクリックしてください。
図S1:(A)遠心ポンプ、(B)左心室、(C)右心室、(D)左右心房の入力信号を、一括パラメータシミュレーション用に示す。 このファイルをダウンロードするには、ここをクリックしてください。
図S2:(A)大動脈および(B)、FEAによって得られたベースライン、狭窄、およびHFpEFプロファイルの僧帽弁流信号。略語: E = 早期緩和段階;A = 心房収縮;FEA = 有限要素解析;HFpEF = 保存された駆出率を有する心不全。 このファイルをダウンロードするには、ここをクリックしてください。
表 S1. ベースラインの一括パラメータシミュレーションの幾何学的パラメータと機械的パラメータ。 このテーブルをダウンロードするには、ここをクリックしてください。
表 S2.ベースラインの一括パラメータシミュレーションのパラメータの広範なセット。 このテーブルをダウンロードするには、ここをクリックしてください。
表 S3. 機械的有限要素解析(FEA)モデル18における流体キャビティ値。 このテーブルをダウンロードするには、ここをクリックしてください。
表 S4. 有限要素解析(FEA)モデル18における流体交換リンクの境界条件。 このテーブルをダウンロードするには、ここをクリックしてください。
表 S5. 有限要素解析 (FEA) モデル18に必要なシミュレーション ファイル。 このテーブルをダウンロードするには、ここをクリックしてください。
表 S6. 大動脈狭窄集中パラメータシミュレーションのパラメータ。 このテーブルをダウンロードするには、ここをクリックしてください。
表 S7. 有限要素解析(FEA)モデル18における流体交換リンク定義。 このテーブルをダウンロードするには、ここをクリックしてください。
表 S8. HFpEF の一括パラメータ シミュレーションのパラメータ。 このテーブルをダウンロードするには、ここをクリックしてください。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
この研究で提案された一括パラメータおよびFEAプラットフォームは、狭窄誘発圧過負荷の急性期および慢性HFpEFの両方において、生理学的状態下での心血管系血液力学を再現した。HFpEFの急性および慢性段階における圧力過負荷が果たす役割を捉えることによって、これらのモデルの結果は、大動脈狭窄による経大動脈圧勾配の発症、左心室圧の増加、壁面硬化による末期拡張期体積の減少を含むHFpEFの臨床文献と一致している。さらに、このFEAモデルは、心臓周期を通してHFpEFにおける心筋ストレスの上昇を捉えることができた。これらのシミュレーションを正しく設定するには、上記のプロトコルセクションで概説した手順を厳格に実行する必要があります。集中パラメータモデルの場合、 図1 に示すように油圧要素のネットワークを正しく再作成し、所定の値を入力パラメータとして提供することが重要です(表S1 および 表S2)。さらに、ソルバ ブロックは、任意のノードでネットワークに対して定義および接続されている必要があります。
FEA モデルを機能させるには、表 S5にリストされているソルバ18と共にパッケージ化されたすべてのシミュレーション ファイルが必要です。前提条件コンポーネントが省略されている場合、シミュレーションが早期終了する可能性があります。どちらのプラットフォームでも、狭窄とHFpEF血動プロファイルを再作成する前に、デフォルトの入力パラメータを使用してベースラインシミュレーションを取得することが重要です。ベースライン シミュレーション15と補足ファイルのシミュレーションにリンクされたドキュメントの概要を示す元の研究記事を参照して、一括パラメータ モデルのトラブルシューティングを行うことができます。同様に、この FEA フレームワークには、トラブルシューティング18のためのソフトウェア ドキュメントとツールボックス フォルダが含まれています。シミュレーションエラーが発生した場合、ユーザはツールボックスフォルダ18内の相対プラグインを実行してシミュレーション診断を呼び出すことができます。一括パラメータモデルからの血行力学結果は、各シミュレートされた条件でFEAを介して計算されたものと類似しており、HFpEFの臨床文献と一致した。高次元FEAプラットフォームは、心臓の複雑な生体力学的挙動を捕捉することを可能にし、高められた計算需要を犠牲にして、心血管の力学の正確な記述を提供する。ただし、一括パラメータモデルでは、ランタイムは数時間から数分に短縮され、シリコモデルの上位に比べ、大きな利点を構成します。
さらに、より多くの心血管コンパートメントをモデル化することにより、この一括パラメータシミュレーションは、心血管樹の様々な部位における血液の流れおよび圧力の検査を可能にし、したがって、心臓室および近位脈管系を超えて広がる研究に適している。しかし、世界的なヘモダイナミクスを再現できる一方で、この記述は構造的相互作用の軽微な影響を捉えることができないため、FEA表現に典型的な精度が欠けている。有限要素アプローチを通じて本研究で得られた心臓力学の分析は、以前の調査からそれらを裏付けた。具体的には、これらの平均応力値は、慢性心不全34,37の間に部分的に支持された心臓の成長モデルによって予測されるものと同じ範囲にある。これらのモデルと比較して、これらの研究で見つかった応力値は、圧力過負荷を誘発するためにシミュレートされた大動脈狭窄の上昇レベルのために適度に高かった。さらに、HFpEFにおける左心室コンプライアンスの喪失は、心内膜ストレスに大きな影響を与えることがわかった。
しかし、この研究では、拡張期の剛性とその感受性はパラメトリックに調査されなかった。実際、このパラメータは、慢性左心室圧過負荷の生理学的に関連する血行性プロファイルを捕捉するように調整された。拡張期のコンプライアンスの低下の影響を十分に特徴付けるために、広範な感度分析を行う必要があります。この計算モデルは、HFpEFにおける心臓構造の生体力学的変化がリモデリングの主要な原動力となり、HFpEFのヘモダイナミクスおよび疾患進行に大きな影響を及ぼす可能性があることをさらに示唆している。FEAシミュレーションの流体構造相互作用との動的成長モデルの統合は、圧力過負荷によって誘発される心臓リモデリングおよび血行力収差のダイナミクスをより包括的に捕捉する今後の研究で考慮される可能性がある。さらに、拡張期機能障害の異なるフェノタイプをシミュレートするために、Kadryらら30 と同様の活性弛緩の効果に関するさらなる研究と電気伝導および収縮性が必要となる可能性がある。
HFpEFの研究に適したシミュレーションプラットフォームの開発は、文献ではほとんど報告されていません。この文脈において、この研究はHFpEF病態生理学の研究のための独特な環境を提供する。解剖学的に導出された一括パラメータモデルは、患者固有の血行力学パラメータ(例えば、血管の発光領域およびコンプライアンス)が健康およびHFpEF状態のための世界的な血行力学で遊ぶ効果の迅速なシミュレーションを可能にする。さらに、FEAモデリングにより、HFpEF中に心臓組織が徐々に変化する際の機械的特性および興奮性の一時的変化の影響を詳細に調査することが可能です。さらに、提案されたモデルは、HFpEFの新規治療法のシミュレーションのための潜在的な有用性を有し、一部は、生体内、インビトロ、 および不十分なデバイス最適化による臨床試験の中断を担当する可能性のあるHFpEFのシリコモデルの欠如に対処する。最後に、FEA アプローチの基礎となる単純化された一括パラメータ記述を数値ソルバー モデルに置き換えることで、これらのモデルを 1 つのシミュレーションに統合することが今後の作業に伴う場合があります。これにより、これらのモデルの精度がさらに向上し、HFpEFおよび他の心血管状態の計算研究をさらにサポートする可能性があります。
要約すると、HFpEFの2つの異なる計算モデルが本研究で説明された。開発されたプラットフォームが生理学的条件下でベースラインの血液力学を記述する能力が最初に実証された。その後、大動脈狭窄から生じる変化を、最終的には左心室リモデリングによるHFpEFから発生する変化を調査し、その結果が文献で報告されたものと一致していることを実証した。最後に、シミュレートされた血行力学的状態は、生理学的状態と比較してHFpEF心臓の心臓壁ストレスの上昇を示した。HFpEFが表す信じられないほど差し迫ったヘルスケアの課題の文脈では、これらの提案されたプラットフォームは、HFpEFのヘモダイナミクスとバイオメカニクスに関する洞察を提供することができるシリコ記述の最初の一つです。これらの計算モデルは、HFpEFの治療法の開発のためのツールとしてさらに使用され、最終的にはこの分野における翻訳研究を支援する。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
この作業に関連する利益相反はありません。
Acknowledgments
ハーバード・マサチューセッツ工科大学の健康科学技術プログラムと、医工学科学研究所のSITA財団賞からの資金提供を認めます。
Materials
Name | Company | Catalog Number | Comments |
Abaqus Software | Dassault Systèmes Simulia Corp. | Version used: 2018; FEA simulation software | |
HETVAL | Dassault Systèmes Simulia Corp. | Version used: 2018 | |
Hydraulic (Isothermal) library | MathWorks | Version used: 2020a | |
Living Heart Human Model | Dassault Systèmes Simulia Corp. | Version used: V2_1, anatomically accurate FEA platform of 4-chamber adult human heart | |
MATLAB | MathWorks | Version used: 2020a, object-oriented numerical solver | |
SIMSCAPE FLUIDS | MathWorks | ||
UAMP | Dassault Systèmes Simulia Corp. | Version used: 2018 | |
VUANISOHYPER | Dassault Systèmes Simulia Corp. | Version used: 2018 |
References
- Borlaug, B. A., Paulus, W. J. Heart failure with preserved ejection fraction: Pathophysiology, diagnosis, and treatment. European Heart Journal. 32 (6), 670-679 (2011).
- Borlaug, B. A., Kane, G. C., Melenovsky, V., Olson, T. P. Abnormal right ventricular-pulmonary artery coupling with exercise in heart failure with preserved ejection fraction. European Heart Journal. 37 (43), 3293-3302 (2016).
- Borlaug, B. A. Evaluation and management of heart failure with preserved ejection fraction. Nature Reviews Cardiology. 17 (9), 1-15 (2020).
- Carabello, B. A., Paulus, W. J.
Aortic stenosis. The Lancet. 373 (9667), 956-966 (2009). - Lam, C. S. P., Donal, E., Kraigher-Krainer, E., Vasan, R. S. Epidemiology and clinical course of heart failure with preserved ejection fraction. European Journal of Heart Failure. 13 (1), 18-28 (2011).
- Omote, K., et al. Left ventricular outflow tract velocity time integral in hospitalized heart failure with preserved ejection fraction. ESC Heart Failure. 7 (1), 167-175 (2020).
- Samson, R., Jaiswal, A., Ennezat, P. V., Cassidy, M., Jemtel, T. H. L. Clinical phenotypes in heart failure with preserved ejection fraction. Journal of the American Heart Association. 5 (1), (2016).
- Weber, K. T., Brilla, C. G., Janicki, J. S. Myocardial fibrosis: Functional significance and regulatory factors. Cardiovascular Research. 27 (3), 341-348 (1993).
- Borbély, A., et al. Cardiomyocyte stiffness in diastolic heart failure. Circulation. 111 (6), 774-781 (2005).
- Borlaug, B. A., Lam, C. S. P., Roger, V. L., Rodeheffer, R. J., Redfield, M. M. Contractility and Ventricular Systolic Stiffening in Hypertensive Heart Disease. Insights Into the Pathogenesis of Heart Failure With Preserved Ejection Fraction. Journal of the American College of Cardiology. 54 (5), 410-418 (2009).
- Penicka, M., et al. Heart Failure With Preserved Ejection Fraction in Outpatients With Unexplained Dyspnea. A Pressure-Volume Loop Analysis. Journal of the American College of Cardiology. 55 (16), 1701-1710 (2010).
- Owen, B., Bojdo, N., Jivkov, A., Keavney, B., Revell, A.
Structural modelling of the cardiovascular system. Biomechanics and Modeling in Mechanobiology. 17 (5), 1217-1242 (2018). - Zhou, S., et al. A review on low-dimensional physics-based models of systemic arteries: Application to estimation of central aortic pressure. BioMedical Engineering Online. 18 (1), 41 (2019).
- Sagawa, K., Lie, R. K., Schaefer, J. Translation of Otto frank's paper "Die Grundform des arteriellen Pulses" zeitschrift für biologie 37. Journal of Molecular and Cellular Cardiology. 22 (1899), 253-254 (1990).
- Rosalia, L., Ozturk, C., Van Story, D., Horvath, M., Roche, E. T. Object-oriented lumped-parameter modeling of the cardiovascular system for physiological and pathophysiological conditions. Advanced theory and simulations. , (2021).
- Lopez-Perez, A., Sebastian, R., Ferrero, J. M. Three-dimensional cardiac computational modelling: METHODS, features and applications. BioMedical Engineering Online. 14, 35 (2015).
- Xie, X., Zheng, M., Wen, D., Li, Y., Xie, S. A new CFD based non-invasive method for functional diagnosis of coronary stenosis. BioMedical Engineering Online. 17 (1), 36 (2018).
- Abaqus Dassault, S. SIMULIA living heart human model user documentation. , (2017).
- 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).
- Moscato, F., et al. Use of continuous flow ventricular assist devices in patients with heart failure and a normal ejection fraction: a computer-simulation study. The Journal of Thoracic and Cardiovascular Surgery. 145 (5), 1352-1358 (2013).
- Fresiello, L., Meyns, B., Di Molfetta, A., Ferrari, G. A Model of the Cardiorespiratory Response to Aerobic Exercise in Healthy and Heart Failure Conditions. Frontiers in Physiology. 7 (189), (2016).
- Moscato, F., et al. Left ventricle afterload impedance control by an axial flow ventricular assist device: a potential tool for ventricular recovery. Artificial Organs. 34 (9), 736-744 (2010).
- Colacino, F. M., Moscato, F., Piedimonte, F., Arabia, M., Danieli, G. A. Left ventricle load impedance control by apical VAD can help heart recovery and patient perfusion: a numerical study. Asaio Journal. 53 (3), 263-277 (2007).
- Gu, K., et al. Lumped parameter model for heart failure with novel regulating mechanisms of peripheral resistance and vascular compliance. Asaio Journal. 58 (3), 223-231 (2012).
- Suga, H., Sagawa, K., Kostiuk, D. P. Controls of ventricular contractility assessed by pressure-volume ratio, Emax. Cardiovascular Research. 10 (5), 582-592 (1976).
- Fernandez de Canete, J., Saz-Orozco, P. d, Moreno-Boza, D., Duran-Venegas, E. Object-oriented modeling and simulation of the closed loop cardiovascular system by using SIMSCAPE. Computers in Biology and Medicine. 43 (4), 323-333 (2013).
- Heldt, T., Shim, E. B., Kamm, R. D., Mark, R. G., et al. Computational modeling of cardiovascular response to orthostatic stress. Journal of Applied Physiology. 92 (3), 1239-1254 (2002).
- Granegger, M., et al. A Valveless Pulsatile Pump for the Treatment of Heart Failure with Preserved Ejection Fraction: A Simulation Study. Cardiovascular Engineering and Technology. 10 (1), 69-79 (2019).
- Hay, I., Rich, J., Ferber, P., Burkhoff, D., Maurer, M. S. Role of impaired myocardial relaxation in the production of elevated left ventricular filling pressure. American Journal of Physiology-Heart and Circulatory Physiology. 288 (3), 1203-1208 (2005).
- Kadry, K., et al. Biomechanics of diastolic dysfunction: a one-dimensional computational modeling approach. American Journal of Physiology-Heart and Circulatory Physiology. 319 (4), 882-892 (2020).
- Luo, C., Ramachandran, D., Ware, D. L., Ma, T. S., Clark, J. W. Modeling left ventricular diastolic dysfunction: classification and key indicators. Theoretical Biology & Medical Modelling. 8, 14 (2011).
- Burkhoff, D., et al. Left atrial decompression pump for severe heart failure with preserved ejection fraction: theoretical and clinical considerations. JACC: Heart Failure. 3 (4), 275-282 (2015).
- Ahmad Bakir, A., Al Abed, A., Stevens, M. C., Lovell, N. H., Dokos, S. A Multiphysics Biventricular Cardiac Model: Simulations With a Left-Ventricular Assist Device. Frontiers in Physiology. 9 (1259), (2018).
- Genet, M., Lee, L. C., Baillargeon, B., Guccione, J. M., Kuhl, E. Modeling pathologies of diastolic and systolic heart failure. Annals of Biomedical Engineering. 44 (1), 112-127 (2016).
- Sack, K. L., et al. Investigating the Role of Interventricular Interdependence in Development of Right Heart Dysfunction During LVAD Support: A Patient-Specific Methods-Based Approach. Frontiers in Physiology. 9 (520), (2018).
- Baillargeon, B., et al. Human cardiac function simulator for the optimal design of a novel annuloplasty ring with a sub-valvular element for correction of ischemic mitral regurgitation. Cardiovascular Engineering and Technology. 6 (2), 105-116 (2015).
- Sack, K. L., et al. Partial LVAD Restores Ventricular Outputs and Normalizes LV but not RV Stress Distributions in the Acutely Failing Heart in Silico. The International Journal of Artificial Organs. 39 (8), 421-430 (2016).
- Baumgartner, H., et al. Echocardiographic assessment of valve stenosis: EAE/ASE recommendations for clinical practice. Journal of the American Society of Echocardiography. 22 (1), 1-23 (2009).
- Rajani, R., Hancock, J., Chambers, J. The art of assessing aortic stenosis. Heart. 98, 14 (2012).
- Vahanian, A., et al. Guidelines on the management of valvular heart disease: The Task Force on the Management of Valvular Heart Disease of the European Society of Cardiology. European Heart Journal. 28 (2), 230-268 (2007).
- Matiwala, S., Margulies, K. B. Mechanical approaches to alter remodeling. Current Heart Failure Reports. 1 (1), 14-18 (2004).
- NIH Clinical Trials Registry. ImCardia for DHF to Treat Diastolic Heart Failure (DHF) Patient a Pilot Study (ImCardia). , (2011).