Numerical Simulation of Excitation-Contraction in Isolated Cardiomyocytes

  • José Alberto Rodrigues Department of Mathematics, Instituto Superior de Engenharia de Lisboa-ISEL, Portugal GIMOSM, ISEL, IPL Polytechnic institute of lisbon
  • Joana Oliveira Biomedical Engineering, Instituto Superior de Engenharia de Lisboa-ISEL, Portugal
Keywords: Biomedical Engineering, Ca2 Oscillations, Mechano-chemical Mathematical Model, Finite Element Methods.


We present and study a mathematical model to simulate the calcium flux and contractile activity of the cardiomyocyte, resorting on the Finite Element Method. This Model of cardiac cells provide the possibility to understand thebiochemistry and biomechanics of cardiac cells and cardiovascular diseases. Heart Failure is considered the ultimate cardiac disease, a condition with no effective cure which is highly related with the function of cardiomyocytes and consequently with the concentration of calcium, affecting the contractile activity of the cardiomyocyte.In order to test our model’s performance, several tests were applied varying the local active cellular tension driven by the intracellular calcium concentration and the localization of the main calcium influx. The results are expressed in the graphics of calcium concentration over time, maximum cardiomyocyte contraction and the gradient of calcium diffusion.Our results show that the behaviour of our models is faithful to what is known to be true with cell’s physiology andpathological conditions.


U. Primessnig, P. Schönleitner, and A. Höll, Novel pathomechanis of cardiomyocyte dysfunction in a model of heart failure with preserved ejection fraction, European Journal of Heart Failure (2016) 18, 987–997.

P. Gorski2015, D. Ceholski, and R. Hajjar Altered Myocardial Calcium Cycling and Energetics in Heart Failure — A Rational Approach for Disease Treatment, Cell Metabolism 21,pp. 183–194, 2015.

N. Sperelakis, Y. Kurachi, A. Terzic, and M. Cohen, Heart Physiology and Pathophysiology - 4th Edition, Academic Press, 2001.

K. Barrett, S. Barman, S. Boitano, and H. Brooks, Ganong’s Review of Medical Physiology, 25th Ed. McGraw-Hill Education

D. Eisner, J. Cadwell, K. Kistamás and W. Trafford, Calcium and Excitation-Contraction Coupling in the Heart, Circulation Research, 121(2):pp. 181–195, 2017.

M. Jafri, Models of Excitation–Contraction Coupling in Cardiac Ventricular Myocytes, Methods Mol Biol. 2012 ; 910: 309–335. doi:10.1007/978 − 1 − 61779 − 965 − 514.

R. Klaubunde, Cardiac Excitation-Contraction Coupling, Cardiovascular Physiology Concepts -

L. Vilac¸a, Mathematical Modeling and Numerical Simulation of Excitation-Contraction Phenomena in the Heart, Departement of Chemistry in EPFL, 2013

G. Ikonnikov, and E. Wong, Action potentials and impulse conduction,

N. El-Sherif, and G. Turitto, Electrolyte disorders and arrhythmogenesis, Cardiology Journal, Vol. 18, No. 3, pp. 233–245, 2011.

A. Goldbeter, G. Dupont, and M. Berridge, Minimal model for signal-induced Ca2+ oscillations and for their frequency encoding through protein phosphorylation., Proceedings of the National Academy

of Sciences of the United States of America, 1990.

A. Beck, and M. Teboulle, A fast iterative shrinkage-thresholding algorithm for linear inverse problems, SIAM Journal on Imaging Sciences, vol. 2, no. 1, pp. 183–202, 2009.

S. Subramanian, S. Viatchenko-Karpinski, V. Lukyanenko,S Györk and T.F. Wiesner, Underlying mechanisms of symmetric calcium wave propagation in rat ventricular myocytes, Biophysical Journal,80(1), pp. 1–11, 2001.

O. Cazorla, A. Lacampagne, J. Fauconnier, and G. Vassort, , SR 33805 a Ca2+antagonist with lengtdependent Ca2+-sensitizing properties in cardiacmyocytes, Brit. J. Pharmacol., 139, pp 99–108, 2003.

Joseph D. Bronzino, eds., Biomedical Engineering Handbook Vol 1, Springer-Verlag Berlin and Heidelberg GmbH & Co. K ,2000.

C. Johnson , Numerical Solution of Partial Differential Equations by the Finite Element Method, Springer-Verlag Berlin and Heidelberg GmbH & Co. K ,2000.

, F. Hecht, New development in FreeFem++, Journal of Numerical Mathematics, 20 (3-4), pp. = 251–265, 2012.

, F. del Monte, and S. Harding, and U. Schmidt, and T. Matsui, and et al, Restoration of Contractile Function in Isolated Cardiomyocytes From Failing Human Hearts by Gene Transfer of SERCA2a, Circulation.,23 (100), pp 2308–2311, 1999.

, M. Høydal, and I. Kirkeby-Garstad, and A. Karevold, and R. Wiseth, and et al., Human cardiomyocyte calcium handling and transverse tubules in mid-stage of post-myocardial-infarction heart failure, ESC Hear. Fail., February, pp 332–342, 2018, doi = 10.1002/ehf2.12271.

, A. Kamgoué, and J. Ohayon, and Y. Usson, and L. Riou, and P. Tracqui, Quantification of cardiomyocyte contraction based on image correlation analysis, Cytom. Part A, 4(75), pp 298–308,2009, doi = 10.1002/cyto.a.20700.

, P. Garcia-Canadilla, and J. Rodriguez, and MJ. Palazzi, and A. Gonzalez-Tendero, and et al., A two dimensional electromechanical model of a cardiomyocyte to assess Intra-cellular regional mechanical heterogeneities PLoS One, 8(12), pp 1–20, 2017, doi = 10.1371/journal.pone.0182915.

How to Cite
Rodrigues, J. A., & Oliveira, J. (2019). Numerical Simulation of Excitation-Contraction in Isolated Cardiomyocytes. Statistics, Optimization & Information Computing, 7(3), 604-616.
Research Articles