×

A condition for setting off ectopic waves in computational models of excitable cells. (English) Zbl 1145.92019

Summary: The purpose of this paper is to study the stability of steady state solutions of the monodomain model equipped with C. H. Luo and Y. Rudy [Circ. Res. 68, 1501 ff (1991)] I kinetics. It is well established that re-entrant arrhythmias can be created in computational models of excitable cells. Such arrhythmias can be initiated by applying an external stimulus that interacts with a partially refractory region, and spawn breaking waves that can eventually generate extremely complex wave patterns commonly referred to as fibrillation. An ectopic wave is one possible stimulus that may initiate fibrillation. Physiologically, it is well known that ectopic waves exist, but the mechanism for initiating ectopic waves in a large collection of cells is poorly understood. We consider computational models of collections of excitable cells in one and two spatial dimensions. The cells are modeled by Luo-Rudy I kinetics, and we assume that the spatial dynamics is governed by the monodomain model. The mathematical analysis is carried out for a reduced model that is known to provide good approximations of the initial phase of solutions of the Luo-Rudy I model. A further simplification is also introduced to motivate and explain the results for the more complicated models.
In the analysis the cells are divided into two regions; one region (N) consists of normal cells as modelled by the standard Luo-Rudy I model, and another region (A) where the cells are automatic in the sense that they would act as pacemaker cells if they where isolated from their surroundings. We let \(\delta \) denote the spatial diffusion and \(a\) denote a characteristic length of the automatic region. It has previously been shown that reducing diffusion or increasing the automatic region enhances ectopic activity. Here we derive a condition for the transition from stable resting state to ectopic wave spread. Under suitable assumptions on the model we provide mathematical and computational arguments indicating that there is a constant \(\eta \) such that a steady state solution of this system is stable whenever \(\delta \gtrsim \eta a^{2}\), and unstable whenever \(\delta \lesssim \eta a^{2}\).

MSC:

92C50 Medical applications (general)
92C37 Cell biology
92-08 Computational methods for problems pertaining to biology
37N25 Dynamical systems in biology
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Luo, C. H.; Rudy, Y., A model of the ventricular cardiac action potential: depolarisation, repolarisation, and their interaction, Circ. Res., 68, 1501 (1991)
[2] Weiss, J. N.; Karma, A.; Shiferaw, Y.; Chen, P.-S.; Garfinkel, A.; Qu, Z., From pulsus to pulseless: the saga of cardiac alternans, Circ. Res., 98, 1244 (2006)
[3] Weiss, J. N.; Qu, Z.; Chen, P.-S.; Lin, S.-F.; Karagueuzian, H. S.; Hayashi, H.; Garfinkel, A.; Karma, A., The dynamics of cardiac fibrillation, Circulation, 112, 1232 (2005)
[4] D.X. Tran, M.-J. Yang, J.N. Weiss, A. Garfinkel, Z. Qua, Vulnerability to re-entry in simulated two-dimensional cardiac tissue: effects of electrical restitution and stimulation sequence, Chaos 17 (2007) 043115-1-11.; D.X. Tran, M.-J. Yang, J.N. Weiss, A. Garfinkel, Z. Qua, Vulnerability to re-entry in simulated two-dimensional cardiac tissue: effects of electrical restitution and stimulation sequence, Chaos 17 (2007) 043115-1-11. · Zbl 1163.37377
[5] Keener, J. P., Model for the onset of fibrillation following coronary artery occlusion, J. Cardiovasc. Electrophysiol., 14, 11, 1225 (2003)
[6] Pumir, A.; Arutunyan, A.; Krinsky, V.; Sarvazyan, N., Genesis of ectopic waves: role of coupling, automaticity, and heterogeneity, Biophys. J., 89, 2332 (2005)
[7] Haïssaguerre, M.; Jaïs, P.; Shah, D. C.; Takahashi, A.; Hocini, M.; Quiniou, G.; Garrigue, S.; Le Mouroux, A.; Le Mtayer, P.; Clémenty, J., Spontaneous initiation of atrial fibrillation by ectopic beats originating in the pulmonary veins, N. Engl. J. Med., 339, 659 (1998)
[8] Hayashi, H.; Wang, C.; Miyauchi, Y.; Omichi, C.; Pak, H.-N.; Zhou, S.; Ohara, T.; Mandel, W. J.; Lin, S.-F.; Fishbein, M. C.; Chen, P.-S.; Karagueuzian, H. S., Aging-related increase to inducible atrial fibrillation in the rat model, J. Cardiovasc. Electrophysiol., 13, 8, 801 (2002)
[9] Joyner, R. W.; Wilders, R.; Wagner, M. B., Propagation of pacemaker activity, Med. Biol. Eng. Comput., 45, 2, 177 (2006)
[10] K.H.W.J. ten Tusscher, A.V. Panfilov, Influence of nonexcitable cells on spiral breakup in two-dimensional and three-dimensional excitable media, Phys. Rev. E 68 (2003) 062902-1-4.; K.H.W.J. ten Tusscher, A.V. Panfilov, Influence of nonexcitable cells on spiral breakup in two-dimensional and three-dimensional excitable media, Phys. Rev. E 68 (2003) 062902-1-4.
[11] V. Jacquemet, Pacemaker activity resulting from the coupling with nonexcitable cells, Phys. Rev. E 74 (2006) 011908-1-7.; V. Jacquemet, Pacemaker activity resulting from the coupling with nonexcitable cells, Phys. Rev. E 74 (2006) 011908-1-7.
[12] Aliev, R. R.; Panfilov, A. V., A simple two-variable model of cardiac excitation, Chaos, Solitons & Fractals, 7, 3, 293 (1996)
[13] Khan, R., Identifying and understanding the role of pulmonary vein activity in atrial fibrillation, Cardiovasc. Res., 64, 387 (2004)
[14] Jacquemet, V., Steady-state solutions in mathematical models of atrial cell electrophysiology and their stability, Math. Biosci., 208, 1, 241 (2007) · Zbl 1117.92018
[15] CellML Model Repository. Available from: <www.cellml.org/models/>; CellML Model Repository. Available from: <www.cellml.org/models/>
[16] ten Tusscher, K. H.W. J.; Bernus, O.; Hren, R.; Panfilov, V., Comparison of electrophysiological models for human ventricular cells and tissues, Prog. Biophys. Mol. Biol., 90, 1-3, 326 (2006)
[17] Hinch, R., An analytical study of the physiology and pathology of the propagation of cardiac action potentials, Prog. Biophys. Mol. Biol., 78, 1, 45 (2002)
[18] S.O. Linge, G.T. Lines, J. Sundnes, A. Tveito, On the frequency of automaticity during ischemia in simulations based on stochastic perturbations of the Luo-Rudy 1 model, Comput. Biol. Med. (2008), submitted for publication.; S.O. Linge, G.T. Lines, J. Sundnes, A. Tveito, On the frequency of automaticity during ischemia in simulations based on stochastic perturbations of the Luo-Rudy 1 model, Comput. Biol. Med. (2008), submitted for publication.
[19] MacLachlan, M. C.; Sundnes, J.; Skavhaug, O.; Lysaker, O. M.; Nielsen, B. F.; Tveito, A., A linear system of partial differential equations modeling the resting potential of a heart with regional ischemia, Math. Biosci., 210, 1, 238 (2007) · Zbl 1129.92043
[20] Logan, J. D., An Introduction to Nonlinear Partial Differential Equations (1994), Wiely Interscience: Wiely Interscience New York · Zbl 0834.35001
[21] Tveito, A.; Winther, R., Introduction to Partial Differential Equation, Volume 29 of Texts in Applied Mathematics (1998), Springer: Springer New York
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.