The Neuroglial Potassium Cycle during Neurotransmission: Role of Kir4.1 Channels

Neuronal excitability relies on inward sodium and outward potassium fluxes during action potentials. To prevent neuronal hyperexcitability, potassium ions have to be taken up quickly. However, the dynamics of the activity-dependent potassium fluxes and the molecular pathways underlying extracellular potassium homeostasis remain elusive. To decipher the specific and acute contribution of astroglial Kir4.1 channels in controlling potassium homeostasis and the moment to moment neurotransmission, we built a tri-compartment model accounting for potassium dynamics between neurons, astrocytes and the extracellular space. We here demonstrate that astroglial Kir4.1 channels are sufficient to account for the slow membrane depolarization of hippocampal astrocytes and crucially contribute to extracellular potassium clearance during basal and high activity. By quantifying the dynamics of potassium levels in neuron-glia-extracellular space compartments, we show that astrocytes buffer within 6 to 9 seconds more than 80% of the potassium released by neurons in response to basal, repetitive and tetanic stimulations. Astroglial Kir4.1 channels directly lead to recovery of basal extracellular potassium levels and neuronal excitability, especially during repetitive stimulation, thereby preventing the generation of epileptiform activity. Remarkably, we also show that Kir4.1 channels strongly regulate neuronal excitability for slow 3 to 10 Hz rhythmic activity resulting from probabilistic firing activity induced by sub-firing stimulation coupled to Brownian noise. Altogether, these data suggest that astroglial Kir4.1 channels are crucially involved in extracellular potassium homeostasis regulating theta rhythmic activity.


Introduction
Astrocytic processes enwrap more than half of CA1 hippocampal synapses to form tripartite synapses [1,2]. Perisynaptic astroglial processes are enriched in ionic channels, neurotransmitter receptors and transporters, enabling astrocytes to detect neuronal activity via calcium signaling [3] and ionic currents with various components, such as glutamate and GABA transporter [4][5][6][7] or potassium (K + ) [8][9][10]. Thus astrocytes regulate neuronal activity through multiple mechanisms, involving signaling or homeostasis of extracellular space volume, glutamate, GABA or K + levels [11]. Interestingly, membrane depolarization was the first activitydependent signal identified in glial cells and was attributed to K + entry across their membrane [10]. Such K + entry was suggested to contribute to K + spatial buffering, consisting in glial uptake of excess extracellular K + ([K + ] o ), redistribution via gap-junction astroglial networks and subsequent release at sites of low [K + ] o [12].
Modeling studies have mostly investigated astroglial regulation of [K + ] o during pathological conditions to clarify its impact on aberrant neuronal activity. In particular astrocytes, by regulating [K + ] o , have been shown to contribute to initiation and maintenance of epileptic seizures [13][14][15], as well as to the severity of ischemia following stroke, with a neuroprotective or neurotoxic role, depending on [K + ] o [16,17]. In addition, experimental data suggest that several K + channels or transporters contribute to astroglial K + clearance, such as inward rectifier 4.1 and two pore K + channels (K ir 4.1 and K 2P , respectively) or Na/K ATPases [18,19]. Remarkably, recent work suggest that K ir 4.1 channels play a prominent role in astroglial regulation of [K + ] o [20][21][22][23]. However, the mouse model used to draw these conclusions, i.e. conditional K ir 4.1 knockout mice directed to glial cells (GFAP-Cre-K ir 4.1fl/fl mice, K ir 4.1 -/-), exhibits several limitations: 1) K ir 4.1 channels are not specifically deleted in astrocytes, but also in other glial cells such as oligodendrocytes or retinal Müller cells [22]; 2) astrocytes are severely depolarized [21,22]; 3) K ir 4.1 -/mice die prematurely (~3 weeks) and display ataxia, seizures, hindleg paralysis, visual placing deficiency, white matter vacuolization and growth retardation [22], highlighting that chronic deletion of K ir 4.1 channels induces multiple brain alterations and possibly compensations. Thus, the specific and acute contribution of astroglial K ir 4.1 channels to [K + ] o and to the moment to moment neurotransmission is still unclear. To decipher the acute role of astrocytes in controlling K + homeostasis and neuronal activity, we built a tricompartment model accounting for K + dynamics between neurons, astrocytes and the extracellular space. We quantified K + neuroglial interactions during basal and high activity, and found that K ir 4.1 channels play a crucial role in K + clearance and astroglial and neuronal membrane potential dynamics, especially during repetitive stimulations, and prominently regulate neuronal excitability for 3 to 10 Hz rhythmic activity.

Modeling potassium dynamics between neuronal, glial and extracellular compartments
To model K + ions dynamics during neuronal activity, we built a biophysical model that includes three compartments: the neuron, the astrocyte and the extracellular space (Fig. 1A). As performed in several studies [13,16,24], the neuron is approximated by a single compartment conductance-based neuron containing Na + and K + voltage-gated channels, enabling action potential discharge. The associated neuronal membrane potential is coupled with the dynamics of intracellular and extracellular Na + and K + levels via the dependence of the neuronal currents to the Nernst equation. The ion concentrations depend also on the activity of neuronal and astroglial Na/K ATPases, which maintain resting [K + ] i by balancing K + and Na + fluxes. Similarly, the astrocyte is approximated by a single compartment conductance-based astrocyte containing K ir 4.1 channels, which are inward rectifier K + channels strongly expressed in astrocytes that generate dynamic K + currents [25]. In the model, neurons and astrocytes are separated by a homogenous extracellular space compartment. The model is based on balancing ionic fluxes between the three compartments (Fig. 1B). The model starts with the induction of a synaptic current (I app , see Materials and Methods). This current is the initial input of a classical Hodgkin-Huxley model, which describes the neuronal membrane potential dynamics (entry of Na + and exit of K + ). Released extracellular K + is taken up by astrocytes through K ir 4.1 channels and Na/K ATPases ( Fig. 1B and Materials and Methods). Because K ir 4.1 channels are strongly involved in K + uptake [22], we fitted the I-V curve of K + ions through K ir 4.1 channels using equation 22 (see materials and methods) and predicted the I-V curve at various values of [K + ] o (Fig. 1C). We obtain that K + fluxes through K ir 4.1 channels vanish around astrocytic resting membrane potential (~-80 mV) and are outward during astrocytic depolarization for a fixed [K + ] o (2.5 mM, Fig. 1C). However, they become inward when [K + ] o increases (5-10 mM, Fig. 1C). Using this model, we shall investigate quantitatively the contribution of K ir 4.1 channels to K + uptake in relation to neuronal activity associated with different [K + ] o .

Astroglial membrane potential dynamics induced by stimulation
To validate our tri-compartment model, we compared simulation results with electrophysiological recordings. To account for the synaptic properties of CA1 pyramidal neurons, we generated a synaptic current (I app ) using the depression-facilitation model (equation 1) (see Materials and Methods with input f(t) = δ(t)) ( Fig. 2A,E,I). We first investigated responses to single stimulation. Using the Hodgkin-Huxley model, this synaptic current induces a firing activity (S1A Fig.), resulting in a~0.9 mM increase of [K + ] o within 300 milliseconds, which slowly decayed back to baseline levels during 10 seconds (S1B Fig.). The extracellular K + dynamics was associated in our model with a small astrocytic depolarization of ΔV = −1.35 mV (equations 22, 23, 25) (Fig. 2C). Using electrophysiological recordings of evoked field excitatory postsynaptic potential (fEPSP) by a single stimulation of Schaffer collaterals in acute hippocampal slices (Fig. 2B), we measured astroglial membrane potential depolarization and found that it reached~1.3 mV (1.3 ± 0.2 mV, n = 6) (Fig. 2C), confirming the result of our simulation. After validating the responses of the tri-compartment model to basal stimulation, we investigated the impact of trains of stimulations on the dynamics of astroglial membrane potential. During tetanic stimulation (100 Hz for 1 second), variations in neuronal membrane potential described by the Hodgkin-Huxley equation show a bursting activity during~1 second (S1C Fig.). This is associated with a depolarization of astrocytic membrane potential of~5 mV, which lasts~6 seconds (Fig. 2G,H) and an increase in [K + ] o that reaches a peak value of 4.4 mM (S1D Fig.).
For repetitive stimulations (10 Hz for 30 seconds), the neuron exhibited firing activity during the whole stimulation (S1E Fig.). This was associated with an astroglial depolarization of 12 mV (Fig. 2K) and an increase in [K + ] o peaking at 6.9 mM after 17.5 seconds of stimulation (S1F Fig.). Although the stimulation lasted 30 seconds, the astrocytic depolarization started to decay after 17 seconds (Fig. 2K).
The kinetics of astroglial membrane potential dynamics obtained with the numerical simulations are comparable to the results obtained with electrophysiological recordings performed in individual astrocytes during single stimulation (rise time: 48.4 ms for numerical stimulation, 42 ± 19 ms n = 6 for experiments; time of peak: 740 ms for numerical simulation, 730 ± 60 ms n = 6 for experiments; decay time: 3.67 s for numerical simulation, 4.50 s ± 0.2 n = 6 for experiments, Fig. 2D Fig. 2L). These data show that the dynamics of astroglial membrane potential changes obtained from numerical simulations and from electrophysiological recordings are similar. Thus our model captures the key players sufficient to mimic the evoked astroglial membrane potential dynamics observed experimentally in different regimes of activity. Tri-compartment model of the potassium cycle between the neuron, the extracellular space and the astrocyte. A, Schematic representation of the tri-compartment model: neuronal activity induces the release of K + in the extracellular space, which is taken up by astrocytes. B, Reduction of the tricompartment model to ionic fluxes exchanges between a generic postsynaptic neuron, astrocyte and extracellular space. The model includes channels and pumps carrying K + and Na + ions. C, Current-Voltage relationship (I-V curve) of K ir 4.1 channels. We identify the free parameters in equation 22 by fitting the simulated IV curve (light blue) to experimental recordings performed in isolated astrocytes (sampling data, black rectangles) [65]. Using equation 22, we plot the I-V curve for different ratios of extracellular to intracellular astrocytic K + concentrations (2.5/135 (light blue), 5/135 (blue) and 10/135 (dark blue). At resting membrane potential (-80 mV) and resting [K + ] o (2.5 mM), the K ir 4.1 current is outward, but as illustrated here, it reverses by increasing [K + ] o . Potassium redistribution in neuronal, astroglial and extracellular space compartments for different regimes of activity We investigated the dynamics of the K + cycle between neurons, extracellular space and astrocytes induced by neuronal activity to decipher the time needed to restore basal extracellular and intra-neuronal K + levels. We studied K + redistribution induced by single, tetanic (100 Hz, 1 s) and repetitive (10 Hz, 30 s) stimulations, and found that the general behavior of K + dynamics was divided into three phases (phases 0, 1 and 2; Fig. 3).
During phase 0 (t = 0 to t1), neuronal K + is released in the extracellular space (peak [K + ] o during phase 0: 0.9 mM at 300 ms for single stimulation; 1.9 mM at 1.3 s for tetanic stimulation; 4.4 mM at 30 s for repetitive stimulation, Fig. 3A,D,G). Compared to basal K + levels in each compartment (at t = 0), the relative transient evoked increase in K + concentration is prominent only in the extracellular space (~+37% for single stimulation, +76% for tetanic stimulation and +168% for repetitive stimulation, Fig. 3B,E,H).
During phases 0 and 1, released K + is then mostly buffered by astrocytes (~80 to 99% at the end of phase 1) during the different regimes of activity (time t2 (at the end of phase 1) for single stimulation: 8.  Finally during phase 2 (t2 to end), astroglial buffered K + is slowly redistributed back to neurons, which ends the K + cycle. The long-lasting phase 2 is marked by an inversion of K + fluxes in astrocytes, suggesting moderate K + release by astrocytes over time. Indeed, K + redistribution to neurons depends on K + release through K ir 4.1 channels, which is limited by the low outward rectification of these channels (Fig. 1C). Altogether, these data suggest a slow, but dynamic and efficient astroglial clearance capacity for the different regimes of activity. For all regimes of activity, neuronal K + (red) is released, increasing K + in the extracellular space (black) during the stimulation initiated at time t = 0 (phase 0, t = 0 to t1), and is then cleared by the astrocyte (blue) (phase 1, t1 to t2). K + levels are illustrated for the different regimes in each compartment (A, D,G) and are normalized to basal [K + ] o (B,E,H) or to the total amount of released K + by neurons (C,F,I). Finally, the buffered K + is slowly redistributed back to neurons, which ends the K + cycle (phase 2, t2 to end). t1 represents the time point where neuronal release of K + stops, whereas t2 is the time point where astroglial K + uptake peaks. doi:10.1371/journal.pcbi.1004137.g003

K ir 4.1 channel contribution to neuronal firing and extracellular K + levels
To study quantitatively the acute and selective role of astroglial K ir 4.1 channels in neuroglial K + dynamics, we inhibited the K ir 4.1 current in our tri-compartment model. Because K ir 4.1 -/mice display altered synaptic plasticity compared to wild type mice [22,26], we recalibrated the synaptic current (I app ) parameters τ rec and τ inact in equations 1,2 (see Table 1) for the facilitation-depression model to get an optimal fit to the recorded postsynaptic responses [26]. Another change in the model consisted in setting at zero both the K ir 4.1 current and the leak term. In addition, to compensate for the loss of K + fluxes through astroglial K ir 4.1 channels, we added in equation 27 a constant K + flux to maintain [K + ] o at an equilibrium value of 2.5 mM. Consequently, the astrocytic membrane potential displayed no change during stimulation, in agreement with electrophysiological recordings [21,22].
The numerical simulations show that inhibition of astroglial K ir 4.1 channels leads to higher transient peak increase in [K + ] o during repetitive and tetanic stimulation compared to control conditions ( Fig. 4E,F,I,J), while no difference is observed for single stimulation (Fig. 4A,B). In  Fig. 4B,F,J). Finally, K ir 4.1 channel inhibition only slightly increased neuronal firing induced by single stimulation (Fig. 4C,D) and tetanic stimulation (Fig. 4G,H), while it had major effect on neuronal excitability during repetitive stimulation (Fig. 4K,L). Indeed, although firing frequency was only slightly increased during the first 8 seconds of repetitive stimulation when [K + ] o reached 10 mM (Fig. 4I), action potential amplitude and firing rate then progressively decreased due to neuronal depolarization (from-33 mV to-19 mV after 14 and 30 seconds of stimulation, respectively), suppressing neuronal firing after 14 seconds of stimulation (Fig. 4K). Altogether, these data show that astroglial K ir 4.1 channels are prominently involved in K + buffering during high level of activity, and thereby have a major impact on neuronal resting membrane potential controlling firing during trains of stimulations.
Astrocytic K ir 4.1 channels modulate firing probability induced by low frequency sub-firing stimulation in noisy neurons To investigate the effect of astroglial K ir 4.1 channels on endogenous physiological rhythmic activity, we generated probabilistic firing induced by sub-firing stimulation coupled to neuronal Brownian noise (Fig. 5A,B). To simulate the firing activity, we generated a sub-firing periodic stimulation (5 ms squared stimulus), which defines the applied synaptic intensity in our tripartite compartment model (Fig. 5A), and added a Brownian noise in the neuronal membrane potential (equation 21, Fig. 5B). Such stimulation induces an increase in [K + ] o (Fig. 5C), and thus firing over time (Fig. 5D-E). We found that astroglial K ir 4.1 channels had no effect on the firing probability (computed over 100 simulations) for basal (0.1 Hz, Fig. 5F), low (1 Hz, Fig. 5G) and high (50 Hz, Fig. 5K) frequency stimulations. However, K ir 4.1 channels directly regulate the firing probability for 3 and 5 Hz stimulations after 7 and 12 s of sub-firing stimulation, respectively (Fig. 5H,I). In contrast, K ir 4.1 channels regulate only transiently the firing probability induced by 10 Hz stimulation (Fig. 5J). These data suggest a prominent and specific involvement of astroglial K ir 4.1 channels in regulation of firing during theta rhythmic activity.

Discussion
[K + ] o modulates neuronal membrane potential, excitability, release probability and synaptic efficacy [27][28][29][30][31][32]. To unravel the acute role of astrocytes in extracellular K + homeostasis and neuronal activity, we used electrophysiological recordings with a tri-compartment model accounting for K + dynamics between neurons, astrocytes and the extracellular space. We found that K ir 4.1 channels play a key role in extracellular K + clearance, astroglial and neuronal membrane potential dynamics, especially during trains of stimulation, and strongly regulate neuronal excitability for slow rhythmic activity (3-10 Hz).
To study seizure discharges and spreading depression, a first tri-compartment model including the neurons, astrocytes and extracellular space was proposed [24], although the astrocytic membrane potential was not taken into account, and K + accumulation in the interstitial volume was controlled by a first-order buffering scheme that simulated an effective glial K + uptake system. With such model, after evoked firing, it took~17 s for the neuronal membrane potential to return to resting values, via activation of Na/K ATPases. The model also  predicted that elevated [K + ] o have a key role in the initiation and maintenance of epileptiform activity. In our study, we accounted for the astroglial modulation of K + buffering capacity regulated by its membrane potential, and found that the biophysical properties of astrocytic membranes including K ir 4.1 channels were sufficient to account for the long-lasting clearance of extracellular K + . Interestingly, we confirm that alteration in K + clearance leading to an extracellular K + accumulation induces epileptiform activity, and show specifically that K ir 4.1 channel acute inhibition leads to such pathological bursting activity during repetitive stimulation.
A similar tri-compartment model has been simplified as a one-dimensional two-layer network model to study how neuronal networks can switch to a persistent state of activity, as well as the stability of the persistent state to perturbations [13]. In this model, Na + and K + affect neuronal excitability, seizure frequency, and stability of activity persistent states. In particular, the quantitative contribution of intrinsic neuronal currents, Na/K ATPases, glia, and extracellular Na + and K + diffusion to slow and large-amplitude oscillations in extracellular and neuronal Na + and K + levels was revealed. In the model, the estimated [K + ] o during epileptiform activity are comparable to the ones observed experimentally [35,36]. Although this model does not account for astroglial K ir 4.1 channels, it shows that a local persistent network activity not only needs balanced excitation and inhibition, but also glial regulation of [K + ] o [15]. Finally, a model accounting for the extracellular space and astroglial compartments has quantified the involvement of several astroglial ionic channels and transporters (Na/K ATPase, NKCC1, NBC, Na + , K + , and aquaporin channels) in the regulation of firing activity [34].
To account for K + dynamics between neurons, astrocytes and the extracellular space, we built for the first time a tri-compartment model, where we included neuronal voltage-gated channels, Na/K pumps and astrocytic K ir 4.1 channels according to their biophysical properties, as well as membrane potential of astrocytes. Because functional expression of voltage-gated calcium channels on hippocampal mature astrocytes in situ in physiological conditions and its impact on astrocytic functions is still a matter of debate [37], such channels were not included in our model. However, many other astroglial K + channels (such as two pore domain K + channels (K 2P ) (TWIK-1, TREK-1, TREK-2 and TASK-1), inward rectifier K + channels (Kir2.1, 2.2, 2.3, 3.1, 6.1, 6.2), delayed rectifier K + channels (Kv1.1, 1.2, 1.5, 1.6), rapidly inactivating A-type K + channels (Kv1.4), calcium-dependent K + channels (K Ca 3.1)), but also other channels, transporters or exchangers (such as Cx hemichannels, Na + /K + /Clco-transporter (NKCC1) K + /Clexchanger, glutamate transporters) [16,38,39] could also play a role in the regulation of activity-dependent changes in [K + ] i or [K + ] o . Functional evidence of the contribution of these channels, transporters or exchangers in astroglial K + clearance is actually scarce, although K 2P channels have been suggested to participate in astroglial K + buffering [40], while NKCC1 were recently shown in hippocampal slices not to be involved in activity-dependent K + clearance [41]. Similarly, adding slower timescale K + dependent conductances in the neuron model could modulate the slow redistribution of K + to neurons, and thus the duration of the neuroglial potassium cycle, and is of interest to implement in future development of the model. In our study, the aim was to simplify the system to capture in the model the minimal set of astroglial channels and pumps accounting for our experimental data related to activity-dependent changes in astroglial membrane potential. In addition our tri-compartment model, as most existing models, did not account for the complex multiscale geometry of astrocytes and neurons. Incorporating in our current model additional astroglial and neuronal channels, as well as complex cell geometry is of particular interest to identify modulatory effects of other specific channels and of microdomain geometry on the neuroglial potassium cycle.
In accordance with previous studies, where K ir 4.1 channels were chronically deleted genetically in glial cells [20,21,23], we found that acute inhibition of K ir 4.1 channels leads to altered regulation of extracellular K + excess and affects the kinetics of [K + ] o (Fig. 4I,J). However, in contrast to these studies, we found that K ir 4.1 channel inhibition also alters significantly [K + ] o peak amplitudes during repetitive stimulation, suggesting that K ir 4.1 -/mice may display some compensatory mechanisms attempting to maintain extracellular K + homeostasis. In addition, our model reveals that specific and acute inhibition of K ir 4.1 channels slows down, but does not abolish, astroglial uptake of excess K + during single, tetanic and repetitive stimulations, confirming that astroglial Na/K ATPases, included in our model, also contribute to K + clearance [41].
The long-lasting astrocytic potassium uptake is due in part to the slow K ir 4.1 conductance dynamics Contrary to action potentials, characterized by a very fast dynamics in the order of a few milliseconds, astroglial K + buffering lasts tens of seconds. As shown in the present study, most of extracellular K + released by neurons is first cleared by astrocytes through K ir 4.1 channels. To determine the factors controlling the slow timescale of astroglial K + clearance, we focused on K ir 4.1 channels. Because the astroglial leak conductance (equation 23) is six times smaller than the K ir 4.1 channel conductance, we neglected it.
The dynamics of astrocytic membrane potential V A is described by equation 23, where the membrane capacitance is C A % 15 pF and the maximal K ir 4.1 channel conductivity isG K ir % 60pS. In that case, using equation 23, the time constant of K ir 4.1 channel-mediated return to equilibrium of astroglial membrane potential τ A is defined as We obtain the following approximation τ A % 0.6s using equation 23 and the parameters of table 1. This time constant is consistent with the fitted exponential decay time obtained in our simulations and experiments for a single stimulation where we obtained τ % 0.7s. However, simulations for stronger stimulations indicate an increase of τ to approximatively 4 seconds (tetanic stimulation) and 9 seconds (repetitive stimulation). This increase in clearance duration is due to the dependence of the K ir 4.1 current to [K + ] o , as illustrated by the IV relation (Fig. 1C) and described in equation 22. The Nernst potential V KA increases for strong stimulations (tetanic and repetitive), which slow down the kinetics of astrocytic membrane potential in equation 22. We conclude that the slow time scale of K + clearance is in part due to the availability of K ir 4.1 channels at low and high [K + ] o . This clearance timescale is much longer than the glutamate clearance rate of τ glu % 15 ms that we previously reported [42]. Moreover, the redistribution of K + released by neurons during the different regimes of activity shows that the higher the activity, the lower the proportion of released K + remains transiently in the extracellular space. This suggests that K ir 4.1 channels have a strong uptake capacity, especially for high regimes of activity ([K + ] o up to 5-6 mM).

Impact of K ir 4.1-mediated potassium buffering on neuronal activity in physiology and pathology
Remarkably, our model reveals that astroglial K ir 4.1 channels strongly regulate neuronal firing induced by high stimulation regime such as repetitive stimulation. K ir 4.1channels are crucially involved in regulation of [K + ] o during this regime of activity, most likely because such stimulation triggered long-lasting neuronal release of K + (20 mM over 30 seconds, Fig. 3G) resulting in a sustained, but moderate increase in [K + ] o (>6 mM for~22 s, Fig. 4I,J), compared to the neuronal release. These data suggest that during repetitive stimulation, astrocytes can buffer up to~14 mM of [K + ] o and thereby preserve neuronal firing. However, astroglial K ir 4.1 channels slightly impact neuronal firing induced by single and tetanic stimulations, probably because these stimulations only triggered transient neuronal K + release (0.9 mM over 300 ms (Fig. 3A) and 1.9 mM over 1.3 s (Fig. 3D), respectively), resulting in a short and small increase in [K + ] o (>2.7 mM for~450 ms for single stimulation (Fig. 4A,B), and >3.5 mM for 1.5 s for tetanic stimulation (Fig. 4E,F)). Nevertheless, we show a prominent and specific involvement of astroglial K ir 4.1 channels in probabilistic firing activity induced by 3 to 10 Hz sub-firing stimulations (Fig. 5), suggesting a key role of these channels in sustained theta rhythmic activity. Interestingly, these data imply that K ir 4.1 channels can contribute to fine tuning of neuronal spiking involving low, but long-lasting, increase in [K + ] o . Thus besides gliotransmission, regulation of [K + ] o by K ir 4.1 channel provides astrocytes with an alternative active and efficient mechanism to regulate neuronal activity. Several studies have reported decreased K ir 4.1 protein levels and K ir functional currents in sclerotic hippocampus from epileptic patients [43][44][45][46]. Whether these changes are the cause or the consequence of epilepsy is still an open question. However, K ir 4.1 -/mice display an epileptic phenotype [22,47] and missense mutations in KCNJ10, the gene encoding K ir 4.1, have been associated with epilepsy in humans [48,49]. These data thus suggest that impairment in K ir 4.1 function leading to alterations in [K + ] o dynamics, as shown in our study, may cause epilepsy. Remarkably, dysfunction of [K + ] o regulation by K ir 4.1 channels is likely involved in other pathologies, since it contributed to neuronal dysfunction in a mouse model of Huntington's disease [50] and the presence of antibodies against K ir 4.1 channels in glial cells was recently found in almost 50% of multiple sclerosis patients [51]. Thus astroglial K ir 4.1 channels may well represent an alternative therapeutic target for several diseases.

Materials and Methods
Ethics statement
Modeling potassium dynamics in the tripartite neuron-astrocyteextracellular compartment We present here the biophysical model we have built to describe K + dynamics during neuronal activity and specifically the role of astroglial K ir 4.1 channels. After Schaffer collateral stimulation, excitatory synapses release glutamate molecules that activate postsynaptic neurons. We modeled this step by classical facilitation/depression model [55]. The resulting postsynaptic activity triggers ionic release in the extracellular space and a change in the astrocytic membrane potential through ion uptake. We used the average neuronal potential and mass conservation equations for ionic concentrations to model changes in astrocytes. We have built a tri-compartment model, which accounts for: 1) the neuron, 2) the astrocyte and 3) the extracellular space. We included voltage gated channels, Na/K pumps and astrocytic K ir 4.1 channels.

Facilitation-depression model
To account for the stimulation of Schaffer collaterals that induce a postsynaptic response in the CA1 stratum radiatum region, we used a facilitation-depression model [55][56][57].
where f is the input function. For a single stimulation generated at time t stim , f(t) = δ(t-t stim ). A stimulation instantaneously activates a fraction U se of synaptic resources r, which then inactivates with a time constant τ inac and recovers with a time constant τ rec In the simulations, at time t = t stim , r and e respectively decreases and increases by the value U se r. The synaptic current I app is proportional to the fraction of synaptic resources in the effective state e and is given by I app = A se e (the parameter A se is defined in table 1). We used the following definitions for the input function f:

Modeling neuronal activity
The dynamics of the neuronal membrane potential, V N , follows the classic Hodgkin Huxley (HH) equations [58].
with rate equations V rest is the resting membrane potential and V KN and V NaN are respectively the K + and Na + equilibrium potentials and are given by the Nernst equations where Na 0 and Na N are respectively the extracellular and neuronal sodium concentrations, and K 0 and K N are respectively the extracellular and neuronal K + concentrations that may vary as we shall describe below. We complete the description of all the neuronal currents with a leak current which stabilizes the membrane potential at its resting value. Finally, the neuronal membrane potential satisfies the equation where I app is the synaptic current derived from equation 1.
Modeling astrocytic potassium uptake by K ir 4.1 channels To account for the K + dynamics in astrocytes, we modeled the K ir 4.1 channel according to its biophysical properties [59] and I-V curve [60]. The total astroglial current I Kir depends on the membrane potential, the extracellular (K 0 ) and the astrocytic (K A ) K + concentrations, and is approximated by where V KA is the Nernst astrocyte K + potential, V A , the astrocyte membrane potential, K 0 is the extracellular K + concentration and V A1 (an equilibrium parameter, which sets K ir current to 0 at-80 mV), V A2 and V A3 are constant parameters calibrated by the I-V curve (Fig. 1C, [60]), as detailed below. The second term of equation 22 describes the dependence of I Kir to the square root of K 0 [60][61][62][63][64] and to the steady state open/close partition function of Kir4.1 channels according to the Boltzmann distribution [59], which includes dynamic variations of potassium Nernst potential during neuronal activity. Adding a leak current I lA = g lA (V A -V lA ), which stabilizes the astrocyte membrane potential at-80 mV, the astrocyte membrane potential V A satisfies the equation where I Kir is defined by relation 22. We fitted the K ir 4.1 channel I-V curve (equation 22) using the experimental recordings for the K ir 4.1 channel (3 mM [K + ] (Fig. 4 in [60,65]). We first obtained that V A1 = (V restA − 26ln(3/ 145)) = −14.83 mV where V restA = −80mV (potential for which the current is zero). We then used the Matlab fitting procedure for a single exponential with formula 22 changed to ðVÀV A1 À26lnð3=145ÞÞ ffi ffi 3 p I with (V from-100 to 20 mV) to get that V A2 = 34 mV and V A3 = 19.23 mV (table 1). Varying [K + ] o by 0.5 mM did not affect significantly the K ir 4.1 channel I-V curve, confirming its robustness.

Na/K pump ionic flux for astrocytes and neurons
The K + resting concentrations in neurons and astrocytes are maintained by Na/K pumps that balance the outward K + and inward Na + fluxes. The associated pump currents i pump,k (index k = N for the neuron, k = A for the astrocyte) depend on the extracellular K + K 0 and intracellular Na + concentrations (Na N for the neuron and Na A for the astrocyte) and follow the same equation as [66], where i maxk is a constant (table 1).

Balance of ionic fluxes
We converted the different electrogenic neuronal and astrocytic channel currents into ionic fluxes [13]. A current I across a membrane induces a flow of charge i equals to δQ = I per unit of time. The corresponding change in extracellular concentration is given by I/(qN A Vol 0 ), where q = 1.6 Ã 10 -19 C is the charge of an electron, N A the Avogadro number and Vol N , Vol A andVol 0 are the neuronal, astrocytic and extracellular volume respectively. To model the ionic concentration dynamics, we converted the currents I Na , I K and I Kir to the corresponding ionic fluxes i Na , i K and i Kir We describe in the following paragraphs the equations for the ionic concentrations in the three compartments (neuron, extracellular space and astrocyte).

Potassium fluxes
To determine the system of equations for the K + fluxes, we use the mass conservation law for the extracellular K 0 , the neuronal K N and the astrocytic K A K + concentrations. The extracellular K + K 0 increases with the neuronal current I K (see equation 8), which is here converted to i K (ion flux), but it is also uptaken back into neurons with a flux 2 i pumpN (the factor 2 is described in [67] and into astrocytes as the sum of the two fluxes 2 i pumpA plus i Kir . Similarly, we obtain the equations for the neuronal and astrocytic K + to balance the various fluxes. Finally, we get To study quantitatively the acute and selective role of astroglial K ir 4.1 channels in neuroglial K + dynamics, we inhibited the K ir 4.1 current in our tri-compartment model. We thus set at zero both the K ir 4.1 current and the leak term. To compensate for the loss of K + fluxes through astroglial K ir 4.1 channels, we added in equation 27 a constant K + flux to maintain [K + ] o at an equilibrium value of 2.5 mM. This constant K + flux in astrocytes could be mediated by various channels or transporters such as two pore domain potassium channels (K2P such as TWIK-1, TREK-1, TREK-2 and TASK-1), delayed rectifier potassium channels (Kv1.1, 1.2, 1.5 and 1.6), rapidly inactivating A-type potassium channels (Kv1.4), glutamate transporters or connexin43 hemichannels. However, since TASK-1 [68] and Cx43 hemichannels [69] are thought to be active in basal conditions, they are more likely to mediate such flux.

Sodium fluxes
Similarly to the K + dynamics, the equations for the Na + fluxes are derived using the balance between the neuronal, astrocytic and extracellular concentrations. However, the main differences are that the pump exchanges 2 K + for 3 Na + ions, leading to the coefficient 3 in front of the pump term. In addition, to stabilize the sodium concentrations, we added two constant leak terms i NalA and i NalN (values given in table 1), as classically used [24], Numerical implementations and fitting procedures Numerical simulations. Simulations, numerical integrations and fitting computations were performed in Matlab. We used Runge Kunta fourth order method for the simulations, which were numerically stable. We used a time step of Δt = 0.1 ms (simulations were repeated with smaller time step to check whether numerical accuracy was affecting results). The leak currents parameters were adjusted to stabilize the model at the resting membrane potentials (-70 mV and -80 mV for neurons and astrocytes respectively) and resting concentrations (neuronal [K + ] and [Na + ]: 135 mM and 12 mM, respectively; extracellular [K + ] and [Na + ]: 2.5 mM and 116 mM, respectively; astrocytic [K + ] and [Na + ]: 135 mM and 12 mM, respectively). The parameters for the Hodgkin Huxley equations were also adjusted to these concentrations.
Approximation of time constants. Time constants τ of simulated extracellular K + transients were fitted to curves using a single exponential ðe À t t Þ (Fig. 4B,F,J). For all the fits obtained on the numerical simulation curves, we obtained an error estimation R-square ! 0.97. Time constants τ of experimental and simulated astroglial membrane potentials were calculated by computing the rise and decay times between 20% and 80% of the maximal peak amplitude responses (Fig. 2D,H,L).
All time constants τ were fitted to curves using a single exponentialðe À t t Þ. For all the fits obtained on the numerical simulation curves, we obtained an error estimation R-square ! 0.97.
Approximation of facilitation/depression model parameters. To account for the synaptic properties of CA1 pyramidal neurons following single, tetanic and repetitive stimulations, we generated a synaptic current using the depression-facilitation model (equation 1) where I app depends on the input functions f s (t) (equation 4), f TT (t) (equation 5) and f Rs (t) (equation 6), respectively ( Fig. 2A,E,I). The synaptic current parameters were fitted to experimental recordings [26] by matching the time of maximal peak amplitude of fEPSP with the one of I app in control conditions (τ = 300 ms, τ inact = 200 ms). The parameters for the K ir 4.1 inhibition condition in the model were extracted from our experimental results on K ir 4.1 glial conditional knockout mice [26] and are given by τ rec = 500 ms, τ inact = 160 ms. When K ir 4.1 channels are inhibited (model) or knockedout (experiment), the maximal peak amplitudes of the applied synaptic currents in the model (I app ) and fEPSPs recorded experimentally are increased compared to control conditions [26].

Simulation of neuronal firing at various frequencies
We imposed an initial input at various frequencies (0.1, 1, 3, 5, 10, 50 Hz). Each input is generated by a sub-firing square current lasting 5 ms (I app ). In addition, we added a Brownian noise of amplitude σ = 0.68 pA 2 ms -1 to induce neuronal membrane potential fluctuation (equation 21), which amplitude (1 mV) was chosen to induce a probabilistic firing of 0.2, matching the CA1 pyramidal cells synaptic release probability p = 0.2 (probability to induce a postsynaptic response in equation 1) [70]. Using the tri-compartment model, we simulated at various frequencies a quantity that we called the observed firing probability defined empirically at time t as the time dependent ratio of the number of spikes observed at time t to the total number of simulations.