Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Mechanistic investigation of Ca2+ alternans in human heart failure and its modulation by fibroblasts

  • Maria T. Mora,

    Roles Formal analysis, Methodology, Visualization, Writing – original draft

    Affiliation Centro de Insvestigación e Innovación en Bioingeniería, Universitat Politècnica de València, Valencia, Spain

  • Juan F. Gomez,

    Roles Methodology, Writing – review & editing

    Affiliation Centro de Insvestigación e Innovación en Bioingeniería, Universitat Politècnica de València, Valencia, Spain

  • Gregory Morley,

    Roles Writing – review & editing

    Affiliation Leon H. Charney Division of Cardiology, New York University School of Medicine, New York, New York, United States of America

  • Jose M. Ferrero ,

    Contributed equally to this work with: Jose M. Ferrero, Beatriz Trenor

    Roles Conceptualization, Formal analysis, Methodology, Supervision, Writing – review & editing

    Affiliation Centro de Insvestigación e Innovación en Bioingeniería, Universitat Politècnica de València, Valencia, Spain

  • Beatriz Trenor

    Contributed equally to this work with: Jose M. Ferrero, Beatriz Trenor

    Roles Conceptualization, Formal analysis, Project administration, Supervision, Writing – review & editing

    btrenor@eln.upv.es

    Affiliation Centro de Insvestigación e Innovación en Bioingeniería, Universitat Politècnica de València, Valencia, Spain

Abstract

Background

Heart failure (HF) is characterized, among other factors, by a progressive loss of contractile function and by the formation of an arrhythmogenic substrate, both aspects partially related to intracellular Ca2+ cycling disorders. In failing hearts both electrophysiological and structural remodeling, including fibroblast proliferation, contribute to changes in Ca2+ handling which promote the appearance of Ca2+ alternans (Ca-alt). Ca-alt in turn give rise to repolarization alternans, which promote dispersion of repolarization and contribute to reentrant activity. The computational analysis of the incidence of Ca2+ and/or repolarization alternans under HF conditions in the presence of fibroblasts could provide a better understanding of the mechanisms leading to HF arrhythmias and contractile function disorders.

Methods and findings

The goal of the present study was to investigate in silico the mechanisms leading to the formation of Ca-alt in failing human ventricular myocytes and tissues with disperse fibroblast distributions. The contribution of ionic currents variability to alternans formation at the cellular level was analyzed and the results show that in normal ventricular tissue, altered Ca2+ dynamics lead to Ca-alt, which precede APD alternans and can be aggravated by the presence of fibroblasts. Electrophysiological remodeling of failing tissue alone is sufficient to develop alternans. The incidence of alternans is reduced when fibroblasts are present in failing tissue due to significantly depressed Ca2+ transients. The analysis of the underlying ionic mechanisms suggests that Ca-alt are driven by Ca2+-handling protein and Ca2+ cycling dysfunctions in the junctional sarcoplasmic reticulum and that their contribution to alternans occurrence depends on the cardiac remodeling conditions and on myocyte-fibroblast interactions.

Conclusion

It can thus be concluded that fibroblasts modulate the formation of Ca-alt in human ventricular tissue subjected to heart failure-related electrophysiological remodeling. Pharmacological therapies should thus consider the extent of both the electrophysiological and structural remodeling present in the failing heart.

Introduction

Heart failure (HF) with reduced left ventricular ejection fraction is characterized by a progressive loss of contractile function and the genesis of malignant arrhythmias, leading to cardiac dysfunction and sudden cardiac death. Altered action potentials (APs), Ca2+ transients (CaTs) in myocytes and abnormal electrical conduction in cardiac tissue may result from the electrophysiological and structural remodeling of the failing heart, in which fibroblast proliferation has been observed [13].

Despite the extensive research on heterocellular connections between cardiac myocytes and non-excitable cells, the existence of this heterocellular coupling in vivo is still controversial [4]. The first experimental evidence was found in isolated cells and cultures in animal models and showed functional gap junctions between fibroblasts and myocytes, allowing electrical propagation [58]. More recent studies have demonstrated electrical coupling between scar tissue and the surrounding myocardium in live injured hearts [911]. These electrotonic interactions affect myocyte impulse conduction and electrophysiological activity [1215] and increase the arrhythmogenic consequences [16,17]. Although numerous studies have used computer simulations of cardiac tissue with fibroblast coupling to investigate the origin of reentrant activity [18,19], including our previous work on failing tissue [20], the role of beat-to-beat alternans has not yet been considered.

Mechanoelectrical alternans, i.e. T-wave and contractile beat-to-beat oscillations, have been observed in HF patients [21]. Electrical alternans are linked to arrhythmogenesis because of dispersion of refractoriness and have become a risk indicator of sudden cardiac death [22], in addition to the ventricular dysfunction resulting from mechanical alternans. Cardiac alternans manifest as AP and Ca2+ cycling fluctuations in myocytes due to the fact that the membrane voltage and Ca2+ processes are coupled. However, the underlying mechanisms leading to alternans have not yet been fully identified and the feedback nature of coupling between voltage and Ca2+ dynamics complicates the task. Among the proposed mechanisms, voltage-driven Ca2+-mediated alternans and the contribution of both systems [2327] have been suggested.

The most extensively developed theory is that alternans are induced by Ca2+ cycling instabilities, by which an altered sarcoplasmic reticulum (SR) Ca2+ release process through the ryanodine receptor (RyR) impairs Ca2+ dynamics and initiates alternans [28,29]. Fluctuations in the SR Ca2+ content [30] have suggested an impairment beween Ca2+ uptake and release mechanisms. Computational studies have analyzed the role of the different Ca2+-handling proteins and their interactions that contribute to Ca2+ alternans [31,32]. In agreement with computational results, SERCA overexpression has been shown to supress cellular alternans [33].

Since the normal heart can develop alternans at rapid heart rates, the vast majority of experimental and computational studies on cellular alternans [28,31,34] have been performed on non-failing myocytes or without considering all the main HF features, such as electrophysiological remodeling, intercellular decoupling and fibroblast coupling. These changes in the cardiac tissue modify myocyte electrophysiology and could alter the susceptibility to alternans. We consider that the passive characteristics of the tissue also play a role in generating alternans. Discordant alternans can produce spatial repolarization gradients that may generate electrical reentries leading to ventricular fibrillation [35], so that it is crucial to determine the causes and susceptability to spatially discordant versus concordant alternans.

The aim of this study was thus to characterize the onset of alternans in HF under the influence of disperse fibroblasts using a virtual two-dimensional human model. In the first part of the study, we used a mathematical model of a single human ventricular myocyte to study its bioelectrical activity, including Ca2+ dynamics in health and during HF, coupled or uncoupled to fibroblasts. Simulations were performed to determine the influence of pacing frequency on the alternans occurrence threshold. In the second part we simulated the activity of a virtual ventricular tissue with normal and failing myocyte populations, with or without inserted fibroblasts. We built populations of models by imposing variability to ion channels and pumps. The biological variability represented in the different models can reveal the most favorable conditions that induce Ca2+ beat-to-beat alternans and, by means of the Ca2+-driven theory, the important role of Ca2+-handling proteins in modifying Ca2+ dynamics.

Methods

Cellular simulations

Cardiac tissue consists of interconnected cardiomyocytes and other non-muscular cells, such as fibroblasts, which might establish electrical connections with myocytes through gap junctions. The role of fibroblast-myocyte interactions at the onset of alternans was first studied in single-cell simulations and compared to two coupled myocytes to identify the key parameters responsible for beat-to-beat oscillations, as well as the impact of HF. The reduction of the variables at the cellular level made it possible to reduce the computational time in the 2D study. The second study was essential in order to understand the origin, features and effects of alternans, such as arrhythmogenic spatial discordance.

The O’Hara et al. [36] human AP model (ORd) was used to simulate the electrophysiological activity of endocardial myocytes. As in our previous work [37], small changes were made to the sodium current (INa) in the model to improve electrical conduction and update the late component to the most recent experimental evidence. The active fibroblast model developed by MacCannell et al. [38] was coupled to myocytes to analyze fibroblast-myocyte interactions. Myocyte-fibroblast coupling was performed taking into account an individual gap junction conductance (Ggap) of 3 nS, within the experimental range [8], which was increased five-fold to establish a 1:5 ratio between cells, accounting for an increased degree of disperse fibroblasts. Although cell-to-cell coupling between myocytes is stronger, the experimental measurements in myocyte pairs span a wide range [39,40], so that the Ggap was set to an intermediate value of 1.3 μS. HF conditions, i.e. electrophysiological myocyte remodeling, were also simulated. In HF, myocyte-myocyte coupling was reduced by 50%.

Electrophysiological remodeling in HF was simulated by applying a scale factor to the maximal conductances of the selected ionic channels (Table 1), which represents an upregulation or downregulation of protein expression and/or function obtained from experimental human myocytes when possible, as detailed in Gomez et al. [41].

thumbnail
Table 1. Electrophysiological remodeling in HF applied to the basic ORd model.

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

Alternans typically occur at fast pacing rates. The pacing cycle length (PCL) threshold was determined by computing a restitution curve starting at a PCL of 1000 ms and reducing the cycle length until loss of capture was observed. Each simulation was run until steady-state was reached. The last two APs and CaTs were used to measure AP duration from maximal upstroke to 90% of repolarization (APD90) and the maximum intracellular [Ca2+], respectively.

After determining pacing threshold in the baseline models, populations of models were built and simulations were carried out at a fixed low PCL in which alternans tend to occur. Populations of models consisted of different cellular models built after varying the selected parameters related to ionic currents of the baseline model, mimicking biological inter-individual variability. We modified 13 maximal ion channel conductances and fluxes simultaneously using random scale factors obtained from a normal distribution of mean equal to 1 and standard deviation (σ) equal to 0.15, so that 95% (±2σ) of the parameters varied in the range ±30% of their baseline values. The selected variables affect the conductances of the fast Na+ current (INa), late Na+ current (INaL), transient outward K+ current (Ito), L-type Ca2+ current (ICaL), rapid delayed rectifier K+ current (IKr), slow delayed rectifier K+ current (IKs), inward rectifier K+ current (IK1), Na+/ K+ ATPase current (INaK), Na+/Ca2+ exchange current (INCX), SR Ca2+ uptake via SERCA pump (JSERCA), SR Ca2+ release flux via RyR (Jrel), SR Ca2+ leak (Jleak), and Na+ background current (INab). Four basic cellular models were considered with different baseline conditions and/or coupling to build the populations: normal myocyte coupled to normal myocyte (N-N), normal myocyte coupled to fibroblast (N-Fb), failing myocyte coupled to failing myocyte (HF-HF), and failing myocyte coupled to fibroblast (HF-Fb). Populations of 300 different human AP models were analyzed and compared to study the role of ion channel variation in alternans generation.

Alternans were measured in AP and CaT waveforms, as shown in Eqs (1) and (2), but the final model classification criterion was the Ca2+ alternans ratio [34], i.e. the relative differences between the small and large CaT amplitude (CaTA), since these were more evident than APD variations. While 0 indicated no alternation, models with Ca-alt > 0.1 were considered as having alternans.

(1)

APD alternans (APD-alt) were also calculated as the APD90 difference between two consecutive beats.

(2)

The models were classified into two groups according to whether or not the CaTs presented alternation. The analysis of the differences in parameter scaling factors between both groups determined the mechanistic contribution of ionic parameters to alternans generation. Statistical differences between groups were determined by the Wilcoxon rank-sum test (p < 0.05). Simulations and data analysis were performed on MATLAB (R2017b, Mathworks).

Tissue simulations

Electrical propagation was simulated in a set of two-dimensional strands (15 mm along the x-axis and 2 mm along the y-axis) formed by myocytes and fibroblasts. Longitudinal fiber orientation was assigned along the x-axis. The tissue was isotropic, i.e. longitudinal and transversal conductivities were identical. The domain was discretized at a spatial resolution of 0.1 mm, which provided 3000 squared elements and 3171 nodes within the tissue, where the cellular model was computed. The time-step for the numerical method was set to 0.01 ms to ensure numerical convergence.

20% of disperse fibroblasts were considered in half of the cases and were randomly intercalated among the myocytes, as shown in Fig 1A, i.e. 20% of the mesh nodes were solved with the fibroblast model. Electrotonic coupling between myocytes and fibroblasts was weaker than myocyte-myocyte coupling. A myocyte-myocyte diffusion coefficient (D) of 0.0024 cm2/ms was used, which yielded a conduction velocity (CV) of 0.7 m/s along the fibers in normal conditions [42]. To simulate intercellular uncoupling, diffusion between myocytes was reduced by 50% in HF (D = 0.0012 cm2/ms) and 3-fold for the fibroblast-myocyte and fibroblast-fibroblast coupling (D = 0.0008 cm2/ms). The assignment of the diffusion coefficient to each element was chosen as in Gomez et al. [20].

thumbnail
Fig 1. Electrophysiology in cardiac tissue simulations.

A) Distribution of fibroblasts (20%; black squares) in a tissue of myocytes. Each square represents a node of the mesh and measurement nodes are marked in red. B) Action potential (AP) and calcium transient (CaT) waveforms obtained in the central node of the tissue without fibroblasts by stimulating at 1.25 Hz (blue traces) and with 20% of fibroblasts (green traces).

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

Stimulation was applied to the left (short) side of the tissue and APs and CaTs were recorded in three different myocytes after 400 beats to ensure steady-state waveforms and alternans. We used the critical PCL obtained in the first part of the study (couple single-cells simulations) at which alternans appeared. Initial simulations at a basal PCL of 800 ms were conducted in tissues with different configurations to ensure that the results were independent of variations in fibroblast distribution (Fig 1B).

To reduce the computational cost of these tissue simulations, instead of considering a large population, as in the coupled single cells, we considered a subpopulation of models by varying the four most crucial factors that led to cellular alternans in the coupled single-cell simulations in the first part of the study, assuming that this subpopulation could be used to elucidate the onset of arrhythmogenic events and also the different behavior of failing and fibrotic tissue.

The most important parameters that contributed to alternans were selected from the coupled single-cell populations of models. The cardiac tissue subpopulation involved varying parameters by ±30% of their maximum value and accounted for all the possible combinations (34 = 81). As in the coupled single-cell populations, we classified models into two groups: alternans (Ca-alt > 0.1) and no alternans, and calculated APD-alt. The magnitude of alternans and the spatial distribution of alternans along the tissue were also taken into account.

The simulations were performed on ELVIRA software [43], which is based on the finite element method to solve the monodomain equation by the operator splitting technique [20,44].

Results

Coupled single-cell simulations

Restitution curve.

The initial simulations in coupled single cells provided a preliminary overview of the electrophysiological behavior in myocytes under pathological conditions and fast pacing rates to force beat-to-beat alternation. The following models were compared: N-N, N-Fb, HF-HF, and HF-Fb.

Fig 2 shows the restitution curves of APD90 and systolic Ca2+ of the four baseline models. These data indicate that an increase in heart rate leads to fluctuations in APD90 in consecutive beats in N-N, which are concomitant with Ca-alt. The large magnitude of Ca-alt is more evident than APD-alt. Fig 2 also shows the different alternans thresholds, according to the simulated conditions of the baseline models. While the CaT of N-N begins to alternate at a PCL of 300 ms, the electrical interactions between myocytes and fibroblasts in N-Fb increase this value to 360 ms. Myocyte-myocyte interactions seemed to slightly affect the bidirectional coupling between APD90 and Ca2+alternans, as there were small Ca-alt not translated into APD-alt.

thumbnail
Fig 2. Restitution curves of baseline coupled single-cell models.

A) Action potential duration (APD90) and B) Systolic [Ca2+]i. Comparison of electrophysiological indicators and alternans onset in myocyte-myocyte (solid circles) and myocyte-fibroblast (1:5 ratio) aggregations (empty squares) under normal (blue) and failing (red) conditions.The four traces are: normal myocyte coupled to normal myocyte (N-N), normal myocyte coupled to fibroblasts (N-Fb), failing myocyte coupled to failing myocyte (HF-HF), and failing myocyte coupled to fibroblasts (HF-Fb).

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

HF-HF myocytes, with changes in APD90 and CaT amplitude, do not show alternans at fast pacing rates and end up in loss of capture at an early PCL equal to 250 ms. The HF-Fb model presents a lower APD90 and loss of capture is delayed. Although it is difficult to compare the slope of the restitution curves for the four models at high frequencies, when the PCL is longer, HF-HF shows a steeper slope than N-N and fibroblast coupling contributes to reducing the slope.

To clarify the origin of alternans, an AP clamp test was simulated (results not shown). When forcing fixed APs, Ca2+ dynamics alternans were not eliminated, while when fixing [Ca2+]i, [Ca2+]JSR or Jrel, all the alternans types disappeared, revealing a Ca2+-driven mechanism. Positive coupling between Ca2+and APD alternans was also observed, as the large (small) CaT was always in phase with the long (short) APD. NCX was the mechanism coupling between Ca-alt and APD-alt, by generating a higher inward current during the large CaT to remove Ca2+, which prolonged APD and suppressed the small reducing effects of an inward Ca2+-inactivated ICaL.

Populations of models.

The four baseline models described above were used to build populations of cells by introducing natural variability into the electrophysiological variables. Small inter-individual differences can have remarkable effects on alternans development and may help to understand their origin. The boxplots in Fig 3 show the distribution of ionic conductances after classifying the models into the non-alternans or alternans groups. With scale factors in the range of 0.3–1.3, significant differences were found in INa, INaL, IKr, NCX, INaK, SERCA, Jrel and Jleak for the N-N population. In N-Fb models, the effect of ICaLwas also significant, as well as that of NCX, INaK, SERCA and Jrel. In the HF setting (HF-HF and HF-Fb), although there was agreement between the main mechanisms that contribute to alternans, unlike healthy myocytes, the failing myocytes in the alternans group presented the inverse modulation of Ca+2-related parameters, regardless of fibroblast content. While high ICaL and SERCA values and small NCX and Jleak values facilitated alternans in the failing populations, normal populations showed the opposite behaviour. As an exception, the impact of Jrel was consistent in all the populations because the higher Ca+2 release induced alternans.

thumbnail
Fig 3. Parameter distribution according to alternans occurrence in coupled single-cell populations.

The four populations are: A) normal myocyte coupled to normal myocyte (N-N), B) normal myocyte coupled to fibroblasts (N-Fb), C) failing myocyte coupled to failing myocyte (HF-HF), and D) failing myocyte coupled to fibroblasts (HF-Fb). Parameter values are represented by the scaling factors applied to the baseline models. A value of p < 0.05 (*) in the medians was considered statistically significant.

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

Cardiac tissue simulations

Two-dimensional populations of models.

Cardiac tissue myocyte simulations provide more reliable information than coupled single-cell simulations, due to electrical propagation. Four tissue configurations were considered: normal myocytes (N-tissue), failing myocytes (HF-tissue), normal myocytes and fibroblasts (N-Fb-tissue), and failing myocytes and fibroblasts (HF-Fb-tissue). Ionic parameters in myocytes (SERCA, ICaL, NCX and INaK) were varied to generate sets of tissue prone to alternans, on the basis of the mechanisms observed in the coupled single-cell populations. The variation of the four most important parameters related to alternans generation yielded 81 tissue models for each configuration (see Methods). The numbers of alternans cases in the different tissues are given in Table 2.

thumbnail
Table 2. Alternans occurrence in four subpopulations of cardiac tissue.

https://doi.org/10.1371/journal.pone.0217993.t002

Although the magnitude of alternans can be different, accounting for the central myocytes exhibiting Ca-alt > 0.1, the results obtained highlight the major incidence of alternans in HF-tissue (including loss of AP capture). The effect of disperse fibroblasts in cardiac tissue increases alternans events in N-Fb-tissue, while it reduces the occurrence in HF-Fb-tissue. In HF-tissue, many simulations end up in one AP for every two consecutive stimulations (loss of capture), which were considered as a special case of alternans, although the origin could be different. This 2:1 block activity was not observed in any of the other cases.

Fig 4 contains a graph showing Ca-alt and APD-alt values in circles and a gray scale of all the tissue sets arranged in 9x9 matrixes. The corresponding SERCA, ICaL, NCX, and INaK scaling factors specific to each myocyte variant classify the models in the matrix squares. In N-tissue, alternans mainly occur with reduced SERCA and ICaL and high NCX. While Ca-alt can rise to values of up to 0.8 in amplitude ratio, the maximum APD90 variation between consecutive beats is only 15 ms. The same mechanisms induce alternans in N-Fb-tissue, but most APD-alt become negligible despite the greater incidence of Ca-alt. When the tissue included failing myocytes, the ionic modulation inducing Ca-alt of similar magnitude is different to the one observed in healthy myocytes. In HF-tissue, the repolarization alternans values are higher (up to 30 ms), although they are still of considerable magnitude in HF-Fb-tissue. The analysis of the loss of capture of AP in HF-tissue (shown by a cross) reveals that increased NCX, together with low SERCA and low INaK activity lead to this outcome.

thumbnail
Fig 4. Ca2+ and APD alternans (Ca-alt and APD-alt) in tissues.

Ca-alt (upper panels) and APD-alt (lower panels) of every model variant is classified according to the parameter values that lead to those results. The four tissue configurations are: normal myocytes (N-tissue), failing myocytes (HF-tissue), normal myocytes and fibroblasts (N-Fb-tissue), and failing myocytes and fibroblasts (HF-Fb-tissue). Alternans magnitude is proportional to circle size and gray scale. Crosses represent 2:1 block.

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

Differences were found in CaT amplitude of the alternans and no alternans groups of the different populations. Since intracellular Ca2+ is linked to the amount of Ca2+ stored and released from the SR, the systolic Ca2+ peak was compared to maximum diastolic [Ca2+]JSR (see Fig 5). The linear correlation confirmed that cytosolic and sarcoplasmic Ca2+ content are related and also indicated that Ca-alt occurs within a specific range of [Ca2+]JSR values. This range (shaded area) was very similar for all four groups, but the normal myocytes (N-tissue and N-Fb-tissue) had the lowest values of [Ca2+], while HF (HF-tissue and HF-Fb-tissue) had the highest Ca2+ levels of the interval covered by the different models.

thumbnail
Fig 5. Relation between SR Ca2+ content, Ca2+ peak and alternans.

The CaT peak (systolic [Ca2+]i) is related to the the maximum junctional sarcoplasmic reticulum Ca2+content (diastolic [Ca2+]JSR) of the myocyte. Comparison of 4 populations obtained in cardiac tissue: A) normal myocytes (N-tissue), B) failing myocytes (HF-tissue), C) normal myocytes and fibroblasts (N-Fb-tissue), and D) failing myocytes and fibroblasts (HF-Fb-tissue). Shaded areas highlight regions where alternans cases prevail.

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

Spatial and mechanical discordance.

As cardiac alternans can become arrhythmogenic and the impact depends on the spatial distribution of the electrical fluctuations, we analyzed the evolution of alternans along the strand. Fig 6 shows APD90 and [Ca2+]i peak for two consecutive beats (blue for even and red for odd beats) in a representative case from each population, quantified in the central nodes of the tissue along the longitudinal axis (x). APD-alt patterns are long-short (L-S) and Ca-alt patterns are large-small (L-S). When a long APD is concomitant with a large CaT (same color) there is positive V-Ca coupling or electromechanical concordance, otherwise there is electromechanical discordance. Consecutive dots with the same color show that alternans between adjacent myocytes are in phase (spatially corcordant), while changes in color indicate a change in alternans patterns between consecutive myocytes (spatially discordant). Spatial concordance in APD alternans is predominant in the four populations because all the myocytes are in phase following the same pattern in APs. This can be seen in both N-tissue (Fig 6A) and HF-tissue (Fig 6C), where spatial concordance is found in APD and Ca alternans.The CV in N-tissue was 0.65 m/s and the characteristic cellular uncoupling in HF reduced it to 0.4 m/s in HF-tissue but the alternans electrical propagation does not seem to be affected. In some cases, discordance can appear in Ca2+ dynamics in normal tissue or in both APD and Ca2+ activity in failing tissue. Of all the N-tissue simulated, only 13% of tissue variants with Ca-alt are spatially discordant, while in N-Fb-tissue all the Ca-alt become discordant. In these tissues with fibroblasts, APD90 alternans are concordant, so that electromechanical discordance occurs (Fig 6B). However, in failing tissues, discordant Ca-alt do not increase significantly from HF-tissue (7.5% alternans) to HF-Fb-tissue (12.5% alternans). What distinguishes HF-Fb-tissue from N-Fb-tissue is that spatial discordance occurs in both Ca-alt and APD-alt and is less chaotic (Fig 6D), which could be related to a slower CV (0.34 m/s vs 0.44 m/s in normal tissue with fibroblasts).

thumbnail
Fig 6. Spatial distribution of Ca2+ peak and APD90 in tissues.

Node values (myocytes) are given for the central fiber of the strand for 2 consecutive beats (red and blue). The four tissue configurations are: A) normal myocytes (N-tissue), B) normal myocytes with fibroblasts (N-Fb-tissue), C) failing myocytes (HF-tissue), and D) failing myocytes with fibroblasts (HF-Fb-tissue). Y-axis has the same scale in Ca2+ and APD graphs, respectively but values can vary for the sake of clarity.

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

In addition to spatial differences, Fig 6 shows maximum and minimum APD90 and systolic [Ca2+]i values. Although these values also depend on the modulation of ionic parameters, the figures show that although there are differences in APD90, small and large CaTs from different myocytes with alternating activity can become similar in magnitude. The APD90 prolongation characteristic of HF appears in HF-tissue and affects APD-alt, but there is also a significant reduction in APD90 in fibrotic tissues (N-Fb-tissue and HF-Fb-tissue) due to myocyte-fibroblast interactions.

Since fibroblasts seem to contribute to a more chaotic spatial distribution of Ca-alt, the whole tissue is represented in Fig 7. Black squares indicate the specific location of disperse fibroblasts and the color scale represents alternans magnitude. Alternans of opposite signs, i.e. L-S and S-L, in both CaT (upper panels) and APD (lower panels) are shown in red and blue. As in Figs 6B and 7A shows that in N-Fb-tissue the spatially discordant Ca-alt affect the entire tissue but are not translated to APD-alt. Discordant areas vary between tissue variants of the N-Fb-tissue population, but all show a similar chaotic Ca-alt distribution (results not shown) with no patterns around the fibroblast locations. When spatial Ca-alt occur in HF-Fb-tissue (Fig 7B), they are confined to certain tissue regions in which APD-alt also appear, enhancing electromechanical concordance. Fig 7B also shows that there is a gradual transition between positive and negative APD-alt where there is no alternation (white areas in upper panel separate blue and red areas), although Ca-alt change abruptly from one myocyte to another (blue nodes next to red nodes).

thumbnail
Fig 7. Spatial concordance of alternans (Ca-alt and APD-alt) along the tissue.

A) Tissue with normal myocytes and fibroblasts (N-Fb-tissue), and B) tissue with failing myocytes and fibroblasts (HF-Fb-tissue). Opposite signs (red vs blue) indicate out-of-phase alternans. Fibroblasts are represented in black.

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

Electrically discordant alternans are the most arrhythmogenic type, but concordant alternans often precede discordant alternans and play an important role. Ca-alt help in the identification of APD-alt, as they are concomitant, but also indicate dysfuntional contraction behavior. Discordant Ca-alt could aggravate the situation by increasing the number of myocytes that do not contract synchronously.

Discussion

Main findings

Mechanical and electrical alternans, which are known to be bidirectionally coupled, may compromise cardiac function. This study shows that alternans have their origin in Ca2+ alterations and reveals how Ca2+ impairment in HF or fibroblast interaction influences the development of cellular alternans and how they appear in cardiac tissue. Our main findings are as follows: (i) Ca2+ mechanisms, involving SERCA, ICaL, and NCX, are responsible for the occurrence of cellular alternans by modulating Ca2+ cycling to threshold levels; these mechanisms act differently in normal and failing mocytes, (ii) the degree of Ca2+ impairment, given by [Ca2+]JSR, determines the myocyte susceptibility to Ca-alt development, (iii) fibroblasts induce discordant Ca-alt in normal cardiac tissues, (iv) the high number of alternans in failing tissue with myocyte-fibroblast interaction is reduced due to exacerbated Ca2+ impairment, and (v) discordant Ca-alt are more prone to be translated into discordant APD-alt in failing tissue with fibroblasts than in normal tissue with fibroblasts. Fibroblasts and HF can therefore promote alternans because they alter Ca2+ dynamics, deteriorate contractility and increase the risk of cardiac arrhythmia. Restoring Ca2+ dynamics could thus potentially be used to prevent alternans and improve cellular contraction and cardiac function, although the differences between normal and failing tissues underline the need to consider HF when developing new pharmacological treatments.

Origin of cellular alternans

Although cardiac alternans usually occur at fast pacing rates, the threshold for the onset of fluctuations can vary between individuals. It is also a characteristic outcome in heart disease. The first cellular theory for cardiac alternans was based on the relationship between APD alternans and the APD restitution slope [26]. Alternation was generated and sustained when the slope was higher than 1, so steeper restitution curves led to alternans at a longer PCL. However, we found that myocyte-fibroblast coupling, which shortens APD, promotes alternans at a longer PCL, despite a flatter restitution curve. These results indicate that there are other factors that contribute to generating alternans and that the analysis of the restitution curve has some limitations. If alternans are voltage-driven, as Xie et al. [45] found in the Phase I of the Luo and Rudy model [46], fibroblasts that shorten APD can induce alternans at a lower PCL. But when the electrophysiological model incorporates detailed Ca2+ cycling (as in our case) myocytes-fibroblast coupling promotes alternans because the alternans onset is then driven by Ca2+ cycling instabilities [45]. In fact, Ca2+ dynamics, which also oscillates from beat to beat, is one of the most important factors in alternans research. The bidirectional coupling between membrane voltage and intracellular Ca2+ explains why APD-alt and Ca-alt occur concomitantly [47], as can be seen in the restitution curves obtained for APD90 and Ca2+. However, there is still controversy about the origin and underlying causes of these fluctuations, due to the Ca2+ dependence of ion channels, which also control Ca2+ cycling [23,4749].

In the present computational study, in which the ORd model was used [36], the alternans found in single-cell pairs and tissue are similar to the alternans in the initial isolated myocytes. The ORd model was validated to reproduce human APD alternans (around 10 ms) with beat-to-beat alternans in the Ca2+ subsystem at high frequencies (PCL< 300 ms). They also occur in a 100-endo-cell strand, and pacing rates faster than 230 ms cause 2:1 block [36], which agrees with our results in tissue. This type of conduction block can also be observed in cellular simulations, usually at fast frequencies, under conditions that prolong APs. The difference is that a second short AP is developed after an AP longer than the cycle length. In previous studies focusing on alternans mechanisms, these cases were considered as alternans because of their long-short pattern [50]. In tissue simulations, we found that the short APs developed in the externally stimulated myocytes did not propagate through the tissue because the myocytes were in a refractory period, leading to 2:1 block. Fig 2A shows that the long APD in failing myocytes causes loss of capture at higher PCLs and this explains why, at a fixed PCL of 300 ms, 2:1 block cases were only observed in HF tissue without fibroblasts. Despite the different mechanisms of alternans and the conduction block, both were included in the study due to their potential arrhythmogenic consequences.

Our first simulations at different PCLs determined that Ca-alt were easier to identify than APD-alt and so this was used as the criterion to detect cellular alternans. In previous studies, AP amplitude alternans have also been used to explain T-wave alternans instead of APD variations and were able to predict malignant arrhythmias in patients with systolic dysfunction [51]. However, the link between mechanical and electrical alternans has also proved that mechanical fluctuations are more visible than T-wave alternans, so that research on these factors could be useful for managing HF patients [21].

Alternans mechanisms

Introducing natural electrophysiological variability into ionic channel conductances and transporters provides a population of coupled cells with different electrophysiological characteristics. The analysis of the model variants prone to alternans revealed the electrophysiological mechanisms that promote these fluctuations, i.e. the Ca2+-handling proteins (SERCA, ICaL, NCX, and Jrel). This agrees with the findings of Zhou et al. [31] and Walmsley et al. [50] for non-failing single myocytes. Walmsley et al. [50] associated the appearance of alternans to APD prolongation in failing myocytes when IKr was low, but not to other parameters, since the alternans were measured as differences in AP and CaT duration. As a novelty, we introduced fibroblast coupling to normal and failing myocytes and then compared them to myocyte-myocyte coupling. The results indicated that although the mechanisms driving alternans in the different populations were similar they had opposite effects on normal conditions and HF.

CaMK appears to play an important role in alternans generation as it controls Ca2+ cycling [36,52]. This could partially explain the fact that HF enhances cellular alternans, since a 50% increase in the CaMK active fraction is a characteristic alteration of failing myocytes (Table 1). Ca2+-handling proteins help to maintain normal Ca2+ cycling in myocytes to ensure contractility and their disruption can lead to impaired Ca2+ dynamics in the different cell compartments. This has been linked to beat-to-beat Ca2+ oscillations [53], including fluctuations in the SR Ca2+ content [30], which are the root cause of cardiac alternans. Recent studies have reported that the SR imbalance between Ca2+ reuptake and Ca2+ release explains the genesis of Ca-alt [30,31]. We found that the SERCA pump, which uptakes Ca2+ into the sarcoplasmic reticulum (SR), is involved in alternans-prone myocytes, while SERCA overexpression may possibly suppress cardiac alternans, as has been shown by Cutler et al. [33]. However, unlike previous studies, in HF we found that alternans were not eliminated by increasing Ca2+ uptake into the SR [5456]. Our results indicate that in HF, when CaT is strongly depressed alternans are not generated. However, if SERCA activity increases, Ca2+ cycling improves although it still remains at failure levels and can lead to alternans. These levels are mainly determined by junctional SR Ca2+ maximal content. Normal myocytes can also reach these failing Ca2+ levels when mechanisms like ICaL, SERCA, NCX, and INaK are altered, after which Ca-alt may occur. With reduced ICaL, Ca2+-induced Ca2+ release (CICR) tends to fail, high NCX contributes to intracellular Ca2+ imbalance by extruding Ca2+ out of the cell, and the higher INaK could indirectly facilitate inward NCX.

Many experimental studies have also highlighted the role of RyR in alternans generation due to RyR refractoriness and a failure in the CICR process [28,29,57]. Computer simulations with a spatially distributed Ca2+ cycling model integrated in the AP model could explain how Ca2+ sparks can lead to cellular alternans due to the randomness, refractoriness and recruitment properties of the couplon network [58]. However, in a less detailed Ca2+ cycling model, such as the one formulated in the ORd model, RyRs cannot show refractoriness, which could explain why we found that, unlike other parameters, changes in RyR hardly affected CaT. However, Tomek et al. [59] showed that slow dynamics and depressed SR Ca2+ release through RyR can also affect vulnerability to alternans. We therefore believe that incomplete JSR recovery is the factor that affects Ca2+ release and generates alternans [52].

Alternans spatial distribution

Regarding the spatiotemporal patterns in cardiac tissue, we found that alternans entail different degrees of risk that depend on individual conditions. Electrical propagation runs smoothly in normal tissue, which explains the tissue’s spatial and electromechanical concordance. However, we found that cellular uncoupling in HF may promote the formation of alternans, in agreement with Hammer et al. [60], who found that myocytes with reduced intercellular coupling had higher alternans occurrence and amplitude, with a concordant pattern between neighboring cells. Fibroblasts also reduce intercellular coupling but can lead to spatially discordant alternans in normal tissue with fibroblasts, because, apart from reducing the CV they alter the myocytes’ electrical activity. Majumder et al. [61] observed discordant APD-alt preceding reentries in cultured monolayers of neonatal rat ventricular myocytes with Cx43 inactivation and myofibroblasts. Both slow CV and interstitial fibrosis were critical to the formation of such arrhythmogenic alternans at fast pacing. This agrees with our simulated HF and HF with fibroblast results, in which CV is depressed and discordant APD-alt are prone to occur. Despite the cellular uncoupling with fibroblasts we did not detect APD discordance in normal fibrotic fibers, but Ca-alt became spatially asynchronous due to the slow diffusion of Ca2+ ions [62] in comparison to the fast diffusion of membrane potential. As a consequence, cardiac tissue can present adjacent myocytes with both in-phase and out-of-phase V-Ca2+ coupling, which reduces APD-alt [63]. Electromechanically discordant alternans were also found by Xie et al. [45] when passive fibroblasts were inserted into a 2D tissue of rabbit ventricular myocytes. According to MacCannell et al. [38], the use of a passive fibroblast model hardly modulates myocyte AP or CaT, whereas the electrophysiological model does have a strong effect. As our objective was to analyze the most critical conditions only active fibroblasts were studied.

Electrophysiological remodeling of Ca2+-handling proteins is the main cause of Ca-alt in HF, which is then translated to APD-alt. Although we found spatially concordant alternans in HF tissue, they could become discordant at lower PCL, as was shown in [64]. There is also experimental evidence to show that HF enhances susceptibility to arrhythmogenic cardiac alternans, since ventricular fibrillation is only inducible after spatial discordance appears, which is more common in failing wedge preparations and occurs at a slower heart rate [54]. It is difficult to say whether or not these tissues included fibroblasts, but as optical mapping can be hampered by fibrotic areas, many experimental measurements avoid tissues with a high proportion of fibroblasts. In any case, of the computational alternans studies, this is the first to consider myocyte-fibroblast coupling in HF tissue. Apparently, fibroblasts do not increase discordance in failing tissues, since we did not find any differences with tissues without fibroblasts. Compared to normal tissue, HF involves more Ca2+ impairment and reduced CV, which may contribute to the positive APD-alt and Ca-alt coupling. Failing myocytes become susceptible to Ca-alt with unstable Ca2+ dynamics provided that CaT is not very depressed, and the resultant repolarization alternans can be high. If the tissue properties were to trigger spatial discordance, these alternans could become arrhythmogenic.

Arrhythmogenic activity in cardiac tissue can have a non-electric origin if Ca-alt is the initial cause. In contrast to normal hearts, the attempt to restore Ca2+ dynamics in HF could not improve cardiac function but did contribute to Ca-alt development and its associated effects. The conclusion is therefore that mechanisms that restore Ca2+ instabilities produced by fibroblasts and/or HF help to recover normal myocyte electrophysiological activity and avoid heart problems only if they are applied to specific heart conditions.

Limitations

Mathematical models can help in the mechanistic analysis of alternans onset and complement experimental findings, although computational results may be limited by a number of uncertainties. Only one fibroblast model and one configuration for myocyte-fibroblast coupling were considered. The electrophysiological fibroblast model is thought to be more accurate, as it represents the ion currents in the fibroblasts, but using a different model (eg. passive fibroblast) or variations in the parameters could change the results. In our simulations in particular the fibroblast model does not take into account the potential Ca2+ cycling pathways in non-myocytes because of the lack of data related to cardiac fibroblasts, and this could contribute to Ca2+ dynamics in myocytes and therefore to developing alternans. Myofibroblasts, the contractile fibroblast activated form, arise in response to injury and present a different phenotype from fibroblasts. However, as no specific electrophysiological models of these contractile cells have been developed to date, in the present study we used the same mathematical fibroblast model in both normal conditions and HF, in spite of the fact that fibroblast electrophysiology may be altered in response to HF. We implemented the most critical of the different configurations for myocyte-fibroblast coupling, in which fibroblasts form conduction pathways and act as an electrotonic load. Representing fibroblasts as the same size as myocytes in a continuous mesh still has certain limitations.

Another limitation was the scarcity of human electrophysiological data, which would have been very useful to validate our findings. Our results were thus compared to animal models, including some species with a different electrophysiology to that of humans.

Despite its limitations, this study provides information that plausibly explains the response of myocytes to HF and fibroblast coupling in humans and can thus be used as a guide in the design of HF treatments.

References

  1. 1. Glukhov A V., Fedorov V V., Kalish PW, Ravikumar VK, Lou Q, Janks D, et al. Conduction remodeling in human end-stage nonischemic left ventricular cardiomyopathy. Circulation. 2012;125: 1835–1847. pmid:22412072
  2. 2. Lou Q, Fedorov V V., Glukhov A V., Moazami N, Fast VG, Efimov IR. Transmural heterogeneity and remodeling of ventricular excitation- contraction coupling in human heart failure. Circulation. 2011;123: 1881–1890. pmid:21502574
  3. 3. Gomez JF, Cardona K, Trenor B. Lessons learned from multi-scale modeling of the failing heart. J Mol Cell Cardiol. 2015;89: 146–159. pmid:26476237
  4. 4. Kohl P, Gourdie RG. Fibroblast-myocyte electrotonic coupling: does it occur in native cardiac tissue? J Mol Cell Cardiol. 2014;70: 37–46. pmid:24412581
  5. 5. Gaudesius G, Miragoli M, Thomas SP, Rohr S. Coupling of cardiac electrical activity over extended distances by fibroblasts of cardiac origin. Circ Res. 2003;93: 421–428. pmid:12893743
  6. 6. Kohl P, Camelliti P, Burton FL, Smith GL. Electrical coupling of fibroblasts and myocytes: relevance for cardiac propagation. J Electrocardiol. 2005;38: 45–50. pmid:16226073
  7. 7. Camelliti P, Green CR, LeGrice I, Kohl P. Fibroblast Network in Rabbit Sinoatrial Node. Circ Res. 2004;94: 828–835. pmid:14976125
  8. 8. Rook MB, van Ginneken AC, de Jonge B, el Aoumari A, Gros D, Jongsma HJ. Differences in gap junction channels between cardiac myocytes, fibroblasts, and heterologous pairs. Am J Physiol Physiol. 1992;263: C959–C977. pmid:1279981
  9. 9. Mahoney VM, Mezzano V, Mirams GR, Maass K, Li Z, Cerrone M, et al. Connexin43 contributes to electrotonic conduction across scar tissue in the intact heart. Sci Rep. 2016;6: 26744. pmid:27244564
  10. 10. Quinn TA, Camelliti P, Rog-Zielinska EA, Siedlecka U, Poggioli T, O’Toole ET, et al. Electrotonic coupling of excitable and nonexcitable cells in the heart revealed by optogenetics. Proc Natl Acad Sci. 2016;113: 14852–14857. pmid:27930302
  11. 11. Rubart M, Tao W, Lu X-L, Conway SJ, Reuter SP, Lin S-F, et al. Electrical coupling between ventricular myocytes and myofibroblasts in the infarcted mouse heart. Cardiovasc Res. 2018;114: 389–400. pmid:29016731
  12. 12. Miragoli M, Gaudesius G, Rohr S. Electrotonic modulation of cardiac impulse conduction by myofibroblasts. Circ Res. 2006;98: 801–810. pmid:16484613
  13. 13. Jacquemet V, Henriquez CS. Loading effect of fibroblast-myocyte coupling on resting potential, impulse propagation, and repolarization: insights from a microstructure model. Am J Physiol Circ Physiol. 2008;294: H2040–H2052. pmid:18310514
  14. 14. Li Y, Asfour H, Bursac N. Age-dependent functional crosstalk between cardiac fibroblasts and cardiomyocytes in a 3D engineered cardiac tissue. Acta Biomater. 2017;55: 120–130. pmid:28455218
  15. 15. Maleckar MM, Greenstein JL, Giles WR, Trayanova NA. Electrotonic coupling between human atrial myocytes and fibroblasts alters myocyte excitability and repolarization. Biophys J. Biophysical Society; 2009;97: 2179–2190. pmid:19843450
  16. 16. Zlochiver S, Muñoz V, Vikstrom KL, Taffet SM, Berenfeld O, Jalife J. Electrotonic myofibroblast-to-myocyte coupling increases propensity to reentrant arrhythmias in two-dimensional cardiac monolayers. Biophys J. 2008;95: 4469–80. pmid:18658226
  17. 17. Nguyen TP, Xie Y, Garfinkel A, Qu Z, Weiss JN. Arrhythmogenic consequences of myofibroblast–myocyte coupling. Cardiovasc Res. 2012;93: 242–251.
  18. 18. Greisas A, Zlochiver S. The Multi-Domain Fibroblast/Myocyte Coupling in the Cardiac Tissue: A Theoretical Study. Cardiovasc Eng Technol. 2016;7: 290–304. pmid:27150222
  19. 19. Sridhar S, Vandersickel N, Panfilov A V. Effect of myocyte-fibroblast coupling on the onset of pathological dynamics in a model of ventricular tissue. Sci Rep. 2017;7: 40985. pmid:28106124
  20. 20. Gomez JF, Cardona K, Martinez L, Saiz J, Trenor B. Electrophysiological and Structural Remodeling in Heart Failure Modulate Arrhythmogenesis. 2D Simulation Study. PLoS One. 2014;9: e103273. pmid:25054335
  21. 21. Kodama M, Kato K, Hirono S, Okura Y, Hanawa H, Yoshida T, et al. Linkage Between Mechanical and Electrical Alternans in Patients with Chronic Heart Failure. J Cardiovasc Electrophysiol. 2004;15: 295–299. pmid:15030419
  22. 22. Rosenbaum DS, Jackson LE, Smith JM, Garan H, Ruskin JN, Cohen RJ. Electrical Alternans and Vulnerability to Ventricular Arrhythmias. N Engl J Med. 1994;330: 235–241. pmid:8272084
  23. 23. Jordan PN, Christini DJ. Action Potential Morphology Influences Intracellular Calcium Handling Stability and the Occurrence of Alternans. Biophys J. 2006;90: 672–680. pmid:16239324
  24. 24. Cherry EM. Distinguishing mechanisms for alternans in cardiac cells using constant-diastolic-interval pacing. Chaos An Interdiscip J Nonlinear Sci. 2017;27: 093902. pmid:28964159
  25. 25. Groenendaal W, Ortega FA, Krogh-Madsen T, Christini DJ. Voltage and calcium dynamics both underlie cellular alternans in cardiac myocytes. Biophys J. 2014;106: 2222–32. pmid:24853751
  26. 26. Nolasco JB, Dahlen RW. A graphic method for the study of alternation in cardiac action potentials. J Appl Physiol. 1968;25: 191–196. pmid:5666097
  27. 27. Qu Z, Liu MB, Nivala M. A unified theory of calcium alternans in ventricular myocytes. Sci Rep. 2016;6: 1–14.
  28. 28. Díaz ME, Eisner DA, O’Neill SC. Depressed ryanodine receptor activity increases variability and duration of the systolic Ca2+ transient in rat ventricular myocytes. Circ Res. 2002;91: 585–593. pmid:12364386
  29. 29. Picht E, DeSantiago J, Blatter LA, Bers DM. Cardiac alternans do not rely on diastolic sarcoplasmic reticulum calcium content fluctuations. Circ Res. 2006;99: 740–748. pmid:16946134
  30. 30. Díaz ME, O’Neill SC, Eisner DA. Sarcoplasmic Reticulum Calcium Content Fluctuation Is the Key to Cardiac Alternans. Circ Res. 2004;94: 650–656. pmid:14752033
  31. 31. Zhou X, Bueno-Orovio A, Orini M, Hanson B, Hayward M, Taggart P, et al. In vivo and in silico investigation into mechanisms of frequency dependence of repolarization alternans in human ventricular cardiomyocytes. Circ Res. 2016;118: 266–278. pmid:26602864
  32. 32. Xie LH, Sato D, Garfinkel A, Qu Z, Weiss JN. Intracellular Ca alternans: Coordinated regulation by sarcoplasmic reticulum release, uptake, and leak. Biophys J. 2008;95: 3100–3110. pmid:18539635
  33. 33. Cutler MJ, Wan X, Laurita KR, Hajjar RJ, Rosenbaum DS. Targeted SERCA2a Gene Expression Identifies Molecular Mechanism and Therapeutic Target for Arrhythmogenic Cardiac Alternans. Circ Arrhythmia Electrophysiol. 2009;2: 686–694. pmid:19948504
  34. 34. Kanaporis G, Blatter LA. The Mechanisms of Calcium Cycling and Action Potential Dynamics in Cardiac Alternans. Circ Res. 2015;116: 846–856. pmid:25532796
  35. 35. Pastore JM, Girouard SD, Laurita KR, Akar FG, Rosenbaum DS. Mechanism Linking T-Wave Alternans to the Genesis of Cardiac Fibrillation. Circulation. 1999;99: 1385–1394. pmid:10077525
  36. 36. O’Hara T, Virág L, Varró A, Rudy Y. Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS Comput Biol. 2011;7: e1002061. pmid:21637795
  37. 37. Mora MT, Ferrero JM, Romero L, Trenor B. Sensitivity analysis revealing the effect of modulating ionic mechanisms on calcium dynamics in simulated human heart failure. PLoS One. 2017;12: e0187739. pmid:29117223
  38. 38. MacCannell KA, Bazzazi H, Chilton L, Shibukawa Y, Clark RB, Giles WR. A mathematical model of electrotonic interactions between ventricular myocytes and fibroblasts. Biophys J. 2007;92: 4121–4132. pmid:17307821
  39. 39. Spach MS, Heidlage JF, Dolber PC, Barr RC. Electrophysiological effects of remodelling cardiac gap junctions and cell size: experimental and model studies of normal cardiac growth. Circ Res. 2000;86: 302–311. pmid:10679482
  40. 40. Kieval RS, Spear JF, Moore EN. Gap junctional conductance in ventricular myocyte pairs isolated from postischemic rabbit myocardium. Circ Res. 1992;71: 127–136. pmid:1606660
  41. 41. Gomez JF, Cardona K, Romero L, Ferrero JM, Trenor B. Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 1D simulation study. PLoS One. 2014;9: e106602. pmid:25191998
  42. 42. Taggart P, Sutton PM, Opthof T, Coronel R, Trimlett R, Pugsley W, et al. Inhomogeneous transmural conduction during early ischaemia in patients with coronary artery disease. J Mol Cell Cardiol. 2000;32: 621–630. pmid:10756118
  43. 43. Heidenreich E. Algoritmos para ecuaciones de reacción difusión aplicados a electrofisiología. Ph.D. Thesis. Universidad de Zaragoza. 2009. https://institutoi4.net/wp-content/uploads/2017/08/TesisEAH.pdf
  44. 44. Heidenreich EA, Ferrero JM, Doblaré M, Rodríguez JF. Adaptive Macro Finite Elements for the Numerical Solution of Monodomain Equations in Cardiac Electrophysiology. Ann Biomed Eng. 2010;38: 2331–2345. pmid:20238165
  45. 45. Xie Y, Garfinkel A, Weiss JN, Qu Z. Cardiac alternans induced by fibroblast-myocyte coupling: mechanistic insights from computational models. Am J Physiol Circ Physiol. 2009;297: H775–H784. pmid:19482965
  46. 46. Luo CH, Rudy Y. A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction. Circ Res. 1991;68: 1501–26. pmid:1709839
  47. 47. Pruvot EJ, Katra RP, Rosenbaum DS, Laurita KR. Role of calcium cycling versus restitution in the mechanism of repolarization alternans. Circ Res. 2004;94: 1083–90. pmid:15016735
  48. 48. Kanaporis G, Blatter LA. Membrane potential determines calcium alternans through modulation of SR Ca 2+ load and L-type Ca 2+ current. J Mol Cell Cardiol. 2017;105: 49–58. pmid:28257761
  49. 49. Goldhaber JI, Xie L-H, Duong T, Motter C, Khuu K, Weiss JN. Action Potential Duration Restitution and Alternans in Rabbit Ventricular Myocytes. Circ Res. 2005;96: 459–466. pmid:15662034
  50. 50. Walmsley J, Rodriguez JF, Mirams GR, Burrage K, Efimov IR, Rodriguez B. mRNA expression levels in failing human hearts predict cellular electrophysiological remodeling: A population-based simulation study. PLoS One. 2013;8: e56359. pmid:23437117
  51. 51. Narayan SM, Bayer JD, Lalani G, Trayanova NA. Action Potential Dynamics Explain Arrhythmic Vulnerability in Human Heart Failure. J Am Coll Cardiol. 2008;52: 1782–1792. pmid:19022157
  52. 52. Livshitz LM, Rudy Y. Regulation of Ca2+ and electrical alternans in cardiac myocytes: role of CAMKII and repolarizing currents. Am J Physiol Circ Physiol. 2007;292: H2854–H2866. pmid:17277017
  53. 53. Wilson LD, Wan X, Rosenbaum DS. Cellular alternans: A mechanism linking calcium cycling proteins to cardiac arrhythmogenesis. Ann N Y Acad Sci. 2006;1080: 216–234. pmid:17132786
  54. 54. Wilson LD, Jeyaraj D, Wan X, Hoeker GS, Said TH, Gittinger M, et al. Heart failure enhances susceptibility to arrhythmogenic cardiac alternans. Hear Rhythm. 2009;6: 251–259. pmid:19187920
  55. 55. Cutler MJ, Wan X, Plummer BN, Liu H, Deschenes I, Laurita KR, et al. Targeted sarcoplasmic reticulum Ca2+ ATPase 2a gene delivery to restore electrical stability in the failing heart. Circulation. 2012;126: 2095–2104. pmid:23019291
  56. 56. Bayer JD, Narayan SM, Lalani GG, Trayanova NA. Rate-dependent action potential alternans in human heart failure implicates abnormal intracellular calcium handling. Hear Rhythm. 2010;7: 1093–1101. pmid:20382266
  57. 57. Wang L, Myles RC, De Jesus NM, Ohlendorf AKP, Bers DM, Ripplinger CM. Optical mapping of sarcoplasmic reticulum Ca2+in the intact heart: Ryanodine receptor refractoriness during alternans and fibrillation. Circ Res. 2014;114: 1410–1421. pmid:24568740
  58. 58. Rovetti R, Cui X, Garfinkel A, Weiss JN, Qu Z. Spark-Induced Sparks As a Mechanism of Intracellular Calcium Alternans in Cardiac Myocytes. Circ Res. 2010;106: 1582–1591. pmid:20378857
  59. 59. Tomek J, Tomková M, Zhou X, Bub G, Rodriguez B. Modulation of Cardiac Alternans by Altered Sarcoplasmic Reticulum Calcium Release: A Simulation Study. Front Physiol. 2018;9: 1306. pmid:30283355
  60. 60. Hammer KP, Ljubojevic S, Ripplinger CM, Pieske BM, Bers DM. Cardiac myocyte alternans in intact heart: Influence of cell-cell coupling and β-adrenergic stimulation. J Mol Cell Cardiol. 2015;84: 1–9.
  61. 61. Majumder R, Engels MC, de Vries AAF, Panfilov A V., Pijnappels DA. Islands of spatially discordant APD alternans underlie arrhythmogenesis by promoting electrotonic dyssynchrony in models of fibrotic rat ventricular myocardium. Sci Rep. 2016;6: 24334. pmid:27072041
  62. 62. Shiferaw Y, Karma A. Turing instability mediated by voltage and calcium diffusion in paced cardiac cells. Proc Natl Acad Sci. 2006;103: 5670–5675. pmid:16574775
  63. 63. Sato D, Shiferaw Y, Garfinkel A, Weiss JN, Qu Z, Karma A. Spatially discordant alternans in cardiac tissue: Role of calcium cycling. Circ Res. 2006;99: 520–527. pmid:16902177
  64. 64. Ponnaluri AVS, Perotti LE, Liu M, Qu Z, Weiss JN, Ennis DB, et al. Electrophysiology of Heart Failure Using a Rabbit Model: From the Failing Myocyte to Ventricular Fibrillation. PLoS Comput Biol. 2016;12: 1–25. pmid:27336310