Impact of slow K+ currents on spike generation can be described by an adaptive threshold model

A neuron that is stimulated by rectangular current injections initially responds with a high firing rate, followed by a decrease in the firing rate. This phenomenon is called spike-frequency adaptation and is usually mediated by slow K+ currents, such as the M-type K+ current (IM) or the Ca2+-activated K+ current (IAHP). It is not clear how the detailed biophysical mechanisms regulate spike generation in a cortical neuron. In this study, we investigated the impact of slow K+ currents on spike generation mechanism by reducing a detailed conductance-based neuron model. We showed that the detailed model can be reduced to a multi-timescale adaptive threshold model, and derived the formulae that describe the relationship between slow K+ current parameters and reduced model parameters. Our analysis of the reduced model suggests that slow K+ currents have a differential effect on the noise tolerance in neural coding.


Introduction
Neuronal adaptation is the change in the responsiveness of a neuron over time. Adaptation may play an important role in the extraction of important information from an ever-changing environment and is the product of several factors, including ion channels, synapses, and network dynamics. In this study, we focus on adaptation at the single neuron level. When a neuron is stimulated by rectangular current injections, it initially responds with a high firing rate, followed by a decrease in the firing rate. This phenomenon is called spike-frequency adaptation and is observed in most pyramidal neurons in various brain areas. The spike-frequency adaptation is usually mediated for by M-type K + current (I M ) (Brown and Adams 1980;Adams et al. 1982), Ca 2+activated K + current (I AHP ) (Brown and Griffith 1983;Madison and Nicoll 1984), Na + -activated K + current (Schwindt et al. 1989), or the slow inactivation of Na + current (Fleidervish et al. 1996;Kim and Rieke 2003). In terms of the spike-frequency adaptation generated by slow K + currents, conductance-based models including slow K + channels have been studied. These models can reproduce the electrophysiological properties of a neuron (see Koch 1999 for a review) and provide insights into the underlying biophysical mechanisms.
Studies using the conductance-based models have suggested that the distinct biophysical mechanisms responsible for the spike-frequency adaptation have different impacts on neural coding (Ermentrout et al. 2001;Prescott and Sejnowski 2008). For example, I M improves spike-timing coding, whereas I AHP improves spike-rate coding (Prescott and Sejnowski 2008). These results indicate that specific biophysical mechanisms underlying adaptation may impact the coding properties of Action Editor: J. Rinzel a neuron. On the other hand, due to the complexity of the detailed models, it remains unclear how the kinetics of slow K + currents influence the spike generation mechanism.
In order to understand the spike generation mechanism, it is essential to reduce the detailed neuron models to simplified models. There have been many attempts to simplify the detailed models (Ermentrout and Kopell 1986;Abbott and Kepler 1990;Destexhe 1997;Kistler et al. 1997;Richardson et al. 2003;Fourcaud-Trocmé et al. 2003 for review Rinzel and Ermentrout 1998;Gerstner and Kistler 2002;Izhikevich 2007). A direct approach to obtain the reduced model is to fit the simplified model to simulated data set generated by the detailed model. This approach has clarified the underlying mechanism of spike generation, such as, integration properties (Kistler et al. 1997;Jolivet et al. 2004), adaptation (Brette and Gerstner 2005), and spike threshold variability (Kobayashi and Shinomoto 2007). However, this approach cannot predict the effect of the detailed model parameters (physiological parameters) on spike generation. Another approach is to develop a mathematical framework to simplify the detailed models. For example, the FitzHugh-Nagumo model and integrate-and-fire models (Gerstner and Kistler 2002) were derived from the Hodgkin-Huxley model (Abbott and Kepler 1990;Richardson et al. 2003).
In this study, we extend the mathematical reduction approach by including the spike history effect that is essential to describe the impact of slow K + currents on spike generation. We show that the detailed conductance-based neuron model can be reduced to a multi-timescale adaptive threshold model Yamauchi et al. 2011), and derive the formulae that describe the relationship between the slow K + current parameters and the reduced model parameters. We evaluate the reduced model by predicting spike trains of the detailed model. Finally, we examine the effect of noise on the coding property of a neuron using the reduced model.

Conductance-based model
We analyzed a single-compartment conductance-based model based on a model for the cerebral cortex and thalamic neurons proposed by Pospischil et al. (2008), that was extended to include Ca 2+ -activated K + (AHP) current (Mainen and Sejnowski 1996;Tsubo et al. 2004). The membrane voltage V of a neuron is described by the following equation: where C m is the membrane capacitance and I ex is the external input current. The ionic currents consist of the leak current I L = g L (V − E L ), Na + current I Na = g Na m 3 h(V − E Na ), delayed rectifier K + current I Kd = g Kd n 4 (V − E K ), muscarinic K + cur- where g x and E x are the maximal ionic conductances and the reversal potentials, respectively. The gating variables w ∈ {m, h, n, p, q, r, s} are described by the Hodgkin − Huxley formalism.
where α w and β w are the activation and inactivation functions, respectively (see Table 1 for details), and [Ca 2+ ] represents the calcium concentration. The Ca 2+ concentration is described by (Mainen and Sejnowski 1996;Tsubo et al. 2004)  (Jolivet et al. 2004). We further confirmed that the results were quantitatively the same for a time step of 0.01 [ms].

Adaptive threshold models
The potential u of a model neuron obeys a linear differential equation, where τ m is the membrane time constant. The neuron generates a spike if the potential u reaches the spike threshold θ u (t) from below, and the threshold is linearly modulated by spikes Yamauchi et al. 2011) If u t ð Þ > θ u t ð Þ→ Emit a spike at time t; where t k is the k-th spike time, H u (t) is the threshold kernel that describes the effect of previous spikes, and the sum is taken up to the most recent spike time. The multi-timescale adaptive threshold (MAT) model ) is a special case of the adaptive threshold model (Eq. (5)). The threshold kernel is given by the sum of exponential functions for each spike in the history, where L is the number of exponential functions and α j and τ j are the weights and the threshold time constants, respectively. It is worth noting that the potential u of the adaptive threshold model is different from the voltage of the leaky integrateand-fire (LIF) model (Gerstner and Kistler 2002). The potential does not reset after a spike and continuously integrates the input current, whereas the voltage in the LIF model is reset after each spike.

Input currents
We used two types of input current I ex (t). The first input is a constant current with a pulse, where I c [μ A/cm 2 ] is the strength of the constant current, q c [nC/cm 2 ] and t p [ms] are the amplitude and timing of the pulse, respectively, and δ(t) is the Dirac's delta function. The constant part I c is tuned to maintain the membrane potential at V c and the pulse amplitude is set to shift the voltage up to −45 [mV], q c = C m (−45 − V c ). The amplitude should be large enough that the neuron always generates a spike. In all simulations, the neuron was stimulated by the pulse after it achieves the steady state. The second input is an in vivo-like current modeled by the Ornstein − Uhlenbeck process (Tuckwell 1988;Kobayashi et al. 2011), where, μ and σ are the mean and standard deviation (SD) of the input, τ syn = 2 [ms] is the synaptic time constant, and η(t) is the Gaussian white noise with zero mean and unit variance.

Calculation of the spike threshold
We evaluated the instantaneous spike threshold of the detailed conductance-based model (Eqs. (1), (2), and (3)). To evaluate the spike threshold at time t 0 , i.e., θ V (t 0 ), we stimulate the model neuron with an impulse, I ex (t) = qδ(t − t 0 ), and observe whether the model neuron generate a spike or not. The spike threshold is defined as is the voltage immediately before the pulse injection and q min is the minimal pulse amplitude for generating a spike (Fig. 1). The minimal amplitude q min can be calculated using the bisection method (Press et al. 2007). Initially, a voltage interval [a, b] is selected such that a (b) is lower (higher) than the spike threshold. The initial interval was set as [−80, 0]. Next, we check whether the midpoint c = (a + b)/2 is larger than the spike threshold by observing the voltage for 50 [ms]. If the neuron emits a spike after the voltage shift, c is higher than the threshold and the subinterval [a, c] is selected. Otherwise the subinterval [c, b] is selected. This procedure is repeated until the interval is sufficiently small (less than 10 − 4 ). Na m The ion channel x, the gating variable w, the activation and inactivation functions α w and β w , the maximal conductance g x , and the reversal potential E x are summarized. α p and β p are given by the equilibrium value p ∞ and the time constant

Reduction of a conductance-based neuron model
We have developed a reduction procedure from a detailed conductance-based model (Koch 1999;Izhikevich 2007) to an adaptive threshold model. We start from a conductancebased neuron model described by where I ion is an ionic current and w ! ¼ w 1 ; ⋯; w d ð Þis a vector of gating variables. Each gating variable w i is described by the kinetic equation (2).
The reduction consists of two approximations. First, we assume that a spike threshold θ V (t) is written as: where θ V ∞ is the spike threshold at the resting state and h ion describes the threshold modulation after a spike by an ionic current. If the voltage V(t) exceeds the threshold θ V (t), the neuron generates a spike. In addition, it is assumed that the previous spikes affect the spike threshold linearly. The validity of the assumption was tested by the comparison with the spike threshold of the detailed neuron model (data not shown). The effect of the spike waveform is incorporated into the reset rule. If the voltage exceeds the threshold, we shift the time and the voltage: where w sp is the spike width and δV is the voltage change during a spike. Specifically, the spike width w sp is approximately 2 ∼ 4 [ms] and the voltage change δV is −20 ∼ −10 [mV]. Second, we assumed that the ionic currents I ion are given by the sum of a spike-triggered ionic current η ion (t) and a leak current in the subthreshold regime (V < θ V ): where g ion , E ion are the average conductance and the reversal potential of an ionic current, respectively. By substituting Eq.
(11) into (9), we obtain where g tot ¼ ∑ ion g ion is the total conductance and E tot ¼ ∑ ion g ion E ion =g tot is the effective reversal potential. The formal solution of Eq. (12) can be written as, where τ m = C m /g tot is the effective membrane time constant. The Eq. (13) is a special case of the Spike Response Model (SRM) (Kistler et al. 1997;Gerstner and Kistler 2002;Jolivet et al. 2004). Here, the SRM is used to interpret the effect of the ionic currents on spike generation in the conductance-based model. Let us consider a new variable u that follows a linear equation without resetting after a spike, As the solution of Eq. (14) is Þe −s=τ m ds, the relationship between the new variable and the voltage is where the last term represents the voltage change during a spike. The spike threshold for u can be written as The effective threshold kernel H(t) is given by the voltage change during a spike, the spike-triggered ionic currents, and the spike threshold variation after a spike.
It should be noted that only the spike-triggered components of the ionic currents are considered in our framework. However, some of these currents, in particular I M , can be activated at voltages lower than the spike threshold (Prescott and Sejnowski 2008). Thus, the accuracy of the approximation may deteriorate if the voltage fluctuations are large.

Evaluation of the reduced model
We evaluated the reduced model by predicting the spike train of the detailed model neuron. The predictive performance was evaluated by injecting six fluctuating currents generated by the Ornstein − Uhlenbeck processes (Eq. (8)). The two types of input currents, i.e., the moderately noisy (σ = μ) input and the highly noisy (σ = 2μ) input, were examined. For each current type, three values of the mean μ were chosen so that the neuron generated spikes with 5, 10, and 20 [Hz]. The input parameters were (μ, σ) = (1.98, 1.98), (2.45, 2.45), Two input-output data sets {I(t), V(t)} were obtained by injecting two independent fluctuating currents for 50 [s], which were characterized by the same parameters (μ, σ, τ s ), into the detailed model.
The performance was evaluated based on the coincidence factor Γ (Kistler et al. 1997;Jolivet et al. 2004 where N d and N m are the number of spikes generated by the detailed model and by the reduced model, N c is the number of coincidences with precision Δ between the two spike trains, N c = 2νN d Δ is the expected number of coincidences using the Poisson process with the same rate ν with which the reduced model generates spikes. 3 Results

Typical behavior of the detailed conductance-based model
We first observed the behavior of a single-compartment conductance-based model with I Na , I Kd , I M , I Ca , and I AHP (Section 2.1). A rectangular current was injected into the three model neurons, i.e., the neuron with no adaptation (g M = g AHP = 0 [mS/cm 2 ]), the neuron with I M (g M = 0.1, g AHP = 0 [mS/cm 2 ]), and the neuron with The neuron with no adaption did not exhibit spikefrequency adaptation, i.e., the firing rate does not decrease during the stimulation (Fig. 2a). By contrast, the neuron with I M or I AHP exhibited spike-frequency adaptation, i.e., the firing rate dropped after the onset of the stimulation (Fig. 2b, c). The firing rate of the neuron with I AHP does not decrease gradually, because the neuron has the s-gate for I AHP . Consistent with previous studies (Benda and Herz 2003;Prescott and Sejnowski 2008), the slow K + currents induced spike-frequency adaptation. Due to the complexity and the nonlinearity of the detailed model, it is not clear how slow K + currents regulate spike generation of a neuron. Thus, we investigated the effects of slow K + currents by mapping the detailed neuron model to a simplified model, and derived a reduced model that clarifies how slow K + currents modulate the effective spike threshold.

Spike triggered ionic current: η ion (t)
A constant current with a pulse (Eq. (7)) was injected into the neuron with I M and the neuron with I AHP , and the spiketriggered ionic currents η ion (t) were calculated. Because Na + , K + , and Ca 2+ currents vanish within a brief period immediately after the spike (typically 4 [ms] after the spike onset), we focused on analyzing slow K + currents, I M and I AHP (Fig. 3a,  b).
First, we examined the spike-triggered current induced by I M , η M (t). By replacing an action potential with a rectangular pulse, similar to the approach of Destexhe (1997), the spiketriggered current can be approximated by the exponential function (Appendix A), where τ p v ð Þ is the p-gate time constant and v is an average voltage after a spike. The formula (19) is in agreement with η M (t) obtained from the detailed neuron model for various values of I M parameters (g M , τ max ) and membrane depolarization V c ( Fig. 3c and data not shown for τ max and V c ). There is a slight discrepancy in η M (t) between the detailed model and Eq. (19) for small t, which may be due to the spike waveform. A more accurate formula can be obtained by incorporating this effect (Appendix A).
Second, we examined the spike-triggered current induced by I AHP , η AHP (t). By replacing the calcium current   Table 1 with an impulse, the spike-triggered current can be approximated by the sum of two exponentials (Appendix A).
where τ Ca is the Ca 2+ time constant,τ s ¼ β −1 s is an approximation of the s-gate time constant, and β s is the inactivation rate of the s-gate. The formula (20) is also in agreement with η AHP (t) obtained from the detailed neuron model for various values of I AHP parameters (g AHP , β s , and τ Ca ) and the membrane depolarization V c (Fig. 3c and data not shown for β s , τ Ca and V c ).

Spike threshold variation by an ionic current: h ion (t)
A constant current with a pulse (Eq. (7)) was injected into the detailed model neurons and the instantaneous spike threshold was calculated (Section 2.3). Again, three neurons were examined, i.e., the neuron with no adaptation, the neuron with I M , and the neuron with I AHP . Whereas the spike threshold decays rapidly after a spike in the neuron with no adaptation, it decays slowly in the neuron with I M or I AHP (Fig. 4b). We can thus conclude that the threshold variation after a spike is mainly caused by the slow K + currents.
The spike threshold variation induced by I M was evaluated by comparing the spike threshold in the neuron with I M to that with no adaptation. The spike threshold variation is approximately proportional to the spike triggered current η M (t) (Appendix B), where the weight b M is proportional to a M in Eq. (19). Equation (21) can accurately describe h M (t) for various values of the I M parameters (g M , τ max ) and of the membrane depolarization V c (Fig. 4c and data not shown for τ max and V c ). Next, the spike threshold variation induced by I AHP was evaluated by comparing the threshold in the neuron with I AHP to that without I AHP . The spike threshold variation is approximately proportional to the spike triggered current η AHP (t) (Appendix B), where the weight b AHP is proportional to a AHP in Eq.
(20). Equation (22) can accurately describe h AHP (t) for various values of the I AHP parameters (g AHP , β s , and τ Ca ) and the membrane depolarization V c (Fig. 4c and data not shown for β s , τ Ca and V c ).  3.4 Reduction of the detailed conductance-based neuron model The conductance-based neuron model can be reduced to an adaptive threshold model (Section 2.4), where τ m is the membrane time constant and θ u is the spike threshold for u (effective spike threshold) written as t k is the k-th spike time, and H u (t) is the effective threshold kernel that describes how the effective spike threshold changes after a spike. We investigated the effect of the slow K + current parameters on the effective threshold kernel. The threshold kernel H u (t) of the neuron with I M can be described by the sum of two exponentials, The threshold kernel is always a monotonically decreasing function in the neuron with I M (Fig. 5a). We can derive a formula that clarifies the relationship between the slow weight α M and I M parameters (Appendix C), where δp is the changes in the p-gate variable during a spike. As predicted by Eq. (26), the slow weight α M increases as g M increases, and decreases as τ max increases (Fig. 6a). Numerical results indicate that I M parameters does not affect on the fast weight α 0 significantly (Fig. 6a). The threshold kernel H u (t) of the neuron with I AHP is described by the sum of three exponentials, Interestingly, the threshold kernel can be a nonmonotonic function in the neuron with I AHP , and a hump was observed in H u (t) (Fig. 5b). We can also derive a formula that clarifies the relation between the slow weight α AHP and I AHP parameters (Appendix C), where δCa is the changes in Ca 2+ concentration during a spike. As predicted by Eq. (28), the slow weight α AHP increases as g AHP increases, and decreases as β s or τ Ca increases (Fig. 6b). Numerical results indicate that I AHP parameters does not affect on the fast weight α 0 significantly (Fig. 6b).  Each parameter (a: g M and τ max ; b: g AHP , β s , and τ Ca ) was tested at three levels. Other parameter values were given in Table 1 3.5 Validation of the reduced model , α 0 , α AHP } for the neuron with I AHP , were determined by maximizing the coincidence factor Γ (Section 2.5) using the simplex downhill method ). Then, the predictive performance was evaluated by calculating the coincidence factor from the other data set (test data) that was not used for parameter optimization. We found that the reduced model can accurately predict spike trains of the detailed model (Fig. 7). The predictive performance Γ for the input currents was 0.854 ± 0.01 (means ± standard errors, unless stated otherwise) for the neuron with I M , and 0.903 ± 0.01 for the neuron with I AHP , and the results are summarized in Table 2

Coding property of the reduced model
We analyzed the impact of slow K + currents on the coding property of a neuron using the reduced model. Here, we focused on the effect of the input noise on firing irregularity. First, we considered that the reduced model neuron (Eqs. (23) and (24)) is stimulated by a constant current, I ex (t) = I 0 . An asymptotic periodic solution with period T is written as where t f is the most recent spike time, and η T ∞ (t) describes the threshold variation between the spikes; η T ∞ (t) for the neuron with I M is and η T ∞ (t) for the neuron with I AHP is The spike condition at the next spike, t = t f + T, leads to We can analytically evaluate the firing rate f = T − 1 by solving Eq. (32), and the analytical results are in agreement with f-I curves calculated from simulated spike trains (Fig. 8a). The f-I curve of the neuron without adaptation (g M = g AHP = 0) can be explicitly written as which is similar to the f-I curve of the LIF neuron. Note that the response of the reduced model with I M to the constant current (Eqs. (29) and (30)) is equivalent to the response of the time-dependent threshold model (Tuckwell 1978;Lindner and Longtin 2005;Tamborrino 2016). Next, we examined the effect of the input noise on spiking irregularity. We have not been able to derive a full analytical result for this effect; however, it is possible to predict the effect of the input noise with the following argument. Let us consider a situation in which a neuron is stimulated by the constant current before the N-th spike (N ≫ 1) and stimulated by the constant current with small noise after the N-th spike. We can evaluate how the input noise changes the timing of the  The performance of spike prediction for the detailed neuron models using the reduced model is summarized. Each neuron was injected with six fluctuating input currents. Current BM^denotes moderately noisy input (σ = μ) and current BH^denotes highly noisy input (σ = 2μ) subsequent spike. At the (N + 1)-th spike time, the threshold should cross the potential where t N is the N-th spike time, δT and δu are perturbations due to the small noise. By Taylor-expanding θ u assuming that δT is small, we obtain Equation (35) indicates that slow K + currents improve the robustness against noise in a different manner (Fig. 8b), i.e., I M suppresses the spike interval dispersion δT/T for a broad firing range, whereas I AHP suppress the dispersion only at a low firing range (~3 [Hz]). The dispersion δT/T is not identical to the coefficient of variation (Cv) of interspike intervals (ISIs), however it has a close relation to Cv. We found that this differential effect was also observed in Cv of the detailed model with slow K + currents (Fig. 8c). Finally, we examined how the slow K + currents modulate autocorrelation of a spike train that was defined as ρ 1 = 〈ISI i ISI i + 1 − 〈ISI i 〉 2 〉/〈ISI i 2 − 〈ISI i 〉 2 〉, where ISI i is the i-th ISI and 〈 … 〉 is the averaging over index i. The autocorrelation quantifies how often a long ISI is followed by a short ISI and vice versa. The neuron model with the slow K + currents can reproduce the negative ISI correlation, which was commonly observed in sensory periphery and cortical neurons (Farkhooi et al. 2009). As shown in Fig. 8d, the effect of I AHP on the autocorrelation is stronger than that of I M in the low firing rate regime (<15 [Hz]), whereas the effect of I AHP is similar to that of I M in the high firing rate regime (>15 [Hz]). A previous work  Table 1. d: Effect of the slow K + currents on autocorrelation. The autocorrelation was calculated from spike trains generated from the reduced model. The parameters were same as a (Chacron et al. 2001) showed that the negative ISI correlation can improve the capacity for encoding time-varying stimulus. Our result implies that the slow K + currents improve the encoding of time-varying stimulus in a different way.

Discussion
We have shown that the detailed conductance-based neuron model with slow K + currents (I M and I AHP ) can be reduced to an adaptive threshold model. The reduced model is written as where θ u is the spike threshold for u, and H u (t) is the threshold kernel. We have also derived formulae that describe the relationship between slow K + current parameters and reduced model parameters (Eqs. (26) and (28)), which provide a physiological interpretation of the reduced model. The reduced model can accurately predict spike trains of the detailed model (Fig. 7). Our analysis of the reduced model revealed that slow K + currents have differential effects on noise tolerance of a neuron, i.e., I M suppresses firing irregularity regardless of the firing rate, whereas I AHP suppresses firing irregularity only at a low firing range (Fig. 8b, c). The slow K + currents induce negative interspike interval correlations, and the effect of I AHP is stronger than that of I M in the low firing regime (Fig. 8d).

Mapping a detailed conductance-based neuron model to a simplified model
As noted in the Introduction, one approach of obtaining a reduced model is to develop a mathematical framework from detailed neuron models to simplified models. This approach has clarified the relationship between these models. For example, the FitzHugh-Nagumo model was derived from the Hodgkin-Huxley model by assuming that Na + activation (m) is instantaneous and that Na + inactivation (h) and K + activation (n) change with a similar time constant (Abbott and Kepler 1990;Rinzel and Ermentrout 1998). A generalized integrate-and-fire model can also be derived from the Hodgkin-Huxley model by linearization (Destexhe 1997;Koch 1999;Richardson et al. 2003).
In this study, we have extended the linearization approach by including the spike history effect, which is essential for describing the effect of slow K + currents on spike generation. The linearized model is a simple linear equation with the effective threshold θ u (t) (Eq. (36)) that incorporates the effect of ionic currents and spike threshold variation on neuronal excitability. We have shown that the effective threshold obtained from the detailed model with slow K + current can be approximated by a modified multi-timescale adaptive threshold (MAT) model ).

Reduced neuron model
Spike-frequency adaptation can be described by simplified models with adaptation, which is modeled by adaptive current (Liu and Wang 2001;Brette and Gerstner 2005;Izhikevich 2007) or adaptive threshold (Chacron et al. 2000(Chacron et al. , 2007Liu and Wang 2001;Jolivet et al. 2004Jolivet et al. , 2006Jolivet et al. , 2008. The adaptive threshold models can reproduce the interspike interval statistics (Chacron et al. 2000), f-I curve (Rauch et al. 2003;Kobayashi 2009), and spike timings (Jolivet et al. 2006(Jolivet et al. , 2008 of a neuron recorded in experiments. On the other hand, the adaptive threshold model was criticized because, unlike the adaptive current model, it cannot reproduce the lateral shift of f-I curves observed in experiments (Benda et al. 2010). Note that the derived MAT model (Eq. (36)) incorporates both effects, i.e., the effect of the adaptive current and threshold. This fact can explain the success of the MAT model in accurately predicting spike times Yamauchi et al. 2011).
The derived model has two advantages. First, the model is essentially linear; the linearity makes mathematical analysis tractable. Indeed, it is possible to examine the effect of noise on firing irregularity, which can predict a qualitative behavior of the detailed model (Fig. 8). In addition, the linearity enables us to efficiently simulate a network of neurons by the exact sub-threshold integration (Morrison et al. 2007;Yamauchi et al. 2011). Second, the reduced model offers a clear relationship between the slow K + parameters and reduced model parameters (Eqs. (26) and (28)). This relationship is important because it enables us to analyze the effect of slow K + currents using the reduced model.

Spike threshold variation in experiments
Conventionally, it was considered that a neuron has a fixed voltage threshold for generating an action potential. However, experimental studies in vivo have suggested that the spike threshold is not constant but is highly variable (Azouz and Gray 2000;Henze and Buzsaki 2001;Chacron et al. 2007). Studies in the rodent hippocampus (Henze and Buzsaki 2001) and fish (Chacron et al. 2007) have demonstrated that the spike threshold increases after each action potential, which is referred to as Bthreshold fatigue.^We found that the spike threshold of the detailed neuron model jumps and decays exponentially after each spike with a time constant of~100 [ms] (Fig. 4), suggesting that slow K + currents may be the possible cellular mechanism underlying threshold fatigue. Other biophysical mechanisms, particularly Na + currents, may underlie the threshold variability. The spike threshold also varies with the voltage derivative preceding a spike (Azouz and Gray 2000). Interestingly, it was shown that Na + inactivation modulates the spike threshold, which varies with the membrane voltage with a small time constant (τ h (v) ≈ 2~10 [ms]) (Platkiewicz and Brette 2010;Fontaine et al. 2014) and Na + inactivation can explain the voltagedependence of the spike threshold observed in experiments (Platkiewicz and Brette 2011). The modulation of the spike threshold differs depending on its factor, in other words, the threshold modulation by slow K + currents is slow and accumulative, whereas that by Na + inactivation is rapid.
The instantaneous spike threshold of layer-5 pyramidal neurons has been estimated from the membrane potential recorded in vitro (Badel et al. 2008). The results suggest that the threshold modulation after a spike is described by the sum of two exponentials, where θ V (t) is the spike threshold and t f is the most recent spike time. The fast components were A 1 ≈ 10 [mV] and τ 1 ≈ 20 [ms], whereas the slow components were A 2 ≈ 1 [mV] and τ 2 ≈ 100 [ms]. The detailed model used here reproduces the slow component in the threshold modulation; however, it does not reproduce the fast component. This is presumably due to the difference in Na + current kinetics that describes the shape of an action potential. Indeed, it was reported that the spike waveform recorded from experiments is much shaper than that of Hodgkin-Huxley models (Badel et al. 2008).

Functional implications of the slow K + currents
It is well known that slow K + currents induce the spikefrequency adaptation, which acts as a spike-triggered self-inhibition ( Fig. 2; Benda and Herz 2003;Prescott and Sejnowski 2008). Several studies have proposed functional consequences of spike-frequency adaptation. For instance, the adaptation generates the Bforward masking^effect, which suppresses the neuronal response under a prolonged stimulus (Liu and Wang 2001), improve signal transmission for low frequency stimulus (Chacron et al. 2007), and contributes to sparse and reliable coding (Farkhooi et al. 2013). Here, we have derived a simplified model that can reproduce the differential effects of slow K + currents. The reduced model can accurately predict spike trains of the detailed neuron model (Fig. 7) and reproduce the f-I curve and spike train power spectrum (Data not shown).
Previous studies have suggested that slow K + currents have differential effects on the coding property of a single neuron. For instance, I M facilitates spike-timing coding because it improves the robustness of spike pattern against the input noise. In contrast, I AHP enhances spike-rate coding, because it regularizes the spike train elicited by slow inputs (Prescott and Sejnowski 2008). It has also been suggested that I M increases, whereas I AHP decreases, the response to low-frequency input signals (Deemyad et al. 2012). Our analysis revealed a new differential effect underlying slow K + currents (Fig. 8b, c), i.e., I M suppresses firing irregularity regardless of the firing rate, whereas I AHP suppresses the irregularity only at a low firing range (~3 [Hz]). This result suggests that neurons with I AHP can contribute to the generation of rhythmical activity at a low firing rate. We hope that the reduced model will be useful for analyzing how the slow K + currents impact on the coding properties of single neurons and neural populations. membrane time constant is much smaller than the slow K + time constants: τ m ≪τ p v ð Þ; τ Ca ;τ s , the formula can be simplified further, where the weights are α 0 ¼ −δV e w sp =τ m −a M =g tot , α M = a M / g tot + b M , α AHP = a AHP /g tot + b AHP . The slow weights can be related to the slow K + parameters by using Eqs. (40), (46), , where δp :¼p w sp À Á −p ∞ v ð Þ and δCa are the changes in p and [Ca 2+ ] during a spike.
Open Access This article is distributed under the terms of the Creative Comm ons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.