Intracellular calcium (Ca i
IntroductionIntracellular calcium handling is a key process in the electrophysiology of cardiac myocytes [1,2]. On the one hand, it determines cell contractility [1], because intracellular calcium is the main regulator of myocardial excitation-contraction coupling. On the other hand, it is intimately related to arrhythmogenesis, because calcium overload may trigger abnormal cell depolarization [2,3]. The dynamic changes in intracellular calcium concentration ([Ca 2+ ] i ) during the action potential are complex. Upon cell depolarization at the onset of an action potential (AP), calcium enters the cell through the L-type calcium channels. This triggers a massive release of calcium from the sarcoplasmic reticulum (SR) to the cytosol through the ryanodine receptors, a process that may elevate [Ca 2+ ] i almost an order of magnitude compared to its basal value (the diastolic calcium level). After [Ca 2+ ] i peaks (reaching its systolic level), it slowly declines due to (a) extrusion to the extracellular medium via the sodiumcalcium exchanger (NCX) in the sarcolemma, and (b) uptake to the SR via the SERCA pump [4].In the past two decades, mathematical modelling and computational simulations have been a powerful tool used to better understand the intricate mechanisms related to cellular electrophysiology of cardiac myocytes (see [5] for review). These computational studies rely on detailed and comprehensive mathematical models of action potentials for different animal species [6][7][8], including human [9]. Using this models, different studies have attempted to characterize the influence of the different sarcolemmal and SR ionic currents in AP and intracellular Ca 2+ preclinical biomarkers [10,11], but none of them has been carried out in dog, an animal that has widely been used in the past for drug and mutation effect studies [12,13].In this study, we have used a well-known existing model of the action potential and ionic currents of a dog ventricular myocyte [7] to perform a sensitivity analysis to elucidate the ionic mechanisms that affect intracellular Ca 2+ regulation.
MethodsThe action potential, underlying ionic currents and dynamic changes in ionic concentrations of an isolated canine ventricular myocyte were simulated using a modified version of the Decker-Rudy model [7,13], which includes a detailed formulation of calcium dynamics. A sensitivity analysis was performed by applying a one-at-atime variation in the ionic currents through the sarcolemma and the SR membrane. Specifically, each current was multiplied by 0 (complete block), 0.3, 0.7, 1 (no change) and 2 in order to study their impact on selected calciumrelated biomarkers. These changes could represent the effect of drugs or mutations, or simply inter-individual physiological variability [14]. Simulations were performed at three different pacing rates (1 Hz, 2 Hz and 3 Hz) and ran until achieving steady-state. The measured calcium