Intrinsic bursting enhances the robustness of a neural network model of sequence generation by avian brain area HVC

Avian brain area HVC is known to be important for the production of birdsong. In zebra finches, each RA-projecting neuron in HVC emits a single burst of spikes during a song motif. The population of neurons is activated in a precisely timed, stereotyped sequence. We propose a model of these burst sequences that relies on two hypotheses. First, we hypothesize that the sequential order of bursting is reflected in the excitatory synaptic connections between neurons. Second, we propose that the neurons are intrinsically bursting, so that burst duration is set by cellular properties. Our model generates burst sequences similar to those observed in HVC. If intrinsic bursting is removed from the model, burst sequences can also be produced. However, they require more fine-tuning of synaptic strengths, and are therefore less robust. In our model, intrinsic bursting is caused by dendritic calcium spikes, and strong spike frequency adaptation in the soma contributes to burst termination.     

Anti-phase,asymmetric and aperiodic oscillations in excitable cells—I. Coupled bursters

I seek to explain phenomena observed in simulations of populations of gap junction-coupled bursting cells by studying the dynamics of identical pairs. I use a simplified model for pancreatic β-cells and decompose the system into fast (spike-generating) and slow subsystems to show how bifurcations of the fast subsystem affect bursting behavior. When coupling is weak, the spikes are not in phase but rather are anti-phase, asymmetric or quasi-periodic. These solutions all support bursting with smaller amplitude spikes than the in-phase case, leading to increased burst period. A key geometrical feature underlying this is that the in-phase periodic solution branch terminates in a homoclinic orbit. The same mechanism also provides a model for bursting as an emergent property of populations; cells which are not intrinsic bursters can burst when coupled. This phenomenon is enhanced when symmetry is broken by making the cells differ in a parameter.     

Inhibition modifies the effects of slow calcium-activated potassium channels on epileptiform activity in a neuronal network model

Generation of epileptiform activity typically results from a change in the balance between network excitation and inhibition. Experimental evidence indicates that alterations of either synaptic activity or intrinsic membrane properties can produce increased network excitation. The slow Ca²⁺-activated K⁺ currents (sI_AHP) are important modulators of neuronal firing rate and excitability and have important established and potential roles in epileptogenesis. While the effects of changes in sI_AHP on individual neuronal excitability are readily studied and well established, the effects of such changes on network behavior are less well known. The experiments here utilize a defined small network model of multicompartment pyramidal cells and an inhibitory interneuron to study the effects of changes in sI_AHP on network behavior. The benefits of this model system include the ability to observe activity in all cells in a network and the effects of interactions of multiple simultaneous influences. In the model with no inhibitory interneuron, increasing sI_AHP results in progressively decreasing burst activity. Adding an inhibitory interneuron changes the observed effects; at modest inhibitory strengths, increasing sI_AHP in all network neurons actually results in increased network bursting (except at very high values). The duration of the burst activity is influenced by the length of delay in a feedback loop, with longer loops resulting in more prolonged bursting. These observations illustrate that the study of potential antiepileptogenic membrane effects must be extended to realistic networks. Network inhibition can dramatically alter the observations seen in pure excitatory networks.     

Conflicting effects of excitatory synaptic and electric coupling on the dynamics of square-wave bursters

Using two-cell and 50-cell networks of square-wave bursters, we studied how excitatory coupling of individual neurons affects the bursting output of the network. Our results show that the effects of synaptic excitation vs. electrical coupling are distinct. Increasing excitatory synaptic coupling generally increases burst duration. Electrical coupling also increases burst duration for low to moderate values, but at sufficiently strong values promotes a switch to highly synchronous bursts where further increases in electrical or synaptic coupling have a minimal effect on burst duration. These effects are largely mediated by spike synchrony, which is determined by the stability of the in-phase spiking solution during the burst. Even when both coupling mechanisms are strong, one form (in-phase or anti-phase) of spike synchrony will determine the burst dynamics, resulting in a sharp boundary in the space of the coupling parameters. This boundary exists in both two cell and network simulations. We use these results to interpret the effects of gap-junction blockers on the neuronal circuitry that underlies respiration.     

Timing regulation in a network reduced from voltage-gated equations to a one-dimensional map

 We discuss a method by which the dynamics of a network of neurons, coupled by mutual inhibition, can be reduced to a one-dimensional map. This network consists of a pair of neurons, one of which is an endogenous burster, and the other excitable but not bursting in the absence of phasic input. The latter cell has more than one slow process. The reduction uses the standard separation of slow/fast processes; it also uses information about how the dynamics on the slow manifold evolve after a finite amount of slow time. From this reduction we obtain a one-dimensional map dependent on the parameters of the original biophysical equations. In some parameter regimes, one can deduce that the original equations have solutions in which the active phase of the originally excitable cell is constant from burst to burst, while in other parameter regimes it is not. The existence or absence of this kind of regulation corresponds to qualitatively different dynamics in the one-dimensional map. The computations associated with the reduction and the analysis of the dynamics includes the use of coordinates that parameterize by time along trajectories, and “singular Poincaré maps” that combine information about flows along a slow manifold with information about jumps between branches of the slow manifold. Received: 19 May 1997 / Revised version: 6 April 1998     

Co-existent activity patterns in inhibitory neuronal networks with short-term synaptic depression

A network of two neurons mutually coupled through inhibitory synapses that display short-term synaptic depression is considered. We show that synaptic depression expands the number of possible activity patterns that the network can display and allows for co-existence of different patterns. Specifically, the network supports different types of n-m anti-phase firing patterns, where one neuron fires n spikes followed by the other neuron firing m spikes. When maximal synaptic conductances are identical, n-n anti-phase firing patterns are obtained and there are conductance intervals over which different pairs of these solutions co-exist. The multitude of n-m anti-phase patterns and their co-existence are not found when the synapses are non-depressing. Geometric singular perturbation methods for dynamical systems are applied to the original eight-dimensional model system to derive a set of one-dimensional conditions for the existence and co-existence of different anti-phase solutions. The generality and validity of these conditions are demonstrated through numerical simulations utilizing the Hodgkin-Huxley and Morris-Lecar neuronal models.     

Noise-Induced Coherence and Network Oscillations in a Reduced Bursting Model

The dynamics of the Hindmarsh-Rose (HR) model of bursting thalamic neurons is reduced to a system of two linear differential equations that retains the subthreshold resonance properties of the HR model. Introducing a reset mechanism after a threshold crossing, we turn this system into a resonant integrate-and-fire (RIF) model. Using Monte-Carlo simulations and mathematical analysis, we examine the effects of noise and the subthreshold dynamic properties of the RIF model on the occurrence of coherence resonance (CR). Synchronized burst firing occurs in a network of such model neurons with excitatory pulse-coupling. The coherence level of the network oscillations shows a stochastic resonance-like dependence on the noise level. Stochastic analysis of the equations shows that the slow recovery from the spike-induced inhibition is crucial in determining the frequencies of the CR and the subthreshold resonance in the original HR model. In this particular type of CR, the oscillation frequency strongly depends on the intrinsic time scales but changes little with the noise intensity. We give analytical quantities to describe this CR mechanism and illustrate its influence on the emerging network oscillations. We discuss the profound physiological roles this kind of CR may have in information processing in neurons possessing a subthreshold resonant frequency and in generating synchronized network oscillations with a frequency that is determined by intrinsic properties of the neurons. PACS 05.45.-a, 05.40.Ca, 87.18.Sn, 87.19     

Bursting Reverberation as a Multiscale Neuronal Network Process Driven by Synaptic Depression-Facilitation

Neuronal networks can generate complex patterns of activity that depend on membrane properties of individual neurons as well as on functional synapses. To decipher the impact of synaptic properties and connectivity on neuronal network behavior, we investigate the responses of neuronal ensembles from small (5–30 cells in a restricted sphere) and large (acute hippocampal slice) networks to single electrical stimulation: in both cases, a single stimulus generated a synchronous long-lasting bursting activity. While an initial spike triggered a reverberating network activity that lasted 2–5 seconds for small networks, we found here that it lasted only up to 300 milliseconds in slices. To explain this phenomena present at different scales, we generalize the depression-facilitation model and extracted the network time constants. The model predicts that the reverberation time has a bell shaped relation with the synaptic density, revealing that the bursting time cannot exceed a maximum value. Furthermore, before reaching its maximum, the reverberation time increases sub-linearly with the synaptic density of the network. We conclude that synaptic dynamics and connectivity shape the mean burst duration, a property present at various scales of the networks. Thus bursting reverberation is a property of sufficiently connected neural networks, and can be generated by collective depression and facilitation of underlying functional synapses.     

Coherency and connectivity in oscillating neural networks: linear partialization analysis

 This paper studies the relation between the functional synaptic connections between two artificial neural networks and the correlation of their spiking activities. The model neurons had realistic non-oscillatory dynamic properties and the networks showed oscillatory behavior as a result of their internal synaptic connectivity. We found that both excitation and inhibition cause phase locking of the oscillating activities. When the two networks excite each other the oscillations synchronize with zero phase lag, whereas mutual inhibition between the networks resulted in an anti-phase (half period phase difference) synchronization. Correlations between the activities of the two networks can also be caused by correlated external inputs driving the systems (common input). Our analysis shows that when the networks exhibit oscillatory behavior and the rate of the common input is smaller than a characteristic network oscillator frequency, the cross-correlation functions between the activities of two systems still carry information about the mutual synaptic connectivity. This information can be retrieved with linear partialization, removing the influence of the common input. We further explored the network responses to periodic external input. We found that when the input is of a frequency smaller than a certain threshold, the network responds with bursts at the same frequency as the input. Above the threshold, the network responds with a fraction of the input frequency. This frequency threshold, characterizing the oscillatory properties of the network, is also found to determine the limit to which linear partialization works. Received: 20 October 1995 / Accepted in revised form: 20 May 1996     

Feedback Inhibition and Throughput Properties of an Integrate-and-Fire-or-Burst Network Model of Retinogeniculate Transmission

Computational modeling has played an important role in the dissection of the biophysical basis of rhythmic oscillations in thalamus that are associated with sleep and certain forms of epilepsy. In contrast, the dynamic filter properties of thalamic relay nuclei during states of arousal are not well understood. Here we present a modeling and simulation study of the throughput properties of the visually driven dorsal lateral geniculate nucleus (dLGN) in the presence of feedback inhibition from the perigeniculate nucleus (PGN). We employ thalamocortical (TC) and thalamic reticular (RE) versions of a minimal integrate-and-fire-or-burst type model and a one-dimensional, two-layered network architecture. Potassium leakage conductances control the neuromodulatory state of the network and eliminate rhythmic bursting in the presence of spontaneous input (i.e., wake up the network). The aroused dLGN/PGN network model is subsequently stimulated by spatially homogeneous spontaneous retinal input or spatio-temporally patterned input consistent with the activity of X-type retinal ganglion cells during full-field or drifting grating visual stimulation. The throughput properties of this visually-driven dLGN/PGN network model are characterized and quantified as a function of stimulus parameters such as contrast, temporal frequency, and spatial frequency. During low-frequency oscillatory full-field stimulation, feedback inhibition from RE neurons often leads to TC neuron burst responses, while at high frequency tonic responses dominate. Depending on the average rate of stimulation, contrast level, and temporal frequency of modulation, the TC and RE cell bursts may or may not be phase-locked to the visual stimulus. During drifting-grating stimulation, phase-locked bursts often occur for sufficiently high contrast so long as the spatial period of the grating is not small compared to the synaptic footprint length, i.e., the spatial scale of the network connectivity.     

Modeling of Substance P and 5-HT Induced Synaptic Plasticity in the Lamprey Spinal CPG: Consequences for Network Pattern Generation

Consequences of synaptic plasticity in the lamprey spinal CPG are analyzed by means of simulations. This is motivated by the effects substance P (a tachykinin) and serotonin (5-hydroxytryptamin; 5-HT) have on synaptic transmission in the locomotor network. Activity-dependent synaptic depression and potentiation have recently been shown experimentally using paired intracellular recordings. Although normally activity-dependent plasticity presumably does not contribute to the patterning of network activity, this changes in the presence of the neuromodulators substance P and 5-HT, which evoke significant plasticity. Substance P can induce a faster and larger depression of inhibitory connections but potentiation of excitatory inputs, whereas 5-HT induces facilitation of both inhibitory and excitatory inputs. Changes in the amplitude of the first postsynaptic potential are also seen. These changes could thus be a potential mechanism underlying the modulatory role these substances have on the rhythmic network activity.The aim of the present study has been to implement the activity dependent synaptic depression and facilitation induced by substance P and 5-HT into two alternative models of the lamprey spinal locomotor network, one relying on reciprocal inhibition for bursting and one in which each hemicord is capable of oscillations. The consequences of the plasticity of inhibitory and excitatory connections are then explored on the network level.In the intact spinal cord, tachykinins and 5-HT, which can be endogenously released, increase and decrease the frequency of the alternating left-right burst pattern, respectively. The frequency decreasing effect of 5-HT has previously been explained based on its conductance decreasing effect on K _Ca underlying the postspike afterhyperpolarization (AHP). The present simulations show that short-term synaptic plasticity may have strong effects on frequency regulation in the lamprey spinal CPG. In the network model relying on reciprocal inhibition, the observed effects substance P and 5-HT have on network behavior (i.e., a frequency increase and decrease respectively) can to a substantial part be explained by their effects on the total extent and time dynamics of synaptic depression and facilitation. The cellular effects of these substances will in the 5-HT case further contribute to its network effect.     

From simple to complex oscillatory behaviour: analysis of bursting in a multiply regulated biochemical system

We analyze the transition from simple to complex oscillatory behaviour in a three-variable biochemical system that consists of the coupling in series of two autocatalytic enzyme reactions. Complex periodic behaviour occurs in the form of bursting in which clusters of spikes are separated by phases of relative quiescence. The generation of such temporal patterns is investigated by a series of complementary approaches. The dynamics of the system is first cast into two different time-scales, and one of the variables is taken as a slowly-varying parameter influencing the behaviour of the two remaining variables. This analysis shows how complex oscillations develop from simple periodic behaviour and accounts for the existence of various modes of bursting as well as for the dependence of the number of spikes per period on key parameters of the model. We further reduce the number of variables by analyzing bursting by means of one-dimensional return maps obtained from the time evolution of the three-dimensional system. The analysis of a related piecewise linear map allows for a detailed understanding of the complex sequence leading from a bursting pattern with p spikes to a pattern with p + 1 spikes per period. We show that this transition possesses properties of self-similarity associated with the occurrence of more and more complex patterns of bursting. In addition to bursting, period-doubling bifurcations leading to chaos are observed, as in the differential system, when the piecewise-linear map becomes nonlinear.     

Cell Type-Specific Separation of Subicular Principal Neurons during Network Activities

The hippocampal output structure, the subiculum, expresses two major memory relevant network rhythms, sharp wave ripple and gamma frequency oscillations. To this date, it remains unclear how the two distinct types of subicular principal cells, intrinsically bursting and regular spiking neurons, participate in these two network rhythms. Using concomitant local field potential and intracellular recordings in an in vitro mouse model that allows the investigation of both network rhythms, we found a cell type-specific segregation of principal neurons into participating intrinsically bursting and non-participating regular spiking cells. However, if regular spiking cells were kept at a more depolarized level, they did participate in a specific manner, suggesting a potential bimodal working model dependent on the level of excitation. Furthermore, intrinsically bursting and regular spiking cells exhibited divergent intrinsic membrane and synaptic properties in the active network. Thus, our results suggest a cell-type-specific segregation of principal cells into two separate groups during network activities, supporting the idea of two parallel streams of information processing within the subiculum.     

Responses of a bursting pacemaker to excitation reveal spatial segregation between bursting and spiking mechanisms

Central pattern generators (CPGs) frequently include bursting neurons that serve as pacemakers for rhythm generation. Phase resetting curves (PRCs) can provide insight into mechanisms underlying phase locking in such circuits. PRCs were constructed for a pacemaker bursting complex in the pyloric circuit in the stomatogastric ganglion of the lobster and crab. This complex is comprised of the Anterior Burster (AB) neuron and two Pyloric Dilator (PD) neurons that are all electrically coupled. Artificial excitatory synaptic conductance pulses of different strengths and durations were injected into one of the AB or PD somata using the Dynamic Clamp. Previously, we characterized the inhibitory PRCs by assuming a single slow process that enabled synaptic inputs to trigger switches between an up state in which spiking occurs and a down state in which it does not. Excitation produced five different PRC shapes, which could not be explained with such a simple model. A separate dendritic compartment was required to separate the mechanism that generates the up and down phases of the bursting envelope (1) from synaptic inputs applied at the soma, (2) from axonal spike generation and (3) from a slow process with a slower time scale than burst generation. This study reveals that due to the nonlinear properties and compartmentalization of ionic channels, the response to excitation is more complex than inhibition.     

Ionic mechanisms underlying spontaneous CA1 neuronal firing in Ca2+-free solution

Hippocampal CA1 neurons exposed to zero-[Ca(2+)] solutions can generate periodic spontaneous synchronized activity in the absence of synaptic function. Experiments using hippocampal slices showed that, after exposure to zero-[Ca(2+)](0) solution, CA1 pyramidal cells depolarized 5-10 mV and started firing spontaneous action potentials. Spontaneous single neuron activity appeared in singlets or was grouped into bursts of two or three action potentials. A 16-compartment, 23-variable cable model of a CA1 pyramidal neuron was developed to study mechanisms of spontaneous neuronal bursting in a calcium-free extracellular solution. In the model, five active currents (a fast sodium current, a persistent sodium current, an A-type transient potassium current, a delayed rectifier potassium current, and a muscarinic potassium current) are included in the somatic compartment. The model simulates the spontaneous bursting behavior of neurons in calcium-free solutions. The mechanisms underlying several aspects of bursting are studied, including the generation of triplet bursts, spike duration, burst termination, after-depolarization behavior, and the prolonged inactive period between bursts. We show that the small persistent sodium current can play a key role in spontaneous CA1 activity in zero-calcium solutions. In particular, it is necessary for the generation of an after-depolarizing potential and prolongs both individual bursts and the interburst interval.     

A multivariate population density model of the dLGN/PGN relay

Using a population density approach we study the dynamics of two interacting collections of integrate-and-fire-or-burst (IFB) neurons representing thalamocortical (TC) cells from the dorsal lateral geniculate nucleus (dLGN) and thalamic reticular (RE) cells from the perigeniculate nucleus (PGN). Each population of neurons is described by a multivariate probability density function that satisfies a conservation equation with appropriately defined probability fluxes and boundary conditions. The state variables of each neuron are the membrane potential and the inactivation gating variable of the low-threshold Ca²⁺ current I_T. The synaptic coupling of the populations and external excitatory drive are modeled by instantaneous jumps in the membrane potential of postsynaptic neurons. The population density model is validated by comparing its response to time-varying retinal input to Monte Carlo simulations of the corresponding IFB network composed of 100 to 1000 cells per population. In the absence of retinal input, the population density model exhibits rhythmic bursting similar to the 7 to 14 Hz oscillations associated with slow wave sleep that require feedback inhibition from RE to TC cells. When the TC and RE cell potassium leakage conductances are adjusted to represent cholingergic neuromodulation and arousal of the network, rhythmic bursting of the probability density model may either persists or be eliminated depending on the number of excitatory (TC to RE) or inhibitory (RE to TC) connections made by each presynaptic cell. When the probability density model is stimulated with constant retinal input (10–100 spikes/sec), a wide range of responses are observed depending on cellular parameters and network connectivity. These include asynchronous burst and tonic spikes, sleep spindle-like rhythmic bursting, and oscillations in population firing rate that are distinguishable from sleep spindles due to their amplitude, frequency, or the presence of tonic spikes. In this context of dLGN/PGN network modeling, we find the population density approach using 2,500 mesh points and resolving membrane voltage to 0.7 mV is over 30 times more efficient than 1000-cell Monte Carlo simulations. Action Editor: David Golomb     

Computer simulation of the segmental neural network generating locomotion in lamprey by using populations of network interneurons

Realistic computer simulations of the experimentally established local spinal cord neural network generating swimming in the lamprey have been performed. Populations of network interneurons were used in which cellular properties, like cell size and membrane conductance including voltage dependent ion channels were randomly distributed around experimentally obtained mean values, as were synaptic conductances (kainate/AMPA, NMDA, glycine) and delays. This population model displayed more robust burst activity over a wider frequency range than the more simple subsample model used previously, and the pattern of interneuronal activity was appropriate. The strength of the reciprocal inhibition played a very important role in the regulation of burst frequency, and just by changing the inhibitory bias the entire physiological range could be covered. At the lower frequency range of bursting the segmental excitatory interneurons provide stability as does the activation of voltage dependent NMDA receptors. Spike frequency adaptation by means of summation of afterhyperpolarization (AHP) serves as a major burst terminating factor, and at lower rates the membrane properties conferred by the NMDA receptor activation. The lateral interneurons were not of critical importance for the burst termination. They may, however, be of particular importance for inducing a rapid burst termination during for instance steering and righting reactions. Several cellular factors combine to provide a secure and stable motor pattern in the entire frequency range.     

State-Dependent Effects of Na Channel Noise on Neuronal Burst Generation

We explore the effects of stochastic sodium (Na) channel activation on the variability and dynamics of spiking and bursting in a model neuron. The complete model segregates Hodgin-Huxley-type currents into two compartments, and undergoes applied current-dependent bifurcations between regimes of periodic bursting, chaotic bursting, and tonic spiking. Noise is added to simulate variable, finite sizes of the population of Na channels in the fast spiking compartment.During tonic firing, Na channel noise causes variability in interspike intervals (ISIs). The variance, as well as the sensitivity to noise, depend on the model's biophysical complexity. They are smallest in an isolated spiking compartment; increase significantly upon coupling to a passive compartment; and increase again when the second compartment also includes slow-acting currents. In this full model, sufficient noise can convert tonic firing into bursting.During bursting, the actions of Na channel noise are state-dependent. The higher the noise level, the greater the jitter in spike timing within bursts. The noise makes the burst durations of periodic regimes variable, while decreasing burst length duration and variance in a chaotic regime. Na channel noise blurs the sharp transitions of spike time and burst length seen at the bifurcations of the noise-free model. Close to such a bifurcation, the burst behaviors of previously periodic and chaotic regimes become essentially indistinguishable.We discuss biophysical mechanisms, dynamical interpretations and physiological implications. We suggest that noise associated with finite populations of Na channels could evoke very different effects on the intrinsic variability of spiking and bursting discharges, depending on a biological neuron's complexity and applied current-dependent state. We find that simulated channel noise in the model neuron qualitatively replicates the observed variability in burst length and interburst interval in an isolated biological bursting neuron.     

Integrating neuromorphic action-oriented perceptual inputs to generate a navigation behaviour for a robot

We use neural networks with pointer map architectures to provide simple attentional processing in a robotic task. A pointer map comprises a map of neurons that encode a stimulus. Besides global feedback inhibition, the map receives feedback excitation via a small group of pointer neurons that encode the location of a salient stimulus on the map as a vectorial representation. The pointer neurons are able to apply selective processing to a particular region of the network. The robot uses these properties to manoeuver in relation to an attended object. We implemented a controller composed of two pointer maps, and a motor map. The first pointer map reports the direction of a salient obstacle in a one-dimensional map of distance derived from infrared sensors. The second pointer map reports the direction to potential obstacles in a two-dimensional edge-enhanced image derived from a forward looking CCD-camera. These outputs are applied to a motor map, where they bias the motor control signals issued to the robots wheels, according to navigational intentions.     

Long-term Relationships between Cholinergic Tone, Synchronous Bursting and Synaptic Remodeling

Cholinergic neuromodulation plays key roles in the regulation of neuronal excitability, network activity, arousal, and behavior. On longer time scales, cholinergic systems play essential roles in cortical development, maturation, and plasticity. Presumably, these processes are associated with substantial synaptic remodeling, yet to date, long-term relationships between cholinergic tone and synaptic remodeling remain largely unknown. Here we used automated microscopy combined with multielectrode array recordings to study long-term relationships between cholinergic tone, excitatory synapse remodeling, and network activity characteristics in networks of cortical neurons grown on multielectrode array substrates. Experimental elevations of cholinergic tone led to the abrupt suppression of episodic synchronous bursting activity (but not of general activity), followed by a gradual growth of excitatory synapses over hours. Subsequent blockage of cholinergic receptors led to an immediate restoration of synchronous bursting and the gradual reversal of synaptic growth. Neither synaptic growth nor downsizing was governed by multiplicative scaling rules. Instead, these occurred in a subset of synapses, irrespective of initial synaptic size. Synaptic growth seemed to depend on intrinsic network activity, but not on the degree to which bursting was suppressed. Intriguingly, sustained elevations of cholinergic tone were associated with a gradual recovery of synchronous bursting but not with a reversal of synaptic growth. These findings show that cholinergic tone can strongly affect synaptic remodeling and synchronous bursting activity, but do not support a strict coupling between the two. Finally, the reemergence of synchronous bursting in the presence of elevated cholinergic tone indicates that the capacity of cholinergic neuromodulation to indefinitely suppress synchronous bursting might be inherently limited.