Maintaining phase of the crustacean tri-phasic pyloric rhythm

We construct and analyze a model network of the pyloric rhythm of the crustacean stomatogastric ganglion consisting of an oscillator neuron that inhibits two reciprocally inhibitory follower neurons. We derive analytic expressions that determine the phase of firing of the follower neurons with respect to the oscillator. An important aspect of the model is the inclusion of synapses that exhibit short-term synaptic depression. We show that these type of synapses allow there to be a complicated relationship between the intrinsic properties of the neurons and the synapses between them in determining phase relationships. Our analysis reveals the circumstances and ranges of cycle periods under which these properties work in concert with or independently from one another. In particular, we show that phase maintenance over a range of oscillator periods can be enhanced through the interplay of the two follower neurons if the synapses between these neurons are depressing. Since our model represents the core of the oscillatory pyloric network, the results of our analysis can be compared to experimental data and used to make predictions about the biological network.     

The Activity Phase of Postsynaptic Neurons in a Simplified Rhythmic Network

Many inhibitory rhythmic networks produce activity in a range of frequencies. The relative phase of activity between neurons in these networks is often a determinant of the network output. This relative phase is determined by the interaction between synaptic inputs to the neurons and their intrinsic properties. We show, in a simplified network consisting of an oscillator inhibiting a follower neuron, how the interaction between synaptic depression and a transient potassium current in the follower neuron determines the activity phase of this neuron. We derive a mathematical expression to determine at what phase of the oscillation the follower neuron becomes active. This expression can be used to understand which parameters determine the phase of activity of the follower as the frequency of the oscillator is changed. We show that in the presence of synaptic depression, there can be three distinct frequency intervals, in which the phase of the follower neuron is determined by different sets of parameters. Alternatively, when the synapse is not depressing, only one set of parameters determines the phase of activity at all frequencies.     

龙虾胃肠神经系统的数值分析

利用抑制神经系统的WinnerLess Competition(WLC)模型,通过数值方法分析Mulloney型龙虾胃肠神经系统神经元的电位发放,得到胃研磨囊和幽门神经系统中各个神经元的电位发放和系统的节律变化。结果表明,胃研磨系统内神经元的发放规律显示两侧牙齿和中间牙齿出现切断、挤压和研磨食物等状态,幽门系统内神经元的发放规律显示幽门节律出现依次发放的三个部分。两个神经系统的数值结果,不仅解释了龙虾胃肠神经系统中神经元电位发放与肌肉运动的关系,而且理论再现了龙虾胃肠神经系统的节律变化和实验结果。     

Synchronization and sensitivity enhancement of the Hodgkin-Huxley neurons due to inhibitory inputs

Recent experimental results imply that inhibitory postsynaptic potentials can play a functional role in realizing synchronization of neuronal firing in the brain. In order to examine the relation between inhibition and synchronous firing of neurons theoretically, we analyze possible effects of synchronization and sensitivity enhancement caused by inhibitory inputs to neurons with a biologically realistic model of the Hodgkin-Huxley equations. The result shows that, after an inhibitory spike, the firing probability of a single postsynaptic neuron exposed to random excitatory background activity oscillates with time. The oscillation of the firing probability can be related to synchronous firing of neurons receiving an inhibitory spike simultaneously. Further, we show that when an inhibitory spike input precedes an excitatory spike input, the presence of such preceding inhibition raises the firing probability peak of the neuron after the excitatory input. The result indicates that an inhibitory spike input can enhance the sensitivity of the postsynaptic neuron to the following excitatory spike input. Two neural network models based on these effects on postsynaptic neurons caused by inhibitory inputs are proposed to demonstrate possible mechanisms of detecting particular spatiotemporal spike patterns. Received: 15 April 1999 /Accepted in revised form: 25 November 1999     

Using phase resetting to predict 1:1 and 2:2 locking in two neuron networks in which firing order is not always preserved

Our goal is to understand how nearly synchronous modes arise in heterogenous networks of neurons. In heterogenous networks, instead of exact synchrony, nearly synchronous modes arise, which include both 1:1 and 2:2 phase-locked modes. Existence and stability criteria for 2:2 phase-locked modes in reciprocally coupled two neuron circuits were derived based on the open loop phase resetting curve (PRC) without the assumption of weak coupling. The PRC for each component neuron was generated using the change in synaptic conductance produced by a presynaptic action potential as the perturbation. Separate derivations were required for modes in which the firing order is preserved and for those in which it alternates. Networks composed of two model neurons coupled by reciprocal inhibition were examined to test the predictions. The parameter regimes in which both types of nearly synchronous modes are exhibited were accurately predicted both qualitatively and quantitatively provided that the synaptic time constant is short with respect to the period and that the effect of second order resetting is considered. In contrast, PRC methods based on weak coupling could not predict 2:2 modes and did not predict the 1:1 modes with the level of accuracy achieved by the strong coupling methods. The strong coupling prediction methods provide insight into what manipulations promote near-synchrony in a two neuron network and may also have predictive value for larger networks, which can also manifest changes in firing order. We also identify a novel route by which synchrony is lost in mildly heterogenous networks.     

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.     

Co-variation of ionic conductances supports phase maintenance in stomatogastric neurons

Neuronal networks produce reliable functional output throughout the lifespan of an animal despite ceaseless molecular turnover and a constantly changing environment. Central pattern generators, such as those of the crustacean stomatogastric ganglion (STG), are able to robustly maintain their functionality over a wide range of burst periods. Previous experimental work involving extracellular recordings of the pyloric pattern of the STG has demonstrated that as the burst period varies, the inter-neuronal delays are altered proportionally, resulting in burst phases that are roughly invariant. The question whether spike delays within bursts are also proportional to pyloric period has not been explored in detail. The mechanism by which the pyloric neurons accomplish phase maintenance is currently not obvious. Previous studies suggest that the co-regulation of certain ion channel properties may play a role in governing neuronal activity. Here, we observed in long-term recordings of the pyloric rhythm that spike delays can vary proportionally with burst period, so that spike phase is maintained. We then used a conductance-based model neuron to determine whether co-varying ionic membrane conductances results in neural output that emulates the experimentally observed phenomenon of spike phase maintenance. Next, we utilized a model neuron database to determine whether conductance correlations exist in model neuron populations with highly maintained spike phases. We found that co-varying certain conductances, including the sodium and transient calcium conductance pair, causes the model neuron to maintain a specific spike phase pattern. Results indicate a possible relationship between conductance co-regulation and phase maintenance in STG neurons.     

Neuromodulatory changes in short-term synaptic dynamics may be mediated by two distinct mechanisms of presynaptic calcium entry

Although synaptic output is known to be modulated by changes in presynaptic calcium channels, additional pathways for calcium entry into the presynaptic terminal, such as non-selective channels, could contribute to modulation of short term synaptic dynamics. We address this issue using computational modeling. The neuropeptide proctolin modulates the inhibitory synapse from the lateral pyloric (LP) to the pyloric dilator (PD) neuron, two slow-wave bursting neurons in the pyloric network of the crab Cancer borealis. Proctolin enhances the strength of this synapse and also changes its dynamics. Whereas in control saline the synapse shows depression independent of the amplitude of the presynaptic LP signal, in proctolin, with high-amplitude presynaptic LP stimulation the synapse remains depressing while low-amplitude stimulation causes facilitation. We use simple calcium-dependent release models to explore two alternative mechanisms underlying these modulatory effects. In the first model, proctolin directly targets calcium channels by changing their activation kinetics which results in gradual accumulation of calcium with low-amplitude presynaptic stimulation, leading to facilitation. The second model uses the fact that proctolin is known to activate a non-specific cation current I _MI . In this model, we assume that the MI channels have some permeability to calcium, modeled to be a result of slow conformation change after binding calcium. This generates a gradual increase in calcium influx into the presynaptic terminals through the modulatory channel similar to that described in the first model. Each of these models can explain the modulation of the synapse by proctolin but with different consequences for network activity.     

The role of putative glutamatergic neurons and their connections in the locomotor central pattern generator of the mollusk, Clione limacina

In the pteropod mollusk Clione limacina, locomotor rhythm is produced by the central pattern generator (CPG), due mainly to the activity of interneurons of groups 7 (active in the phase of the dorsal flexion of the wings) and 8 (active in the phase of the ventral flexion). Each of these groups excites the neurons active in the same phase of the locomotor cycle, and inhibits the neurons of the opposite phase. In this work, the nature of connections formed by group 7 interneurons was studied. Riluzole (2-amino-6-trifluoro-methoxybenzothiazole), which is known to inhibit the presynaptic release of glutamate, suppressed the action of the type 7 interneurons onto the follower neurons of the same and of the antagonistic phase of the locomotor cycle. The main pattern of rhythmic activity of CPG with alternation of two phases could be maintained after suppression of inhibitory connections from group 7 interneurons to antagonistic neurons. This suggests redundancy of the mechanisms controlling swimming rhythm generation, which ensures the reliable operation of the system.     

Attentional modulation of firing rate and synchrony in a model cortical network

The response of a neuron in the visual cortex to stimuli of different contrast placed in its receptive field is commonly characterized using the contrast response curve. When attention is directed into the receptive field of a V4 neuron, its contrast response curve is shifted to lower contrast values (Reynolds et al., 2000). The neuron will thus be able to respond to weaker stimuli than it responded to without attention. Attention also increases the coherence between neurons responding to the same stimulus (Fries et al., 2001). We studied how the firing rate and synchrony of a densely interconnected cortical network varied with contrast and how they were modulated by attention. The changes in contrast and attention were modeled as changes in driving current to the network neurons. We found that an increased driving current to the excitatory neurons increased the overall firing rate of the network, whereas variation of the driving current to inhibitory neurons modulated the synchrony of the network. We explain the synchrony modulation in terms of a locking phenomenon during which the ratio of excitatory to inhibitory firing rates is approximately constant for a range of driving current values. We explored the hypothesis that contrast is represented primarily as a drive to the excitatory neurons, whereas attention corresponds to a reduction in driving current to the inhibitory neurons. Using this hypothesis, the model reproduces the following experimental observations: (1) the firing rate of the excitatory neurons increases with contrast; (2) for high contrast stimuli, the firing rate saturates and the network synchronizes; (3) attention shifts the contrast response curve to lower contrast values; (4) attention leads to stronger synchronization that starts at a lower value of the contrast compared with the attend-away condition. In addition, it predicts that attention increases the delay between the inhibitory and excitatory synchronous volleys produced by the network, allowing the stimulus to recruit more downstream neurons. Action Editor: David Golomb     

Synchronized oscillations in the visual cortex — a synergetic model

We present an oscillator network model for the synchronization of oscillatory neuronal activity underlying visual processing. The single neuron is modeled by means of a limit cycle oscillator with an eigenfrequency corresponding to visual stimulation. The eigenfrequency may be time dependent. The mutual coupling strengths are unsymmetrical and activity dependent, and they scatter within the network. Synchronized clusters (groups) of neurons emerge in the network due to the visual stimulation. The different clusters correspond to different visual stimuli. There is no limitation of the number of stimuli. Distinct clusters do not perturb each other, although the coupling strength between all model neurons is of the same order of magnitude. Our analysis is not restricted to weak coupling strength. The scatter of the couplings causes shifts of the cluster frequencies. The model's behavior is compared with the experimental findings. The coupling mechanism is extended in order to model the influence of bicucullin upon the neural network. We additionally investigate repulsive couplings, which lead to constant phase differences between clusters of the same frequency. Finally, we consider the problem of selective attention from the viewpoint of our model.     

Analysis of bursting in a thalamic neuron model

We extend a quantitative model for low-voltage, slow-wave excitability based on the T-type calcium current (Wang et al. 1991) by juxtaposing it with a Hodgkin-Huxley-like model for fast sodium spiking in the high voltage regime to account for the distinct firing modes of thalamic neurons. We employ bifurcation analysis to illustrate the stimulus-response behavior of the full model under both voltage regimes. The model neuron shows continuous sodium spiking when depolarized sufficiently from rest. Depending on the parameters of calcium current inactivation, there are two types of low-voltage responses to a hyperpolarizing current step: a single rebound low threshold spike (LTS) upon release of the step and periodic LTSs. Bursting is seen as sodium spikes ride the LTS crest. In both cases, we analyze the LTS burst response by projecting its trajectory into a fast/slow phase plane. We also use phase plane methods to show that a potassium A-current shifts the threshold for sodium spikes, reducing the number of fast sodium spikes in an LTS burst. It can also annihilate periodic bursting. We extend the previous work of Rose and Hindmarsh (1989a–c) for a thalamic neuron and propose a simpler model for thalamic activity. We consider burst modulation by using a neuromodulator-dependent potassium leakage conductance as a control parameter. These results correspond with experiments showing that the application of certain neurotransmitters can switch firing modes. Received: 18 July 1993/Accepted in revised form: 22 January 1994     

Motor pattern generation of the posterior cardiac plate-pyloric system in the stomatogastric ganglion of the mantis shrimp Squilla oratoria

Activity patterns of the constituent neurons of the posterior cardiac plate-pyloric system in the stomatogastric ganglion of the mantis shrimp Squilla oratoria were studied by recording spontaneous burst discharges intracellularly from neuronal somata. These neurons were identified electrophysiologically, and synaptic connections among them were qualitatively analysed. The posterior cardiac plate constrictor, pyloric constrictor, pyloric dilator and ventricular dilator motoneurons, and the pyloric interneuron were involved in the posterior cardiac plate-pyloric system. All the cell types could produce slow burst-forming potentials which led to repetitive spike discharges. These neurons generated sequentially patterned outputs. Most commonly, the posterior cardiac plate neuron activity was followed by the activity of pyloric constrictor neurons, and then by the activity of pyloric dilator/pyloric interneuron, and ventricular dilator neurons. The motoneurons and interneuron in the posterior cardiac plate-pyloric system were connected to each other either by electrical or by inhibitory chemical synapses, and thus constructed the neural circuit characterized by a wiring diagram which was structurally similar to the pyloric circuit of decapods. The circuitry in the stomatogastric ganglion was strongly conserved during evolution between stomatopods and decapods, despite significant changes in the peripheral structure of the foregut. There were more electrical synapses in stomatopods, and more reciprocal inhibitory synapses in decapods.Abbreviations EJP excitatory junctional potential - IPSP inhibitory postsynaptic potential - CoG commissural ganglion - CPG central pattern generator - ion inferior oesophageal nerve - OG oesophageal ganglion - pcp posterior cardiac plate - son superior oesophageal nerve - STG stomatogastric ganglion - stn stomatogastric nerve - PY pyloric constrictor - PD pyloric dilator - VD ventricular dilator - AB pyloric interneuron - lvn lateral ventricular nerves - tcpm transverse cardiac plate muscle     

A neural network model for trace conditioning

We studied the dynamics of a neural network that has both recurrent excitatory and random inhibitory connections. Neurons started to become active when a relatively weak transient excitatory signal was presented and the activity was sustained due to the recurrent excitatory connections. The sustained activity stopped when a strong transient signal was presented or when neurons were disinhibited. The random inhibitory connections modulated the activity patterns of neurons so that the patterns evolved without recurrence with time. Hence, a time passage between the onsets of the two transient signals was represented by the sequence of activity patterns. We then applied this model to represent the trace eye blink conditioning, which is mediated by the hippocampus. We assumed this model as CA3 of the hippocampus and considered an output neuron corresponding to a neuron in CA1. The activity pattern of the output neuron was similar to that of CA1 neurons during trace eye blink conditioning, which was experimentally observed.     

Avalanches in a Stochastic Model of Spiking Neurons

Neuronal avalanches are a form of spontaneous activity widely observed in cortical slices and other types of nervous tissue, both in vivo and in vitro. They are characterized by irregular, isolated population bursts when many neurons fire together, where the number of spikes per burst obeys a power law distribution. We simulate, using the Gillespie algorithm, a model of neuronal avalanches based on stochastic single neurons. The network consists of excitatory and inhibitory neurons, first with all-to-all connectivity and later with random sparse connectivity. Analyzing our model using the system size expansion, we show that the model obeys the standard Wilson-Cowan equations for large network sizes ( neurons). When excitation and inhibition are closely balanced, networks of thousands of neurons exhibit irregular synchronous activity, including the characteristic power law distribution of avalanche size. We show that these avalanches are due to the balanced network having weakly stable functionally feedforward dynamics, which amplifies some small fluctuations into the large population bursts. Balanced networks are thought to underlie a variety of observed network behaviours and have useful computational properties, such as responding quickly to changes in input. Thus, the appearance of avalanches in such functionally feedforward networks indicates that avalanches may be a simple consequence of a widely present network structure, when neuron dynamics are noisy. An important implication is that a network need not be “critical” for the production of avalanches, so experimentally observed power laws in burst size may be a signature of noisy functionally feedforward structure rather than of, for example, self-organized criticality.     

Bifurcation analysis of a neural network model

This paper describes the analysis of the well known neural network model by Wilson and Cowan. The neural network is modeled by a system of two ordinary differential equations that describe the evolution of average activities of excitatory and inhibitory populations of neurons. We analyze the dependence of the model's behavior on two parameters. The parameter plane is partitioned into regions of equivalent behavior bounded by bifurcation curves, and the representative phase diagram is constructed for each region. This allows us to describe qualitatively the behavior of the model in each region and to predict changes in the model dynamics as parameters are varied. In particular, we show that for some parameter values the system can exhibit long-period oscillations. A new type of dynamical behavior is also found when the system settles down either to a stationary state or to a limit cycle depending on the initial point.     

Synchronized oscillations in the visual cortex — a synergetic model

 We present an oscillator network model for the synchronization of oscillatory neuronal activity underlying visual processing. The single neuron is modeled by means of a limit cycle oscillator with an eigenfrequency corresponding to visual stimulation. The eigenfrequency may be time dependent. The mutual coupling strengths are unsymmetrical and activity dependent, and they scatter within the network. Synchronized clusters (groups) of neurons emerge in the network due to the visual stimulation. The different clusters correspond to different visual stimuli. There is no limitation of the number of stimuli. Distinct clusters do not perturb each other, although the coupling strength between all model neurons is of the same order of magnitude. Our analysis is not restricted to weak coupling strength. The scatter of the couplings causes shifts of the cluster frequencies. The model’s behavior is compared with the experimental findings. The coupling mechanism is extended in order to model the influence of bicucullin upon the neural network. We additionally investigate repulsive couplings, which lead to constant phase differences between clusters of the same frequency. Finally, we consider the problem of selective attention from the viewpoint of our model. Received: 15 February 1995/Accepted in revised form: 18 July 1995     

Rare and spatially segregated release sites mediate a synaptic interaction between two identified network neurons.

Laser-scanning confocal microscopy (LSCM), electron microcopy (EM), and cellular electrophysiology were used in combination to study the structural basis of an inhibitory synapse between two identified neurons of the same network. To achieve this, we examined the chemical inhibitory synapse between identified neurons belonging to the lobster (Homarus gammarus) pyloric network: the pyloric dilator (PD) and the lateral pyloric (LP) neurons. In order to visualize simultaneously these two neurons, we used intrasomatic injection of Lucifer Yellow (LY) in one and rhodamine/horseradish peroxydase (HRP) in the other. Under LSCM, we found only two zones of close apposition in a restricted part of the neuritic tree of the two network neurons. Then, within these two zones, the synaptic release sites were searched using EM. To this end, photoconversion of LY with immunogold and development of HRP with DAB were performed on the previously observed preparations. Structural evidence was found for only one release site per zone. To confirm this result, and because the zones of contact were always segregated in a restricted part of the dendrites, we used laser photoablation to selectively delete, either pre- or postsynaptically, the branches on which the release sites were located. In both cases, such restrictive ablation completely abolished the functional interaction between these neurons. Our results therefore demonstrate that an inhibitory synapse that is essential for the operation of a neural network relies on only very few sites of contact localized in a highly restricted part of each neuron's dendritic arbor.     

A burst-feedback model of fast-phase burst generation during nystagmus

 Vestibular and optokinetic nystagmus are characterized by alternating slow-phase eye rotations that stabilize the retinal image, and fast-phase eye rotations that reset eye position. Nystagmus is coordinated in the brainstem by burst neurons that fire an intense, temporally circumscribed burst that terminates the slow phase and drives the fast phase. This paper demonstrates that such a burst can be generated during nystagmus using a simple neural network model containing only known brainstem neurons and their interconnections. These include the feedback connections of the burst neuron (burst feedback). The burst neuron excites itself directly, and disinhibits itself by inhibiting the pause neuron (positive feedback). It also inhibits itself by inhibiting the vestibular neuron (negative feedback). The burst neuron begins to fire after its inhibitory bias is overcome by excitation from the vestibular neuron, and burst neuron positive feedback then produces an intense burst with an abrupt onset. The burst causes the vestibular and pause neurons to pause. The benefit of the pause neuron loop is that it contributes to burst neuron positive feedback when it is needed at burst onset, but that contribution is eliminated when the pause neuron pauses and opens the loop. The burst can then terminate, with an offset duration proportional to burst amplitude, under the control of burst neuron self-excitation and inhibitory bias. Model neuron behavior is comparable to that of real brainstem neurons. Synchronized bursts can be produced over the population of burst neurons in a distributed version of the network. Variability in connection weights in the distributed network results in variability in prelude activity among burst neurons that is similar to the spread in lead observed for real burst neurons during nystagmus. Received: 11 April 1996 / Accepted in revised form: 6 August 1996     

Existence and stability criteria for phase-locked modes in ring neural networks based on the spike time resetting curve method

We developed a systematic and consistent mathematical approach to predicting 1:1 phase-locked modes in ring neural networks of spiking neurons based on the open loop spike time resetting curve (STRC) and its almost equivalent counterpart—the phase resetting curve (PRC). The open loop STRCs/PRCs were obtained by injecting into an isolated model neuron a triangular shaped time-dependent stimulus current closely resembling an actual synaptic input. Among other advantages, the STRC eliminates the confusion regarding the undefined phase for stimuli driving the neuron outside of the unperturbed limit cycle. We derived both open loop PRC and STRC-based existence and stability criteria for 1:1 phase-locked modes developed in ring networks of spiking neurons. Our predictions were in good agreement with the closed loop numerical simulations. Intuitive graphical methods for predicting phase-locked modes were also developed both for half-centers and for larger ring networks.