CMMM Computational and Mathematical Methods in Medicine 1748-6718 1748-670X Hindawi Publishing Corporation 10.1155/2014/761907 761907 Research Article Effects of Maximal Sodium and Potassium Conductance on the Stability of Hodgkin-Huxley Model http://orcid.org/0000-0002-4515-8931 Zhang Yue 1 Wang Kuanquan 1 Yuan Yongfeng 1 http://orcid.org/0000-0002-7887-2111 Sui Dong 1 Zhang Henggui 1, 2 Jiang Rui 1 Biocomputing Research Center School of Computer Science and Technology Harbin Institute of Technology Harbin 150001 China hit.edu.cn 2 School of Physics & Astronomy University of Manchester Manchester UK manchester.ac.uk 2014 372014 2014 13 02 2014 12 06 2014 15 06 2014 3 7 2014 2014 Copyright © 2014 Yue Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Hodgkin-Huxley (HH) equation is the first cell computing model in the world and pioneered the use of model to study electrophysiological problems. The model consists of four differential equations which are based on the experimental data of ion channels. Maximal conductance is an important characteristic of different channels. In this study, mathematical method is used to investigate the importance of maximal sodium conductance g - N a and maximal potassium conductance g - K . Applying stability theory, and taking g - N a and g - K as variables, we analyze the stability and bifurcations of the model. Bifurcations are found when the variables change, and bifurcation points and boundary are also calculated. There is only one bifurcation point when g - N a is the variable, while there are two points when g - K is variable. The ( g - N a ,   g - K ) plane is partitioned into two regions and the upper bifurcation boundary is similar to a line when both g - N a and g - K are variables. Numerical simulations illustrate the validity of the analysis. The results obtained could be helpful in studying relevant diseases caused by maximal conductance anomaly.

1. Introduction

Hodgkin-Huxley (HH) equation is created on the foundation of huge experimental data of sodium and potassium channels by Hodgkin and Huxley who are both excellent biology scientists and had long engaged in nerve conduction research. In about 1952, they took squid giant axon as experiment subject and continuously published four papers describing the electrical excitation of this kind of cell . In their experiment, all the ion channels were divided into three types, sodium channel, potassium channel, and the others. Now we know there are many ion channels on the cell membrane, such as I Na , I Kr , I Ks , I NaCa , I K 1 , I CaL , I Ca , I to , I NaK , I NaL , and I KATP . However the discovery of sodium and potassium channels was marvelous at that time. Experimental data was obtained by voltage-clamp technology, while the patch-clamp technology is widely used at present. On this basis, a four-dimensional ordinary differential equation set, called HH model, was proposed, which was autonomous and contained intricate transcendental equations.

The work of Hodgkin and Huxley was recognized as excellent achievement and with significant contribution to the development of electrophysiology. It is the basis of the subsequent models of ion channels. Not only was the HH model consistent with the obtained experiment data accurately, but also it could precisely simulate the change of action potential. The model discovered the relationship of transmembrane potential and current and maximum conductance of ions. This made it possible to research the character of ion channel with mathematical methods. In 1960, Professor Nobel who pioneered the cardiac electrophysiology simulations applied HH model to myocardial cell and got the famous Purkinje fiber cell model , which was the first computing myocardial cell model. From then on, HH model was broadly applied to almost all kinds of cardiac cells such as atrial muscle cell model  and sinoatrial node cell model . HH model laid the cornerstone of computing electrophysiology. Even today, a large part of electrophysiological models are created on the foundation of HH model. Verkerk’s sinoatrial model , Butters’s atrial model , O’Hara’s ventricular model , and Li’s Purkinje cell model , and so forth, all belong to HH model type.

Because of the importance of HH model, the stability has long attracted the researcher’s attention. Hassard et al. were the earlier researchers caring about the bifurcation phenomenon of HH model. And they indicated that bifurcation would occur at the equilibrium points when the external current i ext changed which was injected into the neuron from microelectrode . Stable and unstable solutions of the model with regard to i ext were analyzed by Rinzel and Miller, and the influence of temperature was also discussed . Two stable steady states were found by Aihara and Matsumoto ; when the two states existed, the bifurcation structure was complex, which included a stable limit cycle, two unstable equilibrium points, and one asymptotically stable equilibrium point. Guckenheimer and Labouriau investigated the influence of i ext and potassium ion potential V K on the bifurcations of the model . Bedrov et al. gave the relationship between the numbers of negative slope regions and presented some results about the possible bifurcation giving rise to maximal sodium conductance g - Na and maximal potassium conductance g - K [19, 20]. Fukai and his fellows examined the structure of the model’s bifurcations produced by i ext and one of the other parameters . Taking leakage conductance g l and sodium channel effective bias voltage V - m as parameters, Terada et al. analyzed bifurcation in Hodgkin-Huxley model for muscles of frogs .

Wang et al.  did a lot of research on the stability of HH model. They studied the bifurcations caused by leakage conductance g - l and sodium ions antielectromotive force when ELF external electric field was considered. The stability and bifurcation control were analyzed and controllers were designed. Bifurcation in HH model exposed to DC electric fields was investigated in detail.

Bifurcation means qualitative changes in the solution structure of a dynamic system when the parameters vary. From analyzing the bifurcation, we can get the effects of the parameters. Further, changing the corresponding parameters, we could make the solution into an ideal condition. Bifurcation is an important branch in mathematics and applied to much different field . In recent years, it is also widely studied in electrophysiology. Indeed, there are many diseases having close relations with bifurcations, such as Parkinson’s, epilepsy, and pathological heart rhythms .

In the past, for HH model, external current i ext and leakage conductance g - l have been most investigated, because they were easily measured. The sodium current is the contributor which leads to depolarization of the neuron while it is potassium current that plays the major role of repolarization. However, g - Na and g - K are seldom taken into consideration to analyze the stability of model, as the relevant data is not abundant. In this study, the effects of   g - Na and g - K on the stability and bifurcations of the model will be discussed, respectively, and collectively. And we will give the critical points of g - Na and g - K when they play the role separately, and the critical boundaries in g - Na - g - K plane will be provided when together. Simulation results demonstrate the validity of the theoretic analysis.

The rest of the paper is organized as follows. The HH equations are introduced in detail in Section 2. In Section 3 we analyze the effects of g - Na and g - K on the model and calculate the bifurcation points (line). Finally, discussion and conclusion are presented in Section 4.

2. Hodgkin-Huxley Equations

HH model was proposed on the foundation of ion channels. The electrophysiological activities of a cell are shown in Figure 1(a). The gray circle is membrane, which ensures orderly biochemical reaction. I Na , I K , and I L are the ion currents corresponding to respective channels on the membrane. When an electrical stimulation makes the sodium channels open, a large number of Na+ flow inward, forming current I Na , resulting in the rise of transmembrane potential. The open of potassium channels makes a large outflux of K+, creating the current I K and the reduction of potential. The model is comprised of four autonomous ordinary differential equations to describe the electrophysiological activities of cell shown in Figure 1(a). In the model, membrane is taken as a constant capacitance and the ion channels are seen as variable resistances. Figure 1(b) shows the equivalent circuit in detail, in which R Na = 1 / g Na , R K = 1 / g K , and R l = 1 / g l . R Na and R K vary with time.

The electrophysiological process and equivalent circuit of neuron.

Electrophysiological activities of cell

Electrical equivalent circuit of HH model

The equations were obtained according to electrical formulas and experimental data, which are shown as follows: (1) d V d t = 1 C M [ i ext - g - Na m 3 h ( V - V Na ) hhhhhh - g - K n 4 ( V - V K ) - g - l ( V - V l ) ] , d m d t = α m ( V ) ( 1 - m ) - β m ( V ) m , d h d t = α h ( V ) ( 1 - h ) - β h ( V ) h , d n d t = α n ( V ) ( 1 - n ) - β n ( V ) n , where (2) α m ( V ) = 0.1 ( V - 25.0 ) 1 - exp [ - ( V - 25.0 ) / 10 ] , β m ( V ) = 4.0 exp ( - V 18.0 ) , α h ( V ) = 0.07 exp ( - V 20.0 ) , β h ( V ) = 1 1 + exp [ - ( V - 30.0 ) / 10 ] , α n ( V ) = 0.01 ( V - 10.0 ) 1 - exp [ - ( V - 10.0 ) / 10 ] , β n ( V ) = 0.125 exp ( - V 80.0 ) .

In these equations, V is the transmembrane potential. 0 m 1 and 0 h 1 are the gating variables indicating activation and inactivation of sodium ion current, respectively. 0 n 1 is the gating variable showing activation of potassium ion current. g - Na , g - K , and g - l represent the maximal conductance of corresponding currents. C m = 1.0 μF/cm2 is membrane capacitance. i ext is the current injected into the neuron. In our paper, we suppose i ext = 0 and g - Na = 120  mS/cm2, g - K = 36  mS/cm2, and g - l = 0.3  mS/cm2, which are the ideal experimental data.

3. Stability Analysis of HH Model

Stability is one of a model’s important properties. If the model is stable, it will reach a rest state at last. Otherwise, periodic phenomenon or chaos may appear. To analyze an ordinary differential system, equilibrium points are one of its most important aspects, which may be the final state of the system. Suppose ( V * , m * , h * , n * ) is the equilibrium points of the model. So it should make the right side of (1) equal to zero. That is, (3) i ext - g - Na m * 3 h * ( V * - V Na ) - g - K n * 4 ( V * - V K ) - g - l ( V * - V l ) = 0 , α m ( V * ) ( 1 - m * ) - β m ( V * ) m * = 0 , α h ( V * ) ( 1 - h * ) - β h ( V * ) h * = 0 , α n ( V * ) ( 1 - n * ) - β n ( V * ) n * = 0 . Then the linearization of (1) around the equilibrium could be obtained as follows: (4) d V d t = J 11 V + J 12 m + J 13 h + J 14 n , d m d t = J 21 V + J 22 m , d h d t = J 31 V + J 33 h , d n d t = J 41 V + J 44 n , where (5) J 11 = - g - Na m * 3 h * + g - K n * 4 + g - l C M , J 12 = - 3 g - Na m * 2 h * ( V * - V Na ) C M , J 13 = - g - Na m * 3 ( V * - V Na ) C M , J 14 = - 4 g - K n * 3 ( V * - V K ) C M , J 21 = 2 m * 9 exp ( V * / 18.0 ) - { 0.1 ( V * - 25.0 ) exp [ - ( V * - 25.0 ) / 10 ] [ 1 - exp ( - ( V * - 25 ) / 10 ) ] 2 hhhhh - 0.1 1 - exp ( - ( V * - 25 ) / 10 ) 0.1 ( V * - 25.0 ) exp [ - ( ( V * - 25.0 ) / 10 ) ] [ 1 - exp ( - ( ( V * - 25 ) / 10 ) ) ] 2 } ( 1 - m * ) , J 22 = 0.1 ( V * - 25.0 ) 1 - exp [ - ( V * - 25.0 ) / 10 ] - 4.0 exp ( - V * 18.0 ) , J 31 = - 7 exp ( V * / 20.0 ) 2000 ( 1 - h * ) - exp ( - ( V * - 30.0 ) / 10 ) h * 10 [ 1 + exp ( - ( V * - 30.0 ) / 10 ) ] 2 , J 33 = 0.07 exp ( - V * 20.0 ) - 1 1 + exp ( - ( V * - 30.0 ) / 10 ) , J 41 = n * exp ( V * / 80.0 ) 640 + { 0.01 exp [ - ( ( V * - 10.0 ) / 10 ) ] ( V * - 10.0 ) 10 [ 1 - exp ( - ( ( V * - 10.0 ) / 10 ) ) ] 2 1 100 [ 1 - exp ( - ( V * - 10.0 ) / 10 ) ] - 0.01 exp [ - ( V * - 10.0 ) / 10 ] ( V * - 10.0 ) 10 [ 1 - exp ( - ( V * - 10.0 ) / 10 ) ] 2 } ( 1 - n * ) , J 44 = 0.01 ( V * - 10.0 ) 1 - exp [ - ( V * - 10.0 ) / 10 ] - 0.125 exp ( - V * 80.0 ) .

We can get the eigenmatrix of (4): (6) J = ( J 11 J 12 J 13 J 14 J 21 J 22 0 0 J 31 0 J 33 0 J 41 0 0 J 44 ) and then the characteristic equation can be obtained: (7) λ 4 + a λ 3 + b λ 2 + c λ + d = 0 , where (8) a = - ( J 11 + J 22 + J 33 + J 44 ) , b = J 11 ( J 22 + J 33 + J 44 ) + J 22 ( J 33 + J 44 ) + J 33 J 44 - J 12 J 21 - J 13 J 31 - J 14 J 41 , c = J 12 J 21 ( J 33 + J 44 ) + J 13 J 31 ( J 22 + J 44 ) + J 14 J 41 ( J 22 + J 33 ) - J 11 J 22 ( J 33 + J 44 ) - ( J 11 + J 22 ) J 33 J 44 , d = J 11 J 22 J 33 J 44 - J 12 J 21 J 33 J 44 - J 13 J 22 J 31 J 44 - J 14 J 22 J 33 J 41 .

According to Routh-Hurwitz criterion, if a > 0 , a b > c , d > 0 ,   a b c > c 2 + a 2 d , the real parts of all the roots are minus. Otherwise, all the real parts are not negative. In the following, we will analyze the stability of the model according to this criterion.

3.1. Effect of <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M82"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mover accent="true"> <mml:mrow> <mml:mi>g</mml:mi></mml:mrow> <mml:mo>-</mml:mo></mml:mover> </mml:mrow> <mml:mrow> <mml:mtext>Na</mml:mtext></mml:mrow> </mml:msub></mml:mrow> </mml:math></inline-formula> on the Stability

In this section, we will investigate the influence of g - Na on the equilibrium, stability, and bifurcations of the model. g - Na is taken as variable, and the other parameters are all kept with desired values. Because the desired value of g - Na is around 120 mS/cm2, g - Na [ 0,500 ] is taken into consideration. When g - Na changes, making the right side of (1) equal to zero, corresponding V * can be acquired. Taking Matlab as a tool, we could obtain the relationship between g - Na and the equilibrium point V * shown in Figure 2.

The relationship between g - Na and V * .

From Figure 2, we can see that V * changes slowly when g - Na [ 0,300 ] and increases rapidly when g - Na [ 350,500 ] . This means that equilibrium points are sensitive to g - Na when g - Na [ 350,500 ] ; a slight change of g - Na may lead the model to a totally different state even though model is still stable.

Applying bifurcation theory and using the method of bisection, we can get one bifurcation point g - Na * = 212.648720656 when g - Na changes.

Substituting g - Na * into the original equation, we get the equilibrium V * , and then substituting both g - Na * and V * into eigenmatrix of (4), we can gain the eigenvalues as follows: (9) λ 1 = - 4.9711711484 , λ 2 = - 0.1259717048 , λ 3 =    1.9 × 10 - 16 - 0.3798402483 i , λ 4 = 1.9 × 10 - 16 + 0.3798402483 i .

Here, we regard 1.9 × 10 - 16 as 0. With the help of computer, we can get that all the real parts of λ i ( i = 1 , 2 , 3 , 4 ) are negative when g - Na < g - Na * . And there exist positive real parts when g - Na > g - Na * . According to the stability theory, the system is stable around equilibrium when g - Na [ 0 , g - Na * ) and it is unstable when g - Na ( g - Na * , 500 ] . The model undergoes Hopf bifurcation at g - Na = g - Na * .

Figure 3 shows the response of V and m , h , and n to different g - Na . As the analysis above, when g - Na = 198 < g - Na * , the system is stable. Figure 3(a) is the potential-time ( V - t ) curve, which shows that the action potential V becomes steady. Figure 3(b) displays the trajectory of gating variables m , h , and n with time. We can see that the electrophysiological activity of cell reaches an equilibrium state at last.

The response of V and m , h , and n to different g - Na .

The V - t curve ( g - Na = 198 < g - Na * )

The trajectory of m , h , and n ( g - Na = 198 < g - Na * )

The V - t curve ( g - Na = 250 > g - Na * )

The trajectory of m , h , and n ( g - Na = 250 > g - Na * )

Figures 3(c) and 3(d) demonstrate that the system is unstable when g - Na = 250 > g - Na * . Figure 3(c) is the V - t graph, from which we can see the potential changes periodically. Figure 3(d) describes the trend of m , h , and n , whose trajectory is a circle finally. Both Figures 3(c) and 3(d) imply that the system is unstable and the electrophysiological activity of cell is in a periodical state at a certain frequency.

3.2. Effect of <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M147"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mover accent="true"> <mml:mrow> <mml:mi>g</mml:mi></mml:mrow> <mml:mo>-</mml:mo></mml:mover> </mml:mrow> <mml:mrow> <mml:mtext>K</mml:mtext></mml:mrow> </mml:msub></mml:mrow> </mml:math></inline-formula> on the Model

In this part, we choose g - K as variable and keep the other parameters with ideal values. The same method with analysis of g - Na is taken to analyze the effect of g - K on the equilibrium, stability, and bifurcation of HH model. g - K [ 0,200 ] is taken into consideration because the desired value of g - K is 36.0. First, the relationship between g - K and the equilibrium V * is obtained in Figure 4.

The relationship between g - K and V * .

From Figure 4, we can see that V * varies rapidly when g - K [ 0,20 ] and decreases slowly when g - K [ 30,200 ] . This means that equilibrium points are sensitive to g - K when g - K [ 0,20 ] . A slight change of g - K may make the final state of model change greatly.

Using the method of bisection to calculate the eigenvalues, we can find two bifurcation points g - K * 1 = 3.843499029 and g - K * 2 = 19.762260771 when g - K varies. Substitute g - K * i ( i = 1 , 2 ) into (1), and obtain the corresponding equilibrium points V * . Both g - K * i and V * are substituted into (7), and then corresponding eigenvalues could be obtained as follows: (10) λ 1 1 = - 5.3218099843 , λ 2 1 = - 0.4223840650 , λ 3 1 = 3.2 × 10 - 16 - 1.1305093754 i , λ 4 1 =    3.2 × 10 - 16 + 1.1305093754 i , λ 1 2 = - 4.5370272278 , λ 2 2 = - 0.1319002182 , λ 3 2 = 4.2 × 10 - 16 - 0.3436440068 i , λ 4 2 = 4.2 × 10 - 16 + 0.3436440068 i .

Here, 3.2 × 10 - 16 and 4.2 × 10 - 16 can be approximately regarded as 0. From computing, all the real parts of eigenvalues are negative when g - K [ 0 , g - K * 1 ) ( g - K * 2 , 200 ] , and all of them are not negative when g - K ( g - K * 1 , g - K * 2 ) . According to the stability theory, the system is stable around equilibrium when g - K [ 0 , g - K * 1 ) ( g - K * 2 , 200 ] and it is unstable when g - K ( g - K * 1 , g - K * 2 ) . The model undergoes Hopf bifurcations at g - K = g - K * i ( i = 1 , 2 ). The system is from locally stable state ( g - K [ 0 , g - K * 1 ) ) to unstable state ( g - K ( g - K * 1 , g - K * 2 ) ) and becomes stable ( g - K ( g - K * 2 , 200 ] ) again. Responses of V and m , h , and n to different g - K are shown in Figure 5.

The response of V and m , h , and n to different g - K .

The V - t curve ( g - K = 2.8 < g - K * 1 )

The trajectory of m , h , and n ( g - K = 2.8 < g - K * 1 )

The V - t curve ( g - K * 1 < g - K = 15 < g - K * 2 )

The trajectory of m , h , and n ( g - K * 1 < g - K = 15 < g - K * 2 )

The V - t curve ( g - K = 21 > g - K * 2 )

The trajectory of m , h , and n ( g - K = 21 > g - K * 2 )

Figure 5 shows the response of V and m , h , and n to different g - K . When g - K = 2.8 < g - K * 1 , the system is stable. Figure 5(a) shows the trend of potential with time, from which we can see that the potential reaches a fixed value. Figure 5(b) is the trajectory of m , h , and n with time. All the gating variables also stay at fixed values (a steady point in Figure 5(b)) at last. These mean that the electrophysiological activity of cell reaches a steady state ultimately.

Figures 5(c) and 5(d) are V - t and m - h - n graphs, respectively, when g - K * 1 < g - K = 15 < g - K * 2 . Figure 5(c) shows that the action potential changes in a certain period. And Figure 5(d) describes the trajectory of m , h , and n with time, from which we can find that the shape of the trajectory is a loop. Figures 5(c) and 5(d) imply all the gating variables and potential change periodically, which means that the electrophysiological activity of cell is in a periodical state.

Figures 5(e) and 5(f) are V - t and m - h - n curves, respectively, when g - K = 21 > g - K * 2 . From Figure 5(e) we can see that the potential reaches the resting state at this occasion. Figure 5(f) describes the trajectory of m , h , and n , which shows that the three variables stay at a fixed point at last. Both Figures 5(e) and 5(f) show that all the potential and gating variables no longer change with time, which implies the cell reaches the resting state finally.

3.3. Effect of <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M232"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mover accent="true"> <mml:mrow> <mml:mi>g</mml:mi></mml:mrow> <mml:mo>-</mml:mo></mml:mover> </mml:mrow> <mml:mrow> <mml:mtext>Na</mml:mtext></mml:mrow> </mml:msub></mml:mrow> </mml:math></inline-formula> and <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M233"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mover accent="true"> <mml:mrow> <mml:mi>g</mml:mi></mml:mrow> <mml:mo>-</mml:mo></mml:mover> </mml:mrow> <mml:mrow> <mml:mtext>K</mml:mtext></mml:mrow> </mml:msub></mml:mrow> </mml:math></inline-formula> on the Model

Both g - Na and g - K are taken as variables in this part to study the stability and bifurcation of the model when g - Na [ 0,400 ] and g - K [ 0,60 ] . Keeping the other parameters with desired values, using Matlab as a tool, we get the equilibrium points first when g - Na and g - K both vary. Then the points are substituted into eigenmatrix of (4) and the eigenvalues of the model can be calculated. At last, Figure 6 is gained, in which all the real parts of eigenmatrix are negative if g - Na and g - K belong to the pink region and positive real parts appear if g - Na and g - K are in white area.

The g - Na - g - K plane and the critical boundary.

From Figure 6, we can find that the upper boundary of the regions is similar to a line. With the least square method applied, the expression of the line can be gotten as g - K = 0.175 × g - Na - 1.675 . However, the lower boundary is not regular. According to stability theory, we can easily know that the model is stable when g - Na and g - K are in pink region and unstable when g - Na and g - K are in white. This means that the electrophysiological activity can reach a steady state when g - Na and g - K are in pink region and it is periodic when g - Na and g - K are in white. The system undergoes bifurcations when g - Na and g - K are on the boundary.

4. Discussion and Conclusion

The effects of g - Na and g - K on the stability and bifurcation of HH model are analyzed in the paper. The critical values and boundary are obtained. When g - Na increases to the critical value, the model will have bifurcation phenomenon, which means system will reach stable state when g - Na is less than the critical value and the cell will have continuous action potential after stimulation when g - Na is greater. However, there are two critical values about g - K . The system will be stable when g - K is less than the smaller critical value and there are periodic solutions when g - K is greater than the value and meanwhile is less than the larger one. The model will reach steady state again when g - K is greater than the larger critical value. From analyzing g - Na and g - K collectively, we can get a critical line which divides the g - Na - g - K plane into two parts. The system will be stable when ( g - Na , g - K ) is in the upper half plane and model will have periodic solutions when ( g - Na , g - K ) is in the lower half.

In our analysis, when g - Na or g - K are taken as the variable(s), all the other parameters are kept with desired values. However, almost all the biological systems are coupled. All the components influence one another and work together forming the overall functionality. Therefore, when the sodium ( g - Na ) and/or potassium ( g - K ) channels vary, do the other parameters remain unchanged? Is it reasonable to keep the other parameters still with desired values? We could not ensure that it must be reasonable. Nevertheless, some evidences may explain a certain rationality of the method. For example, tetrodotoxin (TTX) selectively binds to the outer vestibule voltage-gated sodium channels, preventing channels from opening . Ivabradine is a sinus node I f channel inhibitor, which is selective for the I f current but does not affect other cardiac ionic currents . Acacetin could suppress the ultrarapid delayed rectifier K+ current and the transient outward K+ current and block the acetylcholine-activated K+ current; however, it has no effect on Na+ current, L-type Ca2+ current, or even inward-rectifier K+ current . All of these demonstrate that to an extent when one channel changes, the others may not be affected. That is, when the parameter describing a channel varies, it is reasonable to keep parameters describing the other channels unchanged.

Stable states indicate that the electrophysiological activity of cell will get to corresponding resting state at last, while periodic phenomenon looks like response of pathological cell’s action potentials caused by cardiac arrhythmias . In other words, g - Na and g - K may be the causes of the similar diseases to cardiac arrhythmias. So given appropriate medicine to change g - Na or g - K to reasonable intervals, the corresponding diseases could be abolished or the discomfort can be ameliorated. After all, our research could be a reference to treat relevant diseases. Some diseases led to by abnormal ion channels may be eased by medicine to adjust the conductance into corresponding intervals.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work is supported by the National Natural Science Foundation of China (NSFC) under Grant no. 61173086 and no. 61179009.

Hodgkin A. L. Huxley A. F. The components of membrane conductance in the giant axon of Loligo The Journal of Physiology 1952 116 4 473 496 2-s2.0-76949125485 Hodgkin A. L. Huxley A. F. The components of membrane conductance in the giant axon of Loligo The Journal of Physiology 1952 116 4 473 496 2-s2.0-76949125485 Hodgkin A. L. Huxley A. F. The dual effect of membrane potential on sodium conductance in the giant axon of Loligo The Journal of Physiology 1952 116 4 497 506 Hodgkin A. L. Huxley A. F. A quantitative description of membrane current and its application to conduction and excitation in nerve. The Journal of Physiology 1952 117 4 500 544 2-s2.0-35649001607 Soeller C. Jayasinghe I. D. Li P. Holden A. V. Cannell M. B. Three-dimensional high-resolution imaging of cardiac proteins to construct models of intracellular Ca2+ signalling in rat ventricular myocytes Experimental Physiology 2009 94 5 496 508 10.1113/expphysiol.2008.043976 2-s2.0-64749099928 Wussling M. Szymanski G. Simulation by two calcium store models of myocardial dynamic properties: potentiation, staircase, and biphasic tension development General Physiology and Biophysics 1986 5 2 135 152 2-s2.0-0022702616 Mahajan A. Shiferaw Y. Sato D. Baher A. Olcese R. Xie L. Yang M. Chen P. Restrepo J. G. Karma A. Garfinkel A. Qu Z. Weiss J. N. A rabbit ventricular action potential model replicating cardiac dynamics at rapid heart rates Biophysical Journal 2008 94 2 392 410 10.1529/biophysj.106.98160 2-s2.0-38349013614 Noble D. Cardiac action and pacemaker potentials based on the Hodgkin-Huxley equations Nature 1960 188 4749 495 497 10.1038/188495b0 2-s2.0-36949064203 Kneller J. Ramirez R. J. Chartier D. Courtemanche M. Nattel S. Time-dependent transients in an ionically based mathematical model of the canine atrial action potential American Journal of Physiology—Heart and Circulatory Physiology 2002 282 4 H1437 H1451 2-s2.0-0036084519 Zhang H. Holden A. V. Kodama I. Honjo H. Lei M. Varghese T. Boyett M. R. Mathematical models of action potentials in the periphery and center of the rabbit sinoatrial node American Journal of Physiology: Heart and Circulatory Physiology 2000 279 1 H397 H421 2-s2.0-0033858925 Verkerk A. O. Wilders R. Hyperpolarization-activated current, If, in mathematical models of rabbit sinoatrial node pacemaker cells BioMed Research International 2013 2013 18 872454 10.1155/2013/872454 2-s2.0-84880914083 Butters T. D. Aslanidi O. V. Zhao J. Smaill B. Zhang H. A novel computational sheep atria model for the study of atrial fibrillation Interface Focus 2013 3 2 20120067 10.1098/rsfs.2012.0067 2-s2.0-84875416142 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 Computational Biology 2011 7 5 e1002061 10.1371/journal.pcbi.1002061 2-s2.0-79958158591 Li P. Rudy Y. A model of canine purkinje cell electrophysiology and Ca2+ cycling: rate dependence, triggered activity, and comparison to ventricular myocytes Circulation Research 2011 109 1 71 79 10.1161/CIRCRESAHA.111.246512 2-s2.0-79959984154 Hassard B. Bifurcation of periodic solutions of the Hodgkin-Huxley model for the squid giant axon Journal of Theoretical Biology 1978 71 3 401 420 10.1016/0022-5193(78)90168-6 MR0479443 2-s2.0-0017854654 Rinzel J. Miller R. N. Numerical calculation of stable and unstable periodic solutions to the Hodgkin-Huxley equations Mathematical Biosciences 1980 49 1 27 59 10.1016/0025-5564(80)90109-1 2-s2.0-0019310660 Aihara K. Matsumoto G. Two stable steady states in the Hodgkin-Huxley axons Biophysical Journal 1983 41 1 87 89 10.1016/S0006-3495(83)84408-7 2-s2.0-0020538833 Guckenheimer J. Labouriau J. S. Bifurcation of the Hodgkin and Huxley equations: a new twist Bulletin of Mathematical Biology 1993 55 5 937 952 10.1007/BF02460693 2-s2.0-0000370809 Bedrov Y. A. Akoev G. N. Dick O. E. On the relationship between the number of negative slope regions in the voltage-current curve of the Hodgkin-Huxley model and its parameter values Biological Cybernetics 1995 73 2 149 154 10.1007/BF00204053 ZBLl0826.92003 2-s2.0-0029339680 Bedrov Y. A. Akoev G. N. Dick O. E. Partition of the Hodgkin-Huxley type model parameter space into the regions of qualitatively different solutions Biological Cybernetics 1992 66 5 413 418 10.1007/BF00197721 ZBLl0745.92004 2-s2.0-0026473629 Fukai H. Doi S. Nomura T. Sato S. Hopf bifurcations in multiple-parameter space of the Hodgkin-Huxley equations I. Global organization of bistable periodic solutions Biological Cybernetics 2000 82 3 215 222 10.1007/s004220050021 ZBLl0962.92005 2-s2.0-0034145622 Terada K. Tanaka H. A. Yoshizawa S. Two-parameter bifurcations in the Hodgkin–Huxley equations for muscle fibers Electronics and Communications in Japan 2000 83 6 86 94 Jiang W. Yanqiu C. Xiangyang F. Li L. Multi-parameter Hopf-bifurcation in Hodgkin-Huxley model exposed to ELF external electric field Chaos, Solitons & Fractals 2005 26 4 1221 1229 10.1016/j.chaos.2005.02.042 2-s2.0-19144372247 Wang J. Chen L. Fei X. Analysis and control of the bifurcation of Hodgkin-Huxley model Chaos, Solitons and Fractals 2007 31 1 247 256 10.1016/j.chaos.2005.09.060 ZBLl1140.37370 2-s2.0-33745854291 Wang J. Chen L. Q. Fei X. Y. Bifurcation control of the Hodgkin-Huxley equations Chaos, Solitons & Fractals 2007 33 1 217 224 10.1016/j.chaos.2006.01.035 2-s2.0-33846827228 Che Y. Wang J. Deng B. Wei X. Han C. Bifurcations in the Hodgkin-Huxley model exposed to DC electric fields Neurocomputing 2012 81 41 48 10.1016/j.neucom.2011.11.019 2-s2.0-84856333741 Hao P. M. Fan D. J. Wei J. J. Liu Q. Dynamic behaviors of a delayed HIV model with stage-structure Communications in Nonlinear Science and Numerical Simulation 2012 17 12 4753 4766 10.1016/j.cnsns.2012.04.004 2-s2.0-84864414005 Chang X. Y. Wei J. J. Hopf bifurcation and optimal control in a diffusive predator-prey system with time delay and prey harvesting Nonlinear Analysis: Modelling and Control 2012 17 4 379 409 2-s2.0-84870332531 Chang X. Y. Wei J. J. Stability and Hopf bifurcation in a diffusive predator-prey system incorporating a prey refuge Mathematical Biosciences and Engineering 2013 10 4 979 996 10.3934/mbe.2013.10.979 Durand D. M. Bikson M. Suppression and control of eplleptiform activity by electrical stimulation: a review Proceedings of the IEEE 2001 89 7 1065 1082 10.1109/5.939821 2-s2.0-0000447520 Venkatesh B. Lu S. Q. Dandona N. See S. L. Brenner S. Soong T. W. Genetic basis of tetrodotoxin resistance in pufferfishes Current Biology 2005 15 22 2069 2072 10.1016/j.cub.2005.10.068 2-s2.0-27844446184 Sulfi S. Timmis A. D. Ivabradine—the first selective sinus node If channel inhibitor in the treatment of stable angina International Journal of Clinical Practice 2006 60 2 222 228 10.1111/j.1742-1241.2006.00817.x 2-s2.0-33646674434 Li G. Wang H. Qin G. Jin M. Tang Q. Sun H. Du X. Deng X. Zhang X. Chen J. Chen L. Xu X. Cheng L. Chiu S. Tse H. Vanhoutte P. M. Lau C. Acacetin, a natural flavone, selectively inhibits human atrial repolarization potassium currents and prevents atrial fibrillation in dogs Circulation 2008 117 19 2449 2457 10.1161/CIRCULATIONAHA.108.769554 2-s2.0-43449090419 Wang H. O. Chen D. Bushnell L. G. Control of bifurcations and chaos in heart rhythms 1 Proceedings of the 36th IEEE Conference on Decision and Control December 1997 San Diego, Calif, USA 395 400 2-s2.0-0031386590 10.1109/CDC.1997.650655