Integrative Computational Modeling of Cardiomyocyte Calcium Handling and Cardiac Arrhythmias: Current Status and Future Challenges

Cardiac arrhythmias in fish induced by natural and anthropogenic changes in environmental conditions

A regular heartbeat is essential for maintaining the homeostasis of the vertebrate body. However, environmental pollutants, oxygen deficiency and extreme temperatures can impair heart function in fish. In this Review, we provide an integrative view of the molecular origins of cardiac arrhythmias and their functional consequences, from the level of ion channels to cardiac electrical activity in living fish. First, we describe the current knowledge of the cardiac excitation-contraction coupling of fish, as the electrical activity of the heart and intracellular Ca2+ regulation act as a platform for cardiac arrhythmias. Then, we compile findings on cardiac arrhythmias in fish. Although fish can experience several types of cardiac arrhythmia under stressful conditions, the most typical arrhythmia in fish - both under heat stress and in the presence of toxic substances - is atrioventricular block, which is the inability of the action potential to progress from the atrium to the ventricle. Early and delayed afterdepolarizations are less common in fish hearts than in the hearts of endotherms, perhaps owing to the excitation-contraction coupling properties of the fish heart. In fish hearts, Ca2+-induced Ca2+ release from the sarcoplasmic reticulum plays a smaller role than Ca2+ influx through the sarcolemma. Environmental changes and ion channel toxins can induce arrhythmias in fish and weaken their tolerance to environmental stresses. Although different from endotherm hearts in many respects, fish hearts can serve as a translational model for studying human cardiac arrhythmias, especially for human neonates.

Feed Forward Modeling: an efficient approach for mathematical modeling of the force frequency relationship in the rabbit isolated ventricular myocyte

Background and Objective. This study addresses the Force-Frequency relationship, a fundamental characteristic of cardiac muscle influenced byβ1-adrenergic stimulation. This relationship reveals that heart rate (HR) changes at the sinoatrial node lead to alterations in ventricular cell contractility, increasing the force and decreasing relaxation time for higher beat rates. Traditional models lacking this relationship offer an incomplete physiological depiction, impacting the interpretation of in silico experiment results. To improve this, we propose a new mathematical model for ventricular myocytes, named 'Feed Forward Modeling' (FFM).Methods. FFM adjusts model parameters like channel conductance and Ca2+pump affinity according to stimulation frequency, in contrast to fixed parameter values. An empirical sigmoid curve guided the adaptation of each parameter, integrated into a rabbit ventricular cell electromechanical model. Model validation was achieved by comparing simulated data with experimental current-voltage (I-V) curves for L-type Calcium and slow Potassium currents.Results. FFM-enhanced simulations align more closely with physiological behaviors, accurately reflecting inotropic and lusitropic responses. For instance, action potential duration at 90% repolarization (APD90) decreased from 206 ms at 1 Hz to 173 ms at 4 Hz using FFM, contrary to the conventional model, where APD90 increased, limiting high-frequency heartbeats. Peak force also showed an increase with FFM, from 8.5 mN mm-2at 1 Hz to 11.9 mN mm-2at 4 Hz, while it barely changed without FFM. Relaxation time at 50% of maximum force (t50) similarly improved, dropping from 114 ms at 1 Hz to 75.9 ms at 4 Hz with FFM, a change not observed without the model.Conclusion. The FFM approach offers computational efficiency, bypassing the need to model all beta-adrenergic pathways, thus facilitating large-scale simulations. The study recommends that frequency change experiments include fractional dosing of isoproterenol to better replicate heart conditionsin vivo.

Electrophysiological effects of stretch-activated ion channels: a systematic computational characterization

Cardiac electrophysiology and mechanics are strongly interconnected. Their interaction is, among others, mediated by mechano-electric feedback through stretch-activated ion channels (SACs). The electrophysiological changes induced by SACs may contribute to arrhythmogenesis, but the precise SAC-induced electrophysiological changes remain incompletely understood. Here, we provide a systematic characterization of stretch effects through three distinguished SACs on cardiac electrophysiology using computational modelling. We implemented potassium-selective, calcium-selective and non-selective SACs in the Tomek-Rodriguez-O'Hara-Rudy model of human ventricular electrophysiology. The model was calibrated to experimental data from isolated cardiomyocytes undergoing stretch, considering inter-species differences, and disease-related remodelling of SACs. SAC-mediated effects on the action potential (AP) were analysed by varying stretch amplitude, application timing and/or duration. Afterdepolarizations of different amplitudes were observed with transient 10-ms stretch stimuli of 15-18% applied during phase 4, while stretch ≥18% during phase 4 elicited triggered APs. Longer stimuli shifted the threshold of AP trigger during phase 4 to lower amplitudes, while shorter stimuli increased it. Continuous stretch provoked electrophysiological remodelling. Furthermore, stretch shortened duration or changed morphology of a subsequent electrically evoked AP, and, if applied during a vulnerable time window with sufficient amplitude, prevented its occurrence because of stretch-induced modulation of sodium and L-type calcium channel gating. These effects were more pronounced with disease-related SAC remodelling due to increased stretch sensitivity of diseased hearts. We showed that SACs may induce afterdepolarizations and triggered activities, and prevent subsequent AP generation or change its morphology. These effects depend on cardiomyocyte stretch characteristics and disease-related SACs remodelling and may contribute to cardiac arrhythmogenesis. KEY POINTS: The interplay between cardiac electrophysiology and mechanics is mediated by mechano-electric feedback through stretch-activated ion channels (SACs). These channels may be pro-arrhythmic, but their precise effect on electrophysiology remains unclear. Here we present a systematic in silico characterization of stretch effects through three SACs by implementing inter-species differences as well as disease-related remodelling of SACs in a novel computational model of human ventricular cardiomyocyte electrophysiology. Our simulations showed that, at the cellular level, SACs may provoke electrophysiological remodelling, afterdepolarizations, triggered activities, change the morphology or shorten subsequent electrically evoked action potentials. The model further suggests that a vulnerable window exists in which stretch prevents the following electrically triggered beat occurrence. The pro-arrhythmic effects of stretch strongly depend on disease-related SAC remodelling as well as on stretch characteristics, such as amplitude, time of application and duration. Our study helps in understanding the role of stretch in cardiac arrhythmogenesis and revealing the underlying cellular mechanisms.

A detailed mathematical model of the human atrial cardiomyocyte: integration of electrophysiology and cardiomechanics

Mechano-electric regulations (MER) play an important role in the maintenance of cardiac performance. Mechano-calcium and mechano-electric feedback (MCF and MEF) pathways adjust the cardiomyocyte contractile force according to mechanical perturbations and affects electro-mechanical coupling. MER integrates all these regulations in one unit resulting in a complex phenomenon. Computational modelling is a useful tool to accelerate the mechanistic understanding of complex experimental phenomena. We have developed a novel model that integrates the MER loop for human atrial cardiomyocytes with proper consideration of feedforward and feedback pathways. The model couples a modified version of the action potential (AP) Koivumäki model with the contraction model by Quarteroni group. The model simulates iso-sarcometric and isometric twitches and the feedback effects on AP and Ca2+-handling. The model showed a biphasic response of Ca2+ transient (CaT) peak to increasing pacing rates and highlights the possible mechanisms involved. The model has shown a shift of the threshold for AP and CaT alternans from 4.6 to 4 Hz under post-operative atrial fibrillation, induced by depressed SERCA activity. The alternans incidence was dependent on a chain of mechanisms including RyRs availability time, MCF coupling, CaMKII phosphorylation, and the stretch levels. As a result, the model predicted a 10% slowdown of conduction velocity for a 20% stretch, suggesting a role of stretch in creation of substrate formation for atrial fibrillation. Overall, we conclude that the developed model provides a physiological CaT followed by a physiological twitch. This model can open pathways for the future studies of human atrial electromechanics. KEY POINTS: With the availability of human atrial cellular data, interest in atrial-specific model integration has been enhanced. We have developed a detailed mathematical model of human atrial cardiomyocytes including the mechano-electric regulatory loop. The model has gone through calibration and evaluation phases against a wide collection of available human in-vitro data. The usefulness of the model for analysing clinical problems has been preliminaryly tested by simulating the increased incidence of Ca2+ transient and action potential alternans at high rates in post-operative atrial fibrillation condition. The model determines the possible role of mechano-electric feedback in alternans incidence, which can increase vulnerability to atrial arrhythmias by varying stretch levels. We found that our physiologically accurate description of Ca2+ handling can reproduce many experimental phenomena and can help to gain insights into the underlying pathophysiological mechanisms.

Local Control Model of a Human Ventricular Myocyte: An Exploration of Frequency-Dependent Changes and Calcium Sparks

Calcium (Ca2+) sparks are the elementary events of excitation-contraction coupling, yet they are not explicitly represented in human ventricular myocyte models. A stochastic ventricular cardiomyocyte human model that adapts to intracellular Ca2+ ([Ca2+]i) dynamics, spark regulation, and frequency-dependent changes in the form of locally controlled Ca2+ release was developed. The 20,000 CRUs in this model are composed of 9 individual LCCs and 49 RyRs that function as couplons. The simulated action potential duration at 1 Hz steady-state pacing is ~0.280 s similar to human ventricular cell recordings. Rate-dependence experiments reveal that APD shortening mechanisms are largely contributed by the L-type calcium channel inactivation, RyR open fraction, and [Ca2+]myo concentrations. The dynamic slow-rapid-slow pacing protocol shows that RyR open probability during high pacing frequency (2.5 Hz) switches to an adapted "nonconducting" form of Ca2+-dependent transition state. The predicted force was also observed to be increased in high pacing, but the SR Ca2+ fractional release was lower due to the smaller difference between diastolic and systolic [Ca2+]SR. Restitution analysis through the S1S2 protocol and increased LCC Ca2+-dependent activation rate show that the duration of LCC opening helps modulate its effects on the APD restitution at different diastolic intervals. Ultimately, a longer duration of calcium sparks was observed in relation to the SR Ca2+ loading at high pacing rates. Overall, this study demonstrates the spontaneous Ca2+ release events and ion channel responses throughout various stimuli.

Understanding Calmodulin Variants Affecting Calcium-Dependent Inactivation of L-Type Calcium Channels through Whole-Cell Simulation of the Cardiac Ventricular Myocyte

Mutations in the calcium-sensing protein calmodulin (CaM) have been linked to two cardiac arrhythmia diseases, Long QT Syndrome 14 (LQT14) and Catecholaminergic Polymorphic Ventricular Tachycardia Type 4 (CPVT4), with varying degrees of severity. Functional characterization of the CaM mutants most strongly associated with LQT14 show a clear disruption of the calcium-dependent inactivation (CDI) of the L-Type calcium channel (LCC). CPVT4 mutants on the other hand are associated with changes in their affinity to the ryanodine receptor. In clinical studies, some variants have been associated with both CPVT4 and LQT15. This study uses simulations in a model for excitation-contraction coupling in the rat ventricular myocytes to understand how LQT14 variant might give the functional phenotype similar to CPVT4. Changing the CaM-dependent transition rate by a factor of 0.75 corresponding to the D96V variant and by a factor of 0.90 corresponding to the F142L or N98S variants, in a physiologically based stochastic model of the LCC prolonger, the action potential duration changed by a small amount in a cardiac myocyte but did not disrupt CICR at 1, 2, and 4 Hz. Under beta-adrenergic simulation abnormal excitation-contraction coupling was observed above 2 Hz pacing for the mutant CaM. The same conditions applied under beta-adrenergic stimulation led to the rapid onset of arrhythmia in the mutant CaM simulations. Simulations with the LQT14 mutations under the conditions of rapid pacing with beta-adrenergic stimulation drives the cardiac myocyte toward an arrhythmic state known as Ca2+ overload. These simulations provide a mechanistic link to a disease state for LQT14-associated mutations in CaM to yield a CPVT4 phenotype. The results show that small changes to the CaM-regulated inactivation of LCC promote arrhythmia and underscore the significance of CDI in proper heart function.

Grapefruit Flavonoid Naringenin Sex-Dependently Modulates Action Potential in an In Silico Human Ventricular Cardiomyocyte Model

Recent in vitro studies showed that grapefruit (Citrus × paradisi) flavonoid naringenin alters the function of cardiac ion channels. Here, we explored the effect of naringenin on cardiomyocyte action potentials (APs) using a detailed in silico model of ventricular electrophysiology. Concentration-dependent effects of naringenin on seven major cardiac ion channels were incorporated into the Tomek–Rodriguez modification of O’Hara–Rudy (ToR-ORd) human ventricular endocardium model. To investigate the sex-dependent effect of naringenin, previously reported sex-specific ionic modifications were implemented into the model. Next, populations of 1000 models accommodating intercellular variability were generated. The results show, naringenin at various concentrations prolonged AP duration (APD) in male and female cardiomyocytes. Pacing cells at higher frequencies abbreviated APD differently in males versus females; for example, at 3 Hz, 50 μM naringenin induced AP and calcium alternans only in the female cardiomyocyte. Finally, a population modeling approach corroborated that naringenin significantly prolonged APD in a concentration-dependent manner, with a larger effect in females than in males. In conclusion, our study demonstrates that the APD-prolonging effect of naringenin was larger in females, and that pacing at faster rates induces AP alternation earlier in females, suggesting a potentially higher proarrhythmic risk of naringenin in females than in males.

SGLT2 Inhibitor Empagliflozin Modulates Ion Channels in Adult Zebrafish Heart

Empagliflozin, an inhibitor of sodium-glucose co-transporter 2 (iSGLT2), improves cardiovascular outcomes in patients with and without diabetes and possesses an antiarrhythmic activity. However, the mechanisms of these protective effects have not been fully elucidated. This study aimed to explore the impact of empagliflozin on ion channel activity and electrophysiological characteristics in the ventricular myocardium. The main cardiac ionic currents (I_Na, I_CaL, I_CaT, I_Kr, I_Ks) and action potentials (APs) were studied in zebrafish. Whole-cell currents were measured using the patch clamp method in the isolated ventricular cardiomyocytes. The conventional sharp glass microelectrode technique was applied for the recording of APs from the ventricular myocardium of the excised heart. Empagliflozin pretreatment compared to the control group enhanced potassium I_Kr step current density in the range of testing potentials from 0 to +30 mV, I_Kr tail current density in the range of testing potentials from +10 to +70 mV, and I_Ks current density in the range of testing potentials from −10 to +20 mV. Moreover, in the ventricular myocardium, empagliflozin pretreatment shortened AP duration APD as shown by reduced APD50 and APD90. Empagliflozin had no influence on sodium (I_Na) and L- and T-type calcium currents (I_CaL and I_CaT) in zebrafish ventricular cardiomyocytes. Thus, we conclude that empagliflozin increases the rapid and slow components of delayed rectifier K⁺ current (I_Kr and I_Ks). This mechanism could be favorable for cardiac protection.

Critical Requirements for the Initiation of a Cardiac Arrhythmia in Rat Ventricle: How Many Myocytes?

Cardiovascular disease is the leading cause of death worldwide due in a large part to arrhythmia. In order to understand how calcium dynamics play a role in arrhythmogenesis, normal and dysfunctional Ca²⁺ signaling in a subcellular, cellular, and tissued level is examined using cardiac ventricular myocytes at a high temporal and spatial resolution using multiscale computational modeling. Ca²⁺ sparks underlie normal excitation-contraction coupling. However, under pathological conditions, Ca²⁺ sparks can combine to form Ca²⁺ waves. These propagating elevations of (Ca²⁺)_i can activate an inward Na⁺-Ca²⁺ exchanger current (I_NCX) that contributes to early after-depolarization (EADs) and delayed after-depolarizations (DADs). However, how cellular currents lead to full depolarization of the myocardium and how they initiate extra systoles is still not fully understood. This study explores how many myocytes must be entrained to initiate arrhythmogenic depolarizations in biophysically detailed computational models. The model presented here suggests that only a small number of myocytes must activate in order to trigger an arrhythmogenic propagating action potential. These conditions were examined in 1-D, 2-D, and 3-D considering heart geometry. The depolarization of only a few hundred ventricular myocytes is required to trigger an ectopic depolarization. The number decreases under disease conditions such as heart failure. Furthermore, in geometrically restricted parts of the heart such as the thin muscle strands found in the trabeculae and papillary muscle, the number of cells needed to trigger a propagating depolarization falls even further to less than ten myocytes.