A DISCRETE APPROACH FOR A MODEL OF TEMPORAL LEARNING BY THE CEREBELLUM:IN SILICOCLASSICAL CONDITIONING OF THE EYEBLINK REFLEX

A computational mechanism for unified gain and timing control in the cerebellum

Precise gain and timing control is the goal of cerebellar motor learning. Because the basic neural circuitry of the cerebellum is homogeneous throughout the cerebellar cortex, a single computational mechanism may be used for simultaneous gain and timing control. Although many computational models of the cerebellum have been proposed for either gain or timing control, few models have aimed to unify them. In this paper, we hypothesize that gain and timing control can be unified by learning of the complete waveform of the desired movement profile instructed by climbing fiber signals. To justify our hypothesis, we adopted a large-scale spiking network model of the cerebellum, which was originally developed for cerebellar timing mechanisms to explain the experimental data of Pavlovian delay eyeblink conditioning, to the gain adaptation of optokinetic response (OKR) eye movements. By conducting large-scale computer simulations, we could reproduce some features of OKR adaptation, such as the learning-related change of simple spike firing of model Purkinje cells and vestibular nuclear neurons, simulated gain increase, and frequency-dependent gain increase. These results suggest that the cerebellum may use a single computational mechanism to control gain and timing simultaneously.

Computational models of timing mechanisms in the cerebellar granular layer

A long-standing question in neuroscience is how the brain controls movement that requires precisely timed muscle activations. Studies using Pavlovian delay eyeblink conditioning provide good insight into this question. In delay eyeblink conditioning, which is believed to involve the cerebellum, a subject learns an interstimulus interval (ISI) between the onsets of a conditioned stimulus (CS) such as a tone and an unconditioned stimulus such as an airpuff to the eye. After a conditioning phase, the subject’s eyes automatically close or blink when the ISI time has passed after CS onset. This timing information is thought to be represented in some way in the cerebellum. Several computational models of the cerebellum have been proposed to explain the mechanisms of time representation, and they commonly point to the granular layer network. This article will review these computational models and discuss the possible computational power of the cerebellum.

Development and validation of a neural population model based on the dynamics of a discontinuous membrane potential neuron model

The major goal of this study was to develop a population density based model derived from statistical mechanics based on the dynamics of a discontinuous membrane potential neuron model. A secondary goal was to validate this model by comparing results from a direct simulation approach on the one hand and our population based approach on the other hand. Comparisons between the two approaches in the case of a synaptically uncoupled and a synaptically coupled neural population produced satisfactory qualitative agreement in terms of firing rate and mean membrane potential. Reasonable quantitative agreement was also obtained for these variables in performed simulations. The results of this work based on the dynamics of a discontinuous membrane potential neuron model provide a basis to simulate phenomenologically large-scale neuronal networks with a reasonably short computing time.

New insights offered by a computational model of deep brain stimulation

Deep brain stimulation (DBS) is a standard neurosurgical procedure used to treat motor symptoms in about 5% of patients with Parkinson's disease (PD). Despite the indisputable success of this procedure, the biological mechanisms underlying the clinical benefits of DBS have not yet been fully elucidated. The paper starts with a brief review on the use of DBS to treat PD symptoms. The second section introduces a computational model based on the population density approach and the Izhikevich neuron model. We explain why this model is appropriate for investigating macroscopic network effects and exploring the physiological mechanisms which respond to this treatment strategy (i.e., DBS). Finally, we present new insights into the ways this computational model may help to elucidate the dynamic network effects produced in a cerebral structure when DBS is applied.

Evidence from retractor bulbi EMG for linearized motor control of conditioned nictitating membrane responses

Classical conditioning of nictitating membrane (NM) responses in rabbits is a robust model learning system, and experimental evidence indicates that conditioned responses (CRs) are controlled by the cerebellum. It is unknown whether cerebellar control signals deal directly with the complex nonlinearities of the plant (blink-related muscles and peripheral tissues) or whether the plant is linearized to ensure a simple relation between cerebellar neuronal firing and CR profile. To study this question, the retractor bulbi muscle EMG was recorded with implanted electrodes during NM conditioning. Pooled activity in accessory abducens motoneurons was estimated from spike trains extracted from the EMG traces, and its temporal profile was found to have an approximately Gaussian shape with peak amplitude linearly related to CR amplitude. The relation between motoneuron activity and CR profiles was accurately fitted by a first-order linear filter, with each spike input producing an exponentially decaying impulse response with time constant of order 0.1 s. Application of this first-order plant model to CR data from other laboratories suggested that, in these cases also, motoneuron activity had a Gaussian profile, with time-of-peak close to unconditioned stimulus (US) onset and SD proportional to the interval between conditioned stimulus and US onsets. These results suggest that for conditioned NM responses the cerebellum is presented with a simplified "virtual" plant that is a linearized version of the underlying nonlinear biological system. Analysis of a detailed plant model suggests that one method for linearising the plant would be appropriate recruitment of motor units.

Response linearity determined by recruitment strategy in detailed model of nictitating membrane control

Many models of eyeblink conditioning assume that there is a simple linear relationship between the firing patterns of neurons in the interpositus nucleus and the time course of the conditioned response (CR). However, the complexities of muscle behaviour and plant dynamics call this assumption into question. We investigated the issue by implementing the most detailed model available of the rabbit nictitating membrane response (Bartha and Thompson in Biol Cybern 68:135-143, 1992a and in Biol Cybern 68:145-154, 1992b), in which each motor unit of the retractor bulbi muscle is represented by a Hill-type model, driven by a non-linear activation mechanism designed to reproduce the isometric force measurements of Lennerstrand (J Physiol 236:43-55, 1974). Globe retraction and NM extension are modelled as linked second order systems. We derived versions of the model that used a consistent set of SI units, were based on a physically realisable version of calcium kinetics, and used simulated muscle cross-bridges to produce force. All versions showed similar non-linear responses to two basic control strategies. (1) Rate-coding with no recruitment gave a sigmoidal relation between control signal and amplitude of CR, reflecting the measured relation between isometric muscle force and stimulation frequency. (2) Recruitment of similar strength motor units with no rate coding gave a sublinear relation between control signal and amplitude of CR, reflecting the increase in muscle stiffness produced by recruitment. However, the system response could be linearised by either a suitable combination of rate-coding and recruitment, or by simple recruitment of motor units in order of (exponentially) increasing strength. These plausible control strategies, either alone or in combination, would in effect present the cerebellum with the simplified virtual plant that is assumed in many models of eyeblink conditioning. Future work is therefore needed to determine the extent to which motor neuron firing is in fact linearly related to the nictitating membrane response.