Reliability and stochastic synchronization in type I vs. type II neural oscillators

Biophysical models of intrinsic homeostasis: Firing rates and beyond

In view of ever-changing conditions both in the external world and in intrinsic brain states, maintaining the robustness of computations poses a challenge, adequate solutions to which we are only beginning to understand. At the level of cell-intrinsic properties, biophysical models of neurons permit one to identify relevant physiological substrates that can serve as regulators of neuronal excitability and to test how feedback loops can stabilize crucial variables such as long-term calcium levels and firing rates. Mathematical theory has also revealed a rich set of complementary computational properties arising from distinct cellular dynamics and even shaping processing at the network level. Here, we provide an overview over recently explored homeostatic mechanisms derived from biophysical models and hypothesize how multiple dynamical characteristics of cells, including their intrinsic neuronal excitability classes, can be stably controlled.

Dopamine Neurons Change the Type of Excitability in Response to Stimuli

The dynamics of neuronal excitability determine the neuron's response to stimuli, its synchronization and resonance properties and, ultimately, the computations it performs in the brain. We investigated the dynamical mechanisms underlying the excitability type of dopamine (DA) neurons, using a conductance-based biophysical model, and its regulation by intrinsic and synaptic currents. Calibrating the model to reproduce low frequency tonic firing results in N-methyl-D-aspartate (NMDA) excitation balanced by γ-Aminobutyric acid (GABA)-mediated inhibition and leads to type I excitable behavior characterized by a continuous decrease in firing frequency in response to hyperpolarizing currents. Furthermore, we analyzed how excitability type of the DA neuron model is influenced by changes in the intrinsic current composition. A subthreshold sodium current is necessary for a continuous frequency decrease during application of a negative current, and the low-frequency "balanced" state during simultaneous activation of NMDA and GABA receptors. Blocking this current switches the neuron to type II characterized by the abrupt onset of repetitive firing. Enhancing the anomalous rectifier Ih current also switches the excitability to type II. Key characteristics of synaptic conductances that may be observed in vivo also change the type of excitability: a depolarized γ-Aminobutyric acid receptor (GABAR) reversal potential or co-activation of α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) leads to an abrupt frequency drop to zero, which is typical for type II excitability. Coactivation of N-methyl-D-aspartate receptors (NMDARs) together with AMPARs and GABARs shifts the type I/II boundary toward more hyperpolarized GABAR reversal potentials. To better understand how altering each of the aforementioned currents leads to changes in excitability profile of DA neuron, we provide a thorough dynamical analysis. Collectively, these results imply that type I excitability in dopamine neurons might be important for low firing rates and fine-tuning basal dopamine levels, while switching excitability to type II during NMDAR and AMPAR activation may facilitate a transient increase in dopamine concentration, as type II neurons are more amenable to synchronization by mutual excitation.

Stimulus features, resetting curves, and the dependence on adaptation

We derive a formula that relates the spike-triggered covariance (STC) to the phase resetting curve (PRC) of a neural oscillator. We use this to show how changes in the shape of the PRC alter the sensitivity of the neuron to different stimulus features, which are the eigenvectors of the STC. We compute the PRC and STC for some biophysical models. We compare the STCs and their spectral properties for a two-parameter family of PRCs. Surprisingly, the skew of the PRC has a larger effect on the spectrum and shape of the STC than does the bimodality of the PRC (which plays a large role in synchronization properties). Finally, we relate the STC directly to the spike-triggered average and apply this theory to an olfactory bulb mitral cell recording.

Identifying type I excitability using dynamics of stochastic neural firing patterns

The stochastic firing patterns are simulated near saddle-node bifurcation on an invariant cycle corresponding to type I excitability in stochastic Morris-Lecar model. In absence of external periodic signal, the stochastic firing manifests continuous distribution in ISI histogram (ISIH), whose amplitude at first increases sharply and then decreases exponentially. In presence of the external periodic signal, stochastic firing patterns appear as two cases of integer multiple firing with multiple discrete peaks in ISIH. One manifests perfect exponential decay in all peaks and the other imperfect exponential decay except a lower first peak. These stochastic firing patterns simulated with or without external periodic signal can be demonstrated in the experiments on rat hippocampal CA1 pyramidal neurons. The exponential decay laws in the multiple peaks are also acquired using probability analysis method. The perfect decay law is determined by the independent characteristic within the firing while the imperfect decay law is from the inhibitory effect. In addition, the stochastic firing patterns corresponding to type I excitability are compared to those of type II excitability. The results not only reveal the dynamics of stochastic firing patterns with or without external signal corresponding to type I excitability, but also provide practical indicators to availably identify type I excitability.

Robustness, variability, phase dependence, and longevity of individual synaptic input effects on spike timing during fluctuating synaptic backgrounds: a modeling study of globus pallidus neuron phase response properties

A neuron's phase response curve (PRC) shows how inputs arriving at different times during the spike cycle differentially affect the timing of subsequent spikes. Using a full morphological model of a globus pallidus (GP) neuron, we previously demonstrated that dendritic conductances shape the PRC in a spike frequency-dependent manner, suggesting different functional roles of perisomatic and distal dendritic synapses in the control of patterned network activity. In the present study we extend this analysis to examine the impact of physiologically realistic high conductance states on somatic and dendritic PRCs and the time course of spike train perturbations. First, we found that average somatic and dendritic PRCs preserved their shapes and spike frequency dependence when the model was driven by spatially-distributed, stochastic conductance inputs rather than tonic somatic current. However, responses to inputs during specific synaptic backgrounds often deviated substantially from the average PRC. Therefore, we analyzed the interactions of PRC stimuli with transient fluctuations in the synaptic background on a trial-by-trial basis. We found that the variability in responses to PRC stimuli and the incidence of stimulus-evoked added or skipped spikes were stimulus-phase-dependent and reflected the profile of the average PRC, suggesting commonality in the underlying mechanisms. Clear differences in the relation between the phase of input and variability of spike response between dendritic and somatic inputs indicate that these regions generally represent distinct dynamical subsystems of synaptic integration with respect to influencing the stability of spike time attractors generated by the overall synaptic conductance.

Spike-timing-dependent plasticity and reliability optimization: the role of neuron dynamics

Plastic changes in synaptic efficacy can depend on the time ordering of presynaptic and postsynaptic spikes. This phenomenon is called spike-timing-dependent plasticity (STDP). One of the most striking aspects of this plasticity mechanism is that the STDP windows display a great variety of forms in different parts of the nervous system. We explore this issue from a theoretical point of view. We choose as the optimization principle the minimization of conditional entropy or maximization of reliability in the transmission of information. We apply this principle to two types of postsynaptic dynamics, designated type I and type II. The first is characterized as being an integrator, while the second is a resonator. We find that, depending on the parameters of the models, the optimization principle can give rise to a wide variety of STDP windows, such as antisymmetric Hebbian, predominantly depressing or symmetric with one positive region and two lateral negative regions. We can relate each of these forms to the dynamical behavior of the different models. We also propose experimental tests to assess the validity of the optimization principle.

Phase-response curves and synchronized neural networks

We review the principal assumptions underlying the application of phase-response curves (PRCs) to synchronization in neuronal networks. The PRC measures how much a given synaptic input perturbs spike timing in a neural oscillator. Among other applications, PRCs make explicit predictions about whether a given network of interconnected neurons will synchronize, as is often observed in cortical structures. Regarding the assumptions of the PRC theory, we conclude: (i) The assumption of noise-tolerant cellular oscillations at or near the network frequency holds in some but not all cases. (ii) Reduced models for PRC-based analysis can be formally related to more realistic models. (iii) Spike-rate adaptation limits PRC-based analysis but does not invalidate it. (iv) The dependence of PRCs on synaptic location emphasizes the importance of improving methods of synaptic stimulation. (v) New methods can distinguish between oscillations that derive from mutual connections and those arising from common drive. (vi) It is helpful to assume linear summation of effects of synaptic inputs; experiments with trains of inputs call this assumption into question. (vii) Relatively subtle changes in network structure can invalidate PRC-based predictions. (viii) Heterogeneity in the preferred frequencies of component neurons does not invalidate PRC analysis, but can annihilate synchronous activity.

A new approach for determining phase response curves reveals that Purkinje cells can act as perfect integrators

Cerebellar Purkinje cells display complex intrinsic dynamics. They fire spontaneously, exhibit bistability, and via mutual network interactions are involved in the generation of high frequency oscillations and travelling waves of activity. To probe the dynamical properties of Purkinje cells we measured their phase response curves (PRCs). PRCs quantify the change in spike phase caused by a stimulus as a function of its temporal position within the interspike interval, and are widely used to predict neuronal responses to more complex stimulus patterns. Significant variability in the interspike interval during spontaneous firing can lead to PRCs with a low signal-to-noise ratio, requiring averaging over thousands of trials. We show using electrophysiological experiments and simulations that the PRC calculated in the traditional way by sampling the interspike interval with brief current pulses is biased. We introduce a corrected approach for calculating PRCs which eliminates this bias. Using our new approach, we show that Purkinje cell PRCs change qualitatively depending on the firing frequency of the cell. At high firing rates, Purkinje cells exhibit single-peaked, or monophasic PRCs. Surprisingly, at low firing rates, Purkinje cell PRCs are largely independent of phase, resembling PRCs of ideal non-leaky integrate-and-fire neurons. These results indicate that Purkinje cells can act as perfect integrators at low firing rates, and that the integration mode of Purkinje cells depends on their firing rate.

By observing how brief current pulses injected at different times between spikes change the phase of spiking of a neuron (and thus obtaining the so-called phase response curve), it should be possible to predict a full spike train in response to more complex stimulation patterns. When we applied this traditional protocol to obtain phase response curves in cerebellar Purkinje cells in the presence of noise, we observed a triangular region devoid of data points near the end of the spiking cycle. This “Bermuda Triangle” revealed a flaw in the classical method for constructing phase response curves. We developed a new approach to eliminate this flaw and used it to construct phase response curves of Purkinje cells over a range of spiking rates. Surprisingly, at low firing rates, phase changes were independent of the phase of the injected current pulses, implying that the Purkinje cell is a perfect integrator under these conditions. This mechanism has not yet been described in other cell types and may be crucial for the information processing capabilities of these neurons.

Reliability, synchrony and noise

The brain is noisy. Neurons receive tens of thousands of highly fluctuating inputs and generate spike trains that appear highly irregular. Much of this activity is spontaneous - uncoupled to overt stimuli or motor outputs - leading to questions about the functional impact of this noise. Although noise is most often thought of as disrupting patterned activity and interfering with the encoding of stimuli, recent theoretical and experimental work has shown that noise can play a constructive role - leading to increased reliability or regularity of neuronal firing in single neurons and across populations. These results raise fundamental questions about how noise can influence neural function and computation.

Class-II neurons display a higher degree of stochastic synchronization than class-I neurons

We describe the relationship between the shape of the phase-resetting curve (PRC) and the degree of stochastic synchronization observed between a pair of uncoupled general oscillators receiving partially correlated Poisson inputs in addition to inputs from independent sources. We use perturbation methods to derive an expression relating the shape of the PRC to the probability density function (PDF) of the phase difference between the oscillators. We compute various measures of the degree of synchrony and cross correlation from the PDF's and use the same to compare and contrast differently shaped PRCs, with respect to their ability to undergo stochastic synchronization. Since the shape of the PRC depends on underlying dynamical details of the oscillator system, we utilize the results obtained from the analysis of general oscillator systems to study specific models of neuronal oscillators. It is shown that the degree of stochastic synchronization is controlled both by the firing rate of the neuron and the membership of the PRC (type I or type II). It is also shown that the circular variance for the integrate and fire neuron and the generalized order parameter for a hippocampal interneuron model have a nonlinear relationship to the input correlation.

Stochastic dynamics of uncoupled neural oscillators: Fokker-Planck studies with the finite element method

We have investigated the effect of the phase response curve on the dynamics of oscillators driven by noise in two limit cases that are especially relevant for neuroscience. Using the finite element method to solve the Fokker-Planck equation we have studied (i) the impact of noise on the regularity of the oscillations quantified as the coefficient of variation, (ii) stochastic synchronization of two uncoupled phase oscillators driven by correlated noise, and (iii) their cross-correlation function. We show that, in general, the limit of type II oscillators is more robust to noise and more efficient at synchronizing by correlated noise than type I.