A Computational Model of Neuro-Glio-Vascular Loop Interactions

We present a computational, biophysical model of neuron-astrocyte-vessel interaction. Unlike other cells, neurons convey “hunger” signals to the vascular network via an intervening layer of glial cells (astrocytes); vessels dilate and release glucose which fuels neuronal firing. Existing computational models focus on only parts of this loop (neuron→astrocyte→vessel→neuron), whereas the proposed model describes the entire loop. Neuronal firing causes release of a neurotransmitter like glutamate which triggers release of vasodilator by astrocytes via a cascade of biochemical events. Vasodilators released from astrocytic endfeet cause blood vessels to dilate and release glucose into the interstitium, part of which is taken up by the astrocyticendfeet. Glucose is converted into lactate in the astrocyte and transported into the neuron. Glucose from the interstitium and lactate (produced from glucose) influx from astrocyte are converted into ATP in the neuron. Neuronal ATP is used to drive the Na⁺/K⁺ATPase pumps, which maintain ionic gradients necessary for neuronal firing. When placed in the metabolic loop, the neuron exhibits sustained firing only when the stimulation current is more than a minimum threshold. For various combinations of initial neuronal [ATP] and external current, the neuron exhibits a variety of firing patterns including sustained firing, firing after an initial pause, burst firing etc. Neurovascular interactions under conditions of constricted vessels are also studied. Most models of cerebral circulation describe neurovascular interactions exclusively in the “forward” neuron→vessel direction. The proposed model indicates possibility of “reverse” influence also, with vasomotion rhythms influencing neural firing patterns. Another idea that emerges out of the proposed work is that brain''s computations may be more comprehensively understood in terms of neuro-glial-vascular dynamics and not in terms of neural firing alone.     

Neuronal Chemokines: Versatile Messengers In Central Nervous System Cell Interaction

Whereas chemokines are well known for their ability to induce cell migration, only recently it became evident that chemokines also control a variety of other cell functions and are versatile messengers in the interaction between a diversity of cell types. In the central nervous system (CNS), chemokines are generally found under both physiological and pathological conditions. Whereas many reports describe chemokine expression in astrocytes and microglia and their role in the migration of leukocytes into the CNS, only few studies describe chemokine expression in neurons. Nevertheless, the expression of neuronal chemokines and the corresponding chemokine receptors in CNS cells under physiological and pathological conditions indicates that neuronal chemokines contribute to CNS cell interaction. In this study, we review recent studies describing neuronal chemokine expression and discuss potential roles of neuronal chemokines in neuron–astrocyte, neuron–microglia, and neuron–neuron interaction.     

Astrocyte- neuron interaction as a mechanism responsible for generation of neural synchrony: a study based on modeling and experiments

Neural synchronization is considered as an important mechanism for information processing. In addition, based on recent neurophysiologic findings, it is believed that astrocytes regulate the synaptic transmission of neuronal networks. Therefore, the present study focused on determining the functional contribution of astrocytes in neuronal synchrony using both computer simulations and extracellular field potential recordings. For computer simulations, as a first step, a minimal network model is constructed by connecting two Morris-Lecar neuronal models. In this minimal model, astrocyte-neuron interactions are considered in a functional-based procedure. Next, the minimal network is extended and a biologically plausible neuronal population model is developed which considers functional outcome of astrocyte-neuron interactions too. The employed structure is based on the physiological and anatomical network properties of the hippocampal CA1 area. Utilizing these two different levels of modeling, it is demonstrated that astrocytes are able to change the threshold value of transition from synchronous to asynchronous behavior among neurons. In this way, variations in the interaction between astrocytes and neurons lead to the emergence of synchronous/asynchronous patterns in neural responses. Furthermore, population spikes are recorded from CA1 pyramidal neurons in rat hippocampal slices to validate the modeling results. It demonstrates that astrocytes play a primary role in neuronal firing synchronicity and synaptic coordination. These results may offer a new insight into understanding the mechanism by which astrocytes contribute to stabilizing neural activities.     

Modulation of synchrony without changes in firing rates

It was often reported and suggested that the synchronization of spikes can occur without changes in the firing rate. However, few theoretical studies have tested its mechanistic validity. In the present study, we investigate whether changes in synaptic weights can induce an independent modulation of synchrony while the firing rate remains constant. We study this question at the level of both single neurons and neuronal populations using network simulations of conductance based integrate-and-fire neurons. The network consists of a single layer that includes local excitatory and inhibitory recurrent connections, as well as long-range excitatory projections targeting both classes of neurons. Each neuron in the network receives external input consisting of uncorrelated Poisson spike trains. We find that increasing this external input leads to a linear increase of activity in the network, as well␣as an increase in the peak frequency of oscillation. In␣contrast, balanced changes of the synaptic weight of␣excitatory long-range projections for both classes of postsynaptic neurons modulate the degree of synchronization without altering the firing rate. These results demonstrate that, in a simple network, synchronization and firing rate can be modulated independently, and thus, may be used as independent coding dimensions. Electronic supplementary material  The online version of this article (doi: ) contains supplementary material, which is available to authorized users.     

Learning-Induced Synchronization and Plasticity of a Developing Neural Network

Learning-induced synchronization of a neural network at various developing stages is studied by computer simulations using a pulse-coupled neural network model in which the neuronal activity is simulated by a one-dimensional map. Two types of Hebbian plasticity rules are investigated and their differences are compared. For both models, our simulations show a logarithmic increase in the synchronous firing frequency of the network with the culturing time of the neural network. This result is consistent with recent experimental observations. To investigate how to control the synchronization behavior of a neural network after learning, we compare the occurrence of synchronization for four networks with different designed patterns under the influence of an external signal. The effect of such a signal on the network activity highly depends on the number of connections between neurons. We discuss the synaptic plasticity and enhancement effects for a random network after learning at various developing stages.     

Neural Synchronization at Tonic-to-Bursting Transitions

We studied the synchronous behavior of two electrically-coupled model neurons as a function of the coupling strength when the individual neurons are tuned to different activity patterns that ranged from tonic firing via chaotic activity to burst discharges. We observe asynchronous and various synchronous states such as out-of-phase, in-phase and almost in-phase chaotic synchronization. The highest variety of synchronous states occurs at the transition from tonic firing to chaos where the highest coupling strength is also needed for in-phase synchronization which is, essentially, facilitated towards the bursting range. This demonstrates that tuning of the neuron’s internal dynamics can have significant impact on the synchronous states especially at the physiologically relevant tonic-to-bursting transitions.     

A Markov model for interspike interval distributions of auditory cortical neurons that do not show periodic firings

Spontaneous firing properties of individual auditory cortical neurons are interpreted in terms of local and global order present in functioning brain networks, such as alternating “up” and “down” states. A four-state modulated Markov process is used to model neuronal firings. The system alternates between a bound and an unbound state, both with Poisson-distributed lifetimes. During the unbound state, active and closed states alternate with Poisson-distributed lifetimes. Inside the active state, spikes are generated as a realization of a Poisson process. This combination of processes constitutes a four-state modulated Markov process, determined by five independent parameters. Analytical expressions for the probability density functions (pdfs) that describe the interspike interval (ISI) distribution and autocorrelation function are derived. The pdf for the ISI distribution is shown to be a linear combination of three exponential functions and is expressed through the five system parameters. Through fitting experimental ISI histograms by the theoretical ones, numerical values of the system parameters are obtained for the individual neurons. Both Monte Carlo simulations and goodness-of-fit tests are used to validate the fitting procedure. The values of the estimated system parameters related to the active-closed and bound–unbound processes and their independence on the neurons’ mean firing rate suggest that the underlying quasi-periodic processes reflect properties of the network in which the neurons are embedded. The characteristic times of autocorrelations, determined by the bound–unbound and active-closed processes, are also independent of the neuron’s firing rate. The agreement between experimental and theoretical ISI histograms and autocorrelation functions allows interpretation of the system parameters of the individual neurons in terms of slow and delta waves, and high-frequency oscillations observed in cortical networks. This procedure can identify and track the influence of changing brain states on the single-unit firing patterns in experimental animals.     

Robust emergence of small-world structure in networks of spiking neurons

Spontaneous activity in biological neural networks shows patterns of dynamic synchronization. We propose that these patterns support the formation␣of a small-world structure—network connectivity␣optimal for distributed information processing. We␣present numerical simulations with connected Hindmarsh–Rose neurons in which, starting from random connection distributions, small-world networks evolve as a result of applying an adaptive rewiring rule. The rule connects pairs of neurons that tend fire in synchrony, and disconnects ones that fail to synchronize. Repeated application of the rule leads to small-world structures. This mechanism is robustly observed for bursting and irregular firing regimes.     

Translation-invariant pattern recognition based on Synfire chains

Most of current neural network architectures are not suited to recognize a pattern at various displaced positions. This lack seems due to the prevailing neuron model which reduces a neuron's information transmission to its firing rate. With this information code, a neuronal assembly cannot distinguish between different combinations of its entities and therefore fails to represent the fine structure within a pattern. In our approach, the main idea of the correlation theory is accepted that spatial relationships in a pattern should be coded by temporal relations in the timing of action potentials. However, we do not assume that synchronized spikes are a sign for strong synapses between the neurons concerned. Instead, the synchronization of Synfire chains can be exploited to produce the relevant timing relationships between the neuronal signals. Therefore, we do not require fast synaptic plasticity to account for the precise timing of action potentials. In order to illustrate this claim, we propose a model for translation-invariant pattern recognition which does not depend on any changes in synaptic efficacies. Received: 14 June 1998 / Accepted in revised form: 9 January 1999     

Synchronization study in ring-like and grid-like neuronal networks

In this paper, we study the synchronization status of both two gap-junction coupled neurons and neuronal network with two different network connectivity patterns. One of the network connectivity patterns is a ring-like neuronal network, which only considers nearest-neighbor neurons. The other is a grid-like neuronal network, with all nearest neighbor couplings. We show that by varying some key parameters, such as the coupling strength and the external current injection, the neuronal network will exhibit various patterns of firing synchronization. Different types of firing synchronization are diagnosed by means of a mean field potential, a bifurcation diagram, a correlation coefficient and the ISI-distance method. Numerical simulations demonstrate that the synchronization status of multiple neurons is much dependent on the network patters, when the number of neurons is the same. It is also demonstrated that the synchronization status of two coupled neurons is similar with the grid-like neuronal network, but differs radically from that of the ring-like neuronal network. These results may be instructive in understanding synchronization transitions in neuronal systems.     

Influences of membrane properties on phase response curve and synchronization stability in a model globus pallidus neuron

The activity patterns of the globus pallidus (GPe) and subthalamic nucleus (STN) are closely associated with motor function and dysfunction in the basal ganglia. In the pathological state caused by dopamine depletion, the STN–GPe network exhibits rhythmic synchronous activity accompanied by rebound bursts in the STN. Therefore, the mechanism of activity transition is a key to understand basal ganglia functions. As synchronization in GPe neurons could induce pathological STN rebound bursts, it is important to study how synchrony is generated in the GPe. To clarify this issue, we applied the phase-reduction technique to a conductance-based GPe neuronal model in order to derive the phase response curve (PRC) and interaction function between coupled GPe neurons. Using the PRC and interaction function, we studied how the steady-state activity of the GPe network depends on intrinsic membrane properties, varying ionic conductances on the membrane. We noted that a change in persistent sodium current, fast delayed rectifier Kv3 potassium current, M-type potassium current and small conductance calcium-dependent potassium current influenced the PRC shape and the steady state. The effect of those currents on the PRC shape could be attributed to extension of the firing period and reduction of the phase response immediately after an action potential. In particular, the slow potassium current arising from the M-type potassium and the SK current was responsible for the reduction of the phase response. These results suggest that the membrane property modulation controls synchronization/asynchronization in the GPe and the pathological pattern of STN–GPe activity.     

Firing synchronization and temporal order in noisy neuronal networks

Noise-induced complete synchronization and frequency synchronization in coupled spiking and bursting neurons are studied firstly. The effects of noise and coupling are discussed. It is found that bursting neurons are easier to achieve firing synchronization than spiking ones, which means that bursting activities are more important for information transfer in neuronal networks. Secondly, the effects of noise on firing synchronization in a noisy map neuronal network are presented. Noise-induced synchronization and temporal order are investigated by means of the firing rate function and the order index. Firing synchronization and temporal order of excitatory neurons can be greatly enhanced by subthreshold stimuli with resonance frequency. Finally, it is concluded that random perturbations play an important role in firing activities and temporal order in neuronal networks.     

Leader neurons in leaky integrate and fire neural network simulations

In this paper, we highlight the topological properties of leader neurons whose existence is an experimental fact. Several experimental studies show the existence of leader neurons in population bursts of activity in 2D living neural networks (Eytan and Marom, J Neurosci 26(33):8465–8476, 2006; Eckmann et al., New J Phys 10(015011), 2008). A leader neuron is defined as a neuron which fires at the beginning of a burst (respectively network spike) more often than we expect by chance considering its mean firing rate. This means that leader neurons have some burst triggering power beyond a chance-level statistical effect. In this study, we characterize these leader neuron properties. This naturally leads us to simulate neural 2D networks. To build our simulations, we choose the leaky integrate and fire (lIF) neuron model (Gerstner and Kistler 2002; Cessac, J Math Biol 56(3):311–345, 2008), which allows fast simulations (Izhikevich, IEEE Trans Neural Netw 15(5):1063–1070, 2004; Gerstner and Naud, Science 326:379–380, 2009). The dynamics of our lIF model has got stable leader neurons in the burst population that we simulate. These leader neurons are excitatory neurons and have a low membrane potential firing threshold. Except for these two first properties, the conditions required for a neuron to be a leader neuron are difficult to identify and seem to depend on several parameters involved in the simulations themselves. However, a detailed linear analysis shows a trend of the properties required for a neuron to be a leader neuron. Our main finding is: A leader neuron sends signals to many excitatory neurons as well as to few inhibitory neurons and a leader neuron receives only signals from few other excitatory neurons. Our linear analysis exhibits five essential properties of leader neurons each with different relative importance. This means that considering a given neural network with a fixed mean number of connections per neuron, our analysis gives us a way of predicting which neuron is a good leader neuron and which is not. Our prediction formula correctly assesses leadership for at least ninety percent of neurons.     

A biologically motivated and analytically soluble model of collective oscillations in the cortex

. Feature linking and pattern separation are shown to be performed as simultaneous processes by a highly connected auto-associative network of spiking neurons (spike response model). In principle, many (e.g., with nine) patterns can be separated, but with a biological set of parameters the number is limited to four. The patterns have been learned by an asymmetric hebbian rule that can handle a low activity which may vary from pattern to pattern (in a range between 4% and 7%). Spikes are generated by a threshold process and – with some delay – transmitted to postsynaptic neurons. There they evoke an excitatory or inhibitory postsynaptic potential (EPSP or IPSP). Spike emission is followed by an absolute refractory period (1 ms) and activates an inhibitory delay loop that prevents continuous firing. Three different network topologies are discussed, i.e., a structureless fully connected system, a network composed of two ‘hemispheres’, and finally a hierarchical network with four subsystems that represent different ‘functions’ and interact via feedforward and feedback connections. Functional feedback turns out to be essential for context-sensitive binding. The coherence between the two hemispheres is dependent on the interhemispheric delays. If these are on average too large, the two hemispheres oscillate coherently by themselves but phase-shifted by half a period with respect to each other. Received: 16 June 1993/Accepted in revised form: 24 March 1994     

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     

Cortical network modeling: analytical methods for firing rates and some properties of networks of LIF neurons.

The circuitry of cortical networks involves interacting populations of excitatory (E) and inhibitory (I) neurons whose relationships are now known to a large extent. Inputs to E- and I-cells may have their origins in remote or local cortical areas. We consider a rudimentary model involving E- and I-cells. One of our goals is to test an analytic approach to finding firing rates in neural networks without using a diffusion approximation and to this end we consider in detail networks of excitatory neurons with leaky integrate-and-fire (LIF) dynamics. A simple measure of synchronization, denoted by S(q), where q is between 0 and 100 is introduced. Fully connected E-networks have a large tendency to become dominated by synchronously firing groups of cells, except when inputs are relatively weak. We observed random or asynchronous firing in such networks with diverse sets of parameter values. When such firing patterns were found, the analytical approach was often able to accurately predict average neuronal firing rates. We also considered several properties of E-E networks, distinguishing several kinds of firing pattern. Included were those with silences before or after periods of intense activity or with periodic synchronization. We investigated the occurrence of synchronized firing with respect to changes in the internal excitatory postsynaptic potential (EPSP) magnitude in a network of 100 neurons with fixed values of the remaining parameters. When the internal EPSP size was less than a certain value, synchronization was absent. The amount of synchronization then increased slowly as the EPSP amplitude increased until at a particular EPSP size the amount of synchronization abruptly increased, with S(5) attaining the maximum value of 100%. We also found network frequency transfer characteristics for various network sizes and found a linear dependence of firing frequency over wide ranges of the external afferent frequency, with non-linear effects at lower input frequencies. The theory may also be applied to sparsely connected networks, whose firing behaviour was found to change abruptly as the probability of a connection passed through a critical value. The analytical method was also found to be useful for a feed-forward excitatory network and a network of excitatory and inhibitory neurons.     

Neural network implementation of a three-phase model of respiratory rhythm generation

A mathematical model of the central neural mechanisms of respiratory rhythm generation is developed. This model assumes that the respiratory cycle consists of three phases: inspiration, post-inspiration, and expiration. Five respiratory neuronal groups are included: inspiratory, late-inspiratory, post-inspiratory, expiratory, and early-inspiratory neurons. Proposed interconnections among these groups are based substantially on previous physiological findings. The model produces a stable limit cycle and generally reproduces the features of the firing patterns of the 5 neuronal groups. When simulated feedback from pulmonary stretch receptors is made to excite late-inspiratory neurons and inhibit early-inspiratory neurons, the model quantitatively reproduces previous observations of the expiratory-prolonging effects of pulses and steps of vagal afferent activity presented in expiration. In addition the model reproduces expected respiratory cycle timing and amplitude responses to change of chemical drive both in the absence and in the presence of simulated stretch receptor feedback. These results demonstrate the feasibility of generating the respiratory rhythm with a simple neural network based on observed respiratory neuronal groups. Other neuronal groups not included in the model may be more important for shaping the waveforms than for generating the basic oscillation.     

Delayed reverberation through time windows as a key to cerebellar function

We present a functional model of the cerebellum comprising cerebellar cortex, inferior olive, deep cerebellar nuclei, and brain stem nuclei. The discerning feature of the model being time coding, we consistently describe the system in terms of postsynaptic potentials, synchronous action potentials, and propagation delays. We show by means of detailed single-neuron modeling that (i) Golgi cells can fulfill a gating task in that they form short and well-defined time windows within which granule cells can reach firing threshold, thus organizing neuronal activity in discrete `time slices', and that (ii) rebound firing in cerebellar nuclei cells is a robust mechanism leading to a delayed reverberation of Purkinje cell activity through cerebellar-reticular projections back to the cerebellar cortex. Computer simulations of the whole cerebellar network consisting of several thousand neurons reveal that reverberation in conjunction with long-term plasticity at the parallel fiber-Purkinje cell synapses enables the system to learn, store, and recall spatio-temporal patterns of neuronal activity. Climbing fiber spikes act both as a synchronization and as a teacher signal, not as an error signal. They are due to intrinsic oscillatory properties of inferior olivary neurons and to delayed reverberation within the network. In addition to clear experimental predictions the present theory sheds new light on a number of experimental observation such as the synchronicity of climbing fiber spikes and provides a novel explanation of how the cerebellum solves timing tasks on a time scale of several hundreds of milliseconds. Received: 23 July 1999 / Accepted in revised form: 31 August 1999     

Selective self-synchronization of pulses in neuronal networks of the visual system

In the present study we have checked the hypothesis that the degree of pulse synchrony in neuronal pools is determined by the level of excitation of the neuronal network or its loci and that this relationship does not depend on the factor that causes the excitation. Pulse reactions of neurons in pools (2–4 cells) of cat lateral geniculate nucleus and visual cortex were registered. Neurons were excited using either visual stimuli or glutamic acid microinjections into neuronal pools. The increase of neural pool excitation level (No) regardless of the type of stimulus was shown to increase the pulse synchrony (Ns), with a correlation coefficient of 0.716 ± 0.217. One may suppose that the level of neuronal network excitation “governs” the synchrony of pulses generated by the network, i.e., neuronal networks function in compliance with the principle of self-synchronization.     

The nature of the memory trace and its neurocomputational implications

The brain processes underlying cognitive tasks must be very robust. Disruptions such as the destruction of large numbers of neurons, or the impact of alcohol and lack of sleep do not have negative effects except when they occur in an extreme form. This robustness implies that the parameters determining the functioning of networks of individual neurons must have large ranges or there must exist stabilizing mechanisms that keep the functioning of a network within narrow bounds. The simulation of a minimal neuronal architecture necessary to study cognitive tasks is described, which consists of a loop of three cell-assemblies. A crucial factor in this architecture is the critical threshold of a cell-assembly. When activated at a level above the critical threshold, the activation in a cell-assembly is subject to autonomous growth, which leads to an oscillation in the loop. When activated below the critical threshold, excitation gradually extinguishes. In order to circumvent the large parameter space of spiking neurons, a rate-dependent model of neuronal firing was chosen. The resulting parameter space of 12 parameters was explored by means of a genetic algorithm. The ranges of the parameters for which the architecture produced the required oscillations and extinctions, turned out to be relatively narrow. These ranges remained narrow when a stabilizing mechanism, controlling the total amount of activation, was introduced. The architecture thus shows chaotic behaviour. Given the overall stability of the operation of the brain, it can be concluded that there must exist other mechanisms that make the network robust. Three candidate mechanisms are discussed: synaptic scaling, synaptic homeostasis, and the synchronization of neural spikes.