Boundary effects influence velocity of transverse propagation of simulated cardiac action potentials
© Sperelakis et al; licensee BioMed Central Ltd. 2005
Received: 18 July 2005
Accepted: 06 September 2005
Published: 06 September 2005
We previously demonstrated that transverse propagation of excitation (cardiac action potentials simulated with PSpice) could occur in the absence of low-resistance connections (gap – junction channels) between parallel chains of myocardial cells. The transverse transmission of excitation between the chains was strongly dependent on the longitudinal resistance of the interstitial fluid space between the chains: the higher this resistance, the closer the packing of the parallel chains within the bundle. The earlier experiments were carried out with 2-dimensional sheets of cells: 2 × 3, 3 × 4, and 5 × 5 models (where the first number is the number of parallel chains and the second is the number of cells in each chain). The purpose of the present study was to enlarge the model size to 7 × 7, thus enabling the transverse velocities to be compared in models of different sizes (where all circuit parameters are identical in all models). This procedure should enable the significance of the role of edge (boundary) effects in transverse propagation to be determined.
It was found that transverse velocity increased with increase in model size. This held true whether stimulation was applied to the entire first chain of cells or only to the first cell of the first chain. It also held true for retrograde propagation (stimulation of the last chain). The transverse resistance at the two ends of the bundle had almost no effect on transverse velocity until it was increased to very high values (e.g., 100 or 1,000 megohms).
Because the larger the model size, the smaller the relative edge area, we conclude that the edge effects slow the transverse velocity.
KeywordsPropagation of cardiac action potential transverse propagation velocity PSpice simulations edge/boundary effects electric field transmission of excitation.
Computer simulation of the propagation of impulses in cardiac muscle shows that the electric field generated in the narrow junctional clefts when an action potential occurs at the prejunctional membrane depolarizes the postjunctional membrane to threshold . Thus, the postjunctional cell is excited after a brief delay at the junction and propagation in cardiac muscle is saltatory. We have modeled APs in this tissue using the PSpice program for circuit design and analysis, and we have corroborated earlier reports that the EF developed in the junctional cleft is sufficiently large to allow transfer of excitation to the contiguous cell without the requirement for a gap-junction [2–6]. To date, however, we have only used small-sized models for these simulation studies.
When our paper on transverse propagation of cardiac action potential (APs) simulated by PSpice in a 5 × 5 model [4, 5] was reviewed by the journal, one unanswered question was whether edge (boundary) effects were important. The purpose of the present study was to address this question. To do this, we expanded the model to a 7 × 7 size (7 parallel chains of 7 cells each). Thus, we could compare transverse velocity in 2-dimensional models of 4 sizes: 7 × 7, 5 × 5, 3 × 4, and 2 × 3. It was essential that all circuit parameters were the same in all four models. It was found that the larger the model, the faster the transverse velocity of propagation, up to a presumed saturation point.
Transverse Propagation Velocity (antegrade (A) and retrograde(R)) of Simulated Cardiac Action Potentials in 2-D Sheets at a Rol2 of 500 KΩ.
No. of Chains Responding
Transv. Velocity cm/sec
7 × 7
Entire A Chain
Cell A1 Only
Entire G Chain
Cell G1 Only
6 (A failed)
5 × 5
Entire A Chain
Cell A1 Only
Entire E Chain
Cell E1 Only
3 × 4
Entire A Chain
Cell A1 Only
Entire C Chain
Cell C1 Only
2 × 3
Entire A Chain
Cell A1 Only
Entire B Chain
Cell B1 Only
The myocardial cell was assumed to be a cylinder 150 μm long and 16 μm in diameter. The cell capacitance was assumed to be 100 pF, and the input resistance to be 20 MΩ. A junctional tortuosity (interdigitation) factor of 4 was assumed for the cell junction [1, 2]. The junctional cleft potential (Vjc) is produced across Rjc, the radial resistance of the narrow and tortuous junctional cleft. The junctional cleft contained two longitudinal resistances of 7Ω each and two radial resistances (Rjc) of 50 MΩ each in parallel.
The tortuosity factor does not interact with the packing factor. The tortuosity factor concerns the complex interdigitation of contiguous cells longitudinally (end-to-end), whereas the packing factor deals with how closely the cell chains are packed transversely (or radially) within a tissue bundle. The value assigned to Rol2 reflects the closeness of this packing. The value assigned to Rjc reflects the thickness of the junctional gap (end-to-end) and the tortuosity factor.
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 of the AP. We wanted only to inscribe the rising phase of the APs to study propagation in the 2-dimensional sheet. The RP was -80 mV and the overshoot potential was +30 mV (AP amplitude of 110 mV). Transverse propagation velocity was calculated from the measured total propagation time (TPT) (measured as the difference between when the APs of the first cell and last cell crossed -20 mV) and cell width (number of chains minus one gives the number of transverse junctions traversed).
Because the PSpice program does not have a V-dependent resistance to represent the increase in conductance for Na+ ions in myocardial cells during depolarization and excitation, this function was simulated by a V-controlled current source (our "black-box") in each of the basic circuit units. The current output of the black-box at various membrane voltages was calculated assuming a sigmoidal relationship between membrane voltage and resistance over the range of -60 mV to -30 mV. The V values used in the GTABLE were those recorded directly across the membrane. The excitability of the basic units was the same as in our previous papers, i.e., it was set at the moderate level .
The upper chain of cells was assumed to be bathed in a large volume of Ringer solution connected to ground. The external resistance (Ro) of this fluid was divided into two components: a radial resistance (Ror) and a longitudinal resistance (Rol). The longitudinal resistance value between the chains (Rol2) was increased over a wide range to reflect closer packing of parallel chains into a bundle of fibers (Fig. 1). The transverse resistance of the interstitial fluid space (Ror2) was found to have almost no effect on the transverse velocity. The cells in each chain were not interconnected by low-resistance pathways (gap-junction channels), so that transmission of excitation from one cell to the next had to be by the electric field (EF) developed in the narrow junctional cleft. In our previous papers, we presented a number of references demonstrating that propagation velocity is slowed only slightly in the absence or paucity of gap junctions [e.g., see refs [1, 3] and ]. There were seven parallel chains (chains A-G) of seven cells each in the 7 × 7 model. The block diagrams and detailed circuitry for the other models (5 × 5; 3 × 4; 2 × 3) were previously published [4, 5]. The ends of each chain had a bundle termination resistance (RBT) of 1.0 KΩ to mimic the physiological condition. However, variation of RBT over a wide range had almost no effect, until very high values of about 500 MΩ were inserted.
Electrical stimulation (rectangular current pulses of 0.50 nA and 0.50 ms duration) was applied to the inside of either the first cell of chain A (cell A1) or simultaneously to all cells of the A-chain. For retrograde propagation, stimulation was applied either to cell G1 or to all cells of the G-chain. For some measurements, the V-recording markers were placed on only one chain at a time. To minimize confusion, the voltage was recorded from only one surface unit (upward-facing) in each cell.
The results to be illustrated here will be from the 7 × 7 model only, because this model is new. However, the results from the smaller models (5 × 5, 3 × 4, 2 × 3), previously published, are summarized in Table 1. Thus, Table 1 enables transverse propagation velocities to be compared in four models differing in size but with identical circuit parameters used in the basic units. Through this comparison, it can be ascertained whether edge (boundary) effects are important in transverse propagation.
Table 1 summarizes all these data, not only for the 7 × 7 model, but also for the smaller models of 5 × 5, 3 × 4, and 2 × 3. These data include antegrade (A) and retrograde (R) propagation, with stimulation of the entire chain or single cell only, for the Rol2 value of 500 KΩ. As can be seen, the calculated transverse propagation velocities were highest in the large 7 × 7 model and slower in the smaller models. This was true for all values of Rol2. Transverse propagation velocities were faster at Rol2 of 500 KΩ than at 200 KΩ.
Comparison of the inverse ratios of the relative edge areas of the various-sized cardiac models with the ratio of the transverse propagation velocities (θtr)
Ratio of relative edge area to interior area
7 × 7
28/49 = 0.57
5 × 5
20/25 = 0.80
3 × 4
14/12 = 1.17
2 × 3
10/6 = 1.67
This means that one can predict the transverse propagation velocities in yet-larger models. For example, in a 10 × 10 model, θtr should be approximately 6.70 cm/s (if compared with the 7 × 7 model) or 6.80 cm/s (if compared with the 5 × 5 model). From the equation given in the footnote of Table 2:
Hence, the two calculations are in close agreement. However, this relationship between transverse velocity and the inverse of the relative edge area probably saturates and levels off at some point, i.e., a maximum θtr is reached. In the intact heart, the velocity of transverse propagation is difficult to measure accurately because of the complicated geometry of bundles, but estimates that θtr is about 1/5th to 1/10th that of θlo (longitudinal velocity) have been given (see references given in ref 1). If θlo is taken to be 0.40 m/s, then θtrshould be between 4.0 and 8.0 cm/s. Thus, the values calculated in the present simulations are in good agreement with physiological measurements.
The ratio of propagation velocities, longitudinal (θlo) to transverse (θtr), is almost what is expected based on the cell geometry (cylinder 150 μm long and 16 μm wide). These dimensions would predict a θlo / θtr ratio of 9.4 (150/16), provided that the longitudinal and transverse transfer function are equal (i.e, the delay time at the two types of junctions were equal). If the average cell length were taken to be only 100 μm, then the θlo / θtr ratio would be 6.3. Thus, for a θlo value of 40 cm/s and a θlo / θtr ratio of 7.9 (average of 9.4 and 6.3), then θtr would be 5.1 cm/s, which is close to the value of 4.7 cm/s measured in the 7 × 7 model (for Rol2 of 200 KΩ). However, it has been reported that the anisotropic conduction velocity observed in the heart is not a result of cell geometry .
Another observation in the large 7 × 7 model is that some chains distal to the point of stimulation failed when Rol2 was only 200 KΩ. Such failures did not occur when the model was smaller (e.g., 5 × 5). Failure of distal chains occurred in both the orthodromic and antidromic directions, but was greater in the antidromic direction. However, increasing Rol2 to 500 KΩ allowed all chains to respond, with the exception of failure of one chain (the most distal A-chain) in the retrograde direction (Table 1). Therefore, in the largest model, there is an increase in probability of failure of one or more distal chains.
Although we don't know the mechanism for this effect, we may speculate about two possibilities. First, if some current leaked out at the ends of each chain, then less current would be available for downstream depolarization. Second, if the phenomenon of reflection occurred at the longitudinal edge of the last chain (G), then this would act to slow the transverse velocity. Thus, both of these mechanisms may be involved in explaining why transverse propagation was faster in the larger models.
In summary, the present results using our enlarged 7 × 7 model for cardiac muscle, with comparisons with our prior smaller models, demonstrate that edge effects are important in determining the transverse velocity of propagation, when all circuit parameters are identical. θtr increased with the inverse of the ratio of the relative edge areas in the various-sized models. This relationship likely levels off at some point, such that a maximum velocity is reached. The transverse velocities measured in the largest model (7 × 7), and estimated for a 10 × 10 model, give values in the same range as the physiological values.
- Sperelakis N, McConnell K: Electric field interactions between closely abutting excitable cells. IEEE-Eng Med Biol. 2002, 21: 77-89. 10.1109/51.993199.View ArticleGoogle Scholar
- Sperelakis N, Ramasamy L: Propagation in cardiac muscle and smooth muscle based on electric field transmission at cell junctions: An analysis by PSpice. IEEE-Eng Med Biol. 2002, 21: 130-143. 10.1109/MEMB.2002.1175149.View ArticleGoogle Scholar
- Sperelakis N, Murali KPV: Combined electric field and gap junctions on propagation of action potentials in cardiac muscle and smooth muscle in PSpice simulation. J Electrocardiol. 2003, 36 (4): 279-293. 10.1016/j.jelectrocard.2003.08.001.View ArticlePubMedGoogle Scholar
- Sperelakis N: Propagation of action potentials between parallel chains of cardiac muscle cells in PSpice simulation. Can J Physiol Pharmacol. 2003, 81: 1-112. 10.1139/y03-019.View ArticleGoogle Scholar
- Sperelakis N, Kalloor B: Transverse propagation of action potentials between parallel chains of cardiac muscle and smooth muscle cells in PSpice simulations. Biomed Eng Online. 2004, 3: 5-10.1186/1475-925X-3-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Sperelakis N, Kalloor B: Effect of variation in membrane excitability on propagation velocity of simulated action potentials for cardiac muscle and smooth muscle in the electric field model for cell to cell transmission of excitation. IEEE-Eng Med Biol. 2004, 51: 2216-2219. 10.1109/TBME.2004.836528.View ArticleGoogle Scholar
- Thomas SP, Kucera JP, Bircher-Lehmann L, Rudy Y, Saffitz JE, Kleber AG: Impulse propagation in synthetic strands of neonatal cardiac myocytes with genetically reduced levels of connexin43. Circ Res. 2003, 91: 1209-1216. 10.1161/01.RES.0000074916.41221.EA.View ArticleGoogle Scholar
- Darrow BJ, Fast VG, Kleber AG, Beyer EC, Saffitz JE: Functional and structural assessment of intercellular communication. Increased conduction velocity and enhanced connexin expression in dibutyryl cAMP-treated cultured cardiac myocytes. Circ Res. 1996, 79: 174-183.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.