Minimal requirements for a neuron to coregulate many properties and the implications for ion channel correlations and robustness

Abstract
Editor's evaluation
Introduction
Results
Discussion
Methods
Data availability
References
Article and author information
Metrics

Abstract

Neurons regulate their excitability by adjusting their ion channel levels. Degeneracy – achieving equivalent outcomes (excitability) using different solutions (channel combinations) – facilitates this regulation by enabling a disruptive change in one channel to be offset by compensatory changes in other channels. But neurons must coregulate many properties. Pleiotropy – the impact of one channel on more than one property – complicates regulation because a compensatory ion channel change that restores one property to its target value often disrupts other properties. How then does a neuron simultaneously regulate multiple properties? Here, we demonstrate that of the many channel combinations producing the target value for one property (the single-output solution set), few combinations produce the target value for other properties. Combinations producing the target value for two or more properties (the multioutput solution set) correspond to the intersection between single-output solution sets. Properties can be effectively coregulated only if the number of adjustable channels (n_in) exceeds the number of regulated properties (n_out). Ion channel correlations emerge during homeostatic regulation when the dimensionality of solution space (n_in − n_out) is low. Even if each property can be regulated to its target value when considered in isolation, regulation as a whole fails if single-output solution sets do not intersect. Our results also highlight that ion channels must be coadjusted with different ratios to regulate different properties, which suggests that each error signal drives modulatory changes independently, despite those changes ultimately affecting the same ion channels.

Editor's evaluation

Neurons develop and maintain rich electrophysiological properties that enable nervous systems to function. The question of how different neural properties are regulated by internal ion channel expression mechanisms remains unresolved. In this paper, Yang and colleagues address the question from an abstract perspective by asking how multiple constraints on the physiological properties of neurons, such as firing rate curves and energy efficiency, narrow down the available regulatory possibilities. Their results from a mixture of modelling and dynamic clamp experiments point to the existence of multiple parallel internal feedback loops for controlling ion channel expression in neurons, and derive conditions under which co-regulation mechanisms will fail.

https://doi.org/10.7554/eLife.72875.sa0

Introduction

Neurons maintain their average firing rate near a target value, or set point, by adjusting their intrinsic excitability and synaptic weights (Turrigiano et al., 1994; O’Leary et al., 2010; Aizenman et al., 2003; Desai et al., 1999; Hengen et al., 2013; van Welie et al., 2006). Homeostatic regulation of intrinsic excitability is achieved through feedback control of diverse ion channels (Turrigiano et al., 1994; O’Leary et al., 2010; Desai et al., 1999; Joseph and Turrigiano, 2017). Computational models have successfully employed negative feedback to adjust ion channel densities (O’Leary et al., 2014; O’Leary and Marder, 2016; Liu et al., 1998; Olypher and Prinz, 2010; LeMasson et al., 1993) and control theory provides a valuable framework to conceptualize how this occurs (O’Leary and Wyllie, 2011). But most of the mechanistic details remain unclear (Davis, 2006; Turrigiano, 2011) and are not straightforward; for instance, different perturbations can trigger similar changes in excitability via different signaling pathways affecting different ion channels (Kulik et al., 2019), or via different combinations of excitability changes and synaptic scaling (Maffei and Turrigiano, 2008).

Firing rate homeostasis is facilitated by the ability of different ion channel combinations to produce equivalent excitability (Marder, 2011). The ability of distinct elements to produce the same outcome is known as degeneracy (Edelman and Gally, 2001) and has attracted increasing attention in neuroscience (Marder, 2011; Marder and Goaillard, 2006; Cropper et al., 2016; O’Leary, 2018; Ratté and Prescott, 2016; Tononi et al., 1999; Goaillard and Dufour, 2014; Mason et al., 2015; Rathour and Narayanan, 2019). Many aspects of neural function at the genetic (Klassen et al., 2011; Trojanowski et al., 2014), synaptic (Mukunda and Narayanan, 2017; Anirudhan and Narayanan, 2015), cellular (Kim et al., 2017; Taylor et al., 2009; Drion et al., 2015; Ratté et al., 2014b; Mittal and Narayanan, 2018; Jain and Narayanan, 2020; Günay et al., 2008; Migliore et al., 2018), and network (Grashow et al., 2010; Marder and Taylor, 2011; Prinz et al., 2004; Price and Friston, 2002; Onasch and Gjorgjieva, 2020) levels are now recognized as being degenerate. Ion channel degeneracy facilitates robust homeostatic regulation of neuronal excitability by enabling a disruptive change in one ion channel to be offset by compensatory changes in other ion channels (O’Leary, 2018; Drion et al., 2015; Ratté et al., 2014b; Rho and Prescott, 2012; Swensen and Bean, 2005; Zhao and Golowasch, 2012; Achard and De Schutter, 2006; Bhalla and Bower, 1993; Olypher and Calabrese, 2007).

But to maintain good neural coding, neurons must regulate various aspects of their activity beyond average firing rate (Stemmler and Koch, 1999), and must do so while also regulating their energy usage, osmolarity, pH, protein levels, etc. (Frere and Slutsky, 2018; Fisher et al., 2010; Balch et al., 2008). These other cellular properties affect and are affected by neuronal activity, which is to say that properties are not regulated in isolation from one another. And just as each property depends on multiple ion channels, each ion channel can affect multiple properties (Taylor et al., 2009; Günay et al., 2008). This pleiotropy or functional overlap (Goaillard and Marder, 2021) confounds the independent regulation of each property. In response to a perturbation, a compensatory ion channel change that restores one property to its target value may exacerbate rather than mitigate the disturbance of a second property (Figure 1A). In theory, both properties could be regulated by adjusting two or more ion channels (Figure 1B), but this suggests that the number of properties that can be coregulated is limited by the number of adjustable ion channels, even if the channels are pleiotropic.

Figure 1

Download asset Open asset

Simultaneous regulation of >1 property in a theoretical neuron.

(A) Challenge: If two ion channels, $g_{M}$ (pink) and $g_{A H P}$ (purple), both affect firing rate (orange), then the change in firing rate caused by perturbing $g_{M}$ can be offset by a compensatory change in $g_{A H P}$ . However, if $g_{M}$ and $g_{A H P}$ also affect firing pattern (blue), but not in exactly the same way, then the compensatory change in $g_{A H P}$ that restores firing rate to its target value may exacerbate rather than resolve the change in firing pattern. (B) Solution: Adjusting $g_{A H P}$ and at least one additional conductance, $g_{x}$ (green), that affects each property in a different way than $g_{A H P}$ , may enable firing rate and firing pattern to be regulated back to their target values. This suggests that the number of adjustable ion channels (‘dials’) relative to the number of regulated properties is important. Note that adjusting $g_{x}$ to restore the firing pattern necessitates a small extra increase in $g_{A H P}$ (compare pale and dark purple); in other words, conductances must be coadjusted.

To identify the requirements for coregulating many properties, we began with an experiment to confirm that regulating one property by adjusting a single channel is liable to disrupt other properties. We then proceeded with modeling to unravel how multiple channels must be coadjusted to ensure proper regulation of >1 property. We show that the number of adjustable ion channels (n_in) must exceed the number of regulated properties (n_out). This is consistent with past work (Olypher and Prinz, 2010; Olypher and Calabrese, 2007; Foster et al., 1993) but the implications for homeostatic regulation have not been explored. To that end, we show that the dimensionality of solution space (n_in − n_out) influences the emergence of ion channel correlations in the presence or absence of noise, and that increased correlations presage regulation failure. We also show that regulating different properties requires that ion channels are coadjusted with different ratios, which necessitates separate master regulators.

Results

Adjusting an ion channel to regulate one property risks disrupting other properties – experiments

To test experimentally if regulating one property by adjusting one ion channel is liable to disrupt a second property, we disrupted the firing rate of a CA1 pyramidal neuron by blocking its voltage-gated M-type K⁺ current (I_M) and then we restored firing rate to its baseline value by inserting a virtual calcium-activated AHP-type K⁺ current (I_AHP) using dynamic clamp, while also monitoring the firing pattern. Specifically, the neuron was stimulated by injecting irregularly fluctuating (noisy) current under four conditions (Figure 2A, B): at baseline (blue), after blocking native I_M (black), and again after introducing virtual I_M (cyan) or virtual I_AHP (red). Adding virtual I_M demonstrates that we can replace a native current with an equivalent virtual current; compensation was modeled by inserting a distinct current with functional overlap, namely I_AHP (see Figure 1). Inserting virtual I_M or virtual I_AHP reversed the depolarization (Figure 2C), spike amplitude attenuation (Figure 2D), and firing rate increase (Figure 2E) caused by blocking native I_M. Replacing native I_M with virtual I_M did not affect firing pattern, quantified as the coefficient of variation of the interspike interval (CV_ISI), whereas virtual I_AHP reduced CV_ISI (Figure 2F), often causing the neuron to spike at different times (in response to different stimulus fluctuations) than the neuron with native or virtual I_M (Figure 2G). The results confirm predictions from Figure 1, namely that compensatory ‘upregulation’ of I_AHP restored firing rate but disrupted firing pattern.

Figure 2

Download asset Open asset

A compensatory change that restores firing rate disrupts firing pattern in a CA1 pyramidal neuron.

(A) Rasters show spiking in a single CA1 pyramidal neuron given the same noisy current injection (gray) on five trials under each of four conditions: baseline (blue), after blocking native I_M with 10 μM XE991 (black), and again after inserting virtual I_M (cyan) or virtual I_AHP (red) using dynamic clamp. See Methods for virtual conductance parameters. (B) Sample voltage traces under each condition. (C) Membrane potential differed across conditions (F_3,16 = 298.49, p < 0.001, one-way ANOVA); specifically, it was depolarized by blocking I_M (t = 18.71, p < 0.001, Tukey test) but that effect was reversed by inserting virtual I_M (t = 24.56, p < 0.001) or virtual I_AHP (t = 27.00, p < 0.001). (D) Spike amplitude also differed across conditions (F_3,16 = 154.04, p < 0.001); specifically, it was attenuated by blocking I_M (t = 20.23, p < 0.001) but that effect was reversed by inserting virtual I_M (t = 13.99, p < 0.001) or virtual I_AHP (t = 16.09, p < 0.001). (E) Firing rate also differed across conditions (F_3,16 = 177.61, p < 0.001); specifically, it was increased by blocking I_M (t = 18.38, p < 0.001) but that effect was reversed by inserting virtual I_M (t = 20.85, p < 0.001) or virtual I_AHP (t = 16.12, p < 0.001). (F) Regularity of spiking, reflected in the coefficient of variation of the interspike interval, also differed across conditions (F_3,16 = 188.23, p < 0.001), but whereas blocking I_M had a modest effect (t = 2.69, p = 0.048) and inserting virtual I_M had no effect (t = 0.39, p = 1.00), inserting virtual I_AHP had a large effect (t = 18.23, p < 0.001). Data are summarized as mean ± standard deviation (SD). Each data point represents a different trial from a single neuron. (G) Enlarged view of rasters to highlight spikes that occurred with native or virtual I_M but not with virtual I_AHP (blue shading) or vice versa (red shading) to illustrate the change in spike pattern caused by correcting the change in firing rate caused by blocking I_M with a compensatory change in I_AHP.

Figure 2—source data 1 Numerical values for experimental data plotted in Figure 2. Data include all spikes times (for panel A) as well as trial-averaged values of membrane potential, spike amplitude, firing rate, and coefficient of variation of the interspike interval (for panels C–F, respectively).: https://cdn.elifesciences.org/articles/72875/elife-72875-fig2-data1-v2.xlsx
Download elife-72875-fig2-data1-v2.xlsx

Few of the ion channel combinations producing the target value for one property also produce the target value for a second property

We proceeded with computational modeling to identify the conditions required to coregulate multiple properties. For our investigation, it was critical to account for all ion channels that are changing. Even for well characterized neurons like pyramidal cells, there is no comprehensive list of which properties are directly regulated and which ion channel are involved in each case. Therefore, we built a simple model neuron in which >1 property can be coregulated by adjusting a small number of ion channels. The relative rates at which ion channel densities are updated using the ‘error’ in each property differ between properties, meaning ion channels are coadjusted according to a certain ratio to regulate one property (O’Leary et al., 2014) but are coadjusted according to a different ratio to regulate a second property, which has some notable implications (see Discussion). We focused on regulation of independent properties like firing rate and energy efficiency per spike – instead of energy consumption rate, for example, which depends on firing rate – but our conclusions do not hinge on which properties are considered. Target values were chosen arbitrarily.

Simulations were conducted in a single-compartment model neuron whose spikes are generated by a fast sodium conductance and a delayed rectifier potassium conductance with fixed densities. A spike-dependent adaptation mechanism, g_AHP was included at a fixed density. Densities of all other channels were either systematically varied (and values producing the target output were selected) or they were adjusted via negative feedback to produce the target output (see Methods). The former approach, which amounts to a grid search, identifies all density combinations producing a target output (i.e., the solution set). The latter approach finds a subset of those combinations through a homeostatic regulation mechanism. Ion channels with adjustable densities included a sodium conductance ( $g_{N a}$ ) and potassium conductance ( $g_{K}$ ) that both activate with kinetics similar to the delayed rectifier channel, a slower-activating M-type potassium conductance ( $g_{M}$ ), and a leak conductance ( $g_{l e a k}$ ).

We began by testing whether ion channel density combinations that produce the target value for one property also produce a consistent value for other properties. Figure 3A shows rheobase for different combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ . The contour from point a ( ${\bar{g}}_{N a}$ = 1.54 mS/cm², ${\bar{g}}_{K}$ = 0 mS/cm²) to point b ( ${\bar{g}}_{N a}$ = 4 mS/cm², ${\bar{g}}_{K}$ = 2.95 mS/cm²) represents all combinations (i.e., the solution set) yielding a rheobase of 30 μA/cm². Despite yielding the same rheobase, density combinations along the a–b contour did not yield the same minimum sustainable firing rate (f_min) (Figure 3B). The value of f_min reflects spike initiation dynamics: f_min >> 0 spk/s is consistent with class 2 excitability (Hodgkin, 1948) and operation as a coincidence detector (Ratté et al., 2014a), whereas f_min≈ 0 spk/s is consistent with class 1 excitability and operation as an integrator, with several consequences (Figure 3—figure supplement 1). After adding $g_{A H P}$ to expand the stimulus range over which fluctuation-driven spikes occur, we plotted the firing rate driven by irregularly fluctuating (noisy) I_stim for different combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ (Figure 3C). Density combinations yielding the same rheobase did not yield equivalent stimulation-evoked firing rates (Figure 3D).

Figure 3 with 2 supplements see all

Download asset Open asset

Most ion channel combinations producing the target value for one property produce inconsistent values for other properties.

(A) Color shows the minimum I_stim required to evoke repetitive spiking (rheobase) for different combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ . Contour linking a and b highlights density combinations yielding a rheobase of 30 μA/cm². (B) Minimum sustainable firing rate (f_min) varied along the iso-rheobase contour (see Figure 3—figure supplement 1). (C) Color shows firing rate evoked by noisy I_stim (τ_stim = 5ms, σ_stim = 10 μA/cm², µ_stim = 40 μA/cm²) for different combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ . Adding spike-dependent and -independent forms of adaptation ( ${\bar{g}}_{A H P}$ = 1.75 mS/cm² and ${\bar{g}}_{M}$ = 0.5 mS/cm²) broadened the dynamic range. Gray curve shows density combinations yielding a firing rate of 40 spk/s (i.e., an iso-firing rate contour; shown in red in Figures 4—6). Insets show sample responses to equivalent noisy stimulation. (D) Firing rate varied along the iso-rheobase contour from panel A. (E) Color shows energy consumption rate for different combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ based on firing rates shown in panel C. Iso-firing rate contour *c–d* (gray line in panel C) does not align with energy contours. ATP consumed by the Na⁺/K⁺ pump was calculated from the total Na⁺ influx and K⁺ efflux determined from the corresponding currents. (F) Energy consumption rate varied along the iso-firing rate contour. (G) Color shows energy efficiency per spike for different combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ . Energy efficiency was calculated as the capacitive minimum, CΔV, divided by total Na⁺ influx, where C is capacitance and ΔV is spike amplitude (see Figure 3—figure supplement 2). Density combinations along contour *i–ii* (dashed line) yield energy efficiency of 23.5% (shown in green in Figure 5). (H) Both rheobase and firing rate varied along the iso-energy efficiency contour.

Having demonstrated that ion channel density combinations yielding the target value for one aspect of excitability (e.g., rheobase) yield differing values for other aspects of excitability (e.g., f_min or firing rate), we predicted that the same lack of generalization would extend to other cellular properties (e.g., energy efficiency). Spikes are energetically costly (Attwell and Laughlin, 2001) but vary in their energy efficiency based on the temporal overlap in sodium and potassium channel activation (Sengupta et al., 2010; Figure 3—figure supplement 2). Using responses reported in Figure 3C, we measured energy consumption rate for combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ (Figure 3E). Energy consumption rate increased with firing rate, as expected, but did not increase equivalently across all density combinations, as evident from the variation in energy consumption rate along the iso-firing rate contour (Figure 3F). This is due to differences in energy efficiency (Figure 3G) determined as the energy consumed per spike relative to the theoretical minimum (see Methods). Density combinations yielding equally efficient spikes yielded different stimulation-evoked firing rates and rheobase values (Figure 3H). One might presume that spikes are produced as efficiently as possible rather than being regulated to a specific value, but energy efficiency decreases propagation safety factor (Al-Basha and Prescott, 2019) and a target value likely emerges by balancing these competing interests (Speakman et al., 2011).

Adjusting an ion channel to regulate one property risks disrupting other properties – simulations

Like the experiment in Figure 2, we tested if restoring firing rate to its target value via compensatory changes in either of two ion channels (Figure 4A) after perturbing a third channel would disrupt energy efficiency in our model neuron. Gray circles on Figure 4B show randomly chosen combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{l e a k}$ that yield a firing rate of 40 spk/s when ${\bar{g}}_{K}$ = 2 mS/cm². When ${\bar{g}}_{K}$ was ‘blocked’ (abruptly reset to 0 mS/cm²), firing rate jumped to ~93 spk/s before being restored to 40 spk/s via feedback control of either ${\bar{g}}_{N a}$ (pink) or ${\bar{g}}_{l e a k}$ (cyan) (Figure 4C). As ${\bar{g}}_{N a}$ or ${\bar{g}}_{l e a k}$ converged on new, compensated densities, firing rate returned to its target value but energy efficiency was affected in opposite ways (Figure 4D). If energy efficiency is also regulated, then ${\bar{g}}_{N a}$ and ${\bar{g}}_{l e a k}$ must be coadjusted in a way that restores firing rate without disrupting energy efficiency.

Figure 4

Download asset Open asset

Ion channel changes mediating the same effect on firing rate can oppositely affect energy efficiency.

(A) Schematic shows how a difference in firing rate from its target value creates an error that is reduced by updating ${\bar{g}}_{N a}$ or ${\bar{g}}_{l e a k}$ . (B) Iso-firing rate contours for 40 spk/s are shown for different combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{l e a k}$ with ${\bar{g}}_{K}$ at baseline (2 mS/cm², dashed red curve) and after ${\bar{g}}_{K}$ was ‘knocked out’ (0 mS/cm², solid red curve). When ${\bar{g}}_{K}$ was abruptly reduced, starting models (gray dots) spiked rapidly (~93 spk/s) before firing rate was regulated back to its target value by compensatory changes in either ${\bar{g}}_{N a}$ (pink) or ${\bar{g}}_{l e a k}$ (cyan). Models evolved in different directions and settled at different positions along the solid curve. (C) Trajectories show evolution of ${\bar{g}}_{N a}$ and ${\bar{g}}_{l e a k}$ . Trajectories are terminated once target firing rate is reached. Sample traces show responses before (gray) and immediately after (black) ${\bar{g}}_{K}$ was reduced, and again after compensatory changes in ${\bar{g}}_{N a}$ (pink) or ${\bar{g}}_{l e a k}$ (cyan). (D) Distributions of energy efficiency are shown before (gray) and after firing rate regulation via control of ${\bar{g}}_{N a}$ (pink) or ${\bar{g}}_{l e a k}$ (cyan).

Geometrical explanation for the relationship between single- and multioutput solutions

Next, we sought a geometrical explanation for how to coadjust >1 ion channel to coregulate >1 property. The top panel of Figure 5A shows all combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ (n_in = 2) yielding the target value for firing rate (red) or energy efficiency (green) (n_out = 1). Like in Figures 3 and 4, the solution set for each property corresponds to a curve. The multioutput solution set for firing rate and energy efficiency (n_out = 2) corresponds to where the two single-output solution sets intersect, which occurs at a point – only one density combination yields the target values for both properties, meaning the multioutput solution is unique. But if another conductance like g_leak is adjustable (n_in = 3), the curves in 2D parameter space (top panel) transform into surfaces in 3D parameter space (bottom panel) and those surfaces intersect along a curve – many density combinations yield the target values for both properties, meaning the multioutput solution is degenerate. The same patterns are evident for firing rate and input resistance (Figure 5B) or energy efficiency and input resistance (Figure 5C). Figure 5D shows that if all three properties – firing rate, energy efficiency and input resistance – are regulated (n_out = 3), the multioutput solution set is empty for n_in = 2 (i.e., the three curves do not intersect at a common point [top panel] unless ${\bar{g}}_{l e a k}$ is reset to 1.95 mS/cm² [inset]) and is unique for n_in = 3 (bottom panel). From this, we conclude that robust regulation of n properties requires >n adjustable ion channels even if each ion channel contributes to regulation of >1 property. This is like a system of linear equations, which is said to be underdetermined if unknowns (inputs) outnumber equations (outputs) (see Discussion).

Figure 5 with 1 supplement see all

Download asset Open asset

A degenerate solution for n properties requires at least n + 1 adjustable ion channels.

Curves in top panels depict single-output solutions sets based on control of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ (n_in = 2); ${\bar{g}}_{l e a k}$ = 2 mS/cm². Surfaces in bottom panels depict single-output solution sets based on control of ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{l e a k}$ (n_in = 3). Intersection (yellow) of single-output solutions at a point constitutes a unique multioutput solution, whereas intersection along a curve (or higher-dimensional manifold) constitutes a degenerate multioutput solution. (A) Curves for firing rate (40 spk/s) and energy efficiency (23.5%) intersect at a point whereas the corresponding surfaces intersect along a curve. Solutions for firing rate and input resistance (0.65 kΩ cm²) (B) and for energy efficiency and input resistance (C) follow the same pattern as in panel A. (D) For n_in = 2 (top), curves for firing rate, energy efficiency and input resistance do not intersect at a common point unless ${\bar{g}}_{l e a k}$ is reset to 1.95 mS/cm² (inset). For n_in = 3 (bottom), the three surfaces intersect at the same point as in the inset. See Figure 5—figure supplement 1 for the effects of tolerance on solution sets.

In all simulations involving homeostatic regulation, properties were regulated to within certain bounds of the target rather than to a precise value, despite how solution sets are depicted. When tolerances are depicted, single-output solution sets correspond to thin strips (rather than curves) or shallow volumes (rather than surfaces) in 2D and 3D parameter space, respectively (Figure 5—figure supplement 1). The width, depth, etc. are proportional to the tolerance. How does this impact multioutput solutions? Thin 2D strips intersect at a small patch (unlike curves intersecting at a point) but that patch in 2D parameter space (Figure 5—figure supplement 1A) is unlike the long 1D curve formed by broad 2D surfaces intersecting in 3D parameter space (see bottom panels of Figure 5). Likewise, shallow 3D volumes intersect at a narrow tube (unlike surfaces intersecting at a curve) but that tube in 3D parameter space (see Figure 5—figure supplement 1B) is unlike the broad 2D surface formed by deep 3D volumes intersecting in 4D parameter space. In other words, increasing tolerance does not increase solution space dimensionality the same way as increasing n_in.

The need for ≥n adjustable channels to regulate n properties has been shown before (Olypher and Prinz, 2010; Olypher and Calabrese, 2007; Foster et al., 1993) and may be obvious to some for mathematical reasons, but the relationship has some important implications (e.g., for parameter estimation Sarkar and Sobie, 2010). The implications for homeostatic regulation have not been thoroughly explored, thus prompting the next steps of our study.

The dimensionality of solution space affects ion channel correlations

We predicted that the dimensionality of solution space – point (0D), curve (1D), surface (2D), volume (3D), etc. – affects ion channel correlations by limiting the degrees of freedom. To explore this, we measured ion channel correlations within single- and multioutput solution sets found through homeostatic regulation. Figure 6A shows all combinations of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ (n_in = 2) producing the target firing rate (n_out = 1). Correlation between ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ is high because homeostatically determined solutions are constrained to fall along a curve; under those conditions, variation in one channel is offset solely by variation in the other channel, and all covariance is thus captured in a single pairwise relationship. If ${\bar{g}}_{M}$ is also allowed to vary (n_in = 3), homeostatically determined solutions spread across a surface and pairwise correlations become predictably weaker (Figure 6B) since variation in one channel can be offset by variation in two other channels, and covariance is diluted across >1 pairwise relationship. If additional channels are allowed to vary (n_in ≥ 4), solutions distribute over higher-dimensional manifolds and pairwise correlations further weaken (Figure 6C). However, if firing rate and energy efficiency are both regulated (n_out = 2; Figure 6D), homeostatically determined solutions are once again constrained to fall along a curve when n_in = 3, and pairwise correlations strengthen (Figure 6E). Increasing n_in to four while keeping n_out at two caused correlations to weaken (Figure 6F). These results confirm the predicted impact of solution space dimensionality on ion channel correlations.

Figure 6

Download asset Open asset

Dimensionality of solution space affects ion channel correlations.

(A) Starting from a normally distributed cluster of ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ (white dots) yielding an average firing rate of 73 spk/s, ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ (n_in = 2) were homeostatically adjusted to regulate firing rate (n_out = 1) to its target value of 40 spk/s. Gray lines show trajectories. Because solutions (gray dots) converge on a curve, the pairwise correlation between ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ is predictably strong. Scatter plots show solutions centered on the mean and normalized by the standard deviation (z-scores). Correlation coefficient (R) is shown in the bottom right corner of each scatter plot. (B) Same as panel A (n_out = 1) but via control of ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{M}$ (n_in = 3). Homeostatically found solutions converge on a surface and ion channel correlations are thus weaker. (C) If ${\bar{g}}_{l e a k}$ is also controlled (n_in = 4), solutions converge on a hard-to-visualize volume (not shown) and pairwise correlations are further weakened. (D) Schematic shows how errors for two regulated properties are combined: the error for each property is calculated separately and is scaled by its respective control rate τ to calculate updates, and all updates for a given ion channel (i.e., originating from each error signal) are summed. (E) Same as panel B (n_in = 3), but for regulation of firing rate and energy efficiency (n_out = 2). Homeostatically found solutions once again converge on a curve (yellow), which now corresponds to the intersection of two surfaces; ion channel correlations are thus strong, like in panel A. (F) If n_in is increased to four while n_out remains at 2, solutions converge on a surface (not shown) and ion channel correlations weaken.

O’Leary et al., 2013 demonstrated how the relative rates at which different ion channel densities are controlled impact their correlation. This is reproduced in Figure 7A, where, from the same initial conditions (density combinations), homeostatic control with different relative rates (shown in pink and cyan) produces solutions with different correlations. Relative regulation rates can affect not only the strength of pairwise correlations, but also the sign (compare correlation between ${\bar{g}}_{N a}$ and ${\bar{g}}_{M}$ ). However, if firing rate and energy efficiency are both homeostatically regulated, correlations strengthen (consistent with results from Figure 6) and become independent of the relative regulation rates because solutions are limited to a lower-dimensional solution space (Figure 7B). Recall that dimensionality of the multioutput solution set corresponds to n_in − n_out. For relative regulation rates to influence ion channel correlations, homeostatically determined solutions must fall on a solution manifold with dimensionality >1, but not >>1 lest pairwise correlations be diluted.

Figure 7

Download asset Open asset

Effect of relative regulation rates on ion channel correlations depends on the dimensionality of solution space.

(A) Homeostatic regulation of firing rate via control of ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{M}$ . Same as Figure 6B, but for two new sets of regulation rates (see Supplementary file 1). Solutions found for each set of rates (pink and cyan) approached the surface from different angles and converged on the surface with different patterns, thus producing distinct ion channel correlations, consistent with O’Leary et al., 2013. (B) Same as panel A (n_in = 3), but for homeostatic regulation of firing rate and energy efficiency (n_out = 2). Solutions converge on a curve (yellow), giving rise to virtually identical ion channel correlations regardless of regulation rates.

However, correlations produced through the homeostatic regulation mechanism described by O’Leary et al., 2013, which is essentially the same mechanism used here, were recently shown to be sensitive to noise (Franci et al., 2020). This occurs because this regulation mechanism brings conductance densities to the solution manifold (by minimizing the error signal) but cannot control how solutions drift on the manifold (since the error signal is 0 everywhere on the manifold). Solutions do not drift in the absence of noise but, of course, noise is ubiquitous in biological systems and its effects must be considered. Accordingly, we reran simulations from Figures 6 and 7 with noise added to the conductance densities. After a few hundred noise-update iterations, solutions spread across the available solution space, which corresponds to a surface when only firing rate is regulated (Figure 8A, top). The distribution of solutions can produce correlations (Figure 8A, bottom). Lower and upper bounds on conductance densities shape the solution space and influence correlations (Figure 8—figure supplement 1A). Correlations also depend on how solutions distribute across the solution space, which depends on regulation rates (Figure 8—figure supplement 1B, C), which means correlations can still depend on regulation rates under noisy conditions, but correlations are not the same as in the absence of noise (Figure 8B). When firing rate and energy efficiency were regulated, noise caused solutions to spread across the intersection (Figure 8C) but correlations were relatively unaffected since they are limited by the low dimensionality of this solution space. Spread was bounded by one or another conductance density reaching 0 mS/cm². Franci et al., 2020 proposed another regulatory scheme in which an attractive subspace emerges through cooperative molecular interactions; correlations created by fluctuations along that subspace are robust to noise. That scheme should preclude correlations from arising as in Figure 8A because solutions are prevented from spreading across the solution space; conversely, a low-dimensional solution space like in Figure 8C may prevent correlations unless the attractive subspace arising from molecular interactions aligns with the solution space. Further work is required to explore these regulatory mechanisms and how they interact.

Figure 8 with 1 supplement see all

Download asset Open asset

Noise can affect ion channel correlations depending on the dimensionality of solution space.

(A) Starting from solution sets in Figure 6B, ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{M}$ had noise added to them and were then updated by homeostatic regulation (using regulation rates from Figure 6) to correct the noise-induced disruption of firing rate. Conductance density combinations are depicted before and after 50, 200, 1000, or 3000 noise-update iterations (left to right). Noise caused solutions to spread across the surface (top). Some solutions (number in italics) drifted beyond the illustrated region but this is prevented by imposing an upper bound (Figure 8—figure supplement 1A). Ion channel correlations reflect the available solution space; for example, changes in ${\bar{g}}_{N a}$ and ${\bar{g}}_{K}$ are limited by ${\bar{g}}_{M}$ remaining positive (dashed green line), while changes in ${\bar{g}}_{N a}$ and ${\bar{g}}_{M}$ are limited by ${\bar{g}}_{K}$ remaining positive (dashed purple line). How solutions distribute within that space depends on regulation rates, and also affects correlations (Figure 8—figure supplement 1B, C). Conductance density distributions and pairwise correlations (bottom) were not centered on the mean or normalized by standard deviation (unlike in other figures) in order to visualize how distributions evolve over iterations. (B) Comparison of regression lines from Figure 6 (gray) and 7 (pink, cyan) without noise (top) and with noise (bottom). Correlations are affected by regulation rates in both conditions, but not in the same way. (C) Same as panel A but for regulation of firing rate *and* energy efficiency. Solutions spread along the intersection of the two surfaces, but the spread is bounded by one or another conductance density reaching 0 mS/cm². Correlations slightly increase under noisy conditions (compare with Figure 6).

The dimensionality of solution space affects the success of homeostatic regulation

Beyond affecting the ion channel correlations that emerge through homeostatic regulation, we predicted that the dimensionality of solution space affects whether multiple properties can be successfully regulated to their target values. Figure 9A shows an example in which coadjusting ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{M}$ successfully regulates firing rate. Figure 9B shows successful regulation of energy efficiency by coadjusting the same three channels. Using the same initial conditions and relative regulation rates as above, the system failed to coregulate firing rate and energy efficiency (Figure 9C). Notably, coordinated regulation of both properties was achieved using other relative regulation rates (see Figures 6C and 7B) or using the same relative rates but starting from different initial conditions (not shown), suggesting that low-dimensional solutions are less accessible (i.e., a smaller set of relative regulation rates succeed in finding the solution space). This may be a challenge for the Franci et al. regulation model if molecular interactions create an attractive subspace that does not align with a low-dimensional solution space.

Figure 9

Download asset Open asset

Low-dimensional solutions can be hard for homeostatic regulation to ‘find’.

(A) Homeostatic regulation of firing rate (n_out = 1) via control of ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{M}$ (n_in = 3), like in Figures 6B and 7A but using a different set of regulation rates (see Supplementary file 1). Solutions converged onto the iso-firing rate surface (top panel) and firing rate was regulated to its target value in <30 iterations (bottom panel). (B) Same as A (n_in = 3) but for homeostatic regulation of energy efficiency (n_out = 1). Solutions converged on the iso-energy efficiency surface (top panel) and energy efficiency was regulated to its target value in ~10 iterations (bottom panel). (C) Homeostatic regulation of firing rate and energy efficiency (n_out = 2) via control of the same ion channels (n_in = 3) using the same relative rates and initial values as in panels A and B. Neither firing rate (red trajectories) nor energy efficiency (green trajectories) reached its target value. Conductance densities were capped at 4 mS/cm² but this does not account for trajectories not reaching their target. Noise was not included in these simulations.

Figure 10 shows additional examples of regulating firing rate and energy efficiency by coadjusting ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{M}$ . In these simulations, firing rate is regulated to a precise target value whereas energy efficiency is maintained above a lower bound; the single-output solution set for energy efficiency thus corresponds to a volume rather than a surface for n_in = 3. For energy efficiency ≥22%, homeostatically determined solutions converge on the iso-firing rate surface without regulation of energy efficiency having much effect (Figure 10A). For energy efficiency ≥27%, solutions initially converge onto part of the iso-firing rate surface that sits outside the targeted energy efficiency volume, but without energy efficiency requirements being met, solutions then move across the surface until reaching the intersection with the iso-energy efficiency volume (Figure 10B). By converging on the intersection, which is a curve, ion channel correlations are stronger than in Figure 10A. In other words, ion channel correlations increase as single-output solution sets start to disconnect, meaning increased correlations presage regulation failure; in general, this is consistent with the impact of increasing constraints (see Figure 8), which go hand in hand with decreasing degrees of freedom. For energy efficiency ≥30%, the iso-firing rate surface and targeted energy efficiency volume do not intersect and solutions thus proceed to a point between the two single-output solution sets constrained by ${\bar{g}}_{K}$ and ${\bar{g}}_{M}$ reaching 0 mS/cm² (Figure 10C). In this last example, neither property is optimally regulated but the outcome is a reasonable compromise. Regulation might have failed outright if, in the absence of an upper bound, ion channel densities increased (wound-up) without properties ever reaching their target values. This highlights that coregulation of multiple properties might fail or find suboptimal solutions not because single-output solutions do not exist, but because even large single-output solution sets might not intersect, meaning no multioutput solution exists.

Figure 10

Download asset Open asset

The outcome of homeostatic regulation depends on if and how single-output solution sets intersect.

Homeostatic regulation of firing rate and energy efficiency (n_out = 2) via control of ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{M}$ (n_in = 3). For these simulations, energy efficiency was maintained above a lower bound rather than being regulated to a specific target value; accordingly, the single-output solution set for energy efficiency corresponds to a volume (green) rather than a surface. (A) For energy efficiency ≥22%, homeostatically determined solutions converge on the iso-firing rate surface (red) in a region sitting within the green volume (top panel). The rate of convergence and resulting ion channel correlations are shown in the middle and bottom panels, respectively. (B) For energy efficiency ≥27%, solutions initially converge on the red surface in a region outside the green volume, but trajectories then bend and proceed across the red surface until that surface reaches the green volume. Because solutions converge on a curve, ion channel correlations are stronger than in A, where solutions distributed across a surface. (C) For energy efficiency ≥30%, the red surface and green volume do not intersect. Consequently, solutions settle between the two single-output solution sets (top panel) without either property reaching its target (middle panel). The outcome represents the balance achieved by the opposing pull of control mechanisms regulating different properties, and depends entirely on ${\bar{g}}_{N a}$ since ${\bar{g}}_{K}$ and ${\bar{g}}_{M}$ cannot become negative (bottom panel). Noise was not included in these simulations.

Discussion

This study has identified conditions required for a neuron to coregulate >1 property. Being able to produce the same output using diverse ion channel combinations ensures that single-output solution sets are large. This is critical because of the many ion channel combinations that produce the desired output for one property, few also produce the desired output for other properties (Figure 3), which means ion channel adjustments that regulate one property are liable to disrupt other properties due to ion channel pleiotropy (Figures 2 and 4). Indeed, the multioutput solution set corresponds to the intersection between single-output solution sets and the dimensionality of the multioutput solution corresponds to the difference between the number of adjustable ion channels (n_in) and the number of regulated outputs (n_out) (Figure 5). Coregulation of n properties requires at least n adjustable ion channels for a unique solution, and at least n + 1 channels for a degenerate solution. This constraint is not alleviated by pleiotropy, but n_in would need to exceed n_out by an even wider margin if ion channels were not pleiotropic. Moreover, channels must be coadjusted with different ratios to regulate different properties, which constrains how feedback loops are organized (see below). These important issues get overlooked if regulation of each property is considered in isolation.

With respect to the number of channels required to coregulate a certain number of properties, a direct analogy can be made with a system of linear equations. Each unknown constitutes a degree of freedom and each equation constitutes a constraint that reduces the degrees of freedom by one. The system is said to be overdetermined if equations outnumber unknowns, and underdetermined if unknowns outnumber equations. An underdetermined system can have infinite solutions. Degeneracy is synonymous with underdetermination, which highlights that degeneracy depends not only on ion channel diversity, but also on those channels not having to coregulate ‘too many’ properties. Extra degrees of freedom broaden the range of available solutions, meaning a good solution is liable to exist over a broader range of conditions. On the other hand, if constraints outnumber the degrees of freedom, the system becomes overdetermined and solutions disappear (e.g., Figure 5D), which can cause regulation to fail (see below). One might reasonably speculate that ion channel diversity has been selected for because it enables the addition of new functionality (e.g., excitability) without the cell becoming overdetermined, lest pre-existing functions (e.g., osmoregulation) become compromised.

There are notable similarities and differences between a neuron adjusting its ion channel densities to regulate properties and a scientist trying to infer conductance densities based on the measured values of those properties (i.e., fitting a model to experimental data). Fitting a model with n parameters to many outputs (e.g., firing rate and input resistance and rheobase and spike height) is more difficult but yields better parameter estimates than fitting the same model to just one output (Foster et al., 1993), just as regulating more properties makes the solution set smaller. But how good are those parameter estimates? Assuming there is no measurement noise, can one confidently infer the true conductance densities in a particular neuron by measuring and fitting enough properties? Degeneracy makes solving this inverse problem difficult, if not impossible (Sarkar and Sobie, 2010; Aster et al., 2013). The neuron solves the forward problem, producing a firing rate, input resistance, etc. based on its channel density combination. That combination is determined by negative feedback but the precise densities are unimportant so long as the combination produces the target values for all regulated properties (which happens for every density combination in the multioutput solution set). As conditions change, the multioutput solution will evolve and negative feedback will adjust densities accordingly. A neuron needs to regulate its properties but does not do so by regulating its conductance densities to particular values; in that respect, a neuron does not solve an inverse problem. (Convergence of conductance densities to the same ‘attractive’ set regardless of initial conditions, or after a perturbation, may suggest otherwise [Franci et al., 2020], but may reflect other, unaccounted for constraints like regulation of another property.) Neurons have evolved under selective pressure to be degenerate (see above), not parsimonious, contrary to how most models are constructed.

Previous studies using grid searches to explore degeneracy have tended to apply selection criteria simultaneously, finding the ‘multioutput’ solution set in one fell swoop. In contrast, we considered one criterion (property) at a time in order to find each single-output solution set, which we then combined to find the multioutput solution. The former approach is akin to aggregating several error functions to create a single-objective problem, whereas the latter resembles multiobjective optimization. Nevertheless, past studies have observed that certain parameter changes (in a particular direction through parameter space) dramatically affect some properties but not others (Drion et al., 2015; Goldman et al., 2001), consistent with manifolds representing the single-output solution sets for sensitive and insensitive properties lying orthogonal to one another in parameter space. The effects of solution dimensionality on ion channel correlations and regulation failure (see below) highlight the value of a multiobjective perspective, which implicitly recognizes that there are multiple properties, each with its own feedback loop, even if those feedback loops intersect. Our results show that coregulating many properties using a common set of ion channels is feasible only if enough different channels are involved. Indeed, for coadjustments to work, channels must differ from each other in how they affect each property (see Figure 1), which emphasizes the distinction between degeneracy and redundancy.

Our results also highlight the need for >1 ‘master regulator’ because regulation of each property requires that ion channels are coadjusted with different ratios (Figure 11). O’Leary et al., 2014 argued in favor of a single master regulator, but that was predicated on using a single factor, namely global calcium, to encode the error signal; because a single factor cannot reach two different targets (i.e., simultaneously minimize two errors), it cannot support two master regulators. But local (Bootman et al., 2001) or kinetically distinct (Liu et al., 1998) variations in calcium could encode >1 error signal if the calcium sensors coupling intracellular calcium to gene expression or channel modulation are spatially segregated or operate with different filter properties. Beyond calcium, evidence points to firing rate homeostasis being activity independent (MacLean et al., 2003), using membrane potential as feedback (Santin and Schulz, 2019), or involving dual feedback loops (Kulik et al., 2019). Furthermore, the AMP:ATP ratio is used to track energy level (Hardie, 2014; Garcia and Shaw, 2017) and free amino acids can also be sensed (Efeyan et al., 2012). This is just the tip of the iceberg. Rather than error signals being encoded by a single factor like calcium, we favor the notion of a multi-input/multioutput system (Papin and Palsson, 2004) in which factors combine to encode multiple error signals. Take mTOR (mechanistic target of rapamycin) for example; it is modulated by diverse factors and, in turn, modulates many processes including transcription, translation and protein degradation (Switon et al., 2017). The regulation model proposed by Franci et al., 2020 provides additional insights. These issues require more investigation but, for now, our results argue that coregulation of multiple properties is incompatible with the amalgamation of error signals or certain control steps. The analogy with single- and multiobjective optimization (see above) seems apropos. Critically, n_in exceeding n_out will not guarantee solutions if degrees of freedom are reduced by bottlenecks elsewhere in the feedback loops.

Figure 11

Download asset Open asset

Each error signal must be able to coadjust channels with different ratios.

Cartoons depict regulation of two properties (A and B) via control of three ion channels (X, Y, and Z). Channel levels are homeostatically adjusted to minimize the difference (error) between each property and its target value. The error signal divided by each channel’s regulation time constant (τ) dictates the magnitude of the change in that channel (represented by length of curved arrows in panel B). (A) The set of regulation time constants define a ‘master regulator’ (black dial) that coadjusts ion channels according to a certain ratio. Convergence of error signals at or before the master regulator limits how ion channels are coadjusted. (B) An additional master regulator (red dial) is needed to enable *error B* to coadjust channels with a different ratio than *error A*. Conditions like in panel B were required for all coregulation simulations (see Figure 6C). Different coadjustment ratios are required to regulate each property because of the different way each ion channel affects each property (see Figure 1). Note that master regulators may be more distributed than depicted here, with each regulation time constant *τ_i,j* likely reflecting the net effect of regulating multiple (transcriptional, translational, etc.) processes. The positive molecular regulatory network proposed by Franci et al., 2020 involves additional feedback mechanisms not easily depicted in this sort of cartoon.

Ion channel correlations have been studied in simulations (O’Leary et al., 2014; Mukunda and Narayanan, 2017; Taylor et al., 2009; Jain and Narayanan, 2020; O’Leary et al., 2013; Franci et al., 2020; Soofi et al., 2012; Hudson and Prinz, 2010; Ball et al., 2010) and experiments (Zhao and Golowasch, 2012; Temporal et al., 2014; Schulz et al., 2006; Tobin et al., 2009; Schulz et al., 2007; Khorkova and Golowasch, 2007). These correlations have been ascribed to the coregulation of ion channels. The relative rates with which different conductance densities are controlled dictate the direction in which trajectories move through parameter space, which in turn dictates how trajectories approach and distribute across the solution manifold (O’Leary et al., 2013). Our results are consistent with that explanation (Figure 7), notwithstanding noise effects (see below), but they also highlight the importance of manifold dimensionality: Correlations are stronger but are insensitive to relative regulation rates if trajectories reach a 1D manifold (curve) rather than a higher-dimensional manifold (surface, volume, etc.) (Figure 6). Under noisy conditions, correlations reflect the available solution space and how solutions drift across it, which depends on regulation rates (Figure 8). The cooperative molecular interactions proposed by Franci et al., 2020 can produce correlations in a high-dimensional solution space, but may be prevented from doing so by a low-dimensional solution space. Notwithstanding molecular interactions, the existence of pairwise correlations suggests that the dimensionality of solution space (=n_in − n_out) is relatively low. This may be surprising since n_in is high, but makes sense if n_out is also high; in other words, there are many channels but they are responsible for regulating many properties. Though n_out is constrained by n_in (i.e., the number of regulated properties cannot exceed the number of adjustable ion channels), neurons might take advantage of regulating as many properties as their ion channel diversity safely allows, leading to a relatively low-dimensional solution space (as n_out approaches n_in).

Homeostatic regulation can fail for different reasons. If there are no solutions (i.e., the solution set is empty), negative feedback cannot regulate a property to its target values. A multioutput solution set can be empty because single-output solution sets do not intersect (Figure 10C), meaning regulation of different properties to their respective target values is incompatible. Regulation can also fail because negative feedback fails to converge on available solutions (Figure 9C). Notably, a multioutput solution set may become less accessible (lower dimensional) as single-output solution sets start to separate (Figure 10B), and may foreshadow the eventual disjunction of those solution sets (Figure 10C). Regulation can also fail because feedback signaling is compromised. These failure modes are not mutually exclusive; for example, a system might reach less accessible solutions if regulation rates are normally flexible but might fail to reach those solutions if regulation rate flexibility is reduced. Cooperative molecular interactions (Franci et al., 2020) are notable in this regard, insofar as they would tend to constrain regulation and thus influence failure modes. Failed homeostatic regulation may have similar consequence regardless of how the failure occurs, consistent with the emergence of common disease phenotypes despite vastly different underlying pathologies (Ramocki and Zoghbi, 2008), but subtle differences in exactly how the failure transpires may provide important clues to help pinpoint the underlying mechanism(s).

In conclusion, neurons can coregulate multiple properties by coadjusting ion channels. Despite reusing channels to regulate more than one property, many channels are required to coregulate many properties and feedback loops must allow for coadjustments with different ratios. This may account for why evolution has yielded such diverse channels and why their transcriptional regulation has become so complicated (Gabel et al., 2015). Given how difficult it is to study regulation of any one property, studying the coregulation of multiple properties seems truly daunting, especially if feedback loops intermingle and error signals involve combinatorial codes. But if that is indeed the case, then accounting for the coregulation of multiple properties may be crucial for appreciating design features that might otherwise elude us.

Methods

Slice electrophysiology

All procedures were approved by the Hospital for Sick Children Animal Care Committee. Using the same procedures and equipment previously described (Khubieh et al., 2016), coronal slices of hippocampus were prepared from an adult mouse and a CA1 pyramidal neuron was recorded using whole-cell patch clamp. A junction potential correction of −9 mV was applied to the recorded membrane potential. A noisy stimulus was generated through an Ornstein–Uhlenbeck process (see Equation 7) with τ_stim = 5 ms, μ_stim = 60 pA, and σ_stim = 10 pA; the same stimulus was replayed on each trial. Native I_M was blocked by bath application of 10 µM XE991 (Tocris). The block was continued throughout dynamic clamp experiments. Virtual I_M and I_AHP were modeled as per Equation 6 (with τ_z = 100 ms, β_z = –35 mV, and γ_z = 4 mV for I_M, and τ_z = 600 ms, β_z = 0 mV, and γ_z = 1 mV for I_AHP) and were applied using the dynamic clamp capabilities of Signal v6 (Cambridge Electronic Design). The density of each virtual conductance was adjusted to produce the desired firing rates ( ${\bar{g}}_{M}$ = 10 nS, ${\bar{g}}_{A H P}$ = 120–150 nS).

Neuron model

Our base model includes a fast-activating sodium conductance ( $g_{f a s t}$ ) and a slower-activating potassium conductance ( $g_{s l o w}$ ) which together are sufficient to produce spikes. Other conductances were added to modulate excitability.

C \frac{d V}{d t} = I_{s t i m} - {\bar{g}}_{f a s t} m_{\infty} (V) (V - E_{N a}) - {\bar{g}}_{s l o w} w (V - E_{K}) - {\bar{g}}_{N a} n (V - E_{N a}) - {\bar{g}}_{K} n (V - E_{K}) - {\bar{g}}_{a H P} p (V - E_{K}) - {\bar{g}}_{M} q (V - E_{K}) - {\bar{g}}_{l e a k} (V - E_{l}),

where V is voltage and m changes instantaneously with V whereas other gating variables change more slowly. The spike-generating conductances $g_{f a s t}$ and $g_{s l o w}$ were modeled using a Morris–Lecar formalism (Rho and Prescott, 2012), with C = 2 µF/cm², E_Na = 50 mV, E_K= −100 mV, E_leak = −70 mV, φ_w= 0.15, ${\bar{g}}_{f a s t}$ = 20 mS/cm², ${\bar{g}}_{s l o w}$ = 20 mS/cm², ${\bar{g}}_{l e a k}$ = 2 mS/cm², β_m = −1.2 mV, γ_m = 18 mV, β_w = −10 mV, and γ_w = 10 mV. A generic sodium conductance ( $g_{N a}$ ) and potassium conductance ( $g_{K}$ ) were modeled using Hodgkin–Huxley formalism as described by Ratté et al., 2014b

I_{N a, K} = {\bar{g}}_{N a, K} n (V - E_{N a, K}),

\frac{d n}{d t} = α (1 - n) - β n,

α = \frac{k_{α} (V - V_{α})}{e^{(V - \frac{V_{α}}{s_{α}})} - 1},

β = k_{β} e^{(V - \frac{V_{α}}{s_{α}})},

where V_α,β = −24 mV, s_α,β = −17 mV, and k_α,β = 1 ms⁻¹. Note that $g_{N a}$ and $g_{K}$ differ only in their reversal potentials. Channels approximating a calcium-activated potassium conductance ( $g_{A H P}$ ) and an M-type potassium conductance ( $g_{M}$ ) were modeled as described by Prescott and Sejnowski, 2008

\frac{d z}{d t} = {\frac{1}{1 + e^{(\frac{β_{z} - V}{γ_{z}})}} - z} / τ_{z},

where τ_z= 100 ms, γ_z = 4 mV, and β_z = 0 mV or −35 mV for $g_{A H P}$ and $g_{M}$ , respectively. Maximal conductance densities for $g_{l e a k}$ , $g_{N a}$ , $g_{K}$ , and $g_{M}$ were systematically varied or adjusted by a homeostatic feedback mechanism (see below). Any other parameters that differ from the sources cited above are reported in the relevant figure legends. Injected current I_stim was applied as either a constant step or as noisy fluctuations modeled with an Ornstein–Uhlenbeck process

\frac{d I_{s t i m}}{d t} = \frac{I_{s t i m} - μ_{s t i m}}{τ_{s t i m}} + \frac{S σ_{s t i m}}{\sqrt{d t}} N (t),

where τ_stim is a time constant that controls the rate at which I_stim drifts back toward the mean μ_stim, and N (t) is a random number drawn from a normal distribution with 0 mean and unit variance that is scaled by Sσ_stim, where a scaling factor S = √2/τ_stim makes the standard deviation σ_stim independent of τ_stim. All simulations were conducted in MATLAB using the forward Euler integration method and a time step of 0.05–0.1 ms.

Energy calculations

Energy consumption rate was calculated as in Hasenstaub et al., 2010 Briefly, models were stimulated with a fast fluctuating stimulus (µ_stim = 40 µA/cm², σ_stim = 10 μA/cm²). Sodium and potassium current through all channels was integrated for 1 s to determine the charge for each ion species, which was then converted to ion flux based on the elementary charge 1.602 × 10⁻¹⁹ C. Based on the 3:2 stoichiometry of the Na⁺/K⁺ pump, we divided the number of sodium and potassium ions by 3 and 2, respectively, and used the maximum of those two values as the energy consumption rate (in ATP/cm² s). Energy efficiency was calculated as the ratio between capacitive minimum and total Na⁺ flux during an action potential (Sengupta et al., 2010). Briefly, voltage was reset to −40 mV to evoke a single spike in all models. Models were treated as pure capacitors to calculate the minimum capacitive current as CΔV, where C is the capacitance and ΔV is the difference between the resting membrane potential and the spike peak.

Grid search

Models were tested with conductance density combinations chosen from a 100 × 100 or 30 × 30 × 30 grid for 2D and 3D plots, respectively. To depict each single-output solution set, all models with outputs within a range (tolerance) of the target value (±3 spk/s for firing rate, ±0.25% for energy efficiency, and ±0.003 kΩ cm² for input resistance) were selected and a curve or surface was fit to those successful models. The same tolerances were implemented in the homeostatic learning rule (see below) to minimize ringing. Tolerances are illustrated in Figure 5—figure supplement 1 but were not shown in other figures for sake of clarity.

Feedback control

We used a homeostatic learning rule similar to O’Leary et al., 2014; O’Leary et al., 2013 with two notable differences: (1) we did not use intracellular calcium or other biological signals (e.g., AMP:ATP ratio) as intermediaries for our error signals and (2) the error for each output was determined as the difference between the current and target values at the end of each iteration, and conductance densities were adjusted before the start of the next iteration, rather than updating and feeding back error signals in ‘real time’. The separation between fast and slow (regulatory) timescales justifies the former approach. For each conductance ${\bar{g}}_{i}$ , error was divided by the regulation time constant (τ_i) and added to the conductance (see Figure 4A). For coregulating >1 property, we used the sum of scaled errors to update conductance densities (see Figure 6D). A single run consists of a maximum of 200 iterations, during which a model must reach and maintain its regulated property within the tolerance for five consecutive iterations. All models that reached the target output(s) did so in well under 100 iterations; regulation was deemed to have failed for models not reaching their target output(s) within 200 iterations. Conductance densties during the last five iterations were averaged and reported as the final value. See Supplementary file 1 for the initial conductance densities and regulation time constants used for each figure.

Conductance noise

For simulations in Figure 8, noisy variations in each conductance density, ${\bar{g}}_{N a}$ , ${\bar{g}}_{K}$ , and ${\bar{g}}_{M}$ , was applied by adding a random number drawn from a Gaussian distribution with a mean of 0 and a standard deviation of 0.05 mS/cm². Noise was independent for each conductance. After applying noise, the simulation was run, errors were calculated, and conductance densities were updated according to the feedback control described above; this qualifies as one noise-update iteration. Up to 3000 noise-update iterations were run. If the addition of noise caused a conductance density to become negative, the density was reset to 0 mS/cm² before applying feedback control. Where indicated, conductance density was likewise reset to 4 mS/cm² if noise caused it to increase above that value.

Code availability

All computer code is available at http://modeldb.yale.edu/267309 and at http://prescottlab.ca/code-for-models.

Data availability

All computer code is available at http://modeldb.yale.edu/267309 and at http://prescottlab.ca/code-for-models. Key parameter values are provided in Supplementary file 1. Other parameter values are identified in the Methods. Source data are provided for Figure 2.

The following data sets were generated

1. Yang J
(2022) ModelDB
ID 267309. A modified Morris-Lecar model with gM and gAHP.

http://modeldb.yale.edu/267309

References

1. Achard P
2. De Schutter E
(2006) Complex parameter landscape for a complex neuron model
PLOS Computational Biology 2:e94.

https://doi.org/10.1371/journal.pcbi.0020094
- PubMed
- Google Scholar
(2003) Visually driven regulation of intrinsic neuronal excitability improves stimulus detection in vivo
Neuron 39:831–842.

https://doi.org/10.1016/s0896-6273(03)00527-0
- PubMed
- Google Scholar
1. Al-Basha D
2. Prescott SA
(2019) Intermittent failure of spike propagation in primary afferent neurons during tactile stimulation
The Journal of Neuroscience 39:9927–9939.

https://doi.org/10.1523/JNEUROSCI.0975-19.2019
- PubMed
- Google Scholar
1. Anirudhan A
2. Narayanan R
(2015) Analogous synaptic plasticity profiles emerge from disparate channel combinations
The Journal of Neuroscience 35:4691–4705.

https://doi.org/10.1523/JNEUROSCI.4223-14.2015
- PubMed
- Google Scholar
Book
(2013) Parameter Estimation and Inverse Problems
Waltham, MA: Academic Press.

https://doi.org/10.1016/C2009-0-61134-X
- Google Scholar
1. Attwell D
2. Laughlin SB
(2001) An energy budget for signaling in the grey matter of the brain
Journal of Cerebral Blood Flow and Metabolism 21:1133–1145.

https://doi.org/10.1097/00004647-200110000-00001
- PubMed
- Google Scholar
(2008) Adapting proteostasis for disease intervention
Science (New York, N.Y.) 319:916–919.

https://doi.org/10.1126/science.1141448
- PubMed
- Google Scholar
1. Ball JM
2. Franklin CC
3. Tobin AE
4. Schulz DJ
5. Nair SS
(2010) Coregulation of ion channel conductances preserves output in a computational model of a crustacean cardiac motor neuron
The Journal of Neuroscience 30:8637–8649.

https://doi.org/10.1523/JNEUROSCI.6435-09.2010
- PubMed
- Google Scholar
1. Bhalla US
2. Bower JM
(1993) Exploring parameter space in detailed single neuron models: simulations of the mitral and granule cells of the olfactory bulb
Journal of Neurophysiology 69:1948–1965.

https://doi.org/10.1152/jn.1993.69.6.1948
- PubMed
- Google Scholar
(2001) The organisation and functions of local Ca²⁺ signals
Journal of Cell Science 114:2213–2222.

https://doi.org/10.1242/jcs.114.12.2213
- PubMed
- Google Scholar
(2016) Consequences of degeneracy in network function
Current Opinion in Neurobiology 41:62–67.

https://doi.org/10.1016/j.conb.2016.07.008
- PubMed
- Google Scholar
1. Davis GW
(2006) Homeostatic control of neural activity: from phenomenology to molecular design
Annual Review of Neuroscience 29:307–323.

https://doi.org/10.1146/annurev.neuro.28.061604.135751
- PubMed
- Google Scholar
(1999) Plasticity in the intrinsic excitability of cortical pyramidal neurons
Nature Neuroscience 2:515–520.

https://doi.org/10.1038/9165
- PubMed
- Google Scholar
(2015) Ion channel degeneracy enables robust and tunable neuronal firing rates
PNAS 112:5361–5370.

https://doi.org/10.1073/pnas.1516400112
- PubMed
- Google Scholar
1. Edelman GM
2. Gally JA
(2001) Degeneracy and complexity in biological systems
PNAS 98:13763–13768.

https://doi.org/10.1073/pnas.231499798
- PubMed
- Google Scholar
(2012) Amino acids and mTORC1: from lysosomes to disease
Trends in Molecular Medicine 18:524–533.

https://doi.org/10.1016/j.molmed.2012.05.007
- PubMed
- Google Scholar
(2010) Receptor regulation of osmolyte homeostasis in neural cells
The Journal of Physiology 588:3355–3364.

https://doi.org/10.1113/jphysiol.2010.190777
- PubMed
- Google Scholar
(1993) Significance of conductances in Hodgkin-Huxley models
Journal of Neurophysiology 70:2502–2518.

https://doi.org/10.1152/jn.1993.70.6.2502
- PubMed
- Google Scholar
(2020) Positive dynamical networks in neuronal regulation: How tunable variability coexists with robustness
IEEE Control Systems Letters 4:946–951.

https://doi.org/10.1109/LCSYS.2020.2997214
- Google Scholar
1. Frere S
2. Slutsky I
(2018) Alzheimer’s disease: from firing instability to homeostasis network collapse
Neuron 97:32–58.

https://doi.org/10.1016/j.neuron.2017.11.028
- PubMed
- Google Scholar
1. Gabel HW
2. Kinde B
3. Stroud H
4. Gilbert CS
5. Harmin DA
6. Kastan NR
7. Hemberg M
8. Ebert DH
9. Greenberg ME
(2015) Disruption of DNA-methylation-dependent long gene repression in Rett syndrome
Nature 522:89–93.

https://doi.org/10.1038/nature14319
- PubMed
- Google Scholar
1. Garcia D
2. Shaw RJ
(2017) AMPK: mechanisms of cellular energy sensing and restoration of metabolic balance
Molecular Cell 66:789–800.

https://doi.org/10.1016/j.molcel.2017.05.032
- PubMed
- Google Scholar
1. Goaillard JM
2. Dufour MA
(2014) The pros and cons of degeneracy
eLife 3:e02615.

https://doi.org/10.7554/eLife.02615
- PubMed
- Google Scholar
1. Goaillard JM
2. Marder E
(2021) Ion channel degeneracy, variability, and covariation in neuron and circuit resilience
Annual Review of Neuroscience 44:335–357.

https://doi.org/10.1146/annurev-neuro-092920-121538
- PubMed
- Google Scholar
(2001)
Global structure, robustness, and modulation of neuronal models

The Journal of Neuroscience 21:5229–5238.
- PubMed
- Google Scholar
(2010) Compensation for variable intrinsic neuronal excitability by circuit-synaptic interactions
The Journal of Neuroscience 30:9145–9156.

https://doi.org/10.1523/JNEUROSCI.0980-10.2010
- PubMed
- Google Scholar
(2008) Channel density distributions explain spiking variability in the globus pallidus: a combined physiology and computer simulation database approach
The Journal of Neuroscience 28:7476–7491.

https://doi.org/10.1523/JNEUROSCI.4198-07.2008
- PubMed
- Google Scholar
1. Hardie DG
(2014) AMP-activated protein kinase: maintaining energy homeostasis at the cellular and whole-body levels
Annual Review of Nutrition 34:31–55.

https://doi.org/10.1146/annurev-nutr-071812-161148
- PubMed
- Google Scholar
(2010) Metabolic cost as a unifying principle governing neuronal biophysics
PNAS 107:12329–12334.

https://doi.org/10.1073/pnas.0914886107
- PubMed
- Google Scholar
(2013) Firing rate homeostasis in visual cortex of freely behaving rodents
Neuron 80:335–342.

https://doi.org/10.1016/j.neuron.2013.08.038
- PubMed
- Google Scholar
1. Hodgkin AL
(1948) The local electric changes associated with repetitive action in a non-medullated axon
The Journal of Physiology 107:165–181.

https://doi.org/10.1113/jphysiol.1948.sp004260
- PubMed
- Google Scholar
1. Hudson AE
2. Prinz AA
(2010) Conductance ratios and cellular identity
PLOS Computational Biology 6:e1000838.

https://doi.org/10.1371/journal.pcbi.1000838
- PubMed
- Google Scholar
1. Jain A
2. Narayanan R
(2020) Degeneracy in the emergence of spike-triggered average of hippocampal pyramidal neurons
Scientific Reports 10:374.

https://doi.org/10.1038/s41598-019-57243-8
- PubMed
- Google Scholar
1. Joseph A
2. Turrigiano GG
(2017) All for one but not one for all: Excitatory synaptic scaling and intrinsic excitability are coregulated by CaMKIV, whereas inhibitory synaptic scaling is under independent control
The Journal of Neuroscience 37:6778–6785.

https://doi.org/10.1523/JNEUROSCI.0618-17.2017
- PubMed
- Google Scholar
1. Khorkova O
2. Golowasch J
(2007) Neuromodulators, not activity, control coordinated expression of ionic currents
The Journal of Neuroscience 27:8709–8718.

https://doi.org/10.1523/JNEUROSCI.1274-07.2007
- PubMed
- Google Scholar
(2016) Regulation of cortical dynamic range by background synaptic noise and feedforward inhibition
Cerebral Cortex (New York, N.Y) 26:3357–3369.

https://doi.org/10.1093/cercor/bhv157
- PubMed
- Google Scholar
1. Kim EZ
2. Vienne J
3. Rosbash M
4. Griffith LC
(2017) Nonreciprocal homeostatic compensation in Drosophila potassium channel mutants
Journal of Neurophysiology 117:2125–2136.

https://doi.org/10.1152/jn.00002.2017
- PubMed
- Google Scholar
1. Klassen T
2. Davis C
3. Goldman A
4. Burgess D
5. Chen T
6. Wheeler D
7. McPherson J
8. Bourquin T
9. Lewis L
10. Villasana D
11. Morgan M
12. Muzny D
13. Gibbs R
14. Noebels J
(2011) Exome sequencing of ion channel genes reveals complex profiles confounding personal risk assessment in epilepsy
Cell 145:1036–1048.

https://doi.org/10.1016/j.cell.2011.05.025
- PubMed
- Google Scholar
1. Kulik Y
2. Jones R
3. Moughamian AJ
4. Whippen J
5. Davis GW
(2019) Dual separable feedback systems govern firing rate homeostasis
eLife 8:e45717.

https://doi.org/10.7554/eLife.45717
- PubMed
- Google Scholar
(1993) Activity-dependent regulation of conductances in model neurons
Science (New York, N.Y.) 259:1915–1917.

https://doi.org/10.1126/science.8456317
- PubMed
- Google Scholar
1. Liu Z
2. Golowasch J
3. Marder E
4. Abbott LF
(1998) A model neuron with activity-dependent conductances regulated by multiple calcium sensors
The Journal of Neuroscience 18:2309–2320.

https://doi.org/10.1523/JNEUROSCI.18-07-02309.1998
- PubMed
- Google Scholar
(2003) Activity-independent homeostasis in rhythmically active neurons
Neuron 37:109–120.

https://doi.org/10.1016/s0896-6273(02)01104-2
- PubMed
- Google Scholar
1. Maffei A
2. Turrigiano GG
(2008) Multiple modes of network homeostasis in visual cortical layer 2/3
The Journal of Neuroscience 28:4377–4384.

https://doi.org/10.1523/JNEUROSCI.5298-07.2008
- PubMed
- Google Scholar
1. Marder E
2. Goaillard JM
(2006) Variability, compensation and homeostasis in neuron and network function
Nature Reviews Neuroscience 7:563–574.

https://doi.org/10.1038/nrn1949
- PubMed
- Google Scholar
1. Marder E
(2011) Variability, compensation, and modulation in neurons and circuits
PNAS 108:15542–15548.

https://doi.org/10.1073/pnas.1010674108
- PubMed
- Google Scholar
1. Marder E
2. Taylor AL
(2011) Multiple models to capture the variability in biological neurons and networks
Nature Neuroscience 14:133–138.

https://doi.org/10.1038/nn.2735
- PubMed
- Google Scholar
(2015) Hidden in plain view: degeneracy in complex systems
Bio Systems 128:1–8.

https://doi.org/10.1016/j.biosystems.2014.12.003
- PubMed
- Google Scholar
1. Migliore R
2. Lupascu CA
3. Bologna LL
4. Romani A
5. Courcol JD
6. Antonel S
7. Van Geit WAH
8. Thomson AM
9. Mercer A
10. Lange S
11. Falck J
12. Rössert CA
13. Shi Y
14. Hagens O
15. Pezzoli M
16. Freund TF
17. Kali S
18. Muller EB
19. Schürmann F
20. Markram H
21. Migliore M
(2018) The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow
PLOS Computational Biology 14:e1006423.

https://doi.org/10.1371/journal.pcbi.1006423
- PubMed
- Google Scholar
1. Mittal D
2. Narayanan R
(2018) Degeneracy in the robust expression of spectral selectivity, subthreshold oscillations, and intrinsic excitability of entorhinal stellate cells
Journal of Neurophysiology 120:576–600.

https://doi.org/10.1152/jn.00136.2018
- PubMed
- Google Scholar
1. Mukunda CL
2. Narayanan R
(2017) Degeneracy in the regulation of short-term plasticity and synaptic filtering by presynaptic mechanisms
The Journal of Physiology 595:2611–2637.

https://doi.org/10.1113/JP273482
- PubMed
- Google Scholar
1. Olypher AV
2. Calabrese RL
(2007) Using constraints on neuronal activity to reveal compensatory changes in neuronal parameters
Journal of Neurophysiology 98:3749–3758.

https://doi.org/10.1152/jn.00842.2007
- PubMed
- Google Scholar
1. Olypher AV
2. Prinz AA
(2010) Geometry and dynamics of activity-dependent homeostatic regulation in neurons
Journal of Computational Neuroscience 28:361–374.

https://doi.org/10.1007/s10827-010-0213-z
- PubMed
- Google Scholar
1. Onasch S
2. Gjorgjieva J
(2020) Circuit stability to perturbations reveals hidden variability in the balance of intrinsic and synaptic conductances
The Journal of Neuroscience 40:3186–3202.

https://doi.org/10.1523/JNEUROSCI.0985-19.2020
- PubMed
- Google Scholar
(2010) Homeostasis of intrinsic excitability in hippocampal neurones: dynamics and mechanism of the response to chronic depolarization
The Journal of Physiology 588:157–170.

https://doi.org/10.1113/jphysiol.2009.181024
- PubMed
- Google Scholar
1. O’Leary T
2. Wyllie DJA
(2011) Neuronal homeostasis: time for a change
The Journal of Physiology 589:4811–4826.

https://doi.org/10.1113/jphysiol.2011.210179
- PubMed
- Google Scholar
(2013) Correlations in ion channel expression emerge from homeostatic tuning rules
PNAS 110:2645–2654.

https://doi.org/10.1073/pnas.1309966110
- PubMed
- Google Scholar
(2014) Cell types, network homeostasis, and pathological compensation from a biologically plausible ion channel expression model
Neuron 82:809–821.

https://doi.org/10.1016/j.neuron.2014.04.002
- PubMed
- Google Scholar
1. O’Leary T
2. Marder E
(2016) Temperature-robust neural function from activity-dependent ion channel regulation
Current Biology 26:2935–2941.

https://doi.org/10.1016/j.cub.2016.08.061
- PubMed
- Google Scholar
1. O’Leary T
(2018) Homeostasis, failure of homeostasis and degenerate ion channel regulation
Current Opinion in Physiology 2:129–138.

https://doi.org/10.1016/j.cophys.2018.01.006
- Google Scholar
1. Papin JA
2. Palsson BO
(2004) Topological analysis of mass-balanced signaling networks: a framework to obtain network properties including crosstalk
Journal of Theoretical Biology 227:283–297.

https://doi.org/10.1016/j.jtbi.2003.11.016
- PubMed
- Google Scholar
1. Prescott SA
2. Sejnowski TJ
(2008) Spike-rate coding and spike-time coding are affected oppositely by different adaptation mechanisms
The Journal of Neuroscience 28:13649–13661.

https://doi.org/10.1523/JNEUROSCI.1792-08.2008
- PubMed
- Google Scholar
1. Price CJ
2. Friston KJ
(2002) Degeneracy and cognitive anatomy
Trends in Cognitive Sciences 6:416–421.

https://doi.org/10.1016/s1364-6613(02)01976-9
- PubMed
- Google Scholar
(2004) Similar network activity from disparate circuit parameters
Nature Neuroscience 7:1345–1352.

https://doi.org/10.1038/nn1352
- PubMed
- Google Scholar
1. Ramocki MB
2. Zoghbi HY
(2008) Failure of neuronal homeostasis results in common neuropsychiatric phenotypes
Nature 455:912–918.

https://doi.org/10.1038/nature07457
- PubMed
- Google Scholar
1. Rathour RK
2. Narayanan R
(2019) Degeneracy in hippocampal physiology and plasticity
Hippocampus 29:980–1022.

https://doi.org/10.1002/hipo.23139
- PubMed
- Google Scholar
(2014a) Subthreshold membrane currents confer distinct tuning properties that enable neurons to encode the integral or derivative of their input
Frontiers in Cellular Neuroscience 8:452.

https://doi.org/10.3389/fncel.2014.00452
- PubMed
- Google Scholar
1. Ratté S
2. Zhu Y
3. Lee KY
4. Prescott SA
(2014b) Criticality and degeneracy in injury-induced changes in primary afferent excitability and the implications for neuropathic pain
eLife 3:e02370.

https://doi.org/10.7554/eLife.02370
- PubMed
- Google Scholar
1. Ratté S
2. Prescott SA
(2016) Afferent hyperexcitability in neuropathic pain and the inconvenient truth about its degeneracy
Current Opinion in Neurobiology 36:31–37.

https://doi.org/10.1016/j.conb.2015.08.007
- PubMed
- Google Scholar
1. Rho YA
2. Prescott SA
(2012) Identification of molecular pathologies sufficient to cause neuropathic excitability in primary somatosensory afferents using dynamical systems theory
PLOS Computational Biology 8:e1002524.

https://doi.org/10.1371/journal.pcbi.1002524
- PubMed
- Google Scholar
1. Santin JM
2. Schulz DJ
(2019) Membrane voltage is a direct feedback signal that influences correlated ion channel expression in neurons
Current Biology 29:1683–1688.

https://doi.org/10.1016/j.cub.2019.04.008
- PubMed
- Google Scholar
1. Sarkar AX
2. Sobie EA
(2010) Regression analysis for constraining free parameters in electrophysiological models of cardiac cells
PLOS Computational Biology 6:e1000914.

https://doi.org/10.1371/journal.pcbi.1000914
- PubMed
- Google Scholar
(2006) Variable channel expression in identified single and electrically coupled neurons in different animals
Nature Neuroscience 9:356–362.

https://doi.org/10.1038/nn1639
- PubMed
- Google Scholar
(2007) Quantitative expression profiling of identified neurons reveals cell-specific constraints on highly variable levels of gene expression
PNAS 104:13187–13191.

https://doi.org/10.1073/pnas.0705827104
- PubMed
- Google Scholar
1. Schwanhäusser B
2. Busse D
3. Li N
4. Dittmar G
5. Schuchhardt J
6. Wolf J
7. Chen W
8. Selbach M
(2011) Global quantification of mammalian gene expression control
Nature 473:337–342.

https://doi.org/10.1038/nature10098
- PubMed
- Google Scholar
(2010) Action potential energy efficiency varies among neuron types in vertebrates and invertebrates
PLOS Computational Biology 6:e1000840.

https://doi.org/10.1371/journal.pcbi.1000840
- PubMed
- Google Scholar
(2012) Co-variation of ionic conductances supports phase maintenance in stomatogastric neurons
Journal of Computational Neuroscience 33:77–95.

https://doi.org/10.1007/s10827-011-0375-3
- PubMed
- Google Scholar
1. Speakman JR
2. Levitsky DA
3. Allison DB
4. Bray MS
5. de Castro JM
6. Clegg DJ
7. Clapham JC
8. Dulloo AG
9. Gruer L
10. Haw S
11. Hebebrand J
12. Hetherington MM
13. Higgs S
14. Jebb SA
15. Loos RJF
16. Luckman S
17. Luke A
18. Mohammed-Ali V
19. O’Rahilly S
20. Pereira M
21. Perusse L
22. Robinson TN
23. Rolls B
24. Symonds ME
25. Westerterp-Plantenga MS
(2011) Set points, settling points and some alternative models: theoretical options to understand how genes and environments combine to regulate body adiposity
Disease Models & Mechanisms 4:733–745.

https://doi.org/10.1242/dmm.008698
- PubMed
- Google Scholar
1. Stemmler M
2. Koch C
(1999) How voltage-dependent conductances can adapt to maximize the information encoded by neuronal firing rate
Nature Neuroscience 2:521–527.

https://doi.org/10.1038/9173
- PubMed
- Google Scholar
1. Swensen AM
2. Bean BP
(2005) Robustness of burst firing in dissociated purkinje neurons with acute or long-term reductions in sodium conductance
The Journal of Neuroscience 25:3509–3520.

https://doi.org/10.1523/JNEUROSCI.3929-04.2005
- PubMed
- Google Scholar
(2017) Molecular neurobiology of mTOR
Neuroscience 341:112–153.

https://doi.org/10.1016/j.neuroscience.2016.11.017
- PubMed
- Google Scholar
(2009) How multiple conductances determine electrophysiological properties in a multicompartment model
Journal of Neuroscience 29:5573–5586.

https://doi.org/10.1523/JNEUROSCI.4438-08.2009
- PubMed
- Google Scholar
(2014) Activity-dependent feedback regulates correlated ion channel mRNA levels in single identified motor neurons
Current Biology 24:1899–1904.

https://doi.org/10.1016/j.cub.2014.06.067
- PubMed
- Google Scholar
(2009) Correlations in ion channel mRNA in rhythmically active neurons
PLOS ONE 4:e6742.

https://doi.org/10.1371/journal.pone.0006742
- PubMed
- Google Scholar
(1999) Measures of degeneracy and redundancy in biological networks
PNAS 96:3257–3262.

https://doi.org/10.1073/pnas.96.6.3257
- PubMed
- Google Scholar
(2014) Neural and genetic degeneracy underlies Caenorhabditis elegans feeding behavior
Journal of Neurophysiology 112:951–961.

https://doi.org/10.1152/jn.00150.2014
- PubMed
- Google Scholar
(1994) Activity-dependent changes in the intrinsic properties of cultured neurons
Science (New York, N.Y.) 264:974–977.

https://doi.org/10.1126/science.8178157
- PubMed
- Google Scholar
1. Turrigiano G
(2011) Too many cooks? Intrinsic and synaptic homeostatic mechanisms in cortical circuit refinement
Annual Review of Neuroscience 34:89–103.

https://doi.org/10.1146/annurev-neuro-060909-153238
- PubMed
- Google Scholar
(2006) Background activity regulates excitability of rat hippocampal CA1 pyramidal neurons by adaptation of a K⁺ conductance
Journal of Neurophysiology 95:2007–2012.

https://doi.org/10.1152/jn.00220.2005
- PubMed
- Google Scholar
1. Zhao S
2. Golowasch J
(2012) Ionic current correlations underlie the global tuning of large numbers of neuronal activity attributes
The Journal of Neuroscience 32:13380–13388.

https://doi.org/10.1523/JNEUROSCI.6500-11.2012
- PubMed
- Google Scholar

Article and author information

Author details

Jane Yang
1. Neurosciences and Mental Health, The Hospital for Sick Children, Toronto, Canada
2. Institute of Biomedical Engineering, University of Toronto, Toronto, Canada
Contribution
Conceptualization, Data curation, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-0114-5503
Husain Shakil
1. Neurosciences and Mental Health, The Hospital for Sick Children, Toronto, Canada
2. Department of Physiology, University of Toronto, Toronto, Canada
Present address
Division of Neurosurgery, Department of Surgery, University of Toronto, Toronto, Canada

Contribution
Conceptualization, Data curation, Formal analysis, Investigation, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-3995-6811
Stéphanie Ratté

Neurosciences and Mental Health, The Hospital for Sick Children, Toronto, Canada

Contribution
Data curation, Formal analysis, Investigation, Writing – review and editing

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-7005-6726
Steven A Prescott
1. Neurosciences and Mental Health, The Hospital for Sick Children, Toronto, Canada
2. Department of Physiology, University of Toronto, Toronto, Canada
Contribution
Conceptualization, Funding acquisition, Methodology, Project administration, Supervision, Writing – original draft, Writing – review and editing

For correspondence
steve.prescott@sickkids.ca

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-3827-4512

Funding

Canadian Institutes of Health Research (Foundation Grant 167276)

Steven Alec Prescott

Natural Sciences and Engineering Research Council of Canada (Discovery Grant RGPIN 436168)

Steven Alec Prescott

The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.

Acknowledgements

This research was funded by a Discovery Grant from the Natural Sciences and Engineering Research Council (NSERC) and by a Foundation Grant from the Canadian Institutes of Health Research (CIHR). H.S. was supported by Ontario Graduate Scholarship, the University of Toronto Centre for the Study of Pain Scholarship, and James F Crothers Family Scholarship. J.Y. was supported by a SickKids Restracomp Studentship and Ontario Graduate Scholarship. We thank Eve Marder, Etay Hay, Shreejoy Tripathy, and Simon Hardy for constructive feedback on the manuscript.

Ethics

All experimental procedures were approved by The Hospital for Sick Children Animal Care Committee (protocol #53451) and were conducted in accordance with guidelines from the Canadian Council on Animal Care.

Copyright

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

1,218

views
208

downloads
34

citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

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)

Jane Yang
Husain Shakil
Stéphanie Ratté
Steven A Prescott

(2022)

Minimal requirements for a neuron to coregulate many properties and the implications for ion channel correlations and robustness

eLife 11:e72875.

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

Categories and tags

Research organism

Mouse

Share this article

Cite this article

Simultaneous regulation of >1 property in a theoretical neuron.

A compensatory change that restores firing rate disrupts firing pattern in a CA1 pyramidal neuron.

Figure 2—source data 1

Most ion channel combinations producing the target value for one property produce inconsistent values for other properties.

Ion channel changes mediating the same effect on firing rate can oppositely affect energy efficiency.

A degenerate solution for n properties requires at least n + 1 adjustable ion channels.

Dimensionality of solution space affects ion channel correlations.

Effect of relative regulation rates on ion channel correlations depends on the dimensionality of solution space.

Noise can affect ion channel correlations depending on the dimensionality of solution space.

Low-dimensional solutions can be hard for homeostatic regulation to ‘find’.

The outcome of homeostatic regulation depends on if and how single-output solution sets intersect.

Each error signal must be able to coadjust channels with different ratios.

Author details

Jane Yang

Contribution

Competing interests

Husain Shakil

Present address

Contribution

Competing interests

Stéphanie Ratté

Contribution

Competing interests

Steven A Prescott

Contribution

For correspondence

Competing interests

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)

Categories and tags

Research organism

Further reading