A simple phenomenological approach for myocardial contraction: formulation, parameter sensitivity study and applications in organ level simulations

Abstract

Contraction in myocardial tissue is the result of a complex process through which chemical energy on the cellular level is converted into the mechanical energy needed to circulate blood throughout the body. Due to its vital role for the organism, myocardial contractility is one of the most intensively investigated subjects in medical research. In this contribution, we suggest a novel phenomenological approach for myocardial contraction that is capable of producing realistic intracellular calcium concentration (ICC) and myocyte shortening graphs, can be easily calibrated to capture different ICC and contraction characteristics and, at the same time, is straightforward to implement and ensures efficient computer simulations. This study is inspired by the fact that existing models for myocardial contractility either contain a number of complex equations and material parameters, which reduce their feasibility, or are very simple and cannot accurately mimic reality, which eventually influences the realm of computer simulations. The proposed model in this manuscript considers first the evolution of the ICC through a logarithmic-type ordinary differential equation (ODE) having the normalized transmembrane potential as the argument. The ICC is further put into an exponential-type ODE which determines the shortening of the myocyte (active stretch). The developed approach can be incorporated with phenomenological or biophysically based models of cardiac electrophysiology. Through examples on the material level, we demonstrate that the shape of the ICC and myocardial shortening curves can be easily modified and accurately fitted to experimental data obtained from rat and mouse hearts. Moreover, the performance of the model in organ level simulations is illustrated through several multi-field initial-boundary value problems in which we show variations in volume-time relations, heterogeneous characteristics of myocardial contraction and application of a drug in a virtual left ventricle model.