Variability in State-Dependent Plasticity of Intrinsic Properties during Cell-Autonomous Self-Regulation of Calcium Homeostasis in Hippocampal Model Neurons1,2,3

Abstract How do neurons reconcile the maintenance of calcium homeostasis with perpetual switches in patterns of afferent activity? Here, we assessed state-dependent evolution of calcium homeostasis in a population of hippocampal pyramidal neuron models, through an adaptation of a recent study on stomatogastric ganglion neurons. Calcium homeostasis was set to emerge through cell-autonomous updates to 12 ionic conductances, responding to different types of synaptically driven afferent activity. We first assessed the impact of theta-frequency inputs on the evolution of ionic conductances toward maintenance of calcium homeostasis. Although calcium homeostasis emerged efficaciously across all models in the population, disparate changes in ionic conductances that mediated this emergence resulted in variable plasticity to several intrinsic properties, also manifesting as significant differences in firing responses across models. Assessing the sensitivity of this form of plasticity, we noted that intrinsic neuronal properties and the firing response were sensitive to the target calcium concentration and to the strength and frequency of afferent activity. Next, we studied the evolution of calcium homeostasis when afferent activity was switched, in different temporal sequences, between two behaviorally distinct types of activity: theta-frequency inputs and sharp-wave ripples riding on largely silent periods. We found that the conductance values, intrinsic properties, and firing response of neurons exhibited differential robustness to an intervening switch in the type of afferent activity. These results unveil critical dissociations between different forms of homeostasis, and call for a systematic evaluation of the impact of state-dependent switches in afferent activity on neuronal intrinsic properties during neural coding and homeostasis.

To address this, we first arrived at a population of 78 experimentally constrained (with 7 different physiological measurements) CA1 pyramidal neuron models involving 12 ion channels, derived from a randomized population of 4000 models built from uniform sampling of 48 different model parameters. Next, we adapted a recent study on cell-autonomous self-evolution of calcium homeostasis in neurons of the crab stomatogastric ganglion  to hippocampal neurons. We ensured that the adapted model included ion channels derived from hippocampal pyramidal neurons, was endowed with detailed calcium-handling mechanisms (including pumps, buffers and the endoplasmic reticulum; Ashhad and Narayanan, 2013) and received different types of afferent activity through AMPA and NMDA receptors. The temporal evolution of messenger RNAs (mRNAs) and conductances corresponding to each of the 12 channels was independently monitored in each of the 78 valid models, with the time courses of mRNA evolution controlled by ionic conductances obtained from the corresponding valid model. Within this modeling framework for cell-autonomous evolution of calcium homeostasis, we tested the impact of switches in afferent activity (between theta oscillations and SWR inputs) on neuronal conductances and intrinsic properties.
Our results suggest that neuronal ion-channel conductances and intrinsic properties could undergo significant plasticity in the process of maintaining calcium homeostasis through a regime of behavioral state-dependent changes in afferent activity. The sign and strength of such intrinsic plasticity were dependent on the specific activity pattern, the temporal sequence of switches, and the specific neuronal model. These results call for a significant reassessment of the impact of state-dependent switches in afferent activity on neuronal and network physiology, especially accounting for potentially adaptive changes in intrinsic properties.

Neuronal model and ion channels
To study state-dependent and cell-autonomous calcium homeostasis in hippocampal CA1 pyramidal neurons, we used a single compartmental cylindrical model of diameter (d) ϭ 100 m and length (L) ϭ 100 m. Passive properties were set as specific membrane resistance (R m ) ϭ 35 k⍀.cm 2 and specific membrane capacitance (C m ) ϭ 1 F/cm 2 . These settings ensured that the passive input resistance (R in ) was ϳ111 M⍀ and the passive membrane time constant was 35 ms (Narayanan and . The neuronal compartment consisted of 11 conductance-based models for ion channels (Fig. 1A) namely, fast sodium (NaF), delayed-rectifier potassium (KDR), A-type potassium (KA), M-type potassium (KM), T-type calcium (CaT), R-type calcium (CaR), N-type calcium (CaN), L-type calcium (CaL), hyperpolarizationactivated cyclic nucleotide gated channel (HCN or h), small conductance (SK) and big conductance calciumactivated potassium (BK) channels. The channel kinetics for NaF, KDR, and KA were obtained from Hoffman et al. (1997) and Migliore et al. (1999), for CaT from Shah et al. (2011), for KM from Migliore et al. (2006), for CaR and CaL from Magee and Johnston (1995) and Poirazi et al. (2003), CaN and SK from Migliore et al. (1995), for HCN from Magee (1998) and Poolos et al. (2002), and for BK from Moczydlowski and Latorre (1983). The reversal potentials for K ϩ and Na ϩ ions were set as -90 and ϩ55 mV, respectively, and for the HCN channel as -30 mV. Accounting for the leak conductance (g leak ϭ1/R m ), this configuration meant the presence of 12 ion channels in our model.

Synaptic receptors
Excitatory synapses containing AMPA and NMDA receptors, modeled using the Goldman-Hodgkin-Katz formulation (Goldman, 1943;Hodgkin and Katz, 1949;Ashhad and Narayanan, 2013) were introduced in the model. Spe- showing the 11 ion channels inserted. The red arrows denote inward currents and the green arrows denote outward currents. B, Voltage responses (top) of the base neuronal model to current pulses (bottom) ranging from -50 to 50 pA in steps of 10 pA. C, The steady state voltages from (B) are plotted against the corresponding current injected. The slope of the resulting V-I plot was defined as the input resistance, R in . D, The voltage response of the base neuronal model to a current injection of 250 pA. The amplitude of the last action potential was defined as the action potential (AP) amplitude, V AP . E, The AP firing frequency (f) versus injected current plot showing the frequency of firing with current injections from 0 to 250 pA in steps of 50 pA. The number of APs elicited by the model in response to a 250 pA, 500 ms current pulse was used to compute the firing rate at 250 pA, f 250 . F, The model's voltage response (top) to a chirp current stimulus of peak-to-peak amplitude of 100 pA, with frequency linearly increasing from 0 to 25 Hz in 25 s (bottom). G, The impedance amplitude profile |Z(f)| derived from traces in F. The frequency at which the impedance amplitude is maximum (|Z| max ) was defined as the resonance frequency, f R . The strength of resonance, Q, was taken as the ratio of |Z(f R )| to |Z(0.5)|. H, The impedance phase profile ((f)) with the area under the inductive part of the curve defined as the total inductive phase (⌽ L ). cifically, a canonical synapse consisting of colocalized NMDA receptor (NMDAR) and AMPA receptor (AMPAR) was modeled as in Narayanan and . The NMDAR current was modeled as a combination of three different types of ionic currents namely Ca 2ϩ , Na ϩ and K ϩ : (4) where P NMDAR defined the maximum permeability of the NMDAR; P Ca ϭ 10.6, P Na ϭ 1, P K ϭ 1 (Mayer and Westbrook, 1987;Canavier, 1999). Extracellular and intracellular concentrations of ions were as follows (in mM): These ionic concentrations set the Na ϩ equilibrium potential at ϩ55 mV and K ϩ equilibrium potential at -90 mV. MgB(v) governed the Mg 2ϩ dependence of the NMDAR current (Jahr and Stevens, 1990): with the default value of [Mg] o set at 2 mM. Current through the AMPAR was modeled as the sum of currents carried by sodium and potassium ions: where, (8) where P AMPAR defined the maximum permeability of the AMPAR. P Na was taken to be equal to P K (Dingledine et al., 1999). The relationship between AMPAR and NMDAR permeabilities was defined as follows: where NAR represented the NMDAR-AMPAR ratio, with its default value set at 1.5.

Calcium dynamics
Calcium handling mechanisms to take care of the calcium reactions, radial diffusion, and buffers were adopted from Ashhad and Narayanan (2013). The following partial differential equation to govern the cytosolic calcium dynamics was used (Sneyd et al., 1995;Fink et al., 2000): where D Ca is the diffusion constant for [Ca 2ϩ ] experimentally determined from Allbritton et al. (1992) and Klingauf and Neher (1997); ␤ is the density of leak channels and SERCA pumps on the endoplasmic reticulum (ER) membrane; J VGCC , J SERCA , R buf , J pump , and J leak are the calcium flux due to voltage-gated calcium channels (VGCCs), sarcoendoplasmic reticulum calcium ATPase (SERCA) pumps, static buffers, membrane pumps, and leak channels respectively. Radial diffusion of calcium was taken care of by compartmentalizing the cylinder into four concentric annuli. The calcium concentration on the outermost annulus was considered as the cytosolic calcium, [Ca 2ϩ ] c (Carnevale and Hines, 2006;Ashhad and Narayanan, 2013). The calcium influx into the cytosol through the ER leak channels was modeled as follows (Fink et al., 2000;Ashhad and Narayanan, 2013): where the leak constant L was chosen such that at resting state (-65 mV), there was no net flux of calcium through the leak channels on the ER membrane. The influx of calcium through the VGCCs and NMDARs in our study (L-, T-, R-, and N-type calcium channels) was modeled as follows (Poirazi et al., 2003;Ashhad and Narayanan, 2013): where I Ca represented the calcium current through the VGCCs/NMDARs, diam is the diameter of the compartment, and F is the Faraday constant. The negative sign indicates the inward nature of I Ca , and accounts for the positive flux of calcium ions with increase in I Ca . The uptake of calcium by the SERCA pump was modeled as follows (Fink et al., 2000;Ashhad and Narayanan, 2013): where V max is the average amplitude of uptake by the pump (1 ϫ 10 Ϫ4 mM/ms) and K p is the dissociation constant of calcium binding to the pump (0.27 M). Ca 2ϩ extrusion through plasma membrane pumps was regulated by a threshold on the cytosolic calcium ([Ca 2ϩ ] c ). The pumps were inactive below a critical Ca 2ϩ concentration, [Ca 2ϩ ] crt , above which the extrusion rate depended linearly on [Ca 2ϩ ] c (Fink et al., 2000): where [Ca 2ϩ ] crt was set at 0.2 M, and ␥ (8 m/s) defines the sensitivity of pump extrusion (Herrington et al., 1996;Fink et al., 2000;Ashhad and Narayanan, 2013). The rate of change in calcium due to the stationary buffers was modeled as follows (Ashhad and Narayanan, 2013): where [B buf ] (ϭ450 M) and [Ca 2ϩ B buf ] represented the concentrations of free buffer and calcium bound buffer in the cell. k on and k off denoted the on and off rate constants for calcium binding to the buffer. Note that Eq. (16) constitutes a pseudo steady-state approximation, considering free buffer and calcium bound buffer to be in equilibrium. The value of K buf was set at 10 M (Klingauf and Neher, 1997;Fink et al., 2000;Ashhad and Narayanan, 2013).

Measurements
The excitability of the neuronal model was characterized by measuring its firing rate at 250 pA (f 250 ; Hz), action potential amplitude (V AP ; mV), input resistance (R in ; M⍀) and maximum impedance amplitude (|Z| max ; M⍀). The intrinsic response dynamics of the neuron were characterized by measuring the resonance frequency (f R ; Hz), strength of resonance (Q) and total inductive phase (⌽ L ; rad.Hz). These standard measurements ( Fig. 1) have been previously used to characterize CA1 pyramidal neurons. Firing rate at 250 pA was taken as twice the number of action potentials fired when a current of 250 pA was injected into the neuron for 500 ms (Fig. 1D,E). Action potential amplitude was calculated as the difference between the peak voltage of the action potential and the resting membrane voltage (-65 mV). Input resistance was measured by injecting currents of -50 pA to 50 pA, in steps of 10 pA, for 500 ms (Fig. 1B), recording the corresponding steady-state voltage deflection from -65 mV and taking the slope of the linear fit to the resulting V-I plot (Fig. 1C). To quantify the intrinsic response dynamics of the neuron, we injected a current in the form of a sinusoidal chirp stimulus of constant amplitude (50 pA) and a linearly increasing frequency (0 -25 Hz in 25 s; Fig. 1F). The impedance as a function of frequency (Z(f)) was obtained by dividing the Fourier transform of the voltage response by the Fourier transform of the injected chirp current. The impedance amplitude profile (Fig. 1G) was calculated as the magnitude of this impedance Z(f) which is given as follows: where Re(Z(f)) and Im(Z(f)) are the real and imaginary parts of the complex valued function Z(f). The frequency at which |Z(f)| is maximum is called the resonance frequency, f R . The impedance amplitude at the resonance frequency, |Z(f R )|, denotes the maximum impedance amplitude, |Z| max . The strength of resonance, Q, was computed as the ratio of the maximum impedance amplitude profile and the impedance amplitude at 0.5 Hz. The impedance phase profile ((f); Fig. 1H) was calculated as follows: Total inductive phase was defined as the area under the inductive part of curve (f) (Narayanan and  which is given by the following: (20)

Global sensitivity analysis
We used the global sensitivity analysis (GSA), a random sampling technique similar to previously used approaches (Bhalla and Bower, 1993;Foster et al., 1993;Goldman et al., 2001;Golowasch et al., 2002;Prinz et al., 2003Prinz et al., , 2004Achard and De Schutter, 2006;Tobin et al., 2006;Reid et al., 2007;Hobbs and Hooper, 2008;Weaver and Wearne, 2008;Taylor et al., 2009;Rathour andNarayanan, 2012b, 2014) to study the effect of the variability in and interactions among the passive and active properties of the neuron on activity-dependent calcium homeostasis. Four thousand models were generated by choosing a unique value for each of the 48 parameters (spanning passive and active properties) from uniform distributions around appropriate base values ( Table 1). The base values of parameters were obtained by hand-tuning a base model. As all the 4000 models could not be expected to have biologically realistic measurements, we validated the models by constraining the seven measurements to have values within experimentally determined ranges for CA1 pyramidal neurons (Table 2). Doing so, we found 78 of 4000 (ϳ2%) neuronal models to be physiologically realistic, and called these as valid models. All further analyses were performed on these 78 valid neuronal models.
where [Ca 2ϩ ] tgt (default value was 200 nM) represented the target value of cytosolic calcium concentration ([Ca 2ϩ ] c ) at which homeostasis should be maintained. The evolution of conductances (g i ) of individual channels from their respective mRNA (translation) was governed by: where g is the time constant for the translational process, and was set to be identically equal to 10 ms for all 12 conductances. We noted, with an additional set of simulations, that changes to the specific value of g merely altered the time-course toward reaching steady-state, but not the steady-state values of the conductances. For the mRNAs and the conductances to evolve as functions of the integral of error in calcium with reference to the target calcium (Eqs. 21, 22), we randomized the initial values of m i 's and g i 's to be very low. The time constants ( i ) for the evolution of mRNAs with reference to each channel were selected from the corresponding conductances of the valid models obtained from the GSA. Specifically, the time constants ( i ) for transcription were set to different values as follows : where i and j correspond to the 12 ion channels in our model, and k varies from 1 to 78 and corresponds to the number of valid models obtained from the GSA. As Eq. 23 defines i ratios from g i ratios, one of the i values needed to be set to obtain the other values. We set the i value associated with the sodium conductance at 10 ms, and computed the other i values from the appropriate con-ductance ratios as in Eq. 23. We noted, with an additional set of simulations, that changes to the specific value of the i value for the sodium channel altered the time-course toward reaching steady-state, but not the steady-state values of the conductances. The temporal evolution of transcription and translation were independently assessed for each of the 78 valid models by setting the i values to be dependent on the corresponding conductance ratios in each of these models . The other parameters associated with the model (the half-maximal activation voltages, time constants of ion channels and passive properties) were all derived from the respective valid model from where the i values were derived from.

Assessing state-dependent evolution of calcium homeostasis
In a manner similar to conditions observed under in vivo conditions, changes in afferent input to the neuron were presented as changes to synaptic receptors, specifically to AMPAR and NMDAR permeabilities (Eqs. 2-9). We tested the evolution of calcium homeostasis with two different types of afferent activity patterns that correspond to different behavioral states. The first corresponded to inputs received by hippocampal neurons during awake/REM-sleep state, and were modeled using a theta frequency (8 Hz) sinusoidal modulation Harvey et al., 2009) of AMPAR/NMDAR permeabilities (Fig. 7B). The amplitude of the sinusoidal permeability modulation was such that it was the minimum amplitude required to elicit action potentials in the neuron. The second type of afferent activity pattern reflected that during non-exploratory/non-REM sleep state ( Fig.  7C), where the neuron received SWRs riding on largely silent periods (Buzsáki, 1986). The shape of the SWR waveform was derived from English et al. (2014), and SWR amplitude was set such that the response to an SWR input resulted in a membrane voltage change of ϳ5-10 mV (English et al., 2014). Specifically, the functional form of SWR inputs was modeled as follows: where t represented time in milliseconds, and f ripple was the ripple frequency set at 150 Hz. Each SWR waveform lasted for around 150 ms, and was set to repeat at 3 Hz (Buzsáki, 1986;English et al., 2014). Afferent activity of the same type was continued until steady state of evolution in conductances (Eq. 22) was achieved (typically around 150 s), at which point measurements of intrinsic properties were noted and a switch in afferent activity was effectuated as necessary. All switches in afferent activity from theta oscillations to SWR inputs were initiated after a reset pulse to -65 mV for 1 s to avoid depolarizationinduced block observed in certain neurons. These bounds were extracted from experimental recordings (somatic recordings) presented by Narayanan and  and Narayanan et al. (2010).

Simulation details
All simulations were performed in the NEURON simulation environment (Carnevale and Hines, 2006) at -65 mV and 35°C with an integration time constant of 25 s. Temperature dependence of channel kinetics was appropriately accounted for with experimentally measured Q 10 values for each channel. The computational complexity of these simulations was enormous, as there were several differential equations associated with the activation and inactivation gates of the 11 ion channels (excluding the leak channel), the differential equation for the voltage and calcium (including calcium diffusion; Eqs. 10-17), and the differential equations governing the evolution of 12 mR-NAs and 12 conductances (Eqs. 21-23). The solutions to all these differential equations were computed at every time step (of 25 s) over a period of 150 s in achieving steady state for one type of afferent activity (a double switch in afferent activity is Ͼ3 ϫ 150ϭ450 s, running for several days in terms of simulation time), for each of the 78 valid models. Data analyses were performed using custom-written code with IGOR Pro (Wavemetrics) and statistical analyses were performed with the R Statistical Package (R Core Team, 2014).

Generation of a valid model population through global sensitivity analysis
As a first step in assessing state-dependence of cellautonomous calcium homeostasis in hippocampal neurons, we generated several biophysically realistic models of CA1 pyramidal neurons using the GSA approach. Specifically, we created a cylindrical base neuronal model (100 ϫ 100 m) containing 12 ion channels (Leak, NaF, KDR, KA, KM, HCN, CaL, CaT, CaN, CaR, BK, and SK) and AMPA and NMDA receptors (Fig. 1A). We hand-tuned the base model such that the seven intrinsic measurements namely, f 250 , V AP , R in , |Z| max , f R , Q, and ⌽ L ( Fig.  1B-H) were within experimentally determined ranges (Narayanan and Narayanan et al., 2010). We then uniformly sampled 48 parameters (spanning passive properties and densities/kinetics of channels in the neuron) from a range determined from the corresponding base model values to generate 4000 neurons (Table 1). We obtained seven measurements from each of these 4000 model neurons, and compared the measurements against their experimental counterparts. A model neuron was declared valid if all seven measurements of the model fell within their respective experimental bounds (Table 2). Upon imposing these experimental constraints on measurements from the 4000 models, we found 78 (ϳ2%) models to be valid. To test whether there were correlations between channel expression profiles and their kinetics, we asked whether there were pairwise correlations between the values of the 48 parameters associated with these 78 valid models. Consistent with previous results on hippocampal neurons (Rathour and Narayanan, 2012a, 2014), we found the parametric values to be weakly correlated (Fig 2A) with the range of correlation coefficients ranging from -0.6 to 0.6 ( Fig 2B-C). Importantly, of the 1128 correlation coefficients, 1125 were in the range of -0.4 to 0.4 suggesting weak pairwise relationships between parameters in the model. As the 78 models were valid models (referred to as GSA models in what follows) for hippocampal pyramidal neuron physiology, we used these for our analysis on state-dependence of intrinsic properties in regulating calcium homeostasis.
Under theta-frequency afferent activity, calciumdependent evolution of ionic conductances resulted in variable plasticity of intrinsic properties We set the time constants ( i ) for the evolution of mRNAs with reference to each channel from the corresponding conductances obtained from the valid model (Eq. 23), and implemented the evolution of calcium homeostasis (Eqs. 21, 22) for each of the 78 valid models. The neurons were presented with theta-frequency afferent activity, modeled as an 8 Hz sinusoidal permeability change in synaptic receptors. The mRNAs and conductances were allowed to evolve in time until a steady state was achieved in the cytosolic calcium concentration and the conductance values (Fig. 3A). During the initial phase of the evolution process, the voltage response of the model neuron to the sinusoidal input conductance corresponded to largeamplitude oscillations with small action potential amplitudes (Fig. 3A). We noted that this was consequent to the initial low values of all conductances, implying a large input resistance leading to large-amplitude voltage oscillations. The small action potential amplitudes, on the other hand, were consequent to the lower values of the spike-generating fast sodium conductance. As the calcium-dependent evolution progressed toward achieving the target calcium value, the voltage response corresponded to large-amplitude action potentials within each theta cycle (Fig. 3A).
Although the overall patterns of voltage evolution across the 78 models were similar to the example presented in Fig. 3A, there was significant variability in the final steady-state firing behavior of the model. Specifically, there was variability in the number of action potentials fired per cycle and there were models that did not fire at every cycle of the theta, with some of them exhibiting skipping spikes in alternate cycles and others skipping several cycles (Fig. 3B). These results imply that the maintenance of calcium homeostasis does not require and does not translate to maintenance of firing rate homeostasis, thereby establishing the dissociation between activity homeostasis and calcium homeostasis.
As neuronal firing properties are dictated by ionic conductances, we asked whether changes in conductance values during theta-dependent evolution was also variable. To answer this, we plotted the histogram of sodium conductances obtained after theta-dependent evolution for the 78 valid models, and found that the variability in firing patterns also reflected in these conductance values (Fig. 3C). We plotted the histogram solely for the sodium conductance, and not for all 12 conductances because the changes in all conductances are correlated given that a single transcription factor regulated all conductances Eqs. 21-23; Fig. 3A). These results imply that the maintenance of calcium homeostasis does not require and does not translate to maintenance of individual channels at specific conductance values, thereby establishing the dissociation between individual channelostasis (Rathour and Narayanan, 2014) and calcium homeostasis .
Variability in the percentage change in conductance values would imply that the intrinsic response properties of the neuron should have undergone variable plasticity during the evolution of calcium-dependent homeostasis. Therefore, we calculated the seven intrinsic measurements at steady state of the -dependent evolution, and compared these measurements with their corresponding GSA values for each of the 78 valid models. We found that -dependent evolution introduced variable plasticity in each of these measurements, with significant variability in the strength and sign of these changes (Fig. 3D). Importantly, although the conductance values obtained from GSA were those that satisfied experimental bounds (Table 2), several measurement values at steady state after -dependent evolution were not valid with reference to these experimental bounds (Fig. 3D). Together these results revealed that plasticity in neuronal intrinsic properties and in conductances could exhibit significant vari-ability when ionic conductances were allowed to evolve toward achieving calcium homeostasis with -oscillations as afferent inputs. These results imply that the mere maintenance of calcium homeostasis does not require and does not translate to maintenance of intrinsic measurements within a "valid" range, thereby establishing the dissociation between functional homeostasis and calcium homeostasis.

Sensitivity of intrinsic plasticity driven by calcium homeostasis to target calcium concentration and to the strength and frequency of afferent activity
What was the impact of changing the target calcium concentration, [Ca 2ϩ ] tgt (from the default value of 200 nM), on the temporal evolution of conductances and consequent plasticity in intrinsic properties? To address this, we repeated the -dependent evolution experiments (Fig. 3) with two other values for [Ca 2ϩ ] tgt , set at 100 and 300 nM. We observed neuronal firing at steady state (150 s) of -dependent evolution and found that the neuronal model fired more action potentials per -cycle upon increase in [Ca 2ϩ ] tgt . However, at a higher target value (300 nM), the neuron entered into depolarization-induced block with the membrane potential hovering at suprathreshold voltagelevels (Fig. 4A). Turning to steady-state values of intrinsic properties after -dependent evolution, we noted that variable plasticity in all seven intrinsic measurements was observed across all three values of [Ca 2ϩ ] tgt , with no qualitative differences observed in measurement variability with different values of [Ca 2ϩ ] tgt (Fig. 4B-H).
We next assessed the impact of afferent activity strength on -dependent evolution by altering the peakto-peak amplitude of the sinusoidal modulation in receptor permeability, with the sinusoidal frequency fixed at 8 Hz (Fig. 5). Whereas the neuron did not fire action potentials for lower strengths of afferent input ( Fig. 5A; permeability value P 1 ), at very high values of input permeability neurons entered into depolarization-induced block with average membrane potential around -30 mV ( Fig. 5A; permeability value P 4 ). In the mid-range between these two extremes, the number of action potentials fired per cycle increased with increase in afferent input strength. Assessing steady-state values of intrinsic properties after -dependent evolution, we noted that variable plasticity in all seven intrinsic measurements was observed across all tested values of input strength, with no qualitative differences observed in the measurement variability with different values of the sinusoidal peak-to-peak amplitude (Fig. 5B-H).
Finally, to understand the impact of afferent activity beyond the theta-frequency range on the evolution of calcium homeostasis and intrinsic properties, we picked samples from different frequency bands (delta: 1 Hz, theta: 8 Hz, slow gamma: 40 Hz, and fast gamma: 100 Hz) and repeated our simulations (until steady-state was achieved) with these afferent input frequencies. Neuronal response reflected theta-frequency band-pass structure of hippocampal pyramidal neurons, eliciting maximal firing response at the theta range with reduced response Bottom, Histograms of the percentage change in the seven measurements at steady state of theta-dependent evolution from the corresponding baseline GSA values, plotted for all 78 valid models. Percentage change in subthreshold measurements (R in , |Z| max , f R , Q, ⌽ L ) were computed only for those models that did not fire action potentials in response to the injected stimulus. In addition, models that showed very high percentage changes in ⌽ L were eliminated. The number of models (n) used for each histogram is mentioned in the respective panel.
with little or no firing at other frequency bands (Fig. 6A). Although the measurements at steady state of calciumdependent evolution did exhibit significant variability across the 78 models, their values did not have any specific dependence on the frequency of the input (Fig.  6B-H). For all further analyses, the default input frequency was 8 Hz to take into account the hippocampal theta rhythms, as mentioned above. Together, although neuronal firing pattern was heavily dependent on target calcium values and the strength/frequency of afferent inputs, intrinsic properties measured across all parametric combinations exhibited significant variability in the plasticity consequent to calcium-dependent evolution of conductances.

Neuronal conductances exhibited differential robustness to an intervening switch in the type of afferent activity
To study the state dependence of such a calcium homeostatic mechanism, we used experimentally well established differences in afferent activity to hippocampal neurons during different behavioral states and during different modules of the sleep cycle (Buzsáki, 2006;Miz-useki et al., 2011;Grosmark et al., 2012). Specifically, during REM sleep or exploratory behavior, afferent activity to the hippocampal CA1 neuron is dominated by oscillations in the theta-frequency (Fig. 7B), whereas during non-REM or non-exploratory behavior, the neuron predominantly receives SWRs riding on largely silent background (Fig. 7C). In this context, to study state dependence of the autonomous self-regulating calcium homeostasis mechanism, we switched the afferent activity to the neuron between -frequency oscillations and SWR activity. We analyzed two different temporal sequences of activity: theta-SWR-theta and SWR-theta-SWR sequences (Fig. 7A), with each afferent state of activity lasting for 150 s (to achieve steady-state conductance values with each phase of afferent activity).
Given the formulation of the dynamical evolution, the calcium levels expectedly converged to the target value of calcium, with transient changes occurring during the period that followed the switches. This maintenance of calcium homeostasis was effectuated by alterations in ionic conductances, which accommodated for the switch in afferent activity patterns (Fig. 7D-F). Although the extent of plasticity in ionic conductances was a continuum (see  (Table 2). below), we classified neurons into two distinct classes based on the amount of plasticity in ionic conductances. The first corresponded to neurons that underwent significant plasticity in their conductance values upon receiving the same type of afferent activity after an intervening period of switch to a different activity pattern (plastic neurons, ϳ80%; Fig. 7D,F). A second class of neurons restored their conductance values, exhibiting similar values upon receiving the same type of afferent activity after an intervening period of switch to a different activity pattern (robust neurons, 20%; Fig. 7E,G). Concurrently, neuronal response (firing) patterns with identical afferent activity were similar in robust neurons (Fig. 7E,G), but were significantly different in plastic neurons (Fig. 7D).
We assessed the percentage changes in sodium conductance after each switch in both sequences (across all 78 model neurons), and found that the changes in conductances introduced by the switches were significantly variable ( Fig. 8B-E for the theta-SWR-theta sequence and G-J for the SWR-theta-SWR sequence). We plotted the histogram solely for the sodium conductance, but not for all 12 conductances because the changes in all conductances are correlated given that a single transcription factor regulated all conductances. However, upon reverting back to the same type of afferent activity, a significant percentage of neurons were robust (with reference to conductance values, firing patterns and intrinsic properties) to the intervening switch to another type of afferent activity (Fig. 8E,J). The robustness of the neuron to an intermediate activity switch, however, was dependent on the specific temporal sequence of activity switch. Specifically, the percentage of neurons robust to an intervening Figure 5. The strength of afferent theta inputs critically regulated changes in intrinsic response properties during cell-autonomous self-regulation of calcium homeostasis. A, The steady state voltage response after theta-dependent evolution, for four different amplitudes of the 8 Hz input sinusoid (peak-to-peak value of sinusoidal permeability: black, 100 nm/s; red, 200 nm/s; green, 300 nm/s; blue, 400 nm/s). Note that at high values of sinusoidal amplitudes (e.g., 400 nm/s) the deflections are large along the hyperpolarized direction because of the large driving force for AMPA/NMDA receptors that mediate the sinusoidal oscillations. Along the depolarized direction, an action potential was elicited once the membrane potential crossed threshold, and the amplitude of the action potential did not cross the sodium reversal potential of ϩ55 mV. B-H, Histograms of the steady state measurement values (f 250 , B; V AP , C; R in , D; |Z| max , E; f R , F; Q, G; ⌽ L , H) for the 78 valid models, obtained after theta-dependent evolution with different sinusoidal amplitudes. The dashed lines in B-H represent the lower and upper bounds for the corresponding measurement (in that order) in the GSA model validation procedure (Table 2). period of activity switch was lower for the theta-SWRtheta (ϳ20%; Fig. 8E) compared with the percentage for SWR-theta-SWR sequence (ϳ90%; Fig. 8J). These results also reveal that the maintenance of calcium homeostasis does not necessarily require or translate to maintenance of individual conductances at specific values , and that significant plasticity in ionic conductances need not necessarily translate to significant changes in afferent-driven firing activity (Fig. 7F), thereby revealing a significant dissociation between individual channelostasis, activity/functional homeostasis and calcium homeostasis.
Neuronal intrinsic properties exhibited differential robustness to an intervening switch in the type of afferent activity Finally, we asked how state-dependent evolution of calcium homeostasis altered neuronal intrinsic proper-ties. To do this, we measured seven intrinsic properties (Fig. 1) at steady states of temporal evolution with activity pattern (theta or SWR) with the theta-SWRtheta (Fig. 9) and the SWR-theta-SWR (Fig. 10) sequences. Analysis of the evolution of intrinsic properties with switch in afferent activity revealed several important observations. First, intrinsic properties underwent significant plasticity as a consequence of these switches in types of afferent activity, and the sign and strength of this plasticity varied across different model neurons. Second, in the process of such activitydependent evolution, neuronal intrinsic properties did not necessarily fall into their established experimental bounds (Figs. 9, 10; Table 2). These results reveal a significant dissociation between functional and calcium homeostasis. Third, when neurons were presented with identical (the initial pattern preceding the switch) activity patterns after an intervening period of switch in  (Table 2). activity type, there were models whose intrinsic measurements did not restore to values before the intervening period. However, a significant proportion of neurons were robust to the intervening period of activity switch, where their intrinsic properties were restored when the activity pattern was switched back (Fig. 11). Similar to our observation with conductances, we noted that the percentage of neurons robust (in terms of changes in measurements) to an intervening period of activity switch was lower for the theta-SWR-theta ( Fig.  11A-G) compared with the percentage for SWR-theta-SWR sequence (Fig. 11H-N). We also noted that neurons fired significantly higher during the theta periods and lesser during SWR periods, which was partly due to overall reduction in intrinsic excitability of the neurons during SWR period (Figs. 7-11). . Switch in afferent activity between oscillations and SWRs triggered significant changes in ionic conductances during cell-autonomous self-regulation of calcium homeostasis. A, Experimental design for assessing state-dependence of ionic conductances during cell-autonomous self-regulation of calcium homeostasis. For each valid neuronal model obtained from GSA, ionic conductances were allowed to evolve toward achieving calcium homeostasis when afferent inputs were oscillations (or SWRs). At steady state of this evolution (150 s), inputs were switched to SWR (or oscillations). When this evolution reached steady state (150 s from the first switch), the input was switched back to oscillations (or SWR). B, Input received by a neuron during awake/REM sleep state was modeled as theta-frequency oscillations (8 Hz) injected as AMPAR and NMDAR permeabilities. C, Input received by a neuron during non-exploratory/non-REM sleep state was modeled as SWR inputs (inset, ripple frequency f ripple was 150 Hz; Eq. 24), repeating at a frequency of 3 Hz. These inputs were injected as AMPAR and NMDAR permeabilities into the neuronal model. D, E, Temporal evolution of the 12 ionic conductances in two different model neurons, where afferent activity switched from oscillations to SWR and back to oscillations. The calcium concentration (target [Ca]: 200 nM) and the conductances were allowed to reach steady state before either of the switches. The firing patterns of the neuron at the three steady states are also shown. F, G, Same as D, E, but with afferent activity switching from SWR to oscillations and back to SWR.

Figure 8.
Across different model neurons, switch in afferent activity between oscillations and SWRs triggered variable changes in ionic conductances during cell-autonomous self-regulation of calcium homeostasis. A, Schematic showing the temporal sequence of the experiment, along with notations for the temporal locations at which steady-state values of the sodium conductance (g Na ) were measured in the course of their calcium-dependent evolution. Note that this schematic represents a theta-SWR-theta temporal sequence in afferent activity, and the notations here hold for B-E. All histograms depict statistics across the 78 valid models obtained after GSA. B, Histogram of g Na values at the steady state of the evolution with -frequency oscillations as afferent inputs. Inset, Histogram of base values of g Na obtained from GSA. C, Histogram of percentage changes in g Na measured at steady state of evolution with SWR inputs (after 1 -SWR), computed with reference to the steady state value after evolution with -frequency oscillations ( 1 ). D, Histogram of percentage changes in g Na measured at steady state of evolution with -frequency oscillations (after 1 -SWR-2 ), computed with reference to the steady state value after evolution with SWR inputs (after 1 -SWR). E, Histogram of percentage changes in g Na measured at steady state of evolution with -frequency oscillations (after 1 -SWR-2 ), computed with reference to g Na measured at steady state of evolution after 1 . F-J, Same as A-E, but for a SWR-theta-SWR temporal sequence in afferent activity, with notations for conductance values shown in F.
Together, our results suggest that neuronal conductances and intrinsic properties could undergo significant plasticity toward maintenance of calcium homeostasis in the face of behavioral-state dependent changes in afferent activity. Although a significant population of neurons exhibited robustness to an intervening switch in the type of afferent activity, there were neurons that manifested significantly distinct intrinsic properties upon restoration of the type of afferent activity after the intervening switch. These results suggest that neuronal conductances and intrinsic properties exhibited differential robustness to an intervening switch in the type of afferent activity, apart from demonstrating significant dissociations between functional/channel/firing-rate and calcium homeostasis.

DISCUSSION
The key conclusion of this study is that neuronal ionic conductances and intrinsic properties could undergo significant plasticity in the process of maintaining calcium homeostasis through a regime of behavioral state-dependent Figure 9. A Theta-SWR-theta switch in afferent activity introduced significant changes to neuronal intrinsic response properties during cell-autonomous self-regulation of calcium homeostasis. A, Experimental design representing a theta-SWR-theta temporal sequence in afferent activity for assessing state-dependence of intrinsic response properties during cell-autonomous self-regulation of calcium homeostasis. The arrows represent time points at which intrinsic measurements were computed, and also associated with different symbols used in B-H. B-H changes in afferent activity. This form of intrinsic plasticity driven by calcium homeostasis was sensitive to the specific value of target calcium and the strength and frequency of afferent activity. We assessed the impact of behavioral state-dependence of afferent activity on calcium homeostasis by switching afferent activity pattern between theta-frequency oscillations (REM sleep/exploratory behavior) and SWR activity (non-REM sleep/nonexploratory behavior). Switches in activity patterns resulted in variable plasticity in ionic conductances and neuronal measurements, with the sign and strength of plasticity dependent on the specific type of activity pat-tern and on the temporal sequence of switch. Additionally, our analysis with temporal sequences of activity switch revealed the presence of two classes of neurons. One that showed significant plasticity in conductance values and in intrinsic properties when presented with identical activity patterns after an intervening period with a different activity pattern, and a second that restored its conductance values and intrinsic properties after the intervening period. The percentage of robust versus plastic neurons was variable in a manner that was critically dependent on the specific sequence of switch in activity. Finally, our results also reveal that the maintenance of calcium homeostasis Figure 10. An SWR-theta-SWR switch in afferent activity introduced significant changes to neuronal intrinsic response properties during cell-autonomous self-regulation of calcium homeostasis. A, Experimental design, representing a SWR-theta-SWR temporal sequence in afferent activity, for assessing state-dependence of intrinsic response properties during cell-autonomous self-regulation of calcium homeostasis. The arrows represent time points at which intrinsic measurements were computed, and also associated with different symbols used in B-H. B-H does not necessarily translate to the emergence of individual channelostasis or of functional homeostasis (across physiological measurements, including firing rate), thereby establishing critical dissociations between these forms of homeostasis. In what follows, we present certain physiological implications for our conclusions, along with detailed analyses on model assumptions and future directions with testable predictions.

Behavioral state-dependent changes in neuronal intrinsic properties
Afferent activity patterns in hippocampal pyramidal neurons exhibit well-established distinctions that reflect the behavioral state of the animal (Anderson et al., 2007). The literature on behavioral state-dependent changes in neuronal activity has largely focused on postulates that involve synaptic normalization and synaptic plasticity during activity switches (Tononi and Cirelli, 2006;Chauvette et al., 2012;Grosmark et al., 2012;Barnes and Wilson, 2014), or on how dendritic nonlinearities could differentially process different forms of activity patterns (Gasparini and Magee, 2006). Our results point to a novel form of behavioral state-dependent plasticity in neuronal intrinsic properties, emerging as a direct consequence of the requirement to maintain calcium homeostasis in the face of changes in afferent activity. In this context, it is important that future studies on behavioral-state dependence of neuronal physiology also consider the role of neuronal intrinsic properties and ionic conductances toward changes in overall firing rate (Tononi and Cirelli, 2006;Grosmark et al., 2012), without limiting the analysis to synaptic changes and neuromodulatory influences (Le-Masson et al., 1993;Siegel et al., 1994;Turrigiano, 1999;Turrigiano andNelson, 2000, 2004;Prinz et al., 2004;Triesch, 2007;Turrigiano, 2007;Marder, 2011;Honnuraiah and Narayanan, 2013;Marder et al., 2014;O'Leary et al., 2014). Additionally, although it is established that the propensity of different activity patterns is higher during specific behavioral states, a constant afferent pattern used in the model is physiologically infeasible under in vivo conditions. Future experiments should therefore explore the relationships between time constants for channel plasticity under in vivo conditions and the temporal extent of specific afferent activity patterns. Although in vitro experiments have shown plasticity in several ion channels to occur in a matter of minutes and have demonstrated that changes in ion channels and intrinsic properties can be concurrent with synaptic changes (Frick et al., 2004;Fan et al., 2005;Kim et al., 2007;Narayanan and Johnston, 2007;Lin et al., 2008;Losonczy et al., 2008;Rosenkranz et al., 2009;Narayanan et al., 2010;Remy et al., 2010;Shah et al., 2010), the temporal aspects of activity-dependent plasticity in ion channels under in vivo conditions needs to be explored in detail, in a manner that accounts for channelostasis individually and collectively Figure 11. Across different model neurons, switch in afferent activity between oscillations and SWRs triggered variable changes in neuronal intrinsic response properties during cell-autonomous self-regulation of calcium homeostasis. A-G, 1 -SWR-2 constitutes the temporal sequence of afferent activity, with notations shown in Figure 9A. Histogram of percentage changes in the seven intrinsic response properties (f 250 , A; V AP , B; R in , C; |Z| max , D; f R , E; Q, F; ⌽ L , G) measured at steady state of evolution with -frequency oscillations (after 1 -SWR-2 ), computed with reference to the corresponding base value of the intrinsic property measured at steady state of evolution after 1 . Percentage changes in subthreshold measurements (R in , |Z| max , f R , Q, ⌽ L ) were computed only for those models that did not fire action potentials in response to the injected stimulus. In addition, models that showed very high percentage changes in ⌽ L were eliminated. The number of models (n) used for each histogram is mentioned in the respective panel. H-N, Same as A-G, but for a SWR 1 -theta-SWR 2 temporal sequence of afferent activity, with notations shown in Figure 10A. (Hanus and Schuman, 2013;Rathour and Narayanan, 2014;Anirudhan and Narayanan, 2015). Additionally, such experiments could test for variability in such statedependent intrinsic plasticity predicted by our model, apart from addressing the impact of such variability on neurophysiology and behavior.
As a specific instance, the existence of plasticity in intrinsic properties (in addition to synaptic plasticity) would expand the putative mechanisms that could be involved in memory consolidation, a postulated function for sleep rhythms (Siegel, 2001;Stickgold et al., 2002;Walker and Stickgold, 2004;Stickgold, 2005;Marshall and Born, 2007). As activity-dependent plasticity in neuronal intrinsic properties is well established (Zhang and Linden, 2003;Kim and Linden, 2007;Johnston and Narayanan, 2008;Shah et al., 2010;Narayanan and Johnston, 2012), the exploration of the postulate that memory consolidation is effectuated through intrinsic plasticity (in conjunction with synaptic changes) is a critical prediction that needs rigorous experimental evaluation. Given this postulate where the possibility of intrinsic changes exists, interpretation of observations from experiments that involve replay or disruption of specific activity pattern (Girardeau et al., 2009;Ego-Stengel and Wilson, 2010;Jadhav et al., 2012;Barnes and Wilson, 2014) should also account for changes in neuronal intrinsic properties that might have been brought about by the specific activity pattern or lack thereof.

Would behavioral state-dependent intrinsic plasticity in a neuronal compartment be dependent on its somatodendritic location?
Although the macroscopic activity patterns recorded in the hippocampus show theta-SWR switches during different behavioral states, it is evident that there are subtle, yet significant, differences in afferent activity at different somatodendritic locations along a hippocampal pyramidal neuron (Colgin et al., 2009;Colgin and Moser, 2010;Bieri et al., 2014;Schomburg et al., 2014). Additionally, there are well-established differences in localization of different ion channels (Johnston et al., 1996;Magee, 2000;Migliore and Shepherd, 2002;Johnston and Narayanan, 2008;Spruston, 2008;Narayanan and Johnston, 2012), in the locus of plasticity in these channels (Frick et al., 2004;Narayanan and Johnston, 2007;Losonczy et al., 2008;Narayanan et al., 2010;Shah et al., 2010), and in calcium source localization and calcium propagation (Magee and Johnston, 1995;Augustine et al., 2003;Hertle and Yeckel, 2007;Ross, 2012) across the somatodendritic arbor. Consistent with this, and given our results with a single-compartmental model (necessitated by the tremendous computational complexity of calciumdependent evolution), we postulate that the behavioralstate dependent intrinsic plasticity reported here would be dependent on the somatodendritic location of the neuronal compartment.
Testing this postulate would require development of specific experimental techniques and computational models to assess the impact of the self-regulating evolution of calcium homeostasis on changes in localization profiles of ion channels across the somatodendritic arbor. Experimental procedures would require direct measurement of somatodendritic channel properties under in vivo conditions, during different stages of sleep or behavior, with location-dependent afferent activity monitored in parallel (Agarwal et al., 2014;Schomburg et al., 2014). These experiments would provide direct answers to questions on whether dendritic channel localization profiles change as a function of activity switch (during specific stages of sleep or behavior), and if neurons implement an efficient form of neural coding that accounts the statistics of their afferent activity (Stemmler and Koch, 1999;Simoncelli and Olshausen, 2001;Simoncelli, 2003;Narayanan and Johnston, 2012).
Computational models, on the other hand, would have to explicitly account for somatodendritic differences in ion channel profiles, physiological measurements and calcium source localization (Berridge, 2006;Vacher et al., 2008;Kotaleski and Blackwell, 2010;Nusser, 2012;Ashhad and Narayanan, 2013;Hanus and Schuman, 2013;Rathour and Narayanan, 2014), apart from ensuring that spatial compartmentalization of the neuronal model is based on the calcium space constant rather than the electrical space constant (Koch and Zador, 1993;Zador and Koch, 1994;Ashhad and Narayanan, 2013). Second, these models would have to address the question on whether calcium homeostasis is maintained globally or locally, and ask if localization profiles of different channels were emergent properties consequent to the cellautonomous calcium homeostasis process (Siegel et al., 1994;Rabinowitch andSegev, 2006, 2008). Third, consistent with the existence of several enzymes that act as calcium sensors (Liu et al., 1998), and the existence of several activity-dependent transcription factors (Dolmetsch, 2003), future models should extend beyond the single transcription factor-based analysis used in our model (see below). Finally, in our study, we have not accounted for neuromodulatory influences, and have resorted to a simplistic classification of afferent activity as theta versus SWR. Future experimental and computational studies on state-dependent calcium homeostasis should also account for differences in neuromodulatory activity during REM versus exploratory behavior and non-REM versus non-exploratory behavior (Vertes, 1984;Steriade, 2004;McCarley, 2007;Lee and Dan, 2012), and the impact of each of such differential neuromodulatory activity on ion channels and their plasticity (Fisher and Hoffman and Johnston, 1999;Cantrell and Catterall, 2001;Marder and Thirumalai, 2002;Perez-Reyes, 2003;Biel et al., 2009;Marder, 2012;Marder et al., 2014).
What is the impact of such location-dependent plasticity in somato-dendritic channel properties on neuronal physiology and behaviorally relevant neural computation? First, such plasticity would alter the following physiological characteristics of these neuronal structures, each of which is known to exhibit location-dependence in a manner dependent on specific ion channel combinations: spectral selectivity (Narayanan and Hu et al., 2009;Das and Narayanan, 2014), coincidence detection London and Häusser, 2005;Sjöström et al., 2008;Das and Narayanan, 2015), impedance phase (Narayanan and Vaidya and Johnston, 2013), and supralinear summation (Losonczy and Losonczy et al., 2008;Sjöström et al., 2008;Spruston, 2008;Takahashi and Magee, 2009). Importantly, such state-dependent plasticity in channel conductances would not just reflect as changes in intrinsic neuronal physiology, but also express as changes in the amplitude and phase of local field potentials (LFP) and associated neuronal spike phases Schomburg et al., 2012;Einevoll et al., 2013;Reimann et al., 2013;Sinha and Narayanan, 2015). Future studies should test if such intrinsically-driven changes in LFP and spike phase could potentially form a cellular substrate for REM-shifting neurons (Mizuseki et al., 2011), a scenario where the phase shift is a consequence of intrinsic plasticity that occurred during an intervening switch to non-REM activity (Figs. 7-9, 11). Finally, extrapolating from recent studies that have demonstrated the importance of dendritic nonlinearities to place cell formation (Bittner et al., 2015;Sheffield and Dombeck, 2015), changes in dendritic sodium/calcium/potassium/HCN channels would alter the propensity for generating dendritic plateau potentials (Golding et al., 1999;Gasparini and Magee, 2002;Gasparini et al., 2004;Losonczy and Magee, 2006;Tsay et al., 2007;Losonczy et al., 2008), potentially resulting in changes in the place cell properties of the associated hippocampal neurons (Bittner et al., 2015).

Implications for the assumption on a single transcription factor
In our model, we have assumed that the channel conductances are regulated by a single transcription factor , an assumption that significantly oversimplifies the complexities of neuronal transcription, where multiple transcription factors coexist (Bading et al., 1993;Dolmetsch, 2003;Lein et al., 2007;Alberini, 2009). This assumption implies that proportions of changes in channel conductances are correlated (Fig. 3A), resulting in correlated channel expression profiles . Although this assumption was motivated by correlated expression profiles of ion channels in certain neuronal subtypes (Toledo-Rodriguez et al., 2004;Schulz et al., , 2007Tobin et al., 2009;Amendola et al., 2012), detailed quantitative analysis of channel conductances and mRNA expression has not been performed in single hippocampal neurons. In the absence of such experimental data, not just at the cell body, but across the somatodendritic arbor (Hanus and Schuman, 2013;Rathour and Narayanan, 2014), model-based extrapolations about correlations in expression profiles of hippocampal channels or mRNAs would be incorrect, because the model outcome is just a direct consequence of the assumption involving a single transcription factor. Therefore, we restrict our inferences from this simple model to: (1) variable state-dependent plasticity of ionic conductances and intrinsic properties toward cell-autonomous maintenance of calcium homeostasis and (2) significant dissociation between different forms of homeostasis (see below), which also form testable predictions from our analysis. The question on whether the expression profiles of different channels/mRNAs are correlated or lack significant correlation in the presence of multiple transcription factors needs to be rigorously addressed both from experimental as well as theoretical standpoints.
Incorporation of multiple transcription factors into a model for cell-autonomous calcium homeostasis has been reported to result in unbounded production of mR-NAs and channels ("windup"), leading to eventual loss of regulatory control . Although this constitutes a significant impediment to the incorporation of multiple transcription factors into models, this analysis was performed in a manner where the different transcription factors were independent of each other . Future theoretical studies should explore the possibility of avoiding such windup by coupling the multiple calcium sensors and multiple transcription factors through established signaling motifs, including negative feedback mechanisms (Thattai and van Oudenaarden, 2001;Losick and Desplan, 2008;Yu et al., 2008;Kotaleski and Blackwell, 2010;Cheong et al., 2011). Experimental studies should explore the relationships between the different transcription factors, mRNAs and channel conductances across the somatodendritic arbor of single hippocampal neurons (Dolmetsch, 2003;Hanus and Schuman, 2013).

Dissociations between different forms of homeostasis
It is clear from our analyses here, and from several others in the literature, that there are significant dissociations between different forms of homeostasis. First, homeostasis in functional properties, including synaptic plasticity profiles (Anirudhan and Narayanan, 2015), could emerge with disparate conductance values for the constituent ion channels and synaptic conductances (Golowasch and Marder, 1992;Golowasch et al., 2002;MacLean et al., 2005;Schulz et al., , 2007Taylor et al., 2009;Rathour andNarayanan, 2012a, 2014), suggesting that homeostatic maintenance of single channels at specific conductance values is not essential for maintaining functional or plasticity profile homeostasis. Second, for maintenance of calcium homeostasis across neurons in a network (O'Leary et al., 2013), or in a neuron that receives state-dependent switch in afferent activity (Figs. 3-11), it is not essential that functional homeostasis across different measurements is maintained. Specifically, despite maintenance of calcium homeostasis across models receiving identical temporal evolution of afferent activity, we noted that the conductance values (Figs. 3, 8) and physiological measurements (9)(10)(11) were significantly variable across these models, with some models manifesting measurements beyond what is expected from CA1 pyramidal neurons. Additionally, a significant proportion of neurons did not restore their intrinsic properties despite restoration of specific type of activity after an intervening switch to a different type of activity (Figs. 9-11). Finally, although calcium homeostasis was achieved across all neuronal models, there was significant variability (across models) in the firing rate and in pattern of firing in response to the same afferent activity (Figs. 3,7;O'Leary et al., 2013O'Leary et al., , 2014. Together, these results clearly establish that maintenance of calcium homeostasis neither translates to nor follows from any of channel/ functional/firing-rate forms of homeostasis, outlining critical dissociations between these forms of homeostasis.
Future experiments should therefore explore the specific form of homeostasis maintained by individual neurons and their dendritic arbor when subjected to behavioral state-dependent afferent activity during different stages of sleep and behavior. Do channel properties/ localization and intrinsic functional properties change in the process of maintaining calcium homeostasis? Or, do neurons implement a mechanism where all forms of homeostasis-including that in spatially distributed channel properties, intrinsic response properties, plasticity profiles, and calcium levels-are concurrently maintained across all somatodendritic locations of the neuron, with firing rate homeostasis emerging as an overall consequence? Should a homeostatic mechanism account not just for the average calcium level in a neuron, but also make provisions for the homeostasis in input-output profiles, intrinsic response properties and synaptic/intrinsic plasticity profiles of the neuron through synergistic interactions between synaptic and intrinsic neuronal properties (LeMasson et al., 1993;Liu et al., 1998;Triesch, 2007;Turrigiano, 2011;Honnuraiah and Narayanan, 2013;O'Leary et al., 2014;Anirudhan and Narayanan, 2015)? Finally, whereas homeostasis covers only one aspect of neuronal function, the other core function (especially of hippocampal neurons), is encoding of new information. Juxtaposed against questions on various forms of homeostasis is the fundamental issue of how neurons change their properties toward encoding new information, without jeopardizing any or some forms of homeostasis. Therefore, further exploration into behavioral statedependent evolution of homeostasis should account for encoding as a critical aspect of neuronal function that depends on changes in intrinsic and/or synaptic properties (Narayanan and , apart from exploring the relationships between different forms of homeostasis.