A Mathematical Model of the Mouse Atrial Myocyte With Inter-Atrial Electrophysiological Heterogeneity

Pro-Arrhythmic Effects of Discontinuous Conduction at the Purkinje Fiber-Ventricle Junction Arising From Heart Failure-Induced Ionic Remodeling - Insights From Computational Modelling

Heart failure is associated with electrical remodeling of the electrical properties and kinetics of the ion channels and transporters that are responsible for cardiac action potentials. However, it is still unclear whether heart failure-induced ionic remodeling can affect the conduction of excitation waves at the Purkinje fiber-ventricle junction contributing to pro-arrhythmic effects of heart failure, as the complexity of the heart impedes a detailed experimental analysis. The aim of this study was to employ computational models to investigate the pro-arrhythmic effects of heart failure-induced ionic remodeling on the cardiac action potentials and excitation wave conduction at the Purkinje fiber-ventricle junction. Single cell models of canine Purkinje fiber and ventricular myocytes were developed for control and heart failure. These single cell models were then incorporated into one-dimensional strand and three-dimensional wedge models to investigate the effects of heart failure-induced remodeling on propagation of action potentials in Purkinje fiber and ventricular tissue and at the Purkinje fiber-ventricle junction. This revealed that heart failure-induced ionic remodeling of Purkinje fiber and ventricular tissue reduced conduction safety and increased tissue vulnerability to the genesis of the unidirectional conduction block. This was marked at the Purkinje fiber-ventricle junction, forming a potential substrate for the genesis of conduction failure that led to re-entry. This study provides new insights into proarrhythmic consequences of heart failure-induced ionic remodeling.

Immediate and Delayed Response of Simulated Human Atrial Myocytes to Clinically-Relevant Hypokalemia

Although plasma electrolyte levels are quickly and precisely regulated in the mammalian cardiovascular system, even small transient changes in K⁺, Na⁺, Ca²⁺, and/or Mg²⁺ can significantly alter physiological responses in the heart, blood vessels, and intrinsic (intracardiac) autonomic nervous system. We have used mathematical models of the human atrial action potential (AP) to explore the electrophysiological mechanisms that underlie changes in resting potential (V_r) and the AP following decreases in plasma K⁺, [K⁺]_o, that were selected to mimic clinical hypokalemia. Such changes may be associated with arrhythmias and are commonly encountered in patients (i) in therapy for hypertension and heart failure; (ii) undergoing renal dialysis; (iii) with any disease with acid-base imbalance; or (iv) post-operatively. Our study emphasizes clinically-relevant hypokalemic conditions, corresponding to [K⁺]_o reductions of approximately 1.5 mM from the normal value of 4 to 4.5 mM. We show how the resulting electrophysiological responses in human atrial myocytes progress within two distinct time frames: (i) Immediately after [K⁺]_o is reduced, the K⁺-sensing mechanism of the background inward rectifier current (I_K1) responds. Specifically, its highly non-linear current-voltage relationship changes significantly as judged by the voltage dependence of its region of outward current. This rapidly alters, and sometimes even depolarizes, V_r and can also markedly prolong the final repolarization phase of the AP, thus modulating excitability and refractoriness. (ii) A second much slower electrophysiological response (developing 5-10 minutes after [K⁺]_o is reduced) results from alterations in the intracellular electrolyte balance. A progressive shift in intracellular [Na⁺]_i causes a change in the outward electrogenic current generated by the Na⁺/K⁺ pump, thereby modifying V_r and AP repolarization and changing the human atrial electrophysiological substrate. In this study, these two effects were investigated quantitatively, using seven published models of the human atrial AP. This highlighted the important role of I_K1 rectification when analyzing both the mechanisms by which [K⁺]_o regulates V_r and how the AP waveform may contribute to "trigger" mechanisms within the proarrhythmic substrate. Our simulations complement and extend previous studies aimed at understanding key factors by which decreases in [K⁺]_o can produce effects that are known to promote atrial arrhythmias in human hearts.

Evolution of mathematical models of cardiomyocyte electrophysiology

Advanced computational techniques and mathematical modeling have become more and more important to the study of cardiac electrophysiology. In this review, we provide a brief history of the evolution of cardiomyocyte electrophysiology models and highlight some of the most important ones that had a major impact on our understanding of the electrical activity of the myocardium and associated transmembrane ion fluxes in normal and pathological states. We also present the use of these models in the study of various arrhythmogenesis mechanisms, particularly the integration of experimental pharmacology data into advanced humanized models for in silico proarrhythmogenic risk prediction as an essential component of the Comprehensive in vitro Proarrhythmia Assay (CiPA) drug safety paradigm.

Computational modeling approaches to cAMP/PKA signaling in cardiomyocytes

The cAMP/PKA pathway is a fundamental regulator of excitation-contraction coupling in cardiomyocytes. Activation of cAMP has a variety of downstream effects on cardiac function including enhanced contraction, accelerated relaxation, adaptive stress response, mitochondrial regulation, and gene transcription. Experimental advances have shed light on the compartmentation of cAMP and PKA, which allow for control over the varied targets of these second messengers and is disrupted in heart failure conditions. Computational modeling is an important tool for understanding the spatial and temporal complexities of this system. In this review article, we outline the advances in computational modeling that have allowed for deeper understanding of cAMP/PKA dynamics in the cardiomyocyte in health and disease, and explore new modeling frameworks that may bring us closer to a more complete understanding of this system. We outline various compartmental and spatial signaling models that have been used to understand how β-adrenergic signaling pathways function in a variety of simulation conditions. We also discuss newer subcellular models of cardiovascular function that may be used as templates for the next phase of computational study of cAMP and PKA in the heart, and outline open challenges which are important to consider in future models.

The Role of CaMKII Overexpression and Oxidation in Atrial Fibrillation-A Simulation Study

This simulation study aims to investigate how the Calcium/calmodulin-dependent protein kinase II (CaMKII) overexpression and oxidation would influence the cardiac electrophysiological behavior and its arrhythmogenic mechanism in atria. A new-built CaMKII oxidation module and a refitted CaMKII overexpression module were integrated into a mouse atrial cell model for analyzing cardiac electrophysiological variations in action potential (AP) characteristics and intracellular Ca²⁺ cycling under different conditions. Simulation results showed that CaMKII overexpression significantly increased the phosphorylation level of its downstream target proteins, resulting in prolonged AP and smaller calcium transient amplitude, and impaired the Ca²⁺ cycling stability. These effects were exacerbated by extra reactive oxygen species, which oxidized CaMKII and led to continuous high CaMKII activation in both systolic and diastolic phases. Intracellular Ca²⁺ depletion and sustained delayed afterdepolarizations (DADs) were observed under co-existing CaMKII overexpression and oxidation, which could be effectively reversed by clamping the phosphorylation level of ryanodine receptor (RyR). We also found that the stability of RyR release highly depended on a delicate balance between the level of RyR phosphorylation and sarcoplasmic reticulum Ca²⁺ concentration, which was closely related to the genesis of DADs. We concluded that the CaMKII overexpression and oxidation have a synergistic role in increasing the activity of CaMKII, and the unstable RyR may be the key downstream target in the CaMKII arrhythmogenic mechanism. Our simulation provides detailed mechanistic insights into the arrhythmogenic effect of CaMKII overexpression and oxidation, which suggests CaMKII as a promising target in the therapy of atrial fibrillation.