Abstract
The fast hybrid operator splitting (HOS) and stable uniformization (UNI) methods have been proposed to save computation cost and enhance stability for Markov chain model in cardiac cell simulations. Moreover, Chen-Chen-Luo’s quadratic adaptive algorithm (CCL) combined with HOS or UNI was used to improve the tradeoff between speedup and stability, but without considering accuracy. To compromise among stability, acceleration, and accuracy, we propose a generalized Trotter operator splitting (GTOS) method combined with CCL independent of the asymptotic property of a particular ion-channel model. Due to the accuracy underestimation of the mixed root mean square error (MRMSE) method, threshold root mean square error (TRMSE) is proposed to evaluate computation accuracy. With the fixed time-step RK4 as a reference, the second-order GTOS combined with CCL (30.8-fold speedup) for the wild-type Markov chain model with nine states (WT-9 model) or (7.4-fold) for the wild-type Markov chain model with eight states (WT-8 model) is faster than UNI combined with CCL (15.6-fold) for WT-9 model or (1.2-fold) for WT-8 model, separately. Besides, the second-order GTOS combined with CCL has 3.81% TRMSE for WT-9 model or 4.32% TRMSE for WT-8 model more accurate than 72.43% TRMSE for WT-9 model or 136.17% TRMSE for WT-8 model of HOS combined with CCL. To compromise speedup and accuracy, low-order GTOS combined with CCL is suggested to have the advantages of high precision and low computation cost. For high-accuracy requirements, high-order GTOS combined with CCL is recommended.
Similar content being viewed by others
References
Lopez-Perez A, Sebastian R, Ferrero JM (2015) Three-dimensional cardiac computational modelling: methods, features and applications. Biomed Eng Online 14(1):35
Pathmanathan P, Gray RA (2018) Validation and trustworthiness of multiscale models of cardiac electrophysiology. Front Physiol 9:106
Fallahi A, Khorram HG, Kokabi A (2019) Electrocardiogram signal generation using electrical model of cardiac cell: application in cardiac ischemia. J Med Eng Technol 43(4):207–216
Gao Y, Xia L, Gong YL, Zheng DC (2018) Electrocardiogram (ECG) patterns of left anterior fascicular block and conduction impairment in ventricular myocardium: a whole-heart model-based simulation study. J Zhejiang Univ-Sci B 19(1):49–56
Prajapati C, Aalto-Setälä K (2019) Modelling of genetic cardiac diseases. In Visions of Cardiomyocyte-Fundamental Concepts of Heart Life and Disease. IntechOpen.
Rudy Y, Silva JR (2006) Computational biology in the study of cardiac ion channels and cell electrophysiology. Q Rev Biophys 39(1):57–116
Mangold KE, Brumback BD, Angsutararux P, Voelker TL, Zhu W, Kang PW, Silva JR (2017) Mechanisms and models of cardiac sodium channel inactivation. Channels 11(6):517–533
Silva JR, Pan H, Wu D, Nekouzadeh A, Decker KF, Cui J, Rudy Y (2009) A multiscale model linking ion-channel molecular dynamics and electrostatics to the cardiac action potential. Proc Natl Acad Sci 106(27):11102–11106
Carbonell-Pascual B, Godoy E, Ferrer A, Romero L, Ferrero JM (2016) Comparison between Hodgkin–Huxley and Markov formulations of cardiac ion channels. J Theor Biol 399:92–102
Asfaw TN, Bondarenko VE (2019) A mathematical model of the human cardiac Na+ channel. J Membr Biol 252(1):77–103
Campana C, Sobie EA (2019) A comparative analysis of parameter estimation strategies for mathematical modeling of ion channel gating. Biophys J 116(3):231a
Niederer SA, Lumens J, Trayanova NA (2019) Computational models in cardiology. Nat Rev Cardiol 16(2):100–111
Fink M, Noble D (2009) Markov models for ion channels: versatility versus identifiability and speed. Philos Trans R Soc A: Math Phys Eng Sci 367(1896):2161–2179
Greer-Short A, George SA, Poelzing S, Weinberg SH (2017) Revealing the concealed nature of long-QT type 3 syndrome. Circ Arrhythm Electrophysiol 10(2):e004400
Campana C, Gando I, Tan RB, Cecchin F, Coetzee WA, Sobie EA (2018) Population-based mathematical modeling to deduce disease-causing cardiac Na+ channel gating defects. Biophys J 114(3):634a–635a
Moreno JD, Lewis TJ, Clancy CE (2016) Parameterization for in-silico modeling of ion channel interactions with drugs. PLoS One 11(3):e0150761
Clancy CE, An G, Cannon WR, Liu Y, May EE, Ortoleva P, Zhou X (2016) Multiscale modeling in the clinic: drug design and development. Ann Biomed Eng 44(9):2591–2610
Tveito A, Maleckar MM, Lines GT (2018) Computing optimal properties of drugs using mathematical models of single channel dynamics. Comput Math Biophys 6(1):41–64
Spiteri RJ, Dean RC (2010) Stiffness analysis of cardiac electrophysiological models. Ann Biomed Eng 38(12):3592–3604
Sidje RB, Burrage K, MacNamara S (2007) Inexact uniformization method for computing transient distributions of Markov chains. SIAM J Sci Comput 29(6):2562–2580
Gomes JM, Alvarenga A, Campos RS, Rocha BM, da Silva APC, dos Santos RW (2014) Uniformization method for solving cardiac electrophysiology models based on the Markov-chain formulation. IEEE Trans Biomed Eng 62(2):600–608
Starý T, Biktashev VN (2014) Exponential integrators for a Markov chain model of the fast sodium channel of cardiomyocytes. IEEE Trans Biomed Eng 62(4):1070–1076
Korhonen T, Tavi P (2008) Automatic time-step adaptation of the forward Euler method in simulation of models of ion channels and excitable cells and tissue. Simul Model Pract Theory 16(6):639–644
Campos RS, Lobosco M, dos Santos RW (2011) Adaptive time step for cardiac myocyte models. Procedia Comput Sci 4:1092–1100
Gomes JM, dos Santos RW (2015) A time adaptive scheme based on the Rush-Larsen method for solving cardiac cell models. In VI Latin American Congress on Biomedical Engineering CLAIB 2014, Paraná, Argentina 29, 30 & 31 October 2014 (pp. 896-899)
García-Molla VM, Liberos A, Vidal A, Guillem MS, Millet J, Gonzalez A, Climent AM (2014) Adaptive step ODE algorithms for the 3D simulation of electric heart activity with graphics processing units. Comput Biol Med 44:15–26
Chen MH, Chen PY, Luo CH (2016) Quadratic adaptive algorithm for solving cardiac action potential models. Comput Biol Med 77:261–273
Chen XJ, Luo CH, Chen MH, Zhou X (2019) Combination of “quadratic adaptive algorithm” and “hybrid operator splitting” or uniformization algorithms for stability against acceleration in the Markov model of sodium ion channels in the ventricular cell model. Med Biol Eng Comput 57(6):1367–1379
Marsh ME, Ziaratgahi ST, Spiteri RJ (2012) The secrets to the success of the Rush–Larsen method and its generalizations. IEEE Trans Biomed Eng 59(9):2506–2515
Suzuki M (1985) Decomposition formulas of exponential operators and Lie exponentials with some applications to quantum mechanics and statistical physics. J Math Phys 26(4):601–612
Suzuki M (1992) Fractal path integrals with applications to quantum many-body systems. Phys A: Stat Mech Appl 191(1–4):501–515
Sundnes J, Artebrant R, Skavhaug O, Tveito A (2009) A second-order algorithm for solving dynamic cell membrane equations. IEEE Trans Biomed Eng 56(10):2546–2548
Spiteri RJ, Dean RC (2008) On the performance of an implicit–explicit Runge--Kutta method in models of cardiac electrical activity. IEEE Trans Biomed Eng 55(5):1488–1495
Clancy CE, Rudy Y (2002) Na+ channel mutation that causes both Brugada and long-QT syndrome phenotypes. Circulation 105:1208–1213
Luo CH, Rudy Y (1994) A dynamic model of the cardiac ventricular action potential. I Simulations of ionic currents and concentration changes. Circ Res 74(6):1071–1096
Hille B (2001) Ion channels of excitable membranes, vol 507. Sinauer, Sunderland
Suzuki M (1992) General nonsymmetric higher-order decomposition of exponential operators and symplectic integrators. J Phys Soc Jpn 61:3015–3019
Funding
This work was supported by Sun Yat-sen University, China, under Scientific Initiation Project (No.67000-18821109) for High-level Experts. M.-H. Chen was supported by the Ministry of Science and Technology of Taiwan under grant MOST 108-2115-M-006-019-.
Author information
Authors and Affiliations
Corresponding authors
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Xing-Ji Chen and Ching-Hsing Luo are first author
Electronic supplementary material
ESM 1
(PDF 305 kb)
Rights and permissions
About this article
Cite this article
Chen, XJ., Luo, CH. & Chen, MH. Combination of “generalized Trotter operator splitting” and “quadratic adaptive algorithm” method for tradeoff among speedup, stability, and accuracy in the Markov chain model of sodium ion channels in the ventricular cell model. Med Biol Eng Comput 58, 2131–2141 (2020). https://doi.org/10.1007/s11517-020-02220-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11517-020-02220-x