Abstract

The ratio of densities of Na-Ca exchanger current () in the t-tubular and surface membranes (-ratio) computed from the values of and membrane capacitances () measured in adult rat ventricular cardiomyocytes before and after detubulation ranges between 1.7 and 25 (potentially even 40). Variations of action potential waveform and of calcium turnover within this span of the -ratio were simulated employing previously developed model of rat ventricular cell incorporating separate description of ion transport systems in the t-tubular and surface membranes. The increase of -ratio from 1.7 to 25 caused a prolongation of APD (duration of action potential at 90% repolarisation) by 12, 9, and 6% and an increase of peak intracellular Ca2+ transient by 45, 19, and 6% at 0.1, 1, and 5 Hz, respectively. The prolonged APD resulted from the increase of due to the exposure of a larger fraction of Na-Ca exchangers to higher Ca2+ transients under the t-tubular membrane. The accompanying rise of Ca2+ transient was a consequence of a higher Ca2+ load in sarcoplasmic reticulum induced by the increased Ca2+ cycling between the surface and t-tubular membranes. However, the reason for large differences in the -ratio assessed from measurements in adult rat cardiomyocytes remains to be explained.

1. Introduction

In cardiac ventricular myocytes, the Na-Ca exchanger forms a major pathway enabling Ca2+ extrusion from the cells [13] and takes part in the process of excitation-contraction (EC) coupling [46]. Hence, the exchanger plays an important role in the regulation of cellular Ca2+ content, Ca2+ transient, and thereby contractility; see [79] for review. The majority of immunolocalization studies found the Na-Ca exchanger to be predominantly localized in the t-tubular membrane of adult ventricular cardiomyocytes [1013]. While these studies reported the presence of Na-Ca exchanger proteins, they were not able to evaluate their activity in the t-tubular and surface membrane pools.

The recent method based on detubulation of cardiac cells [14] made it possible to assess the proportion of ionic current components and the proportion of membrane capacitances located at the t-tubular and surface membranes. Since the Na-Ca exchanger is electrogenic, this method allowed evaluating the membrane distribution of its activity. Application of this method to adult rat cardiomyocytes showed that the density of Na-Ca exchanger current () is higher in the t-tubular membrane [1518]. However, the ratio of Na-Ca current density in t-tubular to surface membrane (-ratio for short) assessed from data reported in these studies ranges between 1.7 (derived from [16]) and 25 (the lower estimate based on the data measured by Gadeberg et al. [15]). The aim of this paper is to explore the physiological consequences of such different values of the -ratio and to discuss the possible factors that may contribute to the dispersion of these estimations.

2. Methods

2.1. Assessment of -Ratio

The assessment of the ratio of Na-Ca current density in t-tubular membrane to that in surface membrane from experimentally measured values of and membrane capacitances () before and after detubulation was done using the following set of relations: where , , , (measured quantities) stand for the average densities of the Na-Ca exchange current (pA/pF) and the membrane capacitances (pF) in intact and detubulated cells and , , , (unknown quantities) stand for the densities of the Na-Ca exchange current and membrane capacitances at the surface and t-tubular membranes. The symbol denotes the fraction of t-tubules that resisted detubulation as estimated from confocal images of intact and detubulated myocytes.

Solution of the set of equations (1) for the unknown quantities , , , and yields:Hence, the -ratio can be expressed as

2.2. Model Evaluation of the Impact of -Ratio on Cellular Activity

To evaluate the impact of the -ratio on cellular activity, we used our model of rat ventricular cardiomyocyte [19] incorporating separate description of ion transport system in the t-tubular and surface membranes, separate dyadic, and subsarcolemmal spaces adjacent to the t-tubular and surface membranes and separate junctional sarcoplasmic reticulum (SR) compartments adjacent to the t-tubular and surface dyadic spaces (for the scheme of the model, see Figure  1 in [19]). The distribution of ion transporters between the t-tubular and surface membranes in the model is determined by the fractions of particular transporters in the t-tubular membrane (Table 1). The t-tubular fraction of Na-Ca exchangers is specified by the parameter . To explore the responses of the model to different values of -ratio, as resulted from (3) after insertion of experimental data, we set the values of according to the following relation:The simulations were performed using the computational software MATLAB v.7.2 (MathWorks, Natick, MA, USA). To attain a dynamic steady state at each stimulation frequency (0.1, 1, and 5 Hz), the model was paced for 600 s of equivalent cell lifetime by 1 ms pulses (the magnitude was set at twice the respective threshold value) under all specified conditions. The MATLAB code of the model is available at http://www.it.cas.cz/en/d3/l033/biophysics-cardiac-cells.

3. Results

3.1. -Ratio in Adult Rat Cardiomyocytes

The -ratio in adult rat ventricular cardiomyocytes was computed according to (3) using the values from different studies that are specified in Table 2. The experimental data published by Thomas et al. [16], Despa et al. [17], and Yang et al. [18] result in average ratios of 1.7, 4.3, and 6, respectively. However, the ratio computed from experimental data measured by Gadeberg et al. [15] is substantially higher and highly dependent on how the average values of and were assessed. If these values were assessed as the ratio of the average of amplitudes to the average of membrane capacitances for a given set of cells (Table 2, values denoted by asterisk), the resulting ratio is 25, as reported by Gadeberg et al. [15]. However, the same ratio, computed from the averages of current densities directly assessed in individual cells (see values in Table 2), amounts to 39.5. This suggests that the extent of the -ratio in adult rat ventricular myocytes might range between 1.7 and nearly 40.

3.2. Impact of -Ratio on Cellular Activity

To investigate the impact of the distribution of Na-Ca exchangers between the t-tubular and surface membranes on cellular electrophysiological activity, we used our quantitative model of the rat ventricular cardiomyocyte [19] with and of 65.02 pF and 34.98 pF, respectively. To meet the -ratio of 1.7 that resulted from experimental data by Thomas et al. [16] and of 25 as reported by Gadeberg et al. [15], the original value of t-tubular fraction of transporters in the model ( = 0.78 [19]) had to be changed to 0.48 and 0.93, respectively (see (4)).

The effect of such changes of on resting concentrations of Ca2+ in subsarcolemmal spaces (, ), cytosol (), and network compartment of sarcoplasmic reticulum () in unstimulated cell is illustrated in Figure 1. The increase of (that corresponds to a reduced number of Na-Ca exchangers assigned to the surface membrane) induced a rise of and consequently also of . This resulted in a higher Ca2+ load in SR. However, was only slightly affected, because the increase of Ca2+ extrusion from this space, induced by higher , was compensated by increased Ca2+ flux from the cytosol to the t-tubular subsarcolemmal space. Thus, the model showed that a redistribution of Na-Ca exchangers would lead to marked changes in net Ca2+ flux from the surface membrane to the t-tubular membrane through both subsarcolemmal spaces and cytosol (defined as Ca2+ cycling and analysed in detail in [19]).

To further explore the impact of changing on cellular electrophysiological activity, we used the model to simulate action potentials, membrane currents, and dynamic changes of Ca2+ concentration in individual intracellular compartments at set to 0.78, 0.48, and 0.93 at stimulation frequencies of 0.1, 1, and 5 Hz. The results obtained at 0.1 and 1 Hz (stimulation frequencies that are used in experimental works) and 5 Hz (stimulation frequency corresponding to resting heart rate in rats) are illustrated in Figures 2 and 3.

A redistribution of Na-Ca exchange proteins caused rather smaller but still apparent frequency-dependent change of action potential (AP). A decrease of from 0.78 to 0.48 led to a shortening of APD90 (duration of AP at 90% repolarisation) by approximately 7, 4, and 4% at 0.1, 1, and 5 Hz whereas increasing to 0.93 resulted in a prolongation of APD90 by 5, 4, and 2%, respectively. However, marked changes were evident in the level of and in the peak value of transient changes of . A decrease of to 0.48 led to a reduction of at the end of the cycle by 13, 7, and 3% and to a decrease of peak by 15, 8, and 4% at 0.1, 1, and 5 Hz, respectively. Conversely, an increase of to 0.93 caused a rise of by 14, 7, and 2% at the end of the cycle and an increase of peak by 23, 9, and 2%, respectively. To clarify the mechanisms underlying the described changes induced by redistribution of Na-Ca exchangers, the amounts of Ca2+ moving across sarcolemma, through cytosol, and across SR membrane within a steady-state cycle at 1 Hz stimulation are shown in Figure 4 (see Discussion for details).

These data show that a redistribution of Na-Ca exchangers causes relative changes in that decrease with the increase of stimulation rate (reflecting a decrease of relative changes of SR Ca2+ uptake at shorter stimulation cycle). This resulted in a reduction of the effect on at higher stimulation frequencies; the increase of from 0.48 to 0.93 caused an increase of peak by 45, 19, and 6% at 0.1, 1, and 5 Hz, respectively. Hence, the model indicates that the changes in SR Ca2+ load, cytosolic Ca2+ transient, and AP configuration caused by redistribution of Na-Ca exchangers between the t-tubular and surface membranes are potentially important in the whole range of the simulated frequencies explored.

4. Discussion

The assessment of the -ratio in adult rat ventricular myocytes from electrophysiological data obtained in intact and detubulated cells [1518] results in values ranging from 1.7 to nearly 40. The simulations on our model of rat ventricular cardiomyocytes incorporating separate description of ion transport in the t-tubular and surface membranes showed that such marked differences in the -ratio would have important consequences for intracellular Ca+ cycling, Ca2+ transient, and thereby inotropic state of cardiomyocytes.

4.1. Effect of -Ratio on Action Potential and Intracellular Transient

In 2014, we have shown that a change in the distribution of Na-Ca exchanger and other Ca2+ removal proteins (SERCA and sarcolemmal Ca ATPase) between membrane parts adjacent to dyadic and extradyadic spaces alters the amount of Ca2+ removed by each pathway from the cytoplasm and causes complex changes in intracellular Ca2+ dynamics and cellular Ca2+ cycling [23]. Despite that, a single relocation of Na-Ca exchangers (30%) from extradyadic to dyadic parts of membrane, explored at 5 Hz stimulation, had relatively small effects on Ca2+ transient in the cytosol (increase by ~3% [23]) due to the feedback effect of these changes on other Ca2+ transporting pathways.

The present study shows that changes in the distribution of Na-Ca exchanger between the t-tubular and surface membranes have substantially higher effects on cellular activity. An increase of from 0.48 to 0.93 caused a prolongation of APD90 by 12, 9, and 6% and a rise of steady-state cytosolic Ca2+ transient by 45, 19, and 6% at 0.1, 1, and 5 Hz, respectively. Thus, although this effect decreased at higher stimulation frequencies, it was still evident even at 5 Hz.

The analysis of simulated results showed that the observed changes in APD90 induced by different values of were caused predominantly by the related changes in (Figures 2 and 3). The higher the fraction of Na-Ca exchangers located in the t-tubular membrane, the higher (negative component) that was induced during AP due to the larger transient changes of Ca2+ concentration under the t-tubular membrane (compare the peak values of and in Figures 2 and 3). This explains the prolongation of APD90 with the increase of .

The increased extrusion of Ca2+ from the cell during AP would be expected to decrease the Ca2+ load in SR. Surprisingly, the simulations indicated that both quantities were increased or decreased simultaneously (Figures 2 and 3). To explain this apparent discrepancy, the Ca2+ turnover during stimulation cycle was explored in detail. As demonstrated in Figure 4, the increase in was accompanied by a decrease in the amount of Ca2+ extruded from the cell across the surface membrane (see the columns NaCa,s), which increased (see insets in Figure 2). The related increase of Ca2+ diffusion from the surface subsarcolemmal space into the cytosol (see in Figure 4(b)) promoted a higher Ca2+ uptake by the network SR () and an equivalent increase of Ca2+ release from junctional SR (). Note, however, that because of higher fraction of junctional SR at the t-tubules (0.8 [19]), the increase of Ca2+ release from SR was higher at its t-tubular side (see the differences in the rise of and with in Figure 4(c)). Simultaneously, the amount of Ca2+ cycling from the surface to t-tubular membrane increased (see the equal rise of and with in Figure 4), which, together with higher , prevented the reduction of (see insets in Figure 2) due to the increased Ca2+ extrusion through the t-tubular membrane. Consequently, the increase of Ca2+ diffusion from the surface subsarcolemmal space into the cytosol () with was not fully compensated by the decrease of Ca2+ diffusion from the t-tubular subsarcolemmal space into the cytosol (). Thus, the model shows that the relocation of Na-Ca exchangers from the surface to the t-tubular membrane changes the intracellular gradients of Ca2+ concentration (mainly because of the rise of ), which results in an increase of Ca2+ entering the cytosol and SR.

It is worth mentioning that a change of affected only minutely the total Ca2+ transfers through individual Ca transporters in the model (, , , and are nearly identical at different values of in Figure 4). The increased Ca2+ extrusion via during AP at higher values of was compensated by its reduction in the later phase of the stimulation cycle (see insets in Figure 2) due to lower (compared with ). This implicated that the total transsarcolemmal Ca2+ transfer was nearly unaffected by higher . Nevertheless, due to the increased Ca2+ cycling, an increased fraction of Na-Ca exchangers at the t-tubules may have a sizeable positive inotropic effect and may play an important role in the adaptation of the heart to the increase in hemodynamic demand during development. Besides, the inverse process could, at least partly, explain the relation between the reduction of cytosolic Ca2+ transient observed at 1 Hz stimulation in cells from failing hearts [21] and a decrease of -ratio in the same cells [15]. The frequency dependence of this effect, as revealed by the model, reflected a decrease in relative changes of SR Ca2+ uptake at shorter stimulation cycle; when was increased from 0.48 to 0.93, the amount of Ca2+ entering SR during a stimulation cycle increased by 69 and 28% at 0.1 and 1 Hz (see the increase in in Figure 4) and by 6% at 5 Hz.

To sum up, changes in the -ratio affect more the intracellular Ca2+ turnover and related excitation-contraction coupling than the electrical activity of rat ventricular cell in our model. However, the reason of sizeable differences in the -ratio as evaluated in cardiomyocytes from adult healthy rats remains to be explained.

4.2. Which Factors Might Be Responsible for Different -Ratio in Rat Cardiomyocytes?

The question about the factors affecting the -ratio in rat ventricular cardiomyocytes deserves attention. The basic characteristics of rats used by Gadeberg et al. [15], Thomas et al. [16], Despa et al. [17], and Yang et al. [18] are summarised in Table 3. As follows from the table, the largest differences between the animals appear to be in the age of rats. The rats used by Gadeberg et al. [15] were more than twice older than those used in all other studies. The studies published by Dan et al. [10] and Chen et al. [12] clearly showed that the distribution of Na-Ca exchangers in rabbit ventricular cardiomyocytes changes during development. While less than 10% of Na-Ca exchange proteins were present at the cell interior (t-tubular membrane) in myocytes from newborn rabbits, about 67% of these proteins were located there at their maturity (8 weeks postpartum) [10]. This distribution of Na-Ca exchangers results in an -ratio of 3.7 in our rat model, which is slightly smaller than the average ratio of 4.3 computed from data by Despa et al. [17] obtained from approximately eleven-week-old rats. However, the rats used by Thomas et al. [16], although being older than those used by Yang et al. [18], exhibited a smaller ratio (1.7 versus 6). This inconsistency in age-dependency of experimental data suggests that the -ratio in adult rat ventricular cardiomyocytes is likely affected by other factors. Besides the age of animals, a potentially important role may also be assigned to different experimental conditions and precision in measurement/evaluation of experimental data. The need to explain this physiologically important point strongly calls for new electrophysiological and immunolabelling data from adult rats of different ages.

5. Conclusion

The -ratio in adult rat ventricular cells, as assessed from electrophysiological data published to date, yields strikingly different values. The simulations performed on a model of rat cardiomyocyte showed that such differences in the -ratio would significantly affect the intracellular gradients of Ca2+ concentration, SR Ca2+ load, and thus the cellular inotropic state. The reason for these differences is however unclear. More studies focused on location and precise assessment of the fraction of Na-Ca exchangers within t-tubules would be of considerable value and would help to explain the differences resulting from available experimental data.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

This study was performed with the Institutional Support RVO: 61388998 and with the support from the Ministry of Health of the Czech Republic, Grant Project 16-30571A. Georges Christé’s work was supported by Université Lyon 1 Research Allowances to the EA4612 Neurocardiology Unit. The authors thank Professor C. Orchard for useful discussion.