Astrocytic Kir4.1 channels and gap junctions account for spontaneous epileptic seizure

Experimental recordings in hippocampal slices indicate that astrocytic dysfunction may cause neuronal hyper-excitation or seizures. Considering that astrocytes play important roles in mediating local uptake and spatial buffering of K+ in the extracellular space of the cortical circuit, we constructed a novel model of an astrocyte-neuron network module consisting of a single compartment neuron and 4 surrounding connected astrocytes and including extracellular potassium dynamics. Next, we developed a new model function for the astrocyte gap junctions, connecting two astrocyte-neuron network modules. The function form and parameters of the gap junction were based on nonlinear regression fitting of a set of experimental data published in previous studies. Moreover, we have created numerical simulations using the above single astrocyte-neuron network module and the coupled astrocyte-neuron network modules. Our model validates previous experimental observations that both Kir4.1 channels and gap junctions play important roles in regulating the concentration of extracellular potassium. In addition, we also observe that changes in Kir4.1 channel conductance and gap junction strength induce spontaneous epileptic activity in the absence of external stimuli.

Introduction Temporal lobe epilepsy, which has a specific clinical presentation of seizures arising from the hippocampus [1], is a serious health risk to affected individuals. Extensive experimental and theoretical studies have shown that elevated K + concentration in the extracellular microenvironment ([K + ] o ) may be linked to spontaneous epileptic seizure activity in the absence of external stimuli [2][3][4][5][6]. Astrocyte-mediated K + buffering generally serves to maintain [K + ] o , supporting normal neuronal electrical activities [7][8][9][10]. Mechanisms for K + removal, regulated by astrocytes, can be broadly categorized as K + uptake and K + spatial buffering [11,12]. In the case of K + uptake, excess K + is absorbed primarily by astrocytes through Na + /K + -ATPase pumps and inwardly rectifying Kir4.1 channels or through Na + -K + -Cl cotransporters. K + is spatially buffered via gap junctions, which redistribute excess K + taken up from astrocytes in areas of excessive neuronal activity to astrocytes with lower K + concentrations within the astrocytic network. Experimental data suggest that Kir4.1 channels play a prominent role in local uptake of extracellular K + [8,13,14], and many experimental studies have observed that dysfunction of local uptake by astrocytes or inactivation of astrocytic gap junction protein expression cause the generation or spread of seizure activity [15][16][17][18][19][20]. It was recently reported that down-regulation or dysfunction of K + uptake channels (Kir4.1) in astrocytes induces tonic-clonic seizures [7]. For example, low-intensity Schaffer collateral stimulation generates epileptiform activity in gap junction-deficient mice. Another experiment found that Tsc1 (TSC: Tuberous Sclerosis Complex) inactivation in astrocytes caused defects in astrocytic gap junction coupling and potassium clearance, leading to epilepsy in Tsc1GFAPCKO mice [21]. Furthermore, expression of the astrocytic gap junction proteins connexin 43 (Cx43) and connexin 30 (Cx30) is altered in epilepsy, and changes in gap junction communication have been observed in sclerotic hippocampal tissue of epileptic animal models [22]. Therefore, we investigated the dynamic mechanisms by which dysfunction of K + uptake or gap junctions leads to changes in extracellular K + concentration and resultant pathological epileptic discharges using a theoretical study with a computational astrocytic-neural network model.
Recent studies have applied theoretical analysis approaches to explore the relationship between astrocytic modulation of [K + ] o and the initiation or maintenance of epileptic seizure activity [23][24][25][26][27][28][29][30][31][32]. Cressman J. R. et al. provided a simplified astrocytic K + uptake function that incorporates both passive and active local K + uptake into a single sigmoidal response function that depends solely on extracellular K + concentration [23][24][25][26][27]. Hübel studied the mechanism of extracellular K + uptake in astrocytes using a phenomenological equation for the rate of K + uptake, which depends entirely on extracellular K + concentration [30]. Øyehaug L et al. modeled the mechanism of K+ uptake from the perspective of volume dynamics [28,29], modeling the cotransporters (Na + /K + /2Cl and K + /Clcotransporters) and the Na + /K + -ATPase pump of astrocytes. However, these model systems did not adequately reflect the complexity of astrocytes. For example, in addition to cotransporters and the Na + /K + -ATPase pump, the inwardly rectifying K + channels (Kir4.1) in astrocyte membranes and inter-astrocyte spatial buffering also contribute to extracellular K + removal. Sibille J et al. developed mathematical models of Kir4.1 channels that depend on both the intra-astrocyte K + concentration and the astrocytic membrane potential [31,32]. In addition, they verified the prominent role of Kir4.1 channels in regulating [K + ] o during repetitive stimulation using their mathematical model.
It is worth mentioning that none of the aforementioned studies considered the relationship between Kir4.1 channel conductance and the occurrence of spontaneous epileptic seizures. Morphologically, previous studies of the extracellular K + local buffer mechanisms in astrocyteneuron networks have shown that astrocytes have a 1:1 quantitative relationship with neurons; however, experimental data show that one cortical neuron is typically surrounded by multiple astrocytes [20,33,34]. In support of this, in 2009, Azevedo et al. reported that the ratio of neurons to glia is approximately 1:4 in the human cerebrum [33].
In light of these findings, we constructed a model of an astrocyte-neuron network consisting of a single compartment neuron connected to four surrounding astrocytes with Kir4.1 channels and Na + /K + -ATPase pumps. Extracellular potassium ions diffused in and out of the space between the neuron and astrocytes. In addition, we developed a novel configuration for astrocytic gap junctions, connecting two astrocyte-neuron network modules. The model function parameters of the gap junction are based on a nonlinear regression fit of a set of experimental data previously published [17]. The model simulation results were validated by comparison with published experimental data [17]. We conducted a series of simulations using the single astrocyte-neuron network models and the coupled astrocyte-neuron module network and showed that changes in Kir4.1 channel conductance and gap junction strength induce spontaneous epileptic activity in the absence of external stimuli.

Astrocytic K + uptake by Kir4.1 channels
To validate the role of Kir4.1 channels in regulating astrocytic and neuronal extracellular K + concentration, a model was constructed representing astrocyte-neuron network modules consisting of a single compartment neuron connected to four surrounding astrocytes with Kir4.1 channels and Na + /K + -ATPase pumps (Fig 1(a)). Based on evidence from previous studies [23,25,26,30,31], there exist several key factors mediating K + concentration in the extracellular space: K + currents across the neuronal membrane, Na + /K + -ATPase pumps for neurons and astrocytes, Kir4.1 channels in astrocytes, spatial diffusion of K + in the extracellular space, and others (illustrated in Fig 1(a)). During an action potential, K + is released by neurons at the extracellular space, followed by K + flowing into adjacent astrocytes via the Na + /K + -ATPase pump and Kir4.1 channels. It is important to note that the Kir4.1 current is larger in magnitude than the outward Na + current from the Na + /K + -ATPase pump, resulting in overall depolarization of the astrocytic membrane (Fig 1(b)).
Effects of Kir4.1 channel blockage on extracellular K + regulation were first tested by setting the Kir4.1 conductance to 45.0 pS for the normal (control) condition. Sine stimulus inputs (the stimulus amplitude and frequency are 1μA/cm 2 and 10 Hz, respectively; here, 10 Hz was used to verify our results in comparison with previous experimental recordings [35]) with different durations, e.g., 30 s, 10 s and 100 ms, were applied to the neuron when Kir4.1 channels were in the control state. Note here that a 100 ms stimulus resulted in a small undershoot and a very short recovery to the [K + ] o baseline. As the duration increased, both the undershoot area and recovery time increased. Compared with 100 ms, the 10 s and 30 s stimulations required significantly longer times to recover to the baseline, as shown in Fig 1(c). This result is consistent with an experimental trend observed in vivo in the mouse hippocampus [35].
We next tested our model with Kir4.1 channel inhibition (g kirA = 0.01 pS). We observed that extracellular K + concentration displayed a more pronounced undershoot after long-  [36]. For other parameter values, see Table 1. In addition, ε = 1.2/s -1 and K bath = 4 mM. lasting sine stimulation (10 Hz, 10 s). The maximum amplitude of extracellular K + concentration for g kirA = 0.01 pS was higher than that for g kirA = 45.0 pS during a 10 s external stimulation input (see lower panel in Fig 1(d)). During the 10 s post-stimulus recovery period, extracellular K + initially dropped to below the baseline levels before returning gradually to the resting-state equilibrium concentration, taking even longer to recover to the K + baseline concentration, as shown in the lower panel of Fig 1(d). Furthermore, it was observed that in response to stimulus under control conditions, astrocytes responded with a rapid increase in intra-astrocyte K + concentration with an accompanying high amplitude, while the baseline K + concentration was low. On the contrary, the K + concentration in astrocytes rose much more slowly to a lower peak concentration when Kir4.1 was blocked, as shown in the top panel of Fig 1(d). These computational results were all in agreement with previously recorded experimental results [35,36].
Additionally, the critical value of g kir was examined for the generation of spontaneous epileptic seizures in the absence of external stimuli. Model simulations suggested that spontaneous periodic epileptic activity is induced when g kir < = 7.0 pS without external stimulation ( shows that spontaneous periodic epileptic firing is accompanied by higher amplitude  oscillations of extracellular K + concentrations and lower peak intra-astrocyte K + concentrations compared with normal neuronal discharges. Fig 3 presents a new model in which astrocyte gap junctions connect two astrocyte-neuron network modules. This network model was constructed to examine the K + spatial buffering mechanisms dominated by astrocytic gap junctions. In this model, individual neurons are surrounded by 4 astrocytes, and they share the same bath K + concentration. During action potential generation, neurons release K + into the extracellular space, and K + enters adjacent astrocytes via the Na + /K + -ATPase pump and Kir4.1 channels. Subsequently, in astrocytes with higher concentrations, K + may flow into adjacent astrocytes via inter-astrocytic gap junctions due to the electrochemical driving force.

Gap junction deficiency causes a delay in extracellular K + buffering
To validate whether this type of inter-astrocytic gap junction behaves the same way as those in experimental observations, the decay dynamics of extracellular K + concentration ([K + ] o ) was examined in response to the sine stimulus trains (10 s, 20 Hz). We chose 20 Hz rather than the 10 Hz used above in order to compare the results in our model with previously published experimental results ( Fig. 5B in [17]). The results demonstrated that after the stimulus trains, [K + ] o initially decayed very rapidly, followed by a prolonged phase of slower decay, often associated with an undershoot [37] (Fig 4(a)). To quantify the faster K + decay phase, the time at which the [K + ] o amplitude has decayed to 1/e of its initial maximal amplitude is defined as the decay time constant (t 1/e ), providing a descriptive parameter for the decay process. Anke Wallraff [17] observed an inverse correlation between t 1/e and [K + ] o , which is well fitted by a power-law function (Fig 4(b)).   Fig. 5B from [17]). Black spots and gray spots are experimental data obtained from connexin30 -/rats and control gap junction rats, respectively. (c) The relationship between extracellular K + concentration and decay factor t 1/e in control conditions or the gap junction-deficient state. The dashed lines and solid lines denote experimental fitted curves (see Fig. 5B from [17]) and fitted curves from numerical simulation data, respectively. for deficient gap junctions obtained from rat hippocampal slice recordings. Clearly, deficient gap junctions are characterized by a larger decay time constant, suggesting slower removal of K + load [17]. Here, the proposed model simulated this process by applying the sine stimulus trains (10 s, 20 Hz, the stimulus amplitude is randomly chosen from 1 to 60μA/cm 2 to induce different spiking patterns and potassium diffusion dynamics) to neuron 1 (N1). The K + concentration K bath , the spatial diffusion coefficient ε and the Na + /K + -ATP pump strength ρ were set to different values to study their effect on the correlations between the diffusion time constant t 1/e and extracellular K + concentration, as shown in Fig 4c. The Kir4.1 channel conductance was set to 45.0 pS for Fig 1c (for values of the other parameters, see Table 1). By setting the strength of the gap junction to F = 20.0 and 0.01, the model simulation reproduced the wild-type control (Fig 4c) and dko deficient (Fig 4d) gap junctions for neuron 1 (N1). We then used the nonlinear least-squares method to fit the relation curve of t 1/e and [K + ] o that was obtained by the coupled astrocyte-neuron network modules. The fitted curve function for normal gap junctions is t 1/e = 3.8179exp(-0.0623[K + ] o ) (left panel, Fig 4(c)). The fitted curve function for deficient gap junctions is t 1/e = 4.2055exp(-0.0569[K + ] o ) (right panel, Fig 4(c)). The simulated fitted curve functions are quantitatively similar to experimental observations. In addition, to validate the physiological significance of the model constructed in this section and the reliability of our numerical results, statistical analysis was conducted on the error between the experimental and numerical fits for different gap junction strengths, as shown in Fig 4(d). Here, the fit error is the mean absolute error (calculations shown in S1 Text). The fit error between the numerical and experimental data is far less than 10%. Considering the experimental accuracy and other environmental factors, we believe that the gap junctionmediated coupled model presented herein is within the permissible error range. In addition, different strengths of gap junctions were applied between astrocytes, e.g., F = 20 and 0.01 (Fig  4(e)). The results show that decreased gap junction strength induces higher [K + ] o peaks, and if the strength of gap junctions experiences a strong decrease, [K + ] o will have a larger undershoot and an increased baseline recovery time.
The above simulations applied the exponential diffusion function observed in experimental data reported by Anke Wallraff [17]. While this is a reasonable method to validate our overall model, it does not necessarily justify modeling gap junction diffusion as an exponential function of intracellular K + . Therefore, we performed a new set of numerical simulation experiments to address this concern. We chose several different types of diffusion functions for our model network. Numerous experiments have observed that the coupling-deficient astrocytic network (with depression of connexin43 expression or connexin30 expression) in the nervous system is related to generation of spontaneous seizures [17,22,[38][39][40]. To verify the relationship between the dysfunction of the gap junctions and spontaneous epileptic discharges in the absence of external stimuli, we used the coupled two astrocyte-neuron network modules. Model simulations suggested that spontaneous epileptic activity is induced when F< = 10.0 pS without external stimulation (Fig 5(a)). Fig 5(a) shows the time sequences of the neuronal membrane potentials (V N1 ) and the K + concentration in the extracellular space ([K + ] o1 ) for F = 20.0, 8.0 and 0.001. The neuron (N1) exhibited normal spontaneous slow oscillation activity when the gap junction strength F was equal to 20.0. However, the activity changed into spontaneous epileptic seizures after certain time periods when the value of F decreased. For example, the neuron (N1) changed into spontaneous seizures firing at 50 s when F was 8.0. Specifically, neurons rapidly enter spontaneous epileptic discharge with smaller gap junction strengths. As such, the neuron (N1) exhibited spontaneous seizure firing at 50 s when F was 8.0, while entering spontaneous seizure firing at 24 s when F was 0.001 (Fig 5(a)). This type of spontaneous epileptic seizure has been observed in various experimental recordings [2,41,42]. Fig 5(b) summarizes the schema of astrocytic gap junction strength F and [K + ] o , or t 1/e when there is no external stimuli to the neurons, the Kir4.1 channel conductance is 45.0 pS and other parameters are fixed, as in Table 1. [K + ] o and t 1/e for N1 gradually decrease as F increases. These relationships can be used to explain the extracellular K + concentration requiring more time to recover to the baseline when the gap junction is deficient. The above results confirm that decreasing astrocytic gap junction strength induces spontaneous epileptic seizures in the astrocyte-neuron network in the absence of external stimulation. Furthermore, neurons will rapidly enter into spontaneous epileptic seizures with smaller gap junction strength.

Discussion
Potassium clearance is widely accepted as a primary function of astrocytes. Several computational models have investigated the buffering mechanism for extracellular K + accompanied by neural activities [23,25,26,28,30,31,43]. A model accounting for K + concentration in the extracellular space and astrocytic compartments has been used to quantify involved astrocytic ion channels and transporters (Na + /K + -ATPase, NKCC1, etc.) [28,29]. While NKCCl has been recently shown to not be involved in activity-dependent K + clearance in hippocampal slices [15], astrocytic Kir4.1 channels are crucial for the recovery of basal extracellular K + levels and neuronal excitability during external stimulation [31]. To unravel how the acute role of astrocytes in [K + ] o homeostasis induces spontaneous epileptic discharge and maintains normal electrophysiological activity in the astrocyte-neuron system, we constructed an astrocyte-neuron network module consisting of a single compartment neuron connected to four surrounding astrocytes with Kir4.1 channels and Na + /K + -ATPase pumps. Extracellular potassium ions were allowed to diffuse in and out of the space between the neuron and astrocytes. In addition, we have built a new model for astrocytic gap junctions, enabling connection of two astrocyte-  neuron network modules. The model function parameters were based on a nonlinear regression fit based on a set of experimental data published previously [17]. The model simulation results were validated by comparison with the published results [17]. We first verified that spontaneous periodic seizures can be induced directly in the absence of stimuli when the astrocytic Kir4.1 channel conductance or the strength of gap junction-mediated K + spatial buffering decreases to a certain level.
We first simulated the specific dynamic characteristics of astrocytic Kir4.1 channel conductance inducing spontaneous epileptic activity in the absence of external stimulus input. To mathematically examine K + dynamics in neurons, astrocytes and extracellular spaces, we constructed an astrocyte-neuron network module consisting of a single compartment neuron connected to 4 surrounding astrocytes. The results demonstrated that the K + undershoot of the extracellular space increases after a period of stimulus in Kir4.1 knockouts [15,35,44]. We verified these results using the proposed computational model compared to controls and found that K + concentration in the extracellular space exhibits an increased undershoot after a 10 s sine stimulus sequence, taking a longer time to return to the baseline in Kir4.1 channel-deficient conditions. Moreover, the baseline K + concentration in astrocytes is higher and rises more slowly to a decreased peak value after an external stimulus train compared to controls. The high baseline K + concentration in astrocytes is due to the normal function of the Na + /K + -ATPase pump. In particular, spontaneous epileptic activity was directly induced in the single astrocyte-neuron network module without external stimuli when the conductance of Kir4.1 channels was lower than a certain threshold, further demonstrating that down-regulation of astrocytic Kir4.1 channels is closely related to neuropathological hyper-excitability [8,45,46]. At the same time, our single astrocyte-neuron network correctly predicts that when there is a partial blockage of Kir4.1 channels in astrocytes, rapid instantaneous neuronal seizure activity is induced. Moreover, it is possible that astrocytic Kir4.1 channels modulate extracellular K + as a means of regulating synaptic activity.
The dynamic characteristics of astrocytic gap junction strength pertaining to spontaneous epileptic activity in the absence of external stimuli were first simulated in this study. Gap junctions in astrocytes appear to play a dual role: on the one hand, they counteract the generation of hyperactivity by facilitating the decrease of extracellular K + levels; on the other hand, they constitute a pathway for energetic substrate delivery that fuels neuronal hyper-excitability [16]. Thus far, research involving the effects of astrocytic networks coupled with gap junctions on extracellular K + spatial buffering is still in the biological experimentation stage. Previous studies have examined the role of astrocytic gap junctions in potassium buffering with respect to epileptic seizures [17,22]. We presented a new exponential function model for astrocyte gap junctions that connect two astrocyte-neuron network modules. To verify the validity of the proposed model of gap junctions, we compared two groups of data of the fitted exponential function curves of t 1/e vs [K + ] o obtained from numerical simulation data based on the coupled astrocytic-neural network modules and previous experimental recordings [17]. Both control conditions and deficient gap junction cases were compared while sine stimulus trains were applied to neuron 1 (N1). We discovered that the curve fitting error of numerical and experimental data was within a reasonable range. In addition, the numerical results validated the finding that gap junction-mediated K + spatial buffering is slower in the absence of astrocytic gap junctions compared to wild-type rats. At the same time, the neuron entered a spontaneous epileptic state more quickly at weaker gap junction strengths. Finally, we summarized the schema of the astrocytic gap junction strength F and [K + ] o or t 1/e when there is no external stimulation to either of the neurons, and this relationship can be used to explain the observation that extracellular K + requires more time to return to the baseline when the gap junction is dysfunctional. Indeed, our results reproduced experimental observations wherein expression-deficient astrocytic gap junction proteins (Cx43 or Cx30) distorted neuronal information processing and the generation of spontaneous epileptiform events [16,47]. Our coupled astrocyteneuron network modules predict that with partially blocked astrocytic gap junctions, neuronal seizure activity will gradually develop after a delay of tens of seconds. This finding could represent a protective mechanism to prevent K + accumulation from reaching neurotoxic levels. These results can be used as a basis for performing further analysis of the characteristics of sodium, calcium, the power factor and other factors in the network, as well as the mutual influences among them.
In this study, we considered the gap junction between one of two adjacent astrocytes. In vivo many glial cells surrounding the firing neuron will buffer potassium into adjacent glial cells that have lower concentrations. We have chosen several other different diffusion functions (including an exponential function, a linear function and a threshold-nonlinear function) to simulate experimental data (Fig 4c and S1 Fig). The simulation results suggest that the dynamic properties of potassium ion channels faithfully follow an exponential decay function. Therefore, we chose the exponential function for our model and all other simulations. Moreover, we examined several new model network configurations to test inverse correlations between t 1/e and [K + ] o . First, we assessed gap junctions between astrocytes in the same neuron-astrocyte module (S2A Fig). In the simulation, we reduced the gap junction weight to 6.5 in order to match previous experimental parameters. The simulation results of our novel model (S2C Fig) are indeed quite similar to our original model with only one gap junction in a two module network (S2B Fig). Therefore, either there are more gap junctions with a weak weight between glial cells or fewer gap junctions with a strong gap junction weight, and the network model reproduced established experimental results. Choosing 1:4 as the ratio of neurons to glial cells might be closer to actual physiology [33]. Second, we scaled up the network to a relatively large network with more neurons and glia, as well as increased gap junctions (same gap junction weight as the original model in S3A and S3B Fig). The simulation results show that scaled up large networks produced similar experimental decay curves as the two module small network (S3D and S3E Fig). Hence, these results suggest that the gap junction form and basic network configuration are robustly represented in this model.
It is well established that astrocyte dysfunction causes hyper-excitation and the generation or spread of seizure activity; dysfunctional astrocytes should be considered promising targets for new therapies. This study aimed to provide an in-depth understanding of Kir4.1 channels and astrocytic gap junctions in regulating the extracellular K + microenvironment during epileptic seizures in order to facilitate construction of more accurate and dynamic models of neuron-astrocyte networks to improve recognition, forecasting and control of epileptic seizures.

Neuron model
It is assumed that neurons have different local architectures. Each excitatory neuron model has I Na , I K , leak current I L [43], after-hyperpolarization current I AHP . The H-H-type dynamic equations for the two neurons are as follows: where, C is the membrane capacitance. The voltage-gated currents I Na and I K , the leak current I L , and the after-hyperpolarization current I AHP in Eq (1) are: Where g Na and g K denote the conductances corresponding to the sodium and potassium currents, respectively. g AHP is the conductance corresponding to the after-hyperpolarization current. g NaL and g KL are the conductances corresponding to the sodium leak current and potassium leak current, respectively. V Na , V K and V Cl denote the Na + , K + and Clchannel reversal potentials, respectively. n, m and h are gating variables for sodium and potassium currents. g Ca and V Ca are the conductance and the reversal potential of calcium, respectively.
[Ca] i corresponds to the intracellular calcium concentration, and the dynamics equation is: The equations for the gating variables in Eq (2) are: The reversal potentials of Na + and K + and Clare given by the Nernst equation:

Balance of ions fluxes
The [K] o value is continuously updated by the K + currents across the neuronal membrane, Na + /K + -ATP pumps of neuron, K + spatial diffusion [6,23,48], Na + /K + -ATP pumps and Kir4.1 channels of four astrocytes. The electrical current through cell membrane can cause ionic concentration of the inside and outside cell changing. A electrical current I across a membrane is equal to ion flow per unit of time. Thus, the K + concentration dynamics for the neurons and astrocytes and the extra-neuronal space are described as follows: Similar to the K + dynamics, The [Na] o value is continuously updated by the Na + currents across the neuronal membrane, Na + /K + -ATPase pumps of neuron [6,23], Na + /K + -ATP pumps of four astrocytes. However, the main differences are that the pump exchanges two K + for three Na + ions, leading to the coefficient 3 in front of the pump term. In addition to sodium concentrations was added two constant leak terms J NaLA and J NaLN , Thus, the Na + concentration dynamics for neurons and astrocytes and the extraneuronal space are modeled as follows: Also, different terms in the expressions (6-11) are described as follows: where, ρ is the pump strength of Na + /K + -ATP, [Na] Ni and [Na] Ai are the sodium concentration for neuron and astrocyte, respectively, ε is the spatial diffusion coefficient of K + , and K bath is the K + concentration in the largest nearby reservoir.

Model of membrane potential and Kir4.1 channels of astrocytes
The equation of the astrocyte membrane potential V A is: Where C A is the astrocytic capacitance, the leak current I IA =g IA (V A -V IA ). The inward rectifier Kir4.1 channels current I kir in astrocyte depends on the membrane potential and the extracellular potassium ion concentration, the expression can be described as [31,32]: Where g kir and E kir are the conductance and the Nernst constant for the astrocyte Kir 4.1 channels, respectively.

Astrocytic coupling
The K + dynamics in two astrocytes coupled by gap junction is described as follows: The electrical current through the cell membrane can cause the ionic concentrations inside and outside of the cell to change. An electrical current I across a membrane is equal to the ion flow per unit of time. We convert the electrical current I into the ionic flux J, which is computed from J = I/(C Ã γ), where C is the astrocytic membrane capacitance. Here, γ is the scaling factor that relates the ion flux to the membrane potential.
where J gap,k,ij is the potassium ions flow model mediated by astrocytic gap junction: where F is the strength of gap junction between astrocytes. θΔ ij K is the K + concentration gradient for two adjacent astrocytes. The threshold function θ(x) = x if x > 0, and 0 otherwise. The astrocytic gap junction mediated K + buffer allows astrocytic network to respond to changes of the extracellular K + concentration. The values of the parameters used in the model are listed in Table 1.
In simulating the numerical results of gap junction dysfunction that induces extracellular K + buffering delay, we assumed that neuron 1 (N1) and neuron 2 (N2) share the same potassium bath concentration (in other words, neuron 1 (N1) and neuron 2 (N2) have the same value of K bath in Eq (12)). The 4th-order Runge-Kutta method was used for numerical simulation with a time step of h = 0.01 ms. Additionally, parameter values used in the numerical simulation are shown in Table 1, assuming no special emphasis.