Type-A medial vestibular nucleus neurons (mVNn) can show spontaneous repetitive spiking without stimulation in vitro. A mathematical model, which reproduces the spiking, has been proposed based on the Hodgkin-Huxley formulation [1] . This model, which is described by a system of nonlinear ordinary differential equations (ODEs), contains several ionic conductances such as the sodium conductance, calcium conductance, delayed-rectifier potassium conductance, transient potassium conductance, and calcium-dependent potassium conductance. The study that proposed this model also indicated that the spontaneous spiking activity of the type-A mVNn model can be suppressed by hyperpolarizing stimulation [1] . However, it is important to conduct detailed analysis of the ionic conductances of excitable cells [2] . Thus, it is important to clarify how the different ionic conductances in the type-A mVNn model each affect the suppression of spiking. However, this issue was not been investigated in the previous study [1] .

Several studies have investigated the characteristics of the ionic conductances of various excitable cell models. For example, various studies have investigated the effect of co-variation of ionic conductances on the dynamics of excitable cell models [3] - [12] . In addition, studies have investigated the effect of variation of ionic conductances on the electrical stimulation-induced changes in the dynamical states [13] - [16] . The above studies have contributed to an understanding of the differences between the characteristics among different ionic conductances. Although analysis of ionic conductances has been extensively performed as described above, analysis of the ionic conductances of the type-A mVNn model has not been performed. Investigation of the type-A mVNn model will contribute to an understanding of the difference in the characteristics among ionic conductances of the type-A mVNn model in more detail. Therefore, the present study performed numerical simulation of the type-A mVNn model to clarify how the suppression of the spiking is dependent on the different ionic conductances.

The present study performed numerical simulations of a mathematical model of type-A mVNn, which was developed in a previous study [1] . The model is described by a system of nonlinear ODEs, which consists of the following five state variables: the membrane potential of the mVNn [V (mV)], three gating variables of ionic currents (n, x, and b), and calcium concentration {[Ca] (mM)}. The time evolution of the state variables is described as follows:

where C_m(1 mF/cm²) is the membrane capacitance; I_app is the externally injected current of constant amplitude; I_Na(V, n), I_Ca(V, x, [Ca]), I_K(V, n), I_KCa(V, [Ca]), I_A(V, b), and I_L(V) are the sodium current, calcium current, delayed-rectifier potassium current, calcium-dependent potassium current, transient potassium current, and leak current, respectively, which are defined in Equations (6)-(11) below;t_X(V) (ms) (X = n, x, b) and X_∞(V) (X = n, x, b) are the time constants of activation/inactivation and the steady-state activation/inactivation functions, respectively, which are defined in Equations (12)-(17) below; K_p (=1) is the calcium influx parameter; and R (=5) is the parameter that regulates calcium removal through calcium transport and mitochondrial uptake.

where g_Na, g_Ca, g_K, g_KCa, g_A, and g_L (=0.3 mS/cm²) are the maximal conductances of I_Na(V, n), I_Ca(V, x, [Ca]), I_K(V, n), I_KCa(V, [Ca]), I_A(V, b), and I_L(V), respectively, and V_Na(=55 mV), V_Ca(=124 mV), V_K(=−80 mV), and V_L(=−50 mV) are the reversal potentials of I_Na(V, n), I_Ca(V, x, [Ca]), three types of potassium currents {I_K(V, n), I_KCa(V, [Ca]), and I_A(V, b)}, and I_L(V), respectively. Refer to [1] for detailed explanations of the equations.

The free and open source software, Scilab (http://www.scilab.org/), was used to numerically solve the above ODEs (initial conditions: V = −60 mV, n = 0.1, x = 0.1, b = 0.9, and [Ca] = 0.1 mM). The total simulation time was 0.6 s in all the simulations. System parameters whose values were varied were I_app, g_Na, g_Ca, g_K, g_KCa, and g_A. I_app between 0.2 s and 0.6 s was varied from −2.0 mA/cm² to 2.5 mA/cm² at an interval of 0.5 mA/cm², whereas I_app between 0.0 s and 0.2 s was set to zero. Default values of g_Na, g_Ca, g_K, g_KCa, and g_A were 20 mS/cm², 1 mS/cm², 2 mS/cm², 1 mS/cm², and 4 mS/cm², respectively. Each of g_Na, g_Ca, g_K, g_KCa, and g_A was varied to 50% or 150% of each default value.

The previous study had indicated that a hyperpolarizing stimulation of −3 mA/cm² suppressed spontaneous repetitive spiking [1] . To determine the effect of the hyperpolarizing stimulation on spiking activity in more detail, the present study first investigated the response of the model to hyperpolarizing stimulations of smaller amplitudes under conditions in which all the ionic conductances were set to default values. When a stimulation of −0.5 mA/cm² was applied to the model, it showed a repetitive spiking behavior, but the frequency of spiking became smaller than that before the stimulation (Figure 1(a)). In contrast, when a stimulation of −1.0 mA/cm² was applied to the model, it did not show a repetitive spiking behavior, but showed a quiescent state (Figure 1(b)). Based on these results, the present study defined −1.0 mA/cm² as the stimulation threshold for suppressing repetitive spiking (we hereafter referred to this threshold as the “suppression threshold”). To reveal the relationship between the ionic conductances and the dynamics of the model, the present study next investigated how the suppression threshold changed in response to variations of each ionic conductance.

The effects of variations in the ionic conductance values on the dynamics of the model are shown in Figure 2. Under conditions in which g_Na was 50% of the default value with other conductance values being set at the de-

fault values, the model showed a quiescent state in response to stimulation from −2.0 mA/cm² to 2.0 mA/cm², whereas it showed a repetitive spiking state in response to stimulation of 2.5 mA/cm² (Figure 2(a)). Therefore, the suppression threshold at 50% g_Na was 2.0 mA/cm². Figure 2(a) also shows that the suppression thresholds at 100% g_Na and 150% g_Na were −1.0 mA/cm² and −2.0 mA/cm², respectively. Similarly, the suppression thresholds under conditions in which other ionic conductance values were varied were calculated. The suppression threshold at 50% g_Ca, 100% g_Ca, and 150% g_Ca was −1.5 mA/cm², −1.0 mA/cm², and −0.5 mA/cm², respectively (Figure 2(b)). The suppression threshold at 50% g_A, 100% g_A, and 150% g_A was −1.5 mA/cm², −1.0 mA/cm², and −0.5 mA/cm², respectively (Figure 2(c)). The suppression threshold at 50% g_KCa, 100% g_KCa, and 150% g_KCa was −1.5 mA/cm², −1.0 mA/cm², and 0.0 mA/cm², respectively (Figure 2(d)). Finally, the suppression threshold at 50%, 100%, and 150% g_K was −1.0 mA/cm² for all three conditions (Figure 2(e)).

The present study performed numerical simulation of the type-A mVNn model and revealed the differential sensitivity of the threshold for the transition from a repetitive spiking state to a quiescent state to variations in the different ionic conductance values. The present study revealed the difference in modulation of the threshold among the ionic conductances. An increase in g_Na lowered the threshold, whereas an increase in g_Ca, g_A, or g_KCaraised the threshold. Finally, an increase in g_K had no effect on the threshold. The sensitivity of the threshold to variations in the ionic conductance values decreased in the following order: g_Na > g_KCa > g_Ca = g_A > g_K. Although the previous study had investigated the response of the type-A mVNn model to hyperpolarizing and depolarizing stimulation [1] , the relationship between the ionic conductances and the threshold for the transition from a repetitive spiking state to a quiescent state had not been clarified prior to the present study.

It is important to determine whether or not variations in system parameters affect mathematical models of neurons linearly. For example, analysis of the ghostbursting model has reported that changes in the dendritic potassium conductance kinetics induce anonlinear effect on the dynamical states in a two-dimensional parameter space of the ghostbursting model [17] . In contrast, in the ghostbursting model, the threshold for the transition from a repetitive spiking state to a bursting state is influenced linearly by several ionic conductances [16] . In addition, one report has indicated that whether or not variations in ionic conductance values affect a certain type of neuron model linearly depends on the type of ionic conductance: in a vibrissa motoneuron model, the threshold for the transition between a subthreshold oscillation state and a tonic firing state changes nonlinearly in response to variations in the M-type potassium and the persistent sodium conductances [15] , whereas the threshold for the transition between a subthreshold oscillation state and a mixed-mode oscillation state changes linearly in response to variations in the cationic h conductance [15] . The study that developed the type-A mVNn model had not investigated the effect of variations in the ionic conductance values on the dynamics of the type-A mVNn model in detail [1] ; thus, whether or not variations in the ionic conductance values affect the type-A mVNn model linearly was unclear. The present study resolved this problem: in the type-A mVNn model, the threshold for the transition from a repetitive spiking state to a quiescent state changed linearly in response to variations in g_Ca and g_A, whereas the threshold changed highly nonlinearly in response to variations in g_Na (the threshold changed slightly nonlinearly in response to variations in g_KCa). In particular, the present study clarified that the dynamics of the type-A mVNn model was affected linearly in response to variations in the calcium conductance. However, because neither the ghostbursting model nor the vibrissa motoneuron model contains the calcium conductance variable [15] [16] , whether or not variations in the calcium conductance affects the dynamics of these two models linearly cannot be studied.

The present study performed numerical simulation of the type-A mVNn model to reveal that the threshold for the transition from a repetitive spiking state to a quiescent state was differentially sensitive to variations in the ionic conductance values among the different types of the ionic conductances, which had not been reported in the previous study of the type-A mVNn model. The present study contributes to a more detailed understanding of the difference among the ionic conductances of the type-A mVNn model.

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

Takaaki Shirahata, (2016) The Effect of Variations in Ionic Conductance Values on the Suppression of Repetitive Spiking in a Mathematical Model of Type-A Medial Vestibular Nucleus Neurons. Applied Mathematics,07,1134-1139. doi: 10.4236/am.2016.710101