Next Article in Journal
Model Systems to Study the Mechanism of Vascular Aging
Next Article in Special Issue
Facilitation of hERG Activation by Its Blocker: A Mechanism to Reduce Drug-Induced Proarrhythmic Risk
Previous Article in Journal
Ubiquitin Proteasome System Role in Diabetes-Induced Cardiomyopathy
Previous Article in Special Issue
The W101C KCNJ5 Mutation Induces Slower Pacing by Constitutively Active GIRK Channels in hiPSC-Derived Cardiomyocytes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model

Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Shiga 525-8577, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(20), 15378; https://doi.org/10.3390/ijms242015378
Submission received: 30 September 2023 / Accepted: 13 October 2023 / Published: 19 October 2023
(This article belongs to the Special Issue Cardiac Arrhythmia: Molecular Mechanisms and Therapeutic Strategies)

Abstract

:
Although repolarization has been suggested to propagate in cardiac tissue both theoretically and experimentally, it has been challenging to estimate how and to what extent the propagation of repolarization contributes to relaxation because repolarization only occurs in the course of membrane excitation in normal hearts. We established a mathematical model of a 1D strand of 600 myocytes stabilized at an equilibrium potential near the plateau potential level by introducing a sustained component of the late sodium current (INaL). By applying a hyperpolarizing stimulus to a small part of the strand, we succeeded in inducing repolarization which propagated along the strand at a velocity of 1~2 cm/s. The ionic mechanisms responsible for repolarization at the myocyte level, i.e., the deactivation of both the INaL and the L-type calcium current (ICaL), and the activation of the rapid component of delayed rectifier potassium current (IKr) and the inward rectifier potassium channel (IK1), were found to be important for the propagation of repolarization in the myocyte strand. Using an analogy with progressive activation of the sodium current (INa) in the propagation of excitation, regenerative activation of the predominant magnitude of IK1 makes the myocytes at the wave front start repolarization in succession through the electrical coupling via gap junction channels.

Graphical Abstract

1. Introduction

It has long been recognized that the repolarization of cardiac action potential (AP) can be regulated by both intrinsic ion channel properties, as well as cable properties conferred by the three-dimensional syncytium of cardiac myocytes [1,2,3,4]. However, the mechanisms of the propagation of repolarization in multicellular preparations of the cardiac muscle remain to be clarified based on ionic current systems. Initially, Weidman (1951) [5] observed all-or-none repolarizations in sheep’s Purkinje fiber when he applied a hyperpolarizing current beyond a certain strength during the plateau phase. He observed that there was a definite ‘threshold’ for the all-or-none repolarization, and suggested that the wave of repolarization could propagate without decrement as in the case of excitation. Essentially the same threshold phenomenon was observed by Cranefield and Hoffman (1958) [6] when inducing repolarization in the papillary muscle preparation. However, attempts to determine the velocity of the propagation of repolarization were not entirely successful. In part, this were due to the fact that the velocity of the repolarization wave was much slower than that of the excitation, so the propagation of repolarization was expected to be overtaken by the normal process of repolarization after the beginning of the AP.
The first theoretical study of cardiac electrophysiology using detailed mathematical descriptors was presented by Noble in 1962 [1]. He modified the Hodgkin–Huxley equations of squid giant axon and combined these with the cable theory to quite accurately reproduce many of the electrical properties of Purkinje fiber APs and pacemaker potentials. Purkinje fiber taken from the conduction system of juvenile sheep hearts had an exceptionally large diameter and length, which allowed extrapolation of the theory developed for the squid nerve to the cardiac muscle preparation after some modifications on parameters. After the epoch-making experimental milestones of the introduction of both dissociated cardiomyocyte preparations [7,8,9] and the patch clamp technique [10] in the next decades, the characteristics of ion channels and transporters were extensively clarified in the isolated cardiac myocytes from various cardiac tissues. These milestones made possible the development of various cardiac myocyte models based on the experimental data on isolated single myocytes [11]. Since then, cardiac myocyte models have become increasingly detailed and have been validated based on experimental data (for human ventricular myocyte models, ORd model [12], GPB model [13], Asakura et al., [14], Himeno et al. [15], etc.).
The aims of the present simulation study were (1) to clarify the ionic mechanisms of induced regenerative repolarization and the threshold phenomenon for the all-or-none repolarization in dissociated ventricular myocytes, (2) to evaluate the impacts of intercellular currents on the propagation of repolarization induced in the myocyte strand, and (3) to determine the velocity of the propagation of the repolarization wave. Accordingly, we first determined whether our comprehensive human ventricular myocyte (hVC) model [15] could reproduce the all-or-none repolarization as demonstrated by Trenor et al. (2017) [2] using the ORd model [12]. After confirming the fundamental behavior of all-or-none repolarization in the hVC model, we developed an in silico 1D strand of myocytes to analyze the effects of the intercellular currents on the regenerative repolarization as well as on the propagation of the repolarization wave. The clinical implications for the antiarrhythmic treatment of the present study will be described in the Section 3.

2. Results

2.1. Threshold Potentials for the All-or-None Repolarization in the hVC Model

Before examining the propagation of repolarization, we examined whether the all-or-none repolarization, i.e., the key experimental basis for suggesting the propagation of repolarization, can be observed in a single hVC model as has been proved in the ORd model by Trenor et al. (2017) [2]. In the simulation in Figure 1, the membrane potential (Vm) was clamped at various test potential (Vt) levels with 0.5 mV intervals for a duration of 5 ms.
After the membrane was released from the temporal voltage clamp, Vm changed automatically before progressive depolarization or hyperpolarization became obvious as seen in Figure 1A. The slope of this change reversed between −55 and −53.5 mV of the clamp pulse applied at 50 ms. By adjusting the Vt bit-by-bit, a relatively flat phase of Vm change, reflecting an unstable equilibrium, was obtained at around −50 mV. Above that level, delayed activation of an early after-depolarization (EAD)-like event started after the release time. A Vt far more positive than this threshold potential level induced the rising foot of an EAD-like event directly, while a more negative Vt induced monotonic repolarization. The threshold potentials, determined at an interval of 2 ms during the plateau phase, were plotted in Figure 1B (red dots). The threshold potential was time- and Vm-dependent; it shifted to more negative levels when the clamp pulse was applied earlier during the AP plateau. There were two apparent ‘jumps’ in the traces of the thresholds: at Vt = −42 and −23 mV. The green dots on the control time course of AP (black trace) indicated that the hyperpolarizing clamp pulse only accelerated the monotonic repolarization. All these findings were very similar to those reported in previous studies [2,3].

2.2. Ionic Mechanisms Responsible for Determining the Threshold Potentials for All-or-None Repolarization

As demonstrated in Figure 1A, the fate of all-or-none repolarization was determined by the initial slow Vm changes after switching off the clamp pulse (Figure 1A). To clarify the ionic mechanisms underlying the initial slow change, the changes in p(O) were plotted in Figure 2B. In both columns I and II, hyperpolarizing clamp pulses were applied at 50 ms. When the pulse induced repolarization to −55 mV, the p(O) of inward-going late sodium current (INaL), p(O)INaL, gradually decreased, and p(O) of the rapid component of delayed rectifier potassium current (IKr), p(O)IKr, also decreased. The magnitude of p(O) of the inward rectifier potassium channel (IK1), p(O)IK1, was slightly but significantly increased to enlarge IK1 (Figure 2C) by the negative jump in Vm from the plateau potential to −55 mV. The sum of the K+ currents was larger, by about 0.3 pA/pF, than the sum of the major inward currents, which were mainly given by INaL and the background non-selective cation current, IbNSC, to drive the gradual repolarization (sum of major Iout and total Iin were represented by gray lines in Figure 2C and 2D, respectively). Finally, the membrane quickly hyperpolarized to the resting potential via a positive feedback cycle; the rapid increase in p(O)IK1 drives repolarization and more repolarization induces an increase in p(O)IK1. Through this cycle, the amplitude of IK1 showed a sharp peak as large as 1.5 pA/pF at about −70 mV. From this fact, it may be concluded that the IK1 takes the central role in inducing regenerative repolarization toward the resting potential.
When the hyperpolarization from the plateau potential was increased only by 1 mV (to −54 mV), the slope of the slow depolarization was inverted from a negative (at −55 mV) to a positive direction (column II). As expected, no significant changes were obvious in the initial levels of p(O) or the amplitudes of INaL, IKr, IbNSC, and IK1 immediately after the break in the voltage clamp pulse in comparison to the pulse to −54 mV. However, an evident exponential rise in the L-type calcium current (ICaL) amplitude was observed 200 ms after the pulse (Figure 3D); the more the p(O)ICaL increased, the more depolarized the membrane became, and vice versa. The accompanying changes in INaL, IKr, IbNSC, and IK1 may occur secondarily due to the membrane depolarization caused by ICaL. The final discharge of an EAD-like event is caused by the much-enhanced positive feedback cycle of ICaL activation.
In column III, the Vt was −19 mV, which gave a stable p(O)ICaL of about 0.02 to supply an ICaL amplitude comparable to INaL at Vt = −54 mV (column II). On the other hand, the INaL amplitude was smaller if compared with that obtained by Vt at −54 mV because of Vm-dependent inactivation and the smaller driving force for the carrier ion Na+. The ICaL was also partially inactivated to give the decelerated rising phase of the EAD-like event. In column IV, the Vt (at 6 mV) was very close to the intrinsic plateau potential at 170 ms, and thereby the time courses of each current were almost the same as in the natural time course of final repolarization of the plateau. There was no threshold potential observed thereafter. It is important to note that the IKto is activated by the rapid depolarization, but is almost instantaneously deactivated by the hyperpolarization induced by the temporal voltage clamp pulse, thereby, IKto plays little role in initiating the full repolarization of the plateau (see Näbauer et al., 1996 [17]).
In summary, the IK1 takes the primary role in generating regenerative repolarization and the ICaL plays a role in regenerative depolarization. Both INaL and ICaL play modulatory roles in determining the threshold potential of repolarization.

2.3. Induction of a Stable Equilibrium Potential in the Current-Voltage Relationship of the hVC Model

In principle, if there is a second stable equilibrium potential in a more positive potential range than a resting membrane potential, as well as a negative conductance region in the steady-state I-V curve of a single myocyte model, it would be possible to induce a wave of repolarization in a multicellular preparation of myocyte model as predicted in the previous studies [1,2,3,4]. It is expected that the propagation of repolarization might be simply visualized in a 1D strand of such myocyte models whose Vm is set to another stable equilibrium potential. In brief, graded hyperpolarizing stimulus, which is strong enough to induce an all-or-none repolarization, is applied, it is expected that the repolarization would propagate as a wave based on the ionic mechanisms revealed for the threshold phenomenon at the single myocyte level (Figure 2): progressive increase in p(O)IK1 supplemented by deactivation of INaL and ICaL and activation of IKr.
A study by Chandler and Meves (1970) [18] used sodium fluoride to generate a long-lasting plateau in the AP of squid axons and reproduced the prolonged plateau potential mathematically by modifying the inactivation parameter of the sodium current (INa). In human ventricular myocytes, a slow inactivating sodium current component (INaL) has been described [16]. This INaL component in the human ventricle inactivates with a sufficiently fast time course to shape APD normally. However, it has been reported that the rate of slow inactivation of INaL is further delayed in heart failure, or by modulation of phospholipid associated with cardiovascular disease, to prolong the AP duration [19]. Based on these findings, we removed the inactivation of INaL to obtain a stable equilibrium potential near the plateau potential [14]. Figure 3D shows the removal of the slow inactivation steps to states (Is + I2), namely kIsf, kIsb, and kI1I2 were all set to zero. The fast inactivated state I1 allows flickering openings of the channel to reconstruct the repetitive brief openings demonstrated in the single channel recordings [16] because the magnitude of kI1O was set close to kOI1.
The whole-cell current–voltage (I-V) relationships in Figure 3A were obtained by using the full set of current systems (Equation (3)) except for the slow inactivated states of INaL fixed at 0.45 for Is and 0.001 for I2, respectively. The I-V curve represented by black points connected by lines clearly demonstrated the stable equilibrium (zero-current) potential at approximately −5 mV. Indeed, the slope of the whole-cell I-V curve was positive on both sides of the stable equilibrium potential. It should be noted that a similar positive slope was also evident around the stable equilibrium (resting) potential at −90.5 mV. In contrast, the slope of the steady-state I-V curve is negative around the unstable zero-current potential at −59 mV (vertical line).
To get deeper insights into the ionic mechanisms of the abolition of the AP plateau phase as well as the propagation of repolarization, we examined the composition of the whole-cell I-V relationship. In Figure 3A, the major determinants of I-V relations, such as IK1 (red), IKr (steel blue), INaL (blue), and ICaL (chocolate) were superimposed on the whole-cell current. Toward the negative potential side of the threshold potential (vertical line), a sharp voltage-dependent increase in both p(O)IK1 and a decrease in p(O)INaL (blue dots in Figure 3B) were evident, indicating that resultant increase in outward IK1 and decrease in inward INaL were mainly responsible for the acceleration of repolarization. On the positive side of the threshold potential, the sharp Vm-dependent increase in p(O)ICaL and the resultant increase in ICaL amplitude with depolarization promotes the depolarization. This influence of ICaL activation was antagonized by the voltage-dependent increase in p(O)IKr to give the stable equilibrium potential at about −5 mV.
The INaL took a key role in inducing the second stable equilibrium potential at the voltage range of AP plateau as demonstrated in Figure 3C1,C2, where the amplitude of INaL was decreased from one to zero in steps of 0.2. In Figure 3C1, it is evident that the unstable equilibrium potential was largely shifted toward negative potentials from −40 mV to −60 mV with increasing magnitude of INaL. If the limiting conductance of IK1 was increased or decreased by +20%, the unstable equilibrium potential was shifted only by +2 mV or −3 mV, respectively (not shown). In Figure 3C2, the ‘window component’ of ICaL was removed to reveal the role of INaL in inducing the unstable equilibrium potential. In the absence of ICaL, two intersections of the I-V relation with the zero-current level disappeared when INaL was intact, shown in in Figure 3C2, leaving the stable equilibrium potential at the resting potential intact. The addition of ICaL scarcely affected the location of unstable equilibrium potential when the relative magnitude of INaL remained in the physiological range (>0.2). If the gating kinetics of INaL (illustrated in Figure 3D) are intact under the physiological condition, the time-dependent slow inactivation accumulates the p(Is)INaL to promote membrane repolarization to the resting potential.

2.4. Propagation of Repolarization from the Second Stable Equilibrium Potential at Around 0 mV

To separate the propagation of repolarization from the natural time course of AP repolarization, the Vm of the myocyte strand was stabilized to the equilibrium potential near the plateau potential by fixing 45.1% of INaL channels to the inactivated states (Is and I2). The stimulus current pulse (−45 pA/pF, 30 ms in duration) was applied to cell No. 1 (Figure 4A) at 10 ms after starting the Vm recording. The first AP of cell No. 1 (blue trace) was conducted sequentially from myocyte to myocyte in the 1D strand model of hVCs (the AP depicted in red is recorded from cell No. 17). With the modified INaL (shown in Figure 3A,D), a train of EADs was initiated spontaneously during the repolarizing phase of the AP, but the EAD subsequently transformed into decremental oscillations in the myocyte strand to level off at a stable depolarization (Figure 4A) as predicted by the steady-state I-V curve.
The rate of AP conduction was measured from the slope of the profile of Vo along the 1D strand of in silico myocytes (Figure 4B). An example snapshot of the original Vo profile along cell No. on the x axis from cell No. 300~330 is demonstrated in the inset of Figure 4B. The Vo of a given myocyte within the myocyte strand is positive (here depicted in red) when the whole-cell membrane current is outward and vice versa (depicted in blue). It was evident that the direction of Io, which corresponded to the sign of Vo, was inverted from negative to positive between cell No. 317 and 318. The wavefront of excitation propagation (in the right-ward direction) was defined by the location of this inversion of the sign of the Io (defined as O/I center). In Figure 4B, the Vo of an individual myocyte was represented by the color code indicated at the bottom of Figure 4C, in which the ordinate is time, and the abscissa represents cell No. of the myocyte strand. The propagation rate of membrane excitation was determined to be 48.0 cm/s from the slope of the linear relationship between the propagation time and the number of the myocytes, whose length was 120 mm in the hVC model (Figure 4B).
The propagation of repolarization could be evoked in the 1D strand model with 600 hVCs which was depolarized continuously, as in Figure 4A. If compared with the propagation of excitation in panel B, the propagation of repolarization was characterized by a much smaller current density (by two orders of magnitude) of the major repolarizing current IK1 than that of the major excitatory current INa. To trigger the propagation of repolarization in the myocyte strand, the hyperpolarizing current should be applied either to a large fraction of myocytes within the strand, or at a large amplitude, or for a longer duration of current pulse to overcome the buffering action of the rest fraction of myocytes in the strand. This was indicated by the findings demonstrated in Figure 4C1, where a trigger current pulse of 400 pA/pF and 40 ms in duration was applied simultaneously to 100 myocytes (cell No. 50~150). Immediately after switching off the test pulse, the distance between the two O/I centers on the left and right sides at the borders of stimulated and unstimulated myocytes slightly decreased for ~10 ms before starting the propagation of repolarization. If the triggering pulse was smaller (<310 pA/pF), and/or shorter in duration, or applied to a smaller number of myocytes, the distance between the two O/I centers was shortened to push back the repolarization, which is consistent with the presence of a threshold phenomenon observed in the multicellular preparations.
The repolarizing electric stimulus applied to the 1D strand model of myocytes could be replaced by increasing the GK1, for example. Increasing the p(O)IK1 beyond 0.03 for 40 ms in the central 100-myocyte segment successfully triggered the propagation of repolarization. Similarly, the propagation of repolarization could be evoked by activation of the ATP-sensitive K+ channel (IKATP), which has a large maximum conductance (GKATP). Since the maximum conductance of IKr (GKr) is much smaller than that of IK1, increasing the p(O)IKr failed to evoke hyperpolarization to trigger the propagation of repolarization. Hyperpolarizing stimulus applied by nullifying the p(O) of the major inward currents, INaL and ICaL, during the AP plateau phase, was not effective in inducing propagation of repolarization.
The rate of the propagation of repolarization (the slope of the movement of the O/I center) slightly decreased with time and reached a steady rate of 1.7 cm/s after ~400 ms. The ionic mechanisms of repolarization propagation were analyzed in panels D1–D5 in Figure 4, which shows snapshots taken 800 ms after the initiation of the triggering pulse. Surprisingly, the Vo profiles in Figure 4D2 were distributed over a much wider range of ~100 myocytes, in which Vm continues to repolarize (Figure 4D1), in contrast to the narrower width of Vo profile over 10 myocytes in the case of propagation of membrane excitation (inset in Figure 4B). The location of the O/I center is indicated by the vertical line in Figure 4D1–D5. On the left side of the O/I center, the p(O)IK1 increased and p(O)INaL decreased with a steep slope (Figure 4D5), and the amplitude of IK1 increased while the amplitude of INaL decreased (Figure 4D4). The profile of all variables shown in panels D1–D5 remained the same as far as the rate of propagation was constant, indicating that the propagation proceeded without decrement.
To compare the contribution of membrane ionic currents (Im) and extracellular current (Io) to varying Vm of individual myocytes during the propagation of repolarization, profiles of all related parameters are illustrated on a common x-axis of cell No. in the strand.
I c = I o I m
where Ic (or Icap in Figure 4D3) is the current through the membrane capacitor to determine dVm/dt, and Im represents the sum of whole-cell ionic currents. Major ionic currents, such as IK1, IKr, INaL, and ICaL, and transporter currents INaK and INCX are depicted in different colors in Figure 4D4. It is evident that the peak amplitudes of these major currents are roughly in the same order of magnitude, indicating that each one of these major currents has the potential to affect the rate of propagation of+ repolarization. However, if the current amplitudes in Figure 4D4 were compared to the Vo or Io profiles in Figure 4D2 or Figure 4D3, it would seem to indicate that IK1 and INaL took the central role in driving repolarization over a potential range more negative than −60 mV. On the right side of the O/I center, the increase in p(O)IKr and decrease in p(O)ICaL over the segment of myocytes from No. 365 to 420 were mainly responsible for the initial repolarization over the Vm range from 0 to −60 mV. The currents generated by the Na/K pump (pink) and Na/Ca exchange (green) were relatively small and were mainly determined by the Vm change.
It is of note that the large increasing outward current component of IK1 and decreasing INaL in the myocytes on the left side of the O/I center with Vm below −60 mV (red and blue in Figure 4D4) not only contribute to repolarize themselves but also repolarize the myocytes on the right side of the O/I center (see Icap in Figure 4D3) through gap junction channels (see Io in Figure 4D3). The area of the positive and negative Io in Figure 4D3 was always equal and was distributed electrotonically along the strand. This electrotonic current started to repolarize the myocytes at the propagation wavefront where the inward and outward currents were still almost balanced at the potential near the quasi-equilibrium around 0 mV. These findings on current dynamics observed near the unstable equilibrium point around −60 mV were quite comparable, in qualitative aspects, with the ionic mechanisms observed in both the steady-state I-V curve (Figure 3) and the abolition experiment of the AP plateau (Figure 1 and Figure 2).

2.5. Effects of Varying the Conductance of INaL, IK1, and Gap Junction Channel (Gg) on the Rate of Propagation of Repolarization

The velocities of propagation for the repolarization were measured at various conductances of INaL and IK1 as listed in Table 1. The limiting conductance of INaL (GNaL) was increased by decreasing the magnitude of fixed (p(I2) + p(Is)), and vice versa when decreasing GNaL. The rate of propagation of repolarization increased with decreasing GNaL, suggesting that the influence of INaL to interfere with the hyperpolarizing influence of Io (= Icap + Im) in myocytes located in front of the O/I center is greater than the accelerating influence of INaL deactivation in myocytes located behind the O/I center. As expected, an increase in the conductance of outward IK1 (GK1) increased the propagation of repolarization to support the central role of IK1 in promoting the propagation.
In Table 1, the letter ‘F’ indicates that the triggering hyperpolarizing pulse failed to evoke the propagation of repolarization because of smaller GK1 and/or larger GNaL. Under these conditions, EAD-like events occurred in myocytes located before the wavefront after switching off the hyperpolarizing pulse. Then, the EAD-like event subsequently triggered an AP in myocytes behind the wavefront that stayed at the resting potential, resulting in a retrograde propagation of AP in the 1D strand.
These findings, obtained by varying the gap junction conductance (Gg), are consistent with the theory that the magnitude of the intercellular current is one of the determinants for promoting the rate of AP propagation (Table 2).

3. Discussion

3.1. Brief Summary of the Results

(1)
The ionic mechanisms underlying the induced all-or-none repolarization at the myocyte level (Figure 2), and the propagation of repolarization in a 1D linear strand of myocytes (Figure 4), were examined. Unambiguous propagation of repolarization that did not overlap with the natural time course of the AP plateau was observed only when the 1D strand of myocytes was depolarized at the second stable equilibrium potential, which was introduced by freezing the fraction of slow inactivation of INaL to 0.451 in Figure 4.
(2)
The key ionic currents were IK1, IKr, INaL, and ICaL in both the abolition of the AP plateau in a single ventricular myocyte and the propagation of repolarization in the in silico 1D strand of the myocyte model. Since the voltage range of the kinetic dynamics for IK1 and INaL are more negative compared to that for IKr and ICaL, as demonstrated in Figure 3A, IK1 and INaL play their role during the late phase of repolarization, while IKr and ICaL work in the early phase of repolarization. The positive feedback mechanisms accelerate the rate of Vm change caused by the kinetics dynamics of these currents.
(3)
In the abolition experiments in Figure 1 and Figure 2, the triggering pulse of hyperpolarization was supplied by the external current source, while the propagation of repolarization in the 1D strand was initiated by local current Io and evoked by the myocytes located behind the wavefront of repolarization before the O/I center (Figure 4D2).
(4)
The kinetic mechanisms of the four currents were most probably responsible for the stable propagation of repolarization.
(5)
The O/I profile of Vo around the O/I center (Figure 4D2) is a mirror image of that in the propagation of excitation (Figure 4B, inset). The width of the O/I profile was much more elongated in the propagation of repolarization (Figure 4D1–D5) compared with that for the propagation of AP (Figure 4B), reflecting the relatively small current amplitude of IK1, INaL, IKr, and ICaL in repolarization. The merit of using multiple channel species may be attributed to its flexibility in adjusting the configuration of the AP plateau phase more precisely and stably.
(6)
The rate of propagation of repolarization was dependent mainly on the amplitude of IK1 and INaL (Table 1). The rate increased with increasing the Gg within a limited range of Gg.

3.2. Clinical Implications

Noble et al. [1,3] demonstrated, for the first time, the principal ionic mechanisms underlying the induction of repolarization as well as the propagation of repolarization in cardiac myocytes, based on the computer model of cardiac membrane excitation. The present study used one of the detailed contemporary computer models (human ventricular myocyte models [12,13,14,15]) to clarify the ionic mechanisms of the propagation of repolarization in terms of individual ion channels as well as ion transporters to reconstruct the key experimental findings. The model was modified only by removing the very slow inactivation of INaL, leaving other fast mechanisms of INaT, ICaL, IKr, the transient outward potassium current (IKto), the slow component of delayed rectifier potassium current (IKs), and IK1 intact to create the stable equilibrium potential near 0 mV. Since the kinetic time constant of the slow inactivation of INaL is in the order of several hundreds ms, fixing the INaL slow inactivation hardly interferes with the repolarizing process from the stable equilibrium potential. Thus, the ionic mechanisms of the repolarization demonstrated in the present study are all relevant to repolarization mechanisms under physiological conditions.
The ionic mechanisms underlying the propagation of repolarization might be examined using an in silico model of a simple 1D strand, in which hVCs were connected by gap junctions. Here, let us assume that the 1D strand is depolarized by preventing its repolarization by removing the slow inactivation step from the state transition scheme of INaL (Figure 3D). This strand will be connected to an intact 1D strand of a similar size. At the junction of the two strands, the IK1 channel in the latter strand will provide a hyperpolarizing current through gap junctions to the depolarized myocytes to repolarize toward the resting potential. This repolarization through the cell-to-cell interaction should propagate into the depolarized strand as demonstrated in Figure 4. Through this propagation of repolarization, the whole strand will recover to its resting potential. Thus, the propagation of repolarization might be relevant in clinical anti-arrhythmic treatment. Indeed, the present study found that the propagation of repolarization could be initiated by amplifying GK1 above 0.03, which is at the same level of GK1 at the resting potential, in a fraction of the in silico 1D strand of hVC model. As expected, it was found in our preliminary study that magnifying GKATP was also useful in initiating the propagation of repolarization.
For the relevance of the present study on cardiac AP repolarization as well as the influences of these theoretical studies on clinical sciences, see the excellent review by Trenor et al. (2017) [2].

3.3. Ionic Mechanisms of Propagation of Repolarization in Comparison to the Propagation of Excitation

Here, a comprehensive view of the propagation of repolarization is discussed in comparison to that of excitation, including the threshold phenomena for the abolition of the AP plateau phase, the steady-state I-V relations, and positive feedback mechanisms. The mechanisms of the propagation of membrane excitation in the cardiac tissue were well established, as described in many textbooks. In essence, the depolarization-dependent activation of INa and the INa-dependent membrane depolarization cause the positive feedback cycle to promote the rapidly rising phase of AP. Depolarization of a myocyte within the 1D strand generates outward-going Io in neighbor myocytes through gap junctions to cause the positive cycle of INa activation in the next myocyte. In contrast, the propagation of membrane repolarization takes place through the positive feedback cycle of repolarization-dependent activation of IK1 and IK1-dependent repolarization in the presence of supplemental ionic currents. Repolarization of myocytes in the strand induces inward-going Io to repolarize myocytes which are still depolarized at the repolarization wavefront of the strand through gap junctions.
In the intact hVC model, all the parameters of INaL, including the slow inactivation states (I2 and Is) of the model, continuously change throughout the time course of the AP plateau phase driven by the time- and voltage-dependent gating parameters to promote the membrane repolarization. On the other hand, the state transitions among the C, O, and I1 states of INaL are much faster than the Vm change during the relatively slow AP plateau, and thereby these fast parameters roughly follow the steady-state values during the time course of the AP plateau phase. Thus, it is quite possible that the model may take nearly the same parameter set as that used in the steady-state I-V curve (Figure 3A) at an instant during the natural time course of the plateau phase.
It seems difficult to find any ionic current that can substitute the role of INaL in the simulation of propagated repolarization in the multicellular preparations. Especially as the unstable equilibrium (threshold) potential at around −60 mV allowed the sharp inversion of Io from inward to outward current at the O/I center, as demonstrated in Figure 4D1–D5. Note that these equilibrium potentials are common for both excitation and repolarization under physiological conditions (Figure 3). The window current of ICaL, which was depicted in chocolate in Figure 3A, can also provide the unstable equilibrium potential at about −40 mV in addition to the inducible stable equilibrium potential at −8.5 mV, even in the absence of INaL (Figure 3C1). However, around this threshold potential (−40 mV), the amplitude of the outward current (IK1) is much smaller than that at −60 mV, so the rate of repolarization may be too slow to evoke the regenerative propagation of repolarization against the buffering effect of the multicellular strand. In contrast, the threshold potential at ~−60 mV allowed the smooth and swift transition from the initial half of repolarization (ahead of the O/I center) to the final (behind the O/I center) repolarization due to the massive magnification of IK1 as shown in Figure 4D1–D5. The rapid upstroke of AP also starts at about −60 mV when the AP is triggered from the resting potential.
Qu et al. [20] in their review paper, introduced a general framework to explain the mechanisms of EADs using the quasi-stable I-V relationships. The same framework is also useful in understanding how the second stable equilibrium potential was introduced in the present study. The description of a given current, IX, is described using three parameters; the limiting conductance GX, the open probability of the channel (p(O)X), and the driving force (VmEX), where EX is the reversal potential of the ionic species of the current X.
I X = G X · p O · V m E X
The steady-state I-V curve is measured in experiments (or is computed in simulation in Figure 4) by clamping Vm until the Im reaches a stable level at different voltages along the Vm axis. While the quasi-steady-state I-V curve in the present study was obtained by the same protocol, the magnitude p(O) of a selected slow parameter was fixed during the natural time course of the AP plateau; for example, the p(O)INaL in the present study (Figure 4C). Thus, at a depolarized Vm during AP plateau, the quasi-steady-state whole-cell I-V relationship shown in Figure 4C moves upward smoothly with the continuous reduction of p(O)INaL accompanied by increasing inactivated fraction (p(Is)INaL + p(I2)INaL). At the plateau potential of the AP, the Vm traces the negative shift of the second stable equilibrium potential. Finally, the quasi-steady-state I-V relation fails to cross the X-axis over the depolarized potential range, and the membrane quickly repolarizes to the resting potential. Thus, the stable equilibrium potential could be induced in the 1D myocyte strand when the decrease of the inactivated fraction of INaL was frozen at an optimum level during the progress of the AP plateau with all the other kinetic parameters left intact.
Lastly, some of the reported values of the velocity of repolarization in previous studies are compared. In the classic study of Cranefield and Hoffmann [6], using papillary muscle fiber, they stated, “no value was given for conduction velocity under Results because of the uncertainty of measurement” and described a velocity of 20 cm/s estimated from the records of a propagated repolarization wave in their Discussion. In addition, they commented on the value as “this value may not be correct, but it is certain that the velocity is low compared with that of propagated depolarization” and referred to the original study of Weidemann [5]. Our simulation using a strand composed of intact hVC demonstrated that it was impossible to determine a unique value for the velocity of propagated repolarization because the velocity increased continuously and approached infinity as the repolarization propagated [21]. If the slope of the O/I center was roughly calculated to be t = 205~255 ms after the hyperpolarizing stimulation, a value of 17 cm/s was obtained, which was in the same order of magnitude as that reported by Cranefield and Hoffmann. However, it should be noted that this velocity is not a pure rate of the propagation of repolarization, instead it is a rate of the propagation of repolarization overlapping the intrinsic spontaneous repolarization of each myocyte. Previous work in this area is of interest. There has been a study on the conduction across the interface between two different cell groups using micropatterned cocultures of genetically engineered (GE) excitable cells and neonatal rat ventricular myocytes (NRVMs) [22]. The GE cell was a monoclonal HEK293 cell line stably expressing human INa (NaV 1.5), IK1 (Kir2.1), and rat connexin-43 (Cx43) gap junctions and had shorter APD (31.9 ms) than NRVMs (153.2 ms). They stimulated the co-cultured strand from the GE cell side and mapped the conduction of activation and repolarization using voltage-sensitive dye optically. In this strand, intrinsic repolarization occurred spontaneously in individual cells, but the GE cells of the former half triggered repolarization of the NRVMs composing the latter half of the strand because the APD of the GE cells was much shorter than NRVMs. If we calculate the rate of propagation of repolarization focusing on the donor–host interface from their result, the rate was about 1.8 cm/s which was much slower than that of activation (22 cm/s). In addition, Sperelakis et al. [23] measured the propagation of repolarization by inducing repolarization from another stable potential level to a resting potential as we did in the present study. They use a five-cell array of an extremely simplified model of electrical activity (2005). They reported that the rate of repolarization propagation is considerably (2.3~15 times) slower than that of excitation propagation, and that it varies with the number of gap junction channels (50~200 cm/s for repolarization propagation compared to 115~3000 cm/s for excitation propagation in the strand with 100~10,000 gap junction channels of 100 pS). Although the values previously reported in various experimental settings varied widely, the conclusion that the velocity of propagation of repolarization was much smaller than that of excitation was consistent in the present study.

4. Materials and Methods

4.1. The hVC Model

A single ventricular myocyte having a cuboid shape (120 µm on the longitudinal axis, 37.6 µm on the transverse axis, and 8.4 µm in depth) was assumed to calculate the propagation rate of excitation and repolarization. The membrane ionic current (Im) was given as a sum of 15 components as described in Equation (3).
I m = I NaT + I NaL + I CaL + I Kto + I Kr + I Ks + I K 1 + I Kpl + I bNSC + I Cab + I L Ca + I KATP + I NaK + I NaCa + I PMCA
As will be described in the present study, the IK1 takes a pivotal role in both the repolarization of the AP plateau as well as in the propagation of repolarization. In this simulation study, we used the IK1 model developed by the Ishihara group [24,25]. Fabbri et al. (2019) [26] compared six computer formulations of IK1 to overcome the deficiency of IK1 conductance in the human-induced pluripotent stem cells (hiPSC-CMs) through the dynamic clamp (DC) technique. They found that the optimum compensation of the AP waveform was satisfactory when the Ishihara IK1 model was used. The subendocardial type of IKto was used (Ito1 in Näbauer et al., 1996 [17]). The open probabilities of individual currents are plotted in Figure 2, Figure 3 and Figure 4. The open probability of whole cell ICaL (p(O)ICaL,X) was a weighted sum of those distributed in three Ca2+ spaces (Xs) in jnc, iz, and blk in the hVC model (Equation (4)).
p O I CaL = f r c C a L j n c · p O I CaL , j n c + f r c C a L i z · p O I CaL , i z + f r c C a L b l k · p O I CaL , b l k
where frcCaLX represents the fraction of ICaL channel distribution in X space (X: jnc, iz, and blk).
For calculating each current component in the hVC, abbreviations, physical constants, ionic composition of external solution, parameters, and initial values of the model components, see model equations in Supporting Materials and Himeno et al. (2015) [15] or Asakura et al. (2014) [14]. The dimensions of variables in the model equation were usually given as time ms, length µm, concentration mM (except [Ca2+]i µM for contraction), electric potential mV, conductance nS, and current pA, unless otherwise described. In the figures for membrane ionic currents in the present paper, the standard pA/pF is used.
The ordinary differential equations (ODE) of the hVC model [15] were time-integrated using the Heun method or Euler method. Computer programs were coded in Visual Basic and run on Visual studio. The time step of integration was set at 0.01 ms, approximately half of the maximum interval, which gives results consistent with those obtained using smaller step sizes. In some preliminary simulations, the Adams–Bashforth and Adams–Moulton integration methods were used to validate the numerical integration using simple methods.

4.2. One-Dimensional Strand of the hVC Model

Figure 5 shows the electrical coupling between two neighboring electric circuits in the myocyte model in a one-dimensional strand of Nc (=400) myocytes. The intracellular potential (Vi) was defined as the sum of the extracellular potential (Vo) and the membrane potential (Vm) of each myocyte (Equation (5)).
V i , n = V o , n + V m , n
The velocity of propagation along the fiber orientation was determined from the slope of the relationship between time and location of the point where the sign of the Vo was inverted.

4.3. Calculating Vo Outside the Myocyte Membrane in the Strand

The equivalent electric circuit of the myocyte strand is shown in Figure 5. The electrically equivalent circuit of the myocyte. Note that this myocyte strand model is different from the bidomain model because all myocytes within the strand have direct contact with the ground through an extracellular medium that has high conductance. This circuit represented experimental conditions of small multicellular preparations soaked in the bath solution, as in the classical experimental studies [5], rather than myocyte sheets tightly packed in the intact heart. In the presented model of the myocyte strand, the Vo of individual myocyte models was determined from individual Vm, Go, and Gg according to the algorithm as described in Algorithm 1.
Algorithm 1. Summary of the calculation for the 1D myocyte strand model.
Step1: Calculate conductance matrix of A, B, C in Equation (9);
Step2: Calculate the inverse matrix of the conductance matrix of A, B, C;
Step3: Calculate the right-hand side of Equation (9) from V m , n t ;
Step4 :   Calculate   V o , n t in Equation (9) by using the inverse matrix obtained in Step 2;
Step5 :   Calculate   V i , n t   from   V o , n t   and   V m , n t using Equation (5);
Step6 :   Calculate   I i , n t   from   V i , n t   and   V i , n + 1 t using Equation (6);
Step7 :   Calculate   I o , n t   from   I i , n t   and   I i , n 1 t using Equation (7);
Step8 :   Determine   d V m , n t d t   using   total   membrane   current   I m , n t   and   I o , n t
                          d V m , n t d t = I c , n t C m = I m , n t I o , n t C m ;
Step9: Integration of parameters for the next time step
                          V m , n t + d t = V m , n t + d V m , n t d t · d t ;
Step10: Go to Step3.
The current along the longitudinal axis (across the gap junction, Ii,n) was determined by Equation (6), which deals with the Gg and the voltage difference across the Gg.
I i , n = G g · V i , n V i , n + 1 = G g · V o , n + V m , n V o , n + 1 + V m , n + 1
The intercellular current along the transverse axis (Io) is driven by the difference between Ii,n−1 and Ii,n through Go.
I o , n = I i , n 1 I i , n = G o · V o , n
Replacing I in Equation (7) with the potentials (Vo, Vm) provides an equation composed of V and G,
G g · V o , n 1 + G o + 2 · G g · V o , n G g · V o , n + 1 = G g · V m , n 1 V m , n V m , n V m , n + 1
which is then rearranged for Vo of three sequential myocyte models (Vo,n−1, Vo,n, and Vo,n+1) as a function of voltages (Vo and Vm) as follows;
A · V o , n 1 + B · V o , n + C · V o , n + 1 = D · V m , n 1 + E · V m , n + F · V m , n + 1
here, the Vm of each unit is obtained beforehand by integrating the whole family of membrane currents to calculate the charge across the membrane capacitance. In brief, values of Vo for all myocytes are calculated by the parallel linear relations of all myocytes using a homemade solver every time after the numerical integration of the variables of individual myocytes. Vo was assigned a positive sign when Io was outward whereas Vo was assigned a negative sign when Io was inward and vice versa. As the Go (= 10 µS) was larger than the Gg (= 1.5 µS), the conduction velocity was largely dependent on Gg. The Go should be a scalable parameter, adjusted to the real size and dimension of the external electrode used for measurements. For a tentative Go magnitude, the recordings of the QRS complex using the epicardial electrogram [27] were referred to. With a Go of 5 µS, full-sized experimental Vo deflections of 5~10 mV were obtained with a standard maximum rate of rise of 250 V/s.
The gap junction current (Igap) was calculated on the right side of the nth equivalent circuit.
I gap = I i , n = G g · V i , n V i , n + 1 + E K , n
To calculate changes in [K+]i, Igap was assumed to be carried purely by K+. Thus, the Nernst equilibrium potential EK was determined across the gap junction.
E K , n = R T F · ln K + n K + n + 1

4.4. Electrical Stimulation to the Myocyte Model and the Myocyte Strand Model

The conditioning AP was evoked by applying a 5 ms current pulse of −12 pA/pF (K+-mediated current) in the single myocyte experiment. In the 1D strand model experiment, the stimulus current was magnified to −35 pA/pF in amplitude and 30 ms in duration to evoke the conditioning AP, which propagated throughout the axis to the right end of the strand.
To confirm whether the threshold phenomenon is observed in the single hVC model, hyperpolarizing voltage clamp pulses of 5 ms in duration at various test potential (Vt) levels in hyperpolarization were applied during the AP plateau phase. The compensation current given by Equation (12) was applied at various selected time points during the AP plateau phase.
I compensation = V t V m · 1 f b gain  
where fbgain was 1/(10 × dt), dt = 0.01 ms.
To trigger the propagation of repolarization in the 1D strand of myocytes, a trigger current pulse of 400 pA/pF and 40 ms in duration was applied simultaneously to 100 myocytes (cell No. 50~150).
These stimulus currents as well as Icompensation were added to the total membrane current Im,n of the myocyte model, which is used to calculate dVm/dt in each myocyte. The propagation of repolarization started from the border of the segment of the 100 myocytes. The movement of the wave of propagation was monitored by plotting the Vo in the 2D display (the myocyte number abscissa and the time ordinate).

4.5. Induction of the Second Stable Equilibrium Vm in Myocyte Model within the In Silico 1D Strand

The propagation of repolarization could be observed in a straightforward manner if the myocyte strand was depolarized to a stable equilibrium potential near the plateau potential around −5 mV. This was achieved by removing the slow Vm-dependent inactivation of INaL from the state transition scheme, described later in Figure 3D. The equilibrium potential was examined by recording the steady-state current–voltage relationship (I-V relationship) by applying long voltage clamp pulses of 4.9 s in duration. After confirming that the time-dependent changes were virtually saturated near the end of clamp pulses, the current amplitude was measured. Influences of varying the INaL were examined by changing the fractional amplitude of the fixed Is state of INaL.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms242015378/s1. Supplementary materials include model equations and computer program codes. References [14,15,24,28,29,30,31,32,33,34,35,36,37] are cited in the supplementary materials.

Author Contributions

Conceptualization, A.N. and Y.H.; methodology, A.N. and N.Y.; software, A.N., N.Y., M.U., H.N., Y.Z. and Y.M.; validation, S.E., Y.H., S.K., M.U., Y.Z. and A.A.; formal analysis, S.K. and S.E.; investigation, N.Y. and Y.Z.; writing—original draft preparation, A.N., Y.M. and Y.H.; writing—review and editing, Y.H.; visualization, N.Y., S.E. and H.N.; supervision, A.A. and Y.H.; project administration, Y.H.; funding acquisition, Y.H. and A.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by The Research Promotion Program for Acquiring Grants in-Aid for Scientific Research (KAKENHI) at Ritsumeikan University to Y.H. and A.A.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that supports the findings of this study are available in the Supplementary Material (Tables S1–S8 and Schemes S1–S6) of this article.

Acknowledgments

We would like to thank Haruka Hirata for her secretarial support.

Conflicts of Interest

The authors declare no conflict of interest.

References and Note

  1. Noble, D. A modification of the Hodgkin-Huxley equations applicable to Purkinje fibre action and pacemaker potentials. J. Physiol. 1962, 160, 317–352. [Google Scholar] [CrossRef] [PubMed]
  2. Trenor, B.; Cardona, K.; Saiz, J.; Noble, D.; Giles, W. Cardiac action potential repolarization revisited: Early repolarization shows all-or-none behaviour. J. Physiol. 2017, 595, 6599–6612. [Google Scholar] [CrossRef] [PubMed]
  3. Noble, D.; Hall, A.E. The Conditions for Initiating “All-or-Nothing” Repolarization in Cardiac Muscle. Biophys. J. 1963, 3, 261–274. [Google Scholar] [CrossRef] [PubMed]
  4. Noble, D. The voltage dependence of the cardiac membrane conductance. Biophys. J. 1962, 2, 381–393. [Google Scholar] [CrossRef] [PubMed]
  5. Weidmann, S. Effect of current flow on the membrane potential of cardiac muscle. J. Physiol. 1951, 115, 227–236. [Google Scholar] [CrossRef] [PubMed]
  6. Cranefield, P.F.; Hoffman, B.F. Propagated repolarization in heart muscle. J. Gen. Physiol. 1958, 41, 633–649. [Google Scholar] [CrossRef] [PubMed]
  7. Powell, T.; Steen, E.M.; Twist, V.W.; Woolf, N. Surface characteristics of cells isolated from adult rat myocardium. J. Mol. Cell. Cardiol. 1978, 10, 287–292. [Google Scholar] [CrossRef]
  8. Powell, T.; Twist, V.W. A rapid technique for the isolation and purification of adult cardiac muscle cells having respiratory control and a tolerance to calcium. Biochem. Biophys. Res. Commun. 1976, 72, 327–333. [Google Scholar] [CrossRef]
  9. Powell, T.; Terrar, D.A.; Twist, V.W. Electrical properties of individual cells isolated from adult rat ventricular myocardium. J. Physiol. 1980, 302, 131–153. [Google Scholar] [CrossRef]
  10. Sakmann, B.; Neher, E. Patch clamp techniques for studying ionic channels in excitable membranes. Annu. Rev. Physiol. 1984, 46, 455–472. [Google Scholar] [CrossRef]
  11. Fink, M.; Niederer, S.A.; Cherry, E.M.; Fenton, F.H.; Koivumäki, J.T.; Seemann, G.; Thul, R.; Zhang, H.; Sachse, F.B.; Beard, D.; et al. Cardiac cell modelling: Observations from the heart of the cardiac physiome project. Prog. Biophys. Mol. Biol. 2011, 104, 2–21. [Google Scholar] [CrossRef] [PubMed]
  12. 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. [Google Scholar] [CrossRef]
  13. Grandi, E.; Pasqualini, F.S.; Bers, D.M. A novel computational model of the human ventricular action potential and Ca transient. J. Mol. Cell. Cardiol. 2010, 48, 112–121. [Google Scholar] [CrossRef] [PubMed]
  14. Asakura, K.; Cha, C.Y.; Yamaoka, H.; Horikawa, Y.; Memida, H.; Powell, T.; Amano, A.; Noma, A. EAD and DAD mechanisms analyzed by developing a new human ventricular cell model. Prog. Biophys. Mol. Biol. 2014, 116, 11–24. [Google Scholar] [CrossRef]
  15. Himeno, Y.; Asakura, K.; Cha, C.Y.; Memida, H.; Powell, T.; Amano, A.; Noma, A. A human ventricular myocyte model with a refined representation of excitation-contraction coupling. Biophys. J. 2015, 109, 415–427. [Google Scholar] [CrossRef] [PubMed]
  16. Undrovinas, A.I.; Maltsev, V.A.; Kyle, J.W.; Silverman, N.; Sabbah, H.N. Gating of the late Na+ channel in normal and failing human myocardium. J. Mol. Cell. Cardiol. 2002, 34, 1477–1489. [Google Scholar] [CrossRef]
  17. Näbauer, M.; Beuckelmann, D.J.; Uberfuhr, P.; Steinbeck, G. Regional differences in current density and rate-dependent properties of the transient outward current in subepicardial and subendocardial myocytes of human left ventricle. Circulation 1996, 93, 168–177. [Google Scholar] [CrossRef]
  18. Chandler, W.K.; Meves, H. Slow changes in membrane permeability and long-lasting action potentials in axons perfused with fluoride solutions. J. Physiol. 1970, 211, 707–728. [Google Scholar] [CrossRef]
  19. Undrovinas, A.I.; Fleidervish, I.A.; Makielski, J.C. Inward sodium current at resting potentials in single cardiac myocytes induced by the ischemic metabolite lysophosphatidylcholine. Circ. Res. 1992, 71, 1231–1241. [Google Scholar] [CrossRef] [PubMed]
  20. Qu, Z.; Xie, L.H.; Olcese, R.; Karagueuzian, H.S.; Chen, P.S.; Garfinkel, A.; Weiss, J.N. Early afterdepolarizations in cardiac myocytes: Beyond reduced repolarization reserve. Cardiovasc. Res. 2013, 99, 6–15. [Google Scholar] [CrossRef]
  21. Himeno, Y. Ritsumeikan University, Shiga, Japan. Our simulation using a strand composed of intact hVC demonstrated that it was impossible to determine a unique value for the velocity of propagated repolarization because the velocity increased continuously and approached infinity as the repolarization propagated. 2023.
  22. Kirkton, R.D.; Badie, N.; Bursac, N. Spatial Profiles of Electrical Mismatch Determine Vulnerability to Conduction Failure Across a Host–Donor Cell Interface. Circ. Arrhythm. Electrophysiol. 2013, 6, 1200–1207. [Google Scholar] [CrossRef] [PubMed]
  23. Sperelakis, N.; Ramasamy, L.; Kalloor, B. Propagated repolarization of simulated action potentials in cardiac muscle and smooth muscle. Theor. Biol. Med. Model 2005, 2, 5. [Google Scholar] [CrossRef]
  24. Yan, D.H.; Ishihara, K. Two Kir2.1 channel populations with different sensitivities to Mg2+ and polyamine block: A model for the cardiac strong inward rectifier K+ channel. J. Physiol. 2005, 563, 725–744. [Google Scholar] [CrossRef] [PubMed]
  25. Ishihara, K.; Sarai, N.; Asakura, K.; Noma, A.; Matsuoka, S. Role of Mg2+ block of the inward rectifier K+ current in cardiac repolarization reserve: A quantitative simulation. J. Mol. Cell. Cardiol. 2009, 47, 76–84. [Google Scholar] [CrossRef]
  26. Fabbri, A.; Goversen, B.; Vos, M.A.; van Veen, T.A.B.; de Boer, T.P. Required GK1 to Suppress Automaticity of iPSC-CMs Depends Strongly on IK1 Model Structure. Biophys. J. 2019, 117, 2303–2315. [Google Scholar] [CrossRef] [PubMed]
  27. Taggart, P.; Sutton, P.; John, R.; Hayward, R.; Swanton, H. The epicardial electrogram: A quantitative assessment during balloon angioplasty incorporating monophasic action potential recordings. Br. Heart J. 1989, 62, 342–352. [Google Scholar] [CrossRef]
  28. Shannon, T.R.; Wang, F.; Puglisi, J.; Weber, C.; Bers, D.M. A mathematical treatment of integrated Ca dynamics within the ventricular myocyte. Biophys. J. 2004, 87, 3351–3371. [Google Scholar] [CrossRef]
  29. Negroni, J.A.; Lascano, E.C. Simulation of steady state and transient cardiac muscle response experiments with a Huxley-based contraction model. J. Mol. Cell. Cardiol. 2008, 45, 300–312. [Google Scholar] [CrossRef]
  30. Rush, S.; Larsen, H. A practical algorithm for solving dynamic membrane equations. IEEE Trans. Biomed. Eng. 1978, 25, 389–392. [Google Scholar] [CrossRef]
  31. Ishihara, K.; Yan, D.H. Low-affinity spermine block mediating outward currents through Kir2.1 and Kir2.2 inward rectifier potassium channels. J. Physiol. 2007, 583, 891–908. [Google Scholar] [CrossRef]
  32. Takeuchi, A.; Tatsumi, S.; Sara, N.; Terashima, K.; Matsuoka, S.; Noma, A. Ionic mechanisms of cardiac cell swelling induced by blocking Na+/K+ pump as revealed by experiments and simulation. J. Gen. Physiol. 2006, 128, 495–507. [Google Scholar] [CrossRef]
  33. Oka, C.; Cha, C.Y.; Noma, A. Characterization of the cardiac Na+/K+ pump by development of a comprehensive and mechanistic model. J. Theor. Biol. 2010, 265, 68–77. [Google Scholar] [CrossRef] [PubMed]
  34. Smith, N.P.; Crampin, E.J. Development of models of active ion transport for whole-cell modelling: Cardiac sodium-potassium pump as a case study. Prog. Biophys. Mol. Biol. 2004, 85, 387–405. [Google Scholar] [CrossRef] [PubMed]
  35. Grandi, E.; Pasqualini, F.S.; Pes, C.; Corsi, C.; Zaza, A.; Severi, S. Theoretical investigation of action potential duration dependence on extracellular Ca2+ in human cardiomyocytes. J. Mol. Cell. Cardiol. 2009, 46, 332–342. [Google Scholar] [CrossRef] [PubMed]
  36. Hinch, R.; Greenstein, J.L.; Tanskanen, A.J.; Xu, L.; Winslow, R.L. A simplified local control model of calcium-induced calcium release in cardiac ventricular myocytes. Biophys. J. 2004, 87, 3723–3736. [Google Scholar] [CrossRef]
  37. Tran, K.; Smith, N.P.; Loiselle, D.S.; Crampin, E.J. A thermodynamic model of the cardiac sarcoplasmic/endoplasmic Ca2+ (SERCA) pump. Biophys. J. 2009, 96, 2029–2042. [Google Scholar] [CrossRef]
Figure 1. Threshold potentials detected by applying Vt of 5 ms in duration to different levels in the hVC model. (A) The application of Vt to −53.5, −54, −54.5, and −55 mV induced dynamic changes in the time course of plateau potential as displayed from top to bottom (four traces), respectively. The black vertical line in panel A indicates the off time (t = 55 ms) of the 5 ms clamp pulse. (B) The threshold potentials at which the given clamp pulse induced abolition of the plateau potential (red) were plotted against the off time of various test pulses. Green dots indicated there were no threshold potential. In this range colored in green, Vt did not induce all-or-none repolarization but led to monotonic repolarization.
Figure 1. Threshold potentials detected by applying Vt of 5 ms in duration to different levels in the hVC model. (A) The application of Vt to −53.5, −54, −54.5, and −55 mV induced dynamic changes in the time course of plateau potential as displayed from top to bottom (four traces), respectively. The black vertical line in panel A indicates the off time (t = 55 ms) of the 5 ms clamp pulse. (B) The threshold potentials at which the given clamp pulse induced abolition of the plateau potential (red) were plotted against the off time of various test pulses. Green dots indicated there were no threshold potential. In this range colored in green, Vt did not induce all-or-none repolarization but led to monotonic repolarization.
Ijms 24 15378 g001
Figure 2. Ionic mechanisms underlying changes in Vm induced by applying 5 ms hyperpolarizing clamp pulses. Each time and voltage of the brief hyperpolarizing voltage clamp pulse applied were indicated at the top of each column. (A) indicates Vm, (B) indicates the open probabilities p(O) of currents, IK1 (red), IKr (steel blue), IKto (black), ICaL multiplied by 2 (chocolate), and INaL multiplied by 100 (blue), (C,D) indicate amplitudes of outward (Iout) and inward (Iin) currents, respectively. The same colors were used as in (B). The sum of major outward and inward currents was calculated and illustrated in gray. The amplitude of INaK (pink) and IbNSC (black) are also shown in (C) and (D), respectively. The vertical lines indicate the end of the clamp pulse.
Figure 2. Ionic mechanisms underlying changes in Vm induced by applying 5 ms hyperpolarizing clamp pulses. Each time and voltage of the brief hyperpolarizing voltage clamp pulse applied were indicated at the top of each column. (A) indicates Vm, (B) indicates the open probabilities p(O) of currents, IK1 (red), IKr (steel blue), IKto (black), ICaL multiplied by 2 (chocolate), and INaL multiplied by 100 (blue), (C,D) indicate amplitudes of outward (Iout) and inward (Iin) currents, respectively. The same colors were used as in (B). The sum of major outward and inward currents was calculated and illustrated in gray. The amplitude of INaK (pink) and IbNSC (black) are also shown in (C) and (D), respectively. The vertical lines indicate the end of the clamp pulse.
Ijms 24 15378 g002
Figure 3. The N-shaped I-V curve with two stable equilibrium potentials and one unstable equilibrium potential of the hVC model. The I-V curves in (A,C1,C2) were measured by applying the voltage clamp test pulses of 4.9 s in duration given every 2 mV step from the holding potential of −80 mV. Current magnitudes saturated within the duration of test pulses at the end of individual test pulses, and were plotted in these figures. Additionally, (A) shows total whole cell current (black) and current components IK1, IKr, INaL, and ICaL depicted in different colors as indicated. The steady-state p(O)s of major ion channels were plotted in (B). Note that p(O)INaL and p(O)ICaL were plotted after multiplication by 100 and 10, respectively. (C1,C2) shows the influence of increasing the amplitude of INaL on the appearance of equilibrium points in the presence and the absence of ICaL, respectively. Increasing the amplitude of INaL shifted the N-shaped I-V curve negatively. All I-V and p(O)-V relationships share the same abscissa of Vm. In (D), the state of slow inactivation, Is and the transitional inactivated state, I2 of INaL, were shaded by gray color to indicate that these states were frozen, namely, the state transition from other states was totally prevented as indicated by cross marks. Thereby, repetitive state transitions between I1 and O produce continuous bursts of brief openings of the channel, and ‘late scattered mode’ of Na channel [16] (kI1O is close to kOI1) during a clamp pulse to depolarized potentials. If integrated over bursts with numerous numbers of brief openings, the whole-cell INaL maintains a significant amplitude during depolarization. With an intact transition scheme, INaL gradually decreases in amplitude through the state transition from O to (Is + I2) when the membrane is depolarized. The state transition between the two closed states (C1,C2) was assumed to be instantaneous on the Vm jump; the probability of (C2) increases with membrane depolarization. In the state transition of transient sodium current (INaT) (see Supplementary Materials) the I1 state is absent, and the inactivation occurs through a state transition from O to I2 at transition rate kOI2, which has the same magnitude as kOI1.
Figure 3. The N-shaped I-V curve with two stable equilibrium potentials and one unstable equilibrium potential of the hVC model. The I-V curves in (A,C1,C2) were measured by applying the voltage clamp test pulses of 4.9 s in duration given every 2 mV step from the holding potential of −80 mV. Current magnitudes saturated within the duration of test pulses at the end of individual test pulses, and were plotted in these figures. Additionally, (A) shows total whole cell current (black) and current components IK1, IKr, INaL, and ICaL depicted in different colors as indicated. The steady-state p(O)s of major ion channels were plotted in (B). Note that p(O)INaL and p(O)ICaL were plotted after multiplication by 100 and 10, respectively. (C1,C2) shows the influence of increasing the amplitude of INaL on the appearance of equilibrium points in the presence and the absence of ICaL, respectively. Increasing the amplitude of INaL shifted the N-shaped I-V curve negatively. All I-V and p(O)-V relationships share the same abscissa of Vm. In (D), the state of slow inactivation, Is and the transitional inactivated state, I2 of INaL, were shaded by gray color to indicate that these states were frozen, namely, the state transition from other states was totally prevented as indicated by cross marks. Thereby, repetitive state transitions between I1 and O produce continuous bursts of brief openings of the channel, and ‘late scattered mode’ of Na channel [16] (kI1O is close to kOI1) during a clamp pulse to depolarized potentials. If integrated over bursts with numerous numbers of brief openings, the whole-cell INaL maintains a significant amplitude during depolarization. With an intact transition scheme, INaL gradually decreases in amplitude through the state transition from O to (Is + I2) when the membrane is depolarized. The state transition between the two closed states (C1,C2) was assumed to be instantaneous on the Vm jump; the probability of (C2) increases with membrane depolarization. In the state transition of transient sodium current (INaT) (see Supplementary Materials) the I1 state is absent, and the inactivation occurs through a state transition from O to I2 at transition rate kOI2, which has the same magnitude as kOI1.
Ijms 24 15378 g003
Figure 4. Propagated repolarization in a linear strand of the EAD-prone myocytes in comparison to that of excitation. Panel A shows the time course of Vm, moving from the resting membrane potential (−91 mV) to the second stable equilibrium potential (−2 mV) in the myocyte strand, in which the deeply inactivated states (Is + I2) of INaL were removed. Panel (A) indicates Vm recordings in cell No. 1 (blue) and 17 (red). The AP was evoked by applying a depolarizing pulse to cell No. 1 in the 1D strand of the hVC model. The Vm stabilized at about −3 mV at the end of the record in all myocytes within the strand. Panel (B) indicates the propagation of the O/I pattern (propagation of excitation) along the myocyte strand. The inset shows the stable O/I pattern of Vo; positive Vo evoked by outward whole-cell current, Io (red), or negative Vo due to inward Io (blue). Note, that the cell number on the abscissa of the inset is the same as in panel B, but the scale of the Vo is reduced by 100 times. Panels (B,C) indicate the movement of O/I patterns (propagation of repolarization) along the myocyte strand by using the color code of Vo shown at the bottom of panel (C). The O/I patterns illustrated in (C1,C2) were obtained at different recording times. See text for the explanation of panels (D1D5). Note that p(O)INaL and p(O)ICaL were plotted after multiplication by 100 and 10, respectively, in panel (D5). The time scale shown in (D5) is common for all (D1D5). The arrow head in (D1) indicates the direction of the propagation of repolarization. The scaling factor of IK1 = 1.2 was used in this simulation. Note that the scale of the Vo is reduced by 100 times in panel (C) if compared to that in panel (B).
Figure 4. Propagated repolarization in a linear strand of the EAD-prone myocytes in comparison to that of excitation. Panel A shows the time course of Vm, moving from the resting membrane potential (−91 mV) to the second stable equilibrium potential (−2 mV) in the myocyte strand, in which the deeply inactivated states (Is + I2) of INaL were removed. Panel (A) indicates Vm recordings in cell No. 1 (blue) and 17 (red). The AP was evoked by applying a depolarizing pulse to cell No. 1 in the 1D strand of the hVC model. The Vm stabilized at about −3 mV at the end of the record in all myocytes within the strand. Panel (B) indicates the propagation of the O/I pattern (propagation of excitation) along the myocyte strand. The inset shows the stable O/I pattern of Vo; positive Vo evoked by outward whole-cell current, Io (red), or negative Vo due to inward Io (blue). Note, that the cell number on the abscissa of the inset is the same as in panel B, but the scale of the Vo is reduced by 100 times. Panels (B,C) indicate the movement of O/I patterns (propagation of repolarization) along the myocyte strand by using the color code of Vo shown at the bottom of panel (C). The O/I patterns illustrated in (C1,C2) were obtained at different recording times. See text for the explanation of panels (D1D5). Note that p(O)INaL and p(O)ICaL were plotted after multiplication by 100 and 10, respectively, in panel (D5). The time scale shown in (D5) is common for all (D1D5). The arrow head in (D1) indicates the direction of the propagation of repolarization. The scaling factor of IK1 = 1.2 was used in this simulation. Note that the scale of the Vo is reduced by 100 times in panel (C) if compared to that in panel (B).
Ijms 24 15378 g004
Figure 5. The electrically equivalent circuit of the myocyte strand. Go; the extracellular conductance of 10 µS, Vo; extracellular potential near the cell surface, Cm; membrane capacitance of 192.5 pF, Vm; membrane potential, Vi; intracellular potential, Io; extracellular current, Ic; capacitor current, Im; membrane current, Ii; gap junction current, Gg; the gap junction conductance of 1.5 µS assumed for the end-to-end junction along the longitudinal axis of the cuboid shape myocyte.
Figure 5. The electrically equivalent circuit of the myocyte strand. Go; the extracellular conductance of 10 µS, Vo; extracellular potential near the cell surface, Cm; membrane capacitance of 192.5 pF, Vm; membrane potential, Vi; intracellular potential, Io; extracellular current, Ic; capacitor current, Im; membrane current, Ii; gap junction current, Gg; the gap junction conductance of 1.5 µS assumed for the end-to-end junction along the longitudinal axis of the cuboid shape myocyte.
Ijms 24 15378 g005
Table 1. The rate of propagation of repolarization in the strand model. The sign ‘F’ indicates that the stimulus failed to evoke the propagation of repolarization. Note that increasing the fraction of I2 in (I2 + Is) by 10 or 20% decreases GNaL.
Table 1. The rate of propagation of repolarization in the strand model. The sign ‘F’ indicates that the stimulus failed to evoke the propagation of repolarization. Note that increasing the fraction of I2 in (I2 + Is) by 10 or 20% decreases GNaL.
0.639 (−20%)0.594 (−10%)Relative GNaL 0.549 after Fixing (I2 + Is)0.504 (+10%)0.459 (+20%)
−10%FFFFF
control GK1FFF1.692.11
+10%F1.521.942.382.94
+20%1.792.1142.402.693.04
Table 2. The rate of propagation of excitation and repolarization calculated by varying gap junction conductances (Gg) in the strand model.
Table 2. The rate of propagation of excitation and repolarization calculated by varying gap junction conductances (Gg) in the strand model.
Gg (nS)1000150020002500
rate excitation propagation (cm/s)36.548.051.7958.5
rate repolarization propagation (cm/s)1.932.432.602.94
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Himeno, Y.; Zhang, Y.; Enomoto, S.; Nomura, H.; Yamamoto, N.; Kiyokawa, S.; Ujihara, M.; Muangkram, Y.; Noma, A.; Amano, A. Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model. Int. J. Mol. Sci. 2023, 24, 15378. https://doi.org/10.3390/ijms242015378

AMA Style

Himeno Y, Zhang Y, Enomoto S, Nomura H, Yamamoto N, Kiyokawa S, Ujihara M, Muangkram Y, Noma A, Amano A. Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model. International Journal of Molecular Sciences. 2023; 24(20):15378. https://doi.org/10.3390/ijms242015378

Chicago/Turabian Style

Himeno, Yukiko, Yixin Zhang, Suzuka Enomoto, Hiroto Nomura, Natsuki Yamamoto, Shotaro Kiyokawa, Mirei Ujihara, Yuttamol Muangkram, Akinori Noma, and Akira Amano. 2023. "Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model" International Journal of Molecular Sciences 24, no. 20: 15378. https://doi.org/10.3390/ijms242015378

APA Style

Himeno, Y., Zhang, Y., Enomoto, S., Nomura, H., Yamamoto, N., Kiyokawa, S., Ujihara, M., Muangkram, Y., Noma, A., & Amano, A. (2023). Ionic Mechanisms of Propagated Repolarization in a One-Dimensional Strand of Human Ventricular Myocyte Model. International Journal of Molecular Sciences, 24(20), 15378. https://doi.org/10.3390/ijms242015378

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop