Simulation and mechanistic investigation of the arrhythmogenic role of the late sodium current in human heart failure

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Optimization of a cardiomyocyte model illuminates role of increased INa,L in repolarization reserve

Cardiac ion currents may compensate for each other when one is compromised by a congenital or drug-induced defect. Such redundancy contributes to a robust repolarization reserve that can prevent the development of lethal arrhythmias. Most efforts made to describe this phenomenon have quantified contributions by individual ion currents. However, it is important to understand the interplay between all major ion-channel conductances, as repolarization reserve is dependent on the balance between all ion currents in a cardiomyocyte. Here, a genetic algorithm was designed to derive profiles of nine ion-channel conductances that optimize repolarization reserve in a mathematical cardiomyocyte model. Repolarization reserve was quantified using a previously defined metric, repolarization reserve current, i.e., the minimum constant current to prevent normal action potential repolarization in a cell. The optimization improved repolarization reserve current up to 84% compared to baseline in a human adult ventricular myocyte model and increased resistance to arrhythmogenic insult. The optimized conductance profiles were not only characterized by increased repolarizing current conductances but also uncovered a previously unreported behavior by the late sodium current. Simulations demonstrated that upregulated late sodium increased action potential duration, without compromising repolarization reserve current. The finding was generalized to multiple models. Ultimately, this computational approach, in which multiple currents were studied simultaneously, illuminated mechanistic insights into how the metric's magnitude could be increased and allowed for the unexpected role of late sodium to be elucidated.NEW & NOTEWORTHY Genetic algorithms are typically used to fit models or extract desired parameters from data. Here, we use the tool to produce a ventricular cardiomyocyte model with increased repolarization reserve. Since arrhythmia mitigation is dependent on multiple cardiac ion-channel conductances, study using a comprehensive, unbiased, and systems-level approach is important. The use of this optimization strategy allowed us to find robust profiles that illuminated unexpected mechanistic determinants of key ion-channel conductances in repolarization reserve.

Discussion on Repolarization Reserve between Patients with Coronary Heart Disease and Normal Controls

Objective

To investigate the repolarization reserve of normal controls (NCs) and patients with coronary heart disease (CHD).

Methods

From January 1st, 2010 to December 31st, 2018, 200 age- and gender-matched inpatients in the Second Hospital of Shanxi Medical University and Heji Hospital Affiliated to Changzhi Medical College were selected for treadmill exercise test (TET), including 67 patients in the myocardial ischemia group, 66 patients in the suspected myocardial ischemia group, and 67 patients in the normal control group. Coronary angiography (CAG) was performed on 49 of 133 patients in the myocardial ischemia group and the suspected myocardial ischemia group, and 9 positives and 40 negatives were identified. The heart rate (HR) and QT interval of TET examiners before exercise, during exercise (90 beats/min, 120 beats/min, maximum HR), and in the recovery period (1 minute and 3 minutes after exercise) were reviewed, and QTc values were calculated after being corrected by BaZett's.

Results

The mean QTc values in NCs were all below 452 ms, before exercise, during exercise (90 beats/min, 120 beats/min and maximum HR), and during the recovery period (1 minute and 3 minutes after exercise). The comparison results of the RR interval between the two groups revealed P > 0.05, indicating no statistical significance. Significant differences were present when comparing the QT intervals when the HRs were 90 beats/minute and 120 beats/minute during exercise (P < 0.05). And comparing the QTc values, it was found that the QTc values during different exercise periods were statistically different between groups (P < 0.05).

Conclusions

NCs have good repolarization reserve. CAG can confirm true myocardial ischemia patients (i.e., patients with CHD) among myocardial ischemia and suspected myocardial ischemia patients screened by TET. Patients with positive CAG have poor repolarization reserve as QT interval represents ventricular repolarization adaptability.

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.

Late Sodium Current of the Heart: Where Do We Stand and Where Are We Going?

Late sodium current has long been linked to dysrhythmia and contractile malfunction in the heart. Despite the increasing body of accumulating information on the subject, our understanding of its role in normal or pathologic states is not complete. Even though the role of late sodium current in shaping action potential under physiologic circumstances is debated, it’s unquestioned role in arrhythmogenesis keeps it in the focus of research. Transgenic mouse models and isoform-specific pharmacological tools have proved useful in understanding the mechanism of late sodium current in health and disease. This review will outline the mechanism and function of cardiac late sodium current with special focus on the recent advances of the area.

Protective Effect of Trimetazidine on Potassium Ion Homeostasis in Myocardial Tissue in Mice with Heart Failure

The occurrence of heart failure (HF) is closely correlated with the disturbance of mitochondrial energy metabolism, and trimetazidine (TMZ) has been regarded as an effective agent in treating HF. Intracellular potassium ion (K⁺) homeostasis, which is modulated by K⁺ channels and transporters, is crucial for maintaining normal myocardial function and can be disrupted by HF. This study is aimed at exploring the protective effect of TMZ on K⁺ homeostasis within myocardial tissue in mice with HF. We observed the pathological changes of myocardial tissue under microscopes and further measured the content of adenosine triphosphate (ATP), the activity of Na⁺-K⁺ ATPase, and the expression of ATP1α1 at the mRNA and protein levels. Moreover, we also analyzed the changes in K⁺ flux across the myocardial tissue in mice. As a result, we found that there was a large amount of myocardial fiber lysis and fracture in HF myocardial tissue. Meanwhile, the potassium flux of mice with HF was reduced, and the expression of ATP1α1, the activity of Na⁺-K⁺ ATPase, and the supply and delivery of ATP were also decreased. In contrast, TMZ can effectively treat HF by restoring K⁺ homeostasis in the local microenvironment of myocardial tissues.

Effect of ivabradine in heart failure: a meta-analysis of heart failure patients with reduced versus preserved ejection fraction

In clinical trials of heart failure reduced ejection fraction (HFrEF), ivabradine seemed to be an effective heart rate lowering agent associated with lower risk of cardiovascular death. In contrast, ivabradine failed to improve cardiovascular outcomes in heart failure preserved ejection fraction (HFpEF) despite the significant effect on heart rate. This meta-analysis is the first to compare the effects of ivabradine on heart rate and mortality parameters in HFpEF versus HFrEF. We screened three databases: PubMed, Embase, and Cochrane Library. The outcomes of these studies were mortality, reduction in heart rate, and left ventricular function improvement. We compared the efficacy of ivabradine treatment in HFpEF versus HFrEF. Heart rate analysis of pooled data showed decrease in both HFrEF (-17.646 beats/min) and HFpEF (-11.434 beats/min), and a tendency to have stronger bradycardic effect in HFrEF (p = 0.094) in randomized clinical trials. Left ventricular ejection fraction analysis revealed significant improvement in HFrEF (5.936, 95% CI: [4.199-7.672], p < 0.001) when compared with placebo (p < 0.001). We found that ivabradine significantly improves left ventricular performance in HFrEF, at the same time it exerts a tendency to have improved bradycardic effect in HFrEF. These disparate effects of ivabradine and the higher prevalence of non-cardiac comorbidities in HFpEF may explain the observed beneficial effects in HFrEF and the unchanged outcomes in HFpEF patients after ivabradine treatment.

Ranolazine as an Alternative Therapy to Flecainide for SCN5A V411M Long QT Syndrome Type 3 Patients

The prolongation of the QT interval represents the main feature of the long QT syndrome (LQTS), a life-threatening genetic disease. The heterozygous SCN5A V411M mutation of the human sodium channel leads to a LQTS type 3 with severe proarrhythmic effects due to an increase in the late component of the sodium current (INaL). The two sodium blockers flecainide and ranolazine are equally recommended by the current 2015 ESC guidelines to treat patients with LQTS type 3 and persistently prolonged QT intervals. However, awareness of pro-arrhythmic effects of flecainide in LQTS type 3 patients arose upon the study of the SCN5A E1784K mutation. Regarding SCN5A V411M individuals, flecainide showed good results albeit in a reduced number of patients and no evidence supporting the use of ranolazine has ever been released. Therefore, we ought to compare the effect of ranolazine and flecainide in a SCN5A V411M model using an in-silico modeling and simulation approach. We collected clinical data of four patients. Then, we fitted four Markovian models of the human sodium current (INa) to experimental and clinical data. Two of them correspond to the wild type and the heterozygous SCN5A V411M scenarios, and the other two mimic the effects of flecainide and ranolazine on INa. Next, we inserted them into three isolated cell action potential (AP) models for endocardial, midmyocardial and epicardial cells and in a one-dimensional tissue model. The SCN5A V411M mutation produced a 15.9% APD90 prolongation in the isolated endocardial cell model, which corresponded to a 14.3% of the QT interval prolongation in a one-dimensional strand model, in keeping with clinical observations. Although with different underlying mechanisms, flecainide and ranolazine partially countered this prolongation at the isolated endocardial model by reducing the APD90 by 8.7 and 4.3%, and the QT interval by 7.2 and 3.2%, respectively. While flecainide specifically targeted the mutation-induced increase in peak INaL, ranolazine reduced it during the entire AP. Our simulations also suggest that ranolazine could prevent early afterdepolarizations triggered by the SCN5A V411M mutation during bradycardia, as flecainide. We conclude that ranolazine could be used to treat SCN5A V411M patients, specifically when flecainide is contraindicated.

Late INa Blocker GS967 Supresses Polymorphic Ventricular Tachycardia in a Transgenic Rabbit Model of Long QT Type 2

BACKGROUND

Long QT syndrome has been associated with sudden cardiac death likely caused by early afterdepolarizations (EADs) and polymorphic ventricular tachycardias (PVTs). Suppressing the late sodium current (INaL) may counterbalance the reduced repolarization reserve in long QT syndrome and prevent EADs and PVTs.

METHODS

We tested the effects of the selective INaL blocker GS967 on PVT induction in a transgenic rabbit model of long QT syndrome type 2 using intact heart optical mapping, cellular electrophysiology and confocal Ca2+ imaging, and computer modeling.

RESULTS

GS967 reduced ventricular fibrillation induction under a rapid pacing protocol (n=7/14 hearts in control versus 1/14 hearts at 100 nmol/L) without altering action potential duration or restitution and dispersion. GS967 suppressed PVT incidences by reducing Ca2+-mediated EADs and focal activity during isoproterenol perfusion (at 30 nmol/L, n=7/12 and 100 nmol/L n=8/12 hearts without EADs and PVTs). Confocal Ca2+ imaging of long QT syndrome type 2 myocytes revealed that GS967 shortened Ca2+ transient duration via accelerating Na+/Ca2+ exchanger (INCX)-mediated Ca2+ efflux from cytosol, thereby reducing EADs. Computer modeling revealed that INaL potentiates EADs in the long QT syndrome type 2 setting through (1) providing additional depolarizing currents during action potential plateau phase, (2) increasing intracellular Na+ (Nai) that decreases the depolarizing INCX thereby suppressing the action potential plateau and delaying the activation of slowly activating delayed rectifier K+ channels (IKs), suggesting important roles of INaL in regulating Nai.

CONCLUSIONS

Selective INaL blockade by GS967 prevents EADs and abolishes PVT in long QT syndrome type 2 rabbits by counterbalancing the reduced repolarization reserve and normalizing Nai. Graphic Abstract: A graphic abstract is available for this article.

Quantitative analysis of variability in an integrated model of human ventricular electrophysiology and β-adrenergic signaling

In ventricular myocytes, stimulation of β-adrenergic receptors activates critical cardiac signaling pathways, leading to shorter action potentials and increased contraction strength during the "fight-or-flight" response. These changes primarily result, at the cellular level, from the coordinated phosphorylation of multiple targets by protein kinase A. Although mathematical models of the intracellular signaling downstream of β-adrenergic receptor activation have previously been described, only a limited number of studies have explored quantitative interactions between intracellular signaling and electrophysiology in human ventricular myocytes. Accordingly, our objective was to develop an integrative mathematical model of β-adrenergic receptor signaling, electrophysiology, and intracellular calcium (Ca²⁺) handling in the healthy human ventricular myocyte. We combined published mathematical models of intracellular signaling and electrophysiology, then calibrated the model results against voltage clamp data and physiological changes occurring after stimulation of β-adrenergic receptors with isoproterenol. We subsequently: (1) explored how molecular variability in different categories of model parameters translated into phenotypic variability; (2) identified the most important parameters determining physiological cellular outputs in the model before and after β-adrenergic receptor stimulation; and (3) investigated which molecular level alterations can produce a phenotype indicative of heart failure with preserved ejection fraction (HFpEF). Major results included: (1) variability in parameters that controlled intracellular signaling caused qualitatively different behavior than variability in parameters controlling ion transport pathways; (2) the most important model parameters determining action potential duration and intracellular Ca²⁺ transient amplitude were generally consistent before and after β-adrenergic receptor stimulation, except for a shift in the importance of K⁺ currents in determining action potential duration; and (3) decreased Ca²⁺ uptake into the sarcoplasmic reticulum, increased Ca²⁺ extrusion through Na⁺/Ca²⁺ exchanger and decreased Ca²⁺ leak from the sarcoplasmic reticulum may contribute to HFpEF. Overall, this study provided novel insight into the phenotypic consequences of molecular variability, and our integrated model may be useful in the design and interpretation of future experimental studies of interactions between β-adrenergic signaling and cellular physiology in human ventricular myocytes.

Late Sodium Current Inhibitors as Potential Antiarrhythmic Agents

Based on recent findings, an increased late sodium current (INa,late) plays an important pathophysiological role in cardiac diseases, including rhythm disorders. The article first describes what is INa,late and how it functions under physiological circumstances. Next, it shows the wide range of cellular mechanisms that can contribute to an increased INa,late in heart diseases, and also discusses how the upregulated INa,late can play a role in the generation of cardiac arrhythmias. The last part of the article is about INa,late inhibiting drugs as potential antiarrhythmic agents, based on experimental and preclinical data as well as in the light of clinical trials.

Impact of disease state on arrhythmic event detection by action potential modelling in cardiac safety pharmacology

INTRODUCTION

The use of in silico cardiac action potential simulations is one of the pillars of the CiPA initiative (Comprehensive in vitro Proarrhythmia Assay) currently under evaluation designed to detect more accurately proarrhythmic liabilities of new drug candidate. In order to take into account the variability of clinical situations, we propose to improve this method by studying the impact of various disease states on arrhythmic events induced by 30 torsadogenic or non-torsadogenic compounds.

METHOD

In silico modelling was done on the human myocytes using the Dutta revised O'Hara-Rudy algorithm. Results were analysed using a new metric based on the compound IC_50s against the seven cardiac ionic currents considered to be the most important by the CiPA initiative (I_Kr, I_Ks, I_Na, I_NaL, I_K1, I_to, I_CaL) and the minimal rate of action potential voltage decrease calculated at the early-afterdepolarization (EAD) take-off membrane voltage (V_min).

RESULTS

The specific threshold at which each torsadogenic compounds induced EAD, was exacerbated by the presence of cardiac risk factors ranked as follows: congestive heart failure > hypertrophic cardiomyopathy > cardiac pause > no risk factor. Non-torsadogenic compounds induced no EAD even in the presence of cardiac risk factors.

DISCUSSION

The present study highlighted the impact of pre-existing cardiovascular disease on arrhythmic event detection suggesting that disease state modelling may need to be incorporated in order to fully realize the goal of the CiPA paradigm in a more accurate predictability of proarrhythmic liabilities of new drug candidate.

Ca2+ Cycling Impairment in Heart Failure Is Exacerbated by Fibrosis: Insights Gained From Mechanistic Simulations

Heart failure (HF) is characterized by altered Ca²⁺ cycling, resulting in cardiac contractile dysfunction. Failing myocytes undergo electrophysiological remodeling, which is known to be the main cause of abnormal Ca²⁺ homeostasis. However, structural remodeling, specifically proliferating fibroblasts coupled to myocytes in the failing heart, could also contribute to Ca²⁺ cycling impairment. The goal of the present study was to systematically analyze the mechanisms by which myocyte–fibroblast coupling could affect Ca²⁺ dynamics in normal conditions and in HF. Simulations of healthy and failing human myocytes were performed using established mathematical models, and cells were either isolated or coupled to fibroblasts. Univariate and multivariate sensitivity analyses were performed to quantify effects of ion transport pathways on biomarkers computed from intracellular [Ca²⁺] waveforms. Variability in ion channels and pumps was imposed and populations of models were analyzed to determine effects on Ca²⁺ dynamics. Our results suggest that both univariate and multivariate sensitivity analyses are valuable methodologies to shed light into the ionic mechanisms underlying Ca²⁺ impairment in HF, although differences between the two methodologies are observed at high parameter variability. These can result from either the fact that multivariate analyses take into account ion channels or non-linear effects of ion transport pathways on Ca²⁺ dynamics. Coupling either healthy or failing myocytes to fibroblasts decreased Ca²⁺ transients due to an indirect sink effect on action potential (AP) and thus on Ca²⁺ related currents. Simulations that investigated restoration of normal physiology in failing myocytes showed that Ca²⁺ cycling can be normalized by increasing SERCA and L-type Ca²⁺ current activity while decreasing Na⁺–Ca²⁺ exchange and SR Ca²⁺ leak. Changes required to normalize APs in failing myocytes depended on whether myocytes were coupled to fibroblasts. In conclusion, univariate and multivariate sensitivity analyses are helpful tools to understand how Ca²⁺ cycling is impaired in HF and how this can be exacerbated by coupling of myocytes to fibroblasts. The design of pharmacological actions to restore normal activity should take into account the degree of fibrosis in the failing heart.

Paradoxical Onset of Arrhythmic Waves from Depolarized Areas in Cardiac Tissue Due to Curvature-Dependent Instability

The generation of abnormal excitations in pathological regions of the heart is a main trigger for lethal cardiac arrhythmias. Such abnormal excitations, also called ectopic activity, often arise from areas with local tissue heterogeneity or damage accompanied by localized depolarization. Finding the conditions that lead to ectopy is important to understand the basic biophysical principles underlying arrhythmia initiation and might further refine clinical procedures. In this study, we are the first to address the question of how geometry of the abnormal region affects the onset of ectopy using a combination of experimental, in silico, and theoretical approaches. We paradoxically find that, for any studied geometry of the depolarized region in optogenetically modified monolayers of cardiac cells, primary ectopic excitation originates at areas of maximal curvature of the boundary, where the stimulating electrotonic currents are minimal. It contradicts the standard critical nucleation theory applied to nonlinear waves in reaction-diffusion systems, where a higher stimulus is expected to produce excitation more easily. Our in silico studies reveal that the nonconventional ectopic activity is caused by an oscillatory instability at the boundary of the damaged region, the occurrence of which depends on the curvature of that boundary. The onset of this instability is confirmed using the Schrödinger equation methodology proposed by Rinzel and Keener [SIAM J. Appl. Math. 43, 907 (1983)]. Overall, we show distinctively novel insight into how the geometry of a heterogeneous cardiac region determines ectopic activity, which can be used in the future to predict the conditions that can trigger cardiac arrhythmias.

Pro-arrhythmic effects of low plasma [K+] in human ventricle: An illustrated review

Potassium levels in the plasma, [K⁺]_o, are regulated precisely under physiological conditions. However, increases (from approx. 4.5 to 8.0mM) can occur as a consequence of, e.g., endurance exercise, ischemic insult or kidney failure. This hyperkalemic modulation of ventricular electrophysiology has been studied extensively. Hypokalemia is also common. It can occur in response to diuretic therapy, following renal dialysis, or during recovery from endurance exercise. In the human ventricle, clinical hypokalemia (e.g., [K⁺]_o levels of approx. 3.0mM) can cause marked changes in both the resting potential and the action potential waveform, and these may promote arrhythmias. Here, we provide essential background information concerning the main K⁺-sensitive ion channel mechanisms that act in concert to produce prominent short-term ventricular electrophysiological changes, and illustrate these by implementing recent mathematical models of the human ventricular action potential. Even small changes (~1mM) in [K⁺]_o result in significant alterations in two different K⁺ currents, I_K1 and HERG. These changes can markedly alter in resting membrane potential and/or action potential waveform in human ventricle. Specifically, a reduction in net outward transmembrane K⁺ currents (repolarization reserve) and an increased substrate input resistance contribute to electrophysiological instability during the plateau of the action potential and may promote pro-arrhythmic early after-depolarizations (EADs). Translational settings where these insights apply include: optimal diuretic therapy, and the interpretation of data from Phase II and III trials for anti-arrhythmic drug candidates.

Sensitivity analysis revealing the effect of modulating ionic mechanisms on calcium dynamics in simulated human heart failure

Abnormal intracellular Ca²⁺ handling is the major contributor to the depressed cardiac contractility observed in heart failure. The electrophysiological remodeling associated with this pathology alters both the action potential and the Ca²⁺ dynamics, leading to a defective excitation-contraction coupling that ends in mechanical dysfunction. The importance of maintaining a correct intracellular Ca²⁺ concentration requires a better understanding of its regulation by ionic mechanisms. To study the electrical activity and ionic homeostasis of failing myocytes, a modified version of the O’Hara et al. human action potential model was used, including electrophysiological remodeling. The impact of the main ionic transport mechanisms was analyzed using single-parameter sensitivity analyses, the first of which explored the modulation of electrophysiological characteristics related to Ca²⁺ exerted by the remodeled parameters. The second sensitivity analysis compared the potential consequences of modulating individual channel conductivities, as one of the main effects of potential drugs, on Ca²⁺ dynamic properties under both normal conditions and in heart failure. The first analysis revealed the important contribution of the sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) dysfunction to the altered Ca²⁺ homeostasis, with the Na⁺/Ca²⁺ exchanger (NCX) and other Ca²⁺ cycling proteins also playing a significant role. Our results highlight the importance of improving the SR uptake function to increase Ca²⁺ content and restore Ca²⁺ homeostasis and contractility. The second sensitivity analysis highlights the different response of the failing myocyte versus the healthy myocyte to potential pharmacological actions on single channels. The result of modifying the conductances of the remodeled proteins such as SERCA and NCX in heart failure has less impact on Ca²⁺ modulation. These differences should be taken into account when designing drug therapies.

Sensitivity analysis of ventricular activation and electrocardiogram in tailored models of heart-failure patients

Cardiac resynchronization therapy is not effective in a variable proportion of heart failure patients. An accurate knowledge of each patient's electroanatomical features could be helpful to determine the most appropriate treatment. The goal of this study was to analyze and quantify the sensitivity of left ventricular (LV) activation and the electrocardiogram (ECG) to changes in 39 parameters used to tune realistic anatomical-electrophysiological models of the heart. Electrical activity in the ventricles was simulated using a reaction-diffusion equation. To simulate cellular electrophysiology, the Ten Tusscher-Panfilov 2006 model was used. Intracardiac electrograms and 12-lead ECGs were computed by solving the bidomain equation. Parameters showing the highest sensitivity values were similar in the six patients studied. QRS complex and LV activation times were modulated by the sodium current, the cell surface-to-volume ratio in the LV, and tissue conductivities. The T-wave was modulated by the calcium and rectifier-potassium currents, and the cell surface-to-volume ratio in both ventricles. We conclude that homogeneous changes in ionic currents entail similar effects in all ECG leads, whereas the effects of changes in tissue properties show larger inter-lead variability. The effects of parameter variations are highly consistent between patients and most of the model tuning could be performed with only ~10 parameters.

Illuminating Myocyte-Fibroblast Homotypic and Heterotypic Gap Junction Dynamics Using Dynamic Clamp

Fibroblasts play a significant role in the development of electrical and mechanical dysfunction of the heart; however, the underlying mechanisms are only partially understood. One widely studied mechanism suggests that fibroblasts produce excess extracellular matrix, resulting in collagenous septa that slow propagation, cause zig-zag conduction paths, and decouple cardiomyocytes, resulting in a substrate for cardiac arrhythmia. An emerging mechanism suggests that fibroblasts promote arrhythmogenesis through direct electrical interactions with cardiomyocytes via gap junction (GJ) channels. In the heart, three major connexin (Cx) isoforms, Cx40, Cx43, and Cx45, form GJ channels in cell-type-specific combinations. Because each Cx is characterized by a unique time- and transjunctional voltage-dependent profile, we investigated whether the electrophysiological contributions of fibroblasts would vary with the specific composition of the myocyte-fibroblast (M-F) GJ channel. Due to the challenges of systematically modifying Cxs in vitro, we coupled native cardiomyocytes with in silico fibroblast and GJ channel electrophysiology models using the dynamic-clamp technique. We found that there is a reduction in the early peak of the junctional current during the upstroke of the action potential (AP) due to GJ channel gating. However, effects on the cardiomyocyte AP morphology were similar regardless of the specific type of GJ channel (homotypic Cx43 and Cx45, and heterotypic Cx43/Cx45 and Cx45/Cx43). To illuminate effects at the tissue level, we performed multiscale simulations of M-F coupling. First, we developed a cell-specific model of our dynamic-clamp experiments and investigated changes in the underlying membrane currents during M-F coupling. Second, we performed two-dimensional tissue sheet simulations of cardiac fibrosis and incorporated GJ channels in a cell type-specific manner. We determined that although GJ channel gating reduces junctional current, it does not significantly alter conduction velocity during cardiac fibrosis relative to static GJ coupling. These findings shed more light on the complex electrophysiological interplay between cardiac fibroblasts and myocytes.

A mathematical model for active contraction in healthy and failing myocytes and left ventricles

Cardiovascular disease is one of the leading causes of death worldwide, in particular myocardial dysfunction, which may lead to heart failure eventually. Understanding the electro-mechanics of the heart will help in developing more effective clinical treatments. In this paper, we present a multi-scale electro-mechanics model of the left ventricle (LV). The Holzapfel-Ogden constitutive law was used to describe the passive myocardial response in tissue level, a modified Grandi-Pasqualini-Bers model was adopted to model calcium dynamics in individual myocytes, and the active tension was described using the Niederer-Hunter-Smith myofilament model. We first studied the electro-mechanics coupling in a single myocyte in the healthy and diseased left ventricle, and then the single cell model was embedded in a dynamic LV model to investigate the compensation mechanism of LV pump function due to myocardial dysfunction caused by abnormality in cellular calcium dynamics. The multi-scale LV model was solved using an in-house developed hybrid immersed boundary method with finite element extension. The predictions of the healthy LV model agreed well with the clinical measurements and other studies, and likewise, the results in the failing states were also consistent with clinical observations. In particular, we found that a low level of intracellular Ca2+ transient in myocytes can result in LV pump function failure even with increased myocardial contractility, decreased systolic blood pressure, and increased diastolic filling pressure, even though they will increase LV stroke volume. Our work suggested that treatments targeted at increased contractility and lowering the systolic blood pressure alone are not sufficient in preventing LV pump dysfunction, restoring a balanced physiological Ca2+ handling mechanism is necessary.

Changes in Intracellular Na+ following Enhancement of Late Na+ Current in Virtual Human Ventricular Myocytes

The slowly inactivating or late Na+ current, INa-L, can contribute to the initiation of both atrial and ventricular rhythm disturbances in the human heart. However, the cellular and molecular mechanisms that underlie these pro-arrhythmic influences are not fully understood. At present, the major working hypothesis is that the Na+ influx corresponding to INa-L significantly increases intracellular Na+, [Na+]i; and the resulting reduction in the electrochemical driving force for Na+ reduces and (may reverse) Na+/Ca2+ exchange. These changes increase intracellular Ca2+, [Ca2+]i; which may further enhance INa-L due to calmodulin-dependent phosphorylation of the Na+ channels. This paper is based on mathematical simulations using the O'Hara et al (2011) model of baseline or healthy human ventricular action potential waveforms(s) and its [Ca2+]i homeostasis mechanisms. Somewhat surprisingly, our results reveal only very small changes (≤ 1.5 mM) in [Na+]i even when INa-L is increased 5-fold and steady-state stimulation rate is approximately 2 times the normal human heart rate (i.e. 2 Hz). Previous work done using well-established models of the rabbit and human ventricular action potential in heart failure settings also reported little or no change in [Na+]i when INa-L was increased. Based on our simulations, the major short-term effect of markedly augmenting INa-L is a significant prolongation of the action potential and an associated increase in the likelihood of reactivation of the L-type Ca2+ current, ICa-L. Furthermore, this action potential prolongation does not contribute to [Na+]i increase.

Current-Voltage Relationship for Late Na(+) Current in Adult Rat Ventricular Myocytes

It is now well established that the slowly inactivating component of the Na(+) current (INa-L) in the mammalian heart is a significant regulator of the action potential waveform. This insight has led to detailed studies of the role of INa-L in a number of important and challenging pathophysiological settings. These include genetically based ventricular arrhythmias (LQT 1, 2, and 3), ventricular arrhythmias arising from progressive cardiomyopathies (including diabetic), and proarrhythmic abnormalities that develop during local or global ventricular ischemia. Inhibition of INa-L may also be a useful strategy for management of atrial flutter and fibrillation. Many important biophysical parameters that characterize INa-L have been identified; and INa-L as an antiarrhythmia drug target has been studied extensively. However, relatively little information is available regarding (1) the ion transfer or current-voltage relationship for INa-L or (2) the time course of its reactivation at membrane potentials similar to the resting or diastolic membrane potential in mammalian ventricle. This chapter is based on our preliminary findings concerning these two very important physiological/biophysical descriptors for INa-L. Our results were obtained using whole-cell voltage clamp methods applied to enzymatically isolated rat ventricular myocytes. A chemical agent, BDF 9148, which was once considered to be a drug candidate in the Na(+)-dependent inotropic agent category has been used to markedly enhance INa-L current. BDF acts in a potent, selective, and reversible fashion. These BDF 9148 effects are compared and contrasted with the prototypical activator of INa-L, a sea anemone toxin, ATX II.

Reduced response to IKr blockade and altered hERG1a/1b stoichiometry in human heart failure

Heart failure (HF) claims 250,000 lives per year in the US, and nearly half of these deaths are sudden and presumably due to ventricular tachyarrhythmias. QT interval and action potential (AP) prolongation are hallmark proarrhythmic changes in the failing myocardium, which potentially result from alterations in repolarizing potassium currents. Thus, we aimed to examine whether decreased expression of the rapid delayed rectifier potassium current, IKr, contributes to repolarization abnormalities in human HF. To map functional IKr expression across the left ventricle (LV), we optically imaged coronary-perfused LV free wall from donor and end-stage failing human hearts. The LV wedge preparation was used to examine transmural AP durations at 80% repolarization (APD80), and treatment with the IKr-blocking drug, E-4031, was utilized to interrogate functional expression. We assessed the percent change in APD80 post-IKr blockade relative to baseline APD80 (∆APD80) and found that ∆APD80s are reduced in failing versus donor hearts in each transmural region, with 0.35-, 0.43-, and 0.41-fold reductions in endo-, mid-, and epicardium, respectively (p=0.008, 0.037, and 0.022). We then assessed hERG1 isoform gene and protein expression levels using qPCR and Western blot. While we did not observe differences in hERG1a or hERG1b gene expression between donor and failing hearts, we found a shift in the hERG1a:hERG1b isoform stoichiometry at the protein level. Computer simulations were then conducted to assess IKr block under E-4031 influence in failing and nonfailing conditions. Our results confirmed the experimental observations and E-4031-induced relative APD80 prolongation was greater in normal conditions than in failing conditions, provided that the cellular model of HF included a significant downregulation of IKr. In human HF, the response to IKr blockade is reduced, suggesting decreased functional IKr expression. This attenuated functional response is associated with altered hERG1a:hERG1b protein stoichiometry in the failing human LV, and failing cardiomyoctye simulations support the experimental findings. Thus, of IKr protein and functional expression may be important determinants of repolarization remodeling in the failing human LV.

An electromechanical left ventricular wedge model to study the effects of deformation on repolarization during heart failure

Heart failure is a major and costly problem in public health, which, in certain cases, may lead to death. The failing heart undergo a series of electrical and structural changes that provide the underlying basis for disturbances like arrhythmias. Computer models of coupled electrical and mechanical activities of the heart can be used to advance our understanding of the complex feedback mechanisms involved. In this context, there is a lack of studies that consider heart failure remodeling using strongly coupled electromechanics. We present a strongly coupled electromechanical model to study the effects of deformation on a human left ventricle wedge considering normal and hypertrophic heart failure conditions. We demonstrate through a series of simulations that when a strongly coupled electromechanical model is used, deformation results in the thickening of the ventricular wall that in turn increases transmural dispersion of repolarization. These effects were analyzed in both normal and failing heart conditions. We also present transmural electrograms obtained from these simulations. Our results suggest that the waveform of electrograms, particularly the T-wave, is influenced by cardiac contraction on both normal and pathological conditions.

Lessons learned from multi-scale modeling of the failing heart

Heart failure constitutes a major public health problem worldwide. Affected patients experience a number of changes in the electrical function of the heart that predispose to potentially lethal cardiac arrhythmias. Due to the multitude of electrophysiological changes that may occur during heart failure, the scientific literature is complex and sometimes ambiguous, perhaps because these findings are highly dependent on the etiology, the stage of heart failure, and the experimental model used to study these changes. Nevertheless, a number of common features of failing hearts have been documented. Prolongation of the action potential (AP) involving ion channel remodeling and alterations in calcium handling have been established as the hallmark characteristics of myocytes isolated from failing hearts. Intercellular uncoupling and fibrosis are identified as major arrhythmogenic factors. Multi-scale computational simulations are a powerful tool that complements experimental and clinical research. The development of biophysically detailed computer models of single myocytes and cardiac tissues has contributed greatly to our understanding of processes underlying excitation and repolarization in the heart. The electrical, structural, and metabolic remodeling that arises in cardiac tissues during heart failure has been addressed from different computational perspectives to further understand the arrhythmogenic substrate. This review summarizes the contributions from computational modeling and simulation to predict the underlying mechanisms of heart failure phenotypes and their implications for arrhythmogenesis, ranging from the cellular level to whole-heart simulations. The main aspects of heart failure are presented in several related sections. An overview of the main electrophysiological and structural changes that have been observed experimentally in failing hearts is followed by the description and discussion of the simulation work in this field at the cellular level, and then in 2D and 3D cardiac structures. The implications for arrhythmogenesis in heart failure are also discussed including therapeutic measures, such as drug effects and cardiac resynchronization therapy. Finally, the future challenges in heart failure modeling and simulation will be discussed.

Representing variability and transmural differences in a model of human heart failure

During heart failure (HF) at the cellular level, the electrophysiological properties of single myocytes get remodeled, which can trigger the occurrence of ventricular arrhythmias that could be manifested in many forms such as early afterdepolarizations (EADs) and alternans (ALTs). In this paper, based on experimentally observed human HF data, specific ionic and exchanger current strengths are modified from a recently developed human ventricular cell model: the O'Hara-Virág-Varró-Rudy (OVVR) model. A new transmural HF-OVVR model is developed that incorporates HF changes and variability of the observed remodeling. This new heterogeneous HF-OVVR model is able to replicate many of the failing action potential (AP) properties and the dynamics of both [Ca(2+)]i and [Na(+)]i in accordance with experimental data. Moreover, it is able to generate EADs for different cell types and exhibits ALTs at modest pacing rate for transmural cell types. We have assessed the HF-OVVR model through the examination of the AP duration and the major ionic currents' rate dependence in single myocytes. The evaluation of the model comes from utilizing the steady-state (S-S) and S1-S2 restitution curves and from probing the accommodation of the HF-OVVR model to an abrupt change in cycle length. In addition, we have investigated the effect of chosen currents on the AP properties, such as blocking the slow sodium current to shorten the AP duration and suppress the EADs, and have found good agreement with experimental observations. This study should help elucidate arrhythmogenic mechanisms at the cellular level and predict unseen properties under HF conditions. In addition, this AP cell model might be useful for modeling and simulating HF at the tissue and organ levels.

In silico assessment of drug safety in human heart applied to late sodium current blockers

Drug-induced action potential (AP) prolongation leading to Torsade de Pointes is a major concern for the development of anti-arrhythmic drugs. Nevertheless the development of improved anti-arrhythmic agents, some of which may block different channels, remains an important opportunity. Partial block of the late sodium current (I(NaL)) has emerged as a novel anti-arrhythmic mechanism. It can be effective in the settings of free radical challenge or hypoxia. In addition, this approach can attenuate pro-arrhythmic effects of blocking the rapid delayed rectifying K(+) current (I(Kr)). The main goal of our computational work was to develop an in-silico tool for preclinical anti-arrhythmic drug safety assessment, by illustrating the impact of I(Kr)/I(NaL) ratio of steady-state block of drug candidates on "torsadogenic" biomarkers. The O'Hara et al. AP model for human ventricular myocytes was used. Biomarkers for arrhythmic risk, i.e., AP duration, triangulation, reverse rate-dependence, transmural dispersion of repolarization and electrocardiogram QT intervals, were calculated using single myocyte and one-dimensional strand simulations. Predetermined amounts of block of I(NaL) and I(Kr) were evaluated. "Safety plots" were developed to illustrate the value of the specific biomarker for selected combinations of IC(50)s for I(Kr) and I(NaL) of potential drugs. The reference biomarkers at baseline changed depending on the "drug" specificity for these two ion channel targets. Ranolazine and GS967 (a novel potent inhibitor of I(NaL)) yielded a biomarker data set that is considered safe by standard regulatory criteria. This novel in-silico approach is useful for evaluating pro-arrhythmic potential of drugs and drug candidates in the human ventricle.

Abnormal calcium homeostasis in heart failure with preserved ejection fraction is related to both reduced contractile function and incomplete relaxation: an electromechanically detailed biophysical modeling study

Heart failure with preserved ejection fraction (HFpEF) accounts for about 50% of heart failure cases. It has features of incomplete relaxation and increased stiffness of the left ventricle. Studies from clinical electrophysiology and animal experiments have found that HFpEF is associated with impaired calcium homeostasis, ion channel remodeling and concentric left ventricle hypertrophy (LVH). However, it is still unclear how the abnormal calcium homeostasis, ion channel and structural remodeling affect the electro-mechanical dynamics of the ventricles. In this study we have developed multiscale models of the human left ventricle from single cells to the 3D organ, which take into consideration HFpEF-induced changes in calcium handling, ion channel remodeling and concentric LVH. Our simulation results suggest that at the cellular level, HFpEF reduces the systolic calcium level resulting in a reduced systolic contractile force, but elevates the diastolic calcium level resulting in an abnormal residual diastolic force. In our simulations, these abnormal electro-mechanical features of the ventricular cells became more pronounced with the increase of the heart rate. However, at the 3D organ level, the ejection fraction of the left ventricle was maintained due to the concentric LVH. The simulation results of this study mirror clinically observed features of HFpEF and provide new insights toward the understanding of the cellular bases of impaired cardiac electromechanical functions in heart failure.

The late sodium current in heart failure: pathophysiology and clinical relevance

Large and growing body of data suggest that an increased late sodium current (I_Na,late ) can have a significant pathophysiological role in heart failure and other heart diseases. The first goal of this article is to describe how I_Na,late functions under physiological circumstances. The second goal is to show the wide range of cellular mechanisms that can increase I_Na,late in cardiac disease, and also to describe how the up-regulated I_Na,late contributes to the pathophysiology of heart failure. The final section of the article discusses the possible use of I_Na,late -modifying drugs in heart failure, on the basis of experimental and preclinical data.

Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 1D simulation study

Background

Heart failure is a final common pathway or descriptor for various cardiac pathologies. It is associated with sudden cardiac death, which is frequently caused by ventricular arrhythmias. Electrophysiological remodeling, intercellular uncoupling, fibrosis and autonomic imbalance have been identified as major arrhythmogenic factors in heart failure etiology and progression.

Objective

In this study we investigate in silico the role of electrophysiological and structural heart failure remodeling on the modulation of key elements of the arrhythmogenic substrate, i.e., electrophysiological gradients and abnormal impulse propagation.

Methods

Two different mathematical models of the human ventricular action potential were used to formulate models of the failing ventricular myocyte. This provided the basis for simulations of the electrical activity within a transmural ventricular strand. Our main goal was to elucidate the roles of electrophysiological and structural remodeling in setting the stage for malignant life-threatening arrhythmias.

Results

Simulation results illustrate how the presence of M cells and heterogeneous electrophysiological remodeling in the human failing ventricle modulate the dispersion of action potential duration and repolarization time. Specifically, selective heterogeneous remodeling of expression levels for the Na⁺/Ca²⁺ exchanger and SERCA pump decrease these heterogeneities. In contrast, fibroblast proliferation and cellular uncoupling both strongly increase repolarization heterogeneities. Conduction velocity and the safety factor for conduction are also reduced by the progressive structural remodeling during heart failure.

Conclusion

An extensive literature now establishes that in human ventricle, as heart failure progresses, gradients for repolarization are changed significantly by protein specific electrophysiological remodeling (either homogeneous or heterogeneous). Our simulations illustrate and provide new insights into this. Furthermore, enhanced fibrosis in failing hearts, as well as reduced intercellular coupling, combine to increase electrophysiological gradients and reduce electrical propagation. In combination these changes set the stage for arrhythmias.

Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 2D simulation study

BACKGROUND

Heart failure is operationally defined as the inability of the heart to maintain blood flow to meet the needs of the body and it is the final common pathway of various cardiac pathologies. Electrophysiological remodeling, intercellular uncoupling and a pro-fibrotic response have been identified as major arrhythmogenic factors in heart failure.

OBJECTIVE

In this study we investigate vulnerability to reentry under heart failure conditions by incorporating established electrophysiological and anatomical remodeling using computer simulations.

METHODS

The electrical activity of human transmural ventricular tissue (5 cm × 5 cm) was simulated using the human ventricular action potential model Grandi et al. under control and heart failure conditions. The MacCannell et al. model was used to model fibroblast electrical activity, and their electrotonic interactions with myocytes. Selected degrees of diffuse fibrosis and variations in intercellular coupling were considered and the vulnerable window (VW) for reentry was evaluated following cross-field stimulation.

RESULTS

No reentry was observed in normal conditions or in the presence of HF ionic remodeling. However, defined amount of fibrosis and/or cellular uncoupling were sufficient to elicit reentrant activity. Under conditions where reentry was generated, HF electrophysiological remodeling did not alter the width of the VW. However, intermediate fibrosis and cellular uncoupling significantly widened the VW. In addition, biphasic behavior was observed, as very high fibrotic content or very low tissue conductivity hampered the development of reentry. Detailed phase analysis of reentry dynamics revealed an increase of phase singularities with progressive fibrotic components.

CONCLUSION

Structural remodeling is a key factor in the genesis of vulnerability to reentry. A range of intermediate levels of fibrosis and intercellular uncoupling can combine to favor reentrant activity.

Carbon monoxide effects on human ventricle action potential assessed by mathematical simulations

Carbon monoxide (CO) that is produced in a number of different mammalian tissues is now known to have significant effects on the cardiovascular system. These include: (i) vasodilation, (ii) changes in heart rate and strength of contractions, and (iii) modulation of autonomic nervous system input to both the primary pacemaker and the working myocardium. Excessive CO in the environment is toxic and can initiate or mediate life threatening cardiac rhythm disturbances. Recent reports link these ventricular arrhythmias to an increase in the slowly inactivating, or “late” component of the Na⁺ current in the mammalian heart. The main goal of this paper is to explore the basis of this pro-arrhythmic capability of CO by incorporating changes in CO-induced ion channel activity with intracellular signaling pathways in the mammalian heart. To do this, a quite well-documented mathematical model of the action potential and intracellular calcium transient in the human ventricular myocyte has been employed. In silico iterations based on this model provide a useful first step in illustrating the cellular electrophysiological consequences of CO that have been reported from mammalian heart experiments. Specifically, when the Grandi et al. model of the human ventricular action potential is utilized, and after the Na⁺ and Ca²⁺ currents in a single myocyte are modified based on the experimental literature, early after-depolarization (EAD) rhythm disturbances appear, and important elements of the underlying causes of these EADs are revealed/illustrated. Our modified mathematical model of the human ventricular action potential also provides a convenient digital platform for designing future experimental work and relating these changes in cellular cardiac electrophysiology to emerging clinical and epidemiological data on CO toxicity.

Dynamics of the late Na(+) current during cardiac action potential and its contribution to afterdepolarizations

The objective of this work is to examine the contribution of late Na(+) current (INa,L) to the cardiac action potential (AP) and arrhythmogenic activities. In spite of the rapidly growing interest toward this current, there is no publication available on experimental recording of the dynamic INa,L current as it flows during AP with Ca(2+) cycling. Also unknown is how the current profile changes when the Ca(2+)-calmodulin dependent protein kinase II (CaMKII) signaling is altered, and how the current contributes to the development of arrhythmias. In this study we use an innovative AP-clamp Sequential Dissection technique to directly record the INa,L current during the AP with Ca(2+) cycling in the guinea pig ventricular myocytes. First, we found that the magnitude of INa,L measured under AP-clamp is substantially larger than earlier studies indicated. CaMKII inhibition using KN-93 significantly reduced the current. Second, we recorded INa,L together with IKs, IKr, and IK1 in the same cell to understand how these currents counterbalance to shape the AP morphology. We found that the amplitude and the total charge carried by INa,L exceed that of IKs. Third, facilitation of INa,L by Anemone toxin II prolonged APD and induced Ca(2+) oscillations that led to early and delayed afterdepolarizations and triggered APs; these arrhythmogenic activities were eliminated by buffering Ca(2+) with BAPTA. In conclusion, INa,L contributes a significantly large inward current that prolongs APD and unbalances the Ca(2+) homeostasis to cause arrhythmogenic APs.

Post-translational modifications of the cardiac Na channel: contribution of CaMKII-dependent phosphorylation to acquired arrhythmias

The voltage-gated Na channel isoform 1.5 (NaV1.5) is the pore forming α-subunit of the voltage-gated cardiac Na channel, which is responsible for the initiation and propagation of cardiac action potentials. Mutations in the SCN5A gene encoding NaV1.5 have been linked to changes in the Na current leading to a variety of arrhythmogenic phenotypes, and alterations in the NaV1.5 expression level, Na current density, and/or gating have been observed in acquired cardiac disorders, including heart failure. The precise mechanisms underlying these abnormalities have not been fully elucidated. However, several recent studies have made it clear that NaV1.5 forms a macromolecular complex with a number of proteins that modulate its expression levels, localization, and gating and is the target of extensive post-translational modifications, which may also influence all these properties. We review here the molecular aspects of cardiac Na channel regulation and their functional consequences. In particular, we focus on the molecular and functional aspects of Na channel phosphorylation by the Ca/calmodulin-dependent protein kinase II, which is hyperactive in heart failure and has been causally linked to cardiac arrhythmia. Understanding the mechanisms of altered NaV1.5 expression and function is crucial for gaining insight into arrhythmogenesis and developing novel therapeutic strategies.

The arrhythmogenic consequences of increasing late INa in the cardiomyocyte

This review presents the roles of cardiac sodium channel NaV1.5 late current (late INa) in generation of arrhythmic activity. The assumption of the authors is that proper Na(+) channel function is necessary to the maintenance of the transmembrane electrochemical gradient of Na(+) and regulation of cardiac electrical activity. Myocyte Na(+) channels' openings during the brief action potential upstroke contribute to peak INa and initiate excitation-contraction coupling. Openings of Na(+) channels outside the upstroke contribute to late INa, a depolarizing current that persists throughout the action potential plateau. The small, physiological late INa does not appear to be critical for normal electrical or contractile function in the heart. Late INa does, however, reduce the net repolarizing current, prolongs action potential duration, and increases cellular Na(+) loading. An increase of late INa, due to acquired conditions (e.g. heart failure) or inherited Na(+) channelopathies, facilitates the formation of early and delayed afterpolarizations and triggered arrhythmias, spontaneous diastolic depolarization, and cellular Ca(2+) loading. These in turn increase the spatial and temporal dispersion of repolarization time and may lead to reentrant arrhythmias.

Redox-dependent regulation of the Na⁺-K⁺ pump: new twists to an old target for treatment of heart failure

By the time it was appreciated that the positive inotropic effect of cardiac glycosides is due to inhibition of the membrane Na(+)-K(+) pump, glycosides had been used for treatment of heart failure on an empiric basis for ~200 years. The subsequent documentation of their lack of clinical efficacy and possible harmful effect largely coincided with the discovery that a raised Na(+) concentration in cardiac myocytes plays an important role in the electromechanical phenotype of heart failure syndromes. Consistent with this, efficacious pharmacological treatments for heart failure have been found to stimulate the Na(+)-K(+) pump, effectively the only export route for intracellular Na(+) in the heart failure. A paradigm has emerged that implicates pump inhibition in the raised Na(+) levels in heart failure. It invokes protein kinase-dependent activation of nicotinamide adenine dinucleotide phosphate oxidase (NADPH oxidase) and glutathionylation, a reversible oxidative modification, of the Na(+)-K(+) pump molecular complex that inhibits its activity. Since treatments of proven efficacy reverse the oxidative Na(+)-K(+) pump inhibition, the pump retains its status as a key pharmacological target in heart failure. Its role as a target is well integrated with the paradigms of neurohormonal abnormalities, raised myocardial oxidative stress and energy deficiency implicated in the pathophysiology of the failing heart. We propose that targeting oxidative inhibition of the pump is useful for the exploration of future treatment strategies. This article is part of a Special Issue entitled "Na(+)Regulation in Cardiac Myocytes".

Na/K pump regulation of cardiac repolarization: insights from a systems biology approach

The sodium-potassium pump is widely recognized as the principal mechanism for active ion transport across the cellular membrane of cardiac tissue, being responsible for the creation and maintenance of the transarcolemmal sodium and potassium gradients, crucial for cardiac cell electrophysiology. Importantly, sodium-potassium pump activity is impaired in a number of major diseased conditions, including ischemia and heart failure. However, its subtle ways of action on cardiac electrophysiology, both directly through its electrogenic nature and indirectly via the regulation of cell homeostasis, make it hard to predict the electrophysiological consequences of reduced sodium-potassium pump activity in cardiac repolarization. In this review, we discuss how recent studies adopting the systems biology approach, through the integration of experimental and modeling methodologies, have identified the sodium-potassium pump as one of the most important ionic mechanisms in regulating key properties of cardiac repolarization and its rate dependence, from subcellular to whole organ levels. These include the role of the pump in the biphasic modulation of cellular repolarization and refractoriness, the rate control of intracellular sodium and calcium dynamics and therefore of the adaptation of repolarization to changes in heart rate, as well as its importance in regulating pro-arrhythmic substrates through modulation of dispersion of repolarization and restitution. Theoretical findings are consistent across a variety of cell types and species including human, and widely in agreement with experimental findings. The novel insights and hypotheses on the role of the pump in cardiac electrophysiology obtained through this integrative approach could eventually lead to novel therapeutic and diagnostic strategies.

mRNA expression levels in failing human hearts predict cellular electrophysiological remodeling: a population-based simulation study

Differences in mRNA expression levels have been observed in failing versus non-failing human hearts for several membrane channel proteins and accessory subunits. These differences may play a causal role in electrophysiological changes observed in human heart failure and atrial fibrillation, such as action potential (AP) prolongation, increased AP triangulation, decreased intracellular calcium transient (CaT) magnitude and decreased CaT triangulation. Our goal is to investigate whether the information contained in mRNA measurements can be used to predict cardiac electrophysiological remodeling in heart failure using computational modeling. Using mRNA data recently obtained from failing and non-failing human hearts, we construct failing and non-failing cell populations incorporating natural variability and up/down regulation of channel conductivities. Six biomarkers are calculated for each cell in each population, at cycle lengths between 1500 ms and 300 ms. Regression analysis is performed to determine which ion channels drive biomarker variability in failing versus non-failing cardiomyocytes. Our models suggest that reported mRNA expression changes are consistent with AP prolongation, increased AP triangulation, increased CaT duration, decreased CaT triangulation and amplitude, and increased delay between AP and CaT upstrokes in the failing population. Regression analysis reveals that changes in AP biomarkers are driven primarily by reduction in I, and changes in CaT biomarkers are driven predominantly by reduction in I and SERCA. In particular, the role of I is pacing rate dependent. Additionally, alternans developed at fast pacing rates for both failing and non-failing cardiomyocytes, but the underlying mechanisms are different in control and heart failure.