- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Computational Intelligence and Neuroscience
Volume 2012 (2012), Article ID 359529, 17 pages
Modeling Spike-Train Processing in the Cerebellum Granular Layer and Changes in Plasticity Reveal Single Neuron Effects in Neural Ensembles
1Amrita School of Biotechnology, Amrita Vishwa Vidyapeetham (Amrita University), Amritapuri, Clappana, Kollam 690525, Kerala, India
2Department of Physiology, University of Pavia, Via Forlanini 7, 21000 Pavia, Italy
3Brain Connectivity Center, IRCCS C. Mondino, Via Mondino 2, 27100 Pavia, Italy
4Department of Mathematics, University of Milan, Via Saldini 50, 27100 Milan, Italy
Received 6 March 2012; Accepted 12 July 2012
Academic Editor: Steven Bressler
Copyright © 2012 Chaitanya Medini et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
The cerebellum input stage has been known to perform combinatorial operations on input signals. In this paper, two types of mathematical models were used to reproduce the role of feed-forward inhibition and computation in the granular layer microcircuitry to investigate spike train processing. A simple spiking model and a biophysically-detailed model of the network were used to study signal recoding in the granular layer and to test observations like center-surround organization and time-window hypothesis in addition to effects of induced plasticity. Simulations suggest that simple neuron models may be used to abstract timing phenomenon in large networks, however detailed models were needed to reconstruct population coding via evoked local field potentials (LFP) and for simulating changes in synaptic plasticity. Our results also indicated that spatio-temporal code of the granular network is mainly controlled by the feed-forward inhibition from the Golgi cell synapses. Spike amplitude and total number of spikes were modulated by LTP and LTD. Reconstructing granular layer evoked-LFP suggests that granular layer propagates the nonlinearities of individual neurons. Simulations indicate that granular layer network operates a robust population code for a wide range of intervals, controlled by the Golgi cell inhibition and is regulated by the post-synaptic excitability.
Decoding neural activity is the key to understand spatiotemporal patterns that the brain receives as sensory information regarding the world. Time-scale of operation is closely correlated to the activity of the neural circuit and decoding such activity reveals principles regarding the function. One of the main circuits in the cerebellum is the large input layer circuit formed of granule and Golgi cells. Spatiotemporal information is one of the unique functional characteristics observed in the cerebellar input layer network [1, 2]. Cerebellar granular layer forms the input stage of the cerebellum in which information coming from the peripheral and central systems converge through the mossy fibers. The granular layer has by far the smallest (~5 μm) and the most numerous neurons (~1011) in humans. Understanding how the granular layer process information appears critical to understand the cerebellar function, since signals coming into upper cortical layers are provided by the granular layer. The granule cells form the largest neuronal population in the mammalian brain and regulate information transfer along the major afferent systems to the cerebellum. The granule layer receives excitatory input primarily from mossy fibers and inhibition via synapses from interneurons like Golgi cell. Mossy fiber input excites both the granule cell and inhibitory interneurons like the Golgi cell. The granule cell is a small neuron with 3–5 dendrites. Timing in the cerebellar granular layer plays a key role via passage-of-time representation (POT), learning or adaptation to movements [3, 4], modulation of information transfer to Purkinje cells (activation of granule cell subsets with respect to time). Knockout and lesion studies revealed that disruption of the granular layer leads to abnormal functioning of the cerebellar mossy fiber-granule cell relay , affects the learning-dependent timing of conditioning eyelid responses , loss of rapid spike processing in the cerebellum (results in ataxia)  which in turn affects the plasticity of granule-Purkinje cell synapses. Prion protein (PrP) knockout mice showed a large proportion of granule cells (~40%) with slow non-overshooting nonrepetitive action potential, slow EPSPs, and no inward rectification . Likewise, FHF1-FHF4 mutants showed impaired granule cell excitability which prevented rapid burst transmission in the cerebellum .
The objective of the current paper is to study how excitation operates in the granular layer network [1, 2, 8–11]. The focus is also on understanding the modulatory role of inhibition [1, 2] in the granule cells  and underlying ensemble activity in terms of combinatorial operations on granular layer network .
In order to estimate spiking behaviour and to test reliability in modeling, we tried using detailed and simple models of neurons in our network models. The use of detailed multicompartmental models was focused towards reproducing the spatiotemporal dynamics of normal cerebellar activity. The detailed models allow reconstructing information including local field potentials , which were not seen while using less detailed models. We also used a spiking neuron model (modified from ) for reconstructing network activity in order to understand the contributions of individual spikes in the cerebellar cortex. The necessity of these spiking models was to retain the biologically plausibility of Hodgkin-Huxley-type dynamics, while maintaining the low computational cost. It was also an attempt to validate whether such simple spiking models also allow us to create computationally simpler yet large-scale models of cerebellum. Using properties of the granule cell  and Golgi neuron [15, 16], we developed simple spiking models  to represent the spiking behavior in a network. Estimates of spiking and reproducibility of spiking could be very useful for computationally efficient and large network models. Spike modulation due to the effect of feed-forward inhibition has been known to play an important role in time-windows hypothesis . The simulations quantify how spikes pass through the granular layer network and the role of feed-forward inhibition in the neuronal microcircuit. This gives a new paradigm on the functional relevance of patterns in the cerebellar granular layer circuitry. The main objectives of the paper were also to reconstruct the center-surround excitation patterns and observe role of combined inhibition and excitation geometry in frequency-dependent transmission of spike information.
The objective was also to understand the effect of combinatorial operations on the granular layer network . Combinatorial operations included combined excitation and inhibition which forms the spatiotemporal pattern in granular layer network in vitro and in vivo. Another objective behind the simulations was to understand information flow in granule neurons via burst-burst transmission and feed-forward inhibition . Together they suggest the role of granule cells in expansion recoding and sparse activation via mossy fiber-granule cell synapses. The paper reports potential reconstruction of network activity in the form of center-surround structures , spike properties of underlying cells, and modulation of spikes due to changes induced by synaptic plasticity and due to inhibition. A special case of NMDA receptor blocking was simulated since GABA (γ-aminobutyric acid)-ergic inhibition is especially effective in controlling NMDA (N-methyl-D-aspartate) receptor-dependent depolarization in the granular layer . Network simulations predict specific computational roles of granule cells in processing bursts and overall spike processing in the cerebellar granular layer.
Understanding population code via comparisons of spatiotemporal properties of simulated neural activity and with experimental measurements using multielectrode recordings  is useful to identify how information encoding happens in microcircuits. Therefore, evoked local field potential (LFP) responses from granular layer in vitro  were reconstructed computationally . The main intention of local field reconstruction on the network in vitro was to study the role of inhibition in generating the wave [4, 20]. The paper also investigates the combined role of excitation-inhibition affecting the granular layer clusters.
The study carried out in this paper involved the use of computational models of neurons based on experimental data from p17–23 Wistar rat cerebellum . Mathematical neuron models of granule cell (GrC) [9, 19] and Golgi cell model (GoC) [15, 16] were used in this network study. Modeling reliability for spiking models was based on the extensive characterization of membrane currents and the compact electrotonic structure of cerebellar granule cells [9, 19]. The models used AMPA (2-amino-3-(5-methyl-3-oxo-1,2-oxazol-4-yl) propanoic acid) and NMDA receptor components as excitatory mossy fiber (MF)-GrC synapses and GABAergic synapses for the GoC-GrC relay [1, 15, 19, 21]. On an average, each granule cell receives excitatory connections from 4-5 mossy fibers .
2.1. Simple Neuronal Models
The objective of using simple models in the study was to understand how spatiotemporal patterns integrated over time to produce responses that are selective to specific patterns and to reconstruct the representation of spiking behavior in networks especially to study how inhibition affects intrinsic electroresponsiveness. We also wanted to see how the high variability in MF spike trains affects firing behavior. A simple spiking model  was used to study the neuronal spiking activity. A good model should be feasible with Hodgkin-Huxley dynamics and be computationally efficient  to reproduce the firing behavior of biorealistic model. The simple spiking model  of neuron primarily used two equations one regulating the membrane potential () and the other regulating adaptation current () where is the membrane capacitance, gL is the leak conductance, is the voltage threshold, is the slope factor which quantifies the sharpness of the spikes, is an adaptation current, and is the injected current. The membrane time constant is The reset values were activated when the membrane potential reached the desired peak voltage: The change in parameters (as shown in (4)) allowed the simple spiking model to replicate amplitude of granule cell firing and Golgi cell firing (see Table 1). Capacitance values of Golgi and granule cells were modelled to simulate similar frequency and firing patterns matching experimentally observed data, as observed in biophysical models [9, 15].
While replicating the basic firing behavior, the Adaptive Exponential Integrate-and-fire (AdEx, see ) granule cell model was tested with various combinations of synaptic connections. AMPA [17, 22, 23] synaptic kinetics was used as the excitatory synaptic dynamics and GABA [17, 24, 25] as the inhibitory synaptic kinetics as observed via experiments : Likewise, GABA synaptic kinetics was modeled using the GABA-A equation : The maximal conductance of AMPA and GABA was adjusted to suit the observed biophysical firing pattern (see Table 2). For simulating the case without inhibition, maximal conductance value of GABA synapses is set to zero (see Table 2). To model inhibition, varying values of maximal conductance for different number of inhibitory synapses were used. With these two receptor kinetics, we were able to match the number of spikes and amplitude to that of biophysical model. The models also reproduced a close correlation to the biophysical model while varying the intrinsic excitability and release probability (data not shown). By changing conductance and dynamic parameters, we simulated LTP and LTD, matching to experimental and computational values .
To understand network dynamics in terms of firing and temporal processing, we used a simple spiking network with 1680 AdEx granule cell models and 1 AdEx Golgi cell model. The properties of the network model were matched (see Tables 1 and 2 and [9, 15, 16, 27, 28]) to experimental data.
2.2. Multicompartmental Models
A detailed multicompartmental GrC model  was used, and simulations were performed by varying the excitatory () and inhibitory () synaptic inputs. The model of the granule cell was based on multicompartmental cable theory and included soma, axon, hillock and dendritic compartments. The model consisted of 52 active compartments connected to each other via the -power law . For each of the compartments, membrane voltage had to be estimated separately:
where is the conductance corresponding to (ion channel), syn (synaptic dynamics), br (neighboring attached branch), and tot (total). Here, which is the time constant of oscillation of the membrane based on its membrance resistance, and membrane capacitance, . The calcium current in the model was included as where is the depth of a shell adjacent to the cell surface of area , determines the loss of calcium ions from the shell approximating the effect of fluxes, ionic pumps, diffusion, and buffers, is resting [Ca], and is the Faraday’s constant. [Ca] is the calcium channel dynamics as reported in .
The model GrC has 1–4 excitatory (one for each dendrite) and 0 (no inhibition)–4 inhibitory connections (one for each dendrite) . The detailed explanations of ionic channel dynamics, compartmental localization of ion channels and electronic structure of this granule neuron model are described elsewhere [9, 19, 21]. Since granule cell is one of the rarest neurons where the ionic channel densities can be accurately determined using whole-cell patch clamp, the ion channel dynamics that was modeled previously [7, 9, 21] is not repeated here. Also, excitatory and inhibitory synaptic inputs to the dendrites were located in dendritic tips although in neighboring dendritic compartments. Presynaptic dynamics for the MF-GrC was modeled separately as in [19, 21] due to components such as facilitation and depression. Excitatory postsynaptic mechanisms were shown as AMPA and NMDA postsynaptic receptor components as seen in granule neurons. AMPA receptor dynamics was modeled using a three-state scheme and a 2D diffusion model, whereas the NMDA receptors used Boltzmann equation as seen in . Both the excitatory presynaptic and excitatory postsynaptic mechanisms are described in detail elsewhere . The GoC-GrC inhibitory synapse model was based on the following presynaptic dynamics: release probability = 0.35, = 36 ms, = 58.5 ms, and = 0.1 ms, respectively and as described in . Effects of blocking inhibition by adding gabazine were also simulated by setting GABAergic conductance in inhibitory fibers to zero.
The Golgi neuron was adapted from [15, 16]. All simulations were performed with NEURON environment  running on HP Blade C3000 node. Timing and initial time-window modulations are mainly affected by the role of feed-forward inhibition as it happens with only a slight delay from the mossy-fiber input and hence the role of feed-back inhibition was not simulated.
2.3. Granular Layer Network
Granular layer spiking network model consisted of 140 homogenous mossy fibers (MF) rosettes, 1680 granule cells (GrC), and 1 Golgi cell (GoC). In this network, about 48 GrC receive 1 excitatory input from the same mossy fiber, and each granule cell receives four excitatory connections from four different mossy fibers. Along with these excitatory inputs given to GrC, mossy fibers also provide excitatory input to GoC whose ratio was set in this model to about 78 : 1 (see [15, 16], each dendrite had 26 synapses in the GoC model and assuming a total of 3 dendrites, we approximated to 78 synapses) providing an overall glomeruli connectivity pattern . The network topology is illustrated in Figure 1(a).
Modeling responses in brain slices in vitro were simulated by giving single spike as input via mossy fiber (MF) terminals. Anaesthetized rat brain recordings in vivo showed bursts as inputs through mossy fibers . Therefore, in vivo inputs to GrC were simulated as bursts of (5 or 9) spikes via the MF input.
2.4. Center-Surround “Spot” Pattern
Stimulating mossy fibers with an electrode at a particular point activates granule cells in the network in a center-surround activation pattern . Within a “spot,” cells which are in close proximity to the electrode will receive high excitation and the periphery layer cells receive less excitation. In the network model, we simulated the center-surround pattern (see Figure 4A in ) defined as a “spot,” showing decreasing strengths of excitation spreading from the center to the periphery. In each spot, 48 cells in the center received 4 mossy fiber (MF) inputs, 144 cells received 3 MF inputs, 48 cells had 2 MF inputs, and 144 cells received 1 MF input.
2.5. Simulating LTP/LTD
The granule cell model was modeled based on data from Wistar rats [9, 19]. By modifying intrinsic excitability and release probability [34, 35], we simulated plasticity in the GrC. Intrinsic excitability was modified by changing ionic current channel gating dynamics. On-off gating characteristics of sodium channel were altered to modify sodium activation and inactivation parameters [13, 26] for higher and lower intrinsic excitability. Three cases were studied where the intrinsic excitability of the GrC is low (low intrinsic excitability), normal (control), and high (high intrinsic excitability). The release probability () of MFs varied from 0.1 to 0.4 for cells with low intrinsic excitability, and from 0.5 to 0.8 for cells with high intrinsic excitability while the control value remained as 0.416 in simulations for normal cells .
2.6. LFP Reconstruction
The extracellular potential of a single granule neuron (see Figure 1(b)) was estimated using NEURON  extracellular mechanism. The mechanism adds two RC compartments (see Figure 1(c)). To understand population code, we reconstructed network evoked LFP response using Laplace equation (see (1)): where is extracellular potential, at boundary condition is the transmembrane current density and ρ is the extracellular resistivity. Each cell generated an extracellular response corresponding to the activation pattern elicited by the mossy fibers. With the granular layer network, an electrode was assumed to be placed at the center. Temporal and spatial delays due to distance from electrode were assumed to be 0–3 ms . The electrode could measure cells that generated extracellular currents that came with a delay of 0–3 ms (see ). Methodology for modeling the latencies used has been detailed elsewhere . Considering the extracellular activity from each granule cell in the region of interest (number of cells = 700, assuming measurements from a tungsten electrode ), we reconstructed evoked LFP response using (10) and (11). Equation (10) adds the delay by padding zeros to linearly time-shift the signal. Equation (11) denotes the process of summing all shifted extracellular signals for all cells linearly. Total signal obtained is the desired evoked LFP: where is the extracellular potential of th cell in the neuronal population within the region of interest. represents the extracellular potential shifted by time delay (0–3 ms) (see ). Equations (10) and (11) were calculated separately. The detailed methodology for reconstructing evoked local field potential (LFP) has been described elsewhere .
We were able to construct two models of granular layer network microcircuit: one using computationally efficient but physiologically limited spiking neurons and other using biophysically detailed multicompartmental neurons and reproduced activation patterns, burst-burst transmission, role of inhibition, and combinatorial coding.
3.1. Time-Windowing Depends on the Feed-Forward Inhibition-Implications from Simple Spiking Network Model
The objective of using a simple spiking model was to understand input-output relationships in terms of firing dynamics in the cerebellar granular layer.
The simulated single granule neuron responses were modeled based on granule neuron electroresponsiveness [9, 28]. Both in vitro (see Figure 2(A)) like behavior with single spike through mossy fibers and in vivo (See Figure 2(C)) like response with burst inputs through mossy fibers were simulated. The responses matched experimental data [9, 19, 28]. The role of feed-forward inhibition [31, 38–40] was also modeled. With inhibition, the granule neuron model showed suppression (see Figure 2(B1)) of spike doublet (see Figure 2(A1)). Synaptic inhibition, because of its delayed activation, controlled generation of the second spike in the doublet [18, 41]. In the in vivo case, the number of spikes (see Figure 2(C)) was reduced due to inhibition (see Figure 2(D)).
Using the spiking granule neuron model, the 1680 granule cell network was reconstructed. The synaptic input in the mossy fibers were reproduced using either a single pulse to mimic electrical stimulation for in vitro simulations or short high-frequency trains mimicking punctuate sensory stimulation for in vivo simulations The network model showed 720 spikes (without inhibition) and reproduced the synaptic activation of the granular layer.
The spiking neuron network model was activated with a center-surround activation pattern , and the raster of spikes were observed in individual cells. Among 1680 cells, 144 cells with 4 MF active, 432 cells with 3 MF, 144 with 2 MF, and 432 with 1 MF active. The configuration was based on Voltage-Sensitive Dye (VSD) imaging  and results matched our previous findings . With LTP and LTD, the numbers of spikes in the network change significantly (see Table 3). The spiking neuron simulations supported a burst-burst transmission modality (see Figure 3) in which high-frequency spike trains are more reliably transmitted.
The time-windowing  depended on the feed-forward inhibitory loop, regulated by the Golgi synapses impinging on the granule neurons. In this model, inhibition was modeled with a delay of 4 ms to account the MF-GoC-GrC circuit. As expected, the feed-forward inhibition reduced the number of spikes from 720 to 576 (see Table 3). The spike raster in the simulations (see Figure 3) showed selective inhibition of granule firing due to blocking of the second input as reported in [1, 2, 18]. The increase in number of spikes in vivo in the network supports the frequency-modulated transition from LTP to LTD . The model also was computationally efficient in comparison to the biophysical model (see Section 3.2) and took 3 s for a 100 ms simulation.
3.2. Spike-Burst Generation and Bidirectional Plasticity
Although simple spiking models allow reconstruction of frequency and amplitude information in terms of firing of constituent cells, role of plasticity and selective pharmacological effects in population code could not be studied. As reported in Section 2.2, we used detailed multicompartmental models [15, 16, 19] to generate a 1680 granule cell network. A detailed model allows to focus on understanding how specific temporal dynamics and the geometry of connections will eventually determine the circuit output, as indicated by the evident anomalies in network functioning and behavior caused by single-gene mutations altering the physiology of single molecules or neurons . The simulations also attempt to understand and indicate certain functional theories of feed-forward inhibition, sparse recoding via, spikes and long-term plasticity.
LTP in granule cells  comprises of variation in release probability and intrinsic excitability. The network model was modified with higher intrinsic excitability observed by changes to sodium channel properties and release probabilities of MF synapses, thereby simulating granule cell LTP. LTD  was also simulated by combining lower intrinsic excitability and low release probability. With varying amounts of excitation, cells generated different number of spikes. In the simulations, granule cells generated repetitive nonadapting spike discharge in response to a continuous stimulus [27, 28, 45]. Simulations on single granule neurons in vitro suggested that granule neurons allowed spike burst generation and resonance in a low-frequency band (between 4 and 10 Hz) [9, 46]. High-frequency bursting [44, 47–49] was simulated to characterize network properties of granule cells in vivo [26, 33].
3.3. Network Excitability Changes with Varying Excitatory Release Probabilities
During events such as in epileptic seizures, heterogeneous spiking activity is noticed . To understand the nature of spiking, we simulated the role of excitation via mossy fiber. We simulated single spikes (low frequency) and bursts (high frequency) so as to understand spiking behavior in vitro and in vivo.
In the case of simulating in vitro behavior (see Section 2.3) in granular layer network, with release probability 0.416 (control) the cells with normal intrinsic excitability receiving 4 excitatory inputs produced spike doublet and cells receiving 3 excitatory inputs produced single spikes (see granule neuron electroresponsiveness in ). Cells receiving 2 excitatory and 1 excitatory inputs did not produce any spikes .
Decreasing MF synapse release probability from 0.3 to 0.1, many granule cells in the network did not generate spikes. With increased release probability of MF synapses from 0.42 to 0.6, there was an increase in number of spikes (the number of spiking cells increased from 192 to 432) with no significant change in first spike latency and spike amplitude .
With the higher release probabilities like 0.7, 0.8 of MF synapses, the number of spikes saturated and the number of spiking cells remained the same (as seen in 0.6 release probability of MF synapses). Varying synaptic release probabilities, it was possible to generate selective responses.
Increasing intrinsic excitability from normal to higher excitability by modifying sodium gating properties (see Section 2.5) showed a significant increase in the spike amplitude (~6%) for all spiking cells, and an increased number of spikes was observed only for the cells with higher release probability and number of active MF synapses. This change corresponded to long-term potentiation in granule neurons  confirming the mechanisms role in spiking and bursting. The number of spiking cells varied from 192 to 432 in the network of 1680 granule cells.
With in vivo inputs, the number of cells showed a greater sensitivity with LTP (see Table 3) and the number of nonspiking cells decreased. LTD showed decrease in firing-nonfiring  cell ratio. The ratio of firing cells did not change with the length of the burst (see Table 4).
3.4. Inhibition and Spike-Count Modulation
Golgi cells can control both the temporal dynamics and the spatial distribution of information transmitted through the cerebellar granular layer network . The strength of the inhibition depends on the number of inhibitory connections and synaptic release probability. The dynamics of the granule cells-Golgi cell circuit were explained by the simultaneous activation of both neurons through the mossy fibers, followed by activation of the feed-forward and feed-back inhibitory loops [18, 51]. The granule and Golgi cell received excitatory inputs from mossy fiber (MF) at the same time. There are two basic patterns of mossy fiber activity that can activate the Golgi cells, namely, protracted frequency-modulated discharges and short high-frequency bursts [48, 52]. The inhibitory input from Golgi cell reaches the granule cell with a loop delay of approximately 4 ms  compared to the mossy fiber input through GABAergic synapses . The inhibition-based time-windowing in granule cells allow one or more spikes and is seemingly regulated by varying inhibitory inputs.
Golgi cells converging through lateral connections onto some granule cell subsets could generate combined inhibition [2, 8]. The impact of the inhibition on granular layer circuitry differs with respect to two different properties: amount of inhibitory connections and the GABAergic release probability. The variation in the number of spikes with and without inhibition was significant in both cases in vitro (see Figure 4(a)) and in vivo (see Figure 4(b)). As expected, LTP showed increased number of spikes compared to control, while LTD showed reduced number of spikes. In vitro LTD suppressed spikes (see Figure 4(a)).
During simulations as seen in vitro (see Figure 4(a)), increased inhibition regulated the spike count rather than affecting the number of spiking cells. Short burst through MF produced ~7 spikes in single neurons, but inhibition showed a sharp modulation by regulating the time-window. A long burst produced slower modulation of spikes in the network.
The increase in inhibitory connections (see Table 5 and Figure 5) to granule cells in the underlying network model decreased number of spikes (see spike count in Figure 5, control refers to release probability being set at normal condition, ), spike amplitude (if the spike rises after the 4 ms time-window when inhibitory inputs reaches the granule cell) and decreased spike latency.
Changing inhibitory (GABAergic) synapse release probability (), spike amplitude, and first spike latency were affected . Spike amplitude decreases whereas spike latency remains unchanged, when varied .
The increase of inhibitory input increases the number of silent cells (visible by the blue plateau in Figure 4), therefore reducing the number of active cells. The simulations indicate that the response of those granule cells that are intensely activated will favor with the generation of a burst, regulated mainly by feed-forward Golgi cell inhibition (see Figure 5 and ).
3.5. Center-Surround Excitation in Populations of Granule Cells
To understand combinatorial effects in the granular network layer and impacts of double mossy fiber bundle stimulation, combined excitation-inhibition in the network was simulated. The “spots” are maps of excitatory activity as seen in the cerebellar granular layer  when MF rosettes were stimulated . In the model configuration (see Table 6), the center of the spot receives stronger excitatory inputs and the consecutive peripheral neurons receive weaker excitatory input, thereby expressing a center-surround configuration (see Figure 6). Both network models could reproduce the firing dynamics  as well as the center-surround structure (see movie, Supplementary Material available online at doi: 10.115/2012/359529). Spiking activity was reconstructed with the morphology (see Section 2.4). A single spike through the mossy fiber activates the center followed by the periphery and the Golgi-granule circuit.
Simulation of LTP and LTD induction in vitro and in vivo on the center-surround spots was modeled by varying release probabilities and intrinsic excitability. The cells in the granular layer network receive GABAergic synaptic inputs equal to the number of excitatory inputs given to the cells in the granular layer network (see Table 5). The high reproducibility indicates that the center-surround organization was a consequence of alternating transitions between burst and silent states at granule cells was not due to the temporal jitter of MFs .
The center-surround structures have complex transmission properties: compared to the surround, the center detects burst on a broader band and emits bursts with shorter lag, higher frequency, and longer duration . Purkinje cells overlaying above these structures may be activated, at the same time, enhance inhibition around them, explaining the spot-like organization of molecular layer responses in vivo .
3.6. Local Field Potential and Selective Blocking of NMDA
Understanding population code through reconstructions was essentially done to suggest how encoding of spike information may happen in cerebellar cortex. Currently, encoding of population activity is explored in microcircuits via comparisons of spatiotemporal properties of simulated neural activity and with experimental measurements using multielectrode recordings [55, 56] or two-photon imaging of activity in blocks of tissue [57, 58]. Evoked responses from granular layer in vitro  have been reconstructed computationally [13, 26]. We used the “spot” to generate and test nature of local field potentials. The postsynaptic evoked LFP response varied as per input pattern and for a combination of 3 MF and 4 MF synaptic activation, spikes were generated. With 4 MF synapses active, a doublet was seen . Correspondingly the responses generated wave and the doublet caused the wave. Inhibition at time = 24 ms via GABAergic synapses suppressed the spike doublet and thereby suppressed the wave  in the evoked LFP response.
Different segments of the network generated varied evoked LFP signals due to the nature of excitatory-inhibitory balance in the network reflecting a relationship different from the extracellular components of a single neuron (see Figure 9 in ). We assumed extracellular space in granular layer to be isopotential  due to close packing of granule neurons. The simulations closely followed experimental results , suggesting that electrotonic compactness of granule neurons contribute to the seemingly linear relationship from granule cell clusters in the granular layer extracellular space. The variations in the nature of spike with number of spiking cells could suggest that sparse coding could be preserved as suggested by Marr  and Albus .
Blocking NMDA receptors  in granule neurons showed reduced excitation. Selective disabling of NMDA receptors, as noticed in mice with NR2A/NR2B  mutations, showed decreased number of spikes which is also seen as a change in amplitude compared to control (Figure 7(a)). In order to predict on the nature of such mutations affecting network computation (in addition to affecting the number of spikes), we randomly disabled (1%, 5%, and 10% of total cells) NMDA receptors in the network (see circles in Figures 7(b), 7(c), and 7(d)) and reconstructed the local field response.
Disabling NMDA receptors in 1–10% of cells showed a 2.5% decrease in number of spikes in a spot (the number of spiking neurons in a 720-cell “spot” changed from 240 to 234). The network model clearly showed a “seemingly linear” outlook in propagating the nonlinearities of individual neurons in population code (evoked LFP, see Figures 7(a) to 7(d)). This “sense of linearity” in population code was observed also when the number of affected neurons was very low (neurons with NMDA disabled were only 1–10% of total cells). NMDA knockout mice show errors in cerebellar motor learning . Plasticity changes were also reflected in the evoked LFP waves. As seen in reconstructions based on single granule neuron simulations , the LFP simulations on the network showed that LTP and LTD were accompanied by changes in the proportion of discharging granule cells (data not shown).
Exploring the geometry of excitation and inhibition in cerebellar granular layer, the simulations highlight the modulatory role of inhibitory inputs on the activities of granule cells. The paper details the effects of combined excitation and combined feed-forward inhibition  on spiking in the granular layer. The study did not simulate feed-back inhibition coming from Golgi cell since it did not affect the modulation of 4 ms time-window that happens because of the early mossy inputs.
The simulations varying inhibition suggest that granule neurons can generate selective responses by varying synaptic strengths. The increase in spikes and modulation during plasticity indicates that the circuit is well adapted to generate enhanced responses such as in theta-burst patterns .
Both in vitro and in vivo simulations indicate that inhibitory input cannot completely block excitation in the network. However, it acts as a modulator that regulates the postsynaptic excitability. Both models support that burst-burst transmission modality in granule neuron and the granular layer through which high-frequency spike trains are more reliably transmitted. The consequences of transformation of spike inputs from mossy fibers to corresponding codes suggest the variable impulse response scheme indicated by previous study  and suggest the granular layer network also operates as an adaptive filter.
The variations of excitatory inputs (without combination of inhibition) showed differences in number of spikes and spike amplitude and did not show variations in first-spike latency . The most promising outcome in variation of spikes and network spiking behavior was with the induction of LTP/LTD where both intrinsic excitability and excitatory release probabilities [21, 43] change the nature of information flow.
Simple spiking neuron models can be tuned to function as network models for accessing timing information. Spatial information in network models  was not seen while using spiking models. Synaptic functions in spiking models are not very reliable. Artificial models have limitations unlike biophysical models for understanding certain population activities like generation of LFP .
The detailed model was used as a test bench to explore the parameter space and induced plasticity. Epileptic seizure-like symptoms seen in voltage-gated sodium channel binding-related knockout mice granule neurons [7, 50] suggests that sparse and asynchronous neuronal activity can evolve into a single hypersynchronous cluster with elevated spiking rates at seizure initiation. The detailed network model suggests that LTP favors burst-burst transmission favoring high-frequency spikes. The presynaptic mechanism coexisted with postsynaptic regulation of ionic channels, which played a major role in determining the granule cell output firing frequency. Intrinsic bursting and modulatory effects of inhibition can be seen by mechanistic control of number of spikes in a granule cell.
With increased excitation, along with an increase in spikes, first-spike latency also decreased. This will also impact the local field potential and could probably explain the observations in vitro . Both in vitro and in vivo simulations indicate that the number of spikes was dependent on the release probability of the synapses, while higher or lower intrinsic excitability caused slight change in spike amplitude.
The key role of local circuit inhibition in determining granular layer combinatorial operations was supported by several model-based predictions. Increasing active inhibitory connections saw lesser number of spikes in the network. In vivo bursts along mossy fibers combined with inhibitory input showed a consistent reduction of at least one spike as inhibition increased. The simulations indicated that the response of those granule cells that are intensely activated will favor with the generation of a burst, whose duration is limited by a brisk feed-forward inhibition in the Golgi cell. Inhibition controlled the number of spikes, thereby modulating spike transmission in the granular layer. The simulations suggest that erratic spikes in the mossy fibers will not be efficiently transmitted so that the burst-burst mechanism would indeed play a role in secure transmission along the mossy fiber pathway . The studies also show that excitation and inhibition may consequently allow complex patterns to be processed .
The paper also shows population signals and effects of mechanism changes on individual neuron affecting population code generated by the network. Reconstructing extracellular properties indicated that plasticity may have similar mechanisms of burst regulation as granule cell burst initiation and may implement an adaptable delay affecting downstream activation into circuitry. The granular layer model indicates a “seemingly linear” tendency to propagate the nonlinearities of individual neurons via the population code even when the variations are little (affected cells 1–10% of total). The simulations suggest that a combined mechanism of NMDA blocking the After-hyperpolarization (AHP) and role of inhibition can help reconstruct transient suppression of spikes in vitro reported during seizures.
The studies on intensity of mossy fiber synapses and inhibitory synapses help to understand spatiotemporal operations  in the cerebellar granular layer. Combining granule neurons and Golgi cell, this study will help to reveal coincidence detection properties and spatial pattern separation . This work is a preliminary start in modeling to understanding long sought spatiotemporal filtering predicted by the motor learning theory .
Simulations suggest how cerebellum granular layer processes spike information and how afferent information may reach cerebellar cortex and predict how spikes are processed as indicated in the sparse recoding hypothesis . The role of inhibition and plasticity may help fine tune the “sparseness” of the code as indicated in Marr’s theory [3, 4]. To evaluate the exact role of firing, a closer view of cells in the region of interest may be needed. The experimental testing of these predictions will require further electrophysiological and imaging investigations of granular layer activity and computational modelling of the cerebellum  and of the cerebro-cerebellar control loops [62, 63].
This work derives direction and ideas from the chancellor of Amrita University, Sri Mata Amritanandamayi Devi. The authors would like to thank Priyanka James, Nimshitha Abdulmanaph, Harilal Parasuram, and Sergio Solinas for their work and support in making this paper. This work is supported by Grants SR/CSI/49/2010 and SR/CSI/60/2011 from the Department of Science and Technology, Government of India.
- S. Solinas, T. Nieus, and E. D. Angelo, “A realistic large-scale model of the cerebellum granular layer predicts circuit spatio-temporal filtering properties,” Frontiers in Cellular Neuroscience, vol. 4, article 12, 2010.
- E. D'Angelo and C. I. De Zeeuw, “Timing and plasticity in the cerebellum: focus on the granular layer,” Trends in Neurosciences, vol. 32, no. 1, pp. 30–40, 2009.
- J. S. Albus, “A theory of cerebellar function,” Mathematical Biosciences, vol. 10, no. 1-2, pp. 25–61, 1971.
- D. Marr, “A theory of cerebellar cortex,” Journal of Physiology, vol. 202, no. 2, pp. 437–470, 1969.
- F. Prestori, P. Rossi, B. Bearzatto et al., “Altered neuron excitability and synaptic plasticity in the cerebellar granular layer of juvenile prion protein knock-out mice with impaired motor control,” Journal of Neuroscience, vol. 28, no. 28, pp. 7091–7103, 2008.
- S. P. Perrett, B. P. Ruiz, and M. D. Mauk, “Cerebellar cortex lesions disrupt learning-dependent timing of conditioned eyelid responses,” Journal of Neuroscience, vol. 13, no. 4, pp. 1708–1718, 1993.
- M. Goldfarb, J. Schoorlemmer, A. Williams et al., “Fibroblast growth factor homologous factors control neuronal excitability through modulation of voltage-gated sodium channels,” Neuron, vol. 55, no. 3, pp. 449–463, 2007.
- J. Mapelli, D. Gandolfi, and E. D'Angelo, “Combinatorial responses controlled by synaptic inhibition in the cerebellum granular layer,” Journal of Neurophysiology, vol. 103, no. 1, pp. 250–261, 2010.
- E. D'Angelo, T. Nieus, A. Maffei et al., “Theta-frequency bursting and resonance in cerebellar granule cells: experimental evidence and modeling of a slow K+-dependent mechanism,” Journal of Neuroscience, vol. 21, no. 3, pp. 759–770, 2001.
- P. Dean, J. Porrill, C. F. Ekerot, and H. Jörntell, “The cerebellar microcircuit as an adaptive filter: experimental and computational evidence,” Nature Reviews Neuroscience, vol. 11, no. 1, pp. 30–43, 2010.
- J. F. Medina and M. D. Mauk, “Computer simulation of cerebellar information processing,” Nature Neuroscience, vol. 3, pp. 1205–1211, 2000.
- R. Maex and E. De Schutter, “Synchronization of Golgi and granule cell firing in a detailed network model of the cerebellar granule cell layer,” Journal of Neurophysiology, vol. 80, no. 5, pp. 2521–2537, 1998.
- H. Parasuram, B. Nair, G. Naldi, E. D'Angelo, and S. Diwakar, “A modeling based study on the origin and nature of evoked post-synaptic local field potentials in granular layer,” Journal of Physiology Paris, vol. 105, no. 1–3, pp. 71–82, 2011.
- R. Naud, N. Marcille, C. Clopath, and W. Gerstner, “Firing patterns in the adaptive exponential integrate-and-fire model,” Biological Cybernetics, vol. 99, no. 4-5, pp. 335–347, 2008.
- S. Solinas, L. Forti, E. Cesana, J. Mapelli, E. De Schutter, and E. D' Angelo, “Fast-reset of pacemaking and theta-frequency resonance patterns in cerebellar Golgi cells: simulations of their impact in vivo,” Frontiers in Cellular Neuroscience, vol. 1, article 4, 2007.
- S. Solinas, L. Forti, E. Cesana, J. Mapelli, E. De Schutter, and E. D' Angelo, “Computational reconstruction of pacemaking and intrinsic electroresponsiveness in cerebellar Golgi cells,” Frontiers in Cellular Neuroscience, vol. 1, article 2, 2007.
- D. A. McCormick, Z. Wang, and J. Huguenard, “Neurotransmitter control of neocortical neuronal activity and excitability,” Cerebral Cortex, vol. 3, no. 5, pp. 387–398, 1993.
- J. Mapelli and E. D'Angelo, “The spatial organization of long-term synaptic plasticity at the input stage of cerebellum,” Journal of Neuroscience, vol. 27, no. 6, pp. 1285–1296, 2007.
- S. Diwakar, J. Magistretti, M. Goldfarb, G. Naldi, and E. D'Angelo, “Axonal Na+ channels ensure fast spike activation and back-propagation in cerebellar granule cells,” Journal of Neurophysiology, vol. 101, no. 2, pp. 519–532, 2009.
- N. Abdulmanaph, H. Parasuram, B. Nair, and S. Diwakar, “Modeling granular layer local field potential using single neuron and network based approaches to predict LTP/LTD in extracellular recordings,” in Proceedings of Neurocomp, Lyon, France, 2010.
- T. Nieus, E. Sola, J. Mapelli, E. Saftenku, P. Rossi, and E. D'Angelo, “LTP regulates burst initiation and frequency at mossy fiber-granule cell synapses of rat cerebellum: experimental observations and theoretical predictions,” Journal of Neurophysiology, vol. 95, no. 2, pp. 686–699, 2006.
- L. Chen and L. Y. M. Huang, “Protein kinase C reduces Mg2+ block of NMDA-receptor channels as a mechanism of modulation,” Nature, vol. 356, no. 6369, pp. 521–523, 1992.
- C. J. McBain and R. Dingledine, “Dual-component miniature excitatory synaptic currents in rat hippocampal CA3 pyramidal neurons,” Journal of Neurophysiology, vol. 68, no. 1, pp. 16–27, 1992.
- N. Ropert, R. Miles, and H. Korn, “Characteristics of miniature inhibitory postsynaptic currents in CA1 pyramidal neurones of rat hippocampus,” Journal of Physiology, vol. 428, pp. 707–722, 1990.
- T. S. Otis, Y. De Koninck, and I. Mody, “Characterization of synaptically elicited responses using patch-clamp recordings in rat hippocampal slices,” Journal of Physiology, vol. 463, pp. 391–407, 1993.
- S. Diwakar, P. Lombardo, S. Solinas, G. Naldi, and E. D'Angelo, “Local field potential modeling predicts dense activation in cerebellar granule cells clusters under LTP and LTD control,” PLoS ONE, vol. 6, no. 7, Article ID e21928, 2011.
- E. D'Angelo, G. De Filippi, P. Rossi, and V. Taglietti, “Synaptic excitation of individual rat cerebellar granule cells in situ: evidence for the role of NMDA receptors,” Journal of Physiology, vol. 484, no. 2, pp. 397–413, 1995.
- E. D'Angelo, G. De Filippi, P. Rossi, and V. Taglietti, “Ionic mechanism of electroresponsiveness in cerebellar granule cells implicates the action of a persistent sodium current,” Journal of Neurophysiology, vol. 80, no. 2, pp. 493–503, 1998.
- W. Rall, “Branching dendritic trees and motoneuron membrane resistivity,” Experimental Neurology, vol. 1, no. 5, pp. 491–527, 1959.
- P. Rossi, E. Sola, V. Taglietti et al., “NMDA receptor 2 (NR2) C-terminal control of NR open probability regulates synaptic transmission and plasticity at a cerebellar synapse,” Journal of Neuroscience, vol. 22, no. 22, pp. 9687–9697, 2002.
- L. Mapelli, P. Rossi, T. Nieus, and E. D'Angelo, “Tonic activation of receptors reduces release probability at inhibitory connections in the cerebellar glomerulus,” Journal of Neurophysiology, vol. 101, no. 6, pp. 3089–3099, 2009.
- M. L. Hines and N. T. Carnevale, “The NEURON simulation environment,” Neural Computation, vol. 9, no. 6, pp. 1179–1209, 1997.
- L. Roggeri, B. Rivieccio, P. Rossi, and E. D'Angelo, “Tactile stimulation evokes long-term synaptic plasticity in the granular layer of cerebellum,” Journal of Neuroscience, vol. 28, no. 25, pp. 6354–6359, 2008.
- E. Sola, F. Prestori, P. Rossi, V. Taglietti, and E. D'Angelo, “Increased neurotransmitter release during long-term potentiation at mossy fibre-granule cell synapses in rat cerebellum,” Journal of Physiology, vol. 557, no. 3, pp. 843–861, 2004.
- S. Armano, P. Rossi, V. Taglietti, and E. D'Angelo, “Long-term potentiation of intrinsic excitability at the mossy fibergranule cell synapse of rat cerebellum,” Journal of Neuroscience, vol. 20, no. 14, pp. 5208–5216, 2000.
- P. James, N. Abdulmanaph, B. Nair, and S. Diwakar, “Exploring input-output characteristics of the cerebellar granule neuron: role of synaptic inhibition, spike timing and plasticity,” in Proceeding of Neurocomp, Lyon, France, 2010.
- R. A. Silver, S. G. Cull-Candy, and T. Takahashi, “Non-NMDA glutamate receptor occupancy and open probability at a rat cerebellar synapse with single and multiple release sites,” Journal of Physiology, vol. 494, no. 1, pp. 231–250, 1996.
- P. Rossi, L. Mapelli, L. Roggeri et al., “Inhibition of constitutive inward rectifier currents in cerebellar granule cells by pharmacological and synaptic activation of receptors,” European Journal of Neuroscience, vol. 24, no. 2, pp. 419–432, 2006.
- J. C. Eccles, M. Ito, and J. Szentagothai, The Cerebellum as a Neuronal Machine, Springer, Berlin, Germany, 1967.
- C. Palay and S. L. Palay, Cerebellar Cortex: Cytology and Organization, Springer, Berlin, Germany, 1974.
- E. D'Angelo, “The critical role of Golgi cells in regulating spatio-temporal integration and plasticity at the cerebellum input stage,” Frontiers in Neuroscience, vol. 2, pp. 35–46, 2008.
- S. Subramaniyam, C. Medini, B. Nair, and S. Diwakar, “Modeling spatio-temporal processing in cerebellar granular layer and effects of controlled inhbition on plasticity,” in Proceedings of Neurocomp, Lyon, France, 2010.
- A. D'Errico, F. Prestori, and E. D'Angelo, “Differential induction of bidirectional long-term changes in neurotransmitter release by frequency-coded patterns at the cerebellar input,” Journal of Physiology, vol. 587, no. 24, pp. 5843–5857, 2009.
- H. Jörntell and C. F. Ekerot, “Properties of somatosensory synaptic integration in cerebellar granule cells in vivo,” Journal of Neuroscience, vol. 26, no. 45, pp. 11786–11797, 2006.
- S. G. Brickley, S. G. Cull-Candy, and M. Farrant, “Development of a tonic form of synaptic inhibition in rat cerebellar granule cells resulting from persistent activation of receptors,” Journal of Physiology, vol. 497, no. 3, pp. 753–759, 1996.
- J. Magistretti, L. Castelli, L. Forti, and E. D'Angelo, “Kinetic and functional analysis of transient, persistent and resurgent sodium currents in rat cerebellar granule cells in situ: an electrophysiological and modelling study,” Journal of Physiology, vol. 573, no. 1, pp. 83–106, 2006.
- P. Chadderton, T. W. Margie, and M. Häusser, “Integration of quanta in cerebellar granule cells during sensory processing,” Nature, vol. 428, no. 6985, pp. 856–860, 2004.
- E. A. Rancz, T. Ishikawa, I. Duguid, P. Chadderton, S. Mahon, and M. Häusser, “High-fidelity transmission of sensory information by single cerebellar mossy fibre boutons,” Nature, vol. 450, no. 7173, pp. 1245–1248, 2007.
- N. H. Barmack and V. Yakhnitsa, “Functions of interneurons in mouse cerebellum,” Journal of Neuroscience, vol. 28, no. 5, pp. 1140–1152, 2008.
- W. Truccolo, J. A. Donoghue, L. R. Hochberg et al., “Single-neuron dynamics in human focal epilepsy,” Nature Neuroscience, vol. 14, no. 5, pp. 635–641, 2011.
- A. Maffei, F. Prestori, P. Rossi, V. Taglietti, and E. D'Angelo, “Presynaptic current changes at the mossy fiber-granule cell synapse of cerebellum during LTP,” Journal of Neurophysiology, vol. 88, no. 2, pp. 627–638, 2002.
- M. Kase, D. C. Miller, and H. Noda, “Discharges of Purkinje cells and mossy fibres in the cerebellar vermis of the monkey during saccadic eye movements and fixation,” Journal of Physiology, vol. 300, pp. 539–555, 1980.
- A. Volny-Luraghi, R. Maex, B. Vos, and E. De Schutter, “Peripheral stimuli excite coronal beams of Golgi cells in rat cerebellar cortex,” Neuroscience, vol. 113, no. 2, pp. 363–373, 2002.
- T. Yamazaki and S. Tanaka, “A spiking network model for passage-of-time representation in the cerebellum,” European Journal of Neuroscience, vol. 26, no. 8, pp. 2279–2292, 2007.
- G. Buzsáki, “Large-scale recording of neuronal ensembles,” Nature Neuroscience, vol. 7, no. 5, pp. 446–451, 2004.
- M. A. L. Nicolelis and S. Ribeiro, “Multielectrode recordings: the next steps,” Current Opinion in Neurobiology, vol. 12, no. 5, pp. 602–606, 2002.
- K. Ohki, S. Chung, Y. H. Ch'ng, P. Kara, and R. C. Reid, “Functional imaging with cellular resolution reveals precise microarchitecture in visual cortex,” Nature, vol. 433, no. 7026, pp. 597–603, 2005.
- C. Stosiek, O. Garaschuk, K. Holthoff, and A. Konnerth, “In vivo two-photon calcium imaging of neuronal networks,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 12, pp. 7319–7324, 2003.
- C. E. Andreescu, F. Prestori, F. Brandalise et al., “NR2A subunit of the N-methyl d-aspartate receptors are required for potentiation at the mossy fiber to granule cell synapse and vestibulo-cerebellar motor learning,” Neuroscience, vol. 176, pp. 274–283, 2011.
- R. L. Jenison, “Decoding first-spike latency: a likelihood approach,” Neurocomputing, vol. 38–40, pp. 239–248, 2001.
- M. Fujita, “Adaptive filter model of the cerebellum,” Biological Cybernetics, vol. 45, no. 3, pp. 195–206, 1982.
- M. Uusisaari and E. de Schutter, “The mysterious microcircuitry of the cerebellar nuclei,” Journal of Physiology, vol. 589, no. 14, pp. 3441–3457, 2011.
- L. Bao, A. M. Rapin, E. C. Holmstrand, and D. H. Cox, “Elimination of the channel's high-affinity Ca2+ sensitivity,” Journal of General Physiology, vol. 120, no. 2, pp. 173–189, 2002.