Figures and data in Biophysical models reveal the relative importance of transporter proteins and impermeant anions in chloride homeostasis

Figures
Tables
Additional files

8 figures, 1 table and 3 additional files

Figures

Figure 1 with 2 supplements

Download asset Open asset

A biophysical model of ion dynamics based on the pump-leak mechanism demonstrates the importance of the sodium-potassium ATPase for setting transmembrane ion gradients including chloride.

(A) A single-cell compartment was modeled as a cylinder with volume changes equivalent to changes in cylindrical radius. Dynamics of membrane permeable potassium (purple, K⁺), sodium (pink, Na⁺) and chloride (green, Cl^-) ions were included. Impermeant anions (orange, X^z) had a mean intracellular charge z of -0.85. The KCC2 transporter moved Cl^- and K⁺ in equal parts according to the transmembrane gradient for the two ions. The Na⁺/ K⁺ ATPase transported 3 Na⁺ ions out for 2 K⁺ ions moved into the cell. (B) Regardless of intracellular starting concentrations of the permeable ions, the model converged to identical steady state values for all parameters without needing to include any means for ‘sensing’ ion concentration, volume or voltage. We show the result for Cl^- as a time series of [Cl^-]_i (top panel) and volume (bottom panel). (C) The ATPase plays a key role in maintaining steady state ion concentration, membrane voltage (V_m) and volume. Switching off the ATPase results in a continuous increase in cell volume (bottom panel), membrane depolarization (middle panel) and ion concentration dysregulation (top panel with colours per ion as in ‘A’). All cellular parameters recovered when the ATPase was reactivated. (D) The model’s analytic solution showed exact correspondence with steady state values generated by numerical, time series runs (dots) for varying ATPase pump rates. Steady state values for the concentrations of the ions with colours as in ‘A’ (top panel); V_m (middle panel) and volume (bottom panel). The dashed line indicates the default ATPase pump rate used in all simulations unless specified otherwise. The result in ‘B’ was replicated with Hamada et al.’s experimentally validated model of the ATPase in Figure 1—figure supplement 1B, and the models’ respective pump fluxes are compared in Figure 1—figure supplement 1A. Note that although sufficient Na⁺/K⁺-ATPase activity is critical for steady state ionic gradients these variables are relatively stable near the default pump rate (dashed line), Figure 1—figure supplement 2.

https://doi.org/10.7554/eLife.39575.003

Figure 1—figure supplement 1

Download asset Open asset

Different models of the sodium-potassium ATPase have distinct kinetic properties but produce similar qualitative effects.

(A) Comparison of pump rate of different ATPase models against input sodium concentration. The model by Keener and Sneyd (1998) (solid line) is explained in the main paper methods. The Hamada et al. (2003) method (dashed line) is based on experimental evidence fitted by Hill equations (Hamada et al., 2003). The simplified equation follows: $J_{p} = Q_{s c a l e} (\frac{1.62}{1 + {(\frac{6.7 m M}{{[{N a}^{+}]}_{i}})}^{3}} + \frac{1.0}{1 + {(\frac{67.6 m M}{{[{N a}^{+}]}_{i}})}^{3}})$ , where $Q_{s c a l e}$ was set such that the same steady state concentration of Na⁺ ions (and indeed of all ions) was achieved using either model, with units as for P the pump rate constant. (B) Figure 1B was run using the Hamada et al. (2003) model. As for the standard ATPase model, regardless of intracellular starting concentrations of the permeable ions, the model converged to identical steady state values for all parameters without needing to include any means for ‘sensing’ ion concentration, volume or voltage. We show the result for Cl^- as a time series of [Cl^-]_i (top panel) and volume (bottom panel).

https://doi.org/10.7554/eLife.39575.004

Figure 1—figure supplement 2

Download asset Open asset

Changes in Na⁺/K⁺-ATPase activity near the model’s default pump rate produce minimal changes in E_Cl and V_m.

Steady-state E_K (purple), E_Cl (green) and V_m (black, top panel) as well as volume (bottom panel) were calculated for different Na⁺/K⁺-ATPase pump rate constants using the analytical solution. Note that these variables are relatively stable near the default Na⁺/K⁺-ATPase pump rate (dashed line).

https://doi.org/10.7554/eLife.39575.005

Figure 2 with 2 supplements

Download asset Open asset

Membrane conductances affect steady-state intracellular chloride concentration only in the presence of cation-chloride cotransport.

Steady state values for different ionic conductance were calculated using the model’s analytic solution. (A) Steady state E_Cl (green), E_K (purple), V_m (black) and volume w (bottom panel) were calculated at different K⁺conductances (g_K). Increasing g_K resulted in a convergence of steady state E_Cl and V_m. (B) Increasing Na⁺conductance (g_Na) resulted in a progressive increase in steady state E_Cl, V_m and volume with a negligible increase in E_K. Figure 2—figure supplement 2 demonstrates how increased K⁺flux through the Na⁺/K⁺ ATPase matches passive K⁺efflux to maintain E_K. (C) In the presence of active cation-chloride cotransport (g_KCC2 = 20 µS/cm²), increasing Cl^- conductance shifted steady state E_Cl from E_K toward V_m. (D) In the absence of KCC2 activity, g_Cl had no effect on steady state parameters. E_Cl equals V_m in all instances. Dashed lines indicate the default values for g_K, g_Na and g_Cl. In Figure 2—figure supplement 1 the classic dependence of the membrane potential on log([K⁺]_o) is shown.

https://doi.org/10.7554/eLife.39575.006

Figure 2—figure supplement 1

Download asset Open asset

The classic dependence of membrane potential (mV) on the logarithm of extracellular potassium is reproduced using our model.

The analytical solution was used to explore the dependence of membrane potential on extracellular potassium concentration. In agreement with experimental results, the slope is non-linear at low [K⁺]_o, which is not seen when using a pure Nernst-potential modelling approach.

https://doi.org/10.7554/eLife.39575.007

Figure 2—figure supplement 2

Download asset Open asset

K⁺flux is balanced by active and passive fluxes, dependent on Na⁺steady states.

Varying g_Na in the analytical solution caused a reduction in E_Na (top panel, pink) and non-zero increases in E_K (middle panel, purple). The E_K changes with increasing g_Na are small because the increased K⁺ efflux through passive leak membrane channels accompanying membrane depolarisation (bottom panel, solid line) are compensated by increased ATPase K⁺ influx driven by the Na⁺ concentration changes (dashed line). KCC2 plays a smaller role (dotted line).

https://doi.org/10.7554/eLife.39575.008

Figure 3 with 2 supplements

Download asset Open asset

Cation-chloride cotransport sets the chloride reversal and driving force for transmembrane chloride flux.

(A) Increasing KCC2 activity in our model by increasing g_KCC2 from 20 µS/cm² to 370 µS /cm² resulted in a persistent decrease in E_Cl (green), a minimal decrease in E_K (purple), V_m (black) and volume w (bottom panel). This resulted in a permanent increase in the DF for Cl^- via a change in E_Cl from −83.9 mV to −93.2 mV (red). (B) The steady state values for E_Cl, E_K, V_m (top panel), volume (middle panel) and Cl^-driving force (DF, bottom panel, red) at different KCC2 conductances. Increasing KCC2 activity resulted in a decrease in steady state E_Cl and an increase in DF. E_K represented a lower bound on E_Cl at high KCC2 conductances. Similar results were noted for other kinetic models of KCC2 (Figure 3—figure supplements 1 and 2). (C) Schematic showing experimental setup. Gramicidin perforated patch-clamp recordings were performed on CA3 pyramidal cells from rat hippocampal organotypic brain slices. (D) Insets depict GABA_AR currents elicited by somatic application (20 ms) of muscimol (10 μM) at different voltages. Calibration: 500 ms, 500 pA. Holding current (reflecting membrane current) and total current (reflecting membrane current plus the muscimol-evoked current) were measured at the points indicated by the vertical grey and back lines, respectively. Current-voltage (IV) plots were drawn to calculate changes in V_m, E_GABA and DF before (left) and after (right) furosemide (pink) was applied. Voltages were corrected for series resistance error. (E) Top, population data showing significant changes in E_GABA and DF but not V_m five minutes after furosemide application. Bottom, changes in DF over time show a significant decrease from baseline once furosemide was introduced. (F) Meta-analysis of experimental studies demonstrates a correlation between KCC2 activity (% change) and Cl^- DF (mV, top plot, red) but not membrane potential (mV, bottom plot, grey), confirming the role of KCC2 for establishing the neuronal Cl^- gradient in adult tissue. The data and scoring system used to generate the regression can be found in Supplementary file 2 (Tables S2-1 and S2-2 ‘ns’, non-significant; *p<0.05; **p<0.01. The data for ‘C’, ‘D’ and ‘E’ is provided in Figure 3—source data 1.

https://doi.org/10.7554/eLife.39575.009

Figure 3—source data 1 Source data for Figure 3C–E. Resting membrane potential, E_GABA and driving force at time points as indicated relative to addition of 1 mM furosemide (at time 0 min), for each experiment.: https://doi.org/10.7554/eLife.39575.012
Download elife-39575-fig3-data1-v2.xlsx

Figure 3—figure supplement 1

Download asset Open asset

Models of KCC2 have different kinetic properties.

The pump rate velocity for different models of KCC2 were compared for differing intracellular chloride concentration. The equation employed in our model was based on Doyon et al. (2016) (solid line) and is described in the main text. Lewin et al. (2012) (dotted line) used the following equation in their computational model: $J_{K C C 2} = J_{s c a l e} ∙ 0.3 ∙ \frac{{[K^{+}]}_{o} ∙ {[{C l}^{-}]}_{o} - {[K^{+}]}_{i} ∙ {[{C l}^{-}]}_{i}}{k_{K} ∙ k_{C l} ((1 + \frac{{[K^{+}]}_{o} ∙ {[{C l}^{-}]}_{o}}{k_{K} ∙ k_{C l}}) (1 + \frac{{[K^{+}]}_{i}}{k_{K}}) (1 + \frac{{[{C l}^{-}]}_{i}}{k_{C l}}) + (1 + \frac{{[K^{+}]}_{i} ∙ {[{C l}^{-}]}_{i}}{k_{K} ∙ k_{C l}}) (1 + \frac{{[K^{+}]}_{o}}{k_{K}}) (1 + \frac{{[{C l}^{-}]}_{o}}{k_{C l}}))}$ , where $k_{K} = 9.0 m M a n d k_{C l} = 6.0 m M$ . The model was based on a model of the renal distal tubule (Lewin et al., 2012). Here, it had a diminished absolute velocity compared to our model. The models of Raimondo et al. (2012) and Fraser and Huang (2004) (dashed and dash-dotted lines respectively) had similar velocities across our range of concentrations. Raimondo et al.’s model is based on a model for KCC2 derived using Michaelis-Menton kinetics and experiments, and incorporates new experimental values: $J_{K C C 2} = J_{s c a l e} ∙ V_{m a x} ∙ \frac{{[{C l}^{-}]}_{i}}{15 m M + {[{C l}^{-}]}_{i}} (\frac{\ln \frac{{[K^{+}]}_{o} ∙ {[{C l}^{-}]}_{o}}{{[K^{+}]}_{i} ∙ {[{C l}^{-}]}_{i}}}{\ln \frac{{[K^{+}]}_{o} ∙ {[{C l}^{-}]}_{o}}{{[K^{+}]}_{i} ∙ 56 m M}})$ , where $V_{m a x} = 5 m M ∙ s^{- 1}$ (Raimondo et al., 2012). Fraser and Huang based their model on previous experimental evidence (Fraser and Huang, 2004). Their equation for KCC2 follows: $J_{K C C 2} = J_{s c a l e} ∙ P_{K} ∙ P_{K C C 2} ({[K^{+}]}_{o} ∙ {[{C l}^{-}]}_{o} - {[K^{+}]}_{i} ∙ {[{C l}^{-}]}_{i})$ , where $P_{K}$ and $P_{K C C 2}$ are the permeability constants for potassium and the pump respectively, the latter of which is tuned as needed. $J_{s c a l e}$ for each model was set such that the same steady state concentrations were achieved using any model with the other parameters kept constant.

https://doi.org/10.7554/eLife.39575.010

Figure 3—figure supplement 2

Download asset Open asset

KCC2 and ATPase models with different kinetics have similar properties when KCC2 is upregulated.

(A) Increasing KCC2 activity in our model using the Raimondo et al. (2012) model by increasing V_max/J_Scale resulted in a persistent decrease in E_Cl (green), a minimal decrease in E_K (purple), V_m (black) and volume w (bottom panel). This resulted in a permanent increase in the driving force for Cl^-. The change in KCC2 strength was matched so that the final shift is similar to the usual shift in driving force gained using our standard KCC2 pump formulation. (B) When the ATPase model of Hamada et al. (2003) was used in the experimental set up in (A) rather than the usual ATPase model, the resulting driving force differed by only 0.1 mV from the usual value.

https://doi.org/10.7554/eLife.39575.011

Figure 4 with 1 supplement

Download asset Open asset

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

(A) Initiating the model with different starting concentrations of intracellular impermeant anions ([X]_i) with the same mean charge z = −0.85 (orange, top panel), led to compensatory volume changes (bottom panel) which resulted in identical steady state concentrations. (B) Steady state E_Cl (green), E_K (purple) and V_m (black) were identical regardless of initial [X]_i. Final volume, however, was a linear function of initial [X]_i (bottom panel). (C) Addition of impermeant anions of the mean charge (z = −0.85) caused transients shifts in E_Cl (green, top panel), E_K (purple), V_m (black) as well as [X]_i (orange, bottom panel) for the duration of the influx, and sustained increases in volume (black, middle panel). No persistent changes in E_Cl, E_K or V_m were observed. The full time-dependent ionic and water fluxes for the experiment are shown in Figure 4—figure supplement 1, which shows that the inward flux of impermeant anions causes fluxes of all other ions. (D) Similarly, the addition of extracellular impermeant anions in an osmoneutral manner causes transient shifts in the permeable ion gradients (top panel, colours as in ‘C’), and sustained changes in cellular volume (black, middle panel) as well as the extracellular X and extracellular and intracellular Cl^- concentrations.

https://doi.org/10.7554/eLife.39575.013

Figure 4—figure supplement 1

Download asset Open asset

Flux of all ions for the manipulation in Figure 4c.

Detailed fluxes of all ions and changes to volume during addition of impermeant anions of average charge. The flux of ions across the membrane for potassium (purple, K⁺), sodium (pink, Na⁺), chloride (green, Cl^-) and impermeant anions (orange, (X) are shown in the top panel. Concentration and reversal potential changes are depicted in the second and third panels from the top. Whilst there is a large inward (intracellular) flux of impermeant anions as dictated by the experiment, this is not solely mitigated by fluxes of chloride, the other intracellular anion, but also by potassium and sodium fluxes as well as the movement of water (ie volume changes). The rate of change in volume (representing water movement) is shown in the fourth panel from the top while the bottom panel depicts volume itself.

https://doi.org/10.7554/eLife.39575.014

Figure 5 with 1 supplement

Download asset Open asset

Adjusting the mean charge of impermeant anions shifts the chloride reversal potential with negligible effects on the driving force for chloride.

(A) Decreasing the mean charge of impermeant anions from −0.85 (the default value) to −1 (orange, bottom panel), without changing the absolute number of intracellular impermeant anions caused a persistent decrease in E_Cl (green, top panel), E_K (purple) and V_m with moderate increases in volume (middle panel) in our default single compartment model. Negligible changes in Cl^-driving force (∆DF = 0.16 mV, red) were observed. (B) Analytic solution (solid lines) for different impermeant anion mean charge (z) exactly matches the steady state values from numerical, time series runs (dots) based on adjusting z as in ‘A’. Steady state E_Cl, E_K, V_m (top panel) and [X]_i (middle) increased with increasing z, while changes in z resulted in very small changes in Cl^- DF (bottom, red). The vertical dashed line indicates the values at the chosen default z of −0.85. (C) Influx of a species of impermeant anions with a charge of −1.5, that decreased the mean charge z (bottom panel) from −0.85 to −1 and increased the number of impermeant anions also caused persistent decreases of E_Cl, E_K and V_m as in ‘A’, but with larger increases in cell volume (middle panel). Again, very small persistent changes in Cl^- DF were observed. The volume changes for different flux amounts and charges are illustrated in Figure 5—figure supplement 1. (D) Top, schematic of the experimental setup where whole-cell recordings were made from CA3 pyramidal cells in mouse organotypic brain slices. Impermeant anions (orange) were delivered via electroporation of the negatively charged fluorescent dextran Alexa Flour 488 via a pipette positioned near the soma of the recorded cell. GABA_AR currents were elicited via photo-activation (100 ms, 470 nm LED via objective) of ChR2-expressing GAD2+ interneurons (green cells) in the presence of 5 μM CGP-35348 to block GABA_BRs. Lower trace, current clamp recording showing V_m changes during electroporation of anionic dextran. Confocal image demonstrating cell-localized fluorescence of the anionic dextran electroporated in ‘E’. (E) Top, widefield images with electroporation pipette (orange dashed lines) and the recording pipette (white dashed lines). Note increased fluorescence in the soma after electroporation. Below, insets show GABA_AR currents evoked by photo-activation of GAD2+ interneurons at different holding potentials. Calibration: 1 s, 100 pA. Holding current (reflecting membrane current) and total current (reflecting membrane current plus the GABA_AR current) were measured at the points indicated by the vertical grey and black lines, respectively. IV plots were used to calculate V_m, E_GABA and DF before (left) and after (right) electroporation. (F) Top, population data showing significant decreases in mean V_m and E_GABA but not DF five minutes after electroporation. Below: changes in DF over time. Point at which electroporation occurred marked with orange arrow. ‘ns’, non-significant; *p<0.05. The data for ‘D’, ‘E’ and ‘F’ is provided in Figure 5—source data 1.

https://doi.org/10.7554/eLife.39575.015

Figure 5—source data 1 Source data for Figure 5D–F. Resting membrane potential, E_GABA and driving force at time points as indicated relative to electroporation of anionic dextran (at time 0 min), for each experiment.: https://doi.org/10.7554/eLife.39575.017
Download elife-39575-fig5-data1-v2.xlsx

Figure 5—figure supplement 1

Download asset Open asset

Large volume shifts occur when average impermeant charge is changed via temporary anion species flux across the membrane.

Repeating the methodology in Figure 5C for ion species with different charges (z_flux) demonstrates the effect of driving shifts in the average intracellular impermeant charge on cell volume. Generating even small shifts in average z by adding or removing impermeant anions can cause large changes in cell volume.

https://doi.org/10.7554/eLife.39575.016

Figure 6 with 1 supplement

Download asset Open asset

Impermeant anions drive small shifts in chloride driving force by modifying the Na⁺/K⁺-ATPase pump rate under conditions of active chloride cotransport.

(A) E_Cl (green), E_K (purple) and V_m (black) (top panel), [Na⁺]_i (pink, middle panel) and mean impermeant anion charge z (orange, lower panel) over time in the default single compartment model. Changing z from −0.85 to −1 generated small, persistent Cl^- driving force (DF) shifts (arrows, red) only when the effective ATPase pump rate (J_p) was variable (solid line) and not when J_p was kept constant (dashed line). (B) Solving analytically across different values of z with either a variable J_p (solid lines) or a constant J_p (dashed line), demonstrates the direct relationship between Na⁺ (pink), effective pump rate (J_p) (top panel) and DF (lower panel, red). (C) E_Cl (green), E_K (purple) and V_m (black) (top panel), [Na⁺]_i (pink, middle panel) and cell volume (black, lower panel) over time in the default single-compartment model. Impermeant anions of the same charge as the mean charge of the cell were added. A volume constraint was incorporated by adding a hydrostatic force dependent on membrane tension (dashed lines), which resulted in an impermeable anion-induced transmembrane osmotic differential. This caused a small change in DF when J_p was variable (dashed line), but not when J_p was held constant (dotted lines). (D) Solving analytically across osmolarity differences demonstrates the direct relationship between Na⁺ (pink), effective pump rate (J_p) (top panel) and DF (lower panel). Note, the small changes in DF. (E) Schematic explaining the mechanism through which impermeant anion-induced cell swelling in the presence of volume constraints (i.e. membrane tension) result in steady states with equal but non-zero osmotic (osmo_P) and hydrostatic pressures (H_P), causing transmembrane osmotic differences (t₁). This causes small changes in Na⁺, and hence J_p. (F) All Na⁺/K⁺-ATPase pump rate-related shifts in the DF require KCC2 activity; in the absence of activity (dashes), no shifts in driving force can occur. In Figure 6—figure supplement 1, we show that the results in (A) and (C) are similar when an experimentally-matched model of the Na⁺/K⁺-ATPase is used.

https://doi.org/10.7554/eLife.39575.018

Figure 6—figure supplement 1

Download asset Open asset

An ATPase model with different kinetics has similar properties to our model’s ATPase when average impermeant anion charge in the cell is changed, and is also dependent on the ATPase pump rate’s dependence on Na⁺.

(A) The experiment for Figure 6A was repeated using the usual ATPase model (solid line) compared to Hamada et al. (2003)’s model and is consistent with the results in the main text. E_Cl (green), E_K (purple) and V_m (black) (top panel), [Na⁺]_i (pink, middle panel) and average impermeant anion charge z (orange, lower panel) over time in the default single compartment model. Changing the average z from −0.85 to −1 generated small, persistent Cl^- driving force shifts (arrows) only when the effective ATPase pump rate (J_p) was variable (solid and dashed line) for both pump models, and not when J_p was kept constant (dotted line). (B) The experiment for Figure 6C was repeated comparing the usual ATPase model (solid line) with Hamada et al., 2003’s model and is consistent with the results in the main text. E_Cl (green), E_K (purple) and V_m (black) (top panel), [Na⁺]_i (pink, middle panel) and cell volume (black, lower panel) over time in the default single-compartment model. Impermeant anions of the same average charge of the cell were added. A volume constraint on Hamada et al.’s model was incorporated by adding a hydrostatic force dependent on membrane tension (dotted lines), which resulted in an impermeable anion induced transmembrane osmotic differential. This caused a small change in Cl^- driving force when J_p was variable (dotted line), but not when J_p was held constant (dash-dotted lines).

https://doi.org/10.7554/eLife.39575.019

Figure 7

Download asset Open asset

Local changes in KCC2 activity generate local differences in chloride reversal and driving force only under conditions of constrained chloride diffusion.

(A) Schematic depicting the multi-compartment model representing a virtual dendrite of length 100 μm. The virtual dendrite consists of 10 compartments of length 10 μm and initial radius 0.5 μm. Each compartment contains the same mechanisms and default parameterization as the single compartment model. All ions, except impermeant anions, could move between compartments by electrodiffusion. (B) Top panel, E_Cl (green), E_K (purple), V_m (black) and DF (arrows, red) from the second from top compartment (indicated with a white triangle) where the conductance of KCC2 was increased. The insets depict the diameter, and absolute change from baseline of E_Cl, V_m and DF for all compartments before (i), during (ii) and after (iii) the activity of KCC2 was selectively increased. This resulted in E_Cl decreasing in all compartments with minimal changes to E_K and V. Consequently, the Cl^- DF (red) increased. In this case the diffusion constant for Cl^- (D_Cl = 2.03×10⁻⁷ dm².s⁻¹) resulted in E_Cl and DF changes being widespread across the virtual dendrite. (C) Reducing the Cl^- diffusion constant to 0.2 × 10⁻⁷ dm².s⁻¹ resulted in a localized effect of compartment-specific KCC2 activity increases on E_Cl and DF.

https://doi.org/10.7554/eLife.39575.020

Figure 8

Download asset Open asset

Local changes in impermeant anions do not establish the local driving force for chloride.

(A) Multi-compartment model of a 10 compartment 100 μm virtual dendrite as in Figure 7. (B) Top panel, E_Cl (green), E_K (purple), V_m (black) and DF (arrows, red) within the compartment where additional impermeant anions were exclusively added (indicated with a white triangle). The insets depict the diameter, and absolute change from baseline for E_Cl, V_m and DF for each compartment of the virtual dendrite before (i), during (ii) and after (iii) impermeant anions were added. The selective addition of impermeant anions of default charge (z = −0.85) to the 2^nd from top compartment resulted in transient but non-permanent shifts in E_Cl, V_m and the Cl^- DF in all compartments. The volume of the compartment where impermeant anion addition occurred increased permanently. (C) Traces and insets as in ‘B’ showing the addition of impermeant anions of more negative charge in order to decrease z in the second from top compartment specifically. Note that during addition of impermeant anions, E_Cl, E_K, V_m changed. We also observed persistent decreases in E_Cl and V_m in the compartment manipulated, with a negligible change in DF (0.01 mV). Again, impermeant anion addition resulted in an increase in the volume of the specific dendritic compartment manipulated.

https://doi.org/10.7554/eLife.39575.021

Tables

Table 1

Constants, default parameters and usual steady state values for variables used in the biophysical models.

https://doi.org/10.7554/eLife.39575.022

	Value	Description
Constants
F	96485.33 C/mol	Faraday constant
R	8.31446 J/(K.mol)	Universal gas constant
T	310.15 K	Absolute temperature (=37°C)
Parameters
C_m	2 $\times$ 10⁻⁶ F/cm²	Unit membrane capacitance (Qian and Sejnowski, 1989)
g_Na	20 µS/cm²	Na⁺ leak conductance (Kager et al., 2000)
g_K	70 µS/cm²	K⁺ leak conductance (Kager et al., 2000)
g_Cl	20 µS/cm²	Cl^- leak conductance
g_KCC2	20 µS/cm²	KCC2 conductance (Doyon et al., 2016)
v_w	0.018 dm³/mol	Partial molar volume of water (Hernández and Cristina, 1998)
p_w	0.0015 dm/s	Osmotic permeability (Hernández and Cristina, 1998)
k_m	25 N/dm	Variable for membrane tension (higher than reported (Dai et al., 1998): used in this paper to accentuate differences in osmolarity)
p	0.1 C/(dm².s)	Default pump rate constant
[Na⁺]_o	145 mM	Extracellular Na⁺ concentration
[K⁺]_o	3.5 mM	Extracellular K⁺ concentration
[Cl^-]_o	119 mM	Extracellular Cl^- concentration
[X^-]_o	29.5 mM	Extracellular impermeant anion (X) concentration (Lodish et al., 2009; Raimondo et al., 2015)
D_Na	1.33 $\times$ 10⁻⁷ dm²/s	Na⁺ diffusion constant (Hille, 2001)
D_K	1.96 $\times$ 10⁻⁷ dm²/s	K⁺ diffusion constant (Hille, 2001)
D_Cl	2.03 $\times$ 10⁻⁷ dm²/s	Cl^- diffusion constant (Hille, 2001)
Variables (default steady state)
V_m	−72.6 mV	Membrane potential
[Na⁺]_i	14.0 mM	Intracellular Na⁺ concentration
[K⁺]_i	122.9 mM	Intracellular K⁺ concentration
[Cl^-]_i	5.2 mM	Intracellular Cl^- concentration
[X^z]_i	154.9 mM	Intracellular impermeant anion (X) concentration
z	−0.85	Mean charge of intracellular X (Lodish et al., 2009; Raimondo et al., 2015)
w	2.0 pL (single compartment) 0.078 pL (multi-compartment)	Volume

Additional files

Supplementary file 1 Analytical solution derivation.: https://doi.org/10.7554/eLife.39575.023
Download elife-39575-supp1-v2.pdf
Supplementary file 2 Data and scoring methodology for the meta-analysis in Figure 3F.: https://doi.org/10.7554/eLife.39575.024
Download elife-39575-supp2-v2.docx
Transparent reporting form: https://doi.org/10.7554/eLife.39575.025
Download elife-39575-transrepform-v2.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Mendeley

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Kira M Düsterwald
Christopher B Currin
Richard J Burman
Colin J Akerman
Alan R Kay
Joseph V Raimondo

(2018)

Biophysical models reveal the relative importance of transporter proteins and impermeant anions in chloride homeostasis

eLife 7:e39575.

https://doi.org/10.7554/eLife.39575

Share this article

Cite this article

A biophysical model of ion dynamics based on the pump-leak mechanism demonstrates the importance of the sodium-potassium ATPase for setting transmembrane ion gradients including chloride.

Different models of the sodium-potassium ATPase have distinct kinetic properties but produce similar qualitative effects.

Changes in Na+/K+-ATPase activity near the model’s default pump rate produce minimal changes in ECl and Vm.

Membrane conductances affect steady-state intracellular chloride concentration only in the presence of cation-chloride cotransport.

The classic dependence of membrane potential (mV) on the logarithm of extracellular potassium is reproduced using our model.

K+flux is balanced by active and passive fluxes, dependent on Na+steady states.

Cation-chloride cotransport sets the chloride reversal and driving force for transmembrane chloride flux.

Figure 3—source data 1

Models of KCC2 have different kinetic properties.

KCC2 and ATPase models with different kinetics have similar properties when KCC2 is upregulated.

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

Flux of all ions for the manipulation in Figure 4c.

Adjusting the mean charge of impermeant anions shifts the chloride reversal potential with negligible effects on the driving force for chloride.

Figure 5—source data 1

Large volume shifts occur when average impermeant charge is changed via temporary anion species flux across the membrane.

Impermeant anions drive small shifts in chloride driving force by modifying the Na+/K+-ATPase pump rate under conditions of active chloride cotransport.

An ATPase model with different kinetics has similar properties to our model’s ATPase when average impermeant anion charge in the cell is changed, and is also dependent on the ATPase pump rate’s dependence on Na+.

Local changes in KCC2 activity generate local differences in chloride reversal and driving force only under conditions of constrained chloride diffusion.

Local changes in impermeant anions do not establish the local driving force for chloride.

Constants, default parameters and usual steady state values for variables used in the biophysical models.

Supplementary file 1

Supplementary file 2

Transparent reporting form

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

Changes in Na⁺/K⁺-ATPase activity near the model’s default pump rate produce minimal changes in E_Cl and V_m.

K⁺flux is balanced by active and passive fluxes, dependent on Na⁺steady states.

Impermeant anions drive small shifts in chloride driving force by modifying the Na⁺/K⁺-ATPase pump rate under conditions of active chloride cotransport.

An ATPase model with different kinetics has similar properties to our model’s ATPase when average impermeant anion charge in the cell is changed, and is also dependent on the ATPase pump rate’s dependence on Na⁺.