- Open Access
Propagated repolarization of simulated action potentials in cardiac muscle and smooth muscle
Theoretical Biology and Medical Modellingvolume 2, Article number: 5 (2005)
Propagation of repolarization is a phenomenon that occurs in cardiac muscle. We wanted to test whether this phenomenon would also occur in our model of simulated action potentials (APs) of cardiac muscle (CM) and smooth muscle (SM) generated with the PSpice program.
A linear chain of 5 cells was used, with intracellular stimulation of cell #1 for the antegrade propagation and of cell #5 for the retrograde propagation. The hyperpolarizing stimulus parameters applied for termination of the AP in cell #5 were varied over a wide range in order to generate strength / duration (S/D) curves. Because it was not possible to insert a second "black box" (voltage-controlled current source) into the basic units representing segments of excitable membrane that would allow the cells to respond to small hyperpolarizing voltages, gap-junction (g.j.) channels had to be inserted between the cells, represented by inserting a resistor (Rgj) across the four cell junctions.
Application of sufficient hyperpolarizing current to cell #5 to bring its membrane potential (Vm) to within the range of the sigmoidal curve of the Na+ conductance (CM) or Ca++ conductance (SM) terminated the AP in cell #5 in an all-or-none fashion. If there were no g.j. channels (Rgj = ∞), then only cell #5 repolarized to its stable resting potential (RP; -80 mV for CM and -55 mV for SM). The positive junctional cleft potential (VJC) produced only a small hyperpolarization of cell #4. However, if many g.j. channels were inserted, more hyperpolarizing current was required (for a constant duration) to repolarize cell #5, but repolarization then propagated into cells 4, 3, 2, and 1. When duration of the pulses was varied, a typical S/D curve, characteristic of excitable membranes, was produced. The chronaxie measured from the S/D curve was about 1.0 ms, similar to that obtained for muscle membranes.
These experiments demonstrate that normal antegrade propagation of excitation can occur in the complete absence of g.j. channels, and therefore no low-resistance pathways between cells, by the electric field (negative VJC) developed in the narrow junctional clefts. Because it was not possible to insert a second black-box into the basic units that would allow the cells to respond to small hyperpolarizing voltages, only cell #5 (the cell injected with hyperpolarizing pulses) repolarized in an all-or-none manner. But addition of many g.j. channels allowed repolarization to propagate in a retrograde direction over all 5 cells.
There are no low-resistance connections between the cells in several different cardiac muscle and smooth muscle preparations [reviewed in refs.  and ]. In a computer simulation study of propagation in cardiac muscle, it was shown that the electric field (EF) that is generated in the narrow junctional clefts, when the prejunctional membrane fires an action potential (AP), depolarizes the postjunctional membrane to its threshold [3–5]. Propagation by mechanisms not requiring low-resistance connections have also been proposed by others [6–9]. This results in excitation of the postjunctional cell, after a brief junctional delay. The total propagation time consists primarily of the summed junctional delays. This results in a staircase-shaped propagation, the surface sarcolemma of each cell firing almost simultaneously . Propagation has been demonstrated to be discontinuous (or saltatory) in cardiac muscle [10–13]. Fast Na+ channels are localized in the junctional membranes of the intercalated disks of cardiac muscle [5, 14, 15], a requirement for the EF mechanism to work [1–5].
We recently modeled propagation of APs of cardiac muscle and smooth muscle using the PSpice program for circuit design and analysis [16–18]. Like the mathematical simulation published in 1977  and 1991 , the EF developed in the junctional clefts (negative VJC) was large and sufficient to allow transfer of excitation to the contiguous cell, without the requirement of gap-junction (g.j.) channels. Propagation of excitation can occur by the EF mechanism alone, even when the excitability of the cells was made low. In connexin-43 (heterozygous) and Cx40 knockout mice, propagation in the heart still occurs, but it is slowed [19–22] as predicted by our PSpice simulation study .
The present experiments were carried out to study propagated repolarization in this model of simulated action potentials (APs). Propagation of repolarization is a phenomenon that occurs in cardiac muscle . It has been shown that propagation of vasodilation occurs in the microvasculature , and that the endothelial cells are involved in the conduction of hyperpolarization and vasodilation in an artery . Therefore, our hypothesis was that propagated repolarization would also occur in our PSpice model.
The methods used and PSpice program (Cadence Co, Portland) have been described in detail previously, including the circuit [17, 18]. In brief, each cell was represented by four basic excitable units, two for the long surface membrane of the cell (one upward-facing and one downward-facing) and one basic unit for each of the two junctional membranes (left end of cell and right end) (Fig 1). The radial (shunt) resistance of the junctional cleft (RJC) was placed in the junctions between adjoining cells. The basic units were connected internally by the intracellular longitudinal resistance (ri). The basic units were connected externally with the extracellular resistance (RO), broken down into a longitudinal component (Rol) and a transverse (radial) component (Ror). RO was connected to ground as depicted in Figure 1. The circuit used for each unit was kept as simple as possible, using only those ion channels that set the resting potential (RP) and predominate during the rising phase and plateau phase of the AP.
The myocardial cell was assumed to be a cylinder 150 μm long and 16 μm in diameter, and the smooth muscle cell a cylinder 200 μm long and 5 μm diameter. Since in vascular smooth muscle (VSM), the muscle fibers run in a circular direction, if transverse velocity is calculated, the fiber diameter should be used. The values of the capacitive and the resistive elements in each basic unit were set to reflect the input resistance (ca 20 MΩ) and input capacitance (ca 100 pF) of the individual cells, and the junctional units were prorated, with respect to the surface units, based on relative areas represented. At rest, the resistance of K+ compared to Na+ (cardiac muscle) or Ca++ (smooth muscle) were set to give resting potentials (RPs) of -80 mV for cardiac muscle and -55 mV for smooth muscle. During excitation, the action potentials (APs) overshot to +32 mV and +11 mV, respectively.
Electrical stimulations (IS1) were always applied internally to the first cell of the chain (cell A1). Rectangular depolarizing current pulses of 0.25 nA amplitude and 0.50 ms duration were applied. The delay time before the IS1 pulse was applied was usually set to 1.0 ms in SM. A second stimulus (IS2) that was hyperpolarizing was applied to the inside of the last cell (A5) of the chain when the APs of all 5 cells were in their plateau phase. The intensity and duration of the IS2 pulses were varied over a wide range in order to generate strength / duration (S/D) curves.
Because the PSpice program does not have a voltage-dependent resistance (to generate the increase in Na+ or Ca++ conductance during excitation), this function had to be done with a V-controlled current source (our "black-box"). The sigmoidal relationship between conductance and membrane potential (VM), over a relatively narrow VM range, was mimicked by the black-box. The Na+ or Ca++ current required for excitation had to be calculated for several VM values and inserted into the GTABLE function.
Experiments were done with a single chain of 5 cells or 2 cells. There were no gap junctions between the cells of the chain under initial conditions. The presence of gap junction connexons (tunnels) was represented by adding a variable shunt resistance (Rgj) across each cell-to-cell junction. This resistor connected the inner surface of the prejunctional membrane with the inner surface of the postjunctional membrane. This Rgj shunt resistance was varied between 10,000 MΩ (1 tunnel), 1000 MΩ (10 tunnels in parallel), 100 MΩ (100 tunnels), 10 MΩ (1,000 tunnels), and 1.0 MΩ (10,000 tunnels). Each tunnel was assumed to have a conductance of 100 pS.
A. All-or- None Repolarization of Stimulated Cell A5
There was a sharp (all-or-none) repolarization of the stimulated cell (A5) of the 5-cell chain in both cardiac muscle (Fig. 2AB) and smooth muscle (Fig. 2CD). As shown, stimulation of cell A1 with a depolarizing current pulse (IS1) produced propagation of APs down the chain. At the plateau (peak) of the APs, a repolarizing pulse applied intracelluarly to cell A5, if of sufficient intensity (duration constant), produced a sudden repolarization of only cell A5 (Fig. 2B for cardiac muscle and D for smooth muscle). A slightly lower current intensity failed to produce a stable repolarization of cell A5 (Fig. 2A for cardiac muscle and 2C for smooth muscle). Note that the potential change (repolarizing) produced in neighboring call A4 was very small (< 1 mV). This emphasizes that there are indeed no low-resistance connections between the modeled cells under standard conditions. The hyperpolarizing pulse had to bring the Vm of cell A5 into the region of the GTABLE's sigmoidal curve. The transient repolarization is in agreement with the biological case [23–25].
B. Propagation of Repolarization
As indicated in the Methods section, it was not possible to insert a second black-box in the K+ leg of the basic circuit, because the PSpice program became erratic. Therefore, in order to achieve propagation of the repolarization of cell A5 in the retrograde direction, it was necessary to insert gap-junction channels between the cells of the chain (1, 10, 100, 1000, 10000 channels). This corresponded to adding resistive shunts between the cells across the junctions (Rgj) of 10000, 1000, 100, 10, and 1.0 MΩ (assuming each channel has a conductance of 100 pS).
The results of doing such an experiment are shown in Fig 3 for cardiac muscle (A – C) and for smooth muscle (D – F). When there were many channels (e.g. 10,000 in Fig. 3A and 3D or 1000 in Fig. 3B and 3E), the rising phase of the APs of all 5 cells were superimposed. This means that all 5 cells fired nearly simultaneously, as expected because of the high degree of low-resistance coupling. However, when a repolarizing current pulse was applied to cell A5, its repolarization spread to the neighboring cells. But the other cells did not repolarize simultaneously, as can be seen. Instead, there was a propagation of the repolarization at a certain velocity. This repolarization velocity became slower and slower as the number of channels was decreased. For example, with 100 channels (Fig. 3C and 3F), the propagated repolarization velocity was slower than with 1000 channels (Fig. 3B and 3E) or 10,000 channels (3A and 3D). With only 10 channels, the repolarization did not persist in either cardiac muscle or smooth muscle (not illustrated).
C. Strength/Duration Curves
The intensity (strength) and duration of the rectangular hyperpolarizing current pulses (IS2) applied to cell A5 were varied over a wide range in order to generate strength / duration curves. This was done when Rgj was infinite (i.e., 0 channels) and when Rgj was 10 MΩ (1000 channels) for strong coupling. The pulse duration was initially constant at 1.0 ms (near the chronaxie value) and then lowered to 0.5 ms and to 0.25 ms. The current intensity was varied until the sharp endpoint occurred, namely the stable repolarization of all cells in the chain. These results are plotted in Figure 4 for cardiac muscle and smooth muscle. Panel A is the strength / duration (S / D) curve for when Rgj was infinite (0 channels), and Panel B is the S/D curve for when Rgj was 10 MΩ(1000 channels). Note that the IS2 intensity was about 8–10-fold greater when the cells were well-coupled, because the applied hyperpolarizing current had to spread to all 5 cells of the chain. Regardless, the chronaxie values were about the same (ca. 1.0 ms).
In principle, the addition of a second black-box into the K+ leg of the basic circuit would allow the cell to repolarize in an all-or-none fashion to small repolarizing currents. When this was attempted, the program behaved erratically. So in the absence of g.j. channels, only the cell (A5) injected with repolarizing current (IS2) was able to repolarize in an all-or-none manner. The neighboring cell (A4) exhibited only a slight repolarization of <1 mV when cell A5 had repolarized completely back to the RP (-80 mV for CM and -55 mV for SM). This fact emphasized that there were no low-resistance connections between the cells under our initial conditions.
However, addition of 10,000, 1,000, or 100 g.j channels (corresponding to Rgj values of 1.0, 10, and 100 MΩ) did allow propagation of repolarization to occur. The borderline value was 10 g.j. channels (1000 MΩ Rgj), e.g., repolarization propagated part-way down the chain in SM and almost succeeded in CM. Of course, inserting the g.j. channels required that the IS2 repolarizing current applied be much greater. This is because the IS2 current had to spread down the entire chain, with the threshold current required to cause all cells to repolarize being determined by sufficient current entering distal cell A1 to repolarize it to the GTABLE sigmoidal region. Thus the proximal cells, like A5 and A4, became hyperpolarized beyond the level required for their repolarization.
The repolarizing IS2 current intensity required for the all-or-none repolarization was lower when the rectangular pulse duration was increased. This was true for both when only the injected cell A5 was repolarized (Rgj = ∞) and when all 5 cells repolarized (Rgj of 1.0, 10, and 100 MΩ). Plots of current intensity (ordinate) versus current duration (abscissa) gave the typical hyperbolic strength/ duration curve for excitable membranes. The chronaxie values were about 1.0 ms, for which a time constant τm of about 1.44 ms could be calculated. The S/D curves for the two conditions (Rgj = ∞ and Rgj = 10 MΩ) show that the current intensity required was about 8–10-fold greater when there were many gj-channels, in both CM and SM.
The calculated velocity for propagated repolarization (θr) varied with the number of gj-channels (Table 1), as expected. The more channels, the faster the velocity. For cardiac muscle, the θr was about 200 cm/s in the very well coupled case (10,000 channels) and about 50 cm/s in the less coupled case (100 channels) (Table 1). In all cases, the velocity for propagated repolarization (θr) was much lower than the velocity for antegrade propagation (θa.). In the 2-cell chain, the calculated velocities of propagated repolarization were similar to those for the 5-cell chain (Table 1).
The present study provides some new and important information about the PSpice simulations. First, it verifies that propagation (orthodromic) can occur in the complete absence of gap-junction channels, as previously reported [3, 4, 16–18]. Second, it demonstrates for the first time that activation of Na+ (in CM) or Ca++ (in SM) channels is reversible, by bringing Vm back to the level of the sigmoidal activation curve (GTABLE). Third, it shows for the first time that, in the PSpice model, the membranes exhibit the characteristic strength/duration curves. Fourth, it shows that the PSpice program has some serious limitations.
In summary, because of technical difficulties with the PSpice program, it was necessary to insert gj-channels in order to produce propagation of repolarization. Otherwise, only the modeled cell injected (A5) with the repolarizing IS2 current was able to repolarize. Since the potential change in the neighboring cell was only about 1 mV or less, this emphasizes that there were no low-resistance connections between the simulated cells under initial conditions. Propagation in the orthodromic direction occurs by the electric field (EF) discussed in previous papers (1–4, 16–18). The repolarizing IS2 current gave S/D plots that were typical rectangular hyperbolic curves for excitable membranes, with chronaxie values of about 1.0 ms, both for CM and SM. The calculated velocity for propagated repolarization was greater when the number of gj-channels was increased. The antidromic (retrograde) propagation velocity was usually considerably slower than the orthodromic (antegrade) propagation velocity for depolarization. The present findings do not necessarily imply that, in biological tissue, gap junctions are required for propagated repolarization to occur.
Sperelakis N, McConnell K: An electric field mechanism for transmission of excitation from cell to cell in cardiac muscle and smooth muscles. Research Advances in Biomedical Engineering. Edited by: Mohan RM. 2001, Global Research Network, 2: 39-66.
Sperelakis N, McConnell K: Electric field interactions between closely abutting excitable cells. IEEE-EMB. 2002, 21: 77-89.
Sperelakis N, Mann E: Evaluation of electric field changes in the cleft between excitable cells. J Theor Biol. 1997, 64: 71-96. 10.1016/0022-5193(77)90114-X.
Picone J, Sperelakis N, Mann JE: Expanded model of the electric field hypothesis for propagation in cardiac muscle. Math & Computer Modeling. 1991, 15: 17-35. 10.1016/0895-7177(91)90079-M.
Kucera JP, Rohr S, Rudy Y: Localization of sodium channels in intercalated disks modulates cardiac conduction. Circ Res. 2002, 91: 1176-1182. 10.1161/01.RES.0000046237.54156.0A.
Pertsov AM, Medvinski AB: Electric coupling in cells without highly permeable cell contacts. Biofizika. 1976, 21: 698-700.
Hogues H, Leon LJ, Roberge FA: A model for study of electric field interactions between cardiac myocytes. IEEE Trans Biomed Eng. 1992, 30: 1232-1243. 10.1109/10.184699.
Suenson M: Ephaptic impulse transmission between ventricular myocardial cells in vitro. Acta Physiol Scand. 1992, 120: 445-455.
Barr RC, Plonsey R: Electrophysiological interaction through the interstitial space between adjacent unmyelinated parallel fibers. Biophys J. 1992, 61: 1164-1175.
Spach MS, Miller WT, Geselowitz DB, Barr R, Kootsey JM, Johnson EA: The discontinuous nature of propagation in normal canine cardiac muscle. Evidence for recurrent discontinuity of intracellular resistance that affects the membrane currents. 1981, 48: 39-54.
Diaz PJ, Rudy Y, Plonsey R: Intercalated discs as a cause for discontinuous propagation in cardiac muscle. A theoretical simulation. Ann Biomed Eng. 1983, 11: 177-189.
Shaw RM, Rudy Y: Ionic mechanisms of propagation in cardiac tissue. Roles of the sodium and L-type calcium currents during reduced excitability and decreased gap junction coupling. Circ Res. 1997, 81: 727-741.
Henriquez AP, Vogel R, Muller-Borer BJ, Henriquez CS, Weingart R, Cascio WE: Influence of dynamic gap junction resistance on impulse propagation in ventricular myocardium: a computer simulation study. Biophys J. 2001, 81: 2112-2121.
Sperelakis N: Cable properties and propagation of action potentials. CH 18. Cell Physiology Source Book. 1995, Academic Press Publishers. San Diego, 245-254. 1
Cohen SA: Immunocytochemical localization of rH1 sodium channel in adult rat heart atria and ventricle. Presence in terminal intercalated disks. Circulation. 1994, 74: 1071-1096.
Sperelakis N, Editorial: An electric field mechanism for transmission of excitation between myocardial cells. Circ Res. 2002, 91: 985-987. 10.1161/01.RES.0000045656.34731.6D.
Sperelakis N, Ramasamy L: Propagation in cardiac muscle and smooth muscle based on electric field transmission at cell junctions: An analysis by PSpice. IEEE-EMB. 2002, 21: 130-143.
Sperelakis N: Combined electric field and gap junctions on propagation of action potentials in cardiac muscle and smooth muscle in PSpice simulation. J Electrocard. 2003, 36 (4): 279-293. 10.1016/j.jelectrocard.2003.08.001.
Morley GE, Vaidya D, Samie FH, Lo CW, Taffet SM, Delmar M, Jalife J: Characterization of conduction in the ventricles of normal and heterozygous Cx43 knockout mice using optical mapping. J Cardiovasc Electrophysiol. 1999, 10: 1361-1375.
Tamaddon HS, Vaidya D, Simon AM, Paul DL, Jalife J, Morley GE: High resolution optical mapping of the right bundle branch in connexin40 knockout mice reveals low conduction in the specialized conduction system. Circ Res. 2000, 87: 929-936.
Gutstein DE, Morley GE, Tamaddon H, Vaidya D, Schneider MD, Chen J, Chien KR, Stuhlmann H, Fishman GI: Conduction slowing and sudden arrythmic death in mice with cardiac restricted inactivation of connexin43. Circ Res. 2001, 88: 333-339.
Vaidya D, Tamaddon HS, Lo CW, Taffet SM, Delmar M, Morley GE: Nullmutation of connexin43 causes slow propagation of ventricular activation in the late stages of mouse embryonic development. Circ Res. 2001, 88: 1196-1202.
Hoffman BF, Cranefield PF: Electrophysiology of the Heart. 1960, McGraw-Hill Publishers
Figueroa XF, Paul DL, Simon AM: Central role of connexin40 in the propagation of electrically activated vasodilation in mouse cremasteric arterioles in vivo. Circ Res. 2003, 92: 793-800. 10.1161/01.RES.0000065918.90271.9A.
Emerson G, Segal S: Endothelial cell pathway for conduction of hyperpolarization and vasodilation along hamster feed artery. Circ Res. 2000, 86: 94-100.
The authors thank Cara Stevens for typing the manuscript.