Reach Us
+441414719275

^{1}Oakland University, Rochester, MI, USA

^{2}University of Alberta, Edmonton, Alberta, Canada

- *Corresponding Author:
- Piyush Kar

University of Alberta

Edmonton, Alberta, Canada

**Tel:**514-7122130

**E-mail:**[email protected]

**Received date:** November 04, 2015; **Accepted date:** November 25, 2015; **Published date:** November 30, 2015

**Citation: **Dhungana S, Kar P (2015) Simulation of Propagation of Action Potentials in Cardiac Tissue with an Inhomogeneous Distribution of Extracellular Potassium. J Comput Sci Syst Biol 8:6 373-379. doi:10.4172/jcsb.1000212

**Copyright:** © 2015 Dhungana S, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Visit for more related articles at** Journal of Computer Science & Systems Biology

We report simulation of experimental data for heterogeneous distribution of extracellular potassium, by developing and applying a mathematical model formulated by incorporating a Luo-Rudy representation of the cell membrane current into the bidomain model of the cardiac tissue. Our model serves to examine propagation of extra cellular potassium from a region of low potassium (normal) to a region of high potassium (ischemic). Starting with computer simulation of unipolar stimulation of cardiac tissue, we reproduced results of the passive bidomain model, and also simulated the space-clamped Luo Rudy model. Finally, we combined the Luo-Rudy model with the bidomain model. The diastolic threshold for anodal stimulation first decreased and then increased with increasing potassium concentration ([K]0), reaching a minimum value at [K]0=12mM.

Mathematical modeling; Luo-Rudy model; Bidomain model; Extracellular potassium; Cardiac tissue; Action potential

Our mind works and our heart beats as long as there is propagation of action potentials in our body. Action potentials are generated by the inward and outward movements of numerous ions, such as Na^{+}, K^{+}, Ca^{2+}, Cl^{-} across the cell membrane. These action potentials travel in the form of waves through the various chambers of the heart. To fully understand the heart, a detailed knowledge of action potential propagation is essential.

In unhealthy (ischemic) tissue, the extracellular potassium concentration, [K]_{0}, is elevated. High [K]_{0} is believed to predispose cardiac tissue to break excitation (Roth and Patel) [1], which in turn plays an important role in re-entry induction and defibrillation. Because fibrillation raises extracellular [K]_{0} levels, break excitation may play a more important role in defibrillation than is suggested by simulations and experiments using normal [K]_{0} values. In their 2003 paper, Roth and Patel theoretically simulated action potential data in a 1 cm^{2} tissue verifying that the action potential may fail to propagate in an ischemic region. The wave propagates for [K]_{0}=2 to 12.5mM and at higher potassium concentrations the wave dies out.

Sidorov et al. studied propagation in ischemic tissue experimentally. Their optical mapping experiment in a rabbit heart using [K]_{0} values of 4, 6, 8, 10 and 12 mM indicated that spatial heterogeneity of the action potential duration (APD) restitution, created with regional elevation of [K^{+}]_{0}, can lead to action potential (AP) instability, 2:1 block, and reentry induction [2].

The goal of this research is to simulate Sidorov et al. experiment using a mathematical model by incorporating a Luo-Rudy representation of the cell membrane current with the bidomain model of the cardiac tissue, and to examine propagation of potassium from a normal to an ischemic region.

**Luo rudy model**

Hodgkin and Huxley published a landmark paper that paved the way for subsequent electrophysiological models. They successfully explained the propagation of action potentials in the giant squid axon by representing the membrane resistance as resulting from ion flow through three ion channels: sodium, potassium and leakage [3]. Each pathway was accompanied by a battery that had a potential equal to the Nernst potential for that ion. Variable resistances signifying the various ion channels opening and closing were expressed in terms of their conductance. The probability of opening or closing of these ionic channels was governed by three gating parameters m, n and h following the first order kinematics.

(1)

In eq. (1), αi and βi represent the rate of opening and closing of the m, n, and h gates. These rate constants depended on the transmembrane voltage.

Hodgkin and Huxley wrote down the rate of change of membrane potential as

(2)

Following the mathematical formulation of Hodgkin and Huxley, in 1991 Luo and Rudy developed a model of the ventricular cardiac action potential. They introduced and analyzed more ion channels then other contemporary physiological models. They replaced the three ionic current terms in Hodgkin and Huxley’s model with six other currents. They were a fast inward sodium Na^{+} current I_{Na}, a slow inward calcium Ca^{++} current I_{Si}, a time dependent delayed rectifier potassium K^{+} current I_{k}, and time-independent potassium inwardly rectifying current I_{k1}, the plateau K^{+} current I_{kp}, and a time independent background current I_{b}. The membrane current depended on [K^{+}]_{0} as a parameter in this model. Considering a membrane patch of area 1 cm^{2} with capacitance 1 μF/cm^{2} equation (2) can be rewritten as

(3)

All these ionic currents also depend on the gating variables, which are actually probabilistic values indicating the closing and opening of various ion channels. The time variation is given as

(4)

where,

(5)

and,

(6)

Here, y can be any of the gating variables (m,h,….), τ_{y} is the time constant, and y_{∞} is the steady state value of y. α_{y} and β_{y} are voltage dependent rate constants.

The fast inward sodium current in eq. (3) was described using the Ebihara-Johnson model [4] (Ebihara and Johnson) coupled with the slow inactivation gate (j) of the Beeler-Reuter model (Beeler and Reuter) [5],

(7)

where, is the ,…maximum conductance of the sodium channel (23 mS/cm^{2}); m the activation parameter and h the inactivation parameter; and E_{Na} is the reversal potential of sodium.

(8)

(9)

where,

(10)

and, d and f are two other gating parameters.

The delayed outward rectifier current I_{k} was based on a formulation proposed by Shibasaki and included a time–dependent activation gate X, a time independent (fast) inactivation gate X_{i} capturing the inward rectification properties of the delayed rectifier current, and a conductance modulated by extracellular K^{+} concentration [6].

(11)

where,

(12)

The formulation of Ik_{1} was based on patch clamp data of Sakmann and Trube as well as other data from single ventricular myocytes [7].

(13)

where,

(14)

The voltage-dependent and time independent plateau current that activates at depolarized plateau potentials was identified by Yue and Marban [8],

(15)

Where, E_{kp}= E_{k1}

(16)

The background current is formulated as

(17)

where, ,

Then the total time independent current is

(18)

Finally the rate of calcium uptake is given by

(19)

complete listing of all the governing equations, parameters and their corresponding values can be found in Appendix.

**Bidomain model**

The two or three dimensional cable model, more commonly known as the bidomain model, describes the interaction of electric current with cardiac tissue. It was developed in the late 1970’s (Muler and Markin, Miller and Geselowitz, Tung) [9-11] and was applied to more problems in cardiac electrophysiology in the 1980s (Plonsey and Barr, Barr and Plonsey, Roth and Wikswo) [12-14].

Here we will solve the coupled parabolic and elliptic equations using a finite difference technique (Roth) and discuss some consequences of unipolar stimulation of cardiac tissue [15]. The first analysis of this kind of work was done by Sepulveda et al. [16]. The starting point here is to solve the set of two coupled partial differential equations governing the transmembrane potential, V_{m}, and extracellular potential, V_{e}

(20)

(21)

To solve the parabolic type of equation, we use the simple Euler method. The idea here is to solve eq. (20) for the new value of V_{m}, and then solve the elliptic eq. (21) for the new value of V_{e} using the new value for V_{m} as a source term, by over-relaxation.

By, rearranging, we get,

(22)

In eq. (22), let,

So, we can write an expression for extracellular potential, V_{e} in terms of V_{m}, as follows:

By developing a code using the FORTRAN software, we simulated the distribution of extracellular potassium by using the bidomain model with a LR action potential. Sepulveda et al. used a passive bidomain model to study unipolar stimulation of cardiac tissue [16]. They found that near a unipolar cathode, the tissue is strongly depolarized under the cathode, and more weakly hyperpolarized in two regions near the cathode along the fiber direction. **Figure 1** shows our model generated transmembrane potential produced around the cathode, located at the origin.

In presence of the cathode the tissue is depolarized (positive transmembrane potential) in a dogbone shape. However, about 2 mm from the cathode, in the direction along the fibers, is a region of hyperpolarization (negative), called a “virtual anode”. In Figure 1, only one quarter of the tissue is shown, so two virtual anodes exist one on each side of the cathode. The values in the plots are given in volts. **Figure 2** shows anodic stimulation, which gives rise to two virtual cathodes on either side of the origin. **Figure 3** shows the calculated transmembrane potential during unipolar cathodal stimulation of an anisotropic twodimensional sheet of passive cardiac tissue in a color plot.

A space clamped calculation of the LR action potential was performed in Luo and Rudy’s original paper (Luo and Rudy) [17]. **Figure 4** shows simulated action potentials for different extracellular concentrations of potassium.

As [K]_{0} decreases, action potential duration increases and resting potential becomes more negative. **Figure 5** shows the time course of I_{K1(T)}, I_{K}, and I_{si} during the action potential. A peak in I_{K1(T)} is observed during repolarization. For [K]_{0}=7mM the peak is almost 3 μA/cm^{2} but for [K]_{0}=3mM it is less than 2 μA/cm^{2}. This large difference in I_{K1(T)} brings about a faster rate of repolarization at higher [K]_{0}. Also the peak of I_{K1(T)} occurs earlier (at more positive membrane potential), reflecting an increase of the reversal potential E_{K1}.

Furthermore, **Figures 5-7** indicates that during repolarization the total time independent potassium current depends strongly on [K]_{0}, while the time–dependent potassium current (I_{K}) and slow inward current exhibit only weak dependence on [K]_{0}.

The 1991 Luo-Rudy Phase I model did not describe intracellular Ca^{2+} cycling, sarcolemmal Na^{+}-Ca^{2+} exchange, Na^{+}-K^{+} pumping and temporal variation of intracellular ion concentrations [17]. Luo and Rudy produced another model called Luo-Rudy phase II model, which is one of the most popular cardiac ventricular cell models even today (Luo and Rudy) [18]. However, because focus is made on extracellular potassium concentration and not on intracellular calcium dynamics, it is not relevant here.

Our next analysis is made by combining the bidomain model with an active Luo-Rudy membrane, which bears similarity with simulations performed previously by Roth [15]. **Figure 8** shows stimulation threshold as a function of [K]_{0}.. The threshold initially decreases as [K]_{0} increases, but then rises abruptly for [K]_{0} above 12mM. **Figure 9** shows the resting potential as a function of potassium ion concentration. The resting potential follows approximately the potassium Nernst potential, with a very negative value at low [K]_{0} and then rising to less negative values as [K]_{0} increases.

Having verified that the model gives results similar to previous calculations, we can now apply the model to simulating the experiments of Sidorov et al. [2]. We can consider a tissue that has [K]_{0}=4mM (normal) on the left half of the tissue, and [K]_{0}=10mM (ischemic) on the right (**Figure 10**). The tissue sheet is 10 mm by 10 mm (100 × 100 nodes), with a space step in each direction of 0.1 mm.

We can stimulate the tissue by raising the extracellular potential along the left edge to 20 mV at time t=0. An action potential then propagates across the tissue from left to right. Plots of the transmembrane potential versus time at three locations are shown in **Figures 11 and 12**. Point A (x=-4.92 mm, y=-4.92 mm) is near the left edge of the tissue sheet in the region of normal [K]_{0}, point B (x=-0.025, y=-0.025) is at the center (just on the boundary between normal and ischemic tissue), and point C (x=4.9 mm, y=4.9 mm) is near the right edge in the region of elevated [K]_{0}.

The second stimulus was applied at various times, but never found a case where propagation was successful in the normal tissue but then failed to enter the ischemic tissue. Also attempt was made by applying more than two stimuli (up to five) with various timings, and varied the distribution of extracellular potassium so that in one case the left was 5.4mM and the right was 8mM (**Figure 13**), and in another case the left was 8mM and the right was 12mM (**Figure 14**). Nevertheless, after extensive simulations none of these simulations resulted in propagation on the left but then failure on the right.

We validated our calculations by combining both models by verifying previous calculations, in that we reproduced results reported by Sepulveda et al. [16], Luo-Rudy action potential (Luo and Rudy, 1991) and by Roth and Chen [17,19]. This exercise brought about the conclusion that resting potential rises with [K]_{0}, which tends to make propagation easier (lower threshold) initially as [K]_{0} increases, but makes propagation harder (higher threshold) for very high [K]_{0} (Figures 8 and 9).

Finally the Sidorov et al. experiment was simulated. In their case, the wave front propagated from a region of [K]_{0}=4mM to a region of [K]_{0}=10mM in a Langendorff-perfused rabbit heart [2]. They often found cases in which propagation was successful in the normal (4mM) tissue but then failed at the boundary between the normal and ischemic (10mM) tissue. Here this behavior was never observed. The reason appears to be the dependence of action potential duration on [K]_{0}. In the Luo-Rudy model, the action potential duration decreases as [K]_{0} increases, so that the refractory period in the ischemic region is short. Any action potential that can propagate through the normal tissue (with a long refractory period) should have no problem also propagating through the ischemic tissue (with a short refractory period). Interestingly, in Sidorov et al. data, the action potential duration (and therefore the refractory period) appeared insensitive to extracellular potassium ion concentration. This appears to be the primary difference between our simulations (reported herein) and their experiments. Future simulations of Sidorov et al.’s data would need to include a model of the membrane kinetics that properly captures the [K]_{0} dependence of the action potential duration.

We simulated Sidorov’s experiment, in rabit heart, using a mathematical model, wherein we incorporated Luo Rudy’s representation of the cell membrane current with a bidomain model of the cardiac tissue.

In addition to matching experimental results by Sidorov, reliability of our model is also demonstrated by reproduction results of previously reported models that predicted dependence of extracellular potassium on action potential duration.

Using our model, we examined the propagation of potassium from low to high concentration and found that [K]_{0} is dependent on the action potential duration. Further research is necessary to refine our modelling approach by making it more inclusive of physiological parameters.

We thank Professor Bradley Roth (Oakland University) for his important suggestion and careful editing of the manuscript.

- Roth BJ, Patel SG (2003) Effects of elevated extracellular potassium ion concentration on anodal excitation of cardiac tissue.J CardiovascElectrophysiol 14: 1351-1355.
- Sidorov VY, UzelacI, Wikswo JP (2011) Regional increase of extracellular potassium leads to electrical instability and reentry occurrence through the spatial heterogeneity of APD restitution. Am J Physoiol Heart CircPhysiol 301:H209-H220.
- HODGKIN AL, HUXLEY AF (1952) A quantitative description of membrane current and its application to conduction and excitation in nerve.J Physiol 117: 500-544.
- Ebihara L, Johnson EA (1980) Fast sodium current in cardiac muscle. A quantitative description.Biophys J 32: 779-790.
- Beeler GW, Reuter H (1977) Reconstruction of the action potential of ventricular myocardial fibres.J Physiol 268: 177-210.
- Shibasaki T(1987) Conductance and kinetics of delayed rectifier potassium channels in nodal cells of the rabbit heart.J Physiol 387: 227-250.
- Sakmann B, Trube G (1984) Conductance properties of single inwardly rectifying potassium channels in ventricular cells from guinea-pig heart.J Physiol 347: 641-657.
- Yue DT, Marban E (1988) A novel cardiac potassium channel that is active and conductive at depolarized potentials.Pflugers Arch 413: 127-133.
- Muler AL, Markin VS (1977) [Electrical properties of anisotropic neuromuscular syncytia. I. Distribution of the electrotonic potential].Biofizika 22: 307-312.
- Miller WT, Geselowitz DB (1978) Simulation studies of the electrocardiogram. I. The normal heart.Circ Res 43: 301-315.
- Tung L (1978) Abidomain model for describing ischemic myocardial dc potentials. PhD Dissertation, Massachusetts Institute of Technology, Cambridge, MA.
- Plonsey R, Barr RC (1984) Current flow patterns in two-dimensional anisotropic bisyncytia with normal and extreme conductivities.Biophys J 45: 557-571.
- Barr RC, Plonsey R (1984) Propagation of excitation in idealized anisotropic two-dimensional tissue.Biophys J 45: 1191-1202.
- Roth BJ, Wikswo JP Jr (1986) Abidomain model for the extracellular potential and magnetic field of cardiac tissue.IEEE Trans Biomed Eng 33: 467-469.
- Roth BJ (1991) Action potential propagation in a thick strand of cardiac muscle.Circ Res 68: 162-173.
- Sepulveda NG, Roth BJ, Wikswo JP Jr (1989) Current injection into a two-dimensional anisotropic bidomain.Biophys J 55: 987-999.
- Luo CH, Rudy Y (1991) A model of the ventricular cardiac action potential. Depolarization, repolarization, and their interaction.Circ Res 68: 1501-1526.
- Luo CH, Rudy Y (1994) A dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes.Circ Res 74: 1071-1096.
- Roth, BJ, Chen J (1999) Mechanism of anode break excitation in the heart: The relative influence of membrane and electrotonic factors. J Biol Systems7:541-552.

**Appendix**

Governing equations, parameters and their corresponding values for the Luo Rudy Model

Slow inward Current:

Calcium uptake: d([Ca]_{i})/dt=-10^{-4}.I_{si}+0.07(10^{-4}-[Ca]_{i})

Expressions for outward currents are given below.

Time dependent Potassium Current:

Time- independent potassium current:

Plateau Potassium current:

Background current:

Total time independent potassium current:

Select your language of interest to view the total content in your interested language

- Algorithm
- Artificial Intelligence Studies
- Big data
- Bioinformatics Modeling
- Biostatistics: Current Trends
- Cloud Computation
- Computational Sciences
- Computer Science
- Computer-aided design (CAD)
- Data Mining Current Research
- Findings on Machine Learning
- Mathematical Modeling
- Neural Network
- Robotics Research
- Soft Computing
- Studies on Computational Biology
- Systems Biology

- Total views:
**9181** - [From(publication date):

November-2015 - Jul 17, 2019] - Breakdown by view type
- HTML page views :
**9042** - PDF downloads :
**139**

**Make the best use of Scientific Research and information from our 700 + peer reviewed, Open Access Journals**

International Conferences 2019-20