Simulation and Mechanistic Investigation of the Arrhythmogenic Role of the Late Sodium Current in Human Heart Failure

Beatriz Trenor; Karen Cardona; Juan F. Gomez; Sridharan Rajamani; Jose M. Ferrero Jr; Luiz Belardinelli; Javier Saiz

doi:10.1371/journal.pone.0032659

Abstract

Heart failure constitutes a major public health problem worldwide. The electrophysiological remodeling of failing hearts sets the stage for malignant arrhythmias, in which the role of the late Na⁺ current (I_NaL) is relevant and is currently under investigation. In this study we examined the role of I_NaL in the electrophysiological phenotype of ventricular myocytes, and its proarrhythmic effects in the failing heart. A model for cellular heart failure was proposed using a modified version of Grandi et al. model for human ventricular action potential that incorporates the formulation of I_NaL. A sensitivity analysis of the model was performed and simulations of the pathological electrical activity of the cell were conducted. The proposed model for the human I_NaL and the electrophysiological remodeling of myocytes from failing hearts accurately reproduce experimental observations. The sensitivity analysis of the modulation of electrophysiological parameters of myocytes from failing hearts due to ion channels remodeling, revealed a role for I_NaL in the prolongation of action potential duration (APD), triangulation of the shape of the AP, and changes in Ca²⁺ transient. A mechanistic investigation of intracellular Na⁺ accumulation and APD shortening with increasing frequency of stimulation of failing myocytes revealed a role for the Na⁺/K⁺ pump, the Na⁺/Ca²⁺ exchanger and I_NaL. The results of the simulations also showed that in failing myocytes, the enhancement of I_NaL increased the reverse rate-dependent APD prolongation and the probability of initiating early afterdepolarizations. The electrophysiological remodeling of failing hearts and especially the enhancement of the I_NaL prolong APD and alter Ca²⁺ transient facilitating the development of early afterdepolarizations. An enhanced I_NaL appears to be an important contributor to the electrophysiological phenotype and to the dysregulation of [Ca²⁺]_i homeostasis of failing myocytes.

Citation: Trenor B, Cardona K, Gomez JF, Rajamani S, Ferrero JM Jr, Belardinelli L, et al. (2012) Simulation and Mechanistic Investigation of the Arrhythmogenic Role of the Late Sodium Current in Human Heart Failure. PLoS ONE 7(3): e32659. https://doi.org/10.1371/journal.pone.0032659

Editor: Marcello Rota, Brigham & Women's Hospital - Harvard Medical School, United States of America

Received: August 30, 2011; Accepted: February 2, 2012; Published: March 12, 2012

Copyright: © 2012 Trenor et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work was partially supported by the European Commission preDiCT grant (DG-INFSO-224381), by the Plan Nacional de Investigación Científica, Desarrollo e Innovación Tecnológica del Ministerio de Ciencia e Innovación of Spain (TEC2008-02090), by the Programa de Apoyo a la Investigación y Desarrollo (PAID-06-09-2843, PAID-06-11-2002) de la Universidad Politécnica de Valencia, by the Dirección General de Política Científica de la Generalitat Valenciana (GV/2010/078). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. This work was also supported by Gilead Sciences Palo Alto. Dr. Rajamani and Dr. Belardinelli are employees of Gilead Sciences and have collaborated in the study design and manuscript preparation. No additional external funding received for this study.

Competing interests: The authors have read the journal's policy and have the following conflicts: Dr. Belardinelli and Dr. Rajamani are employees of Gilead Sciences, the US company that markets Ranexa, and own stock and stock options. They are reimbursed for travel expenses to lecture on this subject. They have numerous patents related to late INA inhibitors, a list of which can be provided on request. This does not alter the authors' adherence to all the PLoS ONE policies on sharing data and materials.

Introduction

Over 5 million persons in the United States suffer from heart failure (HF) and more than 250,000 die annually [1]. Patients with congestive HF are prone to develop complex ventricular tachyarrhythmias and some die suddenly [2]. Experimental studies conducted using animal models of HF have shown that ventricular arrhythmias are mainly due to non-reentrant mechanisms, most likely triggered activity based on afterdepolarizations [2].

Much attention has been paid to the understanding of the arrhythmogenic mechanisms induced by the structural, electrical, and metabolic remodeling of the failing heart. The electrophysiological remodeling of the failing heart has been well described (see [3], [4] for review). Action potential (AP) prolongation, altered Ca²⁺ handling, as well as intracellular Na⁺ ([Na⁺]_i) accumulation have been established as the hallmark characteristics of myocytes and tissues isolated from failing human and canine hearts [5]–[8]. These alterations are closely related to arrhythmogenic mechanisms, such as early (EADs) and delayed (DADs) afterdepolarizations, observed in HF [9]. Functional remodeling of ion channels and pumps is the underlying cause for AP duration (APD) prolongation and altered intracellular Ca²⁺ ([Ca²⁺]_i) homeostasis. Downregulation of outward K⁺ currents is the most consistent change in animal models and human HF [3], [4], [6], [7], [10]. Major changes in intracellular and sarcoplasmic reticulum (SR) Ca²⁺ homeostasis are also associated with HF in several animal species, included human [9], [11]–[13]. In myocytes from failing hearts [Na⁺]_i concentration and Ca²⁺ handling are closely linked; [Na⁺]_i is increased in failing ventricular myocytes from human and other animal species [3], [4], [9], [14] and a prominent increase of the human late Na⁺ current (I_NaL) has also been documented [15], [16], and has been proposed as a therapeutic target [17]–[19]. Experimental studies have shown that the I_NaL is involved in the generation of EADs in myocytes [17], [18] and life-threatening arrhythmias, such as torsade de pointes (TdP) [20], especially under conditions of reduced repolarization reserve in several animal species [21]–[23]. Other changes in biomarkers for arrhythmic risk such as the increase in the reverse rate-dependent APD prolongation [20], [23], [24] have been attributed to an increase in I_NaL concomitant with inhibition of outward K⁺ currents.

Thus, the goal of our study was to analyze, using computer simulations, the role of the I_NaL in the setting of human HF. A mathematical model of human HF is proposed at cellular level using a modified version of the Grandi et al. model (from herein referred to as the GPB model) [25] of endocardial AP, in which a new formulation of the I_NaL is included. The sensitivity analysis performed for the HF model, as well as the simulations of the rate-dependence of APD and EAD generation revealed that the mechanisms underlying the arrhythmogenic processes taking place in HF conditions, from a theoretical point of view, include an enhanced I_NaL.

Results

Effects of I_NaL on human AP

Figure 1A shows the generation of I_NaL trace using a voltage clamp protocol similar to that of Maltsev et al. [26] in ventricular mid-myocardial myocytes at room temperature. The simulated I_NaL had similar time evolution and amplitude to the experimental findings. Panels B and C show the effects of this current on AP and APD at 90% of repolarization (APD₉₀), respectively. The inclusion of I_NaL in the GPB model slightly prolongs steady-state APD at 1 Hz (see Figure 1B). APD remains within the physiological range for human endocardial myocytes [8], [27], [28]. By increasing the magnitude of I_NaL by a factor of two-, five- or ten-fold, we obtained APD₉₀ prolongations of 11%, 44% and 78%, respectively. Experimental recordings of monophasic APs (MAPs) under the effects of veratridine, an enhancer of I_NaL, show similar effects on APD prolongation in rabbit myocytes [29]. We determined the sensitivity of APD₉₀ to I_NaL amplitude at different stimulation rates. As shown in Figure 1C, I_NaL/I_NaT ratios, where I_NaT is the transient peak I_Na, were varied from 0.0298% to 1.26% and APD₉₀ significantly increased according to the magnitude of I_NaL. Furthermore, the sensitivity of APD₉₀ increased with the slowing of the stimulation rate (compare 1 Hz and 0.5 Hz). Similar simulations were carried out in a rabbit model of AP [30]. As shown in Figure 1D, the difference in sensitivity of APD₉₀ to changes in frequency of stimulation (0.5 and 1 Hz) was also observed. Note that the morphology of the curves in Figure 1D is different to our results, as different models for I_NaL and different AP models were used. The ratio chosen in our model of 0.12% yields an APD₉₀ that is within the physiological range for 1 Hz indicated by the discontinuous lines.

Download:

Figure 1. Time course of late Na⁺ current and its effects on AP.

A: Simulation of late Na⁺ current (I_NaL) using a voltage clamp protocol similar to that of the experimental measurements obtained by Maltsev et al. [26] (shown in the right inset) in ventricular myocytes at room temperature. The left inset shows an expanded view of the current between 480 ms and 2500 ms. B: Simulated action potentials (APs) at 1-Hz pacing rate using the GPB model, the GPB model modified with control I_NaL, and I_NaL enhanced 2-fold, 5-fold and 10-fold. C: APD₉₀ sensitivity to the I_NaL amplitude. APs were simulated at 0.5-Hz (square symbols) and 1-Hz (circle symbols) pacing rate by varying I_NaL/I_NaT ratio from 0.0298% to 1.26%. The range of experimental APD at 90% repolarization (APD₉₀) for human is represented by the two discontinuous lines. D: APD₉₀ sensitivity to the I_NaL amplitude (open symbols) taken from Grandi et al. [30] who used a rabbit model.

https://doi.org/10.1371/journal.pone.0032659.g001

Heart Failure model and sensitivity analysis

APs and [Ca²⁺]_i transients were simulated under conditions of HF. Figure 2 shows APs (panel A), [Ca²⁺]_i transients (panel B) and Na⁺/Ca²⁺ exchanger (NCX) function (I_NCX) (panel C) of failing (light line) and nonfailing (dark line) ventricular myocytes after achieving steady-state conditions for a stimulation rate of 1 Hz. Our results showed an APD₉₀ prolongation of 24% in failing myocytes versus normal ones, as well as a 18% prolongation in APD₅₀, so that triangulation (APD₉₀–APD₅₀) increased by 43% under HF conditions. The experimental observations shown in the inset of Figure 2 (panel A) taken from reference [5], show the variability of the APD in failing myocytes from human hearts, which falls within the simulated values of APD observed.

Download:

Figure 2. Electrophysiological changes in heart failure.

Simulated APs (panel A), [Ca²⁺]_i transients (panel B), and Na⁺/Ca²⁺ exchanger (NCX) activity (panel C) at 1-Hz pacing rate in control (dark line) and in heart failure (HF) conditions (light line). The insets show experimental recordings of Priebe and Beuckelmann et al. [5] (panel A) and Weber et al. [31] (panel B).

https://doi.org/10.1371/journal.pone.0032659.g002

Figure 2B depicts an altered [Ca²⁺]_i transient under HF conditions, as has been documented experimentally and is shown in the inset [31]. Indeed, diastolic [Ca²⁺]_i is slightly increased, whereas peak systolic [Ca²⁺]_i is reduced to 41% of the one observed during normal conditions. An additional reported feature of the [Ca²⁺]_i transient of a failing myocyte is its slow decay. In our simulations Ca²⁺ transient decay, quantified as the time needed from the peak value to reach 10% of the transient amplitude (ô_Ca decay), yielded 630 ms and 380 ms in the failing and nonfailing myocytes, respectively. Finally, Figure 2C illustrates the changes in I_NCX during HF, mainly a shift in the time of the reversal potential for the NCX (t_NCXRP) of 20 ms. Similar shifts in t_NCXRP have been reported in experimental studies [31].

Changes in [Na⁺]_i and [Ca²⁺]_i levels under HF conditions at various stimulation rates were also investigated using a staircase protocol described in the methods section and described previously [32]–[34]. Figure 3A shows that diastolic [Ca²⁺]_i level is higher in HF than in normal conditions and systolic level is always higher in normal conditions than in HF, as reported experimentally [5], [35], [36]. The impact of the variability of ion channel remodeling on these results are presented in Figures S4, S5, S6, S7.

Download:

Figure 3. [Ca²⁺]_i and [Na⁺]_i changes with increasing frequency in HF.

Influence of the stimulation rate on [Ca⁺²]_i (panel A) and [Na⁺]_i in nmol/L (panel B) determined using the staircase protocol. The simulations were performed using the modified GPB model with I_NaL incorporated in normal conditions (blue) and the modified GPB model incorporating HF conditions (red). In panel B, conditions of HF without I_NaL, I_NaK or I_NCX remodeling are also shown.

https://doi.org/10.1371/journal.pone.0032659.g003

Another feature of myocytes from failing hearts is the high concentration of [Na⁺]_i regardless of the stimulation rate, as observed in Figure 3B and experimentally reported by Pieske et al. [33] in failing human myocardial cells. However, as acknowledged by Pieske et al. [33], it should be noted that the absolute [Na⁺]_i values in their experimental work were overestimated in normal and failing hearts. The impact of the variability of ion channel remodeling on these results are presented in Figures S8 and S9.

To analyze the mechanisms responsible for [Na⁺]_i accumulation with frequency in failing myocytes, the above pacing protocol was applied for a) normal conditions (blue), b) HF conditions (red), c) HF with no I_NaL upregulation (purple), d) HF with no Na⁺/K⁺ pump current (I_NaK) downregulation (green) and e) HF with no I_NCX upregulation (black). From the results obtained in Figure 3B, it can be deduced that [Na⁺]_i accumulation is mainly driven in HF by alterations in I_NaK, and much less by the I_NCX and I_NaL.

Our basic HF model, as described in Figure 2 and 3, reproduced the main changes in the electrophysiological (EP) characteristics of myocytes from failing hearts. However, the experimentally reported ranges for EP changes, as well as the magnitude of ion channels, transporters and respective current remodeling, vary substantially for different experimental settings and HF stages. Hence, we performed a sensitivity analysis to assess the impact of the main ionic parameters remodeled in HF on the described EP characteristics. Figure 4A highlights the role of I_NaL and I_NaK in APD₉₀ under conditions of HF. When I_NaL undergoes a 2-fold increase (i.e. doubling) with respect to the basic HF increase, APD is prolonged 22% and no change at all in this current leads to a shortening of the APD₉₀ by 10% with respect to basic HF conditions. This result is in agreement with experimental recordings in failing human myocytes, where I_NaL has a crucial role in APD [19]. Similarly, the downregulation of Na⁺/K⁺ pump in HF has a relevant effect on APD₉₀ shortening, and further reduction of the current leads to a decrease of 22% with respect to HF APD₉₀ value, which is in agreement with experimental observations [37]. Figure 4B reveals the importance of I_NaK, I_NCX, the background Ca²⁺ current (I_Cab), the leak Ca²⁺ current (I_leak), and I_NaL in determining the value of peak systolic [Ca²⁺]_i during HF. Regarding the regulation of [Na⁺]_i value in failing myocytes, Figure 4C shows the important role of I_NaK, I_Cab, the inward rectifier K⁺ current (I_K1), I_NCX and I_NaL. To summarize the sensitivity of some of the EP characteristics (1^st column) during HF to the altered ionic parameters (1^st row), Figure 5 shows the relative sensitivity normalized to the maximum sensitivity for that particular characteristic, as described in the methods section. The positive and negative signs indicate whether the change of the ionic current and the HF EP characteristic follow the same or inverse tendency, respectively. Percentages in each box indicate the maximum absolute sensitivity of the EP parameter correspondent to that row for all ionic properties. From this sensitivity analysis it could be deduced that APD was particularly sensitive to I_NaL and to I_NaK (green and burgundy colors in rows 1 and 2). Furthermore, I_K1, I_NCX and I_NaK have an important effect on AP triangulation. The main features of Ca²⁺ transient in HF (3 medium rows) were mainly influenced in this order by I_NaK, SERCA function, I_NCX, I_leak, I_Cab, and I_NaL. The SR Ca²⁺ concentration ([Ca²⁺]_SR) is also influenced by other currents but in this case I_leak becomes more important than I_NCX. [Na⁺]_i is mainly regulated by I_NaK, I_Cab, I_K1, I_NCX and I_NaL. Finally, t_NCXRP is mainly modulated by the SR Ca²⁺-ATPase activity (I_SERCA). The absolute effects of small changes (±15%) in the ionic remodeling of the basic HF model on the main results of our simulations are shown in Figures S1, S2, S3, S4, S5, S6, S7, S8, S9, S10, and S11. These absolute effects are in agreement with the relative role of the ionic parameters described above. We have especially focused on the sensitivity of APD₉₀ in failing myocytes to I_NaL/I_NaT ratio, as shown in Figure S1. In this figure, the values of I_NaL/I_NaT ratios for 1 Hz in the modified the GPB model and in our basic HF model are indicated by arrows and are within experimental ranges (indicated by blue dotted lines in normal myocytes and by red discontinuous lines in failing myocytes, respectively). Changes in the selected ratios affect the APD₉₀ values of failing and normal myocytes, which stay within experimental ranges provided that the changes of these ratios are not very large.

Download:

Figure 4. Sensitivity of electrophysiological parameters to changes in ionic current properties.

Changes in APD₉₀ (panel A), peak systolic [Ca²⁺]_i (panel B), and [Na⁺]_i (panel C) with changes in I_NaL, I_NaK, I_leak, I_NCX, I_K1 and I_Cab, as labeled next to each curve. Axis x represents the simulation conditions; for “HF Conditions” the remodeling of the basic HF model is considered, for “No Change” the labeled current is unchanged as it is in the GPB model, for “Double Change” the labeled current undergoes a double change with respect to the change exerted in “HF conditions.”

https://doi.org/10.1371/journal.pone.0032659.g004

Download:

Figure 5. Relative sensitivities of the electrophysiological parameters to changes in ionic current properties.

Dark blue color indicates lack of dependency between the ionic property (1^st row) and the EP parameter (1^st column), and dark red color indicates strong direct (+) dependency or inverse dependency (−). Percents in each box indicate the maximum absolute sensitivity of the EP parameter correspondent to that row for all ionic properties.

https://doi.org/10.1371/journal.pone.0032659.g005

APD rate-dependence and the role of I_NaL in HF

The role of I_NaL in steady-state APD₉₀ rate-dependence during HF was analyzed. An enhanced reverse frequency-dependence of APD₉₀ is considered to be a major proarrhythmic risk (see [38] for review). We performed simulations and measured APD₉₀ under normal and HF conditions for different magnitudes of I_NaL activation (see Figure 6A). Also 40% and 60% inhibition of the rapid component of the delayed rectifier K⁺ current (I_Kr) in failing myocytes were simulated with different degrees of I_NaL enhancement (see Figure 6B). Indeed, the concomitant alteration of I_NaL and I_Kr has been reported to have important effects on arrhythmogenicity [20]–[24]. An important result obtained in our simulations, as shown in Figure 6A, is the reverse rate-dependence effect on APD₉₀ exerted by I_NaL enhancement for normal conditions (compare circles with stars). The rate-dependence yielded 110 and 190 ms maximal APD₉₀ prolongation in normal conditions versus I_NaL increased by 4-fold, respectively. Similar observations were reported by others [20], [24] in rabbit ventricular myocytes under the effects of I_NaL activators, as well as in isolated perfused rabbit heart [22]. When considering HF conditions (squares), the rate-dependence is also increased with respect to normal conditions. The ΔAPD value is very similar to the value obtained with I_NaL enhanced (stars). This implies that I_NaL might be the main driver of the increased reverse rate-dependent prolongation of APD. Additionally, greater magnitudes of I_NaL (black triangles) are associated with greater increase in the reverse rate-dependency of APD of failing myocytes. Experimental studies [7] have also documented the increase of APD rate-dependence in human failing epicardial myocytes, as depicted in the inset. Simulations of APD rate-dependence in HF caused by small changes (±15%) in the ionic remodeling of the basic HF model (see Figures S2, S3, S4, S5, S6, S7, S8, S9, S10, and S11) to take into account the experimental variability on electrical remodeling during HF were also performed. The results do not change significantly with respect to the basic HF model, except for I_NaK and I_NaL, which have an important effect on APD₉₀, as predicted by the sensitivity analysis shown in Figure 5.

Download:

Figure 6. Role of I_NaL in APD rate-dependence in HF and reduced repolarization reserve.

Simulated APD dependence on stimulation frequency for normal (circles in panel A) and HF (squares in panel A) conditions and in the presence of enhanced I_NaL in normal (stars in panel A) and HF (triangles in panel A) conditions. The inset shows experimental results of Li et al. [7]. Panel B, depicts the APD rate-dependence for HF combined with different degrees of I_Kr inhibition and I_NaL enhancement.

https://doi.org/10.1371/journal.pone.0032659.g006

Figure 6B summarizes the concomitant contribution of I_NaL enhancement and I_Kr inhibition in the reverse rate-dependency of APD, shown in various experimental studies [20], [24], in failing myocytes.

To analyze the ionic mechanisms responsible for steady-state APD rate-dependence in failing myocytes, we determined the behavior of several ionic currents (I_NaL, I_NCX and I_NaK) and [Na⁺]_i at different frequencies of stimulation, as depicted in Figure 7. Our results show that the shortening of the APD with increased frequency in HF is explained on the one hand by the decrease in I_NaL, indeed faster rates lead to an incomplete recovery from inactivation of this current [39]–[41], and on the other hand by [Na⁺]_i accumulation leading to an increase of I_NaK. The role of I_NCX on APD is not evident as both outward and inward modes are enhanced. Although several experimental studies suggest the direct involvement of this current on APD shortening with frequency, as will be discussed later, its implication in our simulations is not evidenced. Indeed, when we stimulated the cell at high frequency (1.6 Hz) but clamped [Na⁺]_i value to the lower [Na⁺]_i value corresponding to the low frequency (0.5 Hz), we observed a decrease in APD₉₀ with respect to APD₉₀ at low frequency, concomitant with a decrease and increase of the outward and inward modes of NCX activity, respectively (results not shown). These changes in I_NCX should lead to a longer APD₉₀, instead of shorter, suggesting that I_NCX is not determinant in APD rate-dependence in the present model.

Download:

Figure 7. Mechanisms for APD rate-dependence in HF.

Simulated APs and ionic currents at different stimulating frequencies (0.5, 1, 1.25 and 1.6 Hz) under HF conditions.

https://doi.org/10.1371/journal.pone.0032659.g007

Arrhythmogenic effects of I_NaL

Under basic HF conditions, APD prolongation favors the occurrence of EADs. We were therefore interested in determining the role of I_NaL in the repolarization abnormalities of human failing hearts. We stimulated the myocyte at 1 Hz and simulated conditions of low repolarization reserve prone to EADs generation. Under HF conditions using the basic HF model, I_Kr was reduced by 50% and the Ca²⁺ current (I_CaL) was increased by 30%. Figure 8A shows EADs (dark line) when I_NaL was doubled, whereas APs displayed no EADs when I_NaL was normal (light line). The impact of the variability in ion channel remodeling on EADs generation can be observed in Figures S10 and S11. The mechanisms by which EADs arise are dictated by a very delicate balance of ionic currents during AP plateau. Slight changes in this balance can suppress EADs, as occurs in the case when all ionic parameters are reduced in 15% with respect to the basic HF model (see Figure S11).

Download:

Figure 8. Mechanisms for early afterdepolarizations with enhanced I_NaL.

Simulated APs and ionic currents at a stimulating rate of 1 Hz for HF conditions, 50% inhibition of I_Kr and 30% increase of I_CaL. Panel A shows early afterdepolarizations (EADs; dark line) when I_NaL was doubled and APs with no EADs when I_NaL was normal (light line). The temporal evolutions of I_NaL (panel B), I_rel (panel D), I_CaL (panel C), NCX activity (panel E), and activation gate of I_CaL are also depicted when I_NaL was doubled (dark line) and when I_NaL was normal (light line).

https://doi.org/10.1371/journal.pone.0032659.g008

Figure 8 depicts APs and ionic currents to provide insight into the mechanisms by which enhancement of I_NaL indirectly triggers EADs [17]. The APD prolongation (panel A) caused by enhancement of I_NaL (panel B) leads to the reactivation of the activation gate d (panel F) of the I_CaL, which triggers the EAD, as was suggested by January et al. [42]. However, the important contribution of the altered Ca²⁺ transient should also be taken into consideration. Indeed, the APD prolongation caused by I_NaL enhancement alters Ca²⁺ handling. As observed in panel C, [Ca²⁺]_SR remains elevated and provokes the trigger of a small spontaneous release (panel D, arrow) just before the EAD arises (panel A, arrow). NCX (panel E) extrudes Ca²⁺ operating in the inward mode, which contributes to membrane potential increase. This would contribute to I_CaL reactivation which triggers the EAD. The involvement of these currents is also shown in Figures S12, S13, S14, and S15, as a partial block of SR Ca²⁺ release (I_rel) or I_NCX suppresses the EADs.

Discussion

Major Findings

In this article, one of the most recent and detailed AP model for human ventricular myocytes, the GPB model, was modified to simulate HF functional remodeling. The I_NaL was formulated and introduced in the model on the basis of human ventricular voltage-clamp data. An important novelty of the study is the sensitivity analysis of the modulation of the main EP characteristics due to ion channel and transporters remodeling, highlighting the role of I_NaL in the changes of APD, triangulation and Ca²⁺ transient under HF conditions. A mechanistic investigation of [Na⁺]_i accumulation and APD shortening with increasing stimulation frequency in failing myocytes, revealed the roles of I_NaK, I_NCX and I_NaL. Finally, the arrhythmogenic effects of increased I_NaL were investigated. Our results showed the important role of I_NaL in APD reverse rate-dependence and EAD generation in failing myocytes, representing the first realistic simulation study on this issue, which complements experimental studies and sheds light into the ionic mechanisms responsible for HF phenotype.

Simulation of heart failure

Recently, O'Hara et al. [43] published a new AP model for human ventricular myocytes; however, in the present study the GPB model was chosen because its behavior was adequate to analyze the changes exerted in HF in APD, Ca²⁺ transient, the APD rate-dependence and EADs generation. Nevertheless future studies comparing the behavior of both models in pathological situations should be conducted. The GPB model was developed using human voltage-clamp data, and validated against recent experimental results [25]. In the present study, we demonstrate for the first time that the GPB model, subject to ion channel and transporter remodeling, is also appropriate to reproduce the electrical behavior of failing myocytes.

A large body of experimental research has focused on the functional remodeling of failing hearts [3], [4], [6], [7], [10], [11], [14]. However, little has been done in this field from a theoretical perspective. One of the first human ventricular AP models, including HF remodeling, was formulated by Priebe and Beuckelmann in 1998 [5]. However, I_NaL was not considered, and the formulation of Ca²⁺ handling was not as realistic as it is in the GPB model. This represents an important limitation for HF, when Ca²⁺ transient undergoes critical changes. Later, in 2005, Zhang et al. [44] modified the human ten Tusscher model [45] to simulate HF. In the latter model, I_NaL was considered, but with experimental data derived using dog ventricular myocytes [46], and remodeling was derived using experimental data from multiple species. However, controversy about the changes introduced in their model for I_NaT, the slow component of the delayed rectifier K⁺ current (I_Ks) or I_CaL exists for human failing myocytes [2]–[4]. Other authors have simulated ion remodeling of myocytes from canine and rabbit failing hearts [12], [47] that reproduce the APD and Ca²⁺ transient changes associated with HF. Our model of basic HF not only reproduces experimental observations on APD prolongation and Ca²⁺ transients alterations [5], [7], [11], [31], but also the APD, [Ca²⁺]_i and [Na⁺]_i rate-dependence reported in failing myocytes [33], [35], [36]. The ionic model of ventricular HF presented in this study is not only based on recent human experimental data, but also takes into account the variability documented in EP studies using sensitivity analysis. Sensitivity analysis also elucidates the ionic mechanisms responsible for the main changes in the EP characteristics of HF. Previous sensitivity analysis have been performed using the human ten Tusscher et al. [48] AP model [32], a modified version of the GPB model [34] and O'Hara et al. model [43], as well as rabbit [48] AP model. In these studies, the authors analyzed the effects of the physiological variability of the main ionic properties on several electrophysiological characteristics and biomarkers for arrhythmic risk, under normal physiological conditions. The main similarities of our results in HF with their studies on normal myocytes, are the importance of I_NaK [32], [34] and I_NaL [43] in regulating APD₉₀, and I_K1 [32], [34] and I_NaK in regulating AP triangulation. The variability of I_CaL and I_Kr has also an important effect in their studies but was not considered in HF, as these currents are not remodeled. Also the impact of I_NaK and I_NCX was important on Ca²⁺ transients under normal conditions [32], [34].

Regarding the variability of functional remodeling in HF, it is worth noting that experimental measurements reported in literature have been conducted in different stages of HF (see [3], [4] for review). In the present study, ionic currents are within experimental ranges, but do not take into account differences in various stages of HF. Different combinations of changes were tested in the sensitivity analysis in order to determine the relative effects on the EP characteristics in HF (Figure 5) and the absolute effects of small changes (±15%) in ionic parameters on APD-rate-dependence (Figures S2 and S3), [Ca²⁺]_i transients (Figures S4, S5, S6, and S7), [Na⁺]_i values (Figures S8 and S9) and EADs generation (Figures S10 and S11). It should be noted that the results of the present sensitivity analysis are dependent on the values chosen for the ionic parameters in the basic HF model, as the percents of change also depend on these chosen values. The results suggest that the downregulation of the Na⁺/K⁺ pump has an important role in the electrophysiological and Ca²⁺ transient changes in HF, highlighting the need for new experimental data. Electrophysiological studies [49] conducted in failing human hearts revealed increased Ca²⁺ transients under the effects of strophanthidin, a Na⁺/K⁺ pump blocker, which is in agreement with our results. However, little is known about the effects of Na⁺/K⁺ pump activity on APD in HF. Experiments carried out on hypertrophic rat hearts by Levi et al. [37] showed that the application of ouabain, a Na⁺/K⁺ pump blocker, resulted in APD shortening, which was caused by [Na⁺]_i accumulation and the subsequent increase in the reverse mode of the NCX. Similar studies should be performed using human failing hearts to determine the role of the effects of Na⁺/K⁺ pump on APD. Also the [Na⁺]_i accumulation with increasing frequency during HF appears to be mainly due to the downregulation of this pump, and a small contribution of I_NaL and I_NCX upregulation (see Figure 3B). Experimental studies reporting [Na⁺]_i accumulation in human HF [33] or animal species [35], [50] suggest that I_NaK downregulation, I_NaL and I_NCX upregulation, and altered activity of the Na⁺/H⁺ exchanger (not included in the GPB model), might be involved. However, no experiments have been performed to in human failing myocytes to clarify the responsible mechanisms. Other simulation studies [30], [51] using the Shannon et al. rabbit AP model [52] have also reported the small and important influence of I_NaL and I_NaK, respectively, on [Na⁺]_i accumulation with frequency under pathological situations other than HF.

Our results also reveal the role of I_NaK in APD rate-dependence as stated by Carmeliet [53] and by Eisner et al. [54] in their review on AP rate adaptation as well as shown in the theoretical and experimental findings by Pueyo et al. [55]. Additional simulation studies using different human AP models have also demonstrated this fact in normal conditions [32], [43]. Although O'Hara et al. [43] and Carro et al. [34] human AP models exhibit an improved behavior in the fast phase of APD rate adaptation than the GBP model, our study focused on the steady-state APD rate-dependence, where the GPB model is valid except for low frequencies [43]. Also Faber and Rudy [56] showed that APD shortening at high rates in a guinea pig AP model was due to the increase of the outward NCX activity. This observation is also consistent with experiments conducted on myocytes from human failing hearts [49], where increases in [Na⁺]_i led to APD shortening provoked by an increase of the reverse mode of I_NCX. However, the GPB model exhibits a much shorter outward mode of NCX than the Faber and Rudy model. Thus its contribution to AP repolarization is not evident in the present study. Furthermore, in failing human myocytes, Ca²⁺ influx via NCX is prominent during most of the plateau phase [31]. Thus changes in I_NCX formulation of the GPB model should be addressed to better understand the role of this current in HF APD rate-dependence. Also new experiments in failing human myocytes demonstrating the role of I_NaK and I_NCX in APD rate-dependence would be of great interest. This finding also suggests that the functional increase of these outward currents caused by [Na⁺]_i accumulation may limit AP prolongation in HF, counterbalancing the increase of I_NaL.

Our results also highlight the pivotal role of I_NaL in the changes of EP characteristics in HF, which will be discussed in the following section. Finally, another mechanistic effect of relative importance in HF, revealed by the simulations of the present study is the increase of AP triangulation associated with I_K1 downregulation. Similar observations have been reported in ventricular myocytes from normal dogs when I_K1 was reduced with barium [57] and in the sensitivity analysis performed by Romero et al. [32] and Carro et al. [34] in normal physiological conditions.

Important role of I_NaL during heart failure and clinical implications

The pivotal role of I_NaL in the changes of EP characteristics and in arrhythmogenesis in HF has been uncovered in the present work. The formulation of the I_NaL model and its introduction in the GPB model was essential to draw our conclusions. Previous formulations have been proposed to model the behavior of I_NaL for different animal species using the Hodgkin-Huxley or Markov formalisms [19], [20], [30], [46], [58], [59]. Noble and colleagues [59] in 1998 formulated the I_NaL for guinea pig by adding to their AP model a Na⁺ current with reduced conductance and slow inactivation. Similarly, other Hodgkin-Huxley formulations have been proposed for rabbit [20] and dog [19], [46]. To avoid the complexity of Markov models for I_NaL [30], we adapted the Hodgkin-Huxley mathematical description of Hund et al. [46] to human data, as we did in a previous work using ten Tusscher et al. AP model [48]. Our model has been validated against experimental recordings of human APD in normal conditions [8], [27], [28], APD prolongation under the effect of drugs [29] and voltage clamp experiments [58]. Recently, a new AP model for human ventricular myocytes has been published by O'Hara et al. [43] including a formulation for I_NaL as in the present work. However, as their AP model includes CaMK regulation, their formulation of I_NaL also includes the effect of CaMK regulation.

In agreement with previous experimental observations, the upregulation of I_NaL during HF [15] or its enhancement in situations of reduced repolarization reserve [21] prolongs APD, causes AP triangulation and increase the reverse rate-dependent prolongation of APD, which are important harbingers for cardiac arrhythmias. The contribution of I_NaL in the reverse APD rate-dependence was also previously simulated by our group [60] using a different human AP model [48] and conditions of reduced repolarization reserve. In regard to HF conditions, there is background of decreased outward currents, and even a small increase of I_NaL becomes more efficient in prolonging APD and favoring the trigger of EADs. The high probability of formation of EAD, leading to TdP, during HF has been experimentally demonstrated [7] but few studies have addressed the role of I_NaL in such situations [18], [26], [61]. Our results are in keeping with experimental findings, in which the I_NaL blocker ranolazine, effectively suppresses EADs from ventricular myocytes from failing hearts [18], [23]. Furthermore, during conditions of reduced repolarization reserve, in which outward K⁺ currents are inhibited by drugs or by diseases, results of experimental studies [21]–[23], [62] have unmasked the role of endogenous or enhanced I_NaL in exerting proarrhythmic effects. This is in agreement with the results of our simulations, and provides evidence of this effect in HF. Our findings support the hypothesis that I_NaL represents an important target for triggered-arrhythmias treatment [13]. Similarly, inhibition of CaMKII, which is known to be responsible for I_CaL and ryanodine receptor phosphorylation as well as I_NaL regulation, appears to be an important therapeutic target for suppressing arrhythmias in HF [13], [17], [63].

Limitations of the study

Several limitations need to be considered, when drawing conclusions from the present study. The model for I_NaL was formulated on the basis of voltage-clamp human data for mid-myocardial cells [58], as no measurements were available for endocardial or epicardial cells. Indeed, I_NaL is difficult to record due to the very low magnitude (∼30–60 pA) of this current under normal conditions [58]. Maltsev and Undrovinas reported a double exponential decay of I_NaL [58]. However, we chose a single-exponential decay formulation as proposed in previous studies [46], [64], and also in the new AP model by O'Hara et al. [43], as this model was able to reproduce its main effects on the AP. To build the model of HF used in the present study, the ionic remodeling was mainly based on experimental data observed in human hearts. Data from a large number of experimental studies have been taken into account, thus resulting in a high variability not only in the ionic remodeling but also in the stage of HF. In addition to the difficulties associated with gaining access to human hearts, explanted diseased hearts are usually in the end stage of HF. Moreover, there are controversies regarding specific ion channel currents remodeling during HF, such as I_NaT, I_CaL I_Kr, and I_Ks, and changes in these currents were not included in our HF model. Their effects on the EP characteristics would be significant contributors to the phenotype. We did however simulate changes in I_NaL concomitant with inhibition of I_Kr and increase in I_CaL (Figures 6 and 8). Finally, our basic HF model has the inherited limitations described for the GPB model. In this sense, improvements related to the rapid phase of APD rate adaptation have been accomplished in later human AP models [34], [43]. The response of failing myocytes to abrupt changes of the stimulation rate was however out of the scope of the present study because of the lack of availability of experimental data in HF. With regard to steady-state APD rate-dependence, the GPB model is accurate except for very low frequencies [43]. Thus the results obtained for APD₉₀ at low frequencies should be taken with caution.

Conclusion

This study aimed to investigate in silico the role of I_NaL in the electrophysiological and Ca²⁺ homeostasis phenotype of myocytes from failing hearts. Our results showed that the enhancement of this current during HF can lead to important prolongation of APD and triangulation, increases in reverse rate-dependent prolongation of APD, significantly contributes to Ca²⁺ handling changes and has an indirect but pivotal role in the genesis of EADs.

Methods

Model of the human ventricular I_NaL

To simulate the electrical activity of human ventricular myocytes, the AP model formulated by Grandi et al. [25] for endocardial cells was used. This model is one of the latest and most detailed mathematical model for ionic currents and Ca²⁺ handling of the human ventricular AP. A particular strength of the GPB model is its ability to reproduce the rate-dependence of APD upon outward K⁺ currents block and their individual role in repolarization. Thus, this model provides a powerful tool to explore repolarization abnormalities under conditions of disease, such as HF. However, in order to realistically simulate HF, an important issue remains unresolved in this model, namely the role of I_NaL.

We included in the GPB model a new formulation for human I_NaL that is described in equations (1) to (4) and is based on the formulation we included previously [60] in the AP model formulated by ten Tusscher et al. [48].(1)(2)(3)(4)We adopted these equations following Hodgkin Huxley formalism, where the activation gate (m_L) and the Nernst potential for Na⁺ (E_NaL) were unchanged with respect to I_NaT formulation in the GPB model. The steady-state inactivation gate (h_L∞) was taken from Maltsev et al. [64], as did Hund et al. in their model for canine ventricular myocytes [46]. The maximum conductance (g_NaL) and the time constant of inactivation (τ_hL) were fitted to reproduce I_NaL data taken from human mid-myocardial myocytes [58]. In their experiments, Maltsev et al. [58] measured an I_NaL/I_NaT ratio of approximately 0.1% using a voltage clamp protocol at room temperature. In our model g_NaL was fitted accordingly using voltage clamp simulations, yielding 0.015 mS/µF. No correction factor for temperature was applied, as we assumed that I_NaL/I_NaT ratio does not change with temperature in myocytes, as reported in transfected HEK-293 cells [65]. However, because experimental data indicate that ion channel dynamics are altered by temperature, the time constant of inactivation (τ_hL) of I_NaL reported by Maltsev et al. [26] for human ventricular mid-myocardial myocytes was multiplied by a Q₁₀ factor of 2.2 [65], yielding 233 ms at 37°C.

Heart failure cellular model for human

The remodeling of myocyte electrophysiology in HF is well described [3] and can explain for the most, the hallmark characteristics of failing cardiac tissues and myocytes, such as AP prolongation and alteration of Ca²⁺ and Na⁺ handling [2], [6], [66]. On the basis of experimental observations [3], [4], [7], [9], [11], [13]–[15], [67] and previous simulation studies [5], [12], [44], [47], we hereby propose various changes in the formulation of several ionic currents of the GPB model to reproduce the reported experimental changes in AP and intracellular Ca²⁺ and Na⁺ handling in ventricular myocytes from failing human hearts [8], [31], [66]. Our model will be referred to as the basic HF model. Table 1 summarizes the changes we made in the different ionic properties with respect to the GPB model. The ionic remodeling is mainly based on experimental data observed in human hearts.

Download:

Table 1. Heart failure remodeling.

https://doi.org/10.1371/journal.pone.0032659.t001

A prominent increase in I_NaL and slowing of current decay has been described in ventricular mid-myocardial myocytes isolated from failing hearts of dogs and humans [15], [67]. Accordingly, the current density and the time constant of inactivation were increased two-fold compared to non-failing cells. Although Maltsev et al. [15] reported a smaller increase in both values in failing myocytes, their experiments were performed at room temperature and stated that the changes are expected to increase at physiological temperature. Furthermore, the values of the time constant of inactivation could be measured for human myocytes, but the current density was only given for dog myocytes, as the authors state the difficulty to measure it in human myocytes where the variability was very high. Also, in a previous study by Valdivia et al. [16], I_NaL was increased 5-fold in human failing myocytes with respect to normal myocytes. Downregulation of K⁺ currents is the most consistent ionic current change observed in myocytes isolated from failing hearts from animal models and humans [6], [7], [10]. Mainly, I_to is downregulated without a significant change in the voltage dependence or kinetics [10]. We reduced the transient outward K⁺ current (I_to) to 40% of its normal value. Reported changes in I_K1 functional expression are more variable than I_to [3], [10] and have a strong dependence on the etiology. The conductance of this ion channel was multiplied by 0.68 as in [12], a value within the experimental range. In regards to [Na⁺]_i handling, which is also altered in HF, we reduced I_NaK activity by 10%, as the preponderance of experimental data reveal that the expression and function of the Na⁺/K⁺-ATPase are reduced in HF [3], [4], [14]. Similar to that of Priebe and Beuckelmann simulation study [5], the background Na²⁺ current (I_Nab) was not included in our HF model. Finally, the changes in intracellular and SR Ca²⁺ homeostasis were achieved by an increase of 75% of I_NCX [5], [12], [68], and a decrease of 50% of I_SERCA [11]. To reproduce the experimentally observed changes in Ca²⁺ sensitivity of the ryanodine receptor (RyR) [13], [69], I_leak was increased 3-fold and the parameter EC_50SR for SR [Ca²⁺]-dependent activation of SR release (see Grandi [25] supplementary data) was reduced to 0.4 mM. The background Ca²⁺ current was changed as in the study of Priebe and Beuckelmann [5].

Sensitivity analysis of the HF model

A sensitivity analysis was performed to investigate how the reported variability in HF remodeling data might modulate the main EP characteristics in HF. These characteristics “c” include APD at 90% and 50% of repolarization (APD₉₀ and APD₅₀), triangulation, peak systolic and diastolic [Ca²⁺]_i transient, ô_Ca decay, peak systolic and diastolic [Ca²⁺]_SR, [Na⁺]_i, and t_NCXRP. These characteristics were calculated at steady-state HF conditions (HF basic model) and after varying one parameter “p” at a time. The parameters included each of the ionic current properties modified in the HF basic model, and were varied to its normal value, as in the GPB model, and to a value representing twice that observed in the HF basic remodeling. Although an important change has been implemented (double change) to be considered as a high variability, only the sensitivity, i.e. the relative effect was analyzed. The sensitivity analysis performed similar to that of Romero et al. [32], wherein the indexes percentage of change (D_c,p,x) and sensitivities (S_c,p) were calculated as follows:(5)(6)with C_p,x and C_HFbasic being the magnitude of the characteristic c when parameter p undergoes a double change (x = 2) with respect to the HF basic change, or no change at all (x = 1); Δ_a is the total interval of change of parameter p.

Stimulation protocols

Voltage clamp was used to simulate I_NaL behavior. A voltage pulse to −20 mV was applied from a holding potential of −120 mV as depicted in the inset of Figure 1A. To calculate I_NaL/I_NaT ratios (see Figure 1C) I_NaL was measured at 40 ms after the application of the depolarizing pulse and was divided by peak I_NaT.

Cellular simulations were conducted at a stimulation rate of 1 Hz. Measurements were taken after achieving steady-state conditions.

Changes in [Na⁺]_i and [Ca²⁺]_i levels at various stimulation rates were measured using a staircase protocol as previously described [32], [33]. Cardiomyocytes were stimulated at increasing fast frequencies (0.5, 1, 1.5, 2 and 2.5 Hz) for 10 minutes, and [Ca²⁺]_i and [Na⁺]_i levels were recorded for each of the frequencies.

To analyze the APD rate-dependence, simulations were performed at different frequencies of stimulation (0.5, 1, 1.25 and 1.6 Hz). Measurements of APD₉₀ were taken after achieving the steady-state; likewise for APs, [Na⁺]_i and several ionic currents. The rate-dependence was measured as the difference between the maximum APD₉₀ for a stimulation rate of 0.5 Hz and the minimum APD₉₀ corresponding to the highest stimulation frequency.

Supporting Information

Figure S1.

Sensitivity of APD₉₀ to the I_NaL amplitude in HF. Steady-state APs at 1 Hz (square symbols and solid lines) and 0.5 Hz (circle symbols and dashed lines) pacing rates with varying I_NaL/I_NaT for normal conditions using the GPB model (thick lines) and under basic HF conditions (thin lines) where I_NaL is doubled with respect to normal conditions. The range of experimental APD₉₀ for human for normal conditions is represented by the two dotted blue lines. The range of experimental APD₉₀ for human for HF conditions is represented by the two discontinuous red lines.

https://doi.org/10.1371/journal.pone.0032659.s001

(TIF)

Figure S2.

Sensitivity of APD rate-dependence to variations in individual ionic parameters in HF. The steady-state APD₉₀ for different stimulation frequencies is shown for normal conditions using the GPB model (thick line), for basic HF conditions (solid line), and for a 15% increase (long dashed line) and a 15% reduction (short dashed line) of one ionic parameter with respect to its value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s002

(TIF)

Figure S3.

Sensitivity of APD rate-dependence to variations in all ionic parameters in HF. The steady-state APD₉₀ for different stimulation frequencies is shown for normal conditions using the GPB model (thick line), for basic HF conditions (solid line), and for a 15% increase (long dashed line) and a 15% reduction (short dashed line) of all the ionic parameters simultaneously with respect to their value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s003

(TIF)

Figure S4.

Sensitivity of rate-dependent changes in systolic [Ca²⁺]_i to variations in individual ionic parameters in HF. Systolic [Ca²⁺]_i after 10 minutes of stimulation at increasing rates is shown for normal conditions using the GPB model (thick line), for basic HF conditions (solid line), and for a 15% increase (long dashed line) and a 15% reduction (short dashed line) of one ionic parameter with respect to its value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s004

(TIF)

Figure S5.

Sensitivity of rate-dependent changes in systolic [Ca²⁺]_i to variations in all ionic parameters in HF. Systolic [Ca²⁺]_i after 10 minutes of stimulation at increasing rates is shown for normal conditions using the GPB model (thick line), for basic HF conditions (solid line), and for a 15% increase (long dashed line) and a 15% reduction (short dashed line) of all the ionic parameters simultaneously with respect to their value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s005

(TIF)

Figure S6.

Sensitivity of rate-dependent changes in diastolic [Ca²⁺]_i to variations in individual ionic parameters in HF. Diastolic [Ca²⁺]_i after 10 minutes of stimulation at increasing rates is shown for normal conditions using the GPB model (thick line), for basic HF conditions (solid line), and for a 15% increase (long dashed line) and a 15% reduction (short dashed line) of one ionic parameter with respect to its value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s006

(TIF)

Figure S7.

Sensitivity of rate-dependent changes in diastolic [Ca²⁺]_i to variations in all ionic parameters in HF. Diastolic [Ca²⁺]_i after 10 minutes of stimulation at increasing rates is shown for normal conditions using the GPB model (thick line), for basic HF conditions (solid line), and for a 15% increase (long dashed line) and a 15% reduction (short dashed line) of all the ionic parameters simultaneously with respect to their value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s007

(TIF)

Figure S8.

Sensitivity of rate-dependent changes in [Na⁺]_i to variations in individual ionic parameters in HF. [Na⁺]_i after 10 minutes of stimulation at increasing rates is shown for normal conditions using the GPB model (thick line), for basic HF conditions (solid line), and for a 15% increase (long dashed line) and a 15% reduction (short dashed line) of one ionic parameter with respect to its value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s008

(TIF)

Figure S9.

Sensitivity of rate-dependent changes in [Na⁺]_i to variations in all ionic parameters in HF. [Na⁺]_i after 10 minutes of stimulation at increasing rates is shown for normal conditions using the GPB model (thick line), for basic HF conditions (solid line), and for a 15% increase (long dashed line) and a 15% reduction (short dashed line) of all the ionic parameters simultaneously with respect to their value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s009

(TIF)

Figure S10.

Sensitivity of EAD generation to variations in individual ionic parameters in HF. Steady-state APs at 1-Hz pacing rate with 50% inhibition of I_Kr, 30% increase of I_CaL. The simulated results using the basic HF model are shown with a thick line, the solid and dashed lines show the results obtained for a 15% increase and a 15% reduction, respectively, of one ionic parameter with respect to its value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s010

(TIF)

Figure S11.

Sensitivity of EAD generation to variations in all ionic parameters in HF. Steady-state APs at 1-Hz pacing rate with 50% inhibition of I_Kr, 30% increase of I_CaL. The simulated results using the basic HF model are shown with a thick line, the solid and dashed lines show the results obtained for a 15% increase and a 15% reduction, respectively, of all the ionic parameters simultaneously with respect to their value in the basic HF model.

https://doi.org/10.1371/journal.pone.0032659.s011

(TIF)

Figure S12.

Mechanisms for early afterdepolarizations with a 50% continuous block of I_rel. Simulated APs and ionic currents at 1-Hz pacing rate under HF conditions, 50% inhibition of I_Kr, 30% increase of I_CaL. Panel A shows EADs (dark line) with the basic HF value of I_rel and APs with no EADs when I_rel was 50% blocked (light line). The temporal evolutions of I_NaL (panel B), [Ca²⁺]_SR (panel C), I_rel (panel D), NCX activity (panel E), and activation gate of I_CaL (panel F) are also depicted with the basic HF value of I_rel (dark line) and when I_rel was 50% blocked (light line).

https://doi.org/10.1371/journal.pone.0032659.s012

(TIF)

Figure S13.

Mechanisms for early afterdepolarizations with a 50% transitory block of I_rel. Simulated APs and ionic currents at 1-Hz pacing rate under HF conditions, 50% inhibition of I_Kr, 30% increase of I_CaL. Panel A shows EADs (dark line) with the basic HF value of I_rel and APs with no EADs when I_rel was 50% blocked (light line) during the 5 stimulations indicated by a horizontal line. The temporal evolutions of I_NaL (panel B), [Ca²⁺]_SR (panel C), I_rel (panel D), NCX activity (panel E), and activation gate of I_CaL (panel F) are also depicted with the basic HF value of I_rel (dark line) and when I_rel was 50% blocked (light line) during the 5 stimulations indicated by a horizontal line.

https://doi.org/10.1371/journal.pone.0032659.s013

(TIF)

Figure S14.

Mechanisms for early afterdepolarizations with a 50% continuous block of I_NCX. Simulated APs and ionic currents at 1-Hz pacing rate under HF conditions, 50% inhibition of I_Kr, 30% increase of I_CaL. Panel A shows EADs (dark line) with the basic HF value of I_NCX and APs with no EADs when I_NCX was 50% blocked (light line). The temporal evolutions of I_NaL (panel B), [Ca²⁺]_SR (panel C), I_rel (panel D), NCX activity (panel E), and activation gate of I_CaL (panel F) are also depicted with the basic HF value of I_NCX (dark line) and when I_NCX was 50% blocked (light line).

https://doi.org/10.1371/journal.pone.0032659.s014

(TIF)

Figure S15.

Mechanisms for early afterdepolarizations with a 50% transitory block of I_NCX. Simulated APs and ionic currents at 1-Hz pacing rate under HF conditions, 50% inhibition of I_Kr, 30% increase of I_CaL. Panel A shows EADs (dark line) with the basic HF value of I_NCX and APs with no EADs when I_NCX was 50% blocked (light line during the 5 stimulations indicated by a horizontal line. The temporal evolutions of I_NaL (panel B), [Ca²⁺]_SR (panel C), I_rel (panel D), NCX activity (panel E), and activation gate of I_CaL (panel F) are also depicted with the basic HF value of I_NCX (dark line) and when I_NCX was 50% blocked (light line) during the 5 stimulations indicated by a horizontal line.

https://doi.org/10.1371/journal.pone.0032659.s015

(TIF)

Author Contributions

Conceived and designed the experiments: BT. Performed the experiments: KC JFG BT. Analyzed the data: KC JFG BT LB SR JS JMF. Wrote the paper: BT SR LB.

References

1. Rosamond W, Flegal K, Friday G, Furie K, Go A, et al. (2007) Heart disease and stroke statistics–2007 update: a report from the American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Circulation 115: e69–171.
- View Article
- Google Scholar
2. Janse MJ (2004) Electrophysiological changes in heart failure and their relationship to arrhythmogenesis. Cardiovasc Res 61: 208–217.
- View Article
- Google Scholar
3. Tomaselli GF, Marban E (1999) Electrophysiological remodeling in hypertrophy and heart failure. Cardiovasc Res 42: 270–283.
- View Article
- Google Scholar
4. Tomaselli GF, Zipes DP (2004) What causes sudden death in heart failure? Circ Res 95: 754–763.
- View Article
- Google Scholar
5. Priebe L, Beuckelmann DJ (1998) Simulation study of cellular electric properties in heart failure. Circ Res 82: 1206–1223.
- View Article
- Google Scholar
6. Li GR, Lau CP, Ducharme A, Tardif JC, Nattel S (2002) Transmural action potential and ionic current remodeling in ventricles of failing canine hearts. Am J Physiol Heart Circ Physiol 283: H1031–H1041.
- View Article
- Google Scholar
7. Li GR, Lau CP, Leung TK, Nattel S (2004) Ionic current abnormalities associated with prolonged action potentials in cardiomyocytes from diseased human right ventricles. Heart Rhythm 1: 460–468.
- View Article
- Google Scholar
8. Li GR, Feng J, Yue L, Carrier M (1998) Transmural heterogeneity of action potentials and Ito1 in myocytes isolated from the human right ventricle. Am J Physiol 275: H369–H377.
- View Article
- Google Scholar
9. Bers DM, Despa S, Bossuyt J (2006) Regulation of Ca2+ and Na+ in normal and failing cardiac myocytes. Ann N Y Acad Sci 1080: 165–177.
- View Article
- Google Scholar
10. Beuckelmann DJ, Nabauer M, Erdmann E (1993) Alterations of K+ currents in isolated human ventricular myocytes from patients with terminal heart failure. Circ Res 73: 379–385.
- View Article
- Google Scholar
11. Piacentino V , Weber CR, Chen X, Weisser-Thomas J, Margulies KB, et al. (2003) Cellular basis of abnormal calcium transients of failing human ventricular myocytes. Circ Res 92: 651–658.
- View Article
- Google Scholar
12. Winslow RL, Rice J, Jafri S, Marban E, O'Rourke B (1999) Mechanisms of altered excitation-contraction coupling in canine tachycardia-induced heart failure, II: model studies. Circ Res 84: 571–586.
- View Article
- Google Scholar
13. Antoons G, Oros A, Bito V, Sipido KR, Vos MA (2007) Cellular basis for triggered ventricular arrhythmias that occur in the setting of compensated hypertrophy and heart failure: considerations for diagnosis and treatment. J Electrocardiol 40: S8–14.
- View Article
- Google Scholar
14. Bundgaard H, Kjeldsen K (1996) Human myocardial Na,K-ATPase concentration in heart failure. Mol Cell Biochem 163–164: 277–283.
- View Article
- Google Scholar
15. Maltsev VA, Silverman N, Sabbah HN, Undrovinas AI (2007) Chronic heart failure slows late sodium current in human and canine ventricular myocytes: implications for repolarization variability. Eur J Heart Fail 9: 219–227.
- View Article
- Google Scholar
16. Valdivia CR, Chu WW, Pu J, Foell JD, Haworth RA, et al. (2005) Increased late sodium current in myocytes from a canine heart failure model and from failing human heart. J Mol Cell Cardiol 38: 475–483.
- View Article
- Google Scholar
17. Zaza A, Belardinelli L, Shryock JC (2008) Pathophysiology and pharmacology of the cardiac “late sodium current.”. Pharmacol Ther 119: 326–339.
- View Article
- Google Scholar
18. Undrovinas AI, Belardinelli L, Undrovinas NA, Sabbah HN (2006) Ranolazine improves abnormal repolarization and contraction in left ventricular myocytes of dogs with heart failure by inhibiting late sodium current. J Cardiovasc Electrophysiol 17: Suppl 1S169–S177.
- View Article
- Google Scholar
19. Undrovinas NA, Maltsev VA, Belardinelli L, Sabbah HN, Undrovinas A (2010) Late sodium current contributes to diastolic cell Ca2+ accumulation in chronic heart failure. J Physiol Sci 60: 245–257.
- View Article
- Google Scholar
20. Milberg P, Pott C, Fink M, Frommeyer G, Matsuda T, et al. (2008) Inhibition of the Na+/Ca2+ exchanger suppresses torsades de pointes in an intact heart model of long QT syndrome-2 and long QT syndrome-3. Heart Rhythm 5: 1444–1452.
- View Article
- Google Scholar
21. Wu L, Guo D, Li H, Hackett J, Yan GX, et al. (2008) Role of late sodium current in modulating the proarrhythmic and antiarrhythmic effects of quinidine. Heart Rhythm 5: 1726–1734.
- View Article
- Google Scholar
22. Wu L, Rajamani S, Li H, January CT, Shryock JC, et al. (2009) Reduction of repolarization reserve unmasks the proarrhythmic role of endogenous late Na(+) current in the heart. Am J Physiol Heart Circ Physiol 297: H1048–H1057.
- View Article
- Google Scholar
23. Wu L, Ma J, Li H, Wang C, Grandi E, et al. (2011) Late Sodium Current Contributes to the Reverse Rate-Dependent Effect of IKr Inhibition on Ventricular Repolarization. Circulation 123: 1713–1720.
- View Article
- Google Scholar
24. Wu L, Shryock JC, Song Y, Li Y, Antzelevitch C, et al. (2004) Antiarrhythmic effects of ranolazine in a guinea pig in vitro model of long-QT syndrome. J Pharmacol Exp Ther 310: 599–605.
- View Article
- Google Scholar
25. Grandi E, Pasqualini FS, Bers DM (2010) A novel computational model of the human ventricular action potential and Ca transient. J Mol Cell Cardiol 48: 112–121.
- View Article
- Google Scholar
26. Maltsev VA, Sabbah HN, Higgins RS, Silverman N, Lesch M, et al. (1998) Novel, ultraslow inactivating sodium current in human ventricular cardiomyocytes. Circulation 98: 2545–2552.
- View Article
- Google Scholar
27. Drouin E, Charpentier F, Gauthier C, Laurent K, Le Marec H (1995) Electrophysiologic characteristics of cells spanning the left ventricular wall of human heart: evidence for presence of M cells. J Am Coll Cardiol 26: 185–192.
- View Article
- Google Scholar
28. Glukhov AV, Fedorov VV, Lou Q, Ravikumar VK, Kalish PW, et al. (2010) Transmural dispersion of repolarization in failing and nonfailing human ventricle. Circ Res 106: 981–991.
- View Article
- Google Scholar
29. So PP, Backx PH, Dorian P (2008) Slow delayed rectifier K+ current block by HMR 1556 increases dispersion of repolarization and promotes Torsades de Pointes in rabbit ventricles. Br J Pharmacol 155: 1185–1194.
- View Article
- Google Scholar
30. Grandi E, Puglisi JL, Wagner S, Maier LS, Severi S, et al. (2007) Simulation of Ca-calmodulin-dependent protein kinase II on rabbit ventricular myocyte ion currents and action potentials. Biophys J 93: 3835–3847.
- View Article
- Google Scholar
31. Weber CR, Piacentino V , Houser SR, Bers DM (2003) Dynamic regulation of sodium/calcium exchange function in human heart failure. Circulation 108: 2224–2229.
- View Article
- Google Scholar
32. Romero L, Pueyo E, Fink M, Rodriguez B (2009) Impact of ionic current variability on human ventricular cellular electrophysiology. Am J Physiol Heart Circ Physiol.
- View Article
- Google Scholar
33. Pieske B, Maier LS, Piacentino V 3rd, Weisser J, Hasenfuss G, et al. (2002) Rate dependence of [Na+]i and contractility in nonfailing and failing human myocardium. Circulation 106: 447–453.
- View Article
- Google Scholar
34. Carro J, Rodriguez JF, Laguna P, Pueyo E (2011) A human ventricular cell model for investigation of cardiac arrhythmias under hyperkalaemic conditions. Philos Transact A Math Phys Eng Sci 369: 4205–4232.
- View Article
- Google Scholar
35. Baartscheer A, Schumacher CA, Belterman CN, Coronel R, Fiolet JW (2003) [Na+]i and the driving force of the Na+/Ca2+-exchanger in heart failure. Cardiovasc Res 57: 986–995.
- View Article
- Google Scholar
36. Despa S, Islam MA, Weber CR, Pogwizd SM, Bers DM (2002) Intracellular Na(+) concentration is elevated in heart failure but Na/K pump function is unchanged. Circulation 105: 2543–2548.
- View Article
- Google Scholar
37. Levi AJ, Dalton GR, Hancox JC, Mitcheson JS, Issberner J, et al. (1997) Role of intracellular sodium overload in the genesis of cardiac arrhythmias. J Cardiovasc Electrophysiol 8: 700–721.
- View Article
- Google Scholar
38. Dorian P, Newman D (2000) Rate dependence of the effect of antiarrhythmic drugs delaying cardiac repolarization: an overview. Europace 2: 277–285.
- View Article
- Google Scholar
39. Zygmunt AC, Eddlestone GT, Thomas GP, Nesterenko VV, Antzelevitch C (2001) Larger late sodium conductance in M cells contributes to electrical heterogeneity in canine ventricle. Am J Physiol Heart Circ Physiol 281: H689–H697.
- View Article
- Google Scholar
40. Nagatomo T, January CT, Makielski JC (2000) Preferential block of late sodium current in the LQT3 DeltaKPQ mutant by the class I(C) antiarrhythmic flecainide. Mol Pharmacol 57: 101–107.
- View Article
- Google Scholar
41. Rajamani S, El-Bizri , Shryock JC, Makielski JC, Belardinelli L (2009) Use-dependent block of cardiac late Na current by ranolazine. Heart Rhythm.
- View Article
- Google Scholar
42. January CT, Riddle JM, Salata JJ (1988) A model for early afterdepolarizations: induction with the Ca2+ channel agonist Bay K 8644. Circ Res 62: 563–571.
- View Article
- Google Scholar
43. O'Hara T, Virag L, Varro A, Rudy Y (2011) Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS Comput Biol 7: e1002061.
- View Article
- Google Scholar
44. Zhang Y, Shou G, Xia L (2005) Simulation study of transmural cellular electrical properties in failed human heart. Conf Proc IEEE Eng Med Biol Soc 1: 337–340.
- View Article
- Google Scholar
45. Ten Tusscher KH, Noble D, Noble PJ, Panfilov AV (2004) A model for human ventricular tissue. Am J Physiol Heart Circ Physiol 286: H1573–H1589.
- View Article
- Google Scholar
46. Hund TJ, Rudy Y (2004) Rate dependence and regulation of action potential and calcium transient in a canine cardiac ventricular cell model. Circulation 110: 3168–3174.
- View Article
- Google Scholar
47. Shannon TR, Wang F, Bers DM (2005) Regulation of cardiac sarcoplasmic reticulum Ca release by luminal [Ca] and altered gating assessed with a mathematical model. Biophys J 89: 4096–4110.
- View Article
- Google Scholar
48. Ten Tusscher KH, Panfilov AV (2006) Alternans and spiral breakup in a human ventricular tissue model. Am J Physiol Heart Circ Physiol 291: H1088–H1100.
- View Article
- Google Scholar
49. von Lewinski D, Bisping E, Elgner A, Kockskamper J, Pieske B (2007) Mechanistic insight into the functional and toxic effects of Strophanthidin in the failing human myocardium. Eur J Heart Fail 9: 1086–1094.
- View Article
- Google Scholar
50. Pogwizd SM, Sipido KR, Verdonck F, Bers DM (2003) Intracellular Na in animal models of hypertrophy and heart failure: contractile function and arrhythmogenesis. Cardiovasc Res 57: 887–896.
- View Article
- Google Scholar
51. Wagner S, Ruff HM, Weber SL, Bellmann S, Sowa T, et al. (2011) Reactive oxygen species-activated Ca/Calmodulin Kinase IIdelta is required for late I(Na) augmentation leading to cellular Na and Ca overload. Circ Res 108: 555–565.
- View Article
- Google Scholar
52. Shannon TR, Wang F, Puglisi J, Weber C, Bers DM (2004) A mathematical treatment of integrated Ca dynamics within the ventricular myocyte. Biophys J 87: 3351–3371.
- View Article
- Google Scholar
53. Carmeliet E (2004) Intracellular Ca(2+) concentration and rate adaptation of the cardiac action potential. Cell Calcium 35: 557–573.
- View Article
- Google Scholar
54. Eisner DA, Dibb KM, Trafford AW (2009) The mechanism and significance of the slow changes of ventricular action potential duration following a change of heart rate. Exp Physiol 94: 520–528.
- View Article
- Google Scholar
55. Pueyo E, Husti Z, Hornyik T, Baczko I, Laguna P, et al. (2010) Mechanisms of ventricular rate adaptation as a predictor of arrhythmic risk. Am J Physiol Heart Circ Physiol 298: H1577–H1587.
- View Article
- Google Scholar
56. Faber GM, Rudy Y (2000) Action potential and contractility changes in [Na(+)](i) overloaded cardiac myocytes: a simulation study. Biophys J 78: 2392–2404.
- View Article
- Google Scholar
57. Biliczki P, Virag L, Iost N, Papp JG, Varro A (2002) Interaction of different potassium channels in cardiac repolarization in dog ventricular preparations: role of repolarization reserve. Br J Pharmacol 137: 361–368.
- View Article
- Google Scholar
58. Maltsev VA, Undrovinas AI (2006) A multi-modal composition of the late Na+ current in human ventricular cardiomyocytes. Cardiovasc Res 69: 116–127.
- View Article
- Google Scholar
59. Noble D, Varghese A, Kohl P, Noble P (1998) Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes. Can J Cardiol 14: 123–134.
- View Article
- Google Scholar
60. Cardona K, Trenor B, Rajamani S, Romero L, Ferrero JM, et al. (2010) Effects of Late Sodium Current Enhancement During LQT-relatedArrhythmias. A Simulation Study. IEEE-EMBS 3237–3240.
- View Article
- Google Scholar
61. Undrovinas AI, Maltsev VA, Kyle JW, Silverman N, Sabbah HN (2002) Gating of the late Na+ channel in normal and failing human myocardium. J Mol Cell Cardiol 34: 1477–1489.
- View Article
- Google Scholar
62. Antoons G, Oros A, Beekman JD, Engelen MA, Houtman MJ, et al. (2010) Late na(+) current inhibition by ranolazine reduces torsades de pointes in the chronic atrioventricular block dog model. J Am Coll Cardiol 55: 801–809.
- View Article
- Google Scholar
63. Hund TJ, Koval OM, Li J, Wright PJ, Qian L, et al. (2010) A beta(IV)-spectrin/CaMKII signaling complex is essential for membrane excitability in mice. J Clin Invest 120: 3508–3519.
- View Article
- Google Scholar
64. Maltsev VA, Sabbah HN, Undrovinas AI (2001) Late sodium current is a novel target for amiodarone: studies in failing human myocardium. J Mol Cell Cardiol 33: 923–932.
- View Article
- Google Scholar
65. Nagatomo T, Fan Z, Ye B, Tonkovich GS, January CT, et al. (1998) Temperature dependence of early and late currents in human cardiac wild-type and long Q-T DeltaKPQ Na+ channels. Am J Physiol 275: H2016–H2024.
- View Article
- Google Scholar
66. Beuckelmann DJ, Nabauer M, Erdmann E (1992) Intracellular calcium handling in isolated ventricular myocytes from patients with terminal heart failure. Circulation 85: 1046–1055.
- View Article
- Google Scholar
67. Undrovinas AI, Maltsev VA, Sabbah HN (1999) Repolarization abnormalities in cardiomyocytes of dogs with chronic heart failure: role of sustained inward current. Cell Mol Life Sci 55: 494–505.
- View Article
- Google Scholar
68. Reinecke H, Studer R, Vetter R, Holtz J, Drexler H (1996) Cardiac Na+/Ca2+ exchange activity in patients with end-stage heart failure. Cardiovasc Res 31: 48–54.
- View Article
- Google Scholar
69. Curran J, Brown KH, Santiago DJ, Pogwizd S, Bers DM, et al. (2010) Spontaneous Ca waves in ventricular myocytes from failing hearts depend on Ca(2+)-calmodulin-dependent protein kinase II. J Mol Cell Cardiol 49: 25–32.
- View Article
- Google Scholar
70. Wettwer E, Amos GJ, Posival H, Ravens U (1994) Transient outward current in human ventricular myocytes of subepicardial and subendocardial origin. Circ Res 75: 473–482.
- View Article
- Google Scholar
71. Nabauer M, Beuckelmann DJ, Uberfuhr P, Steinbeck G (1996) Regional differences in current density and rate-dependent properties of the transient outward current in subepicardial and subendocardial myocytes of human left ventricle. Circulation 93: 168–177.
- View Article
- Google Scholar
72. Hasenfuss G, Reinecke H, Studer R, Meyer M, Pieske B, et al. (1994) Relation between myocardial function and expression of sarcoplasmic reticulum Ca(2+)-ATPase in failing and nonfailing human myocardium. Circ Res 75: 434–442.
- View Article
- Google Scholar
73. Schwinger RH, Bohm M, Schmidt U, Karczewski P, Bavendiek U, et al. (1995) Unchanged protein levels of SERCA II and phospholamban but reduced Ca2+ uptake and Ca(2+)-ATPase activity of cardiac sarcoplasmic reticulum from dilated cardiomyopathy patients compared with patients with nonfailing hearts. Circulation 92: 3220–3228.
- View Article
- Google Scholar

[ref1] 1. Rosamond W, Flegal K, Friday G, Furie K, Go A, et al. (2007) Heart disease and stroke statistics–2007 update: a report from the American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Circulation 115: e69–171.
View Article
Google Scholar

[2] View Article

[3] Google Scholar

[ref2] 2. Janse MJ (2004) Electrophysiological changes in heart failure and their relationship to arrhythmogenesis. Cardiovasc Res 61: 208–217.
View Article
Google Scholar

[5] View Article

[6] Google Scholar

[ref3] 3. Tomaselli GF, Marban E (1999) Electrophysiological remodeling in hypertrophy and heart failure. Cardiovasc Res 42: 270–283.
View Article
Google Scholar

[8] View Article

[9] Google Scholar

[ref4] 4. Tomaselli GF, Zipes DP (2004) What causes sudden death in heart failure? Circ Res 95: 754–763.
View Article
Google Scholar

[11] View Article

[12] Google Scholar

[ref5] 5. Priebe L, Beuckelmann DJ (1998) Simulation study of cellular electric properties in heart failure. Circ Res 82: 1206–1223.
View Article
Google Scholar

[14] View Article

[15] Google Scholar

[ref6] 6. Li GR, Lau CP, Ducharme A, Tardif JC, Nattel S (2002) Transmural action potential and ionic current remodeling in ventricles of failing canine hearts. Am J Physiol Heart Circ Physiol 283: H1031–H1041.
View Article
Google Scholar

[17] View Article

[18] Google Scholar

[ref7] 7. Li GR, Lau CP, Leung TK, Nattel S (2004) Ionic current abnormalities associated with prolonged action potentials in cardiomyocytes from diseased human right ventricles. Heart Rhythm 1: 460–468.
View Article
Google Scholar

[20] View Article

[21] Google Scholar

[ref8] 8. Li GR, Feng J, Yue L, Carrier M (1998) Transmural heterogeneity of action potentials and Ito1 in myocytes isolated from the human right ventricle. Am J Physiol 275: H369–H377.
View Article
Google Scholar

[23] View Article

[24] Google Scholar

[ref9] 9. Bers DM, Despa S, Bossuyt J (2006) Regulation of Ca2+ and Na+ in normal and failing cardiac myocytes. Ann N Y Acad Sci 1080: 165–177.
View Article
Google Scholar

[26] View Article

[27] Google Scholar

[ref10] 10. Beuckelmann DJ, Nabauer M, Erdmann E (1993) Alterations of K+ currents in isolated human ventricular myocytes from patients with terminal heart failure. Circ Res 73: 379–385.
View Article
Google Scholar

[29] View Article

[30] Google Scholar

[ref11] 11. Piacentino V , Weber CR, Chen X, Weisser-Thomas J, Margulies KB, et al. (2003) Cellular basis of abnormal calcium transients of failing human ventricular myocytes. Circ Res 92: 651–658.
View Article
Google Scholar

[32] View Article

[33] Google Scholar

[ref12] 12. Winslow RL, Rice J, Jafri S, Marban E, O'Rourke B (1999) Mechanisms of altered excitation-contraction coupling in canine tachycardia-induced heart failure, II: model studies. Circ Res 84: 571–586.
View Article
Google Scholar

[35] View Article

[36] Google Scholar

[ref13] 13. Antoons G, Oros A, Bito V, Sipido KR, Vos MA (2007) Cellular basis for triggered ventricular arrhythmias that occur in the setting of compensated hypertrophy and heart failure: considerations for diagnosis and treatment. J Electrocardiol 40: S8–14.
View Article
Google Scholar

[38] View Article

[39] Google Scholar

[ref14] 14. Bundgaard H, Kjeldsen K (1996) Human myocardial Na,K-ATPase concentration in heart failure. Mol Cell Biochem 163–164: 277–283.
View Article
Google Scholar

[41] View Article

[42] Google Scholar

[ref15] 15. Maltsev VA, Silverman N, Sabbah HN, Undrovinas AI (2007) Chronic heart failure slows late sodium current in human and canine ventricular myocytes: implications for repolarization variability. Eur J Heart Fail 9: 219–227.
View Article
Google Scholar

[44] View Article

[45] Google Scholar

[ref16] 16. Valdivia CR, Chu WW, Pu J, Foell JD, Haworth RA, et al. (2005) Increased late sodium current in myocytes from a canine heart failure model and from failing human heart. J Mol Cell Cardiol 38: 475–483.
View Article
Google Scholar

[47] View Article

[48] Google Scholar

[ref17] 17. Zaza A, Belardinelli L, Shryock JC (2008) Pathophysiology and pharmacology of the cardiac “late sodium current.”. Pharmacol Ther 119: 326–339.
View Article
Google Scholar

[50] View Article

[51] Google Scholar

[ref18] 18. Undrovinas AI, Belardinelli L, Undrovinas NA, Sabbah HN (2006) Ranolazine improves abnormal repolarization and contraction in left ventricular myocytes of dogs with heart failure by inhibiting late sodium current. J Cardiovasc Electrophysiol 17: Suppl 1S169–S177.
View Article
Google Scholar

[53] View Article

[54] Google Scholar

[ref19] 19. Undrovinas NA, Maltsev VA, Belardinelli L, Sabbah HN, Undrovinas A (2010) Late sodium current contributes to diastolic cell Ca2+ accumulation in chronic heart failure. J Physiol Sci 60: 245–257.
View Article
Google Scholar

[56] View Article

[57] Google Scholar

[ref20] 20. Milberg P, Pott C, Fink M, Frommeyer G, Matsuda T, et al. (2008) Inhibition of the Na+/Ca2+ exchanger suppresses torsades de pointes in an intact heart model of long QT syndrome-2 and long QT syndrome-3. Heart Rhythm 5: 1444–1452.
View Article
Google Scholar

[59] View Article

[60] Google Scholar

[ref21] 21. Wu L, Guo D, Li H, Hackett J, Yan GX, et al. (2008) Role of late sodium current in modulating the proarrhythmic and antiarrhythmic effects of quinidine. Heart Rhythm 5: 1726–1734.
View Article
Google Scholar

[62] View Article

[63] Google Scholar

[ref22] 22. Wu L, Rajamani S, Li H, January CT, Shryock JC, et al. (2009) Reduction of repolarization reserve unmasks the proarrhythmic role of endogenous late Na(+) current in the heart. Am J Physiol Heart Circ Physiol 297: H1048–H1057.
View Article
Google Scholar

[65] View Article

[66] Google Scholar

[ref23] 23. Wu L, Ma J, Li H, Wang C, Grandi E, et al. (2011) Late Sodium Current Contributes to the Reverse Rate-Dependent Effect of IKr Inhibition on Ventricular Repolarization. Circulation 123: 1713–1720.
View Article
Google Scholar

[68] View Article

[69] Google Scholar

[ref24] 24. Wu L, Shryock JC, Song Y, Li Y, Antzelevitch C, et al. (2004) Antiarrhythmic effects of ranolazine in a guinea pig in vitro model of long-QT syndrome. J Pharmacol Exp Ther 310: 599–605.
View Article
Google Scholar

[71] View Article

[72] Google Scholar

[ref25] 25. Grandi E, Pasqualini FS, Bers DM (2010) A novel computational model of the human ventricular action potential and Ca transient. J Mol Cell Cardiol 48: 112–121.
View Article
Google Scholar

[74] View Article

[75] Google Scholar

[ref26] 26. Maltsev VA, Sabbah HN, Higgins RS, Silverman N, Lesch M, et al. (1998) Novel, ultraslow inactivating sodium current in human ventricular cardiomyocytes. Circulation 98: 2545–2552.
View Article
Google Scholar

[77] View Article

[78] Google Scholar

[ref27] 27. Drouin E, Charpentier F, Gauthier C, Laurent K, Le Marec H (1995) Electrophysiologic characteristics of cells spanning the left ventricular wall of human heart: evidence for presence of M cells. J Am Coll Cardiol 26: 185–192.
View Article
Google Scholar

[80] View Article

[81] Google Scholar

[ref28] 28. Glukhov AV, Fedorov VV, Lou Q, Ravikumar VK, Kalish PW, et al. (2010) Transmural dispersion of repolarization in failing and nonfailing human ventricle. Circ Res 106: 981–991.
View Article
Google Scholar

[83] View Article

[84] Google Scholar

[ref29] 29. So PP, Backx PH, Dorian P (2008) Slow delayed rectifier K+ current block by HMR 1556 increases dispersion of repolarization and promotes Torsades de Pointes in rabbit ventricles. Br J Pharmacol 155: 1185–1194.
View Article
Google Scholar

[86] View Article

[87] Google Scholar

[ref30] 30. Grandi E, Puglisi JL, Wagner S, Maier LS, Severi S, et al. (2007) Simulation of Ca-calmodulin-dependent protein kinase II on rabbit ventricular myocyte ion currents and action potentials. Biophys J 93: 3835–3847.
View Article
Google Scholar

[89] View Article

[90] Google Scholar

[ref31] 31. Weber CR, Piacentino V , Houser SR, Bers DM (2003) Dynamic regulation of sodium/calcium exchange function in human heart failure. Circulation 108: 2224–2229.
View Article
Google Scholar

[92] View Article

[93] Google Scholar

[ref32] 32. Romero L, Pueyo E, Fink M, Rodriguez B (2009) Impact of ionic current variability on human ventricular cellular electrophysiology. Am J Physiol Heart Circ Physiol.
View Article
Google Scholar

[95] View Article

[96] Google Scholar

[ref33] 33. Pieske B, Maier LS, Piacentino V 3rd, Weisser J, Hasenfuss G, et al. (2002) Rate dependence of [Na+]i and contractility in nonfailing and failing human myocardium. Circulation 106: 447–453.
View Article
Google Scholar

[98] View Article

[99] Google Scholar

[ref34] 34. Carro J, Rodriguez JF, Laguna P, Pueyo E (2011) A human ventricular cell model for investigation of cardiac arrhythmias under hyperkalaemic conditions. Philos Transact A Math Phys Eng Sci 369: 4205–4232.
View Article
Google Scholar

[101] View Article

[102] Google Scholar

[ref35] 35. Baartscheer A, Schumacher CA, Belterman CN, Coronel R, Fiolet JW (2003) [Na+]i and the driving force of the Na+/Ca2+-exchanger in heart failure. Cardiovasc Res 57: 986–995.
View Article
Google Scholar

[104] View Article

[105] Google Scholar

[ref36] 36. Despa S, Islam MA, Weber CR, Pogwizd SM, Bers DM (2002) Intracellular Na(+) concentration is elevated in heart failure but Na/K pump function is unchanged. Circulation 105: 2543–2548.
View Article
Google Scholar

[107] View Article

[108] Google Scholar

[ref37] 37. Levi AJ, Dalton GR, Hancox JC, Mitcheson JS, Issberner J, et al. (1997) Role of intracellular sodium overload in the genesis of cardiac arrhythmias. J Cardiovasc Electrophysiol 8: 700–721.
View Article
Google Scholar

[110] View Article

[111] Google Scholar

[ref38] 38. Dorian P, Newman D (2000) Rate dependence of the effect of antiarrhythmic drugs delaying cardiac repolarization: an overview. Europace 2: 277–285.
View Article
Google Scholar

[113] View Article

[114] Google Scholar

[ref39] 39. Zygmunt AC, Eddlestone GT, Thomas GP, Nesterenko VV, Antzelevitch C (2001) Larger late sodium conductance in M cells contributes to electrical heterogeneity in canine ventricle. Am J Physiol Heart Circ Physiol 281: H689–H697.
View Article
Google Scholar

[116] View Article

[117] Google Scholar

[ref40] 40. Nagatomo T, January CT, Makielski JC (2000) Preferential block of late sodium current in the LQT3 DeltaKPQ mutant by the class I(C) antiarrhythmic flecainide. Mol Pharmacol 57: 101–107.
View Article
Google Scholar

[119] View Article

[120] Google Scholar

[ref41] 41. Rajamani S, El-Bizri , Shryock JC, Makielski JC, Belardinelli L (2009) Use-dependent block of cardiac late Na current by ranolazine. Heart Rhythm.
View Article
Google Scholar

[122] View Article

[123] Google Scholar

[ref42] 42. January CT, Riddle JM, Salata JJ (1988) A model for early afterdepolarizations: induction with the Ca2+ channel agonist Bay K 8644. Circ Res 62: 563–571.
View Article
Google Scholar

[125] View Article

[126] Google Scholar

[ref43] 43. O'Hara T, Virag L, Varro A, Rudy Y (2011) Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS Comput Biol 7: e1002061.
View Article
Google Scholar

[128] View Article

[129] Google Scholar

[ref44] 44. Zhang Y, Shou G, Xia L (2005) Simulation study of transmural cellular electrical properties in failed human heart. Conf Proc IEEE Eng Med Biol Soc 1: 337–340.
View Article
Google Scholar

[131] View Article

[132] Google Scholar

[ref45] 45. Ten Tusscher KH, Noble D, Noble PJ, Panfilov AV (2004) A model for human ventricular tissue. Am J Physiol Heart Circ Physiol 286: H1573–H1589.
View Article
Google Scholar

[134] View Article

[135] Google Scholar

[ref46] 46. Hund TJ, Rudy Y (2004) Rate dependence and regulation of action potential and calcium transient in a canine cardiac ventricular cell model. Circulation 110: 3168–3174.
View Article
Google Scholar

[137] View Article

[138] Google Scholar

[ref47] 47. Shannon TR, Wang F, Bers DM (2005) Regulation of cardiac sarcoplasmic reticulum Ca release by luminal [Ca] and altered gating assessed with a mathematical model. Biophys J 89: 4096–4110.
View Article
Google Scholar

[140] View Article

[141] Google Scholar

[ref48] 48. Ten Tusscher KH, Panfilov AV (2006) Alternans and spiral breakup in a human ventricular tissue model. Am J Physiol Heart Circ Physiol 291: H1088–H1100.
View Article
Google Scholar

[143] View Article

[144] Google Scholar

[ref49] 49. von Lewinski D, Bisping E, Elgner A, Kockskamper J, Pieske B (2007) Mechanistic insight into the functional and toxic effects of Strophanthidin in the failing human myocardium. Eur J Heart Fail 9: 1086–1094.
View Article
Google Scholar

[146] View Article

[147] Google Scholar

[ref50] 50. Pogwizd SM, Sipido KR, Verdonck F, Bers DM (2003) Intracellular Na in animal models of hypertrophy and heart failure: contractile function and arrhythmogenesis. Cardiovasc Res 57: 887–896.
View Article
Google Scholar

[149] View Article

[150] Google Scholar

[ref51] 51. Wagner S, Ruff HM, Weber SL, Bellmann S, Sowa T, et al. (2011) Reactive oxygen species-activated Ca/Calmodulin Kinase IIdelta is required for late I(Na) augmentation leading to cellular Na and Ca overload. Circ Res 108: 555–565.
View Article
Google Scholar

[152] View Article

[153] Google Scholar

[ref52] 52. Shannon TR, Wang F, Puglisi J, Weber C, Bers DM (2004) A mathematical treatment of integrated Ca dynamics within the ventricular myocyte. Biophys J 87: 3351–3371.
View Article
Google Scholar

[155] View Article

[156] Google Scholar

[ref53] 53. Carmeliet E (2004) Intracellular Ca(2+) concentration and rate adaptation of the cardiac action potential. Cell Calcium 35: 557–573.
View Article
Google Scholar

[158] View Article

[159] Google Scholar

[ref54] 54. Eisner DA, Dibb KM, Trafford AW (2009) The mechanism and significance of the slow changes of ventricular action potential duration following a change of heart rate. Exp Physiol 94: 520–528.
View Article
Google Scholar

[161] View Article

[162] Google Scholar

[ref55] 55. Pueyo E, Husti Z, Hornyik T, Baczko I, Laguna P, et al. (2010) Mechanisms of ventricular rate adaptation as a predictor of arrhythmic risk. Am J Physiol Heart Circ Physiol 298: H1577–H1587.
View Article
Google Scholar

[164] View Article

[165] Google Scholar

[ref56] 56. Faber GM, Rudy Y (2000) Action potential and contractility changes in [Na(+)](i) overloaded cardiac myocytes: a simulation study. Biophys J 78: 2392–2404.
View Article
Google Scholar

[167] View Article

[168] Google Scholar

[ref57] 57. Biliczki P, Virag L, Iost N, Papp JG, Varro A (2002) Interaction of different potassium channels in cardiac repolarization in dog ventricular preparations: role of repolarization reserve. Br J Pharmacol 137: 361–368.
View Article
Google Scholar

[170] View Article

[171] Google Scholar

[ref58] 58. Maltsev VA, Undrovinas AI (2006) A multi-modal composition of the late Na+ current in human ventricular cardiomyocytes. Cardiovasc Res 69: 116–127.
View Article
Google Scholar

[173] View Article

[174] Google Scholar

[ref59] 59. Noble D, Varghese A, Kohl P, Noble P (1998) Improved guinea-pig ventricular cell model incorporating a diadic space, IKr and IKs, and length- and tension-dependent processes. Can J Cardiol 14: 123–134.
View Article
Google Scholar

[176] View Article

[177] Google Scholar

[ref60] 60. Cardona K, Trenor B, Rajamani S, Romero L, Ferrero JM, et al. (2010) Effects of Late Sodium Current Enhancement During LQT-relatedArrhythmias. A Simulation Study. IEEE-EMBS 3237–3240.
View Article
Google Scholar

[179] View Article

[180] Google Scholar

[ref61] 61. Undrovinas AI, Maltsev VA, Kyle JW, Silverman N, Sabbah HN (2002) Gating of the late Na+ channel in normal and failing human myocardium. J Mol Cell Cardiol 34: 1477–1489.
View Article
Google Scholar

[182] View Article

[183] Google Scholar

[ref62] 62. Antoons G, Oros A, Beekman JD, Engelen MA, Houtman MJ, et al. (2010) Late na(+) current inhibition by ranolazine reduces torsades de pointes in the chronic atrioventricular block dog model. J Am Coll Cardiol 55: 801–809.
View Article
Google Scholar

[185] View Article

[186] Google Scholar

[ref63] 63. Hund TJ, Koval OM, Li J, Wright PJ, Qian L, et al. (2010) A beta(IV)-spectrin/CaMKII signaling complex is essential for membrane excitability in mice. J Clin Invest 120: 3508–3519.
View Article
Google Scholar

[188] View Article

[189] Google Scholar

[ref64] 64. Maltsev VA, Sabbah HN, Undrovinas AI (2001) Late sodium current is a novel target for amiodarone: studies in failing human myocardium. J Mol Cell Cardiol 33: 923–932.
View Article
Google Scholar

[191] View Article

[192] Google Scholar

[ref65] 65. Nagatomo T, Fan Z, Ye B, Tonkovich GS, January CT, et al. (1998) Temperature dependence of early and late currents in human cardiac wild-type and long Q-T DeltaKPQ Na+ channels. Am J Physiol 275: H2016–H2024.
View Article
Google Scholar

[194] View Article

[195] Google Scholar

[ref66] 66. Beuckelmann DJ, Nabauer M, Erdmann E (1992) Intracellular calcium handling in isolated ventricular myocytes from patients with terminal heart failure. Circulation 85: 1046–1055.
View Article
Google Scholar

[197] View Article

[198] Google Scholar

[ref67] 67. Undrovinas AI, Maltsev VA, Sabbah HN (1999) Repolarization abnormalities in cardiomyocytes of dogs with chronic heart failure: role of sustained inward current. Cell Mol Life Sci 55: 494–505.
View Article
Google Scholar

[200] View Article

[201] Google Scholar

[ref68] 68. Reinecke H, Studer R, Vetter R, Holtz J, Drexler H (1996) Cardiac Na+/Ca2+ exchange activity in patients with end-stage heart failure. Cardiovasc Res 31: 48–54.
View Article
Google Scholar

[203] View Article

[204] Google Scholar

[ref69] 69. Curran J, Brown KH, Santiago DJ, Pogwizd S, Bers DM, et al. (2010) Spontaneous Ca waves in ventricular myocytes from failing hearts depend on Ca(2+)-calmodulin-dependent protein kinase II. J Mol Cell Cardiol 49: 25–32.
View Article
Google Scholar

[206] View Article

[207] Google Scholar

[ref70] 70. Wettwer E, Amos GJ, Posival H, Ravens U (1994) Transient outward current in human ventricular myocytes of subepicardial and subendocardial origin. Circ Res 75: 473–482.
View Article
Google Scholar

[209] View Article

[210] Google Scholar

[ref71] 71. Nabauer M, Beuckelmann DJ, Uberfuhr P, Steinbeck G (1996) Regional differences in current density and rate-dependent properties of the transient outward current in subepicardial and subendocardial myocytes of human left ventricle. Circulation 93: 168–177.
View Article
Google Scholar

[212] View Article

[213] Google Scholar

[ref72] 72. Hasenfuss G, Reinecke H, Studer R, Meyer M, Pieske B, et al. (1994) Relation between myocardial function and expression of sarcoplasmic reticulum Ca(2+)-ATPase in failing and nonfailing human myocardium. Circ Res 75: 434–442.
View Article
Google Scholar

[215] View Article

[216] Google Scholar

[ref73] 73. Schwinger RH, Bohm M, Schmidt U, Karczewski P, Bavendiek U, et al. (1995) Unchanged protein levels of SERCA II and phospholamban but reduced Ca2+ uptake and Ca(2+)-ATPase activity of cardiac sarcoplasmic reticulum from dilated cardiomyopathy patients compared with patients with nonfailing hearts. Circulation 92: 3220–3228.
View Article
Google Scholar

[218] View Article

[219] Google Scholar

Figures

Abstract

Introduction

Results

Effects of INaL on human AP

Heart Failure model and sensitivity analysis

APD rate-dependence and the role of INaL in HF

Arrhythmogenic effects of INaL

Discussion

Major Findings

Simulation of heart failure

Important role of INaL during heart failure and clinical implications

Limitations of the study

Conclusion

Methods

Model of the human ventricular INaL

Heart failure cellular model for human

Sensitivity analysis of the HF model

Stimulation protocols

Supporting Information

Author Contributions

References

Effects of I_NaL on human AP

APD rate-dependence and the role of I_NaL in HF

Arrhythmogenic effects of I_NaL

Important role of I_NaL during heart failure and clinical implications

Model of the human ventricular I_NaL