MioLab: Simulator for cardiac myocyte contractile force of rat based on the dynamics of calcium

Abstract

Background: The calcium ion (Ca2+) is essential in cardiac electrical activity since it is the direct activator of myofilaments, which cause contraction. For this reason, irregular Ca2+ flow in the cardiac myocyte is one of the main causes of cardiac arrhythmias and contractile dysfunction. In this sense, the formulation of mathematical and computational models of the mammalian ventricular myocyte has played an important role in understanding cardiac physiology. This paper proposes a biophysical model for the rat cardiac myocyte in order to reduce the number of degrees of freedom and provides a mathematical framework that makes computation simpler for our understanding of the complex process of excitation-contraction-relaxation (ECR), uses low-order lumped parameter of the myofilaments and Ca2+ handling. The model allows the calculation of contractile force mainly based on the dynamics of Ca2+, provides data on transient Ca2+ and enables the analysis of the effects of drugs used in the treatment of cardiac arrhythmias. Methods: The mathematical model which describes the kinetics of cytosolic Ca2+ and dynamic cross-bridges, was implemented through nine differential equations, auxiliary equations, three 47 biophysical parameters. Each was implemented using C++ programming and its user-friendly interface in the programming language Delphi. The results of each simulation were stored in a text file and its analysis was shown in a graphical interface using executable programs in Matlab (R). Additionally, the validation model was performed by comparing both the experimental data of Ca2+ transient and several drug tests. Results: The model satisfactorily reproduced the Ca2+ transient, as well as the effects of drugs that cause beta-adrenergic stimulation and the inhibition of the Ca2+ uptake mechanism (SERCA). Changes in the parameters regulating Ca2+ entry through L-type channels produced the oscillation amplitude of the Ca2+ transient known as the syndrome of pulses alternans. Changes in parameters related to the Na+-Ca2+ exchange has already stabilized the transient Ca2+ and produced a decrease in amplitude of the contractile force. Conclusion: The simulator proved to be a tool to study and understand the mechanisms that involve the kinetics of Ca2+ and the dynamics of cross-bridges in the unit heart muscle, as well as a tool to analyze the possible factors that may cause arrhythmias and study the effects of drugs that are used in the treatment of cardiovascular diseases in general.