Tuning bat LSO neurons to interaural intensity differences through spike-timing dependent plasticity

Bats, like other mammals, are known to use interaural intensity differences (IID) to determine azimuthal position. In the lateral superior olive (LSO) neurons have firing behaviors which vary systematically with IID. Those neurons receive excitatory inputs from the ipsilateral ear and inhibitory inputs from the contralateral one. The IID sensitivity of a LSO neuron is thought to be due to delay differences between the signals coming from both ears, differences due to different synaptic delays and to intensity-dependent delays. In this paper we model the auditory pathway until the LSO. We propose a learning scheme where inputs to LSO neurons start out numerous with different relative delays. Spike timing-dependent plasticity (STDP) is then used to prune those connections. We compare the pruned neuron responses with physiological data and analyse the relationship between IID’s of teacher stimuli and IID sensitivities of trained LSO neurons.     

Vibrational resonance in a randomly connected neural network

A randomly connected network is constructed with similar characteristics (e.g., the ratio of excitatory and inhibitory neurons, the connection probability between neurons, and the axonal conduction delays) as that in the mammalian neocortex and the effects of high-frequency electrical field on the response of the network to a subthreshold low-frequency electrical field are studied in detail. It is found that both the amplitude and frequency of the high-frequency electrical field can modulate the response of the network to the low-frequency electric field. Moreover, vibrational resonance (VR) phenomenon induced by the two types of electrical fields can also be influenced by the network parameters, such as the neuron population, the connection probability between neurons and the synaptic strength. It is interesting that VR is found to be related with the ratio of excitatory neurons that are under high-frequency electrical stimuli. In summary, it is suggested that the interaction of excitatory and inhibitory currents is also an important factor that can influence the performance of VR in neural networks.     

Oscillation,Conduction Delays,and Learning Cooperate to Establish Neural Competition in Recurrent Networks

Specific memory might be stored in a subnetwork consisting of a small population of neurons. To select neurons involved in memory formation, neural competition might be essential. In this paper, we show that excitable neurons are competitive and organize into two assemblies in a recurrent network with spike timing-dependent synaptic plasticity (STDP) and axonal conduction delays. Neural competition is established by the cooperation of spontaneously induced neural oscillation, axonal conduction delays, and STDP. We also suggest that the competition mechanism in this paper is one of the basic functions required to organize memory-storing subnetworks into fine-scale cortical networks.     

Analysis of stability of neural network with inhibitory neurons

Phase coding in a neural network composed of neural oscillators with inhibitory neurons was studied based on the theory of stochastic phase dynamics. We found that with increasing the coupling coefficients of inhibitory neural oscillators, the firing density in excitatory population transits to a critical state. In this case, when we increase the inhibitory coupling, the firing density will come into dynamic balance again and tend to a fixed value gradually. According to the phenomenon, in the paper we found parameter regions to exhibit those different population states, called dividing zones including flat fading zone, rapid fading zone and critical zone. Based on the dividing zones we can choose the number ratio between inhibitory neurons and excitatory neurons in the neural network, and estimate the coupling action of inhibitory population and excitatory population. Our research also shows that the balance value, enabling the firing density to reach the dynamic balance, does not depend on initial conditions. In addition, the critical value in critical state is only related to the number ratio between inhibitory neurons and excitatory neurons, but is independent of inhibitory coupling and excitatory coupling.     

Effects of phase on homeostatic spike rates

Recent experimental results by Talathi et al. (Neurosci Lett 455:145–149, 2009) showed a divergence in the spike rates of two types of population spike events, representing the putative activity of the excitatory and inhibitory neurons in the CA1 area of an animal model for temporal lobe epilepsy. The divergence in the spike rate was accompanied by a shift in the phase of oscillations between these spike rates leading to a spontaneous epileptic seizure. In this study, we propose a model of homeostatic synaptic plasticity which assumes that the target spike rate of populations of excitatory and inhibitory neurons in the brain is a function of the phase difference between the excitatory and inhibitory spike rates. With this model of homeostatic synaptic plasticity, we are able to simulate the spike rate dynamics seen experimentally by Talathi et al. in a large network of interacting excitatory and inhibitory neurons using two different spiking neuron models. A drift analysis of the spike rates resulting from the homeostatic synaptic plasticity update rule allowed us to determine the type of synapse that may be primarily involved in the spike rate imbalance in the experimental observation by Talathi et al. We find excitatory neurons, particularly those in which the excitatory neuron is presynaptic, have the most influence in producing the diverging spike rates and causing the spike rates to be anti-phase. Our analysis suggests that the excitatory neuronal population, more specifically the excitatory to excitatory synaptic connections, could be implicated in a methodology designed to control epileptic seizures.     

Weak electric fields detectability in a noisy neural network

We investigate the detectability of weak electric field in a noisy neural network based on Izhikevich neuron model systematically. The neural network is composed of excitatory and inhibitory neurons with similar ratio as that in the mammalian neocortex, and the axonal conduction delays between neurons are also considered. It is found that the noise intensity can modulate the detectability of weak electric field. Stochastic resonance (SR) phenomenon induced by white noise is observed when the weak electric field is added to the network. It is interesting that SR almost disappeared when the connections between neurons are cancelled, suggesting the amplification effects of the neural coupling on the synchronization of neuronal spiking. Furthermore, the network parameters, such as the connection probability, the synaptic coupling strength, the scale of neuron population and the neuron heterogeneity, can also affect the detectability of the weak electric field. Finally, the model sensitivity is studied in detail, and results show that the neural network model has an optimal region for the detectability of weak electric field signal.     

Longitudinal neuronal organization and coordination in a simple vertebrate: a continuous,semi-quantitative computer model of the central pattern generator for swimming in young frog tadpoles

When frog tadpoles hatch their swimming requires co-ordinated contractions of trunk muscles, driven by motoneurons and controlled by a Central Pattern Generator (CPG). To study this co-ordination we used a 3.5 mm long population model of the young tadpole CPG with continuous distributions of neurons and axon lengths as estimated anatomically. We found that: (1) alternating swimming-type activity fails to self-sustain unless some excitatory interneurons have ascending axons, (2) a rostro-caudal (R-C) gradient in the distribution of excitatory premotor interneurons with short axons is required to obtain the R-C gradient in excitation and resulting progression of motoneuron firing necessary for forward swimming, (3) R-C delays in motoneuron firing decrease if excitatory motoneuron to premotor interneuron synapses are present, (4) these feedback connections and the electrical synapses between motoneurons synchronise motoneuron discharges locally, (5) the above findings are independent of the detailed membrane properties of neurons.     

Impact of adaptation currents on synchronization of coupled exponential integrate-and-fire neurons

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies.     

Short conduction delays cause inhibition rather than excitation to favor synchrony in hybrid neuronal networks of the entorhinal cortex

How stable synchrony in neuronal networks is sustained in the presence of conduction delays is an open question. The Dynamic Clamp was used to measure phase resetting curves (PRCs) for entorhinal cortical cells, and then to construct networks of two such neurons. PRCs were in general Type I (all advances or all delays) or weakly type II with a small region at early phases with the opposite type of resetting. We used previously developed theoretical methods based on PRCs under the assumption of pulsatile coupling to predict the delays that synchronize these hybrid circuits. For excitatory coupling, synchrony was predicted and observed only with no delay and for delays greater than half a network period that cause each neuron to receive an input late in its firing cycle and almost immediately fire an action potential. Synchronization for these long delays was surprisingly tight and robust to the noise and heterogeneity inherent in a biological system. In contrast to excitatory coupling, inhibitory coupling led to antiphase for no delay, very short delays and delays close to a network period, but to near-synchrony for a wide range of relatively short delays. PRC-based methods show that conduction delays can stabilize synchrony in several ways, including neutralizing a discontinuity introduced by strong inhibition, favoring synchrony in the case of noisy bistability, and avoiding an initial destabilizing region of a weakly type II PRC. PRCs can identify optimal conduction delays favoring synchronization at a given frequency, and also predict robustness to noise and heterogeneity.     

Emerging methodologies for the study of hypothalamic gonadotropin-releasing-hormone (GnRH) neurons

Gonadotropin-releasing-hormone (GnRH) neurons form part of a central neural oscillator that controls sexual reproduction through intermittent release of the GnRH peptide. Activity of GnRH neurons, and by extension release of GnRH, has been proposed to reflect intrinsic properties and synaptic input of GnRH neurons. To study GnRH neurons, we used traditional electrophysiology and computational methods. These emerging methodologies enhance the elucidation of processing in GnRH neurons. We used dynamic current-clamping to understand how living GnRH somata process input from glutamate and GABA, two key neurotransmitters in the neuroendocrine hypothalamus. In order to study the impact of synaptic integration in dendrites and neuronal morphology, we have developed full-morphology models of GnRH neurons. Using dynamic clamping, we have demonstrated that small-amplitude glutamatergic currents can drive repetitive firing in GnRH neurons. Furthermore, application of simulated GABAergic synapses with a depolarized reversal potential have revealed two functional subpopulations of GnRH neurons: one population in which GABA chronically depolarizes membrane potential (without inducing action potentials) and a second population in which GABAergic excitation results in slow spiking. Finally, when AMPA-type and GABA-type simulated inputs are applied together, action potentials occur when the AMPA-type conductance occurs during the descending phase of GABAergic excitation and at the nadir of GABAergic inhibition. Compartmental computer models have shown that excitatory synapses at >300 microns from somtata are unable to drive spiking with purely passive dendrites. In models with active dendrites, distal synapses are more efficient at driving spiking than somatic inputs. We then used our models to extend the results from dynamic current clamping at GnRH somata to distribute synaptic inputs along the dendrite. We show that propagation delays for dendritic synapses alter synaptic integration in GnRH neurons by widening the temporal window of interaction for the generation of action potentials. Finally, we have shown that changes in dendrite morphology can modulate the output of GnRH neurons by altering the efficacy of action potential generation in response to after-depolarization potentials (ADPs). Taken together, the methodologies of dynamic current clamping and multi-compartmental modeling can make major contributions to the study of synaptic integration and structure-function relationships in hypothalamic GnRH neurons. Use of these methodological approaches will continue to provide keen insights leading to conceptual advances in our understanding of reproductive hormone secretion in normal and pathological physiology and open the door to understanding whether the mechanisms of pulsatile GnRH release are conserved across species.     

Conductance-based refractory density model of primary visual cortex

A layered continual population model of primary visual cortex has been constructed, which reproduces a set of experimental data, including postsynaptic responses of single neurons on extracellular electric stimulation and spatially distributed activity patterns in response to visual stimulation. In the model, synaptically interacting excitatory and inhibitory neuronal populations are described by a conductance-based refractory density approach. Populations of two-compartment excitatory and inhibitory neurons in cortical layers 2/3 and 4 are distributed in the 2-d cortical space and connected by AMPA, NMDA and GABA type synapses. The external connections are pinwheel-like, according to the orientation of a stimulus. Intracortical connections are isotropic local and patchy between neurons with similar orientations. The model proposes better temporal resolution and more detailed elaboration than conventional mean-field models. In comparison to large network simulations, it excludes a posteriori statistical data manipulation and provides better computational efficiency and minimal parametrization.     

A spiking neuron model for synchronous flashing of fireflies

Certain species of fireflies show a group behavior of synchronous flashing. Their synchronized and rhythmic flashing has received much attention among many researchers, and there has been a study of biological models for their entrainment of flashing. The synchronous behavior of fireflies resembles the firing synchrony of integrate-and-fire neurons with excitatory or inhibitory connections. This paper shows an analysis of spiking neurons specialized for a firefly flashing model, and provides simulation results of multiple neurons with various transmission delays and coupling strengths. It also explains flashing patterns of some firefly species and examines the synchrony conditions depending on transmission delays and coupling strengths.     

Dynamic Excitatory and Inhibitory Gain Modulation Can Produce Flexible,Robust and Optimal Decision-making

Behavioural and neurophysiological studies in primates have increasingly shown the involvement of urgency signals during the temporal integration of sensory evidence in perceptual decision-making. Neuronal correlates of such signals have been found in the parietal cortex, and in separate studies, demonstrated attention-induced gain modulation of both excitatory and inhibitory neurons. Although previous computational models of decision-making have incorporated gain modulation, their abstract forms do not permit an understanding of the contribution of inhibitory gain modulation. Thus, the effects of co-modulating both excitatory and inhibitory neuronal gains on decision-making dynamics and behavioural performance remain unclear. In this work, we incorporate time-dependent co-modulation of the gains of both excitatory and inhibitory neurons into our previous biologically based decision circuit model. We base our computational study in the context of two classic motion-discrimination tasks performed in animals. Our model shows that by simultaneously increasing the gains of both excitatory and inhibitory neurons, a variety of the observed dynamic neuronal firing activities can be replicated. In particular, the model can exhibit winner-take-all decision-making behaviour with higher firing rates and within a significantly more robust model parameter range. It also exhibits short-tailed reaction time distributions even when operating near a dynamical bifurcation point. The model further shows that neuronal gain modulation can compensate for weaker recurrent excitation in a decision neural circuit, and support decision formation and storage. Higher neuronal gain is also suggested in the more cognitively demanding reaction time than in the fixed delay version of the task. Using the exact temporal delays from the animal experiments, fast recruitment of gain co-modulation is shown to maximize reward rate, with a timescale that is surprisingly near the experimentally fitted value. Our work provides insights into the simultaneous and rapid modulation of excitatory and inhibitory neuronal gains, which enables flexible, robust, and optimal decision-making.     

Population rate coding in recurrent neuronal networks with unreliable synapses

Neuron transmits spikes to postsynaptic neurons through synapses. Experimental observations indicated that the communication between neurons is unreliable. However most modelling and computational studies considered deterministic synaptic interaction model. In this paper, we investigate the population rate coding in an all-to-all coupled recurrent neuronal network consisting of both excitatory and inhibitory neurons connected with unreliable synapses. We use a stochastic on-off process to model the unreliable synaptic transmission. We find that synapses with suitable successful transmission probability can enhance the encoding performance in the case of weak noise; while in the case of strong noise, the synaptic interactions reduce the encoding performance. We also show that several important synaptic parameters, such as the excitatory synaptic strength, the relative strength of inhibitory and excitatory synapses, as well as the synaptic time constant, have significant effects on the performance of the population rate coding. Further simulations indicate that the encoding dynamics of our considered network cannot be simply determined by the average amount of received neurotransmitter for each neuron in a time instant. Moreover, we compare our results with those obtained in the corresponding random neuronal networks. Our numerical results demonstrate that the network randomness has the similar qualitative effect as the synaptic unreliability but not completely equivalent in quantity.     

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     

Effect of the topology and delayed interactions in neuronal networks synchronization

As important as the intrinsic properties of an individual nervous cell stands the network of neurons in which it is embedded and by virtue of which it acquires great part of its responsiveness and functionality. In this study we have explored how the topological properties and conduction delays of several classes of neural networks affect the capacity of their constituent cells to establish well-defined temporal relations among firing of their action potentials. This ability of a population of neurons to produce and maintain a millisecond-precise coordinated firing (either evoked by external stimuli or internally generated) is central to neural codes exploiting precise spike timing for the representation and communication of information. Our results, based on extensive simulations of conductance-based type of neurons in an oscillatory regime, indicate that only certain topologies of networks allow for a coordinated firing at a local and long-range scale simultaneously. Besides network architecture, axonal conduction delays are also observed to be another important factor in the generation of coherent spiking. We report that such communication latencies not only set the phase difference between the oscillatory activity of remote neural populations but determine whether the interconnected cells can set in any coherent firing at all. In this context, we have also investigated how the balance between the network synchronizing effects and the dispersive drift caused by inhomogeneities in natural firing frequencies across neurons is resolved. Finally, we show that the observed roles of conduction delays and frequency dispersion are not particular to canonical networks but experimentally measured anatomical networks such as the macaque cortical network can display the same type of behavior.     

Computational Modeling of Seizure Dynamics Using Coupled Neuronal Networks: Factors Shaping Epileptiform Activity

Epileptic seizure dynamics span multiple scales in space and time. Understanding seizure mechanisms requires identifying the relations between seizure components within and across these scales, together with the analysis of their dynamical repertoire. Mathematical models have been developed to reproduce seizure dynamics across scales ranging from the single neuron to the neural population. In this study, we develop a network model of spiking neurons and systematically investigate the conditions, under which the network displays the emergent dynamic behaviors known from the Epileptor, which is a well-investigated abstract model of epileptic neural activity. This approach allows us to study the biophysical parameters and variables leading to epileptiform discharges at cellular and network levels. Our network model is composed of two neuronal populations, characterized by fast excitatory bursting neurons and regular spiking inhibitory neurons, embedded in a common extracellular environment represented by a slow variable. By systematically analyzing the parameter landscape offered by the simulation framework, we reproduce typical sequences of neural activity observed during status epilepticus. We find that exogenous fluctuations from extracellular environment and electro-tonic couplings play a major role in the progression of the seizure, which supports previous studies and further validates our model. We also investigate the influence of chemical synaptic coupling in the generation of spontaneous seizure-like events. Our results argue towards a temporal shift of typical spike waves with fast discharges as synaptic strengths are varied. We demonstrate that spike waves, including interictal spikes, are generated primarily by inhibitory neurons, whereas fast discharges during the wave part are due to excitatory neurons. Simulated traces are compared with in vivo experimental data from rodents at different stages of the disorder. We draw the conclusion that slow variations of global excitability, due to exogenous fluctuations from extracellular environment, and gap junction communication push the system into paroxysmal regimes. We discuss potential mechanisms underlying such machinery and the relevance of our approach, supporting previous detailed modeling studies and reflecting on the limitations of our methodology.     

Cellularly-driven differences in network synchronization propensity are differentially modulated by firing frequency

Spatiotemporal pattern formation in neuronal networks depends on the interplay between cellular and network synchronization properties. The neuronal phase response curve (PRC) is an experimentally obtainable measure that characterizes the cellular response to small perturbations, and can serve as an indicator of cellular propensity for synchronization. Two broad classes of PRCs have been identified for neurons: Type I, in which small excitatory perturbations induce only advances in firing, and Type II, in which small excitatory perturbations can induce both advances and delays in firing. Interestingly, neuronal PRCs are usually attenuated with increased spiking frequency, and Type II PRCs typically exhibit a greater attenuation of the phase delay region than of the phase advance region. We found that this phenomenon arises from an interplay between the time constants of active ionic currents and the interspike interval. As a result, excitatory networks consisting of neurons with Type I PRCs responded very differently to frequency modulation compared to excitatory networks composed of neurons with Type II PRCs. Specifically, increased frequency induced a sharp decrease in synchrony of networks of Type II neurons, while frequency increases only minimally affected synchrony in networks of Type I neurons. These results are demonstrated in networks in which both types of neurons were modeled generically with the Morris-Lecar model, as well as in networks consisting of Hodgkin-Huxley-based model cortical pyramidal cells in which simulated effects of acetylcholine changed PRC type. These results are robust to different network structures, synaptic strengths and modes of driving neuronal activity, and they indicate that Type I and Type II excitatory networks may display two distinct modes of processing information.     

The influence of depolarization block on seizure-like activity in networks of excitatory and inhibitory neurons

The inhibitory restraint necessary to suppress aberrant activity can fail when inhibitory neurons cease to generate action potentials as they enter depolarization block. We investigate possible bifurcation structures that arise at the onset of seizure-like activity resulting from depolarization block in inhibitory neurons. Networks of conductance-based excitatory and inhibitory neurons are simulated to characterize different types of transitions to the seizure state, and a mean field model is developed to verify the generality of the observed phenomena of excitatory-inhibitory dynamics. Specifically, the inhibitory population’s activation function in the Wilson-Cowan model is modified to be non-monotonic to reflect that inhibitory neurons enter depolarization block given strong input. We find that a physiological state and a seizure state can coexist, where the seizure state is characterized by high excitatory and low inhibitory firing rate. Bifurcation analysis of the mean field model reveals that a transition to the seizure state may occur via a saddle-node bifurcation or a homoclinic bifurcation. We explain the hysteresis observed in network simulations using these two bifurcation types. We also demonstrate that extracellular potassium concentration affects the depolarization block threshold; the consequent changes in bifurcation structure enable the network to produce the tonic to clonic phase transition observed in biological epileptic networks.