[論文レビュー] A Map-Based Model of the Cardiac Action Potential

Nikolai F. Rulkov|ArXiv.org|Aug 8, 2007

Functional Brain Connectivity Studies被引用数 3

ひとこと要約

本論文は、心筋細胞のダイナミクスの主要な段階（急速な去分極、プラトー、再分極）を再現する2次元の離散時間写像を用いた、計算的に効率的なマップベースの心電位モデルを提案する。高速なスパiking写像と、制御パラメータによって調整される遅い回復変数を組み合わせることで、従来のODEモデルよりもはるかに大きな時間ステップを用いた大規模なシミュレーションが可能となり、安定な波動伝播および周波数ロックや波動分裂といった複雑なダイナミクスを示す結合ネットワークを再現する。

ABSTRACT

A discrete time model that is capable of replicating the basic features of cardiac cell action potentials is suggested. The paper shows how the map-based approaches can be used to design highly efficient computational models (algorithms) that enable large-scale simulations and analysis of discrete network models of cardiac activity.

研究の動機と目的

• 従来のODEベースのモデルが抱える数値的制約を回避する計算的に効率的な心電位モデルの開発。
• 非パーカイマーカード心筋細胞の5段階の電位ダイナミクス（急速な去分極、プラトー段階など）をマップベースのダイナミクスで再現すること。
• 離散時間写像を用いた大規模な電気的結合心筋ネットワークのシミュレーションを可能にする。
• 結合系において、波動伝播、周波数ロック、波動分裂といった顕在的ネットワーク行動を示す。

提案手法

• モデルは2次元の離散時間写像を用いる：$ x_{n+1} = P(x_n, y_n) $, $ y_{n+1} = Q(x_n, y_n) $、ここで$ x_n $は高速な去分極を表し、$ y_n $は回復および再分極を表す。
• $ x $-ダイナミクスは非線形関数$ f(x_n, x_{n-1}, u) $によって支配され、$ \varepsilon(y_n) = 0 $のときスパイク様の挙動を示し、$ \varepsilon(y_n) \to 1 $にすることでNa+チャネルの不活性化を模倣する。
• $ y $-ダイナミクスは制御パラメータ$ \mu $を有する1次元写像$ Q(x_n, y_n) $でモデル化され、プラトーおよび再分極段階を可能にする。
• ギャップジャンクションを通じた結合メカニズムとして、$ y_{n+1} = Q(x_n, y_n) + \mu_{\text{gap}} \sum_{j \in J} I_{n,j,i}^{\text{gap}} $を導入し、ネットワーク内での波動伝播を可能にする。
• 状態変数を膜電位にマッピングするために線形結合$ V_n = a x_n + b y_n $を用いる。
• モデルは500個の結合細胞からなる1次元リングで検証され、安定な波動伝播、周波数比の切り替え（1:2から1:1）、障害を示す細胞による波動分裂を示した。

実験結果

リサーチクエスチョン

• RQ1 マップベースの離散時間モデルは、高い計算効率を実現しつつ、心電位の主要な段階を再現できるか？
• RQ2 マップベースの手法は、ダイナミクスの忠実性を損なわずに、心電活動のシミュレーションにおいて大きな時間ステップを可能にするか？
• RQ3 結合マップベース心筋モデルにおいて、波動伝播や周波数ロックといった顕在的ネットワーク行動を捉えることができるか？
• RQ4 局所的な細胞異常（例：活動電位持続時間が短縮）は、結合ネットワーク内で波動分裂や新たな波源を引き起こすか？

主な発見

• モデルは5段階の活動電位（静止状態、急速な去分極（第0期）、初期再分極（第1期）、プラトー（第2期）、再分極（第3期））を正常に再現した。
• 制御可能なマップの使用により、従来のODEベースの心筋モデルで一般的に見られる剛性問題を克服し、大きな統合時間ステップを実現した。
• 500個の結合細胞からなるリングにおいて、1:1および1:2の周波数の安定な波動パターンを示し、外部のトリガーパルス1個によってそれらの間の切り替えが誘発された。
• 活動電位持続時間が短縮された30個の細胞（$ y_s = 1.3 $から$ y_s = 1.2 $に変更）は、摂動後に持続的なダブルト波動パターンを引き起こす新たな波源として機能した。
• モデルは波動ダイナミクスにおける記憶効果を捉えており、外部刺激のタイミングが1:2または1:1の周波数レジームに留まるか切り替わるかを決定する。
• ギャップジャンクション結合モデル（$ \mu_{\text{gap}} = 0.0001 $、$ g_{\text{gap}} = 0.004 $）は、ネットワーク内での現実的な波動伝播および回復ダイナミクスを可能にした。