Home | About | Journals | Submit | Contact Us | Français |

**|**Front Comput Neurosci**|**v.4; 2010**|**PMC2856590

Formats

Article sections

- Abstract
- Introduction
- Materials and Methods
- Results
- Discussion
- Conflict of Interest Statement
- Supplementary Material
- References

Authors

Related links

Front Comput Neurosci. 2010; 4: 6.

PMCID: PMC2856590

Edited by: Peter Dayan, University College London, UK

Reviewed by: Mickey London, University of Hertfordshire, UK; Volker Steuber, University of Hertfordshire, UK

*Correspondence: Fidel Santamaria, Department of Biology, University of Texas at San Antonio, San Antonio, TX 78249, USA. e-mail: ude.astu@airamatnas.ledif

Received 2009 July 6; Accepted 2010 March 5.

Copyright © 2010 Coop, Cornelis and Santamaria.

This is an open-access article subject to an exclusive license agreement between the authors and the Frontiers Research Foundation, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are credited.

This article has been cited by other articles in PMC.

Using an electrophysiological compartmental model of a Purkinje cell we quantified the contribution of individual active dendritic currents to processing of synaptic activity from granule cells. We used mutual information as a measure to quantify the information from the total excitatory input current (*I*_{Glu}) encoded in each dendritic current. In this context, each active current was considered an information channel. Our analyses showed that most of the information was encoded by the calcium (*I*_{CaP}) and calcium activated potassium (*I*_{Kc}) currents. Mutual information between *I*_{Glu} and *I*_{CaP} and *I*_{Kc} was sensitive to different levels of excitatory and inhibitory synaptic activity that, at the same time, resulted in the same firing rate at the soma. Since dendritic excitability could be a mechanism to regulate information processing in neurons we quantified the changes in mutual information between *I*_{Glu} and all Purkinje cell currents as a function of the density of dendritic Ca (*g*_{CaP}) and Kca (*g*_{Kc}) conductances. We extended our analysis to determine the window of temporal integration of *I*_{Glu} by *I*_{CaP} and *I*_{Kc} as a function of channel density and synaptic activity. The window of information integration has a stronger dependence on increasing values of *g*_{Kc} than on *g*_{CaP}, but at high levels of synaptic stimulation information integration is reduced to a few milliseconds. Overall, our results show that different dendritic conductances differentially encode synaptic activity and that dendritic excitability and the level of synaptic activity regulate the flow of information in dendrites.

The contribution of the passive dendritic properties of neurons on processing of synaptic activity is well established (Rall, 1962a,b; Marr, 1969; Pellionisz and Szentagothai, 1973, 1974); however, possibly all central nervous system neurons express voltage sensitive dendritic conductances (Llinas and Yarom, 1981; Hirst and McLachlan, 1986; Hockberger et al., 1989; Huguenard et al., 1989; Masukawa et al., 1991; Segev and London, 1999). Active conductances can carry more current than passive conductances and, thus, have a stronger contribution to the processing and transfer of synaptic activity in dendritic trees (Gollo et al., 2009). Furthermore, the multiple types of conductances found in dendrites could differentially filter synaptic activity depending on its temporal and spatial distribution (Poirazi and Mel, 2001; Migliore and Shepherd, 2002). Unfortunately, the lack of experimental access to determine the distribution and kinetics of active membrane conductances in dendrites has not allowed a thorough study of their influence on processing of synaptic activity.

In order to track the transformation of synaptic activity, and thus obtain a representation of the transfer function of the dendrite, it is necessary to determine how synaptic activity is distributed over different dendritic currents. In many cases the variability observed in dendritic currents is driven by synaptic activity and can be regulated by interactions among conductances. Quantifying the amount of variability due to synaptic activity or due to internal interactions can be done using with statistical tools derived from information theory (Shannon, 1948; Shannon and Weaver, 1949; Rieke, 1997). Particularly, determining mutual information (Foffani and Priori, 2007; Maffezzoli et al., 2008; Quian Quiroga and Panzeri, 2009) between synaptic activity and specific dendritic currents could result in a quantitative method to measure the contribution of individual dendritic currents on synaptic processing (Stemmler and Koch, 1999; London et al., 2002). Furthermore, the calculation of the entropy, conditional entropy and mutual information can help elucidate how much variability observed in the dendritic currents is due to the nature of the synaptic input or due to intrinsic properties of the dendrite.

We present the results of calculating mutual information between the net excitatory input and different dendritic currents using a detailed biophysical model of a Purkinje cell (De Schutter and Bower, 1994a,b). The model was stimulated with randomly distributed excitatory and inhibitory synapses that resembled granule cell and interneuron activity. We found that the calcium (CaP) and a calcium activated potassium (Kc) currents carried the largest amount of information. While the firing rate of the Purkinje cell remained constant for different combinations of synaptic input, mutual information was sensitive to such changes, thus disambiguating synaptic activity in dendrites. Our results suggest that dendritic excitability modulated by CaP and Kc channels is most effective in regulating the mutual information, and thus any possible processing of that information. Conceptually, our approach treats each dendritic conductance as an *information channel* that carries a certain amount of information depending on each conductance density. Since plasticity of dendritic excitability might be an important component of general learning and memory mechanisms in neurons (Xu and Kang, 2005; Komendantov and Ascoli, 2009), our results provide a method to quantify the effects of such changes on synaptic information processing in individual dendritic conductances.

We used a previously published Purkinje cell model (De Schutter and Bower, 1994b,c) with updated synaptic kinetics (Figure (Figure1A).1A). Please refer to the cited work for kinetic details of each conductance. In brief, there is one excitatory synaptic contact made with each dendritic spine and 1,695 inhibitory GABAa-type synaptic contacts distributed at random across the dendrites. The model contains the following dendritic conductances (Figure (Figure1B):1B): two types of Ca channels, a P-type, CaP (Llinas et al., 1989) and a T-type, CaT (Kaneda et al., 1990); two types of Ca-activated K^{+} channels, a BK-type (Latorre et al., 1989) referred to here as Kc and a K2-type, K2 (Gruol et al., 1991); and a persistent K^{+} channel or Km-type (Yamada et al., 1989). The soma had two types of sodium channels, a fast current, NaF (Hirano and Hagiwara, 1989) and a slow persistent current, NaP (French et al., 1990); one type of calcium current T-type; and four types of potassium channels, anomalous rectifier, Ih (Spain et al., 1987); delayed rectifier, Kdr (Yamada et al., 1989); persistent potassium, Km; and an A-type potassium channel, KA (Hirano and Hagiwara, 1989). In depth analysis of the model can be found elsewhere (De Schutter and Bower, 1994a,b,c; Jaeger et al., 1997; Gundappa-Sulur et al., 1999; Santamaria et al., 2002, 2007; Santamaria and Bower, 2005; De Schutter and Steuber, 2009). The updated synaptic properties are described in Santamaria et al. (2007).

The simulations consisted in randomly activating all the excitatory and inhibitory synapses at constant Poisson firing rates. As explained in the Results we used four different combinations of excitatory and inhibitory synaptic activity that resulted in the same firing rate at the soma of the Purkinje cell. We ran simulations for up to 400s saving the value of all dendritic and synaptic currents every 100μs. In order to avoid initial condition effects the first 5s of all traces were not used for the analysis. Simulations were run with a pre-release version of the new GENESIS 3 software^{1}. Since each single simulation required about 25h to complete we used a high performance computing cluster operated by the Computational Biology Initiative at UTSA^{2}.

In order to simplify the analysis we monitored the total value of the synaptic or dendritic currents. This approximation to the state of dendritic currents has allowed us to analyze this model in detail in the past with testable experimental predictions (Jaeger and Bower, 1999; Santamaria et al., 2007). We also chose to use the total excitatory current because then the results of our study could be mapped to dynamic current clamp experiments (see Discussion). We also ran our analysis under two other conditions. In the first one we added all currents over all the dendrites except the thick smooth dendrites (primary). The second condition was adding all the current exclusively in the spiny dendrites (tertiary), without the primary and secondary. The results were practically the same as the ones reported for the summed currents over the entire dendrite, indicating that all the integration is done at the level of tertiary dendrites.

For the purposes of comparing the changes due to background activity we normalized the value of all currents from 1–100 and binned the data in 1000 equally spaced bins. All the analyses described in this paper were performed with the normalized current values.

Initial characterization of currents was done by calculating the histograms under all the different combinations of synaptic activity. Further analysis consisted in calculating the cross-correlation between the excitatory synaptic input (*I*_{Glu}) and dendritic currents.

Since we wanted to characterize the amount of activity in dendritic currents due to synaptic stimulation we used information theory measurements. This process required calculating the entropy (*H*) of each current, the conditional entropy between each dendritic current and *I*_{Glu}, and finally the mutual information (*I*).

The entropy was calculated as

$$H(x)={\sum}_{i}p({x}_{i}){\mathrm{log}}_{2}p({x}_{i})$$

(1)

where *p* is the probability of seeing value *x _{i}*.

The conditional entropy is a number that describes the variability of a signal given that the value of a second signal is known. This is a useful measurement when trying to understand the variability of the dendritic currents as a function of the synaptic current. We calculated the conditional entropy as

$$H\left(y|x\right)={\text{\Sigma}}_{j}{\text{\Sigma}}_{i}p\left({x}_{i},{y}_{j}\right){\mathrm{log}}_{2}p\left({y}_{j}|{x}_{i}\right)$$

(2)

In our case, *x* is the input signal (*I*_{Glu}) and *y* any of the dendritic currents. Conditional probability distributions matrices were calculated based on the binned traces. The conditional probability distributions did not result in distributions that allocated values to every matrix entry. This lack of probability space coverage was due to the biophysical restriction of the particular conductances. Coupling and long time constants also contributed to having areas in the probability density matrix that were densely populated while other remained empty.

Finally, the mutual information was calculated with eq. 3:

$$I(y|x)=H(y)-H(y|x)$$

(3)

The mutual information is usually expressed in bits and gives the capacity of a neuron to discriminate between different stimuli that are part of the stimulus, in this case *I*_{Glu}. For example, if the mutual information is equal to 5 bits, this means that a particular dendritic current is able to distinguish 2^{5} or 32 different stimuli.

It is well known that the value of *I* can be biased due to the binning process and finite size of the data being analyzed (Panzeri et al., 2007). We used a recently developed toolbox in Matlab (Natick, MA, USA) that allows the accurate calculation of the different information measurements and compensation for potential biases (Magri et al., 2009). Independently of using such compensations we determined that a 400-s simulation was long enough by plotting the value of *I* for any pair of *I*_{Glu} and dendritic current as a function of the time window of observation from 5–400s. This analysis showed that after 300s the value of *I*, without bias compensation, reached an asymptotic value (**Figure S1** in Supplementary Material). In any case, the value of *I* can be biased if the joint probability distribution of the two traces being analyzed, for example *I*_{Glu} and *I*_{CaP}, is scattered and does not fill out the joint probability space, in our case a matrix of 1000×1000 entries. The ratio *N*/*m* has been shown to determine the strength of such a bias, where *N* is the number of non-zero entries in the joint probability distribution and *m* number of non-zero entries of probability distribution of the stimulus. If *N*/*m* is less than 1 then the value obtained from calculating the mutual information is biased. As shown in the supplementary materials all our simulations had an *N*/*m*>1.

The objective of our work was to quantify the contribution of dendritic excitability on information processing in dendrites. For this purpose we used an available Purkinje cell model since this type of cell has extensive dendritic conductances and this particular model has been well validated with experimental data (Jaeger and Bower, 1994, 1999; Jaeger et al., 1997; Santamaria et al., 2007).

Information in Purkinje cells flows from synapses, to dendrites, to soma (Santamaria et al., 2007). As opposed to cortical pyramidal cells, Purkinje cells have no backpropagating action potentials that carry somatic activity into the dendrites (Vetter et al., 2001). Thus, Purkinje cells are ideal to quantify the incremental contribution of dendritic channels to information processing. Since excitatory synapses from parallel fibers to Purkinje cells carry the activity from outside the cerebellum (Bower, 2002) we assumed this stimulus to be the input signal. We characterized the input signal as the total excitatory current (*I*_{Glu}) because experimental techniques, such as dynamic current clamp, could be used to test modeling predictions based on this assumption (Suter and Jaeger, 2004). Purkinje cells *in vivo* receive excitatory synaptic activity from parallel fiber and inhibitory activity from interneurons; we simulated these types of inputs with random input firing rates drawn from a Poissonian distribution. Firing rates at the soma can be the result of multiple combinations of excitation and inhibition in the dendrite (Santamaria et al., 2002; Santamaria and Bower, 2005), for that reason we selected various pairs of uncorrelated background excitatory and inhibitory activity that evoked similar firing rates at the soma. Using different combinations of excitatory and inhibitory input with the same somatic firing rate allowed us to determine information processing changes in dendrites due to different levels of synaptic activity. The background input firing rates spanned a wide range of rates 12–56Hz and 0.5–2.0Hz for excitation and inhibition, respectively (Santamaria and Bower, 2005).

We initially quantified the effects of the random excitatory input on the histogram of values of each dendritic current. For the different combinations of background excitatory and inhibitory input we found that while the histogram of *I*_{Glu} remained practically unchanged (Figure (Figure2A)2A) the shape of the distribution of the CaP (*I*_{CaP}, Figure Figure2B)2B) and Kc currents (*I*_{Kc}, Figure Figure2C)2C) varied. All other dendritic currents (*I*_{CaT}, *I*_{Km}, and *I*_{K2}) did not show changes in the shape of their distribution as a function of background levels of synaptic stimulation (Figures (Figures2D–F).2D–F). We did not analyze *I*_{NaF}, *I*_{Kdr}, *I*_{A}, *I*_{NaP}, and *I*_{h} because they are only present in the soma and the thick proximal dendrite. The mean firing rate (*F _{M}*), calculated as the number of spikes across the entire simulation divided by the simulation time (400s) shows that the different combinations of synaptic activity had little influence (90.3±2.4 spikes/s), Figure Figure2G,2G, similarly the distributions arising from calculating the inter-spike interval histograms for all combinations of synaptic input (Figure (Figure2H)2H) differ in less than 1% in the value of their expected value and standard deviation.

The analysis presented in Figure Figure22 suggests that most of the changes in background levels of synaptic activity affect the variability of *I*_{CaP} and *I*_{Kc}. Without synaptic activity the Purkinje cell model is quiescent, thus, the source of the variability observed in the dendritic currents comes from synaptic activity and from interactions among channels. In order to quantify the amount of variability due to *I*_{Glu} we made use of information theory (Victor, 2000; Nemenman et al., 2004; Paninski et al., 2004; Kennel et al., 2005). We initially calculated the entropy of *I*_{Glu}, *I*_{CaP} and *I*_{Kc} (Figure (Figure3A).3A). In order to quantify the causal relationship between synaptic activity and dendritic currents we performed the information analysis with a small time lag. Thus, the conditional entropy and mutual information were calculated between the dendritic currents at time *t* and the glutamate input at a time *t*−Δ*t*, we assumed Δ*t* to be 1ms. The analyses show that the entropy of dendritic currents depends on the combination of background activity. The entropy of *I*_{CaP} increases and reaches the same magnitude as the entropy of *I*_{Glu} as a function of the frequency of the background synaptic activity. Note that all of our calculations were biased corrected using the Panzeri and Treves method (Magri et al., 2009). The conditional entropy quantifies the uncertainty in the value of a variable given the known value of another. Calculating the conditional entropy of *I*_{CaP} and *I*_{Kc} with respect to *I*_{Glu} shows that the variability in the currents is high (Figure (Figure3B),3B), suggesting that most of the activity observed is due to intrinsic activity. However, the value of mutual information between *I*_{CaP} and *I*_{Kc} with respect to *I*_{Glu} shows an increase in information as a function of level of background level of stimulation (Figure (Figure3C).3C). The value of mutual information for the other calcium channel (*I*_{CaT}) was larger (**Figure S2** in Supplementary Material) than for *I*_{CaP} but since must of the current flows through *I*_{CaP} (Santamaria et al., 2002) we only analyzed this channel (see Discussion). Overall, the analysis presented in Figures Figures22 and and3;3; **Figure S2** in Supplementary Material show that most dendritic currents, except for CaP and K_{C}, are not sensitive to changes in background synaptic activity. The mutual information of Kc with respect to *I*_{Glu} is larger than for the same measurement done with CaP. However, the ratio of mutual information of these two currents changes as a function of the level of synaptic activity. More interestingly, this differential encoding of the background excitatory information can take place while the somatic firing rate remains constant.

The time constants of activation and inactivation of dendritic conductances imply that information content at one point in time is influenced by past synaptic activity. We characterized these time dependent changes by calculating the cross-correlation of *I*_{Glu} with *I*_{CaP} and *I*_{Kc} for all different combinations of excitatory and inhibitory input (Figure (Figure4).4). This analysis shows that as the excitatory input frequency increases the value of *I*_{CaP} (green) becomes more tightly coupled with *I*_{Glu} (blue, auto-correlation), while the coupling of the *I*_{Cap} and *I*_{Kc} also increases as we have previously reported using phase plane analysis (Santamaria et al., 2002). The amplitude of lag 0 is indicative of the amount of correlation with previous activity. The half-time auto-correlation time between *I*_{Glu}–*I*_{CaP} extracted from these plots ranges between 3–70ms.

We determined how the value of mutual information between *I*_{Glu} and *I*_{CaP} and *I*_{Kc} varied as an effect of previous activity. We calculated the mutual information between *I*_{CaP} and *I*_{Kc} with *I*_{Glu} for lag times up to 1s. Similar to what we did in Figure Figure3,3, at a given time lag of Δ*t* we computed the mutual information of the value of *I*_{CaP} or *I*_{Kc} at time *t* with the value of *I*_{Glu} at time *t*−Δ*t*. This analysis shows that the mutual information between *I*_{CaP} and *I*_{Glu} decays as a function of time lag (Figure (Figure5A).5A). This decay is due to the leaky nature of the conductance and membrane system. However, as opposed to the cross-correlation analysis the mutual information only shows one secondary peak at around 20ms. Similarly, the time lag analysis of *I*_{Glu} and *I*_{Kc} show an initial decay of the mutual information and a peak at around 20ms (Figure (Figure5B).5B). Channel kinetics and level of synaptic activity generate strong correlations at lag 0 which reflect the average level of activity in the dendrite, while the second peak is the mutual information in each channel from the instantaneous value of *I*_{Glu}. As is known from functional and anatomical studies of the cerebellar cortex, excitatory inputs to Purkinje cells are divided in two groups. The first group comes from parallel fibers, which are being studied in this work; the second group is the ascending segment synapses (Llinas, 1982). The ascending segment synapses deliver a more synchronous input that forces the Purkinje cell to fire robustly (Santamaria et al., 2007). The information analysis presented in Figure Figure55 suggests that paired stimulus stimulation from ascending segment synapses would then be influenced in a window of around 20ms by previous activity. This prediction is consistent with our previous analysis of ascending segment processing in the Purkinje cell model (Santamaria and Bower, 2005).

Plastic changes in the expression of active conductances affect the signal processing properties of dendritic trees (Segev and London, 1999; Stuart and Hausser, 2001). Each dendritic channel contributes to the encoding and transformation of synaptic input depending on channel kinetics and expression. Thus, we decided to quantify the dependence of the value of mutual information between *I*_{Glu} with *I*_{CaP} and *I*_{Kc} as a function of dendritic channel density.

We varied the excitability of the dendrite by changing the conductance of the CaP (*g*_{CaP}) or Kc (*g*_{Kc}). The value of *g*_{CaP} was varied from 25–70S/m^{2}, while keeping *g*_{Kc} at its control value (800S/m^{2}), and *g*_{Kc} ranged from 550–1000S/m^{2}, while keeping *g*_{CaP} at the original value of 45S/m^{2}. This wide range of conductances resulted in changes in the somatic firing rate (Figure (Figure6).6). In the case when we varied *g*_{CaP} the firing rate approached zero when *g*_{CaP}<40S/m^{2}. Varying *g*_{Kc} resulted in an almost linear dependence of the firing rate with the firing rate at each conductance level less than 10% from the collective average (Figure (Figure66B).

The information analysis of the conductance manipulations shown in Figure Figure66 result in that the entropy of *I*_{CaP} as a function of *g*_{CaP} decreases as *g*_{CaP} increased (Figure (Figure7A7A left). This dependence is robust to changes in the value of the firing rate (Figure (Figure6A).6A). The conditional entropy between *I*_{CaP} and *I*_{Glu} had a slightly different dependence of *g*_{CaP} (center). The combination of entropy and conditional entropy resulted in a relatively constant value of the mutual information that then increased as a function of background activity for increasing values of *g*_{CaP} (right). We repeated the same analysis to calculate information transmission of *I*_{Glu} through *I*_{Kc} as a function of *g*_{CaP} (Figure (Figure7B).7B). The functional dependence of the mutual information plots shown in Figures Figures7A,B7A,B showed a similar dependence. Note that the mutual information did not reflect the firing rate behavior. Similarly, we calculated mutual information from *I*_{Glu} to *I*_{CaP} and *I*_{Kc} as a function of *g*_{Kc} (Figures (Figures7C,D).7C,D). In these cases the mutual information increases as the *g*_{Kc} decreased, and the firing rate increased. This is due to the lower inhibitory effect of Kc over CaP. In practically all cases the ratio *N/m* was larger than 1 (**Figures S3 and S4** in Supplementary Material).

We further quantified the effect of time lags and dendritic excitability on information processing of *I*_{Glu} in the dendrite. We calculated the mutual information of *I*_{CaP} with *I*_{Glu} as a function of *g*_{CaP} and time lags up to 1s (Figure (Figure8A),8A), for all levels of background activity. This analysis shows that for high values of *g*_{CaP} and low levels of background activity mutual information is maximized (Figure (Figure8A8A right). The decay of the mutual information increases as a function of the value of *g*_{CaP}. A similar set of results is obtained when calculating mutual information content of *I*_{Glu}–*I*_{Kc} (Figure (Figure8B).8B). However, in this case, the effect of the background level of stimulation is not as strong as in *I*_{CaP}. A more subtle influence is shown by a similar analysis of the effects on the mutual information of *I*_{Glu} with *I*_{CaP} and *I*_{Kc} while varying *g*_{Kc} (Figures (Figures8C,D,8C,D, respectively). Overall, this analysis shows that the window of dendritic integration is directly modulated not only by the excitability provided by active conductances but also by the level of synaptic stimulation, especially in the case of *I*_{Kc}. Since both of these properties are dynamic and plastic, they allow the cell to regulate synaptic information processing.

Using a Purkinje cell model we presented a statistical analysis that shows that most of the information encoded in the background excitatory input is carried by the CaP and Kc currents. We have shown that the excitability of the dendritic tree represented by the density of active dendritic conductances modulates the information content of the total excitatory input in conjunction with the amount of synaptic activity. Consequently, the time window of integration of synaptic activity is in itself dependent on the level of synaptic activity which is not necessarily reflected in the spiking activity at the soma. This dynamical range of information processing consolidates the view that Purkinje cell dendrites processes synaptic activity locally before delivering it to the soma and provide a more sophisticated function than to compensate for the electrotonic decay of synaptic signals (Segev and London, 2000; Bekkers and Hausser, 2007).

Neuronal response to synaptic activity is an analog-to-digital transformation that requires the integration and processing of a large number of synaptic inputs distributed across a dendrite (Segev and London, 2000). Traditionally, information theory has been applied to determine how much of the input stimulus can be reconstructed from observing the spiking activity of a neuron (Rieke, 1997). Under this approach, the neuron is lumped into a point source function in which dendritic processing is lost (Harris and Stark, 1972; Sherry and Klemm, 1980; Yamada et al., 1996; Borst and Theunissen, 1999; Levy and Baxter, 2002; Schultz et al., 2009). In this study we refined the use of mutual information and linked it to a biophysical property of dendrites. In this framework, each conductance is treated as an *information channel* that communicates and interacts with other information channels in the dendrite. We found that the mutual information between excitatory synaptic input and dendritic currents varies as a function of dendritic excitability and the level of synaptic activity.

Although the activity of dendritic currents is dependent on *I*_{Glu} stimulation in the past our analysis shows strong mutual information at lag 0. This ‘background’ information stems from the overall level of activity in the dendrite as also shown by the cross-correlation analysis. These results are dependent on the properties of the stimulus that were meant to replicate the continuous activity arriving from parallel fibers and inhibitory synapses. This stimulation resulted in an overall level of activity of dendritic channels that is reflected at zero lag time. However, the second peak shown in the cross-correlation and mutual information is more indicative of the internal dynamics of the dendritic currents and the information content of the instantaneous value of *I*_{Glu}. The differential encoding of synaptic activity underscores the modulatory property of the parallel fiber system on Purkinje cell activity (Santamaria et al., 2007; Walter et al., 2009).

Purkinje cells also receive synaptic activity from the ascending segment part of the granule cell axon and climbing fibers (Llinas, 1982). Both of these synaptic inputs deliver a more instantaneous stimulus lasting a few milliseconds. Under such types of stimuli and in the absence of background parallel fiber stimulation the instantaneous mutual information would be zero and would increase to reach a maximum, possibly at 20ms, as indicated by our cross-correlation and mutual information analyses. Thus, mutual information from different systems could be encoded differently in the Purkinje cell dendritic tree.

The time window in which we define our stimulus affects the value of the mutual information. For example, our window of analysis (1ms) would be too slow to study information at the some where fast channels are present. However, other types of cells could have fast dendritic channels (Regehr and Tank, 1992). Thus, caution should be taken when applying information theory to analyzing dendritic activity.

Bounding the kinetic and density parameters of dendritic current models has been difficult since present experimental techniques are limited in resolving the spatial distribution of dendritic conductances (Achard and De Schutter, 2008; Van Geit et al., 2008). We propose that by treating each dendritic conductance as an information channel it is possible to constrain the value and distribution of dendritic currents in encoding synaptic information. Although here we have presented data regarding encoding the total excitatory current, the same technique can be extended to total inhibitory, total synaptic current, or it can be subdivided into the different parts of the dendritic tree (primary, secondary, and tertiary dendrites). At a first iteration, traditional template matching algorithms could be used to fit a model to data. These algorithms change the density of dendritic conductances (Achard and De Schutter, 2008). Then models could be classified due to their value of mutual information between the stimulus and each dendritic current. Using this approach we could test the hypothesis that by matching the spiking activity at the soma the information content of the synaptic activity in dendritic currents is maximized.

The simulation analysis presented in this work can be used to design experiments based on dynamic clamp technique (Sharp et al., 1993; Kreiner and Jaeger, 2004). A challenging experiment would consist in patching a Purkinje cell while blocking spiking activity with tetrodotoxin. The pipette could deliver artificial synaptic conductances while monitoring the resulting current in the soma. The mutual information between the synaptic input and the total dendritic current value could be calculated from this experiment. Pharmacological manipulation blocking dendritic calcium or calcium activated potassium channels in the dendrite would then uncover the contribution to information processing. Other similar experiments could be to electrically or photochemically stimulate a group of synapses while recording at the soma (Korkotian and Segal, 2006) or to use the recently developed targeted dendrotomy technique (Bekkers and Hausser, 2007).

Dendritic excitability, in Purkinje and other cell types, is the consequence of the different concentration and distribution of active conductances over a dendritic tree (Gunay et al., 2008; Komendantov and Ascoli, 2009). Neuronal stimulation and development change the excitability of the cell and could be a mechanism for storing memories (Womack and Khodakhah, 2003; Kim and Linden, 2007). Similarly, excitability is affected by aging and illnesses that could reduce the information capacity of the dendrite (Landfield and Pitler, 1984; Chan et al., 2007; Kabaso et al., 2009). It is well established that hippocampal pyramidal cells firing rate remains constant through aging (Barnes et al., 1983; Wilson et al., 2005) but not in Purkinje cells (Rogers et al., 1980), however the contribution of dendritic excitability and afferent stimulation on these changes is not well understood. In any case, our analyses have shown that changes in *g*_{CaP} or *g*_{Kc} can regulate the amount of information encoded and the time synaptic activity affects neuronal activity. As the hyperpolarizing Kc current decreases, the information contained in *I*_{CaP} and *I*_{Kc} increases potentially allowing more information to be transmitted by the neuron (Steuber et al., 2007).

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The Supplementary Material for this article can be found online at http://www.frontiersin.org/computationalneuroscience/paper/10.3389/fncom.2010.00006/

UTSA-TRAC and NSF HRD-0932339.

- Achard P., De Schutter E. (2008). Calcium, synaptic plasticity and intrinsic homeostasis in Purkinje neuron models. Front Comput Neurosci 2:8.10.3389/neuro.10.008.2008 [PMC free article] [PubMed] [Cross Ref]
- Barnes C. A., McNaughton B. L., O'Keefe J. (1983). Loss of place specificity in hippocampal complex spike cells of senescent rat. Neurobiol. Aging 4, 113–11910.1016/0197-4580(83)90034-9 [PubMed] [Cross Ref]
- Bekkers J. M., Hausser M. (2007). Targeted dendrotomy reveals active and passive contributions of the dendritic tree to synaptic integration and neuronal output. Proc. Natl. Acad. Sci. U.S.A. 104, 11447–1145210.1073/pnas.0701586104 [PubMed] [Cross Ref]
- Borst A., Theunissen F. E. (1999). Information theory and neural coding. Nat. Neurosci. 2, 947–95710.1038/14731 [PubMed] [Cross Ref]
- Bower J. M. (2002). The organization of cerebellar cortical circuitry revisited: implications for function. Ann. N. Y. Acad. Sci. 978, 135–15510.1111/j.1749-6632.2002.tb07562.x [PubMed] [Cross Ref]
- Chan C. S., Guzman J. N., Ilijic E., Mercer J. N., Rick C., Tkatch T., Meredith G. E., Surmeier D. J. (2007). /‘Rejuvenation/’ protects neurons in mouse models of Parkinson/’s disease. Nature 447, 1081–108610.1038/nature05865 [PubMed] [Cross Ref]
- De Schutter E., Bower J. M. (1994a). Simulated responses of cerebellar Purkinje cells are independent of the dendritic location of granule cell synaptic inputs. Proc. Natl. Acad. Sci. U.S.A. 91, 4736–474010.1073/pnas.91.11.4736 [PubMed] [Cross Ref]
- De Schutter E., Bower J. M. (1994b). An active membrane model of the cerebellar Purkinje cell II. Simulation of synaptic responses. J. Neurophysiol. 71, 401–419 [PubMed]
- De Schutter E., Bower J. M. (1994c). An active membrane model of the cerebellar Purkinje cell. I. Simulation of current clamps in slice. J. Neurophysiol. 71, 375–400 [PubMed]
- De Schutter E., Steuber V. (2009). Patterns and pauses in Purkinje cell simple spike trains: experiments, modeling and theory. Neuroscience 162, 816–82610.1016/j.neuroscience.2009.02.040 [PubMed] [Cross Ref]
- Foffani G., Priori A. (2007). Information theory, single neurons and gamma oscillations in the human subthalamic nucleus. Exp. Neurol. 205, 292–29310.1016/j.expneurol.2007.02.013 [PubMed] [Cross Ref]
- French C. R., Sah P., Buckett K. J., Gage P. W. (1990). A voltage-dependent persistent sodium current in mammalian hippocampal neurons. J. Gen. Physiol. 95, 1139–115710.1085/jgp.95.6.1139 [PMC free article] [PubMed] [Cross Ref]
- Gollo L. L., Kinouchi O., Copelli M. (2009). Active dendrites enhance neuronal dynamic range. PLoS Comput. Biol. 5, e1000402.10.1371/journal.pcbi.1000402 [PMC free article] [PubMed] [Cross Ref]
- Gruol D. L., Jacquin T., Yool A. J. (1991). Single-channel K+ currents recorded from the somatic and dencritic regions of cerebellar Purkinje neurons in culture. J. Neurosci. 11, 1002–1015 [PubMed]
- Gunay C., Edgerton J. R., Jaeger D. (2008). Channel density distributions explain spiking variability in the globus pallidus: a combined physiology and computer simulation database approach. J. Neurosci. 28, 7476–749110.1523/JNEUROSCI.4198-07.2008 [PubMed] [Cross Ref]
- Gundappa-Sulur G., De Schutter E., Bower J. M. (1999). Ascending granule cell axon: an important component of cerebellar cortical circuitry. J. Comp. Neurol. 408, 580–59610.1002/(SICI)1096-9861(19990614)408:4<580::AID-CNE11>3.0.CO;2-O [PubMed] [Cross Ref]
- Harris D. A., Stark L. (1972). Information theoretical analysis of a crayfish behavioral system. Brain Res. 36, 425–43010.1016/0006-8993(72)90749-4 [PubMed] [Cross Ref]
- Hirano T., Hagiwara S. (1989). Kinetics and distribution of voltage-gated Ca, Na and K channels on the somata of rat cerebellar Purkinje cells. Pflugers Arch. Eur. J. Physiol. 413, 463–46910.1007/BF00594174 [PubMed] [Cross Ref]
- Hirst G. D., McLachlan E. M. (1986). Development of dendritic calcium currents in ganglion cells of the rat lower lumbar sympathetic chain. J. Physiol. (Lond.) 377, 349–368 [PubMed]
- Hockberger P. E., Tseng H. Y., Connor J. A. (1989). Fura-2 measurements of cultured rat Purkinje neurons show dendritic localization of Ca2+ influx. J. Neurosci. 9, 2272–2284 [PubMed]
- Huguenard J. R., Hamill O. P., Prince D. A. (1989). Sodium channels in dendrites of rat cortical pyramidal neurons. Proc. Natl. Acad. Sci. U.S.A. 86, 2473–247710.1073/pnas.86.7.2473 [PubMed] [Cross Ref]
- Jaeger D., Bower J. M. (1994). Prolonged responses in rat cerebellar Purkinje cells following activation of the granule cell layer: an intracellular
*in vitro*and*in vivo*investigation. Exp. Brain Res. 100, 200–21410.1007/BF00227191 [PubMed] [Cross Ref] - Jaeger D., Bower J. M. (1999). Synaptic control of spiking in cerebellar Purkinje cells: dynamic current clamp based on model conductances. J. Neurosci. 19, 6090–6101 [PubMed]
- Jaeger D., De Schutter E., Bower J. M. (1997). The role of synaptic and voltage-gated currents in the control of Purkinje cell spiking: a modeling study. J. Neurosci. 17, 91–106 [PubMed]
- Kabaso D., Coskren P. J., Henry B. I., Hof P. R., Wearne S. L. (2009). The electrotonic structure of pyramidal neurons contributing to prefrontal cortical circuits in macaque monkeys is significantly altered in aging. Cereb. Cortex. 19, 2248–226810.1093/cercor/bhn242 [PMC free article] [PubMed] [Cross Ref]
- Kaneda M., Wakamori M., Ito C., Akaike N. (1990). Low-threshold calcium current in isolated Purkinje cell bodies of rat cerebellum. J. Neurophysiol. 63, 1046–1051 [PubMed]
- Kennel M. B., Shlens J., Abarbanel H. D., Chichilnisky E. J. (2005). Estimating entropy rates with Bayesian confidence intervals. Neural. Comput. 17, 1531–157610.1162/0899766053723050 [PubMed] [Cross Ref]
- Kim S. J., Linden D. J. (2007). Ubiquitous plasticity and memory storage. Neuron 56, 582–59210.1016/j.neuron.2007.10.030 [PubMed] [Cross Ref]
- Komendantov A. O., Ascoli G. A. (2009). Dendritic excitability and neuronal morphology as determinants of synaptic efficacy. J. Neurophysiol. 101, 1847–186610.1152/jn.01235.2007 [PubMed] [Cross Ref]
- Korkotian E., Segal M. (2006). Spatially confined diffusion of calcium in dendrites of hippocampal neurons revealed by flash photolysis of caged calcium. Cell Calcium 40, 441–44910.1016/j.ceca.2006.08.008 [PubMed] [Cross Ref]
- Kreiner L., Jaeger D. (2004). Synaptic shunting by a baseline of synaptic conductances modulates responses to inhibitory input volleys in cerebellar Purkinje cells. Cerebellum 3, 112–12510.1080/14734220410031990 [PubMed] [Cross Ref]
- Landfield P., Pitler T. (1984). Prolonged Ca2+-dependent afterhyperpolarizations in hippocampal neurons of aged rats. Science 226, 1089–109210.1126/science.6494926 [PubMed] [Cross Ref]
- Latorre R., OBerhauser A., Labarca P., Alvarez O. (1989). Varieties of calcium-activated potassium channels. Annu. Rev. Physiol. 51, 385–39910.1146/annurev.ph.51.030189.002125 [PubMed] [Cross Ref]
- Levy W. B., Baxter R. A. (2002). Energy-efficient neuronal computation via quantal synaptic failures. J. Neurosci. 22, 4746–4755 [PubMed]
- Llinas R. (1982). “Radial connectivity in the cerebellar cortex: a novel view regarding the functional organization of the molecular layer,” in The cerebellum: new vistas, eds Plalay S., Chan-Palay C., editors. (New York: Elsevier; ), 189–184
- Llinas R. R., Sugimori M., Cherksey B. (1989). Voltage-dependent calcium conductances in mammalian neurons. The P channel. Ann. N.Y. Acad. Sci. 560, 103–11110.1111/j.1749-6632.1989.tb24084.x [PubMed] [Cross Ref]
- Llinas R., Yarom Y. (1981). Electrophysiology of mammalian inferior olivary neurones
*in vitro*. Different types of voltage-dependent ionic conductances. J. Physiol. 315, 549–567 [PubMed] - London M., Schreibman A., Hausser M., Larkum M. E., Segev I. (2002). The information efficacy of a synapse. Nat. Neurosci. 5, 332–34010.1038/nn826 [PubMed] [Cross Ref]
- Maffezzoli A., Signorini M. G., Gullo F., Wanke E. (2008). Modeling spiking activity of
*in vitro*neuronal networks through non linear methods. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2008, 42–45 [PubMed] - Magri C., Whittingstall K., Singh V., Logothetis N. K., Panzeri S. (2009). A toolbox for the fast information analysis of multiple-site LFP, EEG and spike train recordings. BMC Neurosci. 10, 81.10.1186/1471-2202-10-81 [PMC free article] [PubMed] [Cross Ref]
- Marr D. (1969). A theory of cerebellar cortex. J. Physiol. (Lond.) 202, 437–470 [PubMed]
- Masukawa L. M., Hansen A. J., Shepherd G. (1991). Distribution of single-channel conductances in cultured rat hippocampal neurons. Cell. Mol. Neurobiol. 11, 231–24310.1007/BF00769036 [PubMed] [Cross Ref]
- Migliore M., Shepherd G. M. (2002). Emerging rules for the distributions of active dendritic conductances. Nat. Rev. Neurosci. 3, 362–37010.1038/nrn810 [PubMed] [Cross Ref]
- Nemenman I., Bialek W., de Ruyter van Steveninck R. (2004). Entropy and information in neural spike trains: progress on the sampling problem. Phys. Rev. E Stat. Nonlin. Soft Matter Phys. 69, 056111.10.1103/PhysRevE.69.056111 [PubMed] [Cross Ref]
- Paninski L., Pillow J. W., Simoncelli E. P. (2004). Maximum likelihood estimation of a stochastic integrate-and-fire neural encoding model. Neural Comput. 16, 2533–256110.1162/0899766042321797 [PubMed] [Cross Ref]
- Panzeri S., Senatore R., Montemurro M. A., Petersen R. S. (2007). Correcting for the sampling bias problem in spike train information measures. J. Neurophysiol. 98, 1064–107210.1152/jn.00559.2007 [PubMed] [Cross Ref]
- Pellionisz A., Szentagothai J. (1973). Dynamic single unit simulation of a realistic cerebellar network model. Brain Res. 49, 83–9910.1016/0006-8993(73)90403-4 [PubMed] [Cross Ref]
- Pellionisz A., Szentagothai J. (1974). Dynamic single unit simulation of a realistic cerebellar network model. II. Purkinje cell activity within the basic circuit and modified by inhibitory systems. Brain Res. 68, 19–4010.1016/0006-8993(74)90531-9 [PubMed] [Cross Ref]
- Poirazi P., Mel B. W. (2001). Impact of active dendrites and structural plasticity on the memory capacity of neural tissue. Neuron 29, 779–79610.1016/S0896-6273(01)00252-5 [PubMed] [Cross Ref]
- Quian Quiroga R., Panzeri S. (2009). Extracting information from neuronal populations: information theory and decoding approaches. Nat. Rev. Neurosci. 10, 173–18510.1038/nrn2578 [PubMed] [Cross Ref]
- Rall W. (1962a). Electrophysiology of a dendritic neuron model. Biophys. J. 2, 145–16710.1016/S0006-3495(62)86953-7 [PubMed] [Cross Ref]
- Rall W. (1962b). Theory of physiological properties of dendrites. Ann. N. Y. Acad. Sci. 96, 1071–109210.1111/j.1749-6632.1962.tb54120.x [PubMed] [Cross Ref]
- Regehr W. G., Tank D. W. (1992). Calcium concentration dynamics produced by synaptic activation of CA1 hippocampal pyramidal cells. J. Neurosci. 12, 4202–4223 [PubMed]
- Rieke F. (1997). Spikes: Exploring the Neural Code Cambridge, MA: MIT Press; [PubMed]
- Rogers J., Silver M. A., Shoemaker W. J., Bloom F. E. (1980). Senescent changes in a neurobiological model system: cerebellar Purkinje cell electrophysiology and correlative anatomy. Neurobiol. Aging 1, 3–1110.1016/0197-4580(80)90018-4 [PubMed] [Cross Ref]
- Santamaria F., Bower J. M. (2005). Background synaptic activity modulates the response of a modeled Purkinje cell to paired afferent input. J. Neurophysiol. 93, 237–25010.1152/jn.00458.2004 [PubMed] [Cross Ref]
- Santamaria F., Jaeger D., De Schutter E., Bower J. M. (2002). Modulatory effects of parallel fiber and molecular layer interneuron synaptic activity on Purkinje cell responses to ascending segment input: a modeling study. J. Comput. Neurosci. 13, 217–23510.1023/A:1020266315730 [PubMed] [Cross Ref]
- Santamaria F., Tripp P. G., Bower J. M. (2007). Feedforward inhibition controls the spread of granule cell-induced Purkinje cell activity in the cerebellar cortex. J. Neurophysiol. 97, 248–26310.1152/jn.01098.2005 [PubMed] [Cross Ref]
- Schultz S. R., Kitamura K., Post-Uiterweer A., Krupic J., Hausser M. (2009). Spatial pattern coding of sensory information by climbing fiber-evoked calcium signals in networks of neighboring cerebellar Purkinje cells. J. Neurosci. 29, 8005–801510.1523/JNEUROSCI.4919-08.2009 [PubMed] [Cross Ref]
- Segev I., London M. (1999). “A theoretical view of passive and active dendrites,” in Dendrites, eds Stuart G., Nelson S., Hausser M., editors. (Oxford: Oxford university press; ), 205–230
- Segev I., London M. (2000). Untangling dendrites with quantitative models. Science 290, 744–75010.1126/science.290.5492.744 [PubMed] [Cross Ref]
- Shannon C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423, 623–656
- Shannon C. E., Weaver W. (1949). The Mathematical Theory of Communication Urbana: University of Illinois Press
- Sharp A. A., Oneil M. B., Abbott L. F., Marder E. (1993). The dynamic clamp – artificial conductances in biological neurons. Trends Neurosci. 16, 389–39410.1016/0166-2236(93)90004-6 [PubMed] [Cross Ref]
- Sherry C. J., Klemm W. R. (1980). The statistical relationship between the “entropy” of a neuronal signal and its variability. Int. J. Neurosci. 11, 109–11310.3109/00207458009150333 [PubMed] [Cross Ref]
- Spain W. J., Schwindt P. C., Crill W. E. (1987). Anomalous rectification in neurons from cat sensorimotor cortex
*in vitro*. J. Neurophysiol. 57, 1555–1576 [PubMed] - Stemmler M., Koch C. (1999). How voltage-dependent conductances can adapt to maximize the information encoded by neuronal firing rate. Nat. Neurosci. 2, 521–52710.1038/9173 [PubMed] [Cross Ref]
- Steuber V., Mittmann W., Hoebeek F. E., Silver R. A., De Zeeuw C. I., Häusser M., De Schutter E. (2007). Cerebellar LTD and pattern recognition by Purkinje cells. Neuron 54, 121–13610.1016/j.neuron.2007.03.015 [PMC free article] [PubMed] [Cross Ref]
- Stuart G. J., Hausser M. (2001). Dendritic coincidence detection of EPSPs and action potentials. Nat. Neurosci. 4, 63–7110.1038/82910 [PubMed] [Cross Ref]
- Suter K. J., Jaeger D. (2004). Reliable control of spike rate and spike timing by rapid input transients in cerebellar stellate cells. Neuroscience 124, 305–31710.1016/j.neuroscience.2003.11.015 [PubMed] [Cross Ref]
- Van Geit W., De Schutter E., Achard P. (2008). Automated neuron model optimization techniques: a review. Biol. Cybern. 99, 241–25110.1007/s00422-008-0257-6 [PubMed] [Cross Ref]
- Vetter P., Roth A., Hausser M. (2001). Propagation of action potentials in dendrites depends on dendritic morphology. J. Neurophysiol. 85, 926–937 [PubMed]
- Victor J. D. (2000). Asymptotic bias in information estimates and the exponential (Bell) polynomials. Neural. Comput. 12, 2797–280410.1162/089976600300014728 [PubMed] [Cross Ref]
- Walter J. T., Dizon M.-J., Khodakhah K. (2009). The functional equivalence of ascending and parallel fiber inputs in cerebellar computation. J. Neurosci. 29, 8462–847310.1523/JNEUROSCI.5718-08.2009 [PubMed] [Cross Ref]
- Wilson I. A., Ikonen S., Gallagher M., Eichenbaum H., Tanila H. (2005). Age-associated alterations of hippocampal place cells are subregion specific. J. Neurosci. 25, 6877–688610.1523/JNEUROSCI.1744-05.2005 [PubMed] [Cross Ref]
- Womack M. D., Khodakhah K. (2003). Somatic and dendritic small-conductance calcium-activated potassium channels regulate the output of cerebellar Purkinje neurons. J. Neurosci. 23, 2600–2607 [PubMed]
- Xu J., Kang J. (2005). The mechanisms and functions of activity-dependent long-term potentiation of intrinsic excitability. Rev. Neurosci. 16, 311–323 [PubMed]
- Yamada S., Matsumoto K., Nakashima M., Shiono S. (1996). Information theoretic analysis of action potential trains. II. Analysis of correlation among
*n*neurons to deduce connection structure. J. Neurosci. Methods 66, 35–4510.1016/0165-0270(95)00152-2 [PubMed] [Cross Ref] - Yamada W. M., Koch C., Adams P. R. (1989). “Multiple channels and calcium dynamics,” in Methods in Neuronal Modeling: From Synapses to Networks, eds Koch C., Segev I., editors. (Cambridge: MIT Press; ), 97–133 [PubMed]

Articles from Frontiers in Computational Neuroscience are provided here courtesy of **Frontiers Media SA**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |