Minimal Conductance-Based Model of Auditory Coincidence Detector Neurons

Sound localization is a fundamental sensory function of a wide variety of animals. The interaural time difference (ITD), an important cue for sound localization, is computed in the auditory brainstem. In our previous modeling study, we introduced a two-compartment Hodgkin-Huxley type model to investigate how cellular and synaptic specializations may contribute to precise ITD computation of the barn owl''s auditory coincidence detector neuron. Although our model successfully reproduced fundamental physiological properties observed in vivo, it was unsuitable for mathematical analyses and large scale simulations because of a number of nonlinear variables. In the present study, we reduce our former model into three types of conductance-based integrate-and-fire (IF) models. We test their electrophysiological properties using data from published in vivo and in vitro studies. Their robustness to parameter changes and computational efficiencies are also examined. Our numerical results suggest that the single-compartment active IF model is superior to other reduced models in terms of physiological reproducibility and computational performance. This model will allow future theoretical studies that use more rigorous mathematical analysis and network simulations.     

Parameters for a model of an oscillating neuronal network in the cochlear nucleus defined by genetic algorithms

Chopper neurons in the cochlear nucleus are characterized by intrinsic oscillations with short average interspike intervals (ISIs) and relative level independence of their response (Pfeiffer, Exp Brain Res 1:220–235, 1966; Blackburn and Sachs, J Neurophysiol 62:1303–1329, 1989), properties which are unattained by models of single chopper neurons (e.g., Rothman and Manis, J Neurophysiol 89:3070–3082, 2003a). In order to achieve short ISIs, we optimized the time constants of Rothman and Manis single neuron model with genetic algorithms. Some parameters in the optimization, such as the temperature and the capacity of the cell, turned out to be crucial for the required acceleration of their response. In order to achieve the relative level independence, we have simulated an interconnected network consisting of Rothman and Manis neurons. The results indicate that by stabilization of intrinsic oscillations, it is possible to simulate the physiologically observed level independence of ISIs. As previously reviewed and demonstrated (Bahmer and Langner, Biol Cybern 95:371–379, 2006a), chopper neurons show a preference for ISIs which are multiples of 0.4 ms. It was also demonstrated that the network consisting of two optimized Rothman and Manis neurons which activate each other with synaptic delays of 0.4 ms shows a preference for ISIs of 0.8 ms. Oscillations with various multiples of 0.4 ms as ISIs may be derived from neurons in a more complex network that is activated by simultaneous input of an onset neuron and several auditory nerve fibers.     

Firing rate of a retinal neuron are not predictable from interspike interval statistics.

The intervals between successive action potentials (impulses, or "spikes") produced the maintained firing of a neuron (ISIs) are often treated as if they were independent on each other; that is, an impulse train is considered as a stationary renewal process. If this is so, the variability of the mean rate of firing impulses in a sequence of temporal windows should be predictable from the distribution of ISIs. This was found not to be the case for the maintained firing of retinal ganglion cells in goldfish. Although some evident nonstationarity sometimes resulted in greater variability of the observed rate distributions than those predicted (for relatively long temporal windows), as a general rule the observed rate distributions were considerable less dispersed than would be predicted by sampling of the ISI distributions. This was taken as evidence of long-term serial dependency between successive ISIs; however, two standard test for dependency (autocorrelations and serial correlograms failed to to reveal structure of sufficiently long duration to account for the effect noted.     

A mechanistic approach for modeling temperature-dependent consumer-resource dynamics

Paramount to our ability to manage and protect biological communities from impending changes in the environment is an understanding of how communities will respond. General mathematical models of community dynamics are often too simplistic to accurately describe this response, partly to retain mathematical tractability and partly for the lack of biologically pleasing functions representing the model/environment interface. We address these problems of tractability and plausibility in community/environment models by incorporating the Boltzmann factor (temperature dependence) in a bioenergetic consumer-resource framework. Our analysis leads to three predictions for the response of consumer-resource systems to increasing mean temperature (warming). First, mathematical extinctions do not occur with warming; however, stable systems may transition into an unstable (cycling) state. Second, there is a decrease in the biomass density of resources with warming. The biomass density of consumers may increase or decrease depending on their proximity to the feasibility (extinction) boundary. Third, consumer biomass density is more sensitive to warming than resource biomass density (with some exceptions). These predictions are in line with many current observations and experiments. The model presented and analyzed here provides an advancement in the testing framework for global change scenarios and hypotheses of latitudinal and elevational species distributions.     

Gamma oscillations as a mechanism for selective information transmission

In the past decades, many studies have focussed on the relation between the input and output of neurons with the aim to understand information processing by neurons. A particular aspect of neuronal information, which has not received much attention so far, concerns the problem of information transfer when a neuron or a population of neurons receives input from two or more (populations of) neurons, in particular when these (populations of) neurons carry different types of information. The aim of the present study is to investigate the responses of neurons to multiple inputs modulated in the gamma frequency range. By a combination of theoretical approaches and computer simulations, we test the hypothesis that enhanced modulation of synchronized excitatory neuronal activity in the gamma frequency range provides an advantage over a less synchronized input for various types of neurons. The results of this study show that the spike output of various types of neurons [i.e. the leaky integrate and fire neuron, the quadratic integrate and fire neuron and the Hodgkin–Huxley (HH) neuron] and that of excitatory–inhibitory coupled pairs of neurons, like the Pyramidal Interneuronal Network Gamma (PING) model, is highly phase-locked to the larger of two gamma-modulated input signals. This implies that the neuron selectively responds to the input with the larger gamma modulation if the amplitude of the gamma modulation exceeds that of the other signals by a certain amount. In that case, the output of the neuron is entrained by one of multiple inputs and that other inputs are not represented in the output. This mechanism for selective information transmission is enhanced for short membrane time constants of the neuron.     

Sensitive dependence of the coefficient of variation of interspike intervals on the lower boundary of membrane potential for the leaky integrate-and-fire neuron model

After the report of Softky and Koch [Softky, W.R., Koch, C., 1993. The highly irregular firing of cortical cells is inconsistent with temporal integration of random EPSPs. J. Neurosci. 13, 334-350], leaky integrate-and-fire models have been investigated to explain high coefficient of variation (CV) of interspike intervals (ISIs) at high firing rates observed in the cortex. The purpose of this paper is to study the effect of the position of a lower boundary of membrane potential on the possible value of CV of ISIs based on the diffusional leaky integrate-and-fire models with and without reversal potentials. Our result shows that the irregularity of ISIs for the diffusional leaky integrate-and-fire neuron significantly changes by imposing a lower boundary of membrane potential, which suggests the importance of the position of the lower boundary as well as that of the firing threshold when we study the statistical properties of leaky integrate-and-fire neuron models. It is worth pointing out that the mean-CV plot of ISIs for the diffusional leaky integrate-and-fire neuron with reversal potentials shows a close similarity to the experimental result obtained in Softky and Koch [Softky, W.R., Koch, C., 1993. The highly irregular firing of cortical cells is inconsistent with temporal integration of random EPSPs. J. Neurosci. 13, 334-350].     

Bifurcations of large networks of two-dimensional integrate and fire neurons

Recently, a class of two-dimensional integrate and fire models has been used to faithfully model spiking neurons. This class includes the Izhikevich model, the adaptive exponential integrate and fire model, and the quartic integrate and fire model. The bifurcation types for the individual neurons have been thoroughly analyzed by Touboul (SIAM J Appl Math 68(4):1045–1079, 2008). However, when the models are coupled together to form networks, the networks can display bifurcations that an uncoupled oscillator cannot. For example, the networks can transition from firing with a constant rate to burst firing. This paper introduces a technique to reduce a full network of this class of neurons to a mean field model, in the form of a system of switching ordinary differential equations. The reduction uses population density methods and a quasi-steady state approximation to arrive at the mean field system. Reduced models are derived for networks with different topologies and different model neurons with biologically derived parameters. The mean field equations are able to qualitatively and quantitatively describe the bifurcations that the full networks display. Extensions and higher order approximations are discussed.     

An integrate-and-fire model to generate spike trains with long-range dependence

Long-range dependence (LRD) has been observed in a variety of phenomena in nature, and for several years also in the spiking activity of neurons. Often, this is interpreted as originating from a non-Markovian system. Here we show that a purely Markovian integrate-and-fire (IF) model, with a noisy slow adaptation term, can generate interspike intervals (ISIs) that appear as having LRD. However a proper analysis shows that this is not the case asymptotically. For comparison, we also consider a new model of individual IF neuron with fractional (non-Markovian) noise. The correlations of its spike trains are studied and proven to have LRD, unlike classical IF models. On the other hand, to correctly measure long-range dependence, it is usually necessary to know if the data are stationary. Thus, a methodology to evaluate stationarity of the ISIs is presented and applied to the various IF models. We explain that Markovian IF models may seem to have LRD because of non-stationarities.     

On a response characteristic of a mathematical neuron model

A mathematical neuron model in the form of a nonlinear difference equation is proposed and its response characteristic is investigated.If a sequence of pulses with a fixed frequency is applied to the neuron model as an input, and the amplitude of the input pulses is progressively decreased, the firing frequency of the neuron model, regarded as the output, also decreases. The relationship between them is quite complicated, but a mathematical investigation reveals that it takes the form of an extended Cantor's function. This result explains the unusual and unsuspected phenomenon which was found by L. D. Harmon in experimental studies with his transistor neuron models.Besides this, as an analogue of our mathematical neuron model, a very simple circuit composed of a delay line and a negative resistance element is presented and discussed.     

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 non-linear optimization procedure to estimate distances and instantaneous substitution rate matrices under the GTR model

MOTIVATION: The general-time-reversible (GTR) model is one of the most popular models of nucleotide substitution because it constitutes a good trade-off between mathematical tractability and biological reality. However, when it is applied for inferring evolutionary distances and/or instantaneous rate matrices, the GTR model seems more prone to inapplicability than more restrictive time-reversible models. Although it has been previously noted that the causes for intractability are caused by the impossibility of computing the logarithm of a matrix characterised by negative eigenvalues, the issue has not been investigated further. RESULTS: Here, we formally characterize the mathematical conditions, and discuss their biological interpretation, which lead to the inapplicability of the GTR model. We investigate the relations between, on one hand, the occurrence of negative eigenvalues and, on the other hand, both sequence length and sequence divergence. We then propose a possible re-formulation of previous procedures in terms of a non-linear optimization problem. We analytically investigate the effect of our approach on the estimated evolutionary distances and transition probability matrix. Finally, we provide an analysis on the goodness of the solution we propose. A numerical example is discussed.     

Dominant ionic mechanisms explored in spiking and bursting using local low-dimensional reductions of a biophysically realistic model neuron

The large number of variables involved in many biophysical models can conceal potentially simple dynamical mechanisms governing the properties of its solutions and the transitions between them as parameters are varied. To address this issue, we extend a novel model reduction method, based on “scales of dominance,” to multi-compartment models. We use this method to systematically reduce the dimension of a two-compartment conductance-based model of a crustacean pyloric dilator (PD) neuron that exhibits distinct modes of oscillation—tonic spiking, intermediate bursting and strong bursting. We divide trajectories into intervals dominated by a smaller number of variables, resulting in a locally reduced hybrid model whose dimension varies between two and six in different temporal regimes. The reduced model exhibits the same modes of oscillation as the 16 dimensional model over a comparable parameter range, and requires fewer ad hoc simplifications than a more traditional reduction to a single, globally valid model. The hybrid model highlights low-dimensional organizing structure in the dynamics of the PD neuron, and the dependence of its oscillations on parameters such as the maximal conductances of calcium currents. Our technique could be used to build hybrid low-dimensional models from any large multi-compartment conductance-based model in order to analyze the interactions between different modes of activity.     

神经元的确定性与随机性整数倍放电

在大鼠损伤背根节神经元的自发放电中发现了整数倍放电, 为了阐明这种放电所产生的原因, 首先研究神经元模型中确定性混沌所引起的整数倍放电与噪声所诱发的整数倍放电的峰峰间期（ＩＳＩ） 序列,通过分析得到前者的ＩＳＩ序列是非线性可预报的,具有确定的非线性特性,但由噪声所诱发的整数倍放电的ＩＳＩ序列是不可预报的, 这表明这两种机制所产生的整数倍放电具有不同的特点,存在着定性的差别,并且混沌运动所产生的整数倍放电是由混沌中各阶不稳定周期轨道决定的。从这种差别出发,分析了实验中整数倍放电的ＩＳＩ 序列,得到该ＩＳＩ 序列是可非线性预报的,这表明大鼠损伤背根节神经元自发放电中的整数倍放电更可能是由确定性机制所产生的     

Dependence of spontaneous neuronal firing and depolarisation block on astroglial membrane transport mechanisms

Exposed to a sufficiently high extracellular potassium concentration ([K^?+?]_o), the neuron can fire spontaneous discharges or even become inactivated due to membrane depolarisation (??depolarisation block??). Since these phenomena likely are related to the maintenance and propagation of seizure discharges, it is of considerable importance to understand the conditions under which excess [K^?+?]_o causes them. To address the putative effect of glial buffering on neuronal activity under elevated [K^?+?]_o conditions, we combined a recently developed dynamical model of glial membrane ion and water transport with a Hodgkin?CHuxley type neuron model. In this interconnected glia-neuron model we investigated the effects of natural heterogeneity or pathological changes in glial membrane transporter density by considering a large set of models with different, yet empirically plausible, sets of model parameters. We observed both the high [K^?+?]_o-induced duration of spontaneous neuronal firing and the prevalence of depolarisation block to increase when reducing the magnitudes of the glial transport mechanisms. Further, in some parameter regions an oscillatory bursting spiking pattern due to the dynamical coupling of neurons and glia was observed. Bifurcation analyses of the neuron model and of a simplified version of the neuron-glia model revealed further insights about the underlying mechanism behind these phenomena. The above insights emphasise the importance of combining neuron models with detailed astroglial models when addressing phenomena suspected to be influenced by the astroglia-neuron interaction. To facilitate the use of our neuron-glia model, a CellML version of it is made publicly available.     

Complex dynamics of a neuron model with discontinuous magnetic induction and exposed to external radiation

The last two decades have seen many literatures on the mathematical and computational analysis of neuronal activities resulting in many mathematical models to describe neuron. Many of those models have described the membrane potential of a neuron in terms of the leakage current and the synaptic inputs. Only recently researchers have proposed a new neuron model based on the electromagnetic induction theorem, which considers inner magnetic fluctuation and external electromagnetic radiation as a significant missing part that can participate in neural activity. While the flux coupling of the membrane is considered equivalent to a memductance function of a memristor, standard memductance model of \(\alpha + 3\beta \phi^{2}\) has been used in the literatures, but in this paper we propose a new memductance function based on discontinuous flux coupling. Various dynamical properties of the neuron model with discontinuous flux coupling are studied and interestingly the proposed model shows hyperchaotic behavior which was not identified in the literatures. Furthermore, we consider a ring network of the proposed model and investigate whether the chimera state can emerge. The chimera state relates to the state with simultaneously coherence and incoherence in oscillatory networks and has received much attention in recent years.     

Motoneuron membrane potentials follow a time inhomogeneous jump diffusion process

Stochastic leaky integrate-and-fire models are popular due to their simplicity and statistical tractability. They have been widely applied to gain understanding of the underlying mechanisms for spike timing in neurons, and have served as building blocks for more elaborate models. Especially the Ornstein–Uhlenbeck process is popular to describe the stochastic fluctuations in the membrane potential of a neuron, but also other models like the square-root model or models with a non-linear drift are sometimes applied. Data that can be described by such models have to be stationary and thus, the simple models can only be applied over short time windows. However, experimental data show varying time constants, state dependent noise, a graded firing threshold and time-inhomogeneous input. In the present study we build a jump diffusion model that incorporates these features, and introduce a firing mechanism with a state dependent intensity. In addition, we suggest statistical methods to estimate all unknown quantities and apply these to analyze turtle motoneuron membrane potentials. Finally, simulated and real data are compared and discussed. We find that a square-root diffusion describes the data much better than an Ornstein–Uhlenbeck process with constant diffusion coefficient. Further, the membrane time constant decreases with increasing depolarization, as expected from the increase in synaptic conductance. The network activity, which the neuron is exposed to, can be reasonably estimated to be a threshold version of the nerve output from the network. Moreover, the spiking characteristics are well described by a Poisson spike train with an intensity depending exponentially on the membrane potential.     

Two types of independent bursting mechanisms in inspiratory neurons: an integrative model

The network of coupled neurons in the pre-Bötzinger complex (pBC) of the medulla generates a bursting rhythm, which underlies the inspiratory phase of respiration. In some of these neurons, bursting persists even when synaptic coupling in the network is blocked and respiratory rhythmic discharge stops. Bursting in inspiratory neurons has been extensively studied, and two classes of bursting neurons have been identified, with bursting mechanism depends on either persistent sodium current or changes in intracellular Ca²⁺, respectively. Motivated by experimental evidence from these intrinsically bursting neurons, we present a two-compartment mathematical model of an isolated pBC neuron with two independent bursting mechanisms. Bursting in the somatic compartment is modeled via inactivation of a persistent sodium current, whereas bursting in the dendritic compartment relies on Ca²⁺ oscillations, which are determined by the neuromodulatory tone. The model explains a number of conflicting experimental results and is able to generate a robust bursting rhythm, over a large range of parameters, with a frequency adjusted by neuromodulators.     

Effects of polarization induced by non-weak electric fields on the excitability of elongated neurons with active dendrites

An externally-applied electric field can polarize a neuron, especially a neuron with elongated dendrites, and thus modify its excitability. Here we use a computational model to examine, predict, and explain these effects. We use a two-compartment Pinsky-Rinzel model neuron polarized by an electric potential difference imposed between its compartments, and we apply an injected ramp current. We vary three model parameters: the magnitude of the applied potential difference, the extracellular potassium concentration, and the rate of current injection. A study of the Time-To-First-Spike (TTFS) as a function of polarization leads to the identification of three regions of polarization strength that have different effects. In the weak region, the TTFS increases linearly with polarization. In the intermediate region, the TTFS increases either sub- or super-linearly, depending on the current injection rate and the extracellular potassium concentration. In the strong region, the TTFS decreases. Our results in the weak and strong region are consistent with experimental observations, and in the intermediate region, we predict novel effects that depend on experimentally-accessible parameters. We find that active channels in the dendrite play a key role in these effects. Our qualitative results were found to be robust over a wide range of inter-compartment conductances and the ratio of somatic to dendritic membrane areas. In addition, we discuss preliminary results where synaptic inputs replace the ramp injection protocol. The insights and conclusions were found to extend from our polarized PR model to a polarized PR model with I _h dendritic currents. Finally, we discuss the degree to which our results may be generalized.     

Stochastic modeling of the neuronal activity in the subthalamic nucleus and model parameter identification from Parkinson patient data

Several stochastic models, with various degrees of complexity, have been proposed to model the neuronal activity from different parts of the human brain. In this article, we use a simple Ornstein–Uhlenbeck process (OUP) to model the spike activity recorded from the subthalamic nucleus of patients suffering from Parkinson’s disease at the time of implantation of the electrodes for deep brain stimulation. From the recorded data, which contains information about the spike times of a single neuron, we identify and extract the model parameters of the OUP. We then use these parameters to numerically simulate the inter-spike intervals and the voltage across the neuron membrane. We finally assess how well the proposed mathematical model fits to the measured data and compare it with other commonly adopted stochastic models. We show an excellent agreement between the computer-generated data according to the OUP model and the measured one, as well as the superiority of the OUP model when compared to the Poisson process model and the random walk model; thus, establishing the validity of the OUP as a simple yet biologically plausible model of the neuronal activity recorded from the subthalamic nucleus of Parkinson’s disease patients.