Linking dynamics of the inhibitory network to the input structure

Networks of inhibitory interneurons are found in many distinct classes of biological systems. Inhibitory interneurons govern the dynamics of principal cells and are likely to be critically involved in the coding of information. In this theoretical study, we describe the dynamics of a generic inhibitory network in terms of low-dimensional, simplified rate models. We study the relationship between the structure of external input applied to the network and the patterns of activity arising in response to that stimulation. We found that even a minimal inhibitory network can generate a great diversity of spatio-temporal patterning including complex bursting regimes with non-trivial ratios of burst firing. Despite the complexity of these dynamics, the network’s response patterns can be predicted from the rankings of the magnitudes of external inputs to the inhibitory neurons. This type of invariant dynamics is robust to noise and stable in densely connected networks with strong inhibitory coupling. Our study predicts that the response dynamics generated by an inhibitory network may provide critical insights about the temporal structure of the sensory input it receives.     

Firing responses of bursting neurons with delayed feedback

Thalamic neurons, which play important roles in the genesis of rhythmic activities of the brain, show various bursting behaviors, particularly modulated by complex thalamocortical feedback via cortical neurons. As a first step to explore this complex neural system and focus on the effects of the feedback on the bursting behavior, a simple loop structure delayed in time and scaled by a coupling strength is added to a recent mean-field model of bursting neurons. Depending on the coupling strength and delay time, the modeled neurons show two distinct response patterns: one entrained to the unperturbed bursting frequency of the neurons and one entrained to the resonant frequency of the loop structure. Transitions between these two patterns are explored in the model’s parameter space via extensive numerical simulations. It is found that at a fixed loop delay, there is a critical coupling strength at which the dominant response frequency switches from the unperturbed bursting frequency to the loop-induced one. Furthermore, alternating occurrence of these two response frequencies is observed when the delay varies at fixed coupling strength. The results demonstrate that bursting is coupled with feedback to yield new dynamics, which will provide insights into such effects in more complex neural systems.     

Coemergence of regularity and complexity during neural network development

With the growing recognition that rhythmic and oscillatory patterns are widespread in the brain and play important roles in all aspects of the function of our nervous system, there has been a resurgence of interest in neuronal synchronized bursting activity. Here, we were interested in understanding the development of synchronized bursts as information-bearing neuronal activity patterns. For that, we have monitored the morphological organization and spontaneous activity of neuronal networks cultured on multielectrode-arrays during their self-executed evolvement from a mixture of dissociated cells into an active network. Complex collective network electrical activity evolved from sporadic firing patterns of the single neurons. On the system (network) level, the activity was marked by bursting events with interneuronal synchronization and nonarbitrary temporal ordering. We quantified these individual-to-collective activity transitions using newly-developed system level quantitative measures of time series regularity and complexity. We found that individual neuronal activity before synchronization was characterized by high regularity and low complexity. During neuronal wiring, there was a transient period of reorganization marked by low regularity, which then leads to coemergence of elevated regularity and functional (nonstochastic) complexity. We further investigated the morphology-activity interplay by modeling artificial neuronal networks with different topological organizations and connectivity schemes. The simulations support our experimental results by showing increased levels of complexity of neuronal activity patterns when neurons are wired up and organized in clusters (similar to mature real networks), as well as network-level activity regulation once collective activity forms.     

Emergence of Slow-Switching Assemblies in Structured Neuronal Networks

Unraveling the interplay between connectivity and spatio-temporal dynamics in neuronal networks is a key step to advance our understanding of neuronal information processing. Here we investigate how particular features of network connectivity underpin the propensity of neural networks to generate slow-switching assembly (SSA) dynamics, i.e., sustained epochs of increased firing within assemblies of neurons which transition slowly between different assemblies throughout the network. We show that the emergence of SSA activity is linked to spectral properties of the asymmetric synaptic weight matrix. In particular, the leading eigenvalues that dictate the slow dynamics exhibit a gap with respect to the bulk of the spectrum, and the associated Schur vectors exhibit a measure of block-localization on groups of neurons, thus resulting in coherent dynamical activity on those groups. Through simple rate models, we gain analytical understanding of the origin and importance of the spectral gap, and use these insights to develop new network topologies with alternative connectivity paradigms which also display SSA activity. Specifically, SSA dynamics involving excitatory and inhibitory neurons can be achieved by modifying the connectivity patterns between both types of neurons. We also show that SSA activity can occur at multiple timescales reflecting a hierarchy in the connectivity, and demonstrate the emergence of SSA in small-world like networks. Our work provides a step towards understanding how network structure (uncovered through advancements in neuroanatomy and connectomics) can impact on spatio-temporal neural activity and constrain the resulting dynamics.     

Changing excitation and inhibition in simulated neural networks: effects on induced bursting behavior

 The development of synchronous bursting in neuronal ensembles represents an important change in network behavior. To determine the influences on development of such synchronous bursting behavior we study the dynamics of small networks of sparsely connected excitatory and inhibitory neurons using numerical simulations. The synchronized bursting activities in networks evoked by background spikes are investigated. Specifically, patterns of bursting activity are examined when the balance between excitation and inhibition on neuronal inputs is varied and the fraction of inhibitory neurons in the network is changed. For quantitative comparison of bursting activities in networks, measures of the degree of synchrony are used. We demonstrate how changes in the strength of excitation on inputs of neurons can be compensated by changes in the strength of inhibition without changing the degree of synchrony in the network. The effects of changing several network parameters on the network activity are analyzed and discussed. These changes may underlie the transition of network activity from normal to potentially pathologic (e.g., epileptic) states. Received: 21 May 2002 / Accepted in revised form: 3 December 2002 / Published online: 7 March 2003 Correspondence to: P. Kudela (e-mail: pkudela@jhmi.edu) Acknowledgements. This research was supported by NIH grant NS 38958.     

Crustacean motor pattern generator networks

Crustacean motor pattern-generating networks have played central roles in understanding the cellular and network bases of rhythmic motor patterns for over half a century. We review here the four best investigated of these systems: the stomatogastric, ventilatory, cardiac, and swimmeret systems. Generally applicable observations arising from this work include (1) neurons with active, endogenous cell properties (endogenous bursting, postinhibitory rebound, plateau potentials), (2) nonhierarchical (distributed) network synaptic connectivity patterns characterized by high levels of inter-neuronal connections, (3) nonspiking neurons and graded transmitter release, (4) multiple modulatory inputs, (5) networks that produce multiple patterns and have flexible boundaries, and (6) peripheral properties (proprioceptive feedback loops, low-frequency muscle filtering) playing an important role in motor pattern generation or expression.     

Voltage-sensitive ion channels in rhythmic motor systems

Voltage-sensitive ionic currents shape both the firing properties of neurons and their synaptic integration within neural networks that drive rhythmic motor patterns. Persistent sodium currents underlie rhythmic bursting in respiratory neurons. H-type pacemaker currents can act as leak conductances in spinal motoneurons, and also control long-term modulation of synaptic release at the crayfish neuromuscular junction. Calcium currents travel in rostro-caudal waves with motoneuron activity in the spinal cord. Potassium currents control spike width and burst duration in many rhythmic motor systems. We are beginning to identify the genes that underlie these currents.     

Synchronous and asynchronous bursting states: role of intrinsic neural dynamics

Brain signals such as local field potentials often display gamma-band oscillations (30-70 Hz) in a variety of cognitive tasks. These oscillatory activities possibly reflect synchronization of cell assemblies that are engaged in a cognitive function. A type of pyramidal neurons, i.e., chattering neurons, show fast rhythmic bursting (FRB) in the gamma frequency range, and may play an active role in generating the gamma-band oscillations in the cerebral cortex. Our previous phase response analyses have revealed that the synchronization between the coupled bursting neurons significantly depends on the bursting mode that is defined as the number of spikes in each burst. Namely, a network of neurons bursting through a Ca(2+)-dependent mechanism exhibited sharp transitions between synchronous and asynchronous firing states when the neurons exchanged the bursting mode between singlet, doublet and so on. However, whether a broad class of bursting neuron models commonly show such a network behavior remains unclear. Here, we analyze the mechanism underlying this network behavior using a mathematically tractable neuron model. Then we extend our results to a multi-compartment version of the NaP current-based neuron model and prove a similar tight relationship between the bursting mode changes and the network state changes in this model. Thus, the synchronization behavior couples tightly to the bursting mode in a wide class of networks of bursting neurons.     

Pacemaker neurons and neuronal networks: an integrative view

Rhythmically active neuronal networks give rise to rhythmic motor activities but also to seemingly non-rhythmic behaviors such as sleep, arousal, addiction, memory and cognition. Many of these networks contain pacemaker neurons. The ability of these neurons to generate bursts of activity intrinsically lies in voltage- and time-dependent ion fluxes resulting from a dynamic interplay among ion channels, second messenger pathways and intracellular Ca2+ concentrations, and is influenced by neuromodulators and synaptic inputs. This complex intrinsic and extrinsic modulation of pacemaker activity exerts a dynamic effect on network activity. The nonlinearity of bursting activity might enable pacemaker neurons to facilitate the onset of excitatory states or to synchronize neuronal ensembles--an interactive process that is intimately regulated by synaptic and modulatory processes.     

Modelling Feedback Excitation,Pacemaker Properties and Sensory Switching of Electrically Coupled Brainstem Neurons Controlling Rhythmic Activity

What cellular and network properties allow reliable neuronal rhythm generation or firing that can be started and stopped by brief synaptic inputs? We investigate rhythmic activity in an electrically-coupled population of brainstem neurons driving swimming locomotion in young frog tadpoles, and how activity is switched on and off by brief sensory stimulation. We build a computational model of 30 electrically-coupled conditional pacemaker neurons on one side of the tadpole hindbrain and spinal cord. Based on experimental estimates for neuron properties, population sizes, synapse strengths and connections, we show that: long-lasting, mutual, glutamatergic excitation between the neurons allows the network to sustain rhythmic pacemaker firing at swimming frequencies following brief synaptic excitation; activity persists but rhythm breaks down without electrical coupling; NMDA voltage-dependency doubles the range of synaptic feedback strengths generating sustained rhythm. The network can be switched on and off at short latency by brief synaptic excitation and inhibition. We demonstrate that a population of generic Hodgkin-Huxley type neurons coupled by glutamatergic excitatory feedback can generate sustained asynchronous firing switched on and off synaptically. We conclude that networks of neurons with NMDAR mediated feedback excitation can generate self-sustained activity following brief synaptic excitation. The frequency of activity is limited by the kinetics of the neuron membrane channels and can be stopped by brief inhibitory input. Network activity can be rhythmic at lower frequencies if the neurons are electrically coupled. Our key finding is that excitatory synaptic feedback within a population of neurons can produce switchable, stable, sustained firing without synaptic inhibition.     

Dynamics of Time Delay-Induced Multiple Synchronous Behaviors in Inhibitory Coupled Neurons

The inhibitory synapse can induce synchronous behaviors different from the anti-phase synchronous behaviors, which have been reported in recent studies. In the present paper, synchronous behaviors are investigated in the motif model composed of reciprocal inhibitory coupled neurons with endogenous bursting and time delay. When coupling strength is weak, synchronous behavior appears at a single interval of time delay within a bursting period. When coupling strength is strong, multiple synchronous behaviors appear at different intervals of time delay within a bursting period. The different bursting patterns of synchronous behaviors, and time delays and coupling strengths that can induce the synchronous bursting patterns can be well interpreted by the dynamics of the endogenous bursting pattern of isolated neuron, which is acquired by the fast-slow dissection method, combined with the inhibitory coupling current. For an isolated neuron, when a negative impulsive current with suitable strength is applied at different phases of the bursting, multiple different bursting patterns can be induced. For a neuron in the motif, the inhibitory coupling current, of which the application time and strength is modulated by time delay and coupling strength, can cause single or multiple synchronous firing patterns like the negative impulsive current when time delay and coupling strength is suitable. The difference compared to the previously reported multiple synchronous behaviors that appear at time delays wider than a period of the endogenous firing is discussed. The results present novel examples of synchronous behaviors in the neuronal network with inhibitory synapses and provide a reasonable explanation.     

The locust olfactory system as a case study for modeling dynamics of neurobiological networks: from discrete time neurons to continuous time neurons

Both chaotic and periodic activities are observed in networks of the central nervous systems. We choose the locust olfactory system as a good case study to analyze the relationships between networks' structure and the types of dynamics involved in coding mechanisms. In our modeling approach, we first build a fully connected recurrent network of synchronously updated McCulloch and Pitts neurons (MC-P type). In order to measure the use of the temporal dimension in the complex spatio-temporal patterns produced by the networks, we have defined an index the Normalized Euclidian Distance NED. We find that for appropriate parameters of input and connectivity, when adding some strong connections to the initial random synaptic matrices, it was easy to get the emergence of both robust oscillations and distributed synchrony in the spatiotemporal patterns. Then, in order to validate the MC-P model as a tool for analysis for network properties, we examine the dynamic behavior of networks of continuous time model neuron (Izhikevitch Integrate and Fire model -IFI-), implementing the same network characteristics. In both models, similarly to biological PN, the activity of excitatory neurons are phase-locked to different cycles of oscillations which remind the ones of the local field potential (LFP), and nevertheless exhibit dynamic behavior complex enough to be the basis of spatio-temporal codes.     

Capturing the bursting dynamics of a two-cell inhibitory network using a one-dimensional map

Out-of-phase bursting is a functionally important behavior displayed by central pattern generators and other neural circuits. Understanding this complex activity requires the knowledge of the interplay between the intrinsic cell properties and the properties of synaptic coupling between the cells. Here we describe a simple method that allows us to investigate the existence and stability of anti-phase bursting solutions in a network of two spiking neurons, each possessing a T-type calcium current and coupled by reciprocal inhibition. We derive a one-dimensional map which fully characterizes the genesis and regulation of anti-phase bursting arising from the interaction of the T-current properties with the properties of synaptic inhibition. This map is the burst length return map formed as the composition of two distinct one-dimensional maps that are each regulated by a different set of model parameters. Although each map is constructed using the properties of a single isolated model neuron, the composition of the two maps accurately captures the behavior of the full network. We analyze the parameter sensitivity of these maps to determine the influence of both the intrinsic cell properties and the synaptic properties on the burst length, and to find the conditions under which multistability of several bursting solutions is achieved. Although the derivation of the map relies on a number of simplifying assumptions, we discuss how the principle features of this dimensional reduction method could be extended to more realistic model networks. Action Editor: John Rinzel     

Variable neuronal participation in stereotypic motor programs

To what extent are motor networks underlying rhythmic behaviors rigidly hard-wired versus fluid and dynamic entities? Do the members of motor networks change from moment-to-moment or from motor program episode-to-episode? These are questions that can only be addressed in systems where it is possible to monitor the spiking activity of networks of neurons during the production of motor programs. We used large-scale voltage-sensitive dye (VSD) imaging followed by Independent Component Analysis spike-sorting to examine the extent to which the neuronal network underlying the escape swim behavior of Tritonia diomedea is hard-wired versus fluid from a moment-to-moment perspective. We found that while most neurons were dedicated to the swim network, a small but significant proportion of neurons participated in a surprisingly variable manner. These neurons joined the swim motor program late, left early, burst only on some cycles or skipped cycles of the motor program. We confirmed that this variable neuronal participation was not due to effects of the VSD by finding such neurons with intracellular recording in dye-free saline. Further, these neurons markedly varied their level of participation in the network from swim episode-to-episode. The generality of such unreliably bursting neurons was confirmed by their presence in the rhythmic escape networks of two other molluscan species, Tritonia festiva and Aplysia californica. Our observations support a view that neuronal networks, even those underlying rhythmic and stereotyped motor programs, may be more variable in structure than widely appreciated.     

Oscillations and Synchrony in Large-scale Cortical Network Models

Intrinsic neuronal and circuit properties control the responses of large ensembles of neurons by creating spatiotemporal patterns of activity that are used for sensory processing, memory formation, and other cognitive tasks. The modeling of such systems requires computationally efficient single-neuron models capable of displaying realistic response properties. We developed a set of reduced models based on difference equations (map-based models) to simulate the intrinsic dynamics of biological neurons. These phenomenological models were designed to capture the main response properties of specific types of neurons while ensuring realistic model behavior across a sufficient dynamic range of inputs. This approach allows for fast simulations and efficient parameter space analysis of networks containing hundreds of thousands of neurons of different types using a conventional workstation. Drawing on results obtained using large-scale networks of map-based neurons, we discuss spatiotemporal cortical network dynamics as a function of parameters that affect synaptic interactions and intrinsic states of the neurons.     

Transitions between irregular and rhythmic firing patterns in excitatory-inhibitory neuronal networks

Changes in firing patterns are an important hallmark of the functional status of neuronal networks. We apply dynamical systems methods to understand transitions between irregular and rhythmic firing in an excitatory-inhibitory neuronal network model. Using the geometric theory of singular perturbations, we systematically reduce the full model to a simpler set of equations, one that can be studied analytically. The analytic tools are used to understand how an excitatory-inhibitory network with a fixed architecture can generate both activity patterns for possibly different values of the intrinsic and synaptic parameters. These results are applied to a recently developed model for the subthalamopallidal network of the basal ganglia. The results suggest that an increase in correlated activity, corresponding to a pathological state, may be due to an increased level of inhibition from the striatum to the inhibitory GPe cells along with an increased ability of the excitatory STN neurons to generate rebound bursts. Action Editor: Carson Chow     

Perception of odors by a nonlinear model of the olfactory bulb

The behavior of the olfactory bulb is modeled as a network of interconnected cells with nonlinear dynamics. External inputs from sensory neurons are introduced as perturbations to subsets of cells within the network. We describe the attractors of the system and show how they can be classified and ordered according to their varying degrees of symmetry. By studying networks of attractors in the system's phase space, it is shown how different perturbations may evoke specific switches between various patterns of behavior. This ensures that different odors, even if present at extremely low concentrations, are able to evoke a specific spatio-temporal behavior in the olfactory bulb, permitting their unique perception. The model incorporates many of the processes proposed to mediate perception, such as the topographic organisation of sensory systems, destabilization of cortex by sensory input and synchronisation between neurons. It is also consistent with the character of the olfactory electroencephalogram.     

Self-sustained asynchronous irregular states and Up–Down states in thalamic,cortical and thalamocortical networks of nonlinear integrate-and-fire neurons

Randomly-connected networks of integrate-and-fire (IF) neurons are known to display asynchronous irregular (AI) activity states, which resemble the discharge activity recorded in the cerebral cortex of awake animals. However, it is not clear whether such activity states are specific to simple IF models, or if they also exist in networks where neurons are endowed with complex intrinsic properties similar to electrophysiological measurements. Here, we investigate the occurrence of AI states in networks of nonlinear IF neurons, such as the adaptive exponential IF (Brette-Gerstner-Izhikevich) model. This model can display intrinsic properties such as low-threshold spike (LTS), regular spiking (RS) or fast-spiking (FS). We successively investigate the oscillatory and AI dynamics of thalamic, cortical and thalamocortical networks using such models. AI states can be found in each case, sometimes with surprisingly small network size of the order of a few tens of neurons. We show that the presence of LTS neurons in cortex or in thalamus, explains the robust emergence of AI states for relatively small network sizes. Finally, we investigate the role of spike-frequency adaptation (SFA). In cortical networks with strong SFA in RS cells, the AI state is transient, but when SFA is reduced, AI states can be self-sustained for long times. In thalamocortical networks, AI states are found when the cortex is itself in an AI state, but with strong SFA, the thalamocortical network displays Up and Down state transitions, similar to intracellular recordings during slow-wave sleep or anesthesia. Self-sustained Up and Down states could also be generated by two-layer cortical networks with LTS cells. These models suggest that intrinsic properties such as adaptation and low-threshold bursting activity are crucial for the genesis and control of AI states in thalamocortical networks.     

Population behavior and self-organization in the genesis of spontaneous rhythmic activity by developing spinal networks

During development spinal networks generate recurring episodes of rhythmic bursting that can be recorded from motoneurons and interneurons. Optical imaging has identified a set of propriospinal interneurons that may be important in the production of this activity. These neurons are rhythmically active, are recurrently interconnected and have powerful projections to motoneurons. The excitability of this propriospinal network is depressed by activity and recovers in the interval between episodes. These and other observations have been formulated into a qualitative model in which population behavior and self-organization are responsible for the spontaneous activity generated by developing spinal networks.     

大鼠海马癫痫电网络重建中爆发式放电神经元的活动

本文探讨双侧海马(hippoeampus,HPC)神经网络中爆发式放电神经元(bursting-firing neurons,BFN)的活动规律及其与海马癫痫网络重建的关系。实验用雄性SD大鼠140只(150-250 g),急性强直电刺激(60 Hz,2 s,0.4-0.6 mA)右后背HPC CAl区(acute tetanization of the posterior dorsal hippocampus,ATPDH),同步记录同侧或对侧前背HPC单位放电和深部电图;强直电刺激右前背HPC(acute tetanization of the anterior dorsal hippocampus,AT-ADH),同步记录双侧前背HPC单位放电。实验共记录了13.8％(19/138)双侧前背HPC的BFN,其中13个为刺激诱发性BFN,6个为自发性BFN。强直电刺激引起的诱发反应包括:(1)ATPDH明显调制同侧前背HPC的BFN,产生规则的节律性爆发式放电,刺激后串内动作电位间期(bursting interspike interval,BISI)减小(P<0.001);(2)AT-PDH引起对侧前背HPC的BFN出现抑制后轻度调制效应,刺激后动作电位间期(interspike interval,ISI)增大(P<0.001);(3)ATADH后易化对侧前背HPC的自发性BFN节律,增加ISI(P<0.001)和IBI(P=0.01);(4)ATPDH诱导双侧前背HPC的BFN产生规则的节律性爆发式放电,伴有同步或非同步性网络癫痫的形成。上述实验结果提示,ATPDH沿同侧HPC长轴,跨大脑半球诱发前背HPC单个BFN的形成,其节律性爆