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.
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–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 [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–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 (Kir4. 1 and K2P, respectively) or Na/K ATPases [18,19]. Remarkably, recent work suggest that Kir4. 1 channels play a prominent role in astroglial regulation of [K+]o [20–23]. However, the mouse model used to draw these conclusions, i. e. conditional Kir4. 1 knockout mice directed to glial cells (GFAP-Cre-Kir4. 1fl/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 Müller 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.
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 Kir4. 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 (Iapp, 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 Kir4. 1 channels and Na/K ATPases (Fig. 1B and Materials and Methods). Because Kir4. 1 channels are strongly involved in K+ uptake [22], we fitted the I-V curve of K+ ions through Kir4. 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 Kir4. 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 Kir4. 1 channels to K+ uptake in relation to neuronal activity associated with different [K+]o. 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 (Iapp) 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), tetanic stimulation (rise time: 610 ms for numerical simulation, 491 ± 122 ms n = 5 for experiments; time of peak: 1. 07 s for numerical simulation, 1. 05 s ± 0. 25 n = 5 for experiments; decay time: 4. 18 s for numerical simulation, 4. 55 s ± 0. 45 n = 5 for experiments, Fig. 2H) and repetitive stimulation (rise time: 1. 5 s for numerical simulation, 1. 27 s ± 0. 18 n = 5 for experiments; time of peak: 6. 8 s for numerical simulation, 5. 2 s ± 0. 9 n = 5 for experiments; decay time: 7. 95 s for numerical simulation, 8. 3 s ± 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. 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. 2 s; tetanic stimulation: 8. 7 s; repetitive stimulation: 34. 2 s, Fig. 3C, F, I). The astroglial net K+ uptake increases with the activity-dependent [K+]o transient rises (S2A–C Fig.) evoked by the different regimes (S2D Fig.). Neurons slowly re-uptake only ~5–10% of their released K+ at the end of phase 1 (Fig. 3C, F, I). Remarkably, although [K+]o increases with the strength of stimulation (from 0. 9 to 4. 4 mM, Fig. 3A, D, G and S2A–C Fig.), the time needed for astrocytes to buffer the released K+ is not proportional to [K+]o rises (Fig. 3C, F, I), as shown by the phase diagram illustrating astroglial K+ uptake as a dynamic function of activity-dependent changes in [K+]o evoked by the different stimulations (S2D Fig.), but is to the square root of [K+]o (equation 22). In addition, at the end of phase 1, [K+]o is almost back to baseline levels, whereas intra-astroglial K+ levels reach their peak value (Fig. 3C, F, I). 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 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. To study quantitatively the acute and selective role of astroglial Kir4. 1 channels in neuroglial K+ dynamics, we inhibited the Kir4. 1 current in our tri-compartment model. Because Kir4. 1-/- mice display altered synaptic plasticity compared to wild type mice [22,26], we recalibrated the synaptic current (Iapp) 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 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 [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 Kir4. 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 addition, for all regimes of activity, the rise and decay times of the [K+]o 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, J). Finally, Kir4. 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 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 firing during trains of stimulations. To investigate the effect of astroglial Kir4. 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 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. 5J). These data suggest a prominent and specific involvement of astroglial Kir4. 1 channels in regulation of firing during theta rhythmic 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+]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 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. 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 Kir4. 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 Kir4. 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 (K2P) (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 (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 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 K2P 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 Kir4. 1 channels were chronically deleted genetically in glial cells [20,21,23], we found that acute inhibition of Kir4. 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 Kir4. 1 channel inhibition also alters significantly [K+]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]. 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 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 CA ≈ 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 τ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 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 (VA−VKA−V2AV3A) 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 [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 Kir4. 1 channels have a strong uptake capacity, especially for high regimes of activity ([K+]o up to 5–6 mM). Remarkably, our model reveals that astroglial Kir4. 1 channels strongly regulate neuronal firing induced by high stimulation regime such as repetitive stimulation. Kir4. 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 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 [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 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+]o. Thus besides gliotransmission, regulation of [K+]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 KCNJ10, 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+]o dynamics, as shown in our study, may cause epilepsy. Remarkably, dysfunction of [K+]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.
Experiments were carried out according to the guidelines of European Community Council Directives of January 1st 2013 (2010/63/EU) and our local animal committee (Center for Interdisciplinary Research in Biology in Collège de France). 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). Acute transverse hippocampal slices (400 μm) were prepared as previously described [42,52–54] from 19–25 days-old wild type mice. 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 CaCl2,1. 3 MgSO4,1 NaH2PO4,26. 2 NaHCO3 and 11 glucose, saturated with 95% O2 and 5% CO2, prior to recording. Acute slices were placed in a recording chamber mounted on a microscope including infra-red 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 kΩ). Astrocytes from stratum radiatum were recognized by their small soma size (5–10 μm), very low membrane resistance and hyperpolarized resting membrane potentials (≈- 80 mV), passive properties of their membrane (linear I-V), absence of action potentials, and large coupling through gap junctions. Field excitatory postsynaptic potentials (fEPSPs) were obtained in 400 μm slices using pipettes (4–6 MΩ) located in the stratum radiatum region. Stimulus artifacts were suppressed in representative traces. Whole-cell recordings were obtained from CA1 astrocytes, using 4–6 MΩ glass pipettes containing (in mM): 105 K-Gluconate, 30 KCl, 10 HEPES, 10 Phosphocreatine, 4 ATP-Mg, 0. 3 GTP-Tris, 0. 3 EGTA (pH 7. 4,280 mOsm). Prolonged repetitive stimulation was performed for 30 s at 10 Hz. Post-tetanic potentiation was evoked by stimulation at 100 Hz for 1 s in the presence of 10 μM CPP ( (Rs) -3- (2-Carboxypiperazin-4-yl-) propyl-1-phosphonic acid). Recordings were performed with Axopatch-1D amplifiers (Molecular Devices, USA), at 10 kHz, filtered at 2 kHz, and analyzed using Clampex (Molecular Devices, USA), and Matlab (MathWorks, USA) softwares. The data represent mean ± SEM. Picrotoxin was from Sigma and CPP from Tocris. We present here the biophysical model we have built to describe K+ dynamics during neuronal activity and specifically the role of astroglial Kir4. 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 Kir4. 1 channels. 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]. where f is the input function. For a single stimulation generated at time tstim, f (t) = δ (t-tstim). A stimulation instantaneously activates a fraction Use 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 = tstim, r and e respectively decreases and increases by the value User. 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 f: The dynamics of the neuronal membrane potential, VN, follows the classic Hodgkin Huxley (HH) equations [58]. with rate equations αn (VN) =0. 01 (VN+10) exp (0. 1 (VN+10) ) −1 (12) βn (VN) =0. 125exp (VN/80) (13) αm (VN) =0. 1 (VN+25) exp (0. 1 (VN+25) ) −1 (14) βm (VN) =4exp (VN/18) (15) αh (VN) =0. 07exp (VN/20) (16) βh (VN) =1exp (0. 1 (VN+ 30) ) +1 (17) 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 VNaN=RTFln (Na0NaN) (18) VKN=RTFln (K0KN) (19) where Na0 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 IlN=glN (VN−VlN) (20) which stabilizes the membrane potential at its resting value. Finally, the neuronal membrane potential satisfies the equation CNdVNdt=− (INa+IK+IlN+Iapp) (21) where Iapp is the synaptic current derived from equation 1. To account for the K+ dynamics in astrocytes, we modeled the Kir4. 1 channel according to its biophysical properties [59] and I-V curve [60]. The total astroglial current IKir depends on the membrane potential, the extracellular (K0) and the astrocytic (KA) K+ concentrations, and is approximated by IKir=GKir (VA−VKA−VA1) (K01+exp (VA−VKA−VA2VA3) ) (22) where VKA is the Nernst astrocyte K+ potential, VA, the astrocyte membrane potential, K0 is the extracellular K+ concentration and VA1 (an equilibrium parameter, which sets Kir current to 0 at-80 mV), VA2 and VA3 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 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 IlA = glA (VA—VlA), which stabilizes the astrocyte membrane potential at—80 mV, the astrocyte membrane potential VA satisfies the equation CAdVAdt=− (IKir+IlA) (23) where IKir is defined by relation 22. We fitted the Kir4. 1 channel I-V curve (equation 22) using the experimental recordings for the Kir4. 1 channel (3 mM [K+] (Fig. 4 in [60,65]). We first obtained that VA1 = (VrestA − 26ln (3/145) ) = −14. 83 mV where VrestA = −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−VA1−26ln (3/145) ) 3I with (V from-100 to 20 mV) to get that VA2 = 34 mV and VA3 = 19. 23 mV (table 1). Varying [K+]o by 0. 5 mM did not affect significantly the Kir4. 1 channel I-V curve, confirming its robustness. 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 ipump, k (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 NaA for the astrocyte) and follow the same equation as [66], ipump, k=imaxk (1+7. 3K0) −2 (1+10Nak) −3 for k=N, A (24) where imaxk is a constant (table 1). 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/ (qNAVol0), where q = 1. 6 * 10–19C is the charge of an electron, NA the Avogadro number and VolN, VolA andVol0 are the neuronal, astrocytic and extracellular volume respectively. To model the ionic concentration dynamics, we converted the currents INa, IK and IKir to the corresponding ionic fluxes iNa, iK and iKir We describe in the following paragraphs the equations for the ionic concentrations in the three compartments (neuron, extracellular space and astrocyte). To determine the system of equations for the K+ fluxes, we use the mass conservation law for the extracellular K0, the neuronal KN and the astrocytic KA K+ concentrations. The extracellular K+ K0 increases with the neuronal current IK (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 iKir. 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 Kir4. 1 channels in neuroglial K+ dynamics, we inhibited the Kir4. 1 current in our tri-compartment model. 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 [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. 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 iNalA and iNalN (values given in table 1), as classically used [24], dNa0dt=iNa+ iNalN+3ipumpN+3ipumpA+ iNalA (28) dNaNdt= (−iNa−3ipumpN− iNalN) Volo VolN (29) dNaAdt= (−iNalA−3ipumpA) Volo VolA (30) 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τ) (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τ). 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 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 (τ = 300 ms, τinact = 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 τrec = 500 ms, τinact = 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 (Iapp) and fEPSPs recorded experimentally are increased compared to control conditions [26]. 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 (Iapp). In addition, we added a Brownian noise of amplitude σ = 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 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.
Neural excitability relies on precise inward and outward ionic fluxes. 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 moment-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.
The Neuroglial Potassium Cycle during Neurotransmission: Role of Kir4.1 Channels