Biophysical models reveal the relative importance of transporter proteins and impermeant anions in chloride homeostasis
Abstract
Fast synaptic inhibition in the nervous system depends on the transmembrane flux of Cl- ions based on the neuronal Cl- driving force. Established theories regarding the determinants of Cl- driving force have recently been questioned. Here, we present biophysical models of Cl- homeostasis using the pump-leak model. Using numerical and novel analytic solutions, we demonstrate that the Na+/K+-ATPase, ion conductances, impermeant anions, electrodiffusion, water fluxes and cation-chloride cotransporters (CCCs) play roles in setting the Cl- driving force. Our models, together with experimental validation, show that while impermeant anions can contribute to setting [Cl-]i in neurons, they have a negligible effect on the driving force for Cl- locally and cell-wide. In contrast, we demonstrate that CCCs are well-suited for modulating Cl- driving force and hence inhibitory signaling in neurons. Our findings reconcile recent experimental findings and provide a framework for understanding the interplay of different chloride regulatory processes in neurons.
https://doi.org/10.7554/eLife.39575.001eLife digest
Cells called neurons in the brain communicate by triggering or inhibiting electrical activity in other neurons. To inhibit electrical activity, a signal from one neuron usually triggers specific receptors on the second neuron to open, which allows particles called chloride ions to flow into or out of the neuron.
The force that moves chloride ions (the so-called ‘chloride driving force’) depends on two main factors. Firstly, chloride ions, like other particles, tend to move from an area where they are plentiful to areas where they are less abundant. Secondly, chloride ions are negatively charged and are therefore attracted to areas where the net charge (determined by the mix of positively and negatively charged particles) is more positive than their current position.
It was previously believed that a group of proteins known as CCCs, which transport chloride ions and positive ions together across the membranes surrounding cells, sets the chloride driving force. However, it has recently been suggested that negatively charged ions that are unable to cross the membrane (or ‘impermeant anions’ for short) may set the driving force instead by contributing to the net charge across the membrane. Düsterwald et al. used a computational model of the neuron to explore these two possibilities.
In the simulations, altering the activity of the CCCs led to big changes in the chloride driving force. Changing the levels of impermeant anions altered the volume of cells, but did not drive changes in the chloride driving force. This was because the flow of chloride ions across the membrane led to a compensatory change in the net charge across the membrane.
Düsterwald et al. then used an experimental technique called patch-clamping in mice and rats to confirm the model’s predictions. Defects in controlling the chloride driving force in brain cells have been linked with epilepsy, stroke and other neurological diseases. Therefore, a better knowledge of these mechanisms may in future help to identify the best targets for drugs to treat such conditions.
https://doi.org/10.7554/eLife.39575.002Introduction
Fast synaptic inhibition in the nervous system is mediated by type A γ-aminobutyric acid receptors (GABAARs) and glycine receptors (GlyRs), which are primarily permeable to chloride (Cl-) (Farrant and Kaila, 2007). Together with the neuronal membrane potential, the transmembrane gradient for Cl- sets the driving force for Cl- flux across these receptors, controlling the properties of inhibitory signaling. Modification of neuronal intracellular Cl- concentration ([Cl-]i) has been shown to play a causative role in multiple neurological diseases including epilepsy, chronic pain, schizophrenia and autism (Rivera et al., 2004; Huberfeld et al., 2007; Price et al., 2009; Hyde et al., 2011; Tyzio et al., 2014). Similarly, intracellular Cl- is thought to be modulated during brain development so that GABAergic transmission contributes optimally to the construction of neural circuits (Ben-Ari, 2002). Given the importance of Cl- for brain function and dysfunction, the cellular mechanisms that control its transmembrane gradient and driving force are of considerable interest.
Plasmalemmal Cl- transporters, in particular cation-chloride cotransporters (CCCs), are understood to be the major mechanism by which neurons regulate the driving force for Cl- permeable anion channels (Kaila et al., 2014). Recently, it has been suggested that in fact impermeant anions control local [Cl-]i and driving force (Glykys et al., 2014), rather than the CCCs. The majority of intracellular anions are impermeant to the neuronal membrane; these include ribo- and deoxynucleotides, intracellular proteins and metabolites (Burton, 1983). Impermeant anions induce what is known as the Donnan (or Gibbs-Donnan) effect (Hill, 1956; Sperelakis, 2012) – an uneven distribution of impermeant molecules across the membrane which is osmotically unstable. Without active ion transport to counter this effect, neurons would swell and burst (Kay, 2017). Animal cells, including neurons, maintain cell volume in the presence of impermeant anions by using the Na+/K+-ATPase to pump Na+ out of the cell and K+ in, along with the passive movement of water and other ions (Tosteson and Hoffman, 1960; Armstrong, 2003; Liang et al., 2007; Kay, 2017). This pump-leak mechanism, whilst stabilizing cell volume, also establishes the negative resting membrane potential and transmembrane Na+ and K+ gradients, which serve as energy sources for the coupled transport of other molecules, including Cl- by CCCs. In the absence of active Cl- transport, [Cl-]i is set by the membrane potential; i.e. the Nernst potential of Cl- () equals the transmembrane potential (Vm).
The transmembrane Cl- gradient and the driving force (DF = Vm−ECl) for Cl- permeable ion channels are therefore the outcome of multiple, dynamically interacting mechanisms. This makes experimental investigation of the determinants of Cl- driving force difficult, particularly at a local level. Computational models based on established biophysical first principles are a productive means for exploring the roles of cellular mechanisms in generating local Cl- driving force. Here, we establish numerical and novel analytic solutions for an inclusive model of Cl- homeostasis to elucidate the determinants of the neuronal driving force for Cl-.
We demonstrate that baseline [Cl-]i is a product of the interaction of Na+/K+-ATPase activity, the mean charge of impermeant anions, ion conductances, CCCs and water permeability. Consistent with recent experimental reports (Glykys et al., 2014), and our own experimental validation using electroporation of anionic dextrans and optogenetic probing of EGABA, we find that impermeant anions can contribute to setting [Cl-]i. However, we find that they can only affect the Cl- driving force by modifying active transport mechanisms, and then only negligibly. Impermeant anions therefore do not appreciably modify synaptic signaling properties, contrary to the interpretation of recent experiments (Glykys et al., 2014). In contrast, we demonstrate using biophysical models and gramicidin perforated-patch clamp recordings that CCCs selectively regulate substantial changes in the Cl- driving force. This is consistent with a meta-analysis of experimental data from the field, which shows a strong correlation between the activity of the specific CCC, KCC2, and Cl- driving force. The ability of CCCs to specifically modulate Cl- at a local level depends on the characteristics of Cl- electrodiffusion in the structure concerned, demonstrated using multicompartment modeling. Together, our models provide a theoretical framework for understanding the interplay of chloride regulatory processes in neurons and interpreting experimental findings.
Results
A biophysical model based on the pump-leak mechanism demonstrates the importance of the sodium-potassium ATPase for setting transmembrane ion gradients including chloride
To compare the effects of impermeant anions and Cl- cotransport on Cl- homeostasis we first developed a single compartment model based on the pump-leak formulation (Tosteson and Hoffman, 1960; Kay, 2017) (Figure 1A). This model, defined by a set of differential equations, incorporated mathematical representations of the three major permeable ion species Cl-, K+, Na+ as well as impermeant ions (Xz) with mean charge z. Permeable ions could move across the cellular membrane via passive conductances according to each ion’s respective electrochemical gradient. Further, active transport of Na+ and K+ by the Na+/K+-ATPase (with a 3:2 stoichiometry) and cotransport of Cl- and K+ (1:1 stoichiometry) by the cation-chloride cotransporter KCC2 with a non-zero conductance gKCC2 of 20 µS/cm2 unless otherwise stated were included. Finally, our formulation accounted for the dynamics of cell volume (w), intracellular osmolarity () and transmembrane voltage (Vm) (see Materials and methods). Importantly, regardless of initial starting concentrations of permeant or impermeant ions, cell volume or Vm, the model converged to stable fixed points without needing to include any means for ‘sensing’ ion concentration, volume or voltage. Initial permeable ion concentrations also did not influence the final cellular volume. For example, despite initiating the model with different starting concentrations of Cl- (1, 15, 40 and 60 mM, respectively, Figure 1B), [Cl-]i always converged to the same stable concentration of 5.2 mM, a typical baseline [Cl-]i for adult neurons, and volume always converged to 2.0 pL, a typical volume for hippocampal neurons (Ambros-Ingerson and Holmes, 2005). The model is robust in the sense that its convergence to a stable steady state does not depend on a narrow set of parameters and initial values. Trying an alternative model of the Na+/K+-ATPase (Hamada et al., 2003) produced similar results (Figure 1—figure supplement 1A-B).
Consistent with previous results (Xiao et al., 2002; Dierkes et al., 2006; Dijkstra et al., 2016), ‘turning off’ the activity of the Na+/K+-ATPase in our model led to a progressive collapse of transmembrane ion gradients, progressive membrane depolarization and continuous and unstable cell swelling. Such effects could be reversed by reactivation of the pump (Figure 1C). The relative activity of the Na+/K+-ATPase sets the final stable values for ion concentrations (of both permeable and impermeant ions), Vm and volume (Figure 1C). When we increased the activity of the Na+/K+-ATPase in our model, the final steady-state concentration for K+ increased, whilst Na+ and Cl- dropped to levels that approximate those observed in mature neurons (Figure 1D). Note that although sufficient Na+/K+-ATPase activity is critical for steady state ionic gradients including that of Cl-, these are relatively stable near the default pump rate (Figure 1—figure supplement 2). At the same time, the final, stable-state membrane potential and cell volume also decreased. Interestingly, as has been observed previously (Fraser and Huang, 2004), beyond a certain level, further increases in Na+/K+-ATPase activity have negligible effects on cell volume and transmembrane voltage. For subsequent analysis, we chose a ‘default’ effective pump rate for the Na+/K+-ATPase of approximately 1.0 × 10−2 C/(dm2.s), that is a pump rate constant of 10−1 C/(dm2.s) (equation (2)), and mean intracellular impermeant anion charge (z) of −0.85 as extrapolated from reasonable cellular ionic concentrations and osmolarity (Lodish et al., 2009; Raimondo et al., 2015). This resulted in steady-state ion concentrations and membrane potentials that approximate those experimentally observed in mature neurons: Cl- 5 mM; K+ 123 mM; Na+ 14 mM; Xz 155 mM; and a Vm of −72.6 mV (Jiang and Haddad, 1991; Diarra et al., 2001; Tyzio et al., 2008). We were able to corroborate the numerical solutions for final steady-state values by developing a parametric-analytic solution (Supplementary file 1). We observed exact correspondence between the numerical and analytic solutions within our model (Figure 1D). In subsequent analyses, this novel analytic solution allowed us to explore rapidly a large parameter space to determine how various cellular attributes might affect Cl- homeostasis.
Membrane chloride conductance affects steady-state intracellular chloride concentration only in the presence of cation-chloride cotransport
Using the analytic solution, we investigated how changes in baseline ion conductance for the major ions in our model (gK, gNa and gCl) affected Cl- homeostasis. We calculated the steady-state values for the Cl- reversal potential (ECl) and K+ reversal potential (EK), resting membrane potential (Vm) and volume (w) whilst independently manipulating the conductance for each ion (Figure 2). Increasing the baseline K+ conductance (gK) resulted in ECl, EK and Vm converging to similar steady-state values (Figure 2A) without significantly affecting cell volume. We were also able to replicate the classic dependence of membrane potential on log([K+]o) (Figure 2—figure supplement 1). In contrast, increasing the baseline Na+ conductance (gNa) beyond 20 µS/cm2 resulted in a steady increase of ECl, Vm and volume with a minimal increase of EK (Figure 2B). EK is maintained in the face of increased passive K+ efflux accompanying membrane depolarization due to increased active influx of K+ by the Na+-dependent ATPase which increases its effective pump rate due to increased intracellular Na+ concentration with larger gNa (Figure 2—figure supplement 2).
The effect of manipulating Cl- conductance (gCl) depended on the activity of concurrent cation-chloride cotransport by KCC2 (Figure 2C). In the presence of active KCC2 at very low values of gCl, the steady state [Cl-]i is such that ECl approaches EK. This follows because in the absence of alternative Cl- fluxes, KCC2 utilizes the transmembrane K+ gradient to transport Cl- until ECl equals EK. With increasing gCl however, ECl increases, moving away from EK toward Vm, and at very high Cl- conductances ECl and Vm approached similar values in our model. Without the activity of KCC2, any non-zero gCl had no effect on steady state ECl, EK, Vm or volume (Figure 2D). In this instance ECl always equals Vm as the movement of Cl- across the membrane is purely passive. Without the activity of KCC2, there can be no driving force for Cl- flux at steady state (Vm-ECl = 0). Our model therefore behaved in a manner consistent with established theoretical predictions (Kaila et al., 2014).
Cation-chloride cotransport sets the chloride reversal and driving force for transmembrane chloride flux
Next, we used our single-cell unified model to explore how the activity of cation-chloride cotransport affects Cl- homeostasis. In our model, the activity of KCC2 is set by the conductance of KCC2 (gKCC2). Using the numerical formulation with the default values described in Figures 1 and 2, we steadily increased gKCC2 from 20 µS/cm2 to 370 µS/cm2 and tracked changes to ECl, EK, Vm and volume. Increasing KCC2 activity over time caused a steady decrease in [Cl-]i reflected by a hyperpolarization of ECl (Figure 3A). Vm decreased only modestly, resulting in an increase in the driving force for Cl- flux that tracks the increase in gKCC2. This effect saturates as EK constitutes a lower bound on ECl. Importantly, increases in gKCC2 resulted in persistent changes to ECl and the driving force for Cl-. Employing alternate models for KCC2 (Fraser and Huang, 2004; Lewin et al., 2012; Raimondo et al., 2012) and the Na+/K+-ATPase (Hamada et al., 2003) did not change this result when compensation for parameterization was given, although different KCC2 models result in different kinetic rates for Cl- and K+ transport (Figure 3—figure supplements 1 and 2). Using the analytic solution to our model, we calculated how KCC2 activity affects steady state values of ECl, EK, Vm, volume and Cl- driving force (Figure 3B). In confirmation of our findings in Figure 2D, with no KCC2 activity (gKCC2 = 0), ECl equaled Vm and the Cl- driving force was zero. As we increased gKCC2, steady state ECl pulled away from Vm and approached EK. This resulted in an increase in Cl- driving force (Vm-ECl) with steady state values of 11.3 mV at our chosen default value of gKCC2. The results obtained with our model are therefore fully consistent with the view that CCCs, in this case KCC2, establish the driving force for Cl-.
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.39575.012
To test this theoretical finding, we performed gramicidin perforated patch-clamp recordings from CA3 hippocampal neurons in rat organotypic brain slices whilst activating Cl- permeable GABAA receptors with muscimol (10 µM), in order to measure the GABAAR driving force, which approximates Cl- driving force (Figure 3C). We then tested our model predications by applying the CCC blocker, furosemide (1 mM) (Figure 3D,E). We noted that after furosemide was introduced the EGABA became significantly more depolarized (baseline: −78.8 ± 2.8 mV vs furosemide: −71.6 ± 3.0 mV, n = 10, p=0.01, paired t-test), whilst there was no significant difference in Vm (−69.9 ± 1.8 mV vs −71.5 ± 2.6 mV, n = 10, p=0.36, paired t-test) (Figure 3E). This reflects a significant change in the GABAAR driving force (8.9 ± 3.4 mV vs 0.1 ± 4.6 mV, n = 10, p=0.04, paired t-test). Values were stable prior to baseline recordings with EGABA 5 min prior to furosemide application at −78.4 ± 8.1 mV, Vm−70.1 ± 5.6 mV and DF at 9.1 ± 3.2 (see Figure 3E). We then noted that this change in driving force persisted for at least 15 min post-application of furosemide (Figure 3E). These results are consistent with our model predictions and demonstrate how the application of a CCC blocker reduces the GABAAR driving force (and hence Cl- driving force) by selectively depolarizing EGABA with negligible effects on Vm.
In addition, we sought experimental data from the literature to determine whether changes in KCC2 activity correlate with alterations to steady-state [Cl-]i. We focused on changes in KCC2 expression level, as this is likely to be a strong predictor of changes in KCC2 activity. Indeed, in a meta-analysis of seven studies and eight experiments from our review of 26 studies, weighted for methodological biases and data quality, we observed a significant correlation (R2 = 0.796, p<0.001) between the change in KCC2 expression and Cl- driving force (Figure 3F). Absolute changes in Vm were less than 2 mV in all but one study, meaning that the change in driving force could be ascribed to significant shifts in EGABA (R2 = 0.045, p<0.001). The outlier data point (showing a 8.45 mV change) was from a study into the effects of acute stress, where other factors could have transiently influenced Vm (MacKenzie and Maguire, 2015). The meta-analysis supports the prediction that cation-chloride cotransport by KCC2 is an important determinant of [Cl-]i (R2 = 0.83, p<0.001, nine studies) and driving force (see Supplementary file 2 Table S2-1 for raw data, and the scoring table for weighting in Table S2-2).
Altering the concentration of intracellular or extracellular impermeant anions, without changing the mean charge of impermeant anions, does not affect the steady state gradient or driving force for chloride
To determine the effect of impermeant anions on Cl- homeostasis, we first explored whether adjusting the concentration of impermeant anions ([X]i), while maintaining a constant mean impermeant ion charge (z), had any impact on ECl, EK, Vm or volume. The mean charge (z) is the mean of the charge of all the different species of impermeant molecules in the cell, including uncharged ones, where charge is the difference between the number of protons and electrons of a molecule. Impermeant anions are more abundant than impermeant cations, and so in this manuscript we often refer to the group as impermeant anions rather than impermeant ions or impermeant molecules. For example, were there α impermeant molecules of charge −1 and β impermeant molecules of charge 0, then z would be . We initiated the full single-compartment model with different starting concentrations of impermeant anions all with the same mean charge, z = −0.85, and observed that regardless of the initial concentration of impermeant anions, over a period of minutes, the cell adjusted its volume to give an identical steady-state impermeant anion concentration (Figure 4A, [A]i = 155 mM). Analytically, it can be shown that the number of moles of X determines completely the volume of the compartment, while the permeant ions alone cannot be used to predict steady state volume (Kay, 2017). Similarly, all initial impermeant anion concentrations resulted in identical steady state values of ECl (−83.8 mV), EK (−95.1 mV) and Vm (−72.6 mV) (Figure 4B). This shows that simply adjusting the amount of impermeant anions within a cell has no persistent effect on [Cl-]i.
We then tested the effect of dynamically adding impermeant anions with the default mean charge either intracellularly (Figure 4C) or extracellularly (Figure 4D). While impermeant anions are being added to the cell, the membrane potential hyperpolarizes and ECl decreases. However, following the cessation of impermeant anion influx, ECl, EK, Vm and [X]i return to steady state values due to compensatory changes to cell volume (Figure 4C). There are transient transmembrane fluxes of all ions while anions are added into the cell, and in particular the inward flux of the cations Na+ and K+, such that the sum [X]i + [Cl-]i is not necessarily kept constant during impermeant anion addition (Figure 4—figure supplement 1). Impermeant ions (with the default mean charge) were added to the extracellular space, which is effectively an infinite bath in the model, while proportional decreases in [Cl-]o were applied to correct for the changes to charge and osmotic balance. Additions in the extracellular space, similarly, resulted in a temporary depolarization of ECl and Vm, but no persistent shift in these parameters (Figure 4D). The addition of extracellular impermeant anions did however result in a small compensatory decrease in cell volume secondary to the large shifts in [Cl-]i required to maintain the proportion of [Cl-]i to [Cl-]o according to the Nernst potential. In summary, there is no lasting effect on the reversal potential or driving force for Cl- if only the concentration of a neuron’s intracellular or extracellular impermeant anions is altered. This is because concentration changes alone modulate only osmoneutrality, whereas changes to intracellular charge balance affect electroneutrality and therefore the membrane and ionic potentials, which we tested next.
Changing the mean charge of impermeant anions can drive substantial shifts in the reversal potential for chloride, but has negligible effects on chloride driving force
We next sought to determine how changes in the mean charge of the impermeant ions (z) might influence the driving force for Cl-. Such changes in z could be associated with various cellular processes, including post-translational modifications of proteins that decrease their charge without changing the absolute number of protein molecules. To investigate this parameter, we modified the mean charge (z) of intracellular impermeant anions from −0.85 to −1 whilst measuring accompanying changes in ECl, EK, Vm and cell (Figure 5A). We found that this shift to a more negative z resulted in both a transient and persistent decrease in ECl, EK and V. Importantly, the shifts in ECl were accompanied by broadly matching shifts in EK and V, which resulted in a small change in the driving force for Cl- of <0.2 mV. Both numerical and analytic calculation of steady state values for ECl, EK and Vm in our model showed that changing the mean charge of impermeant anions, while substantially affecting ECl, had very small effects on the driving force for Cl- (Figure 5B). By shifting z within reasonable ranges for mammalian neurons (Lodish et al., 2009; Raimondo et al., 2015), and assuming osmo- and electro-neutrality, only shifts of <1 mV could be generated. In addition, although the absolute number of impermeant anions (moles) remained constant throughout the process of modifying z, cell volume shifted, and as a consequence modest alterations to the concentration of impermeant anions occurred as well.
-
Figure 5—source data 1
- https://doi.org/10.7554/eLife.39575.017
Next, instead of adjusting the charge of some of the intracellular impermeant anions as described above, we directly added new impermeant anions to the cell, which had a more negative charge than the previous mean charge (Figure 5C and D). This had the effect of both increasing the absolute quantity of impermeant anions and adjusting the mean charge of impermeant anions. The ‘addition’ of impermeant anions in this way models the de novo synthesis of impermeant anion species, or their active transport into the cell. This process also resulted in both transient and persistent changes to ECl, EK and V, which was dependent on the extent that z was altered. Again, whilst the large additions of impermeant anions could substantially alter the Cl- reversal potential, this had a negligible effect on the driving force for Cl- due to matching shifts in Vm. Driving shifts in ECl in this manner also resulted in changes to cell volume (Figure 5—figure supplement 1).
To experimentally test our biophysical modeling predictions, we used photo-activation of ChR2 expressing GABAergic interneurons and whole-cell patch clamp recordings of mouse organotypic CA3 hippocampal pyramidal neurons to measure Vm, EGABA, and GABA driving force (which approximates Cl- driving force) before and after addition of impermeant anions (Figure 5D). To add impermeant anions to the recorded cell, we used single-cell electroporation of fluorescently tagged anionic dextrans (Alexa Flour 488, see Materials and methods). Successful addition of impermeant anions could be confirmed visually by observing strong and stable fluorescence of the anionic dextran restricted to the recorded cell (Figure 5D,E). To drive impermeant anions into the cell, negative voltage steps (20 ms, 0.5–1 V) were applied to the electroporation pipette necessarily resulting in direct membrane depolarization, which recovered over a period of 1–5 min (Figure 5D). Once this acute perturbation had settled, as predicted by our model, changing the mean charge of impermeant anions, following the addition of highly negatively charged dextrans to the cell, resulted in a stable, mean negative shift in Vm from a baseline of −67.0 ± 4.0 mV to −74.0 ± 3.1 mV (n = 6, p=0.03, Wilcoxon test, Figure 5E,F). Again in line with our predictions, addition of impermeant anions also resulted in a significant reduction of resting EGABA (a proxy for ECl) from baseline values of −72.5 ± 2.0 mV to −77.3 ± 1.4 mV (p=0.03, Wilcoxon test). Importantly, however, similar shifts in Vm and EGABA resulted in an undetectable shift in GABA and hence Cl- driving force (5.5 ± 4.4 mV vs 3.3 ± 3.5 mV, n = 6, p=0.22, Wilcoxon test, Figure 5E,F).
Our single-compartment model of Cl- homeostasis, in conjunction with experimental validation, demonstrates that whilst the adjustment of mean impermeant anion charge can significantly affect ECl, this results in negligible changes to the driving force for Cl-. This contrasts with the results shown earlier, where adjusting the activity of cation-chloride cotransport modulates both ECl and the driving force for Cl- substantially.
Impermeant anions drive small shifts in chloride driving force by modifying the sodium-potassium-ATPase pump rate under conditions of active chloride cotransport
We next set out to determine how, and under what conditions, the modification of impermeant anions could potentially generate the very small persistent shifts in Cl- driving force we observed in our models. Due to their small size (<1 mV) these were not detectable during the experimental validation. First, we repeated the simulation performed in Figure 5A by changing the mean charge of impermeant anions in the cell, but under conditions where the Na+/K+-ATPase effective pump rate (Jp) was either a cubic function of the transmembrane Na+ gradient (default condition) or was fixed at a constant value (Figure 6A). In the case where the pump rate was fixed, adjusting the mean charge of impermeant anions generated no persistent change in Cl- driving force (Figure 6A). Modifying impermeant anions caused a significant change in steady-state intracellular Na+ concentration when Jp was kept constant. However, small shifts in Cl- driving force occurred only when the effective pump rate was variable, in which minor changes to [Na+]i caused significant changes to Jp, which in turn resulted in a small shift in Cl- driving force. There is a direct relationship between the mean charge of impermeant anions (z), [Na+]i, the effective Na+/K+-ATPase pump rate and Cl- driving force. This relationship was abolished when the effective Na+/K+-ATPase pump rate was held constant by removing its dependence on Na+ (Figure 6B). In addition, even large variations in effective pump rate near the default value caused negligible shifts in Cl- driving force of <1 mV. These results were similar when using the experimentally matched ATPase model by Hamada et al. (2003), with slight differences in final values (Figure 6—figure supplement 1A). These small, impermeant anion driven, Na+/K+-ATPase pump-dependent shifts in Cl- driving force are completely dependent on the presence of cation-chloride cotransport. In the absence of KCC2, there is no Cl- driving force as ECl = Vm (Figure 6F).
We then tested whether relaxing the condition of transmembrane osmoneutrality might also alter impermeant anion induced effects on Cl- driving force. We modeled a situation where increases in cell surface area beyond a certain ‘resting’ surface area generated a hydrostatic pressure (membrane tension), which could balance an osmotic pressure difference of 10 mM between the intra- and extracellular compartments (see Figure 6C, schematic in Figure 6E and Materials and methods). In this case, adding impermeant anions of default charge z resulted in constrained increases in cell volume, which were accompanied by persistent transmembrane differences in osmolarity and intracellular Na+ concentration. This was sufficient to generate small differences in driving force for Cl- of <0.2 mV for reasonable increases in cell surface area (Nichol and Hutter, 1996; Dai et al., 1998). Again, this was entirely due to Na+ driven shifts in the Na+/K+-ATPase effective pump rate. By removing the dependence of Na+/K+-ATPase activity on Na+ concentration, addition of impermeant anions no longer generated persistent shifts in Cl- driving force (Figure 6C). Using the experimentally-matched ATPase model by Hamada et al. (2003) generated similar results (Figure 6—figure supplement 1B) because the model is also directly dependent on [Na+]i. We observed a direct relationship between transmembrane osmotic gradient, [Na+]i, the effective Na+/K+-ATPase pump rate and Cl- driving force. This relationship was removed when the effective Na+/K+-ATPase pump rate was held constant, with no changes in Cl- driving force seen despite the generation of the same shift in the transmembrane osmotic gradient (Figure 6D).
In summary, changes in Cl- driving force generated by changing the ionic contributions to cellular charge (by altering the mean charge of impermeant anions) or osmoneutrality (by increasing the contribution of hydrostatic pressure) are due to the alteration of the dynamics of active ion transport mechanisms in the cell. However, these effects are negligible in magnitude and cannot contribute significantly to setting physiologically observed Cl- driving forces. It is worth reiterating that any non-zero Cl- driving force is entirely dependent on the presence of active Cl- cotransport. In our model, in the absence of KCC2, neither the Na+/K+-ATPase nor impermeable anions can shift Cl- out of equilibrium (Figure 6F).
Changes in cation-chloride cotransport activity generate local differences in chloride reversal and driving force, which depend on cytoplasmic diffusion rates
An important functional question is how Cl- driving force might be modified at a local level within a neuron. We considered local persistent changes of Cl- driving force for the case of active transmembrane Cl- fluxes (Figure 7) and impermeant anions (Figure 8) by extending the single-compartment model described above into a multi-compartment model or ‘virtual dendrite.’ This dendrite was 100 μm in length and consisted of 10 compartments, each of 10 μm length and 1 μm diameter. The compartments contained the same mechanisms and default parameterization as the single compartment model described above. Compartmental volume was changed by altering the radius, while holding the length constant. In addition, all ions, except impermeant anions, could move between compartments by electrodiffusion (Figure 7A and Materials and methods).
To explore the local effects of CCC activity, we increased gKCC2 from our default value of 20 µS/cm2 to 600 µS/cm2 in the second distal compartment of the virtual dendrite exclusively. This resulted in a persistent decrease in ECl, concurrent with a modest decrease in Vm, resulting in a permanent increase in Cl- driving force and minimal change in compartment volume (Figure 7B). The spatial precision of this alteration depended strongly on the diffusion constant for Cl-. With a Cl- diffusion constant of 2.03 × 10−7 dm2.s−1, these alterations spread widely through the virtual dendrite. For example, the change in Cl- driving force was 4.8 mV in the furthermost compartment (90 μm apart) as compared to 5.9 mV in the compartment manipulated. When we decreased the Cl- diffusion constant by one order of magnitude, the change in Cl- driving force was 7.3 mV in the compartment in which KCC2 was adjusted, but only 1.8 mV in the furthermost compartment from the site of manipulation (Figure 7C). These findings suggest that local differences in cation-chloride cotransport activity can drive spatially restricted differences in Cl- driving force under conditions of constrained Cl- diffusion; however, under conditions of typical ionic diffusion the effect of Cl- transport by KCC2 is relatively widespread.
Local impermeant anions do not appreciably affect the local driving force for chloride
Following the last result, we considered whether changing impermeant anions in part of a dendrite could create a local area with a different Cl- driving force compared to the rest of the cell. We first added impermeant anions of the default charge (z = −0.85) exclusively to the second-most distal compartment of the virtual dendrite while measuring the Cl- reversal, Vm and Cl- driving force in all compartments. During addition of the impermeant anions, ECl and Vm decreased with an accompanying decrease in Cl- driving force. However, following cessation of impermeant anion influx, all parameters returned to baseline levels, except for the volume of that specific compartment, which showed a modest increase (Figure 8B). This suggests that local addition of impermeant anions of mean charge has no local effect on Cl- homeostasis but can affect the volume of the compartment concerned.
Next, we again added impermeant anions to the second-last compartment of the virtual dendrite, but this time we added impermeant anions with a more negative charge than (z = −1) than the current mean. This resulted in the mean charge of impermeant anions in that compartment becoming more negative (Figure 8C). During the addition of the impermeant anions, ECl and Vm decreased across the dendrite, but with small accompanying shifts in Cl- driving force. Following cessation of local impermeant anion influx, a persistent shift in ECl and Vm was observed specifically in the compartment manipulated. However, this generated a negligible, persistent change in Cl- driving force (<0.01 mV change in driving force for a compartment specific change in z from −0.85 to −0.93), only within that specific compartment of the virtual dendrite. Again, impermeant anion addition resulted in a permanent increase in the volume of the compartment concerned. This finding suggests that local impermeant anions can adjust the Cl- reversal potential locally, but are not well-placed to cause significant, permanent shifts in the driving force for Cl-. Indeed, electrodiffusion may further limit the degree to which local changes in impermeant ion charge can modify driving forces through alterations in active ionic transport: the resulting permanent Cl- driving force changes in the multi-compartment model are many times smaller than the shifts in the single compartment version (as compared to Figure 6B).
Discussion
The driving force for Cl- is a fundamental parameter affecting the excitability of neuronal networks (Raimondo et al., 2017). Recently, impermeant anions, rather than CCCs, have been suggested as the primary determinants of the neuronal driving force for Cl- (Glykys et al., 2014). Here, we have explored the determinants of the Cl- driving force in neurons by deriving theoretical models based on biophysical first principles. We show that the Na+/K+-ATPase, baseline K+, Na+ and Cl- conductances, mean charge of impermeant anions, water permeability and CCCs, likely all play roles in setting neuronal [Cl-]i. However, our findings suggest that while impermeant anions can contribute to setting the [Cl-]i in neurons, they can only affect Cl- driving force by modifying the activity of active transport mechanisms (i.e. the Na+/K+-ATPase). Our modelling and experimental data demonstrate that under physiologically relevant conditions, impermeant anions do not alter the Cl- driving force significantly. In contrast, CCCs are well placed to modulate Cl- driving force and hence inhibitory signaling.
Previous theoretical models, which account for the dynamics of Cl- ions, have been useful in determining how changes to the driving force for Cl- are critical for controlling the effect of synaptic inhibition in the brain (Qian and Sejnowski, 1990; Staley and Proctor, 1999; Doyon et al., 2011; Jedlicka et al., 2011; Lewin et al., 2012; Mohapatra et al., 2016). Whilst these models have included the Na+/K+-ATPase, the interacting dynamics of several ion species, CCCs (Doyon et al., 2011; Krishnan and Bazhenov, 2011), electrodiffusion (Qian and Sejnowski, 1989) and impermeant anions and volume regulation (Dijkstra et al., 2016), none have combined all these mechanisms to explore how their combination determines the local driving force for Cl-. Our theoretical approach is based on the pump-leak formulation (Tosteson and Hoffman, 1960). It suggests that mammalian cells maintain their volume under osmotic stress generated by impermeant anions and the Donnan effect by employing active transport of Na+ and K+ using the Na+/K+-ATPase (Armstrong, 2003; Kay, 2017). A Donnan equilibrium, a true thermodynamic equilibrium requiring no energy to maintain it, is not possible in cells with pliant membranes like neurons (Sperelakis, 2012).
Our model conforms to the pump-leak formulation: abolishing the activity of the Na+/K+-ATPase leads to cell swelling, progressive membrane depolarization and rundown of ionic gradients, including that of Cl-. Therefore, the Na+/K+ ATPase is a fundamental cellular parameter that stabilizes cell volume and determines all ionic gradients including that of Cl- and hence must be considered in any attempt to model ion homeostasis. Interestingly however, we demonstrate that above a certain level of Na+/K+-ATPase activity, even many fold changes in pump rate have minimal effects on volume, ECl and Vm. This might explain recent experimental findings in which periods of Na+/K+ ATPase inhibition using ouabain caused modest changes to cell volume (Glykys et al., 2014). It therefore seems unlikely that neurons adjust the Na+/K+-ATPase as a means for modulating Cl- driving force.
Baseline ion conductances are another important factor that affect Cl- driving force. Our model is consistent with recent experimental results that demonstrate that increased neuronal Na+ conductance (for example by activating NMDA receptors, or preventing closure of voltage-gated Na+ channels), leads to progressive neuronal swelling, membrane depolarization and Cl- accumulation (Rungta et al., 2015) – the primary pathological process in cytotoxic edema (Liang et al., 2007). We also show that tonic neuronal Cl- conductance only affects baseline [Cl-]i and driving force in the presence of CCCs. Without active Cl- flux, which CCCs provide, there is no driving force for passive Cl- flux and hence no mechanism for [Cl-]i changes resulting from selective modification of a Cl- conductance. This is consistent with both classic (Misgeld et al., 1986; Thompson and Gähwiler, 1989) and recent experimental findings (Berglund et al., 2016).
In our model, we find that elevating the activity of KCC2, the most active CCC in mature neurons (Ben-Ari, 2002), increases the driving force for Cl- by shifting the reversal potential for Cl- closer to that of K+. Interestingly, large shifts (~7 mV) in driving force were associated with very minor (1%) changes in volume or membrane potential. As such, modulating KCC2 represents a specific means for manipulating the neuronal Cl- driving force. This is consistent with traditional dogma, recent and previous experimental results (Kaila et al., 2014; Klein et al., 2018) as well as our own experimental validation using furosemide to block the activity of KCC2, which drove significant changes in driving force with little effect on Vm. In further support of this, our meta-analysis of numerous experimental studies showed a strong correlation between change in KCC2 expression and Cl- driving force, but not between KCC2 expression and Vm. There is an ongoing debate as to whether some cotransporters, including CCCs, might also couple water transport to ion transport (Zeuthen, 1994; MacAulay et al., 2002; Gagnon et al., 2004; Charron et al., 2006). Although we did not model the active movement of water by KCC2, this scenario would not alter the central importance of CCCs for setting the Cl- driving force.
Using our multi-compartment model, which incorporated electrodiffusion, we found that local modification of KCC2 activity has a specific local effect on Cl- driving force that is dependent on the characteristics of intracellular Cl- diffusion. Cytoplasmic Cl- diffusion rates had to be reduced substantially before we observed local changes in Cl- driving force driven by KCC2 (Qian and Sejnowski, 1989; Kuner and Augustine, 2000). Whilst differences in KCC2 activity might generate a gradient in Cl- driving force between large subcellular structures (i.e. dendrites versus soma), our modeling results call into question the idea of synapse-specific regulation of Cl- driving force within the same cellular domain (Földy et al., 2010).
Glykys et al. (2014) used Cl- imaging and various experimental manipulations to claim that intracellular and extracellular concentrations of impermeant anions ([X]i and [X]o) set [Cl-]i and the Cl- driving force. From our theoretical analysis, we find that modifying the amount of impermeant anions inside or outside neurons has no persistent effect on [Cl-]i or Cl- driving force, unless we include a mechanism that allows a transmembrane osmotic pressure differential to develop that indirectly affects active transport mechanisms. Even in this case, under transmembrane pressure differentials that do not lyse the membrane (Nichol and Hutter, 1996), Cl- driving force changes are negligible (<1 mV). Recently, it has been suggested that the viscoelastic properties of the cellular cytoskeleton could allow it to take up osmotic shifts created by impermeant anion movement like a sponge (Sachs and Sivaselvan, 2015). This would mean that one would not see as large a volume shifts as predicted by our models. In our model, we have assumed that water can pass through the neuronal membrane to equalize osmotic differences. Although it is thought that some neurons do not express aquaporin channels (Andrew et al., 2007), water can permeate the phospholipid bilayer (Fettiplace and Haydon, 1980). Therefore, whilst differences in neuronal water permeability might affect the time taken to reach steady-state, the steady state values themselves are unlikely to be affected. We conclude that [Cl-]i and the Cl- driving force are not determined by the concentration of impermeant anions.
However, our theoretical findings offer a potential explanation for recent experimental observations. We show that modifying the mean charge of impermeant anions (i.e. z in [Xz]i), rather than their concentration, can affect [Cl-]i and ECl. Relating this to prior experimental observations, Glykys et al. (2014) used SYTO64 staining of nucleic acids and perfusion of weak organic acids in conjunction with Cl- imaging to suggest that [Cl-]i depends upon internal impermeant anions ([X]i). If such a manipulation modifies the mean charge of internal impermeant anions, and not concentration per se, this could account for the observed changes in [Cl-]i. Glykys et al. (2014) did not measure Vm or the Cl- driving force in these experiments. The clear prediction from our model is that any manipulation, which changes the mean charge of impermeant anions would not appreciably affect the Cl- driving force because any impermeant anion driven change on ECl- is matched by an equivalent effect on Vm due to accompanying shifts in cation concentrations. We have provided experimental support for this prediction by showing that whilst EGABA (and ECl) can be shifted by addition of impermeant anions using electroporation of membrane impermeant anionic dextrans, Vm is shifted in a similar direction resulting in an undetectable change in Cl- driving force. Future experiments could further test our model by electroporating positively charged dextrans which would be predicted to depolarize both Vm and ECl, again with minimal effects on Cl- driving force.
Given prior theoretical predictions (Kaila et al., 2014; Voipio et al., 2014; Savtchenko et al., 2017), it is interesting that our model reveals that changing impermeant anions could affect the Cl- driving force at all. We found that the small (<1 mV) impermeant anion-driven changes in Cl- driving force observed in our model were caused by indirect effects on Na+ concentration and hence Na+/K+-ATPase activity. The impermeant anion-driven changes in Cl- driving force are even smaller in the multi-compartment model (<0.1 mV), in which electrodiffusion allows local changes in Na+ to dissipate. When Na+/K+-ATPase activity was decoupled from the transmembrane Na+ gradient, we found that impermeant anions were unable to cause persistent shifts in Cl- driving force as predicted theoretically (Kaila et al., 2014; Voipio et al., 2014; Savtchenko et al., 2017). It is important to note that these small, impermeant anion-Na+/K+-ATPase-driven shifts in Cl- driving force are dependent on the presence of cation-chloride cotransport in the form of KCC2 and would entail changes in energy use by the Na+/K+-ATPase. In other words, active transport mechanisms are again required to drive changes in Cl- homeostasis.
In summary, our theoretical models, which are derived from well-established physical principles, are consistent with our own experimental data and that of others (Glykys et al., 2014; Kaila et al., 2014; Klein et al., 2018), and suggest that impermeant anions alone cannot shift Cl- out of equilibrium across the neuronal membrane. Were neurons to alter impermeant anion concentration or charge, the resting membrane potential would be modified with little effect on the Cl- driving force. Our work confirms the central importance of CCC activity in determining the effects of inhibitory synaptic transmission in the nervous system.
Materials and methods
Single-compartment model
Request a detailed protocolThe single-compartment model consisted of a cylindrical semipermeable membrane separating the extracellular solution from the intracellular milieu with variable volume (Figure 1 ). The extracellular ionic concentrations were assumed constant (Table 1). Permeable ions in the model were K+, Na+ and Cl- with their usual charges, while impermeant anions X were assumed to be a heterogeneous group of impermeant chemical species with a mean intracellular charge z and a mean extracellular charge -1. The default z (-0.85) was chosen on the basis of known resting intracellular ion concentrations (Lodish et al., 2009; Raimondo et al., 2015) and osmolarity (). Bicarbonate ions were not included in our model as a permeant anion as they were assumed to be important for acute depolarizing effects (via GABAARs) rather than the chronic shifts in Cl- driving force, which are the focus of this work (Staley and Proctor, 1999). The model included ionic leak currents for the permeable ions, Na+/K+-ATPase transporters and a CCC, in this case the K+-Cl- cotransporter (KCC2). KCC2 and not NKCC1 is thought to be the most active CCC in mature neurons (Ben-Ari, 2002), therefore, to maintain conceptual simplicity only KCC2 was modeled. Cell volume (w) change was based on osmotic water flux and incorporated a membrane surface area scaling mechanism. An analytical solution to the model at steady state was derived using standard techniques and can be found in Supplementary file 1. The numerical model was initialized assuming conditions close to electroneutrality and an osmotic equilibrium between the intracellular and extracellular compartments. A forward Euler approach was used to update variables at each time step (dt) of 1 ms. Using a smaller dt did not influence the results in Figure 1–5. Code was written in Python 2 and is available on GitHub (Düsterwald and Currin, 2018; copy archived at https://github.com/elifesciences-publications/model-of-neuronal-chloride-homeostasis). The GitHub repository includes a file of Supplementary figures, in which we display transmembrane fluxes of all ions and water for relevant simulations. An example figure displaying ionic flux for all ions is available for Figure 4C in Figure 4—figure supplement 1.
Membrane potential
Request a detailed protocolThe membrane potential Vm was based on the ‘Charge Difference’ approach (Rybak et al., 1997; Fraser and Huang, 2004) as follows:
where F is Faraday’s constant, Cm is the unit membrane capacitance and Am is calculated as the ratio of the surface area (of the cylinder) to cell volume. The term in brackets is the sum of all ionic charges within the cell. This approach has the advantage that the initial voltage can be calculated without needing to assume a steady state as is required for by the Goldman-Hodgkin-Katz (GHK) equation.
Permeable ion concentrations
Request a detailed protocolIntracellular concentrations of the permeable ions Na+, K+ and Cl- were updated individually by summing trans-membrane fluxes. Leak currents were calculated using Ohm's Law, . In addition, Na+ and K+ were transported actively by the Na+/K+-ATPase, with pump rate Jp, which was approximated by a cubic function dependent on the transmembrane sodium gradient, following (Keener and Sneyd, 1998):
where P is the pump rate constant. Because it is a function of the sodium gradient, Jp decreases as [Na+]i depletes. This formulation has been shown to be similar to more accurate kinetic models reliant on both the Na+ gradient and ATP concentration (Keener and Sneyd, 1998). To switch the ATPase pump on or off (Figure 1C), P was decreased/increased exponentially over 10–20 min, consistent with previous reports of the dynamics of inhibition of the ATPase by ouabain and in turn the inhibition of ouabain’s effects by potassium canrenoate (Baker and Willis, 1972; Yeh and Lazzara, 1973). The ATPase pumps 2 K+ ions into the cell for every 3 Na+ ions out and these constants must be multiplied by Jp for each ion, respectively. K+ and Cl- were also modified by flux through the type 2 K-Cl cotransporter (KCC2), which has a stoichiometry of 1:1 and transports both ions in the same direction. Flux though KCC2, JKCC2 (Doyon et al., 2016), was modeled as follows:
where gKCC2 is a fixed conductance and EK and ECl are the Nernst potentials for K+ and Cl- respectively. JKCC2 is 0 when EK=ECl. The rate of change of the intracellular concentration of the three permeant ions was given by the following equations, with the Nernst potentials for each ion given by , w indicating the cell volume, and as described in Equation 7 (Figure 1–5,6A,B,7 and 8 or 9 or Equation 9 (Figure 6C–E):
Volume
Request a detailed protocolIn most calculations, because the osmotic flux of water is expected to be faster than ion fluxes, the volume of the cell (w) was adjusted to reduce the difference between (intracellular osmolarity) and (extracellular osmolarity) at each time step by explicitly modelling water flux, where vw is the partial molar volume of water, pw the osmotic permeability of a biological membrane and SA the surface area (Hernández and Cristina, 1998):
In some calculations (Figure 6C–E) in which we allowed transmembrane differences in osmolarity to develop, we assumed that at rest the cylindrical cell had a radius of ra and zero pressure across the membrane, and that the tension (T) in the membrane followed Hooke’s law such that the tension was proportional to the difference between the dynamic circumference of the cell and that of the resting state. From Laplace’s law the hydrostatic pressure in the cell was given by:
where km is the spring constant of the membrane (Sachs and Sivaselvan, 2015). Equation 7 was thus reformulated:
In order to simulate extreme conditions of constrained volume, a larger km was employed than is realistic (Dai et al., 1998). Intracellular ion concentrations were updated again after volume change at each time step. Volume changes were manifested in the cylindrical compartment as change in the radius. In Figures 1–6, the cell was initialised with diameter 10 μm and length 25 μm.
Anion flux
Request a detailed protocolImpermeant anions were manipulated in the compartment in Figures 4–6 and Figure 8 through several mechanisms. Anions could be added to the compartment at a constant rate and have either the default intracellular X charge z = −0.85 (Figures 4C, 6C and 8B), or a different charge (Figures 5C, 6A and 8C). In these cases, the number of moles of X in the compartment was increased. Alternatively, the charge of a species of intracellular X was slowly changed imitating a charge-carrying trans-membrane reaction (Figures 5A and 6A). In this case, the number of moles of intracellular X did not change and it was assumed charge imbalance was compensated by the extracellular milieu. Finally, extracellular X- was changed in Figure 4D by removing as much Cl- as X- was added, thus maintaining osmolarity and electroneutrality in the extracellular space.
Multi-compartment model
Request a detailed protocolThe single-compartment dendrite model was incorporated in a multi-compartment model by allowing electrodiffusion to occur between individual compartments operating as described above. Compartments were initialised with a radius of 0.5 μm and length of 10 μm. Compartments were linked linearly without branching; 10 connected compartments in total were used. The time step dt was decreased to 10−3 ms for simulations in multiple compartments. Code was written in Python 3 and is available on GitHub (Düsterwald and Currin, 2018).
Electrodiffusion. The Nernst-Planck equation (NPE) was used to model one-dimensional electrodiffusion, based on Qian and Sejnowski (1989). The NPE incorporates fluxes because of diffusion and drift (i.e. the movement of ions driven by an electric field). It has been shown to be more accurate than using Jdiffusion alone in small structures like dendrites (Qian and Sejnowski, 1989). The NPE for J the flux density of ion C is calculated as:
where D is the diffusion constant of ion C (Table 1), z is its charge, [C] is its concentration and x is the distance along the longitudinal axis over which electrodiffusion occurs. The NPE was implemented between compartments i and i + 1, assuming the i→i + 1 direction was positive, using a forward Euler approach. The midpoints of the compartments were used to calculate dx, i.e., and the concentrations of C in each compartment were averaged to obtain Jdrift, ensuring that Ji→i+1 = Ji+1→i, where the fluxes had units of mol/(s.dm2):
The flux was multiplied by the surface area between compartments and then divided by compartment volume to determine the flux in terms of molar concentration (M/s), i.e. and finally implemented numerically with a forward Euler approach. The implementation mirrored that in Qian and Sejnowski (1989) for non-branching dendrites, but was adjusted for compartments whose volumes can change:
Systematic review
Request a detailed protocolA literature search was performed to identify experimental studies that aimed to correlate a change in KCC2 expression with changes in [Cl-]i. The MEDLINE database was used and accessed via the PubMed online platform. Search terms included ‘chloride’, ‘Cl’, ‘intracellular’, ‘KCC2’, ‘cotransporter’, ‘neuronal’, ‘GABA’ using appropriate Boolean operators. All 26 studies that demonstrated changes in KCC2 expression and EGABA were considered for the meta-analysis. As there is a well-described differential expression of KCC2 and NKCC1 at different stages of development, with KCC2 expression increasing and NKCC1 expression decreasing across development, only studies that used tissue older than postnatal day seven were included (eight included data from younger animals). Other exclusion criteria included: reporting a significant change in NKCC1 (five studies); use of non-rodent tissue (two studies); no quantification of the change in KCC2 (two studies). Nine experiments from eight studies met all criteria and were included (Coull et al., 2003; Lagostena et al., 2010; Lee et al., 2011; Ferrini et al., 2013; Campbell et al., 2015; MacKenzie and Maguire, 2015; Mahadevan et al., 2015; Tang et al., 2015). However, one study did not report the change in Vm and hence was excluded from the figure (Mahadevan et al., 2015). Data used in regression can be seen in Supplementary file 2 (Table S2-1) and includes the analysis for regression against change in [Cl-]i. To accommodate varied experimental preparations and techniques influencing data quality and biases, a 34-point scoring system was designed to weight the studies (Table S2-2). A weighted least squares regression model was then used to correlate the percentage (%) change in KCC2 expression versus change in Cl-d driving force.
Slice preparation and electrophysiology
Request a detailed protocolFor all experiments, organotypic slices were prepared from rodent brain tissue. Wistar rats were used for the experiments testing the effects of CCC blockade. However, to allow for optogenetic manipulation, a crossed mouse strain on a C57BL/6 background was used. This mouse strain was a cross between mice expressing cre-recombinase in glutamic acid decarboxylase 2 (GAD2) positive interneurons (GAD2-IRES-cre, Jax Lab strain 010802) and mice with a loxP-flanked STOP cassette preventing transcription of the red-fluorescent protein tdTomato (a cre-reporter strain, Ai4, Jax Lab strain 007914). This created the GAD2-cre-tdTomato strain which resulted in cre-recombinase and tdTomato expression in all GABAergic interneurons (Taniguchi et al., 2011). The use of all animals was approved by either the University of Oxford (rat) or the University of Cape Town (mouse) animal ethics committees.
Organotypic brain slices were prepared using 7 day old Wistar rats (CCC experiment) or crossed GAD2-IRES-cre (RRID:IMSR_JAX:010802) and Ai14 tdtomato reporter (RRID:IMSR_JAX:007914) mice (impermeant anion experiment) and followed the protocol originally described by Stoppini et al., 1991. Briefly, brains were extracted and swiftly placed in cold (4°C) dissection media consisting of Gey’s Balanced Salt Solution (GBSS #G9779, Sigma-Aldrich, USA) supplemented with D-glucose (#G5767, Sigma-Aldrich, USA). The hemispheres were separated and individual hippocampi were removed and immediately cut into 350 μm slices using a Mcllwain tissue chopper (Mickle, UK). Cold dissection media was used to rinse the slices before placing them onto Millicell-CM membranes (#Z354988, Sigma-Aldrich, USA). Slices were maintained in culture medium consisting of 25% (vol/vol) Earle’s balanced salt solution (#E2888, Sigma-Aldrich, USA); 49% (vol/vol) minimum essential medium (#M2279, Sigma-Aldrich, USA); 25% (vol/vol) heat-inactivated horse serum (#H1138, Sigma-Aldrich, USA); 1% (vol/vol) B27 (#17504044, Invitrogen, Life Technologies, USA) and 6.2 g/l D-glucose. Slices were incubated in a 5% carbon dioxide (CO2) humidified incubator at between 35–37°C. Recordings were made after 6–14 days in culture. Previous studies have shown that after 7 days in culture (equivalent to postnatal day 14) GABAergic signaling and EGABA has sufficiently developed to a level resembling mature nervous tissue (Streit et al., 1989; Wright et al., 2017). For the impermeant anion experiment, the mouse organotypic brain slices were injected 1 day after culture with adeno-associated vector serotype 1 (AAV1) containing a double-floxed sequence for channelrhodopsin (ChR2) linked to a yellow fluorescent protein (YFP) tag driven by the elongation factor one promoter (UNC Vector Core, USA). The vector was diffusely injected into slices using a custom-built Openspritzer pressurized ejection device (Forman et al., 2017). Slices were left for 6 days in culture to allow for robust expression of ChR2-YFP in GAD2+ interneurons.
For recordings, slices were transferred to a submerged chamber which was perfused with artificial cerebro-spinal fluid (aCSF) bubbled with carbogen gas (95% oxgen:5% carbon dioxide). aCSF was composed of (in mM): NaCl (120); KCl (3); MgCl2 (2); CaCl2 (2); NaH2PO4 (1.2); NaHCO3 (23); D-Glucose (11) with pH adjusted to be between 7.35–7.40 using 0.1 mM NaOH. Neurons were visualized using a 20x or 60x water-immersion objective (Olympus, Japan) on a BX51WI upright microscope (Olympus, Japan). Widefield images were obtained using a Mightex CCE-B013-U CCD camera. For the optogenetic experiments, we used epifluorescence microscopy to select slices that exhibited strong ChR2-YFP expression throughout the hippocampal region. Micropipettes were prepared from borosilicate glass capillaries with an outer diameter of 1.20 mm and inner diameter of 0.69 mm (Warner Instruments, USA), using a horizontal puller (Sutter, USA). Recordings were made using Axopatch 1D and Axopatch 200B amplifiers and acquired using pClamp9 (Molecular Devices) or WinWCP software (University of Strathclyde, UK).
For the CCC blockade experiment, gramicidin perforated patch recordings (Kyrozis and Reichling, 1995) were performed using glass pipettes containing (in mM): 135 KCl, 4 Na2ATP, 0.3 Na3GTP, 2 MgCl2, 10 HEPES and 80 μg/ml gramicidin (Calbiochem; pH 7.38; osmolarity, 290 mosmol/l). After obtaining a cell-attached patch, the gramicidin perforation process was evaluated by continuously monitoring the decrease in access resistance. Recordings were started once the access resistance had stabilized between 20–80 MΩ, which usually occurred 20–40 min following gigaseal formation. Rupture of the gramicidin patch, referred to as a break through, induces a large influx of the high Cl- internal solution into the cell. This causes a significant and permanent increase in the EGABA at which time recordings would be discarded. GABAAR activation was achieved by pressure application of muscimol (10 μM, Tocris, UK), a selective GABAAR agonist, via a picospritzer. To calculate EGABA, GABAAR currents were elicited at different command voltages. These were a series of 10 mV steps above and below a holding voltage of −60 mV. Reported membrane potentials were corrected for the voltage drop across the series resistance for each neuron. Holding current (reflecting membrane current) and total current (reflecting membrane current plus the GABAAR-evoked current) were plotted against the corrected holding potential to generate a current-voltage (I-V) curve. Using this graph, the EGABA was defined as the potential where the total current equals the holding current. Some of the data used for calculating EGABA values in Figure 3E was also used in a previous study (Wright et al., 2017). Vm was defined as the x-intercept of the holding current, and the driving force as the difference between the two. To block CCC function, the KCC2 blocker furosemide (1 mM, Sigma, USA) was applied.
For the impermeant anion experiments, whole-cell recordings were utilized as electroporation consistently ruptured the gramicidin patches. Pipettes were filled with internal solution composed of (in mM): K-Gluconate (120); KCl (10); Na2ATP (4); NaGTP (0.3); Na2-phosphocreatinine (10) and HEPES (10). To test the effect of introducing impermeant anions, we electroporated the anionic 10 000 MW Dextran Alexa-Flour 488 (Thermo Fisher, USA) using a separate pipette positioned near the soma of the patched cell. This molecule is a hydrophilic polysaccharide, which is both membrane impermeant and highly negatively charged. Pipettes were filled with a 5% dextran solution in phosphate buffered saline and voltage pulses (5–10 of 20 ms duration, 0.5–1 V) were applied using a stimulus isolator. Successful electroporation of the anionic dextran was confirmed visually by observing the cell fill with the fluorescent dye. Electroporation also resulted in immediate membrane depolarization, which recovered over a period of 1–5 mins. EGABA, Vm and driving force were calculated during voltage steps in voltage clamp mode (as above) with GABAARs activated via endogenous synaptic release of GABA using photo-activation of GAD2+ interneurons expressing ChR2-YFP with 100 ms pulses of blue light using a 470 nm LED (Thorlabs), in the presence of 5 μM CGP-35348 (Torcis Bioscience, UK) to block GABABR activation. EGABA, Vm and driving force were calculated before and atleast five mins following electroporation to allow for stabilization of Vm, EGABA and driving force. The image of the dextran filled cell was acquired using a confocal microscope (LSM510 Meta NLO, Car Zeiss, Jena, Germany). Analysis was performed using custom written scripts in the MATLAB environment (MathWorks). Results are presented as mean ± SEM.
Data availability
Code data is available on GitHub (https://github.com/kiradust/model-of-neuronal-chloride-homeostasis; copy archived at https://github.com/elifesciences-publications/model-of-neuronal-chloride-homeostasis). Experimental data in the form of data spreadsheets has been included, and full experimental data is available on Dryad.
-
Dryad Digital RepositoryData from: Biophysical models reveal the relative importance of transporter proteins and impermeant anions in chloride homeostasis.https://doi.org/10.5061/dryad.kj1f3v4
References
-
Excitatory actions of gaba during development: the nature of the nurtureNature Reviews Neuroscience 3:728–739.https://doi.org/10.1038/nrn920
-
Optogenetic visualization of presynaptic tonic inhibition of cerebellar parallel fibersJournal of Neuroscience 36:5709–5723.https://doi.org/10.1523/JNEUROSCI.4366-15.2016
-
The composition of animal cells: solutes contributing to osmotic pressure and charge balanceComparative Biochemistry and Physiology Part B: Comparative Biochemistry 76:663–671.https://doi.org/10.1016/0305-0491(83)90375-9
-
Membrane tension in swelling and shrinking molluscan neuronsThe Journal of Neuroscience 18:6681–6692.https://doi.org/10.1523/JNEUROSCI.18-17-06681.1998
-
In situ calibration and [H+] sensitivity of the fluorescent Na+ indicator SBFIAmerican Journal of Physiology-Cell Physiology 280:C1623–C1633.https://doi.org/10.1152/ajpcell.2001.280.6.C1623
-
Ionic mechanism of ouabain-induced swelling of leech Retzius neuronsPflügers Archiv - European Journal of Physiology 452:25–35.https://doi.org/10.1007/s00424-005-0009-6
-
A biophysical model for cytotoxic cell swellingThe Journal of Neuroscience 36:11881–11890.https://doi.org/10.1523/JNEUROSCI.1934-16.2016
-
The cellular, molecular and ionic basis of GABA(A) receptor signallingProgress in brain research 160:59–87.https://doi.org/10.1016/S0079-6123(06)60005-8
-
Water permeability of lipid membranesPhysiological Reviews 60:510–550.https://doi.org/10.1152/physrev.1980.60.2.510
-
Regulation of fast-spiking basket cell synapses by the chloride channel ClC-2Nature Neuroscience 13:1047–1049.https://doi.org/10.1038/nn.2609
-
Properties of the na+/K+ pump current in small neurons from adult rat dorsal root gangliaBritish Journal of Pharmacology 138:1517–1527.https://doi.org/10.1038/sj.bjp.0705170
-
Modeling cell volume regulation in Nonexcitable cells: the roles of the na+ pump and of cotransport systemsAmerican Journal of Physiology-Cell Physiology 275:C1067–C1080.https://doi.org/10.1152/ajpcell.1998.275.4.C1067
-
A. fundamental studies. On the theory of the Donnan membrane equilibriumDiscussions of the Faraday Society 21:31–45.https://doi.org/10.1039/df9562100031
-
Perturbed chloride homeostasis and GABAergic signaling in human temporal lobe epilepsyJournal of Neuroscience 27:9866–9873.https://doi.org/10.1523/JNEUROSCI.2761-07.2007
-
Expression of GABA signaling molecules KCC2, NKCC1, and GAD1 in cortical development and schizophreniaJournal of Neuroscience 31:11088–11095.https://doi.org/10.1523/JNEUROSCI.1234-11.2011
-
Effect of anoxia on intracellular and extracellular potassium activity in hypoglossal neurons in vitroJournal of Neurophysiology 66:103–111.https://doi.org/10.1152/jn.1991.66.1.103
-
Cation-chloride cotransporters in neuronal development, plasticity and diseaseNature Reviews Neuroscience 15:637–654.https://doi.org/10.1038/nrn3819
-
How cells can control their size by pumping ionsFrontiers in Cell and Developmental Biology 5:1–14.https://doi.org/10.3389/fcell.2017.00041
-
Tenuous inhibitory GABAergic signaling in the reticular thalamusThe Journal of Neuroscience 38:1232–1248.https://doi.org/10.1523/JNEUROSCI.1345-17.2017
-
Ionic dynamics mediate spontaneous termination of seizures and postictal depression stateJournal of Neuroscience 31:8870–8882.https://doi.org/10.1523/JNEUROSCI.6200-10.2011
-
Perforated-patch recording with gramicidin avoids artifactual changes in intracellular chloride concentrationJournal of Neuroscience Methods 57:27–35.https://doi.org/10.1016/0165-0270(94)00116-X
-
Cytotoxic edema: mechanisms of pathological cell swellingNeurosurgical Focus 22:1–9.https://doi.org/10.3171/foc.2007.22.5.3
-
Biology Molecular CellIntracellular Ion Environment and Membrane Electric Potential, Biology Molecular Cell, Fourth edition, New York, WH Freeman and Company.
-
Neto2-null mice have impaired GABAergic inhibition and are susceptible to seizuresFrontiers in Cellular Neuroscience 9:.https://doi.org/10.3389/fncel.2015.00368
-
Tensile strength and dilatational elasticity of giant sarcolemmal vesicles shed from rabbit muscleThe Journal of Physiology 493:187–198.https://doi.org/10.1113/jphysiol.1996.sp021374
-
Chloride regulation in the pain pathwayBrain Research Reviews 60:149–170.https://doi.org/10.1016/j.brainresrev.2008.12.015
-
Optogenetic silencing strategies differ in their effects on inhibitory synaptic transmissionNature Neuroscience 15:1102–1104.https://doi.org/10.1038/nn.3143
-
Ion dynamics during seizuresFrontiers in Cellular Neuroscience 9:1–14.https://doi.org/10.3389/fncel.2015.00419
-
Neuronal chloride and excitability - the big impact of small changesCurrent Opinion in Neurobiology 43:35–42.https://doi.org/10.1016/j.conb.2016.11.012
-
Mechanism of activity-dependent downregulation of the neuron-specific K-Cl cotransporter KCC2Journal of Neuroscience 24:4683–4691.https://doi.org/10.1523/JNEUROSCI.5265-03.2004
-
Cell volume control in three dimensions: Water movement without solute movementThe Journal of General Physiology 145:373–380.https://doi.org/10.1085/jgp.201411297
-
Electrodiffusion phenomena in neuroscience: a neglected companionNature reviews. Neuroscience 18:598–612.https://doi.org/10.1038/nrn.2017.101
-
Cell Physiology Source Book147–151, Chapter 10 – Gibbs–Donnan Equilibrium Potentials, Cell Physiology Source Book, Academic Press, 10.1016/B978-0-12-387738-3.00010-X.
-
Modulation of mammalian dendritic GABA(A) receptor function by the kinetics of cl- and HCO3- transportThe Journal of Physiology 519:693.
-
A simple method for organotypic cultures of nervous tissueJournal of Neuroscience Methods 37:173–182.https://doi.org/10.1016/0165-0270(91)90128-M
-
Role of the potassium chloride cotransporter isoform 2-mediated spinal chloride homeostasis in a rat model of visceral hypersensitivityAmerican Journal of Physiology-Gastrointestinal and Liver Physiology 308:G767–G778.https://doi.org/10.1152/ajpgi.00313.2014
-
Regulation of cell volume by active cation transport in high and low potassium sheep red cellsThe Journal of General Physiology 44:169–194.https://doi.org/10.1085/jgp.44.1.169
-
Postnatal changes in somatic gamma-aminobutyric acid signalling in the rat hippocampusEuropean Journal of Neuroscience 27:2515–2528.https://doi.org/10.1111/j.1460-9568.2008.06234.x
-
Neuronal chloride regulation via KCC2 is modulated through a GABAb>Bb> receptor protein complexThe Journal of Neuroscience 37:5447–5462.https://doi.org/10.1523/JNEUROSCI.2164-16.2017
Article and author information
Author details
Funding
Mandela Rhodes Foundation
- Kira Michaela Düsterwald
- Richard Joseph Burman
Deutscher Akademischer Austauschdienst (DAAD-NRF)
- Christopher Brian Currin
National Research Foundation of South Africa (DAAD-NRF)
- Christopher Brian Currin
H2020 European Research Council (ERC Grant Agreement 617670)
- Colin J Akerman
Newton Fund
- Joseph Valentino Raimondo
École Polytechnique Fédérale de Lausanne (Blue Brain Project)
- Joseph Valentino Raimondo
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Juha Voipio for critically reading the manuscript and providing helpful comments and suggestions.
Ethics
Animal experimentation: No experiments were performed on animals prior to humane killing using cervical dislocation at P7, every effort was made to minimise stress prior to killing. The use of rats in the UK (furosemide experiment) were in accordance with regulations from the United Kingdom Home Office Animals (Scientific Procedures) Act. Use of mice (electroporation experiment) at the University of Cape Town was in accordance with South African national guidelines (South African National Standard: The care and use of animals for scientific purposes, 2008) and was approved by the University of Cape Town Animal Ethics Committee protocol no 014/035.
Copyright
© 2018, Düsterwald et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,377
- views
-
- 391
- downloads
-
- 33
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.