Multiple Dynamical Mechanisms of Phase-2 Early Afterdepolarizations in a Human Ventricular Myocyte Model: Involvement of Spontaneous SR Ca 2+ Release

Yasutaka Kurata; Kunichika Tsumoto; Kenshi Hayashi; Ichiro Hisatome; Yuhichi Kuda; Mamoru Tanida

doi:10.3389/fphys.2019.01545

Multiple Dynamical Mechanisms of Phase-2 Early Afterdepolarizations in a Human Ventricular Myocyte Model: Involvement of Spontaneous SR Ca ²⁺ Release

Front Physiol. 2020 Jan 10:10:1545. doi: 10.3389/fphys.2019.01545. eCollection 2019.

Authors

Yasutaka Kurata¹, Kunichika Tsumoto¹, Kenshi Hayashi², Ichiro Hisatome³, Yuhichi Kuda¹, Mamoru Tanida¹

Affiliations

¹ Department of Physiology II, Kanazawa Medical University, Uchinada, Japan.
² Department of Cardiovascular and Internal Medicine, Graduate School of Medical Sciences, Kanazawa University, Kanazawa, Japan.
³ Department of Genetic Medicine and Regenerative Therapeutics, Graduate School of Medical Sciences, Tottori University, Yonago, Japan.

Abstract

Early afterdepolarization (EAD) is known to cause lethal ventricular arrhythmias in long QT syndrome (LQTS). In this study, dynamical mechanisms of EAD formation in human ventricular myocytes (HVMs) were investigated using the mathematical model developed by ten Tusscher and Panfilov (Am J Physiol Heart Circ Physiol 291, 2006). We explored how the rapid (I_Kr) and slow (I_Ks) components of delayed-rectifier K⁺ channel currents, L-type Ca²⁺ channel current (I_Ca _L), Na⁺/Ca²⁺ exchanger current (I_NCX), and intracellular Ca²⁺ handling via the sarcoplasmic reticulum (SR) contribute to initiation, termination and modulation of phase-2 EADs during pacing in relation to bifurcation phenomena in non-paced model cells. Parameter-dependent dynamical behaviors of the non-paced model cell were determined by calculating stabilities of equilibrium points (EPs) and limit cycles, and bifurcation points to construct bifurcation diagrams. Action potentials (APs) and EADs during pacing were reproduced by numerical simulations for constructing phase diagrams of the paced model cell dynamics. Results are summarized as follows: (1) A modified version of the ten Tusscher-Panfilov model with accelerated I_CaL inactivation could reproduce bradycardia-related EADs in LQTS type 2 and β-adrenergic stimulation-induced EADs in LQTS type 1. (2) Two types of EADs with different initiation mechanisms, I_CaL reactivation-dependent and spontaneous SR Ca²⁺ release-mediated EADs, were detected. (3) Termination of EADs (AP repolarization) during pacing depended on the slow activation of I_Ks. (4) Spontaneous SR Ca²⁺ releases occurred at higher Ca²⁺ uptake rates, attributable to the instability of steady-state intracellular Ca²⁺ concentrations. Dynamical mechanisms of EAD formation and termination in the paced model cell are closely related to stability changes (bifurcations) in dynamical behaviors of the non-paced model cell, but they are model-dependent. Nevertheless, the modified ten Tusscher-Panfilov model would be useful for systematically investigating possible dynamical mechanisms of EAD-related arrhythmias in LQTS.

Keywords: bifurcation analysis; early afterdepolarization; long QT syndrome; mathematical model; spontaneous SR Ca2+ release.