Modulation of the dynamics of cerebellar Purkinje cells through the interaction of excitatory and inhibitory feedforward pathways

The dynamics of cerebellar neuronal networks is controlled by the underlying building blocks of neurons and synapses between them. For which, the computation of Purkinje cells (PCs), the only output cells of the cerebellar cortex, is implemented through various types of neural pathways interactively routing excitation and inhibition converged to PCs. Such tuning of excitation and inhibition, coming from the gating of specific pathways as well as short-term plasticity (STP) of the synapses, plays a dominant role in controlling the PC dynamics in terms of firing rate and spike timing. PCs receive cascade feedforward inputs from two major neural pathways: the first one is the feedforward excitatory pathway from granule cells (GCs) to PCs; the second one is the feedforward inhibition pathway from GCs, via molecular layer interneurons (MLIs), to PCs. The GC-PC pathway, together with short-term dynamics of excitatory synapses, has been a focus over past decades, whereas recent experimental evidence shows that MLIs also greatly contribute to controlling PC activity. Therefore, it is expected that the diversity of excitation gated by STP of GC-PC synapses, modulated by strong inhibition from MLI-PC synapses, can promote the computation performed by PCs. However, it remains unclear how these two neural pathways are interacted to modulate PC dynamics. Here using a computational model of PC network installed with these two neural pathways, we addressed this question to investigate the change of PC firing dynamics at the level of single cell and network. We show that the nonlinear characteristics of excitatory STP dynamics can significantly modulate PC spiking dynamics mediated by inhibition. The changes in PC firing rate, firing phase, and temporal spike pattern, are strongly modulated by these two factors in different ways. MLIs mainly contribute to variable delays in the postsynaptic action potentials of PCs while modulated by excitation STP. Notably, the diversity of synchronization and pause response in the PC network is governed not only by the balance of excitation and inhibition, but also by the synaptic STP, depending on input burst patterns. Especially, the pause response shown in the PC network can only emerge with the interaction of both pathways. Together with other recent findings, our results show that the interaction of feedforward pathways of excitation and inhibition, incorporated with synaptic short-term dynamics, can dramatically regulate the PC activities that consequently change the network dynamics of the cerebellar circuit.     

High-frequency oscillations and sequence generation in two-population models of hippocampal region CA1

Hippocampal sharp wave/ripple oscillations are a prominent pattern of collective activity, which consists of a strong overall increase of activity with superimposed (140 − 200 Hz) ripple oscillations. Despite its prominence and its experimentally demonstrated importance for memory consolidation, the mechanisms underlying its generation are to date not understood. Several models assume that recurrent networks of inhibitory cells alone can explain the generation and main characteristics of the ripple oscillations. Recent experiments, however, indicate that in addition to inhibitory basket cells, the pattern requires in vivo the activity of the local population of excitatory pyramidal cells. Here, we study a model for networks in the hippocampal region CA1 incorporating such a local excitatory population of pyramidal neurons. We start by investigating its ability to generate ripple oscillations using extensive simulations. Using biologically plausible parameters, we find that short pulses of external excitation triggering excitatory cell spiking are required for sharp/wave ripple generation with oscillation patterns similar to in vivo observations. Our model has plausible values for single neuron, synapse and connectivity parameters, random connectivity and no strong feedforward drive to the inhibitory population. Specifically, whereas temporally broad excitation can lead to high-frequency oscillations in the ripple range, sparse pyramidal cell activity is only obtained with pulse-like external CA3 excitation. Further simulations indicate that such short pulses could originate from dendritic spikes in the apical or basal dendrites of CA1 pyramidal cells, which are triggered by coincident spike arrivals from hippocampal region CA3. Finally we show that replay of sequences by pyramidal neurons and ripple oscillations can arise intrinsically in CA1 due to structured connectivity that gives rise to alternating excitatory pulse and inhibitory gap coding; the latter denotes phases of silence in specific basket cell groups, which induce selective disinhibition of groups of pyramidal neurons. This general mechanism for sequence generation leads to sparse pyramidal cell and dense basket cell spiking, does not rely on synfire chain-like feedforward excitation and may be relevant for other brain regions as well.     

Functional consequences of correlated excitatory and inhibitory conductances in cortical networks

Neurons in the neocortex receive a large number of excitatory and inhibitory synaptic inputs. Excitation and inhibition dynamically balance each other, with inhibition lagging excitation by only few milliseconds. To characterize the functional consequences of such correlated excitation and inhibition, we studied models in which this correlation structure is induced by feedforward inhibition (FFI). Simple circuits show that an effective FFI changes the integrative behavior of neurons such that only synchronous inputs can elicit spikes, causing the responses to be sparse and precise. Further, effective FFI increases the selectivity for propagation of synchrony through a feedforward network, thereby increasing the stability to background activity. Last, we show that recurrent random networks with effective inhibition are more likely to exhibit dynamical network activity states as have been observed in vivo. Thus, when a feedforward signal path is embedded in such recurrent network, the stabilizing effect of effective inhibition creates an suitable substrate for signal propagation. In conclusion, correlated excitation and inhibition support the notion that synchronous spiking may be important for cortical processing.     

Mesoscopic Patterns of Neural Activity Support Songbird Cortical Sequences

Time-locked sequences of neural activity can be found throughout the vertebrate forebrain in various species and behavioral contexts. From “time cells” in the hippocampus of rodents to cortical activity controlling movement, temporal sequence generation is integral to many forms of learned behavior. However, the mechanisms underlying sequence generation are not well known. Here, we describe a spatial and temporal organization of the songbird premotor cortical microcircuit that supports sparse sequences of neural activity. Multi-channel electrophysiology and calcium imaging reveal that neural activity in premotor cortex is correlated with a length scale of 100 µm. Within this length scale, basal-ganglia–projecting excitatory neurons, on average, fire at a specific phase of a local 30 Hz network rhythm. These results show that premotor cortical activity is inhomogeneous in time and space, and that a mesoscopic dynamical pattern underlies the generation of the neural sequences controlling song.     

Sparse balance: Excitatory-inhibitory networks with small bias currents and broadly distributed synaptic weights

Cortical circuits generate excitatory currents that must be cancelled by strong inhibition to assure stability. The resulting excitatory-inhibitory (E-I) balance can generate spontaneous irregular activity but, in standard balanced E-I models, this requires that an extremely strong feedforward bias current be included along with the recurrent excitation and inhibition. The absence of experimental evidence for such large bias currents inspired us to examine an alternative regime that exhibits asynchronous activity without requiring unrealistically large feedforward input. In these networks, irregular spontaneous activity is supported by a continually changing sparse set of neurons. To support this activity, synaptic strengths must be drawn from high-variance distributions. Unlike standard balanced networks, these sparse balance networks exhibit robust nonlinear responses to uniform inputs and non-Gaussian input statistics. Interestingly, the speed, not the size, of synaptic fluctuations dictates the degree of sparsity in the model. In addition to simulations, we provide a mean-field analysis to illustrate the properties of these networks.     

Top-Down Modulation on Perceptual Decision with Balanced Inhibition through Feedforward and Feedback Inhibitory Neurons

Recent physiological studies have shown that neurons in various regions of the central nervous systems continuously receive noisy excitatory and inhibitory synaptic inputs in a balanced and covaried fashion. While this balanced synaptic input (BSI) is typically described in terms of maintaining the stability of neural circuits, a number of experimental and theoretical studies have suggested that BSI plays a proactive role in brain functions such as top-down modulation for executive control. Two issues have remained unclear in this picture. First, given the noisy nature of neuronal activities in neural circuits, how do the modulatory effects change if the top-down control implements BSI with different ratios between inhibition and excitation? Second, how is a top-down BSI realized via only excitatory long-range projections in the neocortex? To address the first issue, we systematically tested how the inhibition/excitation ratio affects the accuracy and reaction times of a spiking neural circuit model of perceptual decision. We defined an energy function to characterize the network dynamics, and found that different ratios modulate the energy function of the circuit differently and form two distinct functional modes. To address the second issue, we tested BSI with long-distance projection to inhibitory neurons that are either feedforward or feedback, depending on whether these inhibitory neurons do or do not receive inputs from local excitatory cells, respectively. We found that BSI occurs in both cases. Furthermore, when relying on feedback inhibitory neurons, through the recurrent interactions inside the circuit, BSI dynamically and automatically speeds up the decision by gradually reducing its inhibitory component in the course of a trial when a decision process takes too long.     

Turning On and Off with Excitation: The Role of Spike-Timing Asynchrony and Synchrony in Sustained Neural Activity

Delay-related sustained activity in the prefrontal cortex of primates, a neurological analogue of working memory, has been proposed to arise from synaptic interactions in local cortical circuits. The implication is that memories are coded by spatially localized foci of sustained activity. We investigate the mechanisms by which sustained foci are initiated, maintained, and extinguished by excitation in networks of Hodgkin-Huxley neurons coupled with biophysical spatially structured synaptic connections. For networks with a balance between excitation and inhibition, a localized transient stimulus robustly initiates a localized focus of activity. The activity is then maintained by recurrent excitatory AMPA-like synapses. We find that to maintain the focus, the firing must be asynchronous. Consequently, inducing transient synchrony through an excitatory stimulus extinguishes the sustained activity. Such a monosynaptic excitatory turn-off mechanism is compatible with the working memory being wiped clean by an efferent copy of the motor command. The activity that codes working memories may be structured so that the motor command is both the read-out and a direct clearing signal. We show examples of data that is compatible with our theory.     

Numerical simulation of motility patterns of the small bowel. 1. formulation of a mathematical model

A complete mathematical model of the periodic myoelectrical activity of a functional unit of the small intestine is presented. Based on real morphological and electrophysiological data, the model assumes that: the functional unit is an electromyogenic syncytium; the kinetics of L-type Ca2+, T-type Ca2+, Ca2+-activated K+, voltage dependent K+and Cl-channels determine the electrical activity of the functional unit; the enteric nervous system is satisfactorily represented by an efferent cholinergic neuron that provides an excitatory input to the functional unit through receptor-linked L-type Ca2+channels and by an afferent pathway composed of the primary and secondary sensory neurons; the dynamics of propagation of the wave of depolarization along the unmyelinated nerve axons satisfy the Hodgkin-Huxley model; the electrical activity of the neural soma reflects the interaction of N-type Ca2+channels, Ca2+-activated K+and voltage dependent Na+, K+and Cl-channels; the smooth muscle syncytium of the locus is a null-dimensional contractile system. With the proposed model the dynamics of active force generation are determined entirely by the concentration of cytosolic calcium. The model describes: the mechanical excitation of the free nerve endings of the mechanoreceptor of the receptive field of the pathway; the electrical processes of the propagation of excitation along the afferent and efferent neural circuits; the chemical mechanisms of nerve-pulse transmission at the synaptic zones; the slow wave and bursting type electrical activity; cytosolic calcium concentration; the dynamics of active force generation. Numerical simulations have shown that the model can display different electrical patterns and mechanical responses of the locus. The results show good qualitative and quantitative agreement with the results of experiments conducted on the small intestine.     

Information coding and oscillatory activity in synfire neural networks with and without inhibitory coupling

When a population spike (pulse-packet) propagates through a feedforward network with random excitatory connections, it either evolves to a sustained stable level of synchronous activity or fades away (Diesmann et al. in Nature 402:529-533 1999; Cateau and Fukai Neur Netw 14:675-685 2001). Here I demonstrate that in the presence of noise, the probability of the survival of the pulse-packet (or, equivalently, the firing rate of output neurons) reflects the intensity of the input. Furthermore, inhibitory coupling between layers can result in quasi- periodic alternation between several levels of firing activity. These results are obtained by analyzing the evolution of pulse-packet activity as a Markov chain. For the Markov chain analysis, the output of the chain is a linear mapping of the input into a lower-dimensional space, and the eigenvalues and eigenvectors of the transition matrix determine the dynamics of the evolution. Synchronous propagation of firing activity in successive pools of neurons are simulated in networks of integrate-and-fire and compartmental model neurons, and, consistent with the discrete Markov process, the activation of each pool is observed to be predominantly dependent upon the number of cells that fired in the previous pool. Simulation results agree with the numerical solutions of the Markov model. When inhibitory coupling between layers are included in the Markov model, some eigenvalues become complex numbers, implying oscillatory dynamics. The quasiperiodic dynamics is validated with simulation with leaky integrate-and-fire neurons. The networks demonstrate different modes of quasiperiodic activity as the inhibition or excitation parameters of the network are varied.     

Analysis of an optimal control model of multi-joint arm movements

 In this paper, we propose a model of biological motor control for generation of goal-directed multi-joint arm movements, and study the formation of muscle control inputs and invariant kinematic features of movements. The model has a hierarchical structure that can determine the control inputs for a set of redundant muscles without any inverse computation. Calculation of motor commands is divided into two stages, each of which performs a transformation of motor commands from one coordinate system to another. At the first level, a central controller in the brain accepts instructions from higher centers, which represent the motor goal in the Cartesian space. The controller computes joint equilibrium trajectories and excitation signals according to a minimum effort criterion. At the second level, a neural network in the spinal cord translates the excitation signals and equilibrium trajectories into control commands to three pairs of antagonist muscles which are redundant for a two-joint arm. No inverse computation is required in the determination of individual muscle commands. The minimum effort controller can produce arm movements whose dynamic and kinematic features are similar to those of voluntary arm movements. For fast movements, the hand approaches a target position along a near-straight path with a smooth bell-shaped velocity. The equilibrium trajectories in X and Y show an ‘N’ shape, but the end-point equilibrium path zigzags around the hand path. Joint movements are not always smooth. Joint reversal is found in movements in some directions. The excitation signals have a triphasic (or biphasic) pulse pattern, which leads to stereotyped triphasic (or biphasic) bursts in muscle control inputs, and a dynamically modulated joint stiffness. There is a fixed sequence of muscle activation from proximal muscles to distal muscles. The order is preserved in all movements. For slow movements, it is shown that a constant joint stiffness is necessary to produce a smooth movement with a bell-shaped velocity. Scaled movements can be reproduced by varying the constraints on the maximal level of excitation signals according to the speed of movement. When the inertial parameters of the arm are altered, movement trajectories can be kept invariant by adjusting the pulse height values, showing the ability to adapt to load changes. These results agree with a wide range of experimental observations on human voluntary movements. Received: 4 December 1995 / Accepted in revised form: 17 September 1996     

Perisomatic GABA release and thalamocortical integration onto neocortical excitatory cells are regulated by neuromodulators

Neuromodulators such as acetylcholine, serotonin, and noradrenaline are powerful regulators of neocortical activity. Although it is well established that cortical inhibition is the target of these modulations, little is known about their effects on GABA release from specific interneuron types. This knowledge is necessary to gain a mechanistic understanding of the actions of neuromodulators because different interneuron classes control specific aspects of excitatory cell function. Here, we report that GABA release from fast-spiking (FS) cells, the most prevalent interneuron subtype in neocortex, is robustly inhibited following activation of muscarinic, serotonin, adenosine, and GABA(B) receptors--an effect that regulates FS cell control of excitatory neuron firing. The potent muscarinic inhibition of GABA release from FS cells suppresses thalamocortical feedforward inhibition. This is supplemented by the muscarinic-mediated depolarization of thalamo-recipient excitatory neurons and the nicotinic enhancement of thalamic input onto these neurons to promote thalamocortical excitation.     

Synchronization without oscillatory neurons

Experimental results suggest that neurons in the cortex synchronize their action potentials on the millisecond time scale. More importantly this binding expresses functional relationships between the neurons. A model of neuronal interactions is proposed in which simultaneous discharges of neurons develop through specialized synaptic circuits. As an important prerequisite for this synchronization it is demonstrated that SynFire chains, generating different levels of excitation, propagate their activity waves at distinct velocities. Two chains were coupled by excitatory synapses and their activity was initiated at different times. Due to synaptic interactions, activity in the earlier-initiated chain accelerates propagation in the other chain until the two activity waves are synchronized. Compared with several neural network models with oscillatory units, physiologically more plausible neurons are simulated. It is still under debate whether neurons in the cortex show oscillatory dischargesper se. In particular, a high rate of noise relative to very weak synaptic gains cannot impair our results in the neural network simulations.     

Modified thalamocortical model: A step towards more understanding of the functional contribution of astrocytes to epilepsy

It is evident that the cortex plays a primary role in seizure generation. At the same time, various experimental results clearly confirm that thalamic neurons are also actively involved in seizure generation and spreading. On the other hand, recent neurophysiologic findings suggest that astrocytes regulate dynamically the synaptic activity in neuronal networks. Therefore, in the present study, the thalamocortical neural population model (TCPM) is modified by embedding into the model the functional role of astrocytes in the regulation of synaptic transmission. Using the modified TCPM (MTCPM) we examined the hypothesis that one of the possible causes of neural hypersynchronization is the dysfunction of astrocytes in the regulatory feedback loop. Then, two MTCPMs are coupled via excitatory synapses and the astrocytes are also coupled together through gap junctions. Utilizing the MTCPM and CMTCPM, the transition from normal to malfunctioned states is analyzed using a dynamical system approach. In this way, the hypothesis is investigated and it is demonstrated that the healthy astrocytes provide feedback control to regulate neural activity. That is, the astrocytes compensate to a large extent variations in the coupling between neural populations and maintain the balance between the excitation and inhibition levels. However, the malfunctioned astrocytes are no longer able to regulate and/or compensate the excessive increase of the inter-population coupling strength. As a consequence, disruption of the signaling function of astrocytes could contribute to the neuronal hyperexcitability and generating epileptiform activity. These results suggest that astrocytes might be one of the potential targets for the treatment of epilepsy.     

Determination of the fraction of active inputs required by a neuron to fire

What fraction of the inputs to a neuron in the primary visual cortex (V1) need to be active for that neuron to reach its firing threshold? The paper describes a numerical method for estimating the selectivity of visual neurons, in terms of the required fraction of active excitatory inputs, from standard data produced by intracellular electro-physiological recordings. The method also provides an estimate of the relative strength of the feedforward inhibition in a push-pull model of the inputs to V1 simple cells. The method is tested on two V1 cells described in Carandini and Ferster [Carandini, M., Ferster, D., 2000. Membrane potential and firing rate in cat primary visual cortex. J. Neurosci. 20, 470-484]. The results indicate that the maximum strength of feedforward inhibition is around 30% of the maximum strength of feedforward excitation. The two V1 neurons investigated fire if more than around 40% of their excitatory LGN inputs are active.     

A simple computer model of excitable synaptically connected neurons

 The space-lumped two-variable neuron model is studied. Extension of the neural model by adding a simple synaptic current allows the demonstration of neural interactions. The production of synchronous burst activity in this simple two-neuron excitatory loop is modeled, including the influence of random background excitatory input. The ability of the neuron model to integrate inputs spatially and temporally is shown. Two refractory periods after stimuli were identified and their role in burst cessation is demonstrated. Our findings show that simple neural units without long-lasting membrane processes are capable of generating long lasting patterns of activity. The results of simulation of simple background activity suggest that an increase in background activity tends to cause decreased activity of the network. This phenomenon, as well as the existence of two refractory periods, allows for burst cessation without inhibition in this simple model. Received: 6 December 1996/Accepted in revised form: 7 April 1997     

A minimal mechanistic model for temporal signal processing in the lateral geniculate nucleus

The receptive fields of cells in the lateral geniculate nucleus (LGN) are shaped by their diverse set of impinging inputs: feedforward synaptic inputs stemming from retina, and feedback inputs stemming from the visual cortex and the thalamic reticular nucleus. To probe the possible roles of these feedforward and feedback inputs in shaping the temporal receptive-field structure of LGN relay cells, we here present and investigate a minimal mechanistic firing-rate model tailored to elucidate their disparate features. The model for LGN relay ON cells includes feedforward excitation and inhibition (via interneurons) from retinal ON cells and excitatory and inhibitory (via thalamic reticular nucleus cells and interneurons) feedback from cortical ON and OFF cells. From a general firing-rate model formulated in terms of Volterra integral equations, we derive a single delay differential equation with absolute delay governing the dynamics of the system. A freely available and easy-to-use GUI-based MATLAB version of this minimal mechanistic LGN circuit model is provided. We particularly investigate the LGN relay-cell impulse response and find through thorough explorations of the model’s parameter space that both purely feedforward models and feedback models with feedforward excitation only, can account quantitatively for previously reported experimental results. We find, however, that the purely feedforward model predicts two impulse response measures, the time to first peak and the biphasic index (measuring the relative weight of the rebound phase) to be anticorrelated. In contrast, the models with feedback predict different correlations between these two measures. This suggests an experimental test assessing the relative importance of feedforward and feedback connections in shaping the impulse response of LGN relay cells.     

Generating oscillatory bursts from a network of regular spiking neurons without inhibition

Avian nucleus isthmi pars parvocellularis (Ipc) neurons are reciprocally connected with the layer 10 (L10) neurons in the optic tectum and respond with oscillatory bursts to visual stimulation. Our in vitro experiments show that both neuron types respond with regular spiking to somatic current injection and that the feedforward and feedback synaptic connections are excitatory, but of different strength and time course. To elucidate mechanisms of oscillatory bursting in this network of regularly spiking neurons, we investigated an experimentally constrained model of coupled leaky integrate-and-fire neurons with spike-rate adaptation. The model reproduces the observed Ipc oscillatory bursting in response to simulated visual stimulation. A scan through the model parameter volume reveals that Ipc oscillatory burst generation can be caused by strong and brief feedforward synaptic conductance changes. The mechanism is sensitive to the parameter values of spike-rate adaptation. In conclusion, we show that a network of regular-spiking neurons with feedforward excitation and spike-rate adaptation can generate oscillatory bursting in response to a constant input.     

Epileptiform activity in a neocortical network: a mathematical model

 A simple mathematical model describing the generation and propagation of epileptiform activity in a cerebral cortical network is presented. The model consists of a system of nonlinear delay differential equations. Physiological properties are taken into account as nonlinear transmission of signals at the synapse, temporal and spatial summation of incoming signals at the soma, active membrane characteristics, and dendritic and axonal propagation times. The influence of the connectivity and the temporal parameters on the oscillatory properties of the model is studied. The computer simulations are in agreement with experimental observations in cortical networks: whereas a weak excitatory or strong inhibitory synaptic connection strength produces a stationary status with short-lasting responses to external stimuli, increases in excitation or decreases in inhibition induce spontaneous and stimulus-evoked rhythmic discharges. Synaptic burst-like activity is observed only for an intermediate range of excitatory and inhibitory connection strengths and external inputs. The form and duration of the bursts can also be controlled by the temporal parameters. The results demonstrate that relatively simple mathematical equations are sufficient to model some of the network properties underlying the generation and propagation of epileptiform activity. Received: 2 October 2000 / Accepted in revised form: 4 March 2001     

Activity-dependent matching of excitatory and inhibitory inputs during refinement of visual receptive fields

The receptive field (RF) of single visual neurons undergoes progressive refinement during development. It remains largely unknown how the excitatory and inhibitory inputs on single developing neurons are refined in a coordinated manner to allow the formation of functionally correct circuits. Using whole-cell voltage-clamp recording from Xenopus tectal neurons, we found that RFs determined by excitatory and inhibitory inputs in more mature tectal neurons are spatially matched, with each spot stimulus evoking balanced synaptic excitation and inhibition. This emerges during development through a gradual reduction in the RF size and a transition from disparate to matched topography of excitatory and inhibitory inputs to the tectal neurons. Altering normal spiking activity of tectal neurons by either blocking or elevating GABA(A) receptor activity significantly impeded the developmental reduction and topographic matching of RFs. Thus, appropriate inhibitory activity is essential for the coordinated refinement of excitatory and inhibitory connections.