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

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 | τin | Gmax AMPA | Subtype | Gmax NMDA | Subtype | GABA-A | Subtype | REF | |
| Mf-GrC | 0.43 | 5 | 8 | 1 | 1,200 | GluR2 | 18,800 | NR2A-NR2C | − | − | a |
| 0.46 | No presynaptic | 1,200 | GluR2 | − | − | − | − | b | |||
| 0.416 | 10.8 | 35.1 | 3 | 700 | GluR2 | − | − | − | − | c | |
| 0.416 | 10.8 | 35.1 | 3 | − | - | 10,000 | NR2A-NR2C | − | − | c | |
| 0.42 | 8 | 5 | 1 | 1,200 | GluR2 | 18,800 | NR2A-NR2C | − | − | d | |
| GoC-GrC | 0.35 | 58 | 36 | 1 | 1,050 | GluR2 | − | − | − | α1/6 | e |
| Mf-GoC | 0.43 | 5 | 8 | 1 | 1,200 | GluR2 | 10,000 | NR2B | − | − | f |
| GrC-GoC pf | 0.4 | 54 | 35.1 | 1 | 1,200 | GluR2 | − | − | − | − | f |
| GrC-GoC aa | 0.4 | 54 | 35.1 | 1 | 1,200 | GluR2 | 10,000 | NR2B | − | − | f |
| GrC-GoC (pf-aa) | No presynaptic | 2,000 | GluR2 | g | |||||||
| GrC-PC pf | 0.13 | 54 | 35.1 | 1 | 2,800 | GluR2 | − | − | − | − | g |
| GrC-PC aa | 0.13 | 54 | 35.1 | 1 | 2,800 | GluR2 | − | − | − | − | g |
| GrC-SC | 0.15 | 10.8 | 35.1 | 1 | 2,300 | GluR2 | 10,000 | NR2B | − | − | h |
| GrC-BC | 0.15 | 10.8 | 35.1 | 1 | 2,300 | GluR2 | 10,000 | NR2B | − | − | i |
| SC-PC | 0.35 | 4 | 15 | 1 | − | − | − | − | 2,600 | α1 | h |
| BC-PC | 0.35 | 4 | 15 | 1 | − | − | − | − | 2,600 | α1 | i |
| SC-SC | 0.42 | 4 | 38.7 | 1 | − | − | − | − | 1,600 | α1 | h |
| BC-BC | 0.42 | 4 | 38.7 | 1 | − | − | − | − | 1,600 | α1 | i |
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

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.Nature569413–417. 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), 274–282. 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.264602–4612. 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.365079–5100. 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:E4798–E4807. 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.114608–623. 10.1152/jn.00628.2014
7
AnastassiouC. A.PerinR.MarkramH.KochC. (2011). Ephaptic coupling of cortical neurons.Nat. Neurosci.14217–223. 10.1038/nn.2727
8
AndreaeL. C.BurroneJ. (2014). The role of neuronal activity and transmitter release on synapse formation.Curr. Opin. Neurobiol.2747–52. 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.Cerebellum11681–693. 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.Synapse70269–276. 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.217969–7984. 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.33295–303. 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.1211381–1397. 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.Cerebellum14197–220. 10.1007/s12311-014-0627-7
18
BaumannS. W.BaurR.SigelE. (2001). Subunit arrangement of gamma-aminobutyric acid type A receptors.J. Biol. Chem.27636275–36280. 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.27746020–46025. 10.1074/jbc.M207663200
20
BenderK. J.TrussellL. O. (2012). The physiology of the axon initial segment.Annu. Rev. Neurosci.35249–265. 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.849–76. 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.428333–357. 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.17289–295. 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.1026168–178. 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.401103–125. 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–4862–73. 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.205899–5905. 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.41528–538. 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.1679–93. 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.3312430–12446. 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.18718–727. 10.1038/nn.3974.Synaptic
41
ChaddertonP.MargrieT. W.HäusserM. (2004). Integration of quanta in cerebellar granule cells during sensory processing.Nature428856–860. 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.13593–24. 10.1007/978-3-030-89439-9_1
44
ChowC. C.KopellN. (2000). Dynamics of spiking neurons with electrical coupling.Neural Comput.121643–1678. 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.Science2581498–1501. 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.22366–2382. 10.1038/nprot.2007.324
48
Cull-CandyS.BrickleyS.FarrantM. (2001). NMDA receptor subunits: Diversity, development and disease.Curr. Opin. Neurobiol.11327–335. 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.484397–413. 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.80493–503. 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.3230–40. 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.3512107–114. 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.21759–770. 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.81277–287. 10.1152/jn.1999.81.1.277
59
D’AngeloE.RossiP.GarthwaiteJ. (1990). Dual-component NMDA receptor currents at a single central synapse.Nature346467–470. 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.Neuroscience53121–130. 10.1016/0306-4522(93)90290-v
61
D’AngeloE. D.JirsaV. (2022). Neurosciences The quest for multiscale brain modeling.Trends Neurosci.xx1–14. 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.3117–32. 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.12481–93. 10.1016/s0079-6123(00)24009-0
66
De ZeeuwC. I.LisbergerS. G.RaymondJ. L. (2021). Diversity and dynamism in the cerebellum.Nat. Neurosci.24160–167. 10.1038/s41593-020-00754-9
67
DeanP.PorrillJ. (2011). Evaluating the adaptive-filter model of the cerebellum.J. Physiol.5893459–3470. 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.39722–737. 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.5875843–5857. 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.79999–1016. 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.1195–230. 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.252906–2916. 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.510845–866. 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.Neuron35521–533. 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.278344–8357. 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.101519–532. 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.Neuron61126–139. 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.46463–506. 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.1161–183.
83
EngbersJ. D. T.FernandezF. R.TurnerR. W. (2013). Bistability in purkinje neurons: Ups and downs in cerebellar research.Neural Netw.4718–31. 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.19121–136. 10.1080/09687680210139839
87
FarrantM.FeldmeyerD.TakahashiT.Cull-candyS. G. (1994). NMDA-receptor channel diversity in the developing cerebellum.Nature368335–339.
88
FarrantM.NusserZ. (2005). Variations on an inhibitory theme: Phasic and tonic activation of GABA(A) receptors.Nat. Rev. Neurosci.6215–229. 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.471–11. 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.12648–66. 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.574711–729. 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.832333–2348. 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. IEEE961254–1265. 10.1109/JPROC.2008.925458
95
GabbianiF.MidtgaardJ.KnöpfelT. (1994). Synaptic integration in a model of cerebellar granule cells.J. Neurophysiol.72999–1009. 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.5883639–3655. 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.Cerebellum2290–299. 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.325707–5715. 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:A261–A267.
103
GilmerJ. I.PersonA. L. (2017). Morphological constraints on cerebellar granule cell combinatorial diversity.J. Neurosci.3712153–12166. 10.1523/JNEUROSCI.0588-17.2017
104
GolgiC. (1874). Sulla fina anatomia del cervelletto umano.Arch. Ital. Mal. Nerv.199–111.
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.30407–416. 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.122–13. 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.Bioinformatics38503–512. 10.1093/bioinformatics/btab638
109
HaasJ. S.LandismanC. E. (2012). Bursts modify electrical synaptic strength.Brain Res.1487140–149. 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.36185–194. 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.Neuron100564–578.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.4467–475. 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), 117–175. 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.Neuroscientist7123–135. 10.1177/107385840100700207
117
HinesM. L.CarnevaleN. T. (1997). The NEURON simulation environment.Neural Comput.91179–1209. 10.1162/neco.1997.9.6.1179
118
HinesM. L.CarnevaleN. T. (2008). Tranlating network models to parallel hardware in Neuron.J. Neurosci. Methods169425–455.
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.116449–472.
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.5225–71.
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.Neuron73149–158. 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.2114–122. 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.38289–301. 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.141569–1572. 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.73357–366. 10.1007/BF00199471
129
JonesM. V.WestbrookG. L. (1995). Desensitized states prolong GABAA channel responses to brief agonist pulses.Neuron15181–191. 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.2611786–11797. 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.312025–2034. 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.288955–8967. 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.300539–555. 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.71601–1613. 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.1091174–1181. 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.298530–8538. 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.Neuropharmacology131453–474. 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.712413–2426. 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.Cerebellum119–26. 10.1080/147342202753203050
141
KochC. (1984). Cable theory in neurons with active, linearized membranes.Biol. Cybern.5015–33. 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.202166–2174.
144
KozarevaV.MartinC.OsornoT.RudolphS.GuoC.VanderburgC.et al (2021). A transcriptomic atlas of mouse cerebellar cortex comprehensively defines cell types.Nature598214–219. 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.Cerebellum13151–177. 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.15291–307. 10.1046/j.0953-816x.2001.01855.x
147
KumarN. M.GilulaN. B. (1996). The gap junction communication channel.Cell84381–388. 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.182954–2961.
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.4251–61. 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.Neuroscience86533–545. 10.1016/S0306-4522(98)00007-4
152
LemaG. M. C.AuerbachA. (2006). Modes and models of GABAA receptor gating.J. Physiol.572183–200. 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.591899–918. 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.8202–211. 10.1038/nn1393
157
LuH.-W.BalmerT. S.RomeroG. E.TrussellL. O. (2017). Slow AMPAR synaptic transmission is determined by stargazin and glutamate transporters.Neuron9673–80.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.2231565–1587. 10.1007/s00429-017-1568-y
159
MaconochieD. J.ZempelJ. M.SteinbachJ. H. (1994). How quickly can GABAA receptors open?Neuron1261–71. 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.802521–2537. 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.271285–1296. 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.10611388–11393. 10.1073/pnas.0901290106
164
MarkramH.MullerE.RamaswamyS.ReimannM. W.AbdellahM.SanchezC. A.et al (2015). Reconstruction and simulation of neocortical microcircuitry.Cell163456–492. 10.1016/j.cell.2015.09.029
165
MarrD. (1969). A theory of cerebellar cortex.J. Physiol.202437–470. 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.31205–1211. 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.468593–607. 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.Neuron38433–445. 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.563369–378. 10.1113/jphysiol.2004.075028
176
MiyazakiT.YamasakiM.TanakaK. F.WatanabeM. (2021). Compartmentalized input–output organization of lugaro cells in the cerebellar cortex.Neuroscience46289–105. 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.Neuron12529–540. 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.114131–150. 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.66220–245. 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.961476–1489. 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.21981–2035. 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.95686–699. 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), 1–43. 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.60243–260. 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.25702–713. 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.424348–4361. 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.Nature4501245–1248. 10.1038/nature05995
198
RasminskyM. (1980). Ephaptic transmission between single nerve fibres in the spinal nerve roots of dystrophic mice.J. Physiol.305151–169. 10.1113/jphysiol.1980.sp013356
199
RinaldoL.HanselC. (2010). Ataxias and cerebellar dysfunction: Involvement of synaptic plasticity deficits?Funct. Neurol.25135–139.
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), 261–283. 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.Nature4571015–1018. 10.1038/nature07604
204
RusakovD. A. (2001). The role of perisynaptic glial sheaths in glutamate spillover and extracellular Ca(2+) depletion.Biophys. J.811947–1959.
205
RusakovD. A.KullmannD. M. (1998). Extrasynaptic glutamate diffusion in the hippocampus: Ultrastructural constraints, uptake, and receptor activation.J. Neurosci.183158–3170. 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.Science337749–753. 10.1126/science.1221140
207
SaftenkuE. E. (2005). Modeling of slow glutamate diffusion and AMPA receptor activation in the cerebellar glomerulus.J. Theor. Biol.234363–382. 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.33424–429. 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.97248–263. 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.1041823–1828. 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.20236–260. 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.326878–6893. 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.1131743–1751. 10.1152/jn.00843.2014
217
ShannonP. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks.Genome Res.132498–2504. 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.Nature355163–166. 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.557843–861. 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.Cerebellum7517–522. 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.20114–122. 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.737–68. 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.Neuron54121–136. 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.Neuron60469–476. 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.393353–373.
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.343193–3209. 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.5955341–5357. 10.1113/JP274252
240
TaoL.NicholsonC. (2004). Maximum geometrical hindrance to diffusion in brain extracellular space surrounding uniformly spaced convex cells.J. Theor. Biol.22959–68. 10.1016/j.jtbi.2004.03.003
241
TaylorC. P.DudekF. E. (1982). Synchronous neural afterdischarges in rat hippocampal slices without active chemical synapses.Science218810–812. 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.10821–835. 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.Cerebellum10637–646. 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.97901–911. 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.Neuron67435–451. 10.1016/j.neuron.2010.06.028
253
VoogdJ.GlicksteinM. (1998). The anatomy of the cerebellum.Trends Neurosci.21370–375. 10.1016/s0166-2236(98)01318-6
254
VoŽehF. (2015). Jan evangelista purkyně and the cerebellum then and now.Physiol. Res.64:S567–S584.
255
VyklickýL. (1993). Calcium-mediated modulation of N-methyl-D-aspartate (n.d.) responses in cultured rat hippocampal neurones.J. Physiol.470575–600. 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.Neuron32301–313.
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.61671–682. 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.Neuron87781–796. 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.369070–9083. 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.954584–4596. 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.1132524–2536. 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
© 2022 Masoli, Rizza, Tognolina, Prestori and D’Angelo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Francesca Prestori, francesca.prestori@unipv.itEgidio D’Angelo, egidio.dangelo@unipv.it
†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.