Abstract

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 KiA.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 Ki,4.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 Ki,4.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 Ki,4.1 channels are crucially involved in extracellular potassium homeostasis regulating theta rhythmic activity.

Author Summary

In particular, potassium ions, released by neurons during activity, have to be taken up efficiently to prevent hyperexcitability. Astrocytes, the third element of the synapse, play a prominent role in extracellular potassium homeostasis. Thus unraveling the dynamics of the neuroglial potassium cycle during neurotransmission and the underlying molecular pathways is a key issue. Here, we have developed a tri-compartment model accounting for potassium dynamics between neurons, astrocytes and the extracellular space to quantify the specific and acute contribution of astroglial Kir4.1 channels to extracellular potassium levels and to the mo-ment-to-moment neurotransmission. We demonstrate that astroglial Kir4.1 channels are sufficient to account for the slow membrane depolarization of astrocytes and crucially contribute to extracellular potassium clearance during basal and high activity. We also show that astrocytes buffer in less than 10 seconds more than 80% of the potassium released by neurons, leading to recovery of basal extracellular potassium levels and neuronal excitability. Remarkably, we found that Kir4.1 channels also prominently regulate slow 3 to 10 HZ rhythmic firing activity. Altogether, these data show that Kir4.1 channels acutely regulate extracellular potassium and neuronal excitability during specific patterns of activity.

Introduction

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—7] or potassium (K+) [8—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 activity-dependent 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 [KJ']O [12].

In particular astrocytes, by regulating [K+]o, have been shown to contribute to initiation and maintenance of epileptic seizures [13—15], as well as to the severity of ischemia following stroke, with a neuroprotective or neurotoxic role, depending on [KJF]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 (Kir4.1 and K219, respectively) or Na/K ATPases [18,19]. Remarkably, recent work suggest that Kir4.1 channels play a prominent role in astroglial regulation of [K+] 0 [20—23]. However, the mouse model used to draw these conclusions, i.e. conditional Kir4.1 knockout mice directed to glial cells (GFAP-Cre-KhAJfl/fl mice, Kir4.1'/'), exhibits several limitations: 1) Kir4.1 channels are not specifically deleted in astrocytes, but also in other glial cells such as oligodendrocytes or retinal Muller cells [22]; 2) astrocytes are severely depolarized [21,22]; 3) Kir4.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 Kir4.1 channels induces multiple brain alterations and possibly compensations. Thus, the specific and acute contribution of astroglial Kir4.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 tri-compartment 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 Kir4.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.

Results

Modeling potassium dynamics between neuronal, glial and extracellular compartments

1A). As performed in several studies [13,16,24] , the neuron is approximated by a single compartment conductance-based neuron containing Na+ and KL 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 astro-glial 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 Kir4.1 channels, which are inward rectif1er 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. 13). The model starts with the induction of a synaptic current (1%,, see Materials and Methods). This current is the initial input of a classical Hodg-kin-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 Kir4.1 channels and Na/K ATPases (Fig. IB and Materials and Methods). Because Kir4.1 channels are strongly involved in K+ uptake [22], we fitted the IV curve of K+ ions through Kir4.1 channels using equation 22 (see materials and methods) and predicted the IV curve at various values of [K‘L]0 (Fig. 1C). We obtain that K+ fluxes through Kir4.1 channels vanish around astrocytic resting membrane potential (~-80 mV) and are outward during astrocytic depolarization for a fixed [KJF]O (2.5 mM, Fig. 1C). However, they become inward when [KJF]0 increases (5—10 mM,

1C). Using this model, we shall investigate quantitatively the contribution of Kir4.1 channels to K+ uptake in relation to neuronal activity associated with different [K+]O.

Astroglial membrane potential dynamics induced by stimulation

To account for the synaptic properties of CA1 pyramidal neurons, we generated a synaptic current (lap?) using the depression-facilitation model (equation 1) (see Materials and Methods with input f(t) = 5(t)) (Fig. 2A,E,I). We first investigated responses to single stimulation. Using the Hodgkin-Huxley model, this synaptic current induces a firing activity (SlA Fig), resulting in a ~ 0.9 mM increase of [K‘L]O within 300 milliseconds, which slowly decayed back to baseline levels during 10 seconds (SIB Fig). The extracellular K+ dynamics was associated in our model with a small astrocytic depolarization of AV 2 —1.35 mV (equations 22, 23, 25)

2C). Using electrophysiological recordings of evoked field excitatory postsynaptic potential (fEPSP) by a single stimulation of Schaffer collaterals in acute hippocampal slices (Fig. 23), we measured astroglial membrane potential depolarization and found that it reached ~ 1.3 mV

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-Hquey equation show a bursting activity during ~ 1 second (SIC 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‘L]O that reaches a peak value of 4.4 mM (SID Fig).

This was associated With an astroglial depolarization of ~ 12 mV (Fig. 2K) and an increase in [KJF]o peaking at 6.9 mM after 17.5 seconds of stimulation (81F Fig). Although the stimulation lasted 30 seconds, the astrocytic depolarization started to decay after 17 seconds (Fig. 2K).

2H) and repetitive stimulation (rise time: 1.5 s for numerical simulation, 1.27 s i 0.18 n = 5 for experiments; time of peak: 6.8 s for numerical simulation, 5.2 s i 0.9 n = 5 for experiments; decay time: 7.95 s for numerical simulation, 8.3 s i 0.4 n = 5 for experiments, 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.

Potassium redistribution in neuronal, astroglial and extracellular space compartments for different regimes of activity

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).

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 Kir4.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.

Kir4.1 channel contribution to neuronal firing and extracellular K+ levels

Because Kir4.1'/' mice display altered synaptic plasticity compared to wild type mice [22,26] , we recalibrated the synaptic current (lapp) parameters TreC and Tinact in equations 1,2 (see Table 1) for the facilita-tion-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 Kir4.1 current and the leak term. In addition, to compensate for the loss of K+ fluxes through astroglial Kir4.1 channels, we added in equation 27 a constant K+ flux to maintain [KJF]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].

4E,F,I,I), while no difference is observed for single stimulation (Fig. 4A,B). In Parameters Recovery time constant Inactivation time constant Absolute synaptic strength Utilization of synaptic efficacy Neuronal sodium channel conductance Neuronal potassium channel conductance Neuronal resting membrane potential Avogadro Number Net charge of single monovalent ion Faraday Constant Gaz constant Temperature Neuronal leak conductance Neuronal leak potential Neuronal membrane capacitance Astrocytic single channel Kir conductance Kir current potential constant 1 Kir current potential constant 2 Kir current potential constant 3 Astrocytic capacitance Astrocytic leaking potential Astrocytic leak conductance Astrocytic Na/K pump rate Neuronal Na/K pump rate Extracellular space volume/ neuronal volume Extracellular space volume/astrocytic volume Neuronal sodium leak rate Astrocytic sodium leak rate Value 0.07 nS 60 pS 0.1 nS addition, for all regimes of activity, the rise and decay times of the [K+] 0 were increased when Kir4.1 channels were inhibited (single stimulation, control: rise time 136 ms, decay time 3.4 s; Kir4.1 inhibition: rise time 232 ms, decay time 4.2 s; tetanic stimulation, control: rise time 638 ms, decay time 4 s; Kir4.1 inhibition: rise time 753 ms, decay time 6 s; repetitive stimulation, control: rise time 6.8 s; Kir4.1 inhibition: rise time 20.2 s, Fig. 4B,F,]).

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‘L]o reached 10 mM (Fig. 41), 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 Kir4.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 f1ring during trains of stimulations.

Astrocytic Kir4.1 channels modulate firing probability induced by low frequency sub-firing stimulation in noisy neurons

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. SB). Such stimulation induces an increase in [KJ']o (Fig. 5C), and thus firing over time (Fig. 5D-E). We found that astroglial Kir4.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, Kir4.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, Kir4.1 channels regulate only transiently the firing probability induced by 10 Hz stimulation (Fig. 5]). These data suggest a prominent and specific involvement of astroglial Kir4.1 channels in regulation of firing during theta rhythmic activity.

Discussion

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 Kir4.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).

A novel tri-compartment model accounting for astroglial Kir4.1 channels and membrane potential dynamics in K“ regulation of neuronal activity

Several models have investigated extracellular K+ regulation of neuronal activity, including glial uptake mechanisms [13—17,24,33,34]. 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+]0 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 Kir4.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 Kir4.1 channel acute inhibition leads to such pathological bursting activity during repetitive stimulation.

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 [KJ']o during epileptiform activity are comparable to the ones observed experimentally [35,36]. Although this model does not account for astroglial Kir4.1 channels, it shows that a local persistent network activity not only needs balanced excitation and inhibition, but also glial regulation of [K‘L]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].

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 (sz) (TWIK-I, 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 (KCa3.1)), but also other channels, transporters or exchangers (such as Cx hemichannels, Na+/K+/Cl' co-transporter (NKCC1) K+/Cl' exchanger, glutamate transporters) [16,38,39] could also play a role in the regulation of activi-ty-dependent changes in [K] or [K+]O. Functional evidence of the contribution of these channels, transporters or exchangers in astroglial K+ clearance is actually scarce, although sz 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 neurogli-al 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.

41,1). However, in contrast to these studies, we found that Kir4.1 channel inhibition also alters significantly [K‘L]O peak amplitudes during repetitive stimulation, suggesting that Kir4.1'/' mice may display some compensatory mechanisms attempting to maintain extracellular K+ homeostasis. In addition, our model reveals that specific and acute inhibition of Kir4.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 Kir4.1 conductance dynamics

As shown in the present study, most of extracellular K+ released by neurons is first cleared by astrocytes through Kir4.1 channels. To determine the factors controlling the slow timescale of astroglial K+ clearance, we focused on Kir4.1 channels. Because the astroglial leak conductance (equation 23) is siX times smaller than the Kir4.1 channel conductance, we neglected it. The dynamics of astrocytic membrane potential VA is described by equation 23, where the membrane capacitance is C A z 15 pF and the maximal Kir4.1 channel conductivity isGKir % 60pS. In that case, using equation 23, the time constant of Kir4.1 channel-mediated return to equilibrium of astroglial membrane potential TA is defined as

This time constant is consistent with the fitted exponential decay time obtained in our simulations and experiments for a single stimulation where we obtained 1‘ z 0.75. However, simulations for stronger stimulations indicate an increase of T to approximatively 4 seconds (tetanic stimulation) and 9 seconds (repetitive stimulation). This increase in clearance duration is due to the dependence of the Kir4.1 current to [K+]o, as illustrated by the IV relation (Fig. 1C) and described in equation 22. The Nernst potential VKA increases for strong stimulations (tetanic and repetitive), which slow down the kinetics of astrocytic membrane potential VA through the term 1 —|— exp in equation 22. We conclude that the slow time scale of K+ clearance is in part due to the availability of Kir4.1 channels at low and high [KL] 0. This clearance timescale is much longer than the glutamate clearance rate of Tglu % 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 Kir4.1 channels have a strong uptake capacity, especially for high regimes of activity ([KJF]0 up to 5—6 mM).

Impact of Kir4.1-mediated potassium buffering on neuronal activity in physiology and pathology

Kir4.1channels are crucially involved in regulation of [K‘L]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 [KJ']o (>6 mM for ~22 s, Fig. 4I,]), compared to the neuronal release. These data suggest that during repetitive stimulation, astrocytes can buffer up to ~14 mM of [K‘L]O and thereby preserve neuronal firing. However, astroglial Kir4.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 [KJ']O (>2.7 mM for ~45O ms for single stimulation (Fig. 4A,B), and >35 mM for 1.5 s for tetanic stimulation (Fig. 4E,F)). Nevertheless, we show a prominent and specific involvement of astroglial Kir4.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 Kir4.1 channels can contribute to fine tuning of neuronal spiking involving low, but long-lasting, increase in [K+]0. Thus besides gliotransmission, regulation of [KJ']O by Kir4.1 channel provides astrocytes with an alternative active and efficient mechanism to regulate neuronal activity. Several studies have reported decreased Kir4.1 protein levels and Kir functional currents in sclerotic hippocampus from epileptic patients [43—46]. Whether these changes are the cause or the consequence of epilepsy is still an open question. However, Kir4.1'/' mice display an epileptic phenotype [22,47] and missense mutations in KCN]10, the gene encoding Kir4.1, have been associated with epilepsy in humans [48,49]. These data thus suggest that impairment in Kir4.1 function leading to alterations in [K‘L]O dynamics, as shown in our study, may cause epilepsy. Remarkably, dysfunction of [KJr]o regulation by Kir4.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 Kir4.1 channels in glial cells was recently found in almost 50% of multiple sclerosis patients [51]. Thus astroglial Kir4.1 channels may well represent an alternative therapeutic target for several diseases.

Materials and Methods

Ethics statement

All efforts were made to minimize the number of used animals and their suffering. Experiments were performed on the hippocampus of wild type mice (C57BL6). For all analyses, mice of both genders and littermates were used (PN19—PN25).

Electrophysiological recordings

Slices were kept at room temperature (21—23°C) in a chamber filled with an artificial cerebrospinal fluid (ACSF) composed of (in mM): 119 NaCl, 2.5 KCl, 2.5 CaClz, 1.3 MgSO4, 1 NaH2P04, 26.2 NaHCO3 and 11 glucose, saturated with 95% 02 and 5% C02, prior to recording. Acute slices were placed in a recording chamber mounted on a microscope including infrared differential interference (IR-DIC) equipment, and were bathed in ACSF perfused at 1.5 ml/min. ACSF contained picrotoxin (100 M), and connections between CA1 and CA3 regions were cut to avoid epileptic-like activity propagation. Extracellular field and whole-cell patch-clamp recordings were obtained using glass pipettes made of borosilicate. Astroglial and postsynaptic responses were evoked by Schaffer collateral stimulation (0.05HZ) in the CA1 stratum radiatum region with glass pipettes filled with ACSF (300—700 k0). Astrocytes from stratum radiatum were recognized by their small soma size (5—10 pm), very low membrane resistance and hyperpolarized resting membrane potentials (w- 80 mV), passive properties of their membrane (linear IV), absence of action potentials, and large coupling through gap junctions. Field excitatory postsynaptic potentials (fEPSPs) were obtained in 400 um slices using pipettes (4—6 MQ) located in the stratum radia-tum region. Stimulus artifacts were suppressed in representative traces. Whole-cell recordings

Modeling potassium dynamics in the tripartite neuron-astrocyte-extracellular compartment

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-compart-ment 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 Kir4.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—57].

The synaptic current Iapp is proportional to the fraction of synaptic resources in the effective state e and is given by Iapp = Asee (the parameter Ase is defined in table 1). We used the following definitions for the input function

Modeling neuronal activity

The dynamics of the neuronal membrane potential, VN, follows the classic Hodgkin Huxley (HH) equations [58]. With rate equations Vrest is the resting membrane potential and VKN and VNaN are respectively the K+ and Na+ equilibrium potentials and are given by the Nernst equations Where Nao and NaN are respectively the extracellular and neuronal sodium concentrations, and K0 and KN 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 Iapp is the synaptic current derived from equation 1.

Modeling astrocytic potassium uptake by Kir4.1 channels

To account for the K+ dynamics in astrocytes, we modeled the Kir4.1 channel according to its biophysical properties [59] and IV curve [60]. The total astroglial current I Kir depends on the membrane potential, the extracellular (K0) and the astrocytic (K A) K+ concentrations, and is approximated by

1C, [60]), as detailed below. The second term of equation 22 describes the dependence of IKir to the square root of K0 [60—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 11A 2 g1A(VA—V1A), Which stabilizes the astrocyte membrane potential at—80 mV, the astrocyte membrane potential VA satisfies the equation Where I Kir is defined by relation 22. used the Matlab fitting procedure for a single exponential With formula 22 changed to (V—VA1—261n(3/145 (table 1). Varying [K‘L]0 by 0.5 mM did not affect significantly the Kir4.1 channel IV 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 ipumpk (index k = N for the neuron, k = A for the astrocyte) depend on the extracellular K+ K0 and intracellular Na+ concentrations (NaN for the neuron and Na A for the astrocyte) and follow the same equation as [66],

Balance of ionic fluxes

A current I across a membrane induces a flow of charge 1' equals to 6Q = I per unit of time. The corresponding change in extracellular concentration is given by I/(qNA V010), Where q = 1.6 >|< 10'19C is the charge of an electron, N A the Avogadro number and VolN, Vol A andVolo are the neuronal, astrocytic and extracellular volume respectively. To model the ionic concentration dynamics, we converted the currents I Na, I K and I Ki, to the corresponding ionic fluxes iNa, iK and iKl-r We describe in the following paragraphs the equations for the ionic concentrations in the three compartments (neuron, extracellular space and astrocyte).

Potassium fluxes

The extracellular K+ K0 increases with the neuronal current I K (see equation 8), which is here converted to iK (ion flux), but it is also uptaken back into neurons with a flux 2 ipumpN (the factor 2 is described in [67] and into astrocytes as the sum of the two fluxes 2 ipumPA plus 1' K1,. Similarly, we obtain the equations for the neuronal and astrocytic K+ to balance the various fluxes. Finally, we get

We thus set at zero both the Kir4.1 current and the leak term. To compensate for the loss of K+ fluxes through astroglial Kir4.1 channels, we added in equation 27 a constant K+ flux to maintain [KJF]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-l, 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

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

Numerical implementations and fitting procedures

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 At = 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 T of simulated extracellular K+ transients were fitted to curves using a single exponential (e‘i) (Fig. 4B,F,]). For all the fits obtained on the numerical simulation curves, we obtained an error estimation R-square 2 0.97. Time constants T 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 T were fitted to curves using a single exponential(e‘%). For all the fits obtained on the numerical simulation curves, we obtained an error estimation R-square 2 0.97.

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 Iapp depends on the input functions fs(t) (equation 4), fTT(t) (equation 5) and fRs(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 Iapp in control conditions (1‘ = 300 ms, 7mm 2 200 ms). The parameters for the Kir4.1 inhibition condition in the model were extracted from our experimental results on Kir4.1 glial conditional knockout mice [26] and are given by rm. 2 500 ms, 7mm 2 160 ms. When Kir4.1 channels are inhibited (model) or knocked-out (experiment), the maximal peak amplitudes of the applied synaptic currents in the model (lapp) and fEPSPs recorded experimentally are increased compared to control conditions [26].

Simulation of neuronal firing at various frequencies

Each input is generated by a sub-firing square current lasting 5 ms (lapp). pyramidalIn addition, we added a Brownian noise of amplitude o = 0.68 pA2 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 py-ramidal 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 tto the total number of simulations.

Supporting Information

AC) Astroglial net K+ uptake (equation 27) and activity-dependent changes in [KL] 0 evoked by single (light blue), tetanic (100 HZ, 1s, green) and repetitive (10 HZ, 30 3, dark blue) stimulations are plotted as a function of time. D) Phase diagram illustrating astroglial K+ uptake as a dynamic function of activity-dependent changes in [KJr]O evoked by the different stimulations. Astroglial K+ net uptake is normalized to the maximum value obtained during repetitive stimulation. (EPS) Acknowledgments We thank O. Chever and U. Pannasch for helpful discussions.

Author Contributions

Performed the experiments: IS KDD. Analyzed the data: IS KDD NR DH. Wrote the paper: IS KDD NR DH.