REVIEW article

Front. Comput. Neurosci., 28 October 2022

Volume 16 - 2022 | https://doi.org/10.3389/fncom.2022.1006989

Computational models of neurotransmission at cerebellar synapses unveil the impact on network computation

  • 1. Department of Brain and Behavioral Sciences, University of Pavia, Pavia, Italy

  • 2. IRCCS Mondino Foundation, Brain Connectivity Center, Pavia, Italy

Article metrics

View details

10

Citations

5,6k

Views

1,5k

Downloads

Abstract

The neuroscientific field benefits from the conjoint evolution of experimental and computational techniques, allowing for the reconstruction and simulation of complex models of neurons and synapses. Chemical synapses are characterized by presynaptic vesicle cycling, neurotransmitter diffusion, and postsynaptic receptor activation, which eventually lead to postsynaptic currents and subsequent membrane potential changes. These mechanisms have been accurately modeled for different synapses and receptor types (AMPA, NMDA, and GABA) of the cerebellar cortical network, allowing simulation of their impact on computation. Of special relevance is short-term synaptic plasticity, which generates spatiotemporal filtering in local microcircuits and controls burst transmission and information flow through the network. Here, we present how data-driven computational models recapitulate the properties of neurotransmission at cerebellar synapses. The simulation of microcircuit models is starting to reveal how diverse synaptic mechanisms shape the spatiotemporal profiles of circuit activity and computation.

Introduction

Synapses are specialized structures that regulate inter-neuronal communication through a process called neurotransmission. While the foundations of synaptic biophysics have been developed in the last century (Koch, 1984; Maxwell Cowan et al., 2001), the advent of sophisticated single-cell recording techniques, such as patch-clamping and calcium imaging, has recently allowed the characterization of specific properties of several synapses in the mammalian brain, in most cases using rodents as the experimental model (Hickman et al., 2017). These results have revealed a previously undisclosed complexity and variety of mechanisms, requiring a further understanding of their functional implications through advanced computational models. Cellular and synaptic biophysics can be translated into mathematical equations (Koch, 1984; Segev, 1998; Brodland, 2015) to generate models of neurons and networks, thus adding a new dimension to the computational investigation of neurotransmission (D’Angelo and Jirsa, 2022). A relevant case is the cerebellar circuit, the synapses of which exhibit a variety of properties. These include presynaptic release and postsynaptic receptor mechanisms, which generate a large spectrum of postsynaptic currents, membrane depolarization kinetics, and short-term plasticity patterns (D’Angelo et al., 2016; De Schepper et al., 2021). In this review, we compared the neurotransmission mechanisms at different cerebellar synapses and the corresponding computational models that unveil the impact of synaptic diversity on microcircuit computation.

Experimental data and reconstruction of synaptic models

Experimental techniques used to investigate neurotransmission have evolved considerably over the last century (Sotelo, 2008; VoŽeh, 2015; Cavero et al., 2017; Bentivoglio et al., 2019). The extensive application of 3D morphological reconstructions (Guo et al., 2022), immunochemistry, patch-clamp recordings, calcium imaging, genomics, and proteomics has allowed the creation of atlases, such as the Allen Brain Atlas (Gilbert, 2018), which summarizes the knowledge about brain neurons and circuits along with the distribution of synaptic receptors and molecules, both in rodents and humans. Computational models can integrate this large set of experimental observations to reconstruct a variety of synaptic receptor kinetics (AMPA, NMDA, and GABA) and neurotransmission properties (D’Angelo and Jirsa, 2022) (Box 1).

BOX 1. Simulation environments.

Simulation environments are a crucial part of biophysical research, providing tools for modeling and testing new concepts and algorithms (Bower, 2015). To choose the most appropriate environment, it is essential to evaluate the levels of biological organization, availability of experimental data, and requirements for computing power.

At the molecular level, highly sophisticated simulators can solve biochemical networks (CellDesigner) (Funahashi et al., 2008), reaction-diffusion systems (STEPS) (Hepburn et al., 2012; Chen and de Schutter, 2017), protein folding and protein-protein interactions (NAMD) (Sotomayor, 2015; Nowak, 2016), and gene regulatory networks (CYTOSCAPE) (Shannon, 2003; Cline et al., 2007). Although these modeling types are grounded in biological reality, they require large amounts of heterogeneous data to visualize large-scale biomolecular networks.

Large-scale simulations can be used to reproduce many different types of tissue and/or organ mechanisms (Dos Santos et al., 2015). Specifically, neural mass models (Jansen and Rit, 1995; Spiegler et al., 2011) are among the most popular mathematical models of brain activity, as they enable the rapid simulation of large populations of neurons and synapses at a spatial scale compatible with electrophysiological experiments (Schirner et al., 2018; Cakan and Obermayer, 2020). An open-source simulation platform optimized for realistic brain connectivity/geometry and neural mass models, The Virtual Brain (TVB) (Sanz Leon et al., 2013), has recently been developed for whole-brain dynamic simulations combining a large-scale brain network model with neuroimaging data, including population activity like electroencephalogram (EEG) Magnetoencephalography (MEG), highly resolved metabolic/vascular signals like the Functional magnetic resonance imaging (fMRI) (Alahmadi et al., 2015), and measures of neuronal connectivity like in the case of diffusion tensor imaging (DTI) (Palesi et al., 2021).

The most popular single-neuron models are:

(1) Conductance-based models, which, with a variety of ion channels of the Hodgkin–Huxley (HH) type (Hodgkin and Huxley, 1952a,b; Hille, 2001), can characterize the biophysical mechanisms of cells in great detail.

(2) Integrate-and-fire (IF) models, which are widely used to describe spiking and bursting behavior (Stein, 1967; Tuckwell, 1988). The most well-known are NEST (Eppler, 2008; Nowke et al., 2015), Brain2 (Stimberg et al., 2019) and LEMS (Cannon et al., 2014). NEURON (Hines and Carnevale, 2001, 2008; Hines et al., 2009), CORENEURON (Kumbhar et al., 2019), GENESIS (Bower and Beeman, 2007), and Arbor (Akar et al., 2019) are well-known simulation environments that allow biophysically realistic single-neuron reconstruction and simulation.

(3) Spike-response neuron (SRN) models, which combine the biological plausibility of HH-type dynamics and the computational performance of IF neurons to simulate tens of thousands of spiking neurons in real time (Izhikevich, 2003). Numerous simulators can combine single-neuron models to generate spiking microcircuit models using specific modeling workflows.

Mechanisms and models of synaptic transmission

Chemical synapses

A typical chemical synapse requires juxtaposition between the axon of one neuron and the dendrite of another (Südhof and Malenka, 2008). The two sites, called pre- and postsynaptic, do not touch each other, leaving a space called the synaptic cleft. The fusion of vesicles with the presynaptic membrane releases the neurotransmitter, which diffuses to the postsynaptic receptors, where it can open ionic channels (ionotropic receptors) (Andreae and Burrone, 2014; Batool et al., 2019) or interact with G proteins and various enzymes to produce second messengers (metabotropic receptors) (Marcaggi et al., 2009). Thus, the models of neurotransmission should contain a representation of (1) the vesicle cycle, (2) neurotransmitter diffusion, and (3) receptor activation (Figure 1). The properties of some cerebellar synaptic computational models are listed in Table 1.

FIGURE 1

FIGURE 1

Biophysical representation of synaptic neurotransmission mechanisms. (A) Schematic representation of the processes taking place in a synapse (Left). The figure shows a synaptic vesicle in three different states: (X) refilled (immediately available store), (Y) releasing (fusion state) and (Z) depleted (under recovery). Following the release, the neurotransmitter diffuses through the presynaptic cleft and binds to postsynaptic receptors. The differential equations describing the temporal evolution of release, diffusion, and gating (for a generic C = O = D scheme) are indicated to the right. (B) The kinetic schemes (Markov chains) describe the transitions of AMPA, NMDA, and GABA receptors through C, D, and O states. These receptors function as ligand-operated channels that can visit multiple states (C: closed, O: open, I: inactive) regulating ionic current flow. Examples of traces are shown at the bottom (the NMDA current has two traces at different holding potentials, as indicated).

TABLE 1

UτfacilτrecτinGmax AMPASubtypeGmax NMDASubtypeGABA-ASubtypeREF
Mf-GrC0.435811,200GluR218,800NR2A-NR2Ca
0.46No presynaptic1,200GluR2b
0.41610.835.13700GluR2c
0.41610.835.13-10,000NR2A-NR2Cc
0.428511,200GluR218,800NR2A-NR2Cd
GoC-GrC0.35583611,050GluR2α1/6e
Mf-GoC0.435811,200GluR210,000NR2Bf
GrC-GoC pf0.45435.111,200GluR2f
GrC-GoC aa0.45435.111,200GluR210,000NR2Bf
GrC-GoC (pf-aa)No presynaptic2,000GluR2g
GrC-PC pf0.135435.112,800GluR2g
GrC-PC aa0.135435.112,800GluR2g
GrC-SC0.1510.835.112,300GluR210,000NR2Bh
GrC-BC0.1510.835.112,300GluR210,000NR2Bi
SC-PC0.3541512,600α1h
BC-PC0.3541512,600α1i
SC-SC0.42438.711,600α1h
BC-BC0.42438.711,600α1i

Properties of cerebellar synapses.

The properties of the cerebellar circuit synapses have been parametrized by fitting synaptic current trains (Nieus et al., 2006, 2014). The specific models have been reported in references a (Nieus et al., 2006), b (Saftenku, 2005), c (Diwakar et al., 2011), d (Sudhakar et al., 2017), e (Nieus et al., 2014), f (Masoli et al., 2020a), g (Masoli and D’Angelo, 2017), h (Rizza et al., 2021), and i (Rizza et al., 2021), and in unpublished data. In addition, all the parameters to be used in local microcircuits (refs a-f), are now used in the BSB model of the whole cerebellar cortical network (De Schepper et al., 2021). It should be noted that the parallel fiber synapses are strongly facilitating (low U), while all other synapses are depressing (high U) or neutral (depending on the balance between U, τfacil, and τrec). The NMDA conductance value is higher than AMPA due to low channel opening probability.

Presynaptic site—Vesicle cycle

The presynaptic site contains complex biochemical machinery for the production, storage, and release of vesicles that contain neurotransmitters (Figure 1A). To reduce the complexity into a manageable set of equations, the Tsodyks–Markram model was built to account for the main states (X, Y, Z). When a presynaptic spike arrives, the variable X represents the neurotransmitter available for release, Y the amount of neurotransmitter released, and Z the recovered neurotransmitter (Tsodyks et al., 1998). The entire process is controlled by the time constants of recovery of the releasable neurotransmitter (τrec), facilitation (τfacil), and inactivation (τin), and the probability of release (U). These parameters effectively regulate short-term synaptic plasticity, including facilitation and depression (D’Angelo et al., 2005). Fitting these models to experimental data allows for the investigation of presynaptic dynamics. The Tsodyks–Markram model has been successfully used in different large-scale networks, such as the neocortex (Markram et al., 2015), hippocampus (Romani et al., 2022), striatum, and cerebellum (Solinas et al., 2010; Sudhakar et al., 2017; Florimbi et al., 2021).

Synaptic cleft—Neurotransmitter diffusion

Neurotransmitter homeostasis in and around a synapse involves complex processes, such as diffusion, binding to receptors, and uptake by glial cells (Figure 1A). Various neurotransmitter diffusion-based modeling studies have focused on synaptic receptor activation (Clements et al., 1992; Franks et al., 2002) and evaluated the access of neurotransmitters to synaptic and extra-synaptic locations, considering sequestration by glial cells (Barbour, 2001; Diamond, 2005). These models have been developed and solved using analytical (Kleinle et al., 1996), continuum (Rusakov and Kullmann, 1998; Rusakov, 2001), and stochastic approaches (Franks et al., 2002; Tao and Nicholson, 2004; Savtchenko and Rusakov, 2007; Zheng et al., 2008), maintaining, as a common simplifying feature, the assumption of simple geometrical configurations of the synaptic cleft.

Postsynaptic site—Receptor activation and current generation

Markov kinetic models can accurately reproduce electrophysiological and biophysical aspects of different ligand-gated receptors. The simplest Markov model is the Markov chain, which is routinely used to model transitions between states (e.g., closed, open, and inactivated states) (Yamada and Zucker, 1992; Destexhe et al., 1994).

AMPA receptors are heterotetramers composed of two pairs of dimers, in which GluR2 is the main subunit. All subunits (1–4) are distributed throughout the cerebellum: GluR1 and GluR2 are highly expressed in all layers, while GluR3 and GluR4 are detected at a lower concentration in the granular layer (Day et al., 1995; Greger et al., 2007), on Purkinje cell (PC) dendrites, which are required to induce long-term potentiation (LTP) and on Bergmann glia (Saab et al., 2012). The typical receptor of GrCs is composed of GluR2 and GluR4 (DiGregorio et al., 2007) and their participation in expansion recoding and associative memory has been reported (Kita et al., 2021). The most basic AMPA receptor model has three states (open, closed, and desensitized), which can be extended to account for subconductance states (e.g., multiple open, desensitization, and deactivation states), with transitions between states described by Markov models (Koike et al., 2000; Wadiche and Jahr, 2001; DiGregorio et al., 2007; Bennani et al., 2009).

NMDA receptors are heterotetramers composed of two obligatory NR1 and NR2 subunits, which exist as four distinct subtypes (A-D). In general, Markov models for NMDA receptors contain more states compared to those for AMPA receptors, to consider the different biophysical properties conferred by NR2 subunit variability (Figure 1B; Laube et al., 1998; Cull-Candy et al., 2001; Cull-Candy and Leszkiewicz, 2004). Kinetic analysis of postsynaptic currents induced by glutamate during prolonged exposure to glycine has revealed two open states and five shut states (three closed and two desensitized states) (Iacobucci and Popescu, 2018). If we consider the modulation of NMDA receptor desensitization by glycine (Benveniste et al., 1990) or calcium-induced inactivation (Vyklický, 1993), 12 states can be predicted.

GABAA receptors are heteropentamers formed by various combinations of subunits: α1-6, β1-4, γ1-4, δ, ε, φ, ρ1-3, and θ (Olsen and Sieghart, 2008). The most common stoichiometry in the mammalian CNS is thought to be two α, two β, and one γ subunits (s) (Baumann et al., 2001, 2002; Olsen and Sieghart, 2008). Although this extensive structural diversity is not completely understood, some subunits confer distinct biophysical and pharmacological properties, especially α subunits, which are thought to play a key role in GABAA receptor function. For example, GABAA receptors containing α1–3 are intra-synaptic and mediate phasic inhibition, whereas receptors containing α4–6 are largely extra-synaptic and mediate tonic inhibition (Farrant and Nusser, 2005). Several functional studies have shown that gating of the GABAA receptor is complex, and there is no consensus kinetic model (Gielen et al., 2012; Szczot et al., 2014; Kisiel et al., 2018). Brief GABA pulses to outside-out macro patches excised from neurons produce macroscopic currents characterized by rapid activation and bi-exponential decay, and desensitization kinetics consistent with long closed states (Maconochie et al., 1994; Jones and Westbrook, 1995). Instead, kinetic analysis of single-channel currents from transiently transfected cells provides simple models for GABAA receptor gating containing multiple open and closed components (Lema and Auerbach, 2006).

Electrical synapses

Electrical synapses are widespread in the mammalian nervous system, including in the human brain. Gap junctions (GJs) mediate electrical synaptic transmission and are composed of two hemichannels (connexons), each of which is formed by integral membrane proteins called connexins (for review see Kumar and Gilula, 1996; Evans and Martin, 2002; Goodenough and Paul, 2009; Nielsen et al., 2012). Emerging evidence indicates that electrical synapses are complex, functionally diverse, and highly modifiable structures that govern the synchrony of rhythmic activity and the timing of spikes in coupled neuronal networks (Destexhe et al., 1998; Chow and Kopell, 2000; Haas and Landisman, 2012; Sevetson and Haas, 2015; Pernelle et al., 2018; Pham and Haas, 2018, 2019). Furthermore, electrical synapses can act as low-pass filters and transfer spikes after hyperpolarization, thereby inducing efficient spike-dependent depression (Vaughn and Haas, 2022).

The functional richness of the cerebellum

The cerebellum plays a well-established role in motor coordination, and recent studies have established that it also contributes to higher-order cognitive, emotional, and perceptual processing (Schmahmann, 2010; Bastian and Haslam, 2011; D’Angelo and Casali, 2012; van Beugen et al., 2013; Koziol et al., 2014; Baumann et al., 2015; D’Angelo, 2019). Given its extremely regular cytoarchitecture and the distinctive morphological and electrophysiological profiles of neurons and synapses (Hansel et al., 2001; D’Angelo and De Zeeuw, 2009), computational and mathematical approaches are more tractable for the cerebellum than for other brain areas (De Zeeuw et al., 2021). However, recent studies have suggested that cerebellar organization and properties are more complex than initially thought, indicating an impressive diversity of molecular (Tang et al., 2017), cellular (Nedelescu and Abdelhack, 2013; Nedelescu et al., 2018), and circuit mechanisms (Steuber et al., 2007; Engbers et al., 2013; Cerminara et al., 2015; Zhou et al., 2015; Kozareva et al., 2021).

Realistic multi-compartmental single-cell models

Cerebellar physiology, covering neuronal elements, synaptic connections, and circuit organization, has been intensively studied in the vertebrate brain, allowing the generation of biologically detailed models at multiple levels of complexity (D’Angelo et al., 2016). These models are based on the biophysical and anatomical properties of real neurons, including three-dimensional (3D) spatial morphology, voltage- and ligand-gated ion channels, intracellular calcium dynamics, cell densities, synaptic connectivity patterns, and cytoplasmic processes (Koch, 1998; De Schutter, 2001; D’Angelo and Casali, 2012). The most recent realistic multi-compartmental models of cerebellar neurons and synapses have been developed by integrating knowledge from various sources and capturing a rich variety of neuronal dynamic behaviors during complex neurotransmission patterns (D’Angelo et al., 2001; Solinas et al., 2007a,b; Diwakar et al., 2009; Subramaniyam et al., 2014; Masoli et al., 2015, 2020a,b; Florimbi et al., 2016, 2018, 2021; Masoli and D’Angelo, 2017; Torti et al., 2019; Rizza et al., 2021) (Box 2).

BOX 2. The NEURON simulator and the brain scaffold builder.

NEURON is a simulation environment for creating electrical and chemical signaling models of neurons and networks of neurons that are closely related to experimental data (Hines and Carnevale, 2001, 2008; Hines et al., 2009). NEURON was originally designed to facilitate the development of neuronal models, in which complex intrinsic properties and extensive cell configurations play significant roles. This includes facilities for characterizing the longitudinal distribution of ion concentrations and representing the connection network in a computationally efficient manner. The most well-known example of syntax in NEURON is the concept of continuous cable sections (analogous to an unbranched axon), which can be combined to form a tree-shaped structure (Yihe and Timofeeva, 2020). A section can be endowed with membrane properties (including channels, synapses, and ionic concentrations) that change continually with the position along its length. The primary purpose is to determine the physical properties of a neuron without regard for the numerical issue of the size of spatial compartments and, thus, to help the user focus on neuronal anatomy rather than computational details (e.g., cable equation) (Hines et al., 2009).

NEURON programming has been performed using a high-order calculator (HOC), an interpreted language with C-like syntax that has been expanded to include a domain-specific language called the “NEURON Model Description Language” (NMODL) (Hines and Carnevale, 1997, 2001). NMODL allows the incorporation of a wide range of membrane mechanisms such as voltage- and ligand-gated ion channels and active transport, and point processes such as synapses. More recently, Python was adopted as an alternative language, making available within NEURON a very extensive suite of libraries, such as BluePyOpt for ionic conductance (Van Geit et al., 2016) and Netpyne for biological neuronal networks (Dura-Bernal et al., 2019).

The brain scaffold builder (BSB) (De Schepper et al., 2021) is a framework for reconstructing and simulating neuronal network models with single neurons that have been previously reconstructed using NEURON. The scaffold modeling workflow consists of three main phases: configuration, reconstruction, and simulation. The core concepts of the framework during the reconstruction phase are (i) the network volume; (ii) the cell types, which determine the properties of cell populations, such as their spatial representation (soma radius, geometrical extension, and/or morphologies) and density information; (iii) the placement of said cell types into subspaces of the network volume; and (iv) the connectivity between cell types using certain connection strategies. With this information, the framework places and connects cells and allows the simulation of various functional conditions. This last point directly addresses synapses, which are the object of this review.

The granular layer

The granular layer is the input layer of the cerebellar cortex. It receives mossy fibers, one of the major inputs to the cerebellum, originating from multiple brainstem nuclei, and is composed of different types of cells. The most numerous are small and closely packed granule cells (GrCs), which emit four short dendrites on average and a single axon ascending vertically into the molecular layer, where they bifurcate to form parallel fibers. A few types of interneurons, including excitatory unipolar brush cells (UBCs), inhibitory Golgi cells (GoCs), and Lugaro cells, are located in this layer. GoC axons are profusely arborized and terminate on GrC dendrites, together with mossy fibers, in a peculiar synaptic structure called a “glomerulus” (Figure 2; Mugnaini et al., 1997, 2011; Voogd and Glickstein, 1998).

FIGURE 2

FIGURE 2

Models of synaptic transmission at cerebellar synapses. Schematic representation of the main elements of the cerebellar circuit along with exemplar simulations of repetitive neurotransmission (EPSC and IPSCs). The cerebellar cortex receives two excitatory inputs: mossy fibers (mf) and climbing fibers (cf). The latter originates in the inferior olive (IO). Both inputs send collaterals to DCN before entering the cerebellar cortex. Mfs contact GrC and GoC. Cfs, contact PC in the purkinje cell layer (PCL). The ascending axons (aa) of the GrCs reach the molecular layer (ML) where it bifurcates into parallel fibers (pf), which synapse onto PC and ML interneurons [stellate cells (SC) and basket cell (BCs)]. The activity of PCs is under inhibitory control by SCs and BCs. The only output of the cortex is provided by PCs, which project to the DCN and provide the main output of the cerebellar circuit. The DCNs project to the IO and in turn to PCs generating a loop mediated by the cfs [modified from D’Angelo et al. (2016)]. The simulated traces are taken from the following papers: mf-GrC, (Nieus et al., 2006), pf-SC (Rizza et al., 2021), GoC-GrC (Nieus et al., 2014), pf-PC (Masoli and D’Angelo, 2017), mf-GoC and GrC-GoC (Masoli et al., 2020a). All the traces reproduce 100 Hz activity. Note the short-term depression at the mf-GrC synapse and short-term facilitation at the pf-PC and pf-SC synapses, while GoC-GrC and mf-GoC demonstrate little temporal dynamics.

The properties of mossy fiber-GrC synapses have been the subject of extensive electrophysiological in vitro studies (Delvendahl and Hallermann, 2016). Specifically, synaptic currents (EPSCs) have two distinct kinetic components: an early fast AMPA receptor-mediated component and a later longer-lasting NMDA receptor-mediated component (Silver et al., 1992; D’Angelo et al., 1999; DiGregorio et al., 2002; Sola et al., 2004; Nieus et al., 2006; Baade et al., 2016). The NMDA receptor subunit composition undergoes a developmental switch from GluN2A/B to GluN2C-containing receptors (Farrant et al., 1994; Monyer et al., 1994). The GrC spiking output is enhanced in a non-linear manner by membrane depolarization owing to voltage-dependent NMDA receptor activation (D’Angelo et al., 1995; Cathala et al., 2000; Imamura et al., 2000; Rothman et al., 2009; Schwartz et al., 2012). Other factors, such as glutamate spillover, synaptic plasticity, and inhibition (e.g., feed-forward, feedback, and lateral) also affect the integration of GrC information (Mapelli and D’Angelo, 2007; Kanichay and Silver, 2008; Hallermann and Silver, 2013). By combining data collected from electrophysiological experiments and biophysical analyses in vitro (D’Angelo et al., 2001; Nieus et al., 2006; Solinas et al., 2007a,b, 2010; Diwakar et al., 2009) and in vivo (Kase et al., 1980; Chadderton et al., 2004; Jörntell and Ekerot, 2006; Rancz et al., 2007), large-scale models of the granular layer have been used to predict circuit spatiotemporal coding properties (Maex and De Schutter, 1998; Medina and Mauk, 2000; Solinas et al., 2010; Kalmbach et al., 2011). Interestingly, the realistic computational model of the granular layer network dynamics developed by Solinas et al. (2010) established different roles for the main components of inhibition. Lateral inhibition is crucial for determining center-surface effects [GrCs in the core are more activated than those in the surrounding area (Mapelli and D’Angelo, 2007; D’Angelo and Casali, 2012; Prestori et al., 2019)], feed-forward inhibition to control the temporal window during which GrCs integrate excitatory inputs (D’Angelo and De Zeeuw, 2009; D’Angelo and Casali, 2012; Gandolfi et al., 2014; Eshra et al., 2019; Masoli et al., 2020b), and feedback inhibition to favor coherent oscillations in the theta frequency band (D’Angelo and Casali, 2012; Mapelli et al., 2014).

In contrast to the canonical view describing GrCs as a homogenous population of neurons generating regular firing, our recently published work combining data from simulations and experimental findings on GrC electro responsive properties has identified different subtypes: adapting, non-adapting and accelerating GrCs (Masoli et al., 2020b). Adaptation occurs via Ca2+ influx through high-threshold Ca2+ conductance and subsequent increase in Ca2+-activated K+ currents, causing an after-hyperpolarization potential after a burst of spikes (D’Angelo et al., 1998) whereas acceleration is correlated with TRPM4 channels, which are primarily described in PCs (Kim et al., 2009, 2013; Menigoz et al., 2016). Remarkably, our model, derived from previous models (D’Angelo et al., 2001; Nieus et al., 2006; Diwakar et al., 2009) and upgraded to account for sophisticated mechanisms of spike generation, predicted that TRPM4 currents, first recorded in GrCs, coupled to calmodulin through intracellular Ca2+ changes, could effectively generate firing acceleration observed experimentally (Masoli et al., 2020b). Remarkably, mossy fiber-GrC synapse simulations showed that fine-tuning of adaptation/acceleration and short-term plasticity generated effective mossy fiber-GrC transmission channels, which could enrich the processing of incoming spike trains and enhance spatiotemporal coding at the cerebellar input stage (Gabbiani et al., 1994), according to the theoretical prediction of the adaptive filter model (Marr, 1969; Dean and Porrill, 2011; Rössert et al., 2015; D’Angelo et al., 2016).

Neurotransmission mechanisms at the mossy fiber-GrC synapse have been rigorously investigated, providing a workbench for synaptic investigations in brain slices (D’Angelo et al., 1990; Silver et al., 1992; Mitchell and Silver, 2003) and accurate neurotransmission models, including spillover (DiGregorio et al., 2002) and single-channel opening properties. The AMPA current transient is ultra-fast and controls the high-fidelity millisecond-precise transmission, while a slow spillover-driven persistent component contributes to regulate bursting. The release probability was set to 0.43 to ensure a marked short-term depression, but these synapses were observed to range from 0.4 to 0.6/0.8 (D’Errico et al., 2009), causing adapting discharges in response to mossy fiber bursts. The release probability increases with LTP and decreases with long-term depression (LTD), tuning the first spike delay over a 100 ms window (Nieus et al., 2006). The NMDA current changes during development (Cathala et al., 2000) and can be activated through spillover and non-linearly boosts GrC spikes, generating burst-burst (detonator) retransmission along the mossy fiber-GrC pathway (D’Angelo et al., 1993, 1995). In combination with these specialized neurotransmission properties, the richness of GrC intrinsic electro-responsiveness can produce a variety of transmission patterns, causing spatiotemporal reconfiguration of incoming patterns in the granular layer (D’Angelo and De Zeeuw, 2009; Eshra et al., 2019; Casali et al., 2020) and determining the adaptive filtering properties of the cerebellum.

Since their discovery (Golgi, 1874), GoCs have been intensively studied, both experimentally and theoretically (Dieudonné, 1998; Maex and De Schutter, 1998; De Schutter et al., 2000; Geurts et al., 2003; Forti et al., 2006; Solinas et al., 2007a,b; Hull and Regehr, 2012; Cesana et al., 2013). In recent years, electrophysiological measurements have revealed some important aspects of the cellular functions of these neurons. In particular, GoCs have a rich and complex variety of intrinsic electro-responsive properties, including pacemaker activity, resonance in the theta band frequency, and phase resetting (Dieudonné, 1998; Forti et al., 2006; Solinas et al., 2007a,b). Subsequently, these properties have been modeled using specific ionic channels, excitatory and inhibitory synapses, and GJ (Solinas et al., 2007a,b). We recently developed detailed multi-compartmental models of GoCs that faithfully capture the large repertoire of findings from electrophysiological experiments (Masoli et al., 2020a). The principal prediction of these models is that the synaptic activation of apical dendrites through parallel fibers causes slow Ca2+-dependent depolarizations that trigger and, more importantly, regulate the number of spikes that back-propagate into basal dendrites with concomitant NMDA receptor unblock (when the mossy fiber synapses onto basal dendrites are active). This dendritic processing can be considered a coincidence detector that controls the availability of NMDA currents in basal dendrites, providing the basis for spike-time-dependent plasticity (STDP) anticipated by theory (Garrido et al., 2016).

UBCs form a complex network with mossy fibers (Balmer and Trussell, 2019), GoCs, GrCs, and other UBCs (Mugnaini et al., 1997). A biologically multi-compartmental model reproduces the intrinsic UBC electroresponsiveness, confirming the key role of H- and Ca2+ currents (Subramaniyam et al., 2014). The model accurately predicted the generation of late-onset responses following mossy fiber stimulation and its dependence on the intracellular cAMP cascade and the consequent regulation of TRP currents (Locatelli et al., 2013). Synaptic transmission properties are also controlled by AMPA receptors (Lu et al., 2017), spillover (Balmer et al., 2021), and metabotropic receptor regulation of intrinsic excitable mechanisms (Guo et al., 2021).

The molecular layer

The outer layer of the cerebellar cortex is a molecular layer. This layer contains massive dendritic trees of PCs and parallel fibers, which travel in a parasagittal direction after bifurcating at the end of the GrC ascending axon. The dendritic trees of PCs are penetrated by thousands of parallel fibers, forming glutamatergic synapses on the spines of at least 400 PCs. These neurons are one of the most complex in the CNS and are endowed with several voltage-dependent ionic channels (Akemann and Knöpfel, 2006; Masoli et al., 2015) and numerous proteins that buffer calcium entry (Anwar et al., 2012). The molecular layer also contains two main inhibitory interneuron types, stellate cells (SCs) and basket cells (BCs). BCs are typically located in the lower third of the molecular layer and form GABAergic synapses on PC bodies (Sultan and Bower, 1998), whereas SCs are usually considered to reside in the upper two-thirds of the molecular layer and form GABAergic synapses on PC dendrites. Both BCs and SCs receive glutamatergic input from parallel fibers, as well as GABAergic input from other BCs and SCs (Figure 2; Palay and Chan-Palay, 1974). There is an open discussion regarding the molecular distinction of interneurons into subtypes, and a third class, the candelabrum cell, has recently been discovered (Kozareva et al., 2021; Osorno et al., 2022). Although specific computational roles have been proposed for GoCs, inhibitory interneurons are located in the granular layer (for review see Galliano et al., 2010), suggesting that the functions of MLIs have been confined to the simplistic idea that they provide feed-forward inhibition to PCs (Eccles et al., 1966; Ito, 1984; Mittmann et al., 2005; Santamaria et al., 2007; Bower, 2010).

Following the initial simplified models (Lennon et al., 2014, 2015), recently, the electro-responsive properties and synaptic dynamics of SCs have been characterized in vitro, allowing the identification of a set of electrophysiological properties that have been used to generate multi-compartmental SC models (Rizza et al., 2021). Specifically, SCs have pacemaker activity, spike frequency that increases almost linearly with current intensity, sagging inward rectification in response to hyperpolarizing current injection, and rebound excitation and pause following the depolarizing current step at the end of the hyperpolarization. Moreover, SCs showed marked short-term facilitation during repetitive parallel fiber transmission. Because of synaptic filtering combined with intrinsic electro-responsiveness, spikes were generated by the SC model after a lag and only at high frequencies, making SCs operate as delay-high-pass filters (Rizza et al., 2021). Our models, based on detailed morphologies and precise membrane mechanisms, faithfully reproduced the entire set of available experimental data and allowed us to explore the different functional configurations of the parallel fiber-SC-PC circuit. Interestingly, simulations predicted a new role for SCs, which could provide the molecular layer low-pass and band-pass filtering properties recoding parallel fiber bursts and regulating the PC gain in a frequency-dependent manner (Rizza et al., 2021). These findings expand the adaptive filter view of cerebellar circuit functions (Dean and Porrill, 2011).

Metabotropic receptor

Metabotropic receptors, rather than generating postsynaptic currents, exert a relatively slow effect on ionic channels, since various enzymes transduce the receptor signal through second messengers and cytoplasmic cascades. This multistep process is critical for fine-tuning slow responses, development, and cellular survival (Yamasaki et al., 2021). In the cerebellum, all mGluR subunits are expressed (except for mGluR3 and mGluR6). The most common is mGluR1, which is prominent in PCs but is also expressed in other cerebellar neurons. Some subtypes of GoC are endowed with mGluR2 and mGluR5, while mGluR4 is restricted to parallel fibers (Knöpfel and Grandes, 2002). This receptor mediates slow electrical responses through TRPC channels that, in the case of PCs, are critical for the development and maintenance of the dendritic structure (Wu and Kapfhammer, 2022).

The GABA(B) receptor can be found in specific locations in the cerebellum (Kulik et al., 2002), especially in PC spines and dendritic shafts, and parallel fibers (Luján et al., 2018). The receptor can cluster with Kir3.x ionic channels and the main PC calcium channel Cav2.1. The presence of both GABA(B) and NMDA receptors in PCs has been suggested to support bistability (Sanders et al., 2013) induced by the activation of climbing fibers (Loewenstein et al., 2005). GABA(B) receptors are also expressed in GrCs, where they modulate the input-output relationship (Mapelli et al., 2014; Bassetti, 2022). A computational model of GABA(B) and Kir3 ionic channels has been built to reproduce physiological responses (Destexhe et al., 1994).

Electrical synapses and ephaptic coupling

GoCs are electrically coupled via GJs, which enhance low-frequency oscillatory synchronization and resonance in GoC networks (Dugué et al., 2009). Asynchronous excitatory inputs to synchronized GoCs cause rapid and considerable desynchronization (Vervaeke et al., 2010). Few studies, combining experiments and modeling, have investigated how GJ coupling affects the GoC and granular layer oscillation using isolated GoC networks (Dugué et al., 2009; Vervaeke et al., 2010) or including the parallel fiber feedback loop. A more realistic computational model with GoCs synaptically connected in the granular layer (Simões de Souza and De Schutter, 2011) confirmed previous findings (Dugué et al., 2009; Vervaeke et al., 2010), suggesting that GJs between GoCs increased the power of cerebellar cortex oscillations, which were synaptically driven by the feedback loop between GoCs and GrCs. These results suggest a novel function for GJs between GoCs in determining a change in the timing of oscillatory cycles, which subsequently improves GoC synchronization (Simões de Souza and De Schutter, 2011). GoC-GoC GJs have recently been introduced into a model of the cerebellar cortex (De Schepper et al., 2021).

GJs are electrical synapses that have advantages over chemical synapses, such as being instantaneous and having no transmission failures (Hoehne et al., 2020). They can transmit signals in both directions; however, through synaptic rectification, they can transmit signals in a preferred direction only (Snipas et al., 2017). In the cerebellum, SCs form GJs with just another SC to form couples (Alcami and Marty, 2013), while BCs tend to form GJs with multiple neighboring BCs (Kim et al., 2014). The physiological consequences of this organization can be monitored using the different responses elicited by PCs. SCs form synapses on the vast dendritic tree of PCs and have a modulatory effect on spontaneous firing, generating pauses that reach a maximum of 60–70 ms (Steuber et al., 2007). BCs form GABAergic synapses on the soma of the PC and ephaptic connections on the axon initial segment (AIS) (Blot and Barbour, 2014). A PC is in contact with up to seven BC axons, but only one or two of them form functional synapses (Kim et al., 2014). The consequences of this synaptic organization propagate to the deep cerebellar nuclei (Chu et al., 2012).

Ephaptic coupling, a mechanism by which neurons communicate via electrical signals through the extracellular space, was first observed between spontaneously active fibers in mouse spinal nerve roots (Rasminsky, 1980). Shortly after, ephaptic signaling was shown to control synchronization and time-spiking in hippocampal neurons (Taylor and Dudek, 1982). More recently, distinctive characteristics of ephaptic interactions have been observed in both cerebral and cerebellar cells (Anastassiou et al., 2011, 2015; Blot and Barbour, 2014; Han et al., 2018). In the cerebellum, ephaptic coupling has been shown to occur between the axon terminal of BCs (pinceau) and the PC somata (Blot and Barbour, 2014), and between PCs near the AISs (Bender and Trussell, 2012; Han et al., 2018). Ephaptic coupling between PCs differs from that provided by BC synapses onto PCs. In the latter case, ephaptic transmission has an inhibitory effect. The extracellular field produced by the BC pinceau is characterized by two current components, a capacitive current and an ionic current through voltage-activated potassium conductances (the pinceau is devoid of sodium channels) (Laube et al., 1996; Southan and Robertson, 2000; Bobik et al., 2004) resulting in a rapid inhibition of PC firing. In contrast, for PC AISs, the extracellular field is dominated by sodium current flow, producing a large excitatory signal and subsequently promoting the synchronization of nearby PCs. These examples reveal the importance of the distribution of ionic channels within a cell in determining the strength and sign of ephaptic coupling. To date, no study has reproduced these effects on PCs using computational models.

A summary view

The most critical aspects of synaptic transmission in the cerebellar cortical circuit can be summarized by the parameters of the vesicle release cycle, neurotransmitter diffusion, and postsynaptic receptor activation (Figure 1 and Table 1). A comparative analysis of these parameters, which were extracted by fitting EPSC and IPSC trains at the most important cerebellar cortical synapses, revealed that excitatory synapses formed by ascending axons and parallel fibers emitted by GrCs usually had a lower release probability than others in the circuit. This implies that these synapses mostly show short-term facilitation, while the others show short-term depression. Therefore, the activation of molecular layer interneurons (MLIs) and PCs is subjected to heavy filtering, so that single spikes in mossy fibers do not have remarkable effects, while bursts emitted by GrCs can pass through. This single-synapse filtering is integrated at the cell level by intrinsic excitability and converted into effective circuit filtering by the network neuronal chains. A clear example is that of the MLI-PC sub-circuit in the molecular layer: parallel fiber-MLI synapses facilitate the generation of a high-pass filter in SCs, which then inhibits PCs, causing low-pass filtering of parallel fiber-PC transmissions (Rizza et al., 2021). In the granular layer, the high mossy fiber-GrC release probability combines spillover and non-linear NMDA receptor activation, generating a high-pass filter controlling the phase and composition of burst emission by GrCs. This circuit is further regulated by the inhibitory synapses of the GoCs.

These observations imply that burst-burst retransmission is at the core of the cerebellar coding scheme. Bursts in mossy fibers are initially transformed and retransmitted in the granular layer and are then used to tune the firing frequency of PCs in the molecular layer circuit. Moreover, input bursts prime impulsive or repetitive burst-pause patterns in PCs, SCs, and GoCs, sustained by delays in the feed-forward and feedback circuit inhibitory loops and amplified by the intrinsic electroresponsive properties of neurons. In this system, specific synapse parameters play a key role in tuning transmission and filtering properties.

Another critical factor is neurotransmitter spillover, which sustains bursting at the granular and molecular layer synapses. Moreover, the properties of specific AMPA receptor subtypes (gating kinetics and desensitization) determine the precision of spike emission, whereas the properties of NMDA receptor subunits (regulating gating kinetics, desensitization, and voltage-dependent Mg block) determine the intensity and duration of bursts. Another dimension is introduced by cell-to-cell variability, with differentiated excitable and neurotransmission properties even within the same class of neurons, [e.g., GrCs (Chabrol et al., 2015; Gilmer and Person, 2017; Masoli et al., 2020b; Straub et al., 2020), MLIs (Rizza et al., 2021), and candelabrum cells (Osorno et al., 2022), although not included in this work yet]. Finally, the strategic location of NMDA receptors at certain synapses (mossy fiber-GrC, mossy fiber-GoC, parallel fiber-PC, parallel fiber-MLI), as well as Ca permeation in certain types of AMPA receptors and voltage-gated calcium channels, causes various forms of plasticity that can further tune postsynaptic receptor conductance and presynaptic release probability, regulating the filtering properties of the circuit and local computations.

Conclusion and future challenges

In conclusion, the properties of neurotransmission and intrinsic electroresponsiveness determined experimentally in the cerebellum have been translated into accurate biophysical models, which, in turn, provide a new view of how the cellular and subcellular organization of neurons and synapses contribute to function. Indeed, while a bewildering number of molecular properties await functional explanation and system integration, simulations of detailed computational models can suggest new directions for experimental investigations and provide hints for understanding neuropathologies (D’Angelo and Jirsa, 2022). For example, alterations in synaptic activity in brain diseases affecting the cerebellum, such as ataxia, dystonia, autism, and paroxysmal kinesigenic dyskinesia (LeDoux et al., 1998; Rinaldo and Hansel, 2010; Brown and Loew, 2012; Won et al., 2013; Power and Empson, 2014; Zhang et al., 2015; Zhang and Sudhof, 2016; White and Sillitoe, 2017; Soda et al., 2019; Binda et al., 2021), remain largely to be elucidated; therefore, their translation into models would require investigation of structural changes, receptors, and biochemical pathways. This entails extending the models to include the mechanisms of metabotropic receptor-dependent neuromodulation and calcium-dependent long-term synaptic plasticity (either in the form of LTP, LTD, or STDP). This, in turn, implies updating and expanding the simulation environment, for example, incorporating a voxelized extracellular space into the simulators, as was recently done in NEURON with the reaction-diffusion module (RxD) (Newton et al., 2018). This module can be used to simulate complex biochemical pathways and diffuse substances into a 3D space. This approach allows for the reconstruction of more detailed synapses with substance diffusion between neurons. Unavoidably, this will bring about an increased computational cost, promoting the development of new modeling paradigms and simulation environments (Chen et al., 2022), implying that brain modeling will depend increasingly on the intersection of biology, computation, and informatics.

Here, we indicated that specific properties of neurotransmission and intrinsic electroresponsiveness shape the spatiotemporal profiles of cerebellar cortical circuit activity, as predicted by motor learning theory and adaptive filter theory (Marr, 1969; Dean and Porrill, 2011; Rössert et al., 2015; D’Angelo et al., 2016). In addition to improving single-neuron and synaptic models, the challenge is threefold. First, the remaining neurons and synapses of the cerebellar network, including unipolar brush, Lugaro, and candelabrum cells (Guo et al., 2021; Miyazaki et al., 2021; Osorno et al., 2022), tune cerebellar circuit activity through recurrent inhibitory loops. Second, these new cellular models and synapses are included in an updated version of the cerebellar cortical model (De Schepper et al., 2021) and then extended to the mesoscale by including the deep cerebellar nuclei (Uusisaari et al., 2007; Uusisaari and Knöpfel, 2011) and inferior olives (Torben-Nielsen et al., 2012). Third, we simulated dynamic activity in the mesoscale network models. Simulations should be paralleled by experimental investigation of synaptic and cellular properties to better understand signal integration and the ensuing dynamics, providing essential information for model construction and validation. The existing models contain experimental data available at the time of construction, so they can be expanded with newer information as they become available. At the same time, missing cerebellar models, such as candelabrum cells (Osorno et al., 2022), have recently seen a surge in interest that can induce other labs to further study them. The missing data for both known and less known neurons can be optimally obtained using advanced recording techniques, including high-density Multi electrode array (Buccino et al., 2022), voltage-sensitive dye imaging and multispot 2-photon laser microscopy (Ebner and Chen, 1995; Gandolfi et al., 2014; Miccoli et al., 2019; Gagliano et al., 2022), as well as tools to target and dissect specific circuit components, such as genetic engineering and optogenetics (Yang et al., 2018; Adam et al., 2019; Tognolina et al., 2022). We envisage that this roadmap will eventually shed new light on the computational implications of synaptic mechanisms and contribute to explaining the intrinsic algorithms of the cerebellar circuit.

Statements

Author contributions

SM, MR, MT, and FP organized and wrote the manuscript and prepared the figures. ED’A revised and contributed to the final version of the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work received funding from the European Union’s Horizon 2020 Framework Programme for Research and Innovation (Specific Grant Agreement No. 945539) (Human Brain Project SGA3) and (Specific Grant Agreement No. 785907) (Human Brain Project SGA2).

Acknowledgments

We thank Editage (www.editage.com) for English language editing.

Conflict of interest

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.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AdamY.KimJ. J.LouS.ZhaoY.XieM. E.BrinksD.et al (2019). Voltage imaging and optogenetics reveal behaviour-dependent changes in hippocampal dynamics.Nature569413417. 10.1038/s41586-019-1166-7

  • 2

    AkarN. A.CummingB.KarakasisV.KustersA.KlijnW.PeyserA.et al (2019). “Arbor — A Morphologically-Detailed Neural Network Simulation Library for Contemporary High-Performance Computing Architectures,” in 2019 27th Euromicro International Conference on Parallel, Distributed and Network-Based Processing (PDP), (Piscataway: IEEE), 274282. 10.1109/EMPDP.2019.8671560

  • 3

    AkemannW.KnöpfelT. (2006). Interaction of Kv3 potassium channels and resurgent sodium current influences the rate of spontaneous firing of Purkinje neurons.J. Neurosci.2646024612. 10.1523/JNEUROSCI.5204-05.2006

  • 4

    AlahmadiA. A. S.PardiniM.SamsonR. S.D’AngeloE.FristonK. J.ToosyA. T.et al (2015). Differential involvement of cortical and cerebellar areas using dominant and nondominant hands: An FMRI study.Hum. Brain Mapp.3650795100. 10.1002/hbm.22997

  • 5

    AlcamiP.MartyA. (2013). Estimating functional connectivity in an electrically coupled interneuron network.Proc. Natl. Acad. Sci. U.S.A.110:E4798E4807. 10.1073/pnas.1310983110

  • 6

    AnastassiouC. A.PerinR.BuzsákiG.MarkramH.KochC. (2015). Cell type- and activity-dependent extracellular correlates of intracellular spiking.J. Neurophysiol.114608623. 10.1152/jn.00628.2014

  • 7

    AnastassiouC. A.PerinR.MarkramH.KochC. (2011). Ephaptic coupling of cortical neurons.Nat. Neurosci.14217223. 10.1038/nn.2727

  • 8

    AndreaeL. C.BurroneJ. (2014). The role of neuronal activity and transmitter release on synapse formation.Curr. Opin. Neurobiol.274752. 10.1016/j.conb.2014.02.008

  • 9

    AnwarH.HongS.De SchutterE. (2012). Controlling Ca 2+-activated K + channels with models of Ca 2+ buffering in purkinje cells.Cerebellum11681693. 10.1007/s12311-010-0224-3

  • 10

    BaadeC.ByczkowiczN.HallermannS. (2016). NMDA receptors amplify mossy fiber synaptic inputs at frequencies up to at least 750 Hz in cerebellar granule cells.Synapse70269276. 10.1002/syn.21898

  • 11

    BalmerT. S.Borges-MerjaneC.TrussellL. O. (2021). Incomplete removal of extracellular glutamate controls synaptic transmission and integration at a cerebellar synapse.eLife10:e63819. 10.7554/eLife.63819

  • 12

    BalmerT. S.TrussellL. O. (2019). Selective targeting of unipolar brush cell subtypes by cerebellar mossy fibers.eLife8:e44964. 10.7554/eLife.44964

  • 13

    BarbourB. (2001). An evaluation of synapse independence.J. Neurosci.2179697984. 10.1523/jneurosci.21-20-07969.2001

  • 14

    BassettiD. (2022). Keeping the balance: GABAB receptors in the developing brain and beyond.Brain Sci.12:419. 10.3390/brainsci12040419

  • 15

    BastianB.HaslamN. (2011). Experiencing dehumanization: Cognitive and emotional effects of everyday dehumanization.Basic Appl. Soc. Psychol.33295303. 10.1080/01973533.2011.614132

  • 16

    BatoolS.RazaH.ZaidiJ.RiazS.HasanS.SyedN. I. (2019). Synapse formation: From cellular and molecular mechanisms to neurodevelopmental and neurodegenerative disorders.J. Neurophysiol.12113811397. 10.1152/jn.00833.2018

  • 17

    BaumannO.BorraR. J.BowerJ. M.CullenK. E.HabasC.IvryR. B.et al (2015). Consensus paper: The role of the cerebellum in perceptual processes.Cerebellum14197220. 10.1007/s12311-014-0627-7

  • 18

    BaumannS. W.BaurR.SigelE. (2001). Subunit arrangement of gamma-aminobutyric acid type A receptors.J. Biol. Chem.2763627536280. 10.1074/jbc.M105240200

  • 19

    BaumannS. W.BaurR.SigelE. (2002). Forced subunit assembly in alpha1beta2gamma2 GABAA receptors. Insight into the absolute arrangement.J. Biol. Chem.2774602046025. 10.1074/jbc.M207663200

  • 20

    BenderK. J.TrussellL. O. (2012). The physiology of the axon initial segment.Annu. Rev. Neurosci.35249265. 10.1146/annurev-neuro-062111-150339

  • 21

    BennaniO.ChauvetG.ChauvetP.DupontJ. M.JouenF. (2009). A hierarchical modeling approach of hippocampus local circuit.J. Integr. Neurosci.84976. 10.1142/S0219635209002083

  • 22

    BentivoglioM.CotrufoT.FerrariS.TesorieroC.MariottoS.BertiniG.et al (2019). The original histological slides of camillo golgi and his discoveries on neuronal structure.Front. Neuroanat.13:2. 10.3389/fnana.2019.00003

  • 23

    BenvenisteM.ClementsJ.VyklickýL.MayerM. L. (1990). A kinetic analysis of the modulation of N-methyl-D-aspartic acid receptors by glycine in mouse cultured hippocampal neurones.J. Physiol.428333357. 10.1113/jphysiol.1990.sp018215

  • 24

    BindaF.ValenteP.MarteA.BaldelliP.BenfenatiF. (2021). Increased responsiveness at the cerebellar input stage in the PRRT2 knockout model of paroxysmal kinesigenic dyskinesia.Neurobiol. Dis.152:105275. 10.1016/j.nbd.2021.105275

  • 25

    BlotA.BarbourB. (2014). Ultra-rapid axon-axon ephaptic inhibition of cerebellar Purkinje cells by the pinceau.Nat. Neurosci.17289295. 10.1038/nn.3624

  • 26

    BobikM.EllismanM. H.RudyB.MartoneM. E. (2004). Potassium channel subunit Kv3.2 and the water channel aquaporin-4 are selectively localized to cerebellar pinceau.Brain Res.1026168178. 10.1016/j.brainres.2004.07.088

  • 27

    BowerJ. M. (2010). Model-founded explorations of the roles of molecular layer inhibition in regulating purkinje cell responses in cerebellar cortex: More trouble for the beam hypothesis.Front. Cell. Neurosci.4:27. 10.3389/fncel.2010.00027

  • 28

    BowerJ. M. (2015). The 40-year history of modeling active dendrites in cerebellar Purkinje cells: Emergence of the first single cell “community model”.Front. Comput. Neurosci.9:129. 10.3389/fncom.2015.00129

  • 29

    BowerJ. M.BeemanD. (2007). Constructing realistic neural simulations with GENESIS.Methods Mol. Biol.401103125. 10.1007/978-1-59745-520-6_7

  • 30

    BrodlandG. W. (2015). How computational models can help unlock biological systems.Semin. Cell Dev. Biol.47–486273. 10.1016/j.semcdb.2015.07.001

  • 31

    BrownS.-A.LoewL. M. (2012). Computational analysis of calcium signaling and membrane electrophysiology in cerebellar Purkinje neurons associated with ataxia.BMC Syst. Biol.6:70. 10.1186/1752-0509-6-70

  • 32

    BuccinoA. P.DamartT.BartramJ.MandgeD.XueX.ZbiliM.et al (2022). A multi-modal fitting approach to construct single-neuron models with patch clamp and high-density microelectrode arrays.bioRxiv.[Preprint]. 2022.08.03.502468.

  • 33

    CakanC.ObermayerK. (2020). Biophysically grounded mean-field models of neural populations under electrical stimulation.PLoS Comput. Biol.16:e1007822. 10.1371/journal.pcbi.1007822

  • 34

    CannonR. C.GleesonP.CrookS.GanapathyG.MarinB.PiasiniE.et al (2014). LEMS: A language for expressing complex biological models in concise and hierarchical form and its use in underpinning NeuroML 2.Front. Neuroinform.8:79. 10.3389/fninf.2014.00079

  • 35

    CasaliS.TognolinaM.GandolfiD.MapelliJ.D’AngeloE. (2020). Cellular-resolution mapping uncovers spatial adaptive filtering at the rat cerebellum input stage.Commun. Biol.3:635. 10.1038/s42003-020-01360-y

  • 36

    CathalaL.MisraC.Cull-CandyS. (2000). Developmental profile of the changing properties of NMDA receptors at cerebellar mossy fiber-granule cell synapses.J. Neurosci.2058995905. 10.1523/JNEUROSCI.20-16-05899.2000

  • 37

    CaveroI.GuillonJ. M.HolzgrefeH. H. (2017). Reminiscing about Jan Evangelista Purkinje: A pioneer of modern experimental physiology.Adv. Physiol. Educ.41528538. 10.1152/advan.00068.2017

  • 38

    CerminaraN. L.LangE. J.SillitoeR. V.AppsR. (2015). Redefining the cerebellar cortex as an assembly of non-uniform Purkinje cell microcircuits.Nat. Rev. Neurosci.167993. 10.1038/nrn3886

  • 39

    CesanaE.PietrajtisK.BidoretC.IsopeP.D’AngeloE.DieudonnéS.et al (2013). Granule cell ascending axon excitatory synapses onto Golgi cells implement a potent feedback circuit in the cerebellar granular layer.J. Neurosci.331243012446. 10.1523/JNEUROSCI.4897-11.2013

  • 40

    ChabrolF. P.ArenzA.WiechertM. T.MargrieT. W.DigregorioD. A. (2015). Europe PMC funders group synaptic diversity enables temporal coding of coincident multi- sensory inputs in single neurons.18718727. 10.1038/nn.3974.Synaptic

  • 41

    ChaddertonP.MargrieT. W.HäusserM. (2004). Integration of quanta in cerebellar granule cells during sensory processing.Nature428856860. 10.1038/nature02442

  • 42

    ChenW.de SchutterE. (2017). Parallel STEPS: Large scale stochastic spatial reaction-diffusion simulation with high performance computers.Front. Neuroinform.11:13. 10.3389/fninf.2017.00013

  • 43

    ChenW.HepburnI.MartyushevA.De SchutterE. (2022). “Modeling neurons in 3D at the nanoscale.Adv. Exp. Med. Biol.1359324. 10.1007/978-3-030-89439-9_1

  • 44

    ChowC. C.KopellN. (2000). Dynamics of spiking neurons with electrical coupling.Neural Comput.1216431678. 10.1162/089976600300015295

  • 45

    ChuC. P.BingY. H.LiuH.QiuD. L. (2012). Roles of molecular layer interneurons in sensory information processing in mouse Cerebellar cortex Crus ii in vivo.PLoS One7:e37031. 10.1371/journal.pone.0037031

  • 46

    ClementsJ. D.LesterR. A.TongG.JahrC. E.WestbrookG. L. (1992). The time course of glutamate in the synaptic cleft.Science25814981501. 10.1126/science.1359647

  • 47

    ClineM. S.SmootM.CeramiE.KuchinskyA.LandysN.WorkmanC.et al (2007). Integration of biological networks and gene expression data using Cytoscape.Nat. Protoc.223662382. 10.1038/nprot.2007.324

  • 48

    Cull-CandyS.BrickleyS.FarrantM. (2001). NMDA receptor subunits: Diversity, development and disease.Curr. Opin. Neurobiol.11327335. 10.1016/s0959-4388(00)00215-4

  • 49

    Cull-CandyS. G.LeszkiewiczD. N. (2004). Role of distinct NMDA receptor subtypes at central synapses.Sci. STKE2004:re16. 10.1126/stke.2552004re16

  • 50

    D’AngeloE. (2019). The cerebellum gets social.Science363:229. 10.1126/science.aaw2571

  • 51

    D’AngeloE.AntoniettiA.CasaliS.CasellatoC.GarridoJ. A.LuqueN. R.et al (2016). Modeling the cerebellar microcircuit: New strategies for a long-standing issue.Front. Cell. Neurosci.10:176. 10.3389/fncel.2016.00176

  • 52

    D’AngeloE.CasaliS. (2012). Seeking a unified framework for cerebellar function and dysfunction: From circuit operations to cognition.Front. Neural Circuits6:116. 10.3389/fncir.2012.00116

  • 53

    D’AngeloE.De FilippiG.RossiP.TagliettiV. (1995). Synaptic excitation of individual rat cerebellar granule cells in situ: Evidence for the role of NMDA receptors.J. Physiol.484397413. 10.1113/jphysiol.1995.sp020673

  • 54

    D’AngeloE.De FilippiG.RossiP.TagliettiV. (1998). Ionic mechanism of electroresponsiveness in cerebellar granule cells implicates the action of a persistent sodium current.J. Neurophysiol.80493503. 10.1152/jn.1998.80.2.493

  • 55

    D’AngeloE.De ZeeuwC. I. (2009). Timing and plasticity in the cerebellum: Focus on the granular layer.Trends Neurosci.323040. 10.1016/j.tins.2008.09.007

  • 56

    D’AngeloE.NieusT.BezziM.ArleoA.CoenenO. J. M. D. (2005). Modeling synaptic transmission and quantifying information transfer in the granular layer of the cerebellum.Lect. Notes Comput. Sci.3512107114. 10.1007/11494669_14

  • 57

    D’AngeloE.NieusT.MaffeiA.ArmanoS.RossiP.TagliettiV.et al (2001). Theta-frequency bursting and resonance in cerebellar granule cells: Experimental evidence and modeling of a slow k+-dependent mechanism.J. Neurosci.21759770. 10.1523/JNEUROSCI.21-03-00759.2001

  • 58

    D’AngeloE.RossiP.ArmanoS.TagliettiV. (1999). Evidence for NMDA and mGlu receptor-dependent long-term potentiation of mossy fiber-granule cell transmission in rat cerebellum.J. Neurophysiol.81277287. 10.1152/jn.1999.81.1.277

  • 59

    D’AngeloE.RossiP.GarthwaiteJ. (1990). Dual-component NMDA receptor currents at a single central synapse.Nature346467470. 10.1038/346467a0

  • 60

    D’AngeloE.RossiP.TagliettiV. (1993). Different proportions of N-methyl-D-aspartate and non-N-methyl-D-aspartate receptor currents at the mossy fibre-granule cell synapse of developing rat cerebellum.Neuroscience53121130. 10.1016/0306-4522(93)90290-v

  • 61

    D’AngeloE. D.JirsaV. (2022). Neurosciences The quest for multiscale brain modeling.Trends Neurosci.xx114. 10.1016/j.tins.2022.06.007

  • 62

    DayN. C.WilliamsT. L.InceP. G.KambojR. K.LodgedD.ShawP. J. (1995). Distribution of AMPA-selective glutamate receptor subunits in the human hippocampus and cerebellum.Mol. Brain Res.311732. 10.1016/0169-328X(95)00021-J

  • 63

    De SchepperR.GeminianiA.MasoliS.RizzaM. F.AntoniettiA.CasellatoC.et al (2021). Scaffold modelling captures the structure-function-dynamics relationship in brain microcircuits.bioRxiv[Preprint]. 2021.07.30.454314.

  • 64

    De SchutterE. (2001). Computational Neuroscience: Realistic Modeling for Experimentalists.Boca Raton: CRC Press.

  • 65

    De SchutterE.VosB.MaexR. (2000). The function of cerebellar Golgi cells revisited.Prog. Brain Res.1248193. 10.1016/s0079-6123(00)24009-0

  • 66

    De ZeeuwC. I.LisbergerS. G.RaymondJ. L. (2021). Diversity and dynamism in the cerebellum.Nat. Neurosci.24160167. 10.1038/s41593-020-00754-9

  • 67

    DeanP.PorrillJ. (2011). Evaluating the adaptive-filter model of the cerebellum.J. Physiol.58934593470. 10.1113/jphysiol.2010.201574

  • 68

    DelvendahlI.HallermannS. (2016). The cerebellar mossy fiber synapse as a model for high-frequency transmission in the mammalian CNS.Trends Neurosci.39722737. 10.1016/j.tins.2016.09.006

  • 69

    D’ErricoA.PrestoriF.D’AngeloE. (2009). Differential induction of bidirectional long-term changes in neurotransmitter release by frequency-coded patterns at the cerebellar input.J. Physiol.58758435857. 10.1113/jphysiol.2009.177162

  • 70

    DestexheA.ContrerasD.SteriadeM. (1998). Mechanisms underlying the synchronizing action of corticothalamic feedback through inhibition of thalamic relay cells.J. Neurophysiol.799991016. 10.1152/jn.1998.79.2.999

  • 71

    DestexheA.MainenZ. F.SejnowskiT. J. (1994). Synthesis of models for excitable membranes, synaptic transmission and neuromodulation using a common kinetic formalism.J. Comput. Neurosci.1195230. 10.1007/BF00961734

  • 72

    DiamondJ. S. (2005). Deriving the glutamate clearance time course from transporter currents in CA1 hippocampal astrocytes: Transmitter uptake gets faster during development.J. Neurosci.2529062916. 10.1523/JNEUROSCI.5125-04.2005

  • 73

    DieudonnéS. (1998). Submillisecond kinetics and low efficacy of parallel fibre-Golgi cell synaptic currents in the rat cerebellum.J. Physiol.510845866. 10.1111/j.1469-7793.1998.845bj.x

  • 74

    DiGregorioD. A.NusserZ.SilverR. A. (2002). Spillover of glutamate onto synaptic AMPA receptors enhances fast transmission at a cerebellar synapse.Neuron35521533. 10.1016/s0896-6273(02)00787-0

  • 75

    DiGregorioD. A.RothmanJ. S.NielsenT. A.SilverR. A. (2007). Desensitization properties of AMPA receptors at the cerebellar mossy fiber granule cell synapse.J. Neurosci.2783448357. 10.1523/JNEUROSCI.2399-07.2007

  • 76

    DiwakarS.LombardoP.SolinasS.NaldiG.D’AngeloE. (2011). Local field potential modeling predicts dense activation in cerebellar granule cells clusters under LTP and LTD control.PLoS One6:e21928. 10.1371/journal.pone.0021928

  • 77

    DiwakarS.MagistrettiJ.GoldfarbM.NaldiG.D’AngeloE. (2009). Axonal Na+ channels ensure fast spike activation and back-propagation in cerebellar granule cells.J. Neurophysiol.101519532. 10.1152/jn.90382.2008

  • 78

    Dos SantosR. W.AlonsoS.CherryE. M.SundnesJ. (2015). Simulations of Heart Function.Biomed Res. Int.2015:626378. 10.1155/2015/626378

  • 79

    DuguéG. P.BrunelN.HakimV.SchwartzE.ChatM.LévesqueM.et al (2009). Electrical coupling mediates tunable low-frequency oscillations and resonance in the cerebellar Golgi cell network.Neuron61126139. 10.1016/j.neuron.2008.11.028

  • 80

    Dura-BernalS.SuterB. A.GleesonP.CantarelliM.QuintanaA.RodriguezF.et al (2019). NetPyNE, a tool for data-driven multiscale modeling of brain circuits.eLife8:e44494. 10.7554/eLife.44494

  • 81

    EbnerT. J.ChenG. (1995). Use of voltage-sensitive dyes and optical recordings in the central nervous system.Prog. Neurobiol.46463506. 10.1016/0301-0082(95)00010-S

  • 82

    EcclesJ. C.LlinásR.SasakiK. (1966). Intracellularly recorded responses of the cerebellar Purkinje cells.Exp. Brain Res.1161183.

  • 83

    EngbersJ. D. T.FernandezF. R.TurnerR. W. (2013). Bistability in purkinje neurons: Ups and downs in cerebellar research.Neural Netw.471831. 10.1016/j.neunet.2012.09.006

  • 84

    EpplerJ. M. (2008). PyNEST: A convenient interface to the NEST simulator.Front. Neuroinform.2:12. 10.3389/neuro.11.012.2008

  • 85

    EshraA.HirrlingerP.HallermannS. (2019). Enriched environment shortens the duration of action potentials in cerebellar granule cells.Front. Cell. Neurosci.13:289. 10.3389/fncel.2019.00289

  • 86

    EvansW. H.MartinP. E. M. (2002). Gap junctions: Structure and function (Review).Mol. Membr. Biol.19121136. 10.1080/09687680210139839

  • 87

    FarrantM.FeldmeyerD.TakahashiT.Cull-candyS. G. (1994). NMDA-receptor channel diversity in the developing cerebellum.Nature368335339.

  • 88

    FarrantM.NusserZ. (2005). Variations on an inhibitory theme: Phasic and tonic activation of GABA(A) receptors.Nat. Rev. Neurosci.6215229. 10.1038/nrn1625

  • 89

    FlorimbiG.TortiE.MasoliS.D’AngeloE.DaneseG.LeporatiF. (2016). The Human brain project: Parallel technologies for biologically accurate simulation of granule cells.Microprocess. Microsyst.47111. 10.1016/j.micpro.2016.05.015

  • 90

    FlorimbiG.TortiE.MasoliS.D’AngeloE.DaneseG.LeporatiF. (2018). Exploiting multi-core and many-core architectures for efficient simulation of biologically realistic models of Golgi cells.J. Parallel Distrib. Comput.1264866. 10.1016/j.jpdc.2018.12.004

  • 91

    FlorimbiG.TortiE.MasoliS.D’AngeloE.LeporatiF. (2021). Granular layer simulator: design and Multi-GPU simulation of the cerebellar granular layer.Front. Comput. Neurosci.15:630795. 10.3389/fncom.2021.630795

  • 92

    FortiL.CesanaE.MapelliJ.D’AngeloE. (2006). Ionic mechanisms of autorhythmic firing in rat cerebellar Golgi cells.J. Physiol.574711729. 10.1113/jphysiol.2006.110858

  • 93

    FranksK. M.BartolT. M.SejnowskiT. J. (2002). A monte carlo model reveals independent signaling at central glutamatergic synapses.Biophys. J.8323332348. 10.1016/S0006-3495(02)75248-X

  • 94

    FunahashiA.MatsuokaY.JourakuA.MorohashiM.KikuchiN.KitanoH. (2008). CellDesigner 3.5: A versatile modeling tool for biochemical networks.Proc. IEEE9612541265. 10.1109/JPROC.2008.925458

  • 95

    GabbianiF.MidtgaardJ.KnöpfelT. (1994). Synaptic integration in a model of cerebellar granule cells.J. Neurophysiol.729991009. 10.1152/jn.1994.72.2.999

  • 96

    GaglianoG.MonteverdiA.CasaliS.LaforenzaU.Gandini Wheeler-KingshottC. A. M.D’AngeloE.et al (2022). Non-linear frequency dependence of neurovascular coupling in the cerebellar cortex implies vasodilation-vasoconstriction competition.Cells11:1047. 10.3390/cells11061047

  • 97

    GallianoE.MazzarelloP.D’AngeloE. (2010). Discovery and rediscoveries of Golgi cells.J. Physiol.58836393655. 10.1113/jphysiol.2010.189605

  • 98

    GandolfiD.PozziP.TognolinaM.ChiricoG.MapelliJ.D’AngeloE. (2014). The spatiotemporal organization of cerebellar network activity resolved by two-photon imaging of multiple single neurons.Front. Cell. Neurosci.8:92. 10.3389/fncel.2014.00092

  • 99

    GarridoJ. A.LuqueN. R.ToluS.D’AngeloE. (2016). Oscillation-Driven Spike-Timing Dependent Plasticity Allows Multiple Overlapping Pattern Recognition in Inhibitory Interneuron Networks.Int. J. Neural Syst.26:1650020. 10.1142/S0129065716500209

  • 100

    GeurtsF. J.De SchutterE.DieudonnéS. (2003). Unraveling the cerebellar cortex: Cytology and cellular physiology of large-sized interneurons in the granular layer.Cerebellum2290299. 10.1080/14734220310011948

  • 101

    GielenM. C.LumbM. J.SmartT. G. (2012). Benzodiazepines modulate GABA A receptors by regulating the preactivation step after GABA binding.J. Neurosci.3257075715. 10.1523/JNEUROSCI.5663-11.2012

  • 102

    GilbertT. L. (2018). The allen brain atlas as a resource for teaching undergraduate neuroscience.J. Undergrad. Neurosci. Educ.16:A261A267.

  • 103

    GilmerJ. I.PersonA. L. (2017). Morphological constraints on cerebellar granule cell combinatorial diversity.J. Neurosci.371215312166. 10.1523/JNEUROSCI.0588-17.2017

  • 104

    GolgiC. (1874). Sulla fina anatomia del cervelletto umano.Arch. Ital. Mal. Nerv.199111.

  • 105

    GoodenoughD. A.PaulD. L. (2009). Gap junctions.Cold Spring Harb. Perspect. Biol.1:a002576. 10.1101/cshperspect.a002576

  • 106

    GregerI. H.ZiffE. B.PennA. C. (2007). Molecular determinants of AMPA receptor subunit assembly.Trends Neurosci.30407416. 10.1016/j.tins.2007.06.005

  • 107

    GuoC.HusonV.MacoskoE. Z.RegehrW. G. (2021). Graded heterogeneity of metabotropic signaling underlies a continuum of cell-intrinsic temporal responses in unipolar brush cells.Nat. Commun.12213. 10.1038/s41467-021-22893-8

  • 108

    GuoS.ZhaoX.JiangS.DingL.PengH. (2022). Image enhancement to leverage the 3D morphological reconstruction of single-cell neurons.Bioinformatics38503512. 10.1093/bioinformatics/btab638

  • 109

    HaasJ. S.LandismanC. E. (2012). Bursts modify electrical synaptic strength.Brain Res.1487140149. 10.1016/j.brainres.2012.05.061

  • 110

    HallermannS.SilverR. A. (2013). Sustaining rapid vesicular release at active zones: Potential roles for vesicle tethering.Trends Neurosci.36185194. 10.1016/j.tins.2012.10.001

  • 111

    HanK. S.GuoC.ChenC. H.WitterL.OsornoT.RegehrW. G. (2018). Ephaptic Coupling Promotes Synchronous Firing of Cerebellar Purkinje Cells.Neuron100564578.e3e. 10.1016/j.neuron.2018.09.018

  • 112

    HanselC.LindenD. J.D’AngeloE. (2001). Beyond parallel fiber LTD: The diversity of synaptic and non-synaptic plasticity in the cerebellum.Nat. Neurosci.4467475. 10.1038/87419

  • 113

    HepburnI.ChenW.WilsS.De SchutterE. (2012). STEPS: Efficient simulation of stochastic reaction–diffusion models in realistic morphologies.BMC Syst. Biol.6:36. 10.1186/1752-0509-6-36

  • 114

    HickmanD. L.JohnsonJ.VemulapalliT. H.CrislerJ. R.ShepherdR. (2017). “Commonly used animal models,” in Principles of Animal Research, (Amsterdam: Elsevier), 117175. 10.1016/B978-0-12-802151-4.00007-4

  • 115

    HilleB. (2001). Ionic Channels of Excitable Membranes.Sunderland, MA: Sinauer.

  • 116

    HinesM.CarnevaleN. T. (2001). Neuron: A tool for neuroscientists.Neuroscientist7123135. 10.1177/107385840100700207

  • 117

    HinesM. L.CarnevaleN. T. (1997). The NEURON simulation environment.Neural Comput.911791209. 10.1162/neco.1997.9.6.1179

  • 118

    HinesM. L.CarnevaleN. T. (2008). Tranlating network models to parallel hardware in Neuron.J. Neurosci. Methods169425455.

  • 119

    HinesM. L.DavisonA. P.MullerE. (2009). NEURON and python.Front. Neuroinform.3:1. 10.3389/neuro.11.001.2009

  • 120

    HodgkinA. L.HuxleyA. F. (1952b). Currents carried by sodium and potassium ions through the membrane of the giant axon of Loligo.J. Physiol.116449472.

  • 121

    HodgkinA. L.HuxleyA. F. (1952a). A quantitative description of membrane current and its application to conduction and excitation in nerve. 1952.Bull. Math. Biol.522571.

  • 122

    HoehneA.McFaddenM. H.DigregorioD. A. (2020). Feed-forward recruitment of electrical synapses enhances synchronous spiking in the mouse cerebellar cortex.eLife9:e57344. 10.7554/ELIFE.57344

  • 123

    HullC.RegehrW. G. (2012). Identification of an inhibitory circuit that regulates cerebellar Golgi cell activity.Neuron73149158. 10.1016/j.neuron.2011.10.030

  • 124

    IacobucciG. J.PopescuG. K. (2018). Kinetic models for activation and modulation of NMDA receptor subtypes.Curr. Opin. Physiol.2114122. 10.1016/j.cophys.2018.02.002

  • 125

    ImamuraY.InokawaH.ItoA.KadotaniH.ToyamaK.NodaM.et al (2000). Roles of GABAergic inhibition and NMDA receptor subunits in the spatio-temporal integration in the cerebellar cortex of mice.Neurosci. Res.38289301. 10.1016/s0168-0102(00)00173-5

  • 126

    ItoM. (1984). The Cerebellum and Neural Control.New York, NY: Raven Press.

  • 127

    IzhikevichE. M. (2003). Simple model of spiking neurons.IEEE Trans. Neural Netw.1415691572. 10.1109/TNN.2003.820440

  • 128

    JansenB. H.RitV. G. (1995). Electroencephalogram and visual evoked potential generation in a mathematical model of coupled cortical columns.Biol. Cybern.73357366. 10.1007/BF00199471

  • 129

    JonesM. V.WestbrookG. L. (1995). Desensitized states prolong GABAA channel responses to brief agonist pulses.Neuron15181191. 10.1016/0896-6273(95)90075-6

  • 130

    JörntellH.EkerotC. F. (2006). Properties of somatosensory synaptic integration in cerebellar granule cells in vivo.J. Neurosci.261178611797. 10.1523/JNEUROSCI.2939-06.2006

  • 131

    KalmbachB. E.VoicuH.OhyamaT.MaukM. D. (2011). A subtraction mechanism of temporal coding in cerebellar cortex.J. Neurosci.3120252034. 10.1523/JNEUROSCI.4212-10.2011

  • 132

    KanichayR. T.SilverR. A. (2008). Synaptic and cellular properties of the feedforward inhibitory circuit within the input layer of the cerebellar cortex.J. Neurosci.2889558967. 10.1523/JNEUROSCI.5469-07.2008

  • 133

    KaseM.MillerD. C.NodaH. (1980). Discharges of Purkinje cells and mossy fibres in the cerebellar vermis of the monkey during saccadic eye movements and fixation.J. Physiol.300539555. 10.1113/jphysiol.1980.sp013178

  • 134

    KimJ.LeeS.TsudaS.ZhangX.AsricanB.GlossB.et al (2014). Optogenetic mapping of cerebellar inhibitory circuitry reveals spatially biased coordination of interneurons via electrical synapses.Cell Rep.716011613. 10.1016/j.celrep.2014.04.047

  • 135

    KimY. S.KangE.MakinoY.ParkS.ShinJ. H.SongH.et al (2013). Characterizing the conductance underlying depolarization-induced slow current in cerebellar Purkinje cells.J. Neurophysiol.10911741181. 10.1152/jn.01168.2011

  • 136

    KimY. S.ShinJ. H.HallF. S.LindenD. J. (2009). Dopamine signaling is required for depolarization-induced slow current in cerebellar Purkinje cells.J. Neurosci.2985308538. 10.1523/JNEUROSCI.0468-09.2009

  • 137

    KisielM.JatczakM.BrodzkiM.MozrzymasJ. W. (2018). Spontaneous activity, singly bound states and the impact of alpha1Phe64 mutation on GABAAR gating in the novel kinetic model based on the single-channel recordings.Neuropharmacology131453474. 10.1016/j.neuropharm.2017.11.030

  • 138

    KitaK.AlbergariaC.MachadoA. S.CareyM. R.MüllerM.DelvendahlI. (2021). Glua4 facilitates cerebellar expansion coding and enables associative memory formation.eLife10:e65152. 10.7554/eLife.65152

  • 139

    KleinleJ.VogtK.LüscherH. R.MüllerL.SennW.WylerK.et al (1996). Transmitter concentration profiles in the synaptic cleft: An analytical model of release and diffusion.Biophys. J.7124132426. 10.1016/S0006-3495(96)79435-3

  • 140

    KnöpfelT.GrandesP. (2002). Metabotropic glutamate receptors in the cerebellum with a focus on their function in Purkinje cells.Cerebellum11926. 10.1080/147342202753203050

  • 141

    KochC. (1984). Cable theory in neurons with active, linearized membranes.Biol. Cybern.501533. 10.1007/BF00317936

  • 142

    KochC. (1998). Biophysics of Computation: Information Processing in Single Neurons.New York, NY: Oxford University Press.

  • 143

    KoikeM.TsukadaS.TsuzukiK.KijimaH.OzawaS. (2000). Regulation of kinetic properties of GluR2 AMPA receptor channels by alternative splicing.J. Neurosci.2021662174.

  • 144

    KozarevaV.MartinC.OsornoT.RudolphS.GuoC.VanderburgC.et al (2021). A transcriptomic atlas of mouse cerebellar cortex comprehensively defines cell types.Nature598214219. 10.1038/s41586-021-03220-z

  • 145

    KoziolL. F.BuddingD.AndreasenN.D’ArrigoS.BulgheroniS.ImamizuH.et al (2014). Consensus paper: The cerebellum’s role in movement and cognition.Cerebellum13151177. 10.1007/s12311-013-0511-x

  • 146

    KulikÁNakadateK.NyíriG.NotomiT.MalitschekB.BettlerB.et al (2002). Distinct localization of GABAB receptors relative to synaptic sites in the rat cerebellum and ventrobasal thalamus.Eur. J. Neurosci.15291307. 10.1046/j.0953-816x.2001.01855.x

  • 147

    KumarN. M.GilulaN. B. (1996). The gap junction communication channel.Cell84381388. 10.1016/s0092-8674(00)81282-9

  • 148

    KumbharP.HinesM.FouriauxJ.OvcharenkoA.KingJ.DelalondreF.et al (2019). CoreNEURON: An optimized compute engine for the NEURON simulator.Front. Neuroinform.13:63. 10.3389/fninf.2019.00063

  • 149

    LaubeB.KuhseJ.BetzH. (1998). Evidence for a tetrameric structure of recombinant NMDA receptors.J. Neurosci.1829542961.

  • 150

    LaubeG.RöperJ.PittJ. C.SewingS.KistnerU.GarnerC. C.et al (1996). Ultrastructural localization of Shaker-related potassium channel subunits and synapse-associated protein 90 to septate-like junctions in rat cerebellar Pinceaux.Brain Res. Mol. Brain Res.425161. 10.1016/s0169-328x(96)00120-9

  • 151

    LeDouxM. S.HurstD. C.LordenJ. F. (1998). Single-unit activity of cerebellar nuclear cells in the awake genetically dystonic rat.Neuroscience86533545. 10.1016/S0306-4522(98)00007-4

  • 152

    LemaG. M. C.AuerbachA. (2006). Modes and models of GABAA receptor gating.J. Physiol.572183200. 10.1113/jphysiol.2005.099093

  • 153

    LennonW.Hecht-NielsenR.YamazakiT. (2014). A spiking network model of cerebellar Purkinje cells and molecular layer interneurons exhibiting irregular firing.Front. Comput. Neurosci.8:157. 10.3389/fncom.2014.00157

  • 154

    LennonW.YamazakiT.Hecht-NielsenR. (2015). A Model of In vitro Plasticity at the parallel fiber—molecular layer interneuron synapses.Front. Comput. Neurosci.9:150. 10.3389/fncom.2015.00150

  • 155

    LocatelliF.BottàL.PrestoriF.MasettoS.D’AngeloE. (2013). Late-onset bursts evoked by mossy fibre bundle stimulation in unipolar brush cells: Evidence for the involvement of H- and TRP-currents.J. Physiol.591899918. 10.1113/jphysiol.2012.242180

  • 156

    LoewensteinY.MahonS.ChaddertonP.KitamuraK.SompolinskyH.YaromY.et al (2005). Bistability of cerebellar Purkinje cells modulated by sensory stimulation.Nat. Neurosci.8202211. 10.1038/nn1393

  • 157

    LuH.-W.BalmerT. S.RomeroG. E.TrussellL. O. (2017). Slow AMPAR synaptic transmission is determined by stargazin and glutamate transporters.Neuron967380.e4. 10.1016/j.neuron.2017.08.043

  • 158

    LujánR.AguadoC.CiruelaF.CózarJ.KleindienstD.de la OssaL.et al (2018). Differential association of GABAB receptors with their effector ion channels in Purkinje cells.Brain Struct. Funct.22315651587. 10.1007/s00429-017-1568-y

  • 159

    MaconochieD. J.ZempelJ. M.SteinbachJ. H. (1994). How quickly can GABAA receptors open?Neuron126171. 10.1016/0896-6273(94)90152-X

  • 160

    MaexR.De SchutterE. (1998). Synchronization of golgi and granule cell firing in a detailed network model of the cerebellar granule cell layer.J. Neurophysiol.8025212537. 10.1152/jn.1998.80.5.2521

  • 161

    MapelliJ.D’AngeloE. (2007). The spatial organization of long-term synaptic plasticity at the input stage of cerebellum.J. Neurosci.2712851296. 10.1523/JNEUROSCI.4873-06.2007

  • 162

    MapelliL.SolinasS.D’AngeloE. (2014). Integration and regulation of glomerular inhibition in the cerebellar granular layer circuit.Front. Cell. Neurosci.8:55. 10.3389/fncel.2014.00055

  • 163

    MarcaggiP.MutohH.DimitrovD.BeatoM.KnöpfelT. (2009). Optical measurement of mGluR1 conformational changes reveals fast activation, slow deactivation, and sensitization.Proc. Natl. Acad. Sci. U.S.A.1061138811393. 10.1073/pnas.0901290106

  • 164

    MarkramH.MullerE.RamaswamyS.ReimannM. W.AbdellahM.SanchezC. A.et al (2015). Reconstruction and simulation of neocortical microcircuitry.Cell163456492. 10.1016/j.cell.2015.09.029

  • 165

    MarrD. (1969). A theory of cerebellar cortex.J. Physiol.202437470. 10.2307/1776957

  • 166

    MasoliS.D’AngeloE. (2017). Synaptic activation of a detailed purkinje cell model predicts voltage-dependent control of burst-pause responses in active dendrites.Front. Cell. Neurosci.11:278. 10.3389/fncel.2017.00278

  • 167

    MasoliS.TognolinaM.LaforenzaU.MocciaF.D’AngeloE. (2020b). Parameter tuning differentiates granule cell subtypes enriching transmission properties at the cerebellum input stage.Commun. Biol.3:222. 10.1038/s42003-020-0953-x

  • 168

    MasoliS.OttavianiA.CasaliS.D’AngeloE. (2020a). Cerebellar Golgi cell models predict dendritic processing and mechanisms of synaptic plasticity.PLoS Comput. Biol.16:e1007937.

  • 169

    MasoliS.SolinasS.D’AngeloE. (2015). Action potential processing in a detailed Purkinje cell model reveals a critical role for axonal compartmentalization.Front. Cell. Neurosci.9:47. 10.3389/fncel.2015.00047

  • 170

    Maxwell CowanW.SüdhofT. C.StevensC. F. (2000). Synapses.Baltimore, MA: Johns Hopkins University Press.

  • 171

    MedinaJ. F.MaukM. D. (2000). Computer simulation of cerebellar information processing.Nat. Neurosci.312051211. 10.1038/81486

  • 172

    MenigozA.AhmedT.SabanovV.PhilippaertK.PintoS.KerselaersS.et al (2016). TRPM4-dependent post-synaptic depolarization is essential for the induction of NMDA receptor-dependent LTP in CA1 hippocampal neurons.Pflugers Arch. Eur. J. Physiol.468593607. 10.1007/s00424-015-1764-7

  • 173

    MiccoliB.LopezC. M.GoikoetxeaE.PutzeysJ.SekeriM.KrylychkinaO.et al (2019). High-density electrical recording and impedance imaging with a multi-modal CMOS multi-electrode array chip.Front. Neurosci.13:641. 10.3389/fnins.2019.00641

  • 174

    MitchellS. J.SilverR. A. (2003). Shunting inhibition modulates neuronal gain during synaptic excitation.Neuron38433445. 10.1016/s0896-6273(03)00200-9

  • 175

    MittmannW.KochU.HäusserM. (2005). Feed-forward inhibition shapes the spike output of cerebellar Purkinje cells.J. Physiol.563369378. 10.1113/jphysiol.2004.075028

  • 176

    MiyazakiT.YamasakiM.TanakaK. F.WatanabeM. (2021). Compartmentalized input–output organization of lugaro cells in the cerebellar cortex.Neuroscience46289105. 10.1016/j.neuroscience.2020.05.026

  • 177

    MonyerH.BurnashevN.LaurieD. J.SakmannB.SeeburgP. H. (1994). Developmental and regional expression in the rat brain and functional properties of four NMDA receptors.Neuron12529540. 10.1016/0896-6273(94)90210-0

  • 178

    MugnainiE.DiñoM. R.JaarsmaD. (1997). The unipolar brush cells of the mammalian cerebellum and cochlear nucleus: Cytology and microcircuitry.Prog. Brain Res.114131150. 10.1016/s0079-6123(08)63362-2

  • 179

    MugnainiE.SekerkováG.MartinaM. (2011). The unipolar brush cell: A remarkable neuron finally receiving deserved attention.Brain Res. Rev.66220245. 10.1016/j.brainresrev.2010.10.001

  • 180

    NedelescuH.AbdelhackM. (2013). Comparative morphology of dendritic arbors in populations of purkinje cells in mouse sulcus and apex.Neural Plast.2013:948587. 10.1155/2013/948587

  • 181

    NedelescuH.AbdelhackM.PritchardA. T. (2018). Regional differences in Purkinje cell morphology in the cerebellar vermis of male mice.J. Neurosci. Res.9614761489. 10.1002/jnr.24206

  • 182

    NewtonA. J. H.McDougalR. A.HinesM. L.LyttonW. W. (2018). Using NEURON for reaction-diffusion modeling of extracellular dynamics.Front. Neuroinform.12:41. 10.3389/fninf.2018.00041

  • 183

    NielsenM. S.AxelsenL. N.SorgenP. L.VermaV.DelmarM.Holstein-RathlouN.-H. (2012). Gap junctions.Compr. Physiol.219812035. 10.1002/cphy.c110051

  • 184

    NieusT.SolaE.MapelliJ.SaftenkuE.RossiP.D’AngeloE. (2006). LTP regulates burst initiation and frequency at mossy fiber-granule cell synapses of rat cerebellum: Experimental observations and theoretical predictions.J. Neurophysiol.95686699. 10.1152/jn.00696.2005

  • 185

    NieusT. R.MapelliL.D’AngeloE. (2014). Regulation of output spike patterns by phasic inhibition in cerebellar granule cells.Front. Cell. Neurosci.8:246. 10.3389/fncel.2014.00246

  • 186

    NowakW. (2016). “Applications of Computational Methods to Simulations of Protein Dynamics,” in Handbook of Computational Chemistry, ed.LeszczynskiJ. (Dordrecht: Springer Netherlands), 143. 10.1007/978-94-007-6169-8_31-2

  • 187

    NowkeC.ZielaskoD.WeyersB.PeyserA.HentschelB.KuhlenT. W. (2015). Integrating visualizations into modeling NEST simulations.Front. Neuroinform.9:29. 10.3389/fninf.2015.00029

  • 188

    OlsenR. W.SieghartW. (2008). International union of pharmacology. LXX. Subtypes of gamma-aminobutyric acid(A) receptors: Classification on the basis of subunit composition, pharmacology, and function. Update.Pharmacol. Rev.60243260. 10.1124/pr.108.00505

  • 189

    OsornoT.RudolphS.NguyenT.KozarevaV.NadafN. M.NortonA.et al (2022). Candelabrum cells are ubiquitous cerebellar cortex interneurons with specialized circuit properties.Nat. Neurosci.25702713. 10.1038/s41593-022-01057-x

  • 190

    PalayS. L.Chan-PalayV. (1974). Cerebellar Cortex Cytology and Organization.Berlin: Springer.

  • 191

    PalesiF.FerranteM.GaviraghiM.MisitiA.SaviniG.LascialfariA.et al (2021). Motor and higher-order functions topography of the human dentate nuclei identified with tractography and clustering methods.Hum. Brain Mapp.4243484361. 10.1002/hbm.25551

  • 192

    PernelleG.NicolaW.ClopathC. (2018). Gap junction plasticity as a mechanism to regulate network-wide oscillations.PLoS Comput. Biol.14:e1006025. 10.1371/journal.pcbi.1006025

  • 193

    PhamT.HaasJ. S. (2018). Electrical synapses between inhibitory neurons shape the responses of principal neurons to transient inputs in the thalamus: A modeling study.Sci. Rep.8:7763. 10.1038/s41598-018-25956-x

  • 194

    PhamT.HaasJ. S. (2019). Electrical synapses regulate both subthreshold integration and population activity of principal cells in response to transient inputs within canonical feedforward circuits.PLoS Comput. Biol.15:e1006440. 10.1371/journal.pcbi.1006440

  • 195

    PowerE. M.EmpsonR. M. (2014). Functional contributions of glutamate transporters at the parallel fibre to Purkinje neuron synapse-relevance for the progression of cerebellar ataxia.Cerebellum and Ataxias1:3. 10.1186/2053-8871-1-3

  • 196

    PrestoriF.MocciaF.D’AngeloE. (2019). Disrupted calcium signaling in animal models of human spinocerebellar ataxia (SCA).Int. J. Mol. Sci.21:E216. 10.3390/ijms21010216

  • 197

    RanczE. A.IshikawaT.DuguidI.ChaddertonP.MahonS.HäusserM. (2007). High-fidelity transmission of sensory information by single cerebellar mossy fibre boutons.Nature45012451248. 10.1038/nature05995

  • 198

    RasminskyM. (1980). Ephaptic transmission between single nerve fibres in the spinal nerve roots of dystrophic mice.J. Physiol.305151169. 10.1113/jphysiol.1980.sp013356

  • 199

    RinaldoL.HanselC. (2010). Ataxias and cerebellar dysfunction: Involvement of synaptic plasticity deficits?Funct. Neurol.25135139.

  • 200

    RizzaM. F.LocatelliF.MasoliS.Sánchez-PonceD.MuñozA.PrestoriF.et al (2021). Stellate cell computational modeling predicts signal filtering in the molecular layer circuit of cerebellum.Sci. Rep.11:3873. 10.1038/s41598-021-83209-w

  • 201

    RomaniA.SchürmannF.MarkramH.MiglioreM. (2022). “Reconstruction of the Hippocampus,” in Computational Modelling of the Brain. Advances in Experimental Medicine and Biology, edsGiuglianoM.NegrelloM.LinaroD. (Berlin: Springer), 261283. 10.1007/978-3-030-89439-9_11

  • 202

    RössertC.DeanP.PorrillJ. (2015). At the edge of chaos: How cerebellar granular layer network dynamics can provide the basis for temporal filters.PLoS Comput. Biol.11:e1004515. 10.1371/journal.pcbi.1004515

  • 203

    RothmanJ. S.CathalaL.SteuberV.SilverR. A. (2009). Synaptic depression enables neuronal gain control.Nature45710151018. 10.1038/nature07604

  • 204

    RusakovD. A. (2001). The role of perisynaptic glial sheaths in glutamate spillover and extracellular Ca(2+) depletion.Biophys. J.8119471959.

  • 205

    RusakovD. A.KullmannD. M. (1998). Extrasynaptic glutamate diffusion in the hippocampus: Ultrastructural constraints, uptake, and receptor activation.J. Neurosci.1831583170. 10.1523/jneurosci.18-09-03158.1998

  • 206

    SaabA. S.NeumeyerA.JahnH. M.CupidoA.ŠimekA. A. M.BoeleH.-J.et al (2012). Bergmann glial AMPA receptors are required for fine motor coordination.Science337749753. 10.1126/science.1221140

  • 207

    SaftenkuE. E. (2005). Modeling of slow glutamate diffusion and AMPA receptor activation in the cerebellar glomerulus.J. Theor. Biol.234363382. 10.1016/j.jtbi.2004.11.036

  • 208

    SandersH.BerendsM.MajorG.GoldmanM. S.LismanJ. E. (2013). NMDA and GABAB (KIR) conductances: The “perfect couple” for bistability.J. Neurosci.33424429. 10.1523/JNEUROSCI.1854-12.2013

  • 209

    SantamariaF.TrippP. G.BowerJ. M. (2007). Feedforward inhibition controls the spread of granule cell-induced Purkinje cell activity in the cerebellar cortex.J. Neurophysiol.97248263. 10.1152/jn.01098.2005

  • 210

    Sanz LeonP.KnockS. A.WoodmanM. M.DomideL.MersmannJ.McIntoshA. R.et al (2013). The virtual brain: A simulator of primate brain network dynamics.Front. Neuroinform.7:10. 10.3389/fninf.2013.00010

  • 211

    SavtchenkoL. P.RusakovD. A. (2007). The optimal height of the synaptic cleft.Proc. Natl. Acad. Sci. U.S.A.10418231828. 10.1073/pnas.0606636104

  • 212

    SchirnerM.McIntoshA. R.JirsaV.DecoG.RitterP. (2018). Inferring multi-scale neural mechanisms with brain network modelling.eLife7:e28927. 10.7554/eLife.28927

  • 213

    SchmahmannJ. D. (2010). The role of the cerebellum in cognition and emotion: Personal reflections since 1982 on the dysmetria of thought hypothesis, and its historical evolution from theory to therapy.Neuropsychol. Rev.20236260. 10.1007/s11065-010-9142-x

  • 214

    SchwartzE. J.RothmanJ. S.DuguéG. P.DianaM. (2012). NMDA receptors with incomplete Mg 2 + block enable low- frequency transmission through the cerebellar cortex.J. Neurosci.3268786893. 10.1523/JNEUROSCI.5736-11.2012.NMDA

  • 215

    SegevI. (1998). Methods in neuronal modeling: From ions to networks (second edition).Comput. Math. Appl.10.1016/s0898-1221(98)91179-8

  • 216

    SevetsonJ.HaasJ. S. (2015). Asymmetry and modulation of spike timing in electrically coupled neurons.J. Neurophysiol.11317431751. 10.1152/jn.00843.2014

  • 217

    ShannonP. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks.Genome Res.1324982504. 10.1101/gr.1239303

  • 218

    SilverR. A.TraynelisS. F.Cull-CandyS. G. (1992). Rapid-time-course miniature and evoked excitatory currents at cerebellar synapses in situ.Nature355163166. 10.1038/355163a0

  • 219

    Simões de SouzaF. M.De SchutterE. (2011). Robustness effect of gap junctions between Golgi cells on cerebellar cortex oscillations.Neural Syst. Circuits1:7. 10.1186/2042-1001-1-7

  • 220

    SnipasM.RimkuteL.KraujalisT.MaciunasK. (2017). Functional asymmetry and plasticity of electrical synapses interconnecting neurons through a 36-state model of gap junction channel gating.PLoS Comput. Biol.13:e1005464. 10.1371/journal.pcbi.1005464

  • 221

    SodaT.MapelliL.LocatelliF.BottaL.GoldfarbM.PrestoriF.et al (2019). Erratum: Hyperexcitability and hyperplasticity disrupt cerebellar signal transfer in the ib2 ko mouse model of autism.J. Neurosci.39:7029. 10.1523/JNEUROSCI.1332-19.2019

  • 222

    SolaE.PrestoriF.RossiP.TagliettiV.D’AngeloE. (2004). Increased neurotransmitter release during long-term potentiation at mossy fibre-granule cell synapses in rat cerebellum.J. Physiol.557843861. 10.1113/jphysiol.2003.060285

  • 223

    SolinasS.FortiL.CesanaE.MapelliJ.De SchutterE.D’AngeloE.et al (2007a). Fast-reset of pacemaking and theta-frequency resonance patterns in cerebellar golgi cells: Simulations of their impact in vivo.Front. Cell. Neurosci.1:4. 10.3389/neuro.03.004.2007

  • 224

    SolinasS.FortiL.CesanaE.MapelliJ.De SchutterE.D’AngeloE. (2007b). Computational reconstruction of pacemaking and intrinsic electroresponsiveness in cerebellar Golgi cells.Front. Cell. Neurosci.1:2. 10.3389/neuro.03.002.2007

  • 225

    SolinasS.NieusT.D’AngeloE. (2010). A realistic large-scale model of the cerebellum granular layer predicts circuit spatio-temporal filtering properties.Front. Cell. Neurosci.4:12. 10.3389/fncel.2010.00012

  • 226

    SoteloC. (2008). Viewing the cerebellum through the eyes of Ram??n Y Cajal.Cerebellum7517522. 10.1007/s12311-008-0078-0

  • 227

    SotomayorM. (2015). Computational exploration of single-protein mechanics by steered molecular dynamics.AIP Conf. Proc.1703:030001. 10.1063/1.4939316

  • 228

    SouthanA. P.RobertsonB. (2000). Electrophysiological characterization of voltage-gated K(+) currents in cerebellar basket and purkinje cells: Kv1 and Kv3 channel subfamilies are present in basket cell nerve terminals.J. Neurosci.20114122. 10.1523/JNEUROSCI.20-01-00114.2000

  • 229

    SpieglerA.KnöscheT. R.SchwabK.HaueisenJ.AtayF. M. (2011). Modeling brain resonance phenomena using a neural mass model.PLoS Comput. Biol.7:e1002298. 10.1371/journal.pcbi.1002298

  • 230

    SteinR. B. (1967). Some models of neuronal variability.Biophys. J.73768. 10.1016/S0006-3495(67)86574-3

  • 231

    SteuberV.MittmannW.HoebeekF. E.SilverR. A.De ZeeuwC. I.HäusserM.et al (2007). Cerebellar LTD and pattern recognition by Purkinje cells.Neuron54121136. 10.1016/j.neuron.2007.03.015

  • 232

    StimbergM.BretteR.GoodmanD. F. M. (2019). Brian 2, an intuitive and efficient neural simulator.eLife8:e47314. 10.7554/eLife.47314

  • 233

    StraubI.WitterL.EshraA.HoidisM.ByczkowiczN.MaasS.et al (2020). Gradients in the mammalian cerebellar cortex enable Fourier-like transformation and improve storing capacity.eLife9:e51771. 10.7554/eLife.51771

  • 234

    SubramaniyamS.SolinasS.PerinP.LocatelliF.MasettoS.D’AngeloE. (2014). Computational modeling predicts the ionic mechanism of late-onset responses in unipolar brush cells.Front. Cell. Neurosci.8:237. 10.3389/fncel.2014.00237

  • 235

    SudhakarS. K.HongS.RaikovI.PublioR.LangC.CloseT.et al (2017). Spatiotemporal network coding of physiological mossy fiber inputs by the cerebellar granular layer.PLoS Comput. Biol.13:e1005754. 10.1371/journal.pcbi.1005754

  • 236

    SüdhofT. C.MalenkaR. C. (2008). Understanding synapses: Past, present, and future.Neuron60469476. 10.1016/j.neuron.2008.10.011

  • 237

    SultanF.BowerJ. M. (1998). Quantitative Golgi study of the rat cerebellar molecular layer interneurons using principal component analysis.J. Comp. Neurol.393353373.

  • 238

    SzczotM.KisielM.CzyzewskaM. M.MozrzymasJ. W. (2014). α1F64 residue at GABAA receptor binding site is involved in gating by influencing the receptor flipping transitions.J. Neurosci.3431933209. 10.1523/JNEUROSCI.2533-13.2014

  • 239

    TangT.XiaoJ.SuhC. Y.BurroughsA.CerminaraN. L.JiaL.et al (2017). Heterogeneity of Purkinje cell simple spike–complex spike interactions: Zebrin- and non-zebrin-related variations.J. Physiol.59553415357. 10.1113/JP274252

  • 240

    TaoL.NicholsonC. (2004). Maximum geometrical hindrance to diffusion in brain extracellular space surrounding uniformly spaced convex cells.J. Theor. Biol.2295968. 10.1016/j.jtbi.2004.03.003

  • 241

    TaylorC. P.DudekF. E. (1982). Synchronous neural afterdischarges in rat hippocampal slices without active chemical synapses.Science218810812. 10.1126/science.7134978

  • 242

    TognolinaM.MonteverdiA.D’AngeloE. (2022). Discovering microcircuit secrets with multi-spot imaging and electrophysiological recordings: The example of cerebellar network dynamics.Front. Cell. Neurosci.16:805670. 10.3389/fncel.2022.805670

  • 243

    Torben-NielsenB.SegevI.YaromY. (2012). The generation of phase differences and frequency changes in a network model of inferior olive subthreshold oscillations.PLoS Comput. Biol.8:e1002580. 10.1371/journal.pcbi.1002580

  • 244

    TortiE.MasoliS.FlorimbiG.D’AngeloE.TicliM.LeporatiF. (2019). “GPU Parallelization of Realistic Purkinje Cells with Complex Morphology,” in Proceedings - 27th Euromicro International Conference on Parallel, Distributed and Network-Based Processing, PDP 2019, (Pavia), 10.1109/EMPDP.2019.8671581

  • 245

    TsodyksM.PawelzikK.MarkramH. (1998). Neural networks with dynamic synapses.Neural Comput.10821835. 10.1162/089976698300017502

  • 246

    TuckwellH. C. (1988). Introduction to Theoretical Neurobiology.Cambridge: Cambridge University Press, 10.1017/CBO9780511623202

  • 247

    UusisaariM.KnöpfelT. (2011). Functional classification of neurons in the mouse lateral cerebellar nuclei.Cerebellum10637646. 10.1007/s12311-010-0240-3

  • 248

    UusisaariM.ObataK.KnöpfelT. (2007). Morphological and electrophysiological properties of GABAergic and non-GABAergic cells in the deep cerebellar nuclei.J. Neurophysiol.97901911. 10.1152/jn.00974.2006

  • 249

    van BeugenB. J.GaoZ.BoeleH.-J.HoebeekF.De ZeeuwC. I. (2013). High frequency burst firing of granule cells ensures transmission at the parallel fiber to purkinje cell synapse at the cost of temporal coding.Front. Neural Circuits7:95. 10.3389/fncir.2013.00095

  • 250

    Van GeitW.GevaertM.ChindemiG.RössertC.CourcolJ.-D.MullerE. B.et al (2016). BluePyOpt: Leveraging open source software and cloud infrastructure to optimise model parameters in neuroscience.Front. Neuroinform.10:17. 10.3389/fninf.2016.00017

  • 251

    VaughnM. J.HaasJ. S. (2022). On the diverse functions of electrical synapses.Front. Cell. Neurosci.16:910015. 10.3389/fncel.2022.910015

  • 252

    VervaekeK.LÖrinczA.GleesonP.FarinellaM.NusserZ.SilverR. A. (2010). Rapid desynchronization of an electrically coupled interneuron network with sparse excitatory synaptic input.Neuron67435451. 10.1016/j.neuron.2010.06.028

  • 253

    VoogdJ.GlicksteinM. (1998). The anatomy of the cerebellum.Trends Neurosci.21370375. 10.1016/s0166-2236(98)01318-6

  • 254

    VoŽehF. (2015). Jan evangelista purkyně and the cerebellum then and now.Physiol. Res.64:S567S584.

  • 255

    VyklickýL. (1993). Calcium-mediated modulation of N-methyl-D-aspartate (n.d.) responses in cultured rat hippocampal neurones.J. Physiol.470575600. 10.1113/jphysiol.1993.sp019876

  • 256

    WadicheJ. I.JahrC. E. (2001). Multivesicular release at climbing fiber-purkinje cell synapses evidence for the one-site, one-vesicle hypothesis comes from combined morphological and electrophysi-ological studies suggesting that the number of morpho.Neuron32301313.

  • 257

    WhiteJ. J.SillitoeR. V. (2017). Genetic silencing of olivocerebellar synapses causes dystonia-like behaviour in mice.Nat. Commun.8:14912. 10.1038/ncomms14912

  • 258

    WonH.MahW.KimE. (2013). Autism spectrum disorder causes, mechanisms, and treatments: Focus on neuronal synapses.Front. Mol. Neurosci.6:19. 10.3389/fnmol.2013.00019

  • 259

    WuQ.KapfhammerJ. P. (2022). The emerging key role of the mGluR1-PKC γ signaling pathway in the pathogenesis of spinocerebellar ataxias: A neurodevelopmental viewpoint.Int. J. Mol. Sci.23:9169. 10.3390/ijms23169169

  • 260

    YamadaW. M.ZuckerR. S. (1992). Time course of transmitter release calculated from simulations of a calcium diffusion model.Biophys. J.61671682. 10.1016/S0006-3495(92)81872-6

  • 261

    YamasakiM.AibaA.KanoM.WatanabeM. (2021). mGluR1 signaling in cerebellar Purkinje cells: Subcellular organization and involvement in cerebellar function and disease.Neuropharmacology194:108629. 10.1016/j.neuropharm.2021.108629

  • 262

    YangW.Carrillo-ReidL.BandoY.PeterkaD. S.YusteR. (2018). Simultaneous two-photon imaging and two-photon optogenetics of cortical circuits in three dimensions.eLife7:e32671. 10.7554/eLife.32671

  • 263

    YiheL.TimofeevaY. (2020). Exact solutions to cable equations in branching neurons with tapering dendrites.J. Math. Neurosci.10:1. 10.1186/s13408-020-0078-z

  • 264

    ZhangB.ChenL. Y.LiuX.MaxeinerS.LeeS. J.GokceO.et al (2015). Neuroligins sculpt cerebellar purkinje-cell circuits by differential control of distinct classes of synapses.Neuron87781796. 10.1016/j.neuron.2015.07.020

  • 265

    ZhangB.SudhofT. C. (2016). Neuroligins are selectively essential for NMDAR signaling in cerebellar stellate interneurons.J. Neurosci.3690709083. 10.1523/JNEUROSCI.1356-16.2016

  • 266

    ZhengK.ScimemiA.RusakovD. A. (2008). Receptor actions of synaptically released glutamate: The role of transporters on the scale from nanometers to microns.Biophys. J.9545844596. 10.1529/biophysj.108.129874

  • 267

    ZhouH.VogesK.LinZ.JuC.SchonewilleM. (2015). Differential purkinje cell simple spike activity and pausing behavior related to cerebellar modules.J. Neurophysiol.11325242536. 10.1152/jn.00925.2014

Summary

Keywords

cerebellum, synapses, receptors, computational model, purkinje cell, granule cell

Citation

Masoli S, Rizza MF, Tognolina M, Prestori F and D’Angelo E (2022) Computational models of neurotransmission at cerebellar synapses unveil the impact on network computation. Front. Comput. Neurosci. 16:1006989. doi: 10.3389/fncom.2022.1006989

Received

29 July 2022

Accepted

10 October 2022

Published

28 October 2022

Volume

16 - 2022

Edited by

Vito Di Maio, Institute of Applied Sciences and Intelligent Systems, Italy

Reviewed by

Igor Delvendahl, University of Zurich, Switzerland; De-Lai Qiu, Yanbian University, China

Updates

Copyright

*Correspondence: Francesca Prestori, Egidio D’Angelo,

These authors share last authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics