Ca²⁺-induced delayed afterdepolarizations are triggered by dyadic subspace Ca2²⁺ affirming that increasing SERCA reduces aftercontractions

Stress-Induced Proteasome Sub-Cellular Translocation in Cardiomyocytes Causes Altered Intracellular Calcium Handling and Arrhythmias

The ubiquitin-proteasome system (UPS) is an essential mechanism responsible for the selective degradation of substrate proteins via their conjugation with ubiquitin. Since cardiomyocytes have very limited self-renewal capacity, as they are prone to protein damage due to constant mechanical and metabolic stress, the UPS has a key role in cardiac physiology and pathophysiology. While altered proteasomal activity contributes to a variety of cardiac pathologies, such as heart failure and ischemia/reperfusion injury (IRI), the environmental cues affecting its activity are still unknown, and they are the focus of this work. Following a recent study by Ciechanover's group showing that amino acid (AA) starvation in cultured cancer cell lines modulates proteasome intracellular localization and activity, we tested two hypotheses in human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs, CMs): (i) AA starvation causes proteasome translocation in CMs, similarly to the observation in cultured cancer cell lines; (ii) manipulation of subcellular proteasomal compartmentalization is associated with electrophysiological abnormalities in the form of arrhythmias, mediated via altered intracellular Ca2+ handling. The major findings are: (i) starving CMs to AAs results in proteasome translocation from the nucleus to the cytoplasm, while supplementation with the aromatic amino acids tyrosine (Y), tryptophan (W) and phenylalanine (F) (YWF) inhibits the proteasome recruitment; (ii) AA-deficient treatments cause arrhythmias; (iii) the arrhythmias observed upon nuclear proteasome sequestration(-AA+YWF) are blocked by KB-R7943, an inhibitor of the reverse mode of the sodium-calcium exchanger NCX; (iv) the retrograde perfusion of isolated rat hearts with AA starvation media is associated with arrhythmias. Collectively, our novel findings describe a newly identified mechanism linking the UPS to arrhythmia generation in CMs and whole hearts.

Mechanisms of pulmonary arterial hypertension-induced atrial fibrillation: insights from multi-scale models of the human atria

This study aimed to use multi-scale atrial models to investigate pulmonary arterial hypertension (PAH)-induced atrial fibrillation mechanisms. The results of our computer simulations revealed that, at the single-cell level, PAH-induced remodelling led to a prolonged action potential (AP) (ΔAPD: 49.6 ms in the right atria (RA) versus 41.6 ms in the left atria (LA)) and an increased calcium transient (CaT) (ΔCaT: 7.5 × 10-2 µM in the RA versus 0.9 × 10-3 µM in the LA). Moreover, heterogeneous remodelling increased susceptibility to afterdepolarizations, particularly in the RA. At the tissue level, we observed a significant reduction in conduction velocity (CV) (ΔCV: -0.5 m s-1 in the RA versus -0.05 m s-1 in the LA), leading to a shortened wavelength in the RA, but not in the LA. Additionally, afterdepolarizations in the RA contributed to enhanced repolarization dispersion and facilitated unidirectional conduction block. Furthermore, the increased fibrosis in the RA amplified the likelihood of excitation wave breakdown and the occurrence of sustained re-entries. Our results indicated that the RA is characterized by increased susceptibility to afterdepolarizations, slow conduction, reduced wavelength and upregulated fibrosis. These findings shed light on the underlying factors that may promote atrial fibrillation in patients with PAH.

Estimating the probability of early afterdepolarizations and predicting arrhythmic risk associated with long QT syndrome type 1 mutations

Early afterdepolarizations (EADs) are action potential (AP) repolarization abnormalities that can trigger lethal arrhythmias. Simulations using biophysically detailed cardiac myocyte models can reveal how model parameters influence the probability of these cellular arrhythmias; however, such analyses can pose a huge computational burden. We have previously developed a highly simplified approach in which logistic regression models (LRMs) map parameters of complex cell models to the probability of ectopic beats. Here, we extend this approach to predict the probability of EADs (P(EAD)) as a mechanistic metric of arrhythmic risk. We use the LRM to investigate how changes in parameters of the slow-activating delayed rectifier current (IKs) affect P(EAD) for 17 different long QT syndrome type 1 (LQTS1) mutations. In this LQTS1 clinical arrhythmic risk prediction task, we compared P(EAD) for these 17 mutations with two other recently published model-based arrhythmia risk metrics (AP morphology metric across populations of myocyte models and transmural repolarization prolongation based on a one-dimensional [1D] tissue-level model). These model-based risk metrics yield similar prediction performance; however, each fails to stratify clinical risk for a significant number of the 17 studied LQTS1 mutations. Nevertheless, an interpretable ensemble model using multivariate linear regression built by combining all of these model-based risk metrics successfully predicts the clinical risk of 17 mutations. These results illustrate the potential of computational approaches in arrhythmia risk prediction.

Bacterial toxins and heart function: heat-labile Escherichia coli enterotoxin B promotes changes in cardiac function with possible relevance for sudden cardiac death

Bacterial toxins can cause cardiomyopathy, though it is not its most common cause. Some bacterial toxins can form pores in the membrane of cardiomyocytes, while others can bind to membrane receptors. Enterotoxigenic E. coli can secrete enterotoxins, including heat-resistant (ST) or labile (LT) enterotoxins. LT is an AB5-type toxin that can bind to specific cell receptors and disrupt essential host functions, causing several common conditions, such as certain diarrhea. The pentameric B subunit of LT, without A subunit (LTB), binds specifically to certain plasma membrane ganglioside receptors, found in lipid rafts of cardiomyocytes. Isolated guinea pig hearts and cardiomyocytes were exposed to different concentrations of purified LTB. In isolated hearts, mechanical and electrical alternans and an increment of heart rate variability, with an IC50 of ~0.2 μg/ml LTB, were observed. In isolated cardiomyocytes, LTB promoted significant decreases in the amplitude and the duration of action potentials. Na+ currents were inhibited whereas L-type Ca2+ currents were augmented at their peak and their fast inactivation was promoted. Delayed rectifier K+ currents decreased. Measurements of basal Ca2+ or Ca2+ release events in cells exposed to LTB suggest that LTB impairs Ca2+ homeostasis. Impaired calcium homeostasis is linked to sudden cardiac death. The results are consistent with the recent view that the B subunit is not merely a carrier of the A subunit, having a role explaining sudden cardiac death in children (SIDS) infected with enterotoxigenic E. coli, explaining several epidemiological findings that establish a strong relationship between SIDS and ETEC E. coli.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-023-01100-6.

Using hiPSC-CMs to Examine Mechanisms of Catecholaminergic Polymorphic Ventricular Tachycardia

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is a potentially lethal inherited cardiac arrhythmia condition, triggered by physical or acute emotional stress, that predominantly expresses early in life. Gain-of-function mutations in the cardiac ryanodine receptor gene (RYR2) account for the majority of CPVT cases, causing substantial disruption of intracellular calcium (Ca²⁺ ) homeostasis particularly during the periods of β-adrenergic receptor stimulation. However, the highly variable penetrance, patient outcomes, and drug responses observed in clinical practice remain unexplained, even for patients with well-established founder RyR2 mutations. Therefore, investigation of the electrophysiological consequences of CPVT-causing RyR2 mutations is crucial to better understand the pathophysiology of the disease. The development of strategies for reprogramming human somatic cells to human induced pluripotent stem cells (hiPSCs) has provided a unique opportunity to study inherited arrhythmias, due to the ability of hiPSCs to differentiate down a cardiac lineage. Employment of genome editing enables generation of disease-specific cell lines from healthy and diseased patient-derived hiPSCs, which subsequently can be differentiated into cardiomyocytes. This paper describes the means for establishing an hiPSC-based model of CPVT in order to recapitulate the disease phenotype in vitro and investigate underlying pathophysiological mechanisms. The framework of this approach has the potential to contribute to disease modeling and personalized medicine using hiPSC-derived cardiomyocytes. © 2021 Wiley Periodicals LLC.

Estimating ectopic beat probability with simplified statistical models that account for experimental uncertainty

Ectopic beats (EBs) are cellular arrhythmias that can trigger lethal arrhythmias. Simulations using biophysically-detailed cardiac myocyte models can reveal how model parameters influence the probability of these cellular arrhythmias, however such analyses can pose a huge computational burden. Here, we develop a simplified approach in which logistic regression models (LRMs) are used to define a mapping between the parameters of complex cell models and the probability of EBs (P(EB)). As an example, in this study, we build an LRM for P(EB) as a function of the initial value of diastolic cytosolic Ca²⁺ concentration ([Ca²⁺]_iⁱⁿⁱ), the initial state of sarcoplasmic reticulum (SR) Ca²⁺ load ([Ca²⁺]_SRⁱⁿⁱ), and kinetic parameters of the inward rectifier K⁺ current (I_K1) and ryanodine receptor (RyR). This approach, which we refer to as arrhythmia sensitivity analysis, allows for evaluation of the relationship between these arrhythmic event probabilities and their associated parameters. This LRM is also used to demonstrate how uncertainties in experimentally measured values determine the uncertainty in P(EB). In a study of the role of [Ca²⁺]_SRⁱⁿⁱ uncertainty, we show a special property of the uncertainty in P(EB), where with increasing [Ca²⁺]_SRⁱⁿⁱ uncertainty, P(EB) uncertainty first increases and then decreases. Lastly, we demonstrate that I_K1 suppression, at the level that occurs in heart failure myocytes, increases P(EB).

An ectopic beat is an abnormal cellular electrical event which can trigger dangerous arrhythmias in the heart. Complex biophysical models of the cardiac myocyte can be used to reveal how cell properties affect the probability of ectopic beats. However, such analyses can pose a huge computational burden. We develop a simplified approach that enables a highly complex biophysical model to be reduced to a rather simple statistical model from which the functional relationship between myocyte model parameters and the probability of an ectopic beat is determined. We refer to this approach as arrhythmia sensitivity analysis. Given the efficiency of our approach, we also use it to demonstrate how uncertainties in experimentally measured myocyte model parameters determine the uncertainty in ectopic beat probability. We find that, with increasing model parameter uncertainty, the uncertainty in probability of ectopic beat first increases and then decreases. In general, our approach can efficiently analyze the relationship between cardiac myocyte parameters and the probability of ectopic beats and can be used to study how uncertainty of these cardiac myocyte parameters influences the ectopic beat probability.

Constitutively Activating GNAS Somatic Mutation in Right Ventricular Outflow Tract Tachycardia

[Figure: see text].

Electrophysiological Remodeling: Cardiac T-Tubules and ß-Adrenoceptors

Beta-adrenoceptors (βAR) are often viewed as archetypal G-protein coupled receptors. Over the past fifteen years, investigations in cardiovascular biology have provided remarkable insights into this receptor family. These studies have shifted pharmacological dogma, from one which centralized the receptor to a new focus on structural micro-domains such as caveolae and t-tubules. Important studies have examined, separately, the structural compartmentation of ion channels and βAR. Despite links being assumed, relatively few studies have specifically examined the direct link between structural remodeling and electrical remodeling with a focus on βAR. In this review, we will examine the nature of receptor and ion channel dysfunction on a substrate of cardiomyocyte microdomain remodeling, as well as the likely ramifications for cardiac electrophysiology. We will then discuss the advances in methodologies in this area with a specific focus on super-resolution microscopy, fluorescent imaging, and new approaches involving microdomain specific, polymer-based agonists. The advent of powerful computational modelling approaches has allowed the science to shift from purely empirical work, and may allow future investigations based on prediction. Issues such as the cross-reactivity of receptors and cellular heterogeneity will also be discussed. Finally, we will speculate as to the potential developments within this field over the next ten years.

Implications of SGLT Inhibition on Redox Signalling in Atrial Fibrillation

Atrial fibrillation (AF) is the most common sustained (atrial) arrhythmia, a considerable global health burden and often associated with heart failure. Perturbations of redox signalling in cardiomyocytes provide a cellular substrate for the manifestation and maintenance of atrial arrhythmias. Several clinical trials have shown that treatment with sodium-glucose linked transporter inhibitors (SGLTi) improves mortality and hospitalisation in heart failure patients independent of the presence of diabetes. Post hoc analysis of the DECLARE-TIMI 58 trial showed a 19% reduction in AF in patients with diabetes mellitus (hazard ratio, 0.81 (95% confidence interval: 0.68-0.95), n = 17.160) upon treatment with SGLTi, regardless of pre-existing AF or heart failure and independent from blood pressure or renal function. Accordingly, ongoing experimental work suggests that SGLTi not only positively impact heart failure but also counteract cellular ROS production in cardiomyocytes, thereby potentially altering atrial remodelling and reducing AF burden. In this article, we review recent studies investigating the effect of SGLTi on cellular processes closely interlinked with redox balance and their potential effects on the onset and progression of AF. Despite promising insight into SGLTi effect on Ca²⁺ cycling, Na⁺ balance, inflammatory and fibrotic signalling, mitochondrial function and energy balance and their potential effect on AF, the data are not yet conclusive and the importance of individual pathways for human AF remains to be established. Lastly, an overview of clinical studies investigating SGLTi in the context of AF is provided.

MD1 Depletion Predisposes to Ventricular Arrhythmias in the Setting of Myocardial Infarction

BACKGROUND

Myeloid differentiation protein 1 (MD1) is expressed in the human heart and is a negative regulator of Toll-like receptor 4 (TLR4) signalling. MD1 exerts anti-arrhythmic effects.

AIM

The aim of this study was to determine the role of MD1 in myocardial infarction (MI)-related ventricular arrhythmias (VAs).

METHOD

Myocardial infarction was induced by surgical ligation of the left anterior coronary artery in MD1 knockout (KO) mice and their wild-type littermates. Myocardial infarction-induced vulnerability to VAs and its underlying mechanisms were evaluated.

RESULTS

Myeloid differentiation protein 1 was downregulated in the MI mice. Myeloid differentiation protein 1 deficiency decreased post-MI left ventricular (LV) function and increased the infarct size. The MI mice exhibited prolonged action potential duration (APD), enhanced APD alternans thresholds, and a higher incidence of VAs. Myocardial infarction-induced LV fibrosis and inflammation decreased the expression levels of Kv4.2, Kv4.3, Kv1.5, and Kv2.1, increased Cav1.2 expression, and disturbed Ca²⁺ handling protein expression. These MI-induced adverse effects were further exacerbated in KO mice. Mechanistically, MD1 deletion markedly enhanced the activation of the TLR4/calmodulin-dependent protein kinase II (CaMKII) signalling pathway in post-MI mice.

CONCLUSIONS

Myeloid differentiation protein 1 deletion increases the vulnerability to VAs in post-MI mice. This is mainly caused by the aggravated maladaptive LV fibrosis and inflammation and interference with the expressions of ion channels and Ca²⁺ handling proteins, which is related to enhanced activation of the TLR4/CaMKII signalling pathway.

Targeting Ca2 + Handling Proteins for the Treatment of Heart Failure and Arrhythmias

Diseases of the heart, such as heart failure and cardiac arrhythmias, are a growing socio-economic burden. Calcium (Ca²⁺) dysregulation is key hallmark of the failing myocardium and has long been touted as a potential therapeutic target in the treatment of a variety of cardiovascular diseases (CVD). In the heart, Ca²⁺ is essential for maintaining normal cardiac function through the generation of the cardiac action potential and its involvement in excitation contraction coupling. As such, the proteins which regulate Ca²⁺ cycling and signaling play a vital role in maintaining Ca²⁺ homeostasis. Changes to the expression levels and function of Ca²⁺-channels, pumps and associated intracellular handling proteins contribute to altered Ca²⁺ homeostasis in CVD. The remodeling of Ca²⁺-handling proteins therefore results in impaired Ca²⁺ cycling, Ca²⁺ leak from the sarcoplasmic reticulum and reduced Ca²⁺ clearance, all of which contributes to increased intracellular Ca²⁺. Currently, approved treatments for targeting Ca²⁺ handling dysfunction in CVD are focused on Ca²⁺ channel blockers. However, whilst Ca²⁺ channel blockers have been successful in the treatment of some arrhythmic disorders, they are not universally prescribed to heart failure patients owing to their ability to depress cardiac function. Despite the progress in CVD treatments, there remains a clear need for novel therapeutic approaches which are able to reverse pathophysiology associated with heart failure and arrhythmias. Given that heart failure and cardiac arrhythmias are closely associated with altered Ca²⁺ homeostasis, this review will address the molecular changes to proteins associated with both Ca²⁺-handling and -signaling; their potential as novel therapeutic targets will be discussed in the context of pre-clinical and, where available, clinical data.

Calibration of ionic and cellular cardiac electrophysiology models

Cardiac electrophysiology models are among the most mature and well-studied mathematical models of biological systems. This maturity is bringing new challenges as models are being used increasingly to make quantitative rather than qualitative predictions. As such, calibrating the parameters within ion current and action potential (AP) models to experimental data sets is a crucial step in constructing a predictive model. This review highlights some of the fundamental concepts in cardiac model calibration and is intended to be readily understood by computational and mathematical modelers working in other fields of biology. We discuss the classic and latest approaches to calibration in the electrophysiology field, at both the ion channel and cellular AP scales. We end with a discussion of the many challenges that work to date has raised and the need for reproducible descriptions of the calibration process to enable models to be recalibrated to new data sets and built upon for new studies. This article is categorized under: Analytical and Computational Methods > Computational Methods Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.

Afterdepolarizations and abnormal calcium handling in atrial myocytes with modulated SERCA uptake: a sensitivity analysis of calcium handling channels

Delayed afterdepolarizations (DADs) and spontaneous depolarizations (SDs) are typically triggered by spontaneous diastolic Ca²⁺ release from the sarcoplasmic reticulum (SR) which is caused by an elevated SR Ca²⁺-ATPase (SERCA) uptake and dysfunctional ryanodine receptors. However, recent studies on the T-box transcription factor gene (TBX5) demonstrated that abnormal depolarizations could occur despite a reduced SERCA uptake. Similar findings have also been reported in experimental or clinical studies of diabetes and heart failure. To investigate the sensitivity of SERCA in the genesis of DADs/SDs as well as its dependence on other Ca²⁺ handling channels, we performed systematic analyses using the Maleckar et al. model. Results showed that the modulation of SERCA alone cannot trigger abnormal depolarizations, but can instead affect the interdependency of other Ca²⁺ handling channels in triggering DADs/SDs. Furthermore, we discovered the existence of a threshold value for the intracellular concentration of Ca²⁺ ([Ca²⁺]_i) for abnormal depolarizations, which is modulated by the maximum SERCA uptake and the concentration of Ca²⁺ in the uptake and release compartments in the SR ([Ca²⁺]_up and [Ca²⁺]_rel). For the first time, our modelling study reconciles different mechanisms of abnormal depolarizations in the setting of 'lone' AF, reduced TBX5, diabetes and heart failure, and may lead to more targeted treatment for these patients. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Human Purkinje in silico model enables mechanistic investigations into automaticity and pro-arrhythmic abnormalities

Cardiac Purkinje cells (PCs) are implicated in lethal arrhythmias caused by cardiac diseases, mutations, and drug action. However, the pro-arrhythmic mechanisms in PCs are not entirely understood, particularly in humans, as most investigations are conducted in animals. The aims of this study are to present a novel human PCs electrophysiology biophysically-detailed computational model, and to disentangle ionic mechanisms of human Purkinje-related electrophysiology, pacemaker activity and arrhythmogenicity. The new Trovato2020 model incorporates detailed Purkinje-specific ionic currents and Ca2+ handling, and was developed, calibrated and validated using human experimental data acquired at multiple frequencies, both in control conditions and following drug application. Multiscale investigations were performed in a Purkinje cell, in fibre and using an experimentally-calibrated population of PCs to evaluate biological variability. Simulations demonstrate the human Purkinje Trovato2020 model is the first one to yield: (i) all key AP features consistent with human Purkinje recordings; (ii) Automaticity with funny current up-regulation (iii) EADs at slow pacing and with 85% hERG block; (iv) DADs following fast pacing; (v) conduction velocity of 160 cm/s in a Purkinje fibre, as reported in human. The human in silico PCs population highlights that: (1) EADs are caused by ICaL reactivation in PCs with large inward currents; (2) DADs and triggered APs occur in PCs experiencing Ca2+ accumulation, at fast pacing, caused by large L-type calcium current and small Na+/Ca2+ exchanger. The novel human Purkinje model unlocks further investigations into the role of cardiac Purkinje in ventricular arrhythmias through computer modeling and multiscale simulations.

Simulation of the Effects of Extracellular Calcium Changes Leads to a Novel Computational Model of Human Ventricular Action Potential With a Revised Calcium Handling

The importance of electrolyte concentrations for cardiac function is well established. Electrolyte variations can lead to arrhythmias onset, due to their important role in the action potential (AP) genesis and in maintaining cell homeostasis. However, most of the human AP computer models available in literature were developed with constant electrolyte concentrations, and fail to simulate physiological changes induced by electrolyte variations. This is especially true for Ca2+, even in the O'Hara-Rudy model (ORd), one of the most widely used models in cardiac electrophysiology. Therefore, the present work develops a new human ventricular model (BPS2020), based on ORd, able to simulate the inverse dependence of AP duration (APD) on extracellular Ca2+ concentration ([Ca2+]o), and APD rate dependence at 4 mM extracellular K+. The main changes needed with respect to ORd are: (i) an increased sensitivity of L-type Ca2+ current inactivation to [Ca2+]o; (ii) a single compartment description of the sarcoplasmic reticulum; iii) the replacement of Ca2+ release. BPS2020 is able to simulate the physiological APD-[Ca2+]o relationship, while also retaining the well-reproduced properties of ORd (APD rate dependence, restitution, accommodation and current block effects). We also used BPS2020 to generate an experimentally-calibrated population of models to investigate: (i) the occurrence of repolarization abnormalities in response to hERG current block; (ii) the rate adaptation variability; (iii) the occurrence of alternans and delayed after-depolarizations at fast pacing. Our results indicate that we successfully developed an improved version of ORd, which can be used to investigate electrophysiological changes and pro-arrhythmic abnormalities induced by electrolyte variations and current block at multiple rates and at the population level.

In silico investigation of the mechanisms underlying atrial fibrillation due to impaired Pitx2

Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia and is a major cause of stroke and morbidity. Recent genome-wide association studies have shown that paired-like homeodomain transcription factor 2 (Pitx2) to be strongly associated with AF. However, the mechanisms underlying Pitx2 modulated arrhythmogenesis and variable effectiveness of antiarrhythmic drugs (AADs) in patients in the presence or absence of impaired Pitx2 expression remain unclear. We have developed multi-scale computer models, ranging from a single cell to tissue level, to mimic control and Pitx2-knockout atria by incorporating recent experimental data on Pitx2-induced electrical and structural remodeling in humans, as well as the effects of AADs. The key findings of this study are twofold. We have demonstrated that shortened action potential duration, slow conduction and triggered activity occur due to electrical and structural remodelling under Pitx2 deficiency conditions. Notably, the elevated function of calcium transport ATPase increases sarcoplasmic reticulum Ca2+ concentration, thereby enhancing susceptibility to triggered activity. Furthermore, heterogeneity is further elevated due to Pitx2 deficiency: 1) Electrical heterogeneity between left and right atria increases; and 2) Increased fibrosis and decreased cell-cell coupling due to structural remodelling slow electrical propagation and provide obstacles to attract re-entry, facilitating the initiation of re-entrant circuits. Secondly, our study suggests that flecainide has antiarrhythmic effects on AF due to impaired Pitx2 by preventing spontaneous calcium release and increasing wavelength. Furthermore, our study suggests that Na+ channel effects alone are insufficient to explain the efficacy of flecainide. Our study may provide the mechanisms underlying Pitx2-induced AF and possible explanation behind the AAD effects of flecainide in patients with Pitx2 deficiency.

Inhibition of PI3Kinase-α is pro-arrhythmic and associated with enhanced late Na+ current, contractility, and Ca2+ release in murine hearts

BACKGROUND

Phosphoinositide 3-kinase α (PI3Kα) is a proto-oncogene with high activity in the heart. BYL719 (BYL) is a PI3Kα-selective small molecule inhibitor and a prospective drug for advanced solid tumors. We investigated whether acute pharmacological inhibition of PI3Kα has pro-arrhythmic effects.

METHODS & RESULTS

In isolated wild-type (WT) cardiomyocytes, pharmacological inhibition of PI3Kα (BYL719) increased contractility by 28%, Ca²⁺ release by 20%, and prolonged action potential (AP) repolarization by 10-15%. These effects of BYL719 were abolished by inhibition of reverse-mode Na⁺/Ca²⁺ exchanger (NCX) (KB-R7943) or by inhibition of late Na⁺ current (I_Na-L) (ranolazine). BYL719 had no effect on PI3Kα-deficient cardiomyocytes, suggesting BYL719 effects were PI3Kα-dependent and mediated via NCX and I_Na-L. I_Na-L was suppressed by activation of PI3Kα, application of exogenous intracellular PIP3, or ranolazine. Investigation of AP and Ca²⁺ release in whole heart preparations using epicardial optical mapping showed that inhibition of PI3Kα similarly led to prolongation of AP and enhancement of Ca²⁺ release. In hearts of PI3Kα-deficient mice, β-adrenergic stimulation in the presence of high Ca²⁺ concentrations and 12-Hz burst pacing led to delayed afterdepolarizations and ventricular fibrillation. In vivo, administration of BYL719 prolonged QT interval [QT_cF (Fridericia) increased by 15%] in WT, but not in PI3Kα-deficient mice.

CONCLUSIONS

Pharmacological inhibition of PI3Kα is arrhythmogenic due to activation of I_Na-L leading to increased sarcoplasmic reticulum Ca²⁺ load and prolonged QT interval. Therefore, monitoring of cardiac electrical activity in patients receiving PI3K inhibitors may provide further insights into the arrhythmogenic potential of PI3Ka inhibition.

Caveolin-3 Microdomain: Arrhythmia Implications for Potassium Inward Rectifier and Cardiac Sodium Channel

In human cardiac ventricular myocytes, caveolin-3 functions as a scaffolding and regulatory protein for signaling molecules and compartmentalizes ion channels. Our lab has recently explored this sub-cellular microdomain and found that potassium inward rectifier Kir2.x is found in association with caveolin-3. The three cardiac Kir2.x isoforms (Kir2.1, Kir2.2, and Kir2.3) are the molecular correlates of I_K1 in the heart, of which Kir2.1 is the dominant isoform in the ventricle. Kir2.1 channels assemble with Kir2.2 and Kir2.3 forming hetero-tetramers that modulate I_K1. I_K1 sets the resting membrane potential and assists with terminal phase 3 ventricular repolarization. In our studies using native human ventricular tissue, Kir2.x co-localizes with caveolin-3 and significance of the association between Kir2.x and caveolin-3 is emphasized in relation to mutations in the gene which encodes caveolin-3, CAV3, associated with Long QT Syndrome 9 (LQT9). LQT9-associated CAV3 mutations cause decreased current density in Kir2.1 and Kir2.2 as homomeric and heteromeric channels, which affects repolarization and membrane potential stability. A portion of Kir2.1 cardiac localization parallels that of the cardiac sodium channel (Nav1.5). This may have implications for Long QT9 in which CAV3 mutations cause an increase in the late current of Nav1.5 (I_Na-L) via nNOS mediated nitrosylation of Nav1.5. In iPS-CMs, expression of LQT9 CAV3 mutations resulted in action potential duration (APD) prolongation and early-after depolarizations (EADs), supporting the arrhythmogenicity of LQT9. To evaluate the combined effect of the CAV3 mutants on I_Na-L and I_K1, we studied both ventricular and Purkinje myocyte mathematical modeling. Interestingly, mathematical ventricular myocytes, similar to iPS-CMs, demonstrated EADs but no sustained arrhythmia. In contrast, Purkinje modeling demonstrated delayed-after depolarizations (DADs) driven mechanism for sustained arrhythmia, dependent on the combined loss of I_K1 and gain of I_Na-L. This finding changes the overall assumed arrhythmia phenotype for LQT9. In future studies, we are exploring caveolar micro-domain disruption in heart failure and how this effects Kir2.x and Nav1.5. Here we review the caveolae cardiac microdomain of Kir2.x and Nav1.5 and explore some of the downstream effects of caveolin-3 and caveolae disruption in specific clinical scenarios.

Arrhythmogenic Mechanisms in Heart Failure: Linking β-Adrenergic Stimulation, Stretch, and Calcium

Heart failure (HF) is associated with elevated sympathetic tone and mechanical load. Both systems activate signaling transduction pathways that increase cardiac output, but eventually become part of the disease process itself leading to further worsening of cardiac function. These alterations can adversely contribute to electrical instability, at least in part due to the modulation of Ca²⁺ handling at the level of the single cardiac myocyte. The major aim of this review is to provide a definitive overview of the links and cross talk between β-adrenergic stimulation, mechanical load, and arrhythmogenesis in the setting of HF. We will initially review the role of Ca²⁺ in the induction of both early and delayed afterdepolarizations, the role that β-adrenergic stimulation plays in the initiation of these and how the propensity for these may be altered in HF. We will then go onto reviewing the current data with regards to the link between mechanical load and afterdepolarizations, the associated mechano-sensitivity of the ryanodine receptor and other stretch activated channels that may be associated with HF-associated arrhythmias. Furthermore, we will discuss how alterations in local Ca²⁺ microdomains during the remodeling process associated the HF may contribute to the increased disposition for β-adrenergic or stretch induced arrhythmogenic triggers. Finally, the potential mechanisms linking β-adrenergic stimulation and mechanical stretch will be clarified, with the aim of finding common modalities of arrhythmogenesis that could be targeted by novel therapeutic agents in the setting of HF.

Automatic Optimization of an in Silico Model of Human iPSC Derived Cardiomyocytes Recapitulating Calcium Handling Abnormalities

The growing importance of human induced pluripotent stem cell-derived cardiomyoyctes (hiPSC-CMs), as patient-specific and disease-specific models for studying cellular cardiac electrophysiology or for preliminary cardiotoxicity tests, generated better understanding of hiPSC-CM biophysical mechanisms and great amount of action potential and calcium transient data. In this paper, we propose a new hiPSC-CM in silico model, with particular attention to Ca²⁺ handling. We used (i) the hiPSC-CM Paci2013 model as starting point, (ii) a new dataset of Ca²⁺ transient measurements to tune the parameters of the inward and outward Ca²⁺ fluxes of sarcoplasmic reticulum, and (iii) an automatic parameter optimization to fit action potentials and Ca²⁺ transients. The Paci2018 model simulates, together with the typical hiPSC-CM spontaneous action potentials, more refined Ca²⁺ transients and delayed afterdepolarizations-like abnormalities, which the old Paci2013 was not able to predict due to its mathematical formulation. The Paci2018 model was validated against (i) the same current blocking experiments used to validate the Paci2013 model, and (ii) recently published data about effects of different extracellular ionic concentrations. In conclusion, we present a new and more versatile in silico model, which will provide a platform for modeling the effects of drugs or mutations that affect Ca²⁺ handling in hiPSC-CMs.

Estimating the probabilities of rare arrhythmic events in multiscale computational models of cardiac cells and tissue

Ectopic heartbeats can trigger reentrant arrhythmias, leading to ventricular fibrillation and sudden cardiac death. Such events have been attributed to perturbed Ca2+ handling in cardiac myocytes leading to spontaneous Ca2+ release and delayed afterdepolarizations (DADs). However, the ways in which perturbation of specific molecular mechanisms alters the probability of ectopic beats is not understood. We present a multiscale model of cardiac tissue incorporating a biophysically detailed three-dimensional model of the ventricular myocyte. This model reproduces realistic Ca2+ waves and DADs driven by stochastic Ca2+ release channel (RyR) gating and is used to study mechanisms of DAD variability. In agreement with previous experimental and modeling studies, key factors influencing the distribution of DAD amplitude and timing include cytosolic and sarcoplasmic reticulum Ca2+ concentrations, inwardly rectifying potassium current (IK1) density, and gap junction conductance. The cardiac tissue model is used to investigate how random RyR gating gives rise to probabilistic triggered activity in a one-dimensional myocyte tissue model. A novel spatial-average filtering method for estimating the probability of extreme (i.e. rare, high-amplitude) stochastic events from a limited set of spontaneous Ca2+ release profiles is presented. These events occur when randomly organized clusters of cells exhibit synchronized, high amplitude Ca2+ release flux. It is shown how reduced IK1 density and gap junction coupling, as observed in heart failure, increase the probability of extreme DADs by multiple orders of magnitude. This method enables prediction of arrhythmia likelihood and its modulation by alterations of other cellular mechanisms.

A Parsimonious Model of the Rabbit Action Potential Elucidates the Minimal Physiological Requirements for Alternans and Spiral Wave Breakup

Elucidating the underlying mechanisms of fatal cardiac arrhythmias requires a tight integration of electrophysiological experiments, models, and theory. Existing models of transmembrane action potential (AP) are complex (resulting in over parameterization) and varied (leading to dissimilar predictions). Thus, simpler models are needed to elucidate the "minimal physiological requirements" to reproduce significant observable phenomena using as few parameters as possible. Moreover, models have been derived from experimental studies from a variety of species under a range of environmental conditions (for example, all existing rabbit AP models incorporate a formulation of the rapid sodium current, INa, based on 30 year old data from chick embryo cell aggregates). Here we develop a simple "parsimonious" rabbit AP model that is mathematically identifiable (i.e., not over parameterized) by combining a novel Hodgkin-Huxley formulation of INa with a phenomenological model of repolarization similar to the voltage dependent, time-independent rectifying outward potassium current (IK). The model was calibrated using the following experimental data sets measured from the same species (rabbit) under physiological conditions: dynamic current-voltage (I-V) relationships during the AP upstroke; rapid recovery of AP excitability during the relative refractory period; and steady-state INa inactivation via voltage clamp. Simulations reproduced several important "emergent" phenomena including cellular alternans at rates > 250 bpm as observed in rabbit myocytes, reentrant spiral waves as observed on the surface of the rabbit heart, and spiral wave breakup. Model variants were studied which elucidated the minimal requirements for alternans and spiral wave break up, namely the kinetics of INa inactivation and the non-linear rectification of IK.The simplicity of the model, and the fact that its parameters have physiological meaning, make it ideal for engendering generalizable mechanistic insight and should provide a solid "building-block" to generate more detailed ionic models to represent complex rabbit electrophysiology.

CaMKII as a target for arrhythmia suppression

Calcium/calmodulin-dependent protein kinase II (CaMKII) has emerged as key enzyme in many cardiac pathologies, especially heart failure (HF), myocardial infarction and cardiomyopathies, thus leading to contractile dysfunction and malignant arrhythmias. While many pathways leading to CaMKII activation have been elucidated in recent years, hardly any clinically viable compounds affecting CaMKII activity have progressed from basic in vitro science to in vivo studies. This review focuses on recent advances in anti-arrhythmic strategies involving CaMKII. Specifically, both inhibition of CaMKII itself to prevent arrhythmias, as well as anti-arrhythmic approaches affecting CaMKII activity via alterations in signaling cascades upstream and downstream of CaMKII will be discussed.

Sarcolemmal Ca(2+)-entry through L-type Ca(2+) channels controls the profile of Ca(2+)-activated Cl(-) current in canine ventricular myocytes

Ca(2+)-activated Cl(-) current (ICl(Ca)) mediated by TMEM16A and/or Bestrophin-3 may contribute to cardiac arrhythmias. The true profile of ICl(Ca) during an actual ventricular action potential (AP), however, is poorly understood. We aimed to study the profile of ICl(Ca) systematically under physiological conditions (normal Ca(2+) cycling and AP voltage-clamp) as well as in conditions designed to change [Ca(2+)]i. The expression of TMEM16A and/or Bestrophin-3 in canine and human left ventricular myocytes was examined. The possible spatial distribution of these proteins and their co-localization with Cav1.2 was also studied. The profile of ICl(Ca), identified as a 9-anthracene carboxylic acid-sensitive current under AP voltage-clamp conditions, contained an early fast outward and a late inward component, overlapping early and terminal repolarizations, respectively. Both components were moderately reduced by ryanodine, while fully abolished by BAPTA, but not EGTA. [Ca(2+)]i was monitored using Fura-2-AM. Setting [Ca(2+)]i to the systolic level measured in the bulk cytoplasm (1.1μM) decreased ICl(Ca), while application of Bay K8644, isoproterenol, and faster stimulation rates increased the amplitude of ICl(Ca). Ca(2+)-entry through L-type Ca(2+) channels was essential for activation of ICl(Ca). TMEM16A and Bestrophin-3 showed strong co-localization with one another and also with Cav1.2 channels, when assessed using immunolabeling and confocal microscopy in both canine myocytes and human ventricular myocardium. Activation of ICl(Ca) in canine ventricular cells requires Ca(2+)-entry through neighboring L-type Ca(2+) channels and is only augmented by SR Ca(2+)-release. Substantial activation of ICl(Ca) requires high Ca(2+) concentration in the dyadic clefts which can be effectively buffered by BAPTA, but not EGTA.

Mechanistic Investigation of the Arrhythmogenic Role of Oxidized CaMKII in the Heart

Oxidative stress and calcium (Ca(2+))/calmodulin (CaM)-dependent protein kinase II (CaMKII) both play important roles in the pathogenesis of cardiac disease. Although the pathophysiological relevance of reactive oxygen species (ROS) and CaMKII has been appreciated for some time, recent work has shown that ROS can directly oxidize CaMKII, leading to its persistent activity and an increase of the likelihood of cellular arrhythmias such as early afterdepolarizations (EADs). Because CaMKII modulates the function of many proteins involved in excitation-contraction coupling, elucidation of its role in cardiac function, in both healthy and oxidative stress conditions, is challenging. To investigate this role, we have developed a model of CaMKII activation that includes both the phosphorylation-dependent and the newly identified oxidation-dependent activation pathways. This model is incorporated into our previous local-control model of the cardiac myocyte that describes excitation-contraction coupling via stochastic simulation of individual Ca(2+) release units and CaMKII-mediated phosphorylation of L-type Ca(2+) channels (LCCs), ryanodine receptors and sodium (Na(+)) channels. The model predicts the experimentally measured slow-rate dependence of H2O2-induced EADs. Upon increased H2O2, simulations suggest that selective activation of late Na(+) current (INaL), although it prolongs action potential duration, is not by itself sufficient to produce EADs. Similar results are obtained if CaMKII effects on LCCs and ryanodine receptors are considered separately. However, EADs emerge upon simultaneous activation of both LCCs and Na(+) channels. Further modeling results implicate activation of the Na(+)-Ca(2+) exchanger (NCX) as an important player in the generation of EADs. During bradycardia, the emergence of H2O2-induced EADs was correlated with a shift in the timing of NCX current reversal toward the plateau phase earlier in the action potential. Using the timing of NCX current reversal as an indicator event for EADs, the model identified counterintuitive ionic changes-difficult to experimentally dissect-that have the greatest influence on ROS-related arrhythmia propensity.

Filament Dynamics during Simulated Ventricular Fibrillation in a High-Resolution Rabbit Heart

The mechanisms underlying ventricular fibrillation (VF) are not well understood. The electrical activity on the heart surface during VF has been recorded extensively in the experimental setting and in some cases clinically; however, corresponding transmural activation patterns are prohibitively difficult to measure. In this paper, we use a high-resolution biventricular heart model to study three-dimensional electrical activity during fibrillation, focusing on the driving sources of VF: “filaments,” the organising centres of unstable reentrant scroll waves. We show, for the first time, specific 3D filament dynamics during simulated VF in a whole heart geometry that includes fine-scale anatomical structures. Our results suggest that transmural activity is much more complex than what would be expected from surface observations alone. We present examples of complex intramural activity, including filament breakup and reattachment, anchoring to the thin right ventricular apex; rapid transitions among various filament shapes; and filament lengths much greater than wall thickness. We also present evidence for anatomy playing a major role in VF development and coronary vessels and trabeculae influencing filament dynamics. Overall, our results indicate that intramural activity during simulated VF is extraordinarily complex and suggest that further investigation of 3D filaments is necessary to fully comprehend recorded surface patterns.

OpenCOR: a modular and interoperable approach to computational biology

Computational biologists have been developing standards and formats for nearly two decades, with the aim of easing the description and exchange of experimental data, mathematical models, simulation experiments, etc. One of those efforts is CellML (cellml.org), an XML-based markup language for the encoding of mathematical models. Early CellML-based environments include COR and OpenCell. However, both of those tools have limitations and were eventually replaced with OpenCOR (opencor.ws). OpenCOR is an open source modeling environment that is supported on Windows, Linux and OS X. It relies on a modular approach, which means that all of its features come in the form of plugins. Those plugins can be used to organize, edit, simulate and analyze models encoded in the CellML format. We start with an introduction to CellML and two of its early adopters, which limitations eventually led to the development of OpenCOR. We then go onto describing the general philosophy behind OpenCOR, as well as describing its openness and its development process. Next, we illustrate various aspects of OpenCOR, such as its user interface and some of the plugins that come bundled with it (e.g., its editing and simulation plugins). Finally, we discuss some of the advantages and limitations of OpenCOR before drawing some concluding remarks.

New antiarrhythmic targets to control intracellular calcium handling

Sudden cardiac death due to ventricular arrhythmias is a major problem. Drug therapies to prevent SCD do not provide satisfying results, leading to the demand for new antiarrhythmic strategies. New targets include Ca²⁺/Calmodulin-dependent protein kinase II (CaMKII), the Na/Ca exchanger (NCX), the Ryanodine receptor (RyR, and its associated protein FKBP12.6 (Calstabin)) and the late component of the sodium current (I_Na-Late), all related to intracellular calcium (Ca²⁺) handling. In this review, drugs interfering with these targets (SEA-0400, K201, KN-93, W7, ranolazine, sophocarpine, and GS-967) are evaluated and their future as clinical compounds is considered. These new targets prove to be interesting; however more insight into long-term drug effects is necessary before clinical applicability becomes reality.

Role of CaMKII in post acidosis arrhythmias: a simulation study using a human myocyte model

Postacidotic arrhythmias have been associated to increased sarcoplasmic reticulum (SR) Ca(2+) load and Ca(2+)/calmodulin-dependent protein kinase II (CaMKII) activation. However, the molecular mechanisms underlying these arrhythmias are still unclear. To better understand this process, acidosis produced by CO2 increase from 5% to 30%, resulting in intracellular pH (pHi) change from 7.15 to 6.7, was incorporated into a myocyte model of excitation-contraction coupling and contractility, including acidotic inhibition of L-type Ca(2+) channel (I(CaL)), Na(+)-Ca(2+) exchanger, Ca(2+) release through the SR ryanodine receptor (RyR2) (I(rel)), Ca(2+) reuptake by the SR Ca(2+) ATPase2a (I(up)), Na(+)-K(+) pump, K(+) efflux through the inward rectifier K(+) channel and the transient outward K(+) flow (I(to)) together with increased activity of the Na(+)-H(+) exchanger (I(NHE)). Simulated CaMKII regulation affecting I(rel), I(up), I(CaL), I(NHE) and I(to) was introduced in the model to partially compensate the acidosis outcome. Late Na(+) current increase by CaMKII was also incorporated. Using this scheme and assuming that diastolic Ca(2+) leak through the RyR2 was modulated by the resting state of this channel and the difference between SR and dyadic cleft [Ca(2+)], postacidotic delayed after depolarizations (DADs) were triggered upon returning to normal pHi after 6 min acidosis. The model showed that DADs depend on SR Ca(2+) load and on increased Ca(2+) leak through RyR2. This postacidotic arrhythmogenic pattern relies mainly on CaMKII effect on I(CaL) and I(up), since its individual elimination produced the highest DAD reduction. The model further revealed that during the return to normal pHi, DADs are fully determined by SR Ca(2+) load at the end of acidosis. Thereafter, DADs are maintained by SR Ca(2+) reloading by Ca(2+) influx through the reverse NCX mode during the time period in which [Na(+)]i is elevated.

Atrial fibrillation in the elderly: the potential contribution of reactive oxygen species

Atrial fibrillation (AF) is the most commonly encountered cardiac arrhythmia, and is a significant source of healthcare expenditures throughout the world. It is an arrhythmia with a very clearly defined predisposition for individuals of advanced age, and this fact has led to intense study of the mechanistic links between aging and AF. By promoting oxidative damage to multiple subcellular and cellular structures, reactive oxygen species (ROS) have been shown to induce the intra- and extra-cellular changes necessary to promote the pathogenesis of AF. In addition, the generation and accumulation of ROS have been intimately linked to the cellular processes which underlie aging. This review begins with an overview of AF pathophysiology, and introduces the critical structures which, when damaged, predispose an otherwise healthy atrium to AF. The available evidence that ROS can lead to damage of these critical structures is then reviewed. Finally, the evidence linking the process of aging to the pathogenesis of AF is discussed.

Computational approaches to understand cardiac electrophysiology and arrhythmias

Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.

How the Hodgkin-Huxley equations inspired the Cardiac Physiome Project

Early modelling of cardiac cells (1960-1980) was based on extensions of the Hodgkin-Huxley nerve axon equations with additional channels incorporated, but after 1980 it became clear that processes other than ion channel gating were also critical in generating electrical activity. This article reviews the development of models representing almost all cell types in the heart, many different species, and the software tools that have been created to facilitate the cardiac Physiome Project.