Applied Mathematics
Vol.07 No.09(2016), Article ID:66767,5 pages
10.4236/am.2016.79072

The Relationship of Sodium and Potassium Conductances with Dynamic States of a Mathematical Model of Electrosensory Pyramidal Neurons

Takaaki Shirahata

Kagawa School of Pharmaceutical Sciences, Tokushima Bunri University, Sanuki, Japan

Copyright © 2016 by author and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 30 March 2016; accepted 23 May 2016; published 26 May 2016

ABSTRACT

Electrosensory pyramidal neurons in weakly electric fish can generate burst firing. Based on the Hodgkin-Huxley scheme, a previous study has developed a mathematical model that reproduces this burst firing. This model is called the ghostbursting model and is described by a system of nonlinear ordinary differential equations. Although the dynamic state of this model is a quiescent state during low levels of electrical stimulation, an increase in the level of electrical stimulation transforms the dynamic state first into a repetitive spiking state and finally into a burst firing state. The present study performed computer simulation analysis of the ghostbursting model to evaluate the sensitivity of the three dynamic states of the model (i.e., the quiescent, repetitive spiking, and burst firing states) to variations in sodium and potassium conductance values of the model. The present numerical simulation analysis revealed the sensitivity of the electrical stimulation threshold required for eliciting the burst firing state to variations in the values of four ionic conductances (i.e., somatic sodium, dendritic sodium, somatic potassium, and dendritic potassium conductances) in the ghostbursting model.

Keywords:

Mathematical Model, Computer Simulation, Ghostbursting, Ionic Conductance

1. Introduction

The ghostbursting model is a mathematical model of electrosensory pyramidal neurons in weakly electric fish, which is described by a system of nonlinear Ordinary Differential Equations (ODEs) (see Methods in [1] ). This model is based on the Hodgkin-Huxley formalism and describes the time evolution of the membrane potentials of the somatic and dendritic compartments of the model. This model contains many parameters such as electrical stimulation and ionic conductances (i.e., sodium and potassium conductances), and previous studies have revealed the relationship of the dynamic states of the model with variations in parameter values. An increase in the amplitude of the electrical stimulation to the somatic compartment changes the dynamic state of the model from a quiescent state to a repetitive spiking state and finally to a bursting firing state [1] . The sensitivity of the electrical stimulation thresholds required for inducing the repetitive spiking and bursting firing states to variations in the potassium conductance of the dendritic compartment has been characterized [1] . In addition, the sensitivity of the electrical stimulation thresholds required for inducing the repetitive spiking and bursting firing states to variations in the somatic-dendritic coupling conductance has also been reported [2] . In particular, the relationship between the potassium conductance of the dendritic compartment and the dynamic states of the model has been extensively characterized; the number of spikes per burst in a two-dimensional parameter space [3] and the influence of the kinetics of the potassium conductance of the dendritic compartment on the dynamics of the model [4] have been reported. Results from these previous studies highlight the importance of extending these investigations to studies of the sensitivity of the dynamics of the ghostbursting model to parameter variations, particularly focusing on detailed analysis of the membrane conductance ( [5] and page 26 in [6] ).

Studies of the characteristics of the potassium conductance of the dendritic compartment have previously been carried out, as described above. However, characteristics of the potassium conductance of the somatic compartment and the sodium conductances of the somatic/dendritic compartments have not been investigated in detail. Doiron and coworkers have implied that variations in values of the sodium conductances of the somatic/ dendritic compartments or the potassium conductance of the somatic compartment may affect electrical stimulation thresholds (see the last sentence of Discussion of [1] ), but this relationship has yet to be explored in detail. A systematic evaluation of the sodium and potassium conductances of both somatic and dendritic compartments is necessary for a thorough understanding of the difference between these conductances. Therefore, the present study performed computer simulation analysis of the ghostbursting model to reveal the sensitivity of the electrical stimulation thresholds to variations in the sodium and potassium conductances.

2. Materials and Methods

The ghostbursting model used in this study is described by a system of ODEs, which consists of six state variables: the membrane potential of the somatic compartment [Vs(t) (mV)] [t is time (ms)], the activating variable of the potassium conductance of the somatic compartment [ns(t)], the membrane potential of the dendritic compartment [Vd(t) (mV)], the inactivating variable of the sodium conductance of the dendritic compartment [hd(t)], the activating variable of the potassium conductance of the dendritic compartment [nd(t)], and the inactivating variable of the potassium conductance of the dendritic compartment [pd(t)]. The dynamic states of the ODEs can change depending on the following system parameters: electrical stimulation of the somatic compartment (Is), the maximal sodium conductance of the somatic compartment (gNa,s), the maximal potassium conductance of the somatic compartment (gDr,s), the maximal sodium conductance of the dendritic compartment (gNa,d), and the maximal potassium conductance of the dendritic compartment (gDr,d). The ODEs that describe the ghostbursting model are described as follows:

(1)

(2)

(3)

(4)

(5)

(6)

Is (μA/cm2) is changed from 5.6 to 9.6. The default values of conductances are gNa,s = 55 mS/cm2, gDr,s = 20 mS/cm2, gNa,d = 5 mS/cm2, and gDr,d = 15 mS/cm2. Detailed explanations of the ODEs are described in [1] .

The free and open source software Scilab (http://www.scilab.org/) was used to numerically solve the ODEs (initial conditions: Vs = −70 mV, ns = 0.00005, Vd = −70 mV, hd = 0.973, nd = 0.002, and pd = 0.697). Solving the equations was performed using the lsoda solver implemented in the program Scilab.

3. Results

Under conditions in which all conductance values were default values [see the columns of 100% conductance in Figures 1(a)-(d)], the dynamic states of the ghostbursting model were a quiescent state (× in Figure 1) when Is

Figure 1. The dependence of the dynamic states of the ghostbursting model on various parameters. The sensitivity of the dynamic states to variations in (a) Is and gNa,s, (b) Is and gNa,d, (c) Is and gDr,s, (d) Is and gDr,d. ×, quiescent state; ○, repetitive spiking state; ●, bursting state.

was 5.6, a repetitive spiking state (○ in Figure 1) when Is was between 5.8 and 8.4, and a bursting firing state (● in Figure 1) when Is was between 8.6 and 9.6. The repetitive spiking threshold (Is = 5.8) did not change even if gNa,s decreased to 95% or increased to 105% of the default value [Figure 1(a)]. The bursting threshold (Is = 8.6) was sensitive to variations in gNa,s: it decreased to 8.2 when gNa,s decreased to 95% of the default value, whereas it increased to 9.0 when gNa,s increased to 105% of the default value [Figure 1(a)]. The repetitive spiking threshold (Is = 5.8) increased to 6.0 when gNa,d decreased to 95% but did not change when gNa,d increased to 105% [Figure 1(b)]. The bursting threshold (Is = 8.6) was sensitive to variations in gNa,d; it increased to 9.0 when gNa,d decreased to 95%, whereas it decreased to 8.2 when gNa,d increased to 105% [Figure 1(b)]. The repetitive spiking threshold (Is = 5.8) did not change even if gDr,s decreased to 90% or increased to 110% [Figure 1(c)]. The bursting threshold (Is = 8.6) was sensitive to variations in gDr,s; it increased to 8.8 when gDr,s decreased to 90%, whereas it decreased to 8.4 when gDr,s increased to 110% [Figure 1(c)]. The repetitive spiking threshold (Is = 5.8) did not change even if gDr,d decreased to 95% or increased to 105% [Figure 1(d)]. The bursting threshold (Is = 8.6) was sensitive to variations in gDr,d; it decreased to 7.8 when gDr,d decreased to 95%, whereas it increased to 9.6 when gDr,d increased to 105% [Figure 1(d)].

4. Discussion

The present numerical simulation analysis revealed the sensitivity of the repetitive spiking threshold and the bursting threshold to variations in the values of four ionic conductances (i.e., gNa,s, gNa,d, gDr,s, and gDr,d) in the ghostbursting model. A previous study illustrated the sensitivity of the repetitive spiking threshold and the bursting threshold to variations in gDr,d [1] . In addition, the same study implied that these thresholds are also sensitive to variations in other ionic conductances such as gNa,s, gNa,d, and gDr,s. However, this was not shown explicitly [1] .

The importance of the present findings are that they clearly demonstrate the sensitivity of the repetitive spiking threshold and the bursting threshold to variations in gNa,s, gNa,d, and gDr,s. Specifically, 1) similar to the case of gDr,d, the repetitive spiking threshold is insensitive to variations in the other ionic conductances, except gNa,d; 2) the bursting threshold is sensitive to variations in the other three ionic conductances; 3) similar to the previous case, in which an increase in gDr,d increases the bursting threshold [1] , an increase in gNa,s increases the bursting threshold, whereas an increase in gNa,d and gDr,s decreases the bursting threshold; and 4) the degree of sensitivity occurs as follows: gDr,d > gNa,s = gNa,d > gDr,s.

The effect of variations in ionic conductance values on the behaviors of mathematical models of excitable cells is an important topic of investigation. For example, action potential duration is differentially modulated by variations in the slow-inward calcium conductance (Gsi) versus the delayed rectifier potassium conductance (GK) in the LR1 model [7] . Similarly, the distinct roles of several ionic conductances (GCaS, Gh, GKd, GA, and GKCa) in regulating maximal gain modulation is revealed in a model of lobster somatogastric neurons [8] . The present study advances our understanding of membrane dynamics in pyramidal neurons by revealing differences in bursting threshold changes among gNa,s, gNa,d, gDr,s, and gDr,d. Notably, the present study examines conductances of the same ion type in both somatic and dendritic compartments (i.e., the model incorporates differences between gNa,s and gNa,d, and between gDr,s and gDr,d). Previous studies described above [7] [8] investigated a single- compartment conductance-based model, and thus, could not reveal the differences in ionic conductances between different compartments.

5. Conclusion

Previous work has reported the sensitivity of the repetitive spiking threshold and the bursting threshold to variations in gDr,d but has not characterized the sensitivity of these thresholds to variations in three other ionic conductances (gNa,s, gNa,d, and gDr,s) [1] . The present study addresses this issue and contributes to a thorough understanding of differential modulation of these thresholds by variations in the four ionic conductances (i.e., gNa,s, gNa,d, gDr,s, and gDr,d) in the ghostbursting model.

Acknowledgements

The author would like to thank Enago (www.enago.jp) for the English language review.

Cite this paper

Takaaki Shirahata, (2016) The Relationship of Sodium and Potassium Conductances with Dynamic States of a Mathematical Model of Electrosensory Pyramidal Neurons. Applied Mathematics,07,819-823. doi: 10.4236/am.2016.79072

References

  1. 1. Doiron, B., Laing, C., Longtin, A. and Maler, L. (2002) Ghostbursting: A Novel Neuronal Burst Mechanism. Journal of Computational Neuroscience, 12, 5-25.
    http://dx.doi.org/10.1023/A:1014921628797

  2. 2. Laing, C., Doiron, B., Longtin, A. and Maler, L. (2002) Ghostbursting: The Effects of Dendrites on Spike Patterns. Neurocomputing, 44-46, 127-132.
    http://dx.doi.org/10.1016/S0925-2312(02)00373-9

  3. 3. Shirahata, T. (2012) Analysis of the Electrosensory Pyramidal Cell Bursting Model for Weakly Electric Fish: Model Prediction under Low Levels of Dendritic Potassium Conductance. Acta Biologica Hungarica, 63, 313-320.
    http://dx.doi.org/10.1556/ABiol.63.2012.3.1

  4. 4. Shirahata, T. (2015) Evaluation of Kinetic Properties of Dendritic Potassium Current in Ghostbursting Model of Electrosensory Neurons. Applied Mathematics, 6, 128-135.
    http://dx.doi.org/10.4236/am.2015.61013

  5. 5. Shirahata, T. (2015) Numerical Study of a Mathematical Model of Vibrissa Motoneurons: The Relationship between Repetitive Spiking and Two Types of Sodium Conductance. International Journal of Theoretical and Mathematical Physics, 5, 48-52.

  6. 6. Ashrafuzzaman, M. and Tuszynski, J. (2012) Membrane Biophysics. Springer, Berlin.

  7. 7. Sobie, E.A. (2009) Parameter Sensitivity Analysis in Electrophysiological Models Using Multivariable Regression. Biophysical Journal, 96, 1264-1274.
    http://dx.doi.org/10.1016/j.bpj.2008.10.056

  8. 8. Patel, A.X. and Burdakov, D. (2015) Mechanisms of Gain Control by Voltage-Gated Channels in Intrinsically-Firing Neurons. PLoS ONE, 10, e0115431.
    http://dx.doi.org/10.1371/journal.pone.0115431