Abstract
We developed new numerical methods to optimally adjust the parameters in cardiac electrophysiology models, using optimal control and non-differentiable optimization methods. We define an optimal control problem to adjust parameters in single-cell models so that the trans-membrane potential predicted by a model fits in a least-square (LS) sense the potential recorded over time. To account for restitution properties, this LS function measures the discrepancy between predictions and experiments for a cell paced at various heart rates (HR) of increasing frequency. The methodology is used to adjust parameters in the Mitchell-Schaeffer model to unscaled non-smoothed experimental recording of the trans-membrane potential obtained in pig heart using optical fluorescence imaging based on voltage-sensitive dye, and simultaneously identify scaling factors for the experimental data. The methodology is validated by adjusting the model for multiple heart beats at a single HR. The fit for a single HR is excellent (LS function = 0.0065–0.02). The methodology is applied to adjust the MS model to multiple heart beats at three different HR. It is observed that the fit remains good when the range of HR is moderately large (LS function = 0.052), while a larger HR gap is more challenging (LS function = 0.17).
Supported by the National Science and Engineering Council of Canada.
This is a preview of subscription content, log in via an institution to check access.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
References
Chong, E.K.P., Zak, S.H.: An Introduction to Optimization, 3rd edn. Wiley, Hoboken (2008)
W. Groenendaal et al.: Cell-specific cardiac electrophysiology models. PLOS Comput. Biol. 22 pages (2015)
Kaur, J., Nygren, A., Vigmond, E.J.: Fitting membrane resistance along with action potential shape in cardiac myocytes improves convergence: application of a multi-objective parallel genetic algorithm. PLOS One 9(9), 1–10 (2014)
Kolda, T.G., Lewis, R.M., Torczon, V.: Optimization by direct search: new perspectives on some classical and modern methods. SIAM Rev. 45(3), 385–482 (2003)
Lombardo, D.M., Fenton, F.H., Narayan, S.M., Rappel, W.-J.: Comparison of detailed and simplified models of human atrial myocytes to recapitulate patient specific properties. PLOS Comput. Biol. 15 pages (2016)
Mitchell, C.C., Schaeffer, D.G.: A two-current model for the dynamics of cardiac membrane. Bull. Math. Biol. 65(5), 767–793 (2003)
Pongui Ngoma, D.V., Bourgault, Y., Pop, M., Nkounkou, H.: Adjustment of parameters in ionic models using optimal control problems. In: Pop, M., Wright, G.A. (eds.) FIMH 2017. LNCS, vol. 10263, pp. 322–332. Springer, Cham (2017). https://doi.org/10.1007/978-3-319-59448-4_31
Pongui-Ngoma, D.V., Bourgault, Y., Nkounkou, H.: Parameter identification for a non-differentiable ionic model used in cardiac electrophysiology. Appl. Math. Sci. 9(150), 7483–7507 (2015)
Pop, M., et al.: Fusion of optical imaging and MRI for the evaluation and adjustment of macroscopic models of cardiac electrophysiology: a feasibility study. Med Image Anal. 13(2), 370–80 (2009)
Relan, J., Pop, M., Delingette, H., Wright, G., Ayache, N., Sermesant, M.: Personalization of a cardiac electrophysiology model using optical mapping and MRI for prediction of changes with pacing. IEEE Trans. Biomed. Eng. 10(10), 11 pages (2011)
Relan, J., et al.: Coupled personalization of cardiac electrophysiology models for prediction of ischaemic ventricular tachycardia. Interface Focus 1, 396–407 (2011)
Rioux, M., Bourgault, Y.: A predictive method allowing the use of a single ionic model in numerical cardiac electrophysiology. ESAIM: Math. Modell. Numer. Anal. 47, 987–1016 (2013)
Syed, Z., Vigmond, E., Nattel, S., Leon, L.J.: Atrial cell action potential parameter fitting using genetic algorithms. Med. Biol. Eng. Comput. 43, 561–571 (2005)
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2019 Springer Nature Switzerland AG
About this paper
Cite this paper
Pearce-Lance, J., Pop, M., Bourgault, Y. (2019). A Numerical Method for the Optimal Adjustment of Parameters in Ionic Models Accounting for Restitution Properties. In: Coudière, Y., Ozenne, V., Vigmond, E., Zemzemi, N. (eds) Functional Imaging and Modeling of the Heart. FIMH 2019. Lecture Notes in Computer Science(), vol 11504. Springer, Cham. https://doi.org/10.1007/978-3-030-21949-9_6
Download citation
DOI: https://doi.org/10.1007/978-3-030-21949-9_6
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-030-21948-2
Online ISBN: 978-3-030-21949-9
eBook Packages: Computer ScienceComputer Science (R0)