Phase-response curves give the responses of neurons to transient inputs

Phase response curve analysis of a full morphological globus pallidus neuron model reveals distinct perisomatic and dendritic modes of synaptic integration

Synchronization of globus pallidus (GP) neurons and cortically entrained oscillations between GP and other basal ganglia nuclei are key features of the pathophysiology of Parkinson's disease. Phase response curves (PRCs), which tabulate the effects of phasic inputs within a neuron's spike cycle on output spike timing, are efficient tools for predicting the emergence of synchronization in neuronal networks and entrainment to periodic input. In this study we apply physiologically realistic synaptic conductance inputs to a full morphological GP neuron model to determine the phase response properties of the soma and different regions of the dendritic tree. We find that perisomatic excitatory inputs delivered throughout the interspike interval advance the phase of the spontaneous spike cycle yielding a type I PRC. In contrast, we demonstrate that distal dendritic excitatory inputs can either delay or advance the next spike depending on whether they occur early or late in the spike cycle. We find this latter pattern of responses, summarized by a biphasic (type II) PRC, was a consequence of dendritic activation of the small conductance calcium-activated potassium current, SK. We also evaluate the spike-frequency dependence of somatic and dendritic PRC shapes, and we demonstrate the robustness of our results to variations of conductance densities, distributions, and kinetic parameters. We conclude that the distal dendrite of GP neurons embodies a distinct dynamical subsystem that could promote synchronization of pallidal networks to excitatory inputs. These results highlight the need to consider different effects of perisomatic and dendritic inputs in the control of network behavior.

Control of the firing patterns of vibrissa motoneurons by modulatory and phasic synaptic inputs: a modeling study

Vibrissa motoneurons (vMNs) generate rhythmic firing that controls whisker movements, even without cortical, cerebellar, or sensory inputs. vMNs receive serotonergic modulation from brain stem areas, which mainly increases their persistent sodium conductance (g(NaP)) and, possibly, phasic input from a putative central pattern generator (CPG). In response to serotonergic modulation or just-suprathreshold current steps, vMNs fire at low rates, below the firing frequency of exploratory whisking. In response to periodic inputs, vMNs exhibit nonlinear suprathreshold resonance in frequency ranges of exploratory whisking. To determine how firing patterns of vMNs are determined by their 1) intrinsic ionic conductances and 2) responses to periodic input from a putative CPG and to serotonergic modulation, we construct and analyze a single-compartment, conductance-based model of vMNs. Low firing rates are supported in extended regimes by adaptation currents and the minimal firing rate decreases with g(NaP) and increases with M-potassium and h-cation conductances. Suprathreshold resonance results from the locking properties of vMN firing to stimuli and from reduction of firing rates at low frequencies by slow M and afterhyperpolarization potassium conductances. h conductance only slightly affects the suprathreshold resonance. When a vMN is subjected to a small periodic CPG input, serotonergically induced g(NaP) elevation may transfer the system from quiescence to a firing state that is highly locked to the CPG input. Thus we conclude that for vMNs, the CPG controls firing frequency and phase and enables bursting, whereas serotonergic modulation controls transitions from quiescence to firing unless the CPG input is sufficiently strong.

Synchrony with shunting inhibition in a feedforward inhibitory network

Recent experiments have shown that GABA(A) receptor mediated inhibition in adult hippocampus is shunting rather than hyperpolarizing. Simulation studies of realistic interneuron networks with strong shunting inhibition have been demonstrated to exhibit robust gamma band (20-80 Hz) synchrony in the presence of heterogeneity in the intrinsic firing rates of individual neurons in the network. In order to begin to understand how shunting can contribute to network synchrony in the presence of heterogeneity, we develop a general theoretical framework using spike time response curves (STRC's) to study patterns of synchrony in a simple network of two unidirectionally coupled interneurons (UCI network) interacting through a shunting synapse in the presence of heterogeneity. We derive an approximate discrete map to analyze the dynamics of synchronous states in the UCI network by taking into account the nonlinear contributions of the higher order STRC terms. We show how the approximate discrete map can be used to successfully predict the domain of synchronous 1:1 phase locked state in the UCI network. The discrete map also allows us to determine the conditions under which the two interneurons can exhibit in-phase synchrony. We conclude by demonstrating how the information from the study of the discrete map for the dynamics of the UCI network can give us valuable insight into the degree of synchrony in a larger feed-forward network of heterogeneous interneurons.

Fluctuation of gamma-band phase synchronization within the auditory cortex in schizophrenia

OBJECTIVE

To study the phase stability of the 40Hz auditory steady-state response (ASSR) in Sz, and in addition, to investigate inter-hemispheric phase synchronization using ipsilateral and contralateral hemisphere gamma band ASSRs.

METHODS

Whole head magnetoencephalography (MEG) was used to detect ASSR from both hemispheres in Sz patients and their control counterparts. Source localization, spatial and temporal filtering were performed to infer gamma band activity from the neural generators of the ASSR. The response gamma band phase stability relative to a reference signal was quantified using the phase synchronization index (PSI).

RESULTS

Results indicated reduced phase synchronization of the ASSR and the stimulus reference signal in Sz patients compared to control subjects, in addition to reduced inter-hemispheric phase synchronization between contralateral and ipsilateral hemispheric responses in Sz patients.

CONCLUSIONS

Greater intra and inter hemispheric fluctuations of ASSR gamma band phase synchronization in Sz add to previous studies suggesting timing deficiencies within neural populations, possibly caused by impairments of neural network parameters.

SIGNIFICANCE

This study provides experimental support that may aid in understanding the dynamics of neural phase synchrony caused by modifications of underlying neurotransmitter systems, as reflected in disease states such as schizophrenia.

Spike timing-dependent plasticity is affected by the interplay of intrinsic and network oscillations

Spike timing-dependent plasticity (STDP) is a form of Hebbian learning which is thought to underlie structure formation during development, and learning and memory in later life. In this paper we show that the intrinsic properties of the postsynaptic neuron might have a deep influence on STDP dynamics by shaping the causal correlation between the pre- and the postsynaptic spike trains. The cell-specific effect of STDP is particularly evident in the presence of an oscillatory component in a cell input. In this case, the cell-specific phase response to an oscillatory modulation biases the oscillating afferents towards potentiation or depression, depending upon the intrinsic dynamics of the postsynaptic neuron and the period of the modulation.

Specific entrainment of mitral cells during gamma oscillation in the rat olfactory bulb

Local field potential (LFP) oscillations are often accompanied by synchronization of activity within a widespread cerebral area. Thus, the LFP and neuronal coherence appear to be the result of a common mechanism that underlies neuronal assembly formation. We used the olfactory bulb as a model to investigate: (1) the extent to which unitary dynamics and LFP oscillations can be correlated and (2) the precision with which a model of the hypothesized underlying mechanisms can accurately explain the experimental data. For this purpose, we analyzed simultaneous recordings of mitral cell (MC) activity and LFPs in anesthetized and freely breathing rats in response to odorant stimulation. Spike trains were found to be phase-locked to the gamma oscillation at specific firing rates and to form odor-specific temporal patterns. The use of a conductance-based MC model driven by an approximately balanced excitatory-inhibitory input conductance and a relatively small inhibitory conductance that oscillated at the gamma frequency allowed us to provide one explanation of the experimental data via a mode-locking mechanism. This work sheds light on the way network and intrinsic MC properties participate in the locking of MCs to the gamma oscillation in a realistic physiological context and may result in a particular time-locked assembly. Finally, we discuss how a self-synchronization process with such entrainment properties can explain, under experimental conditions: (1) why the gamma bursts emerge transiently with a maximal amplitude position relative to the stimulus time course; (2) why the oscillations are prominent at a specific gamma frequency; and (3) why the oscillation amplitude depends on specific stimulus properties. We also discuss information processing and functional consequences derived from this mechanism.

Olfactory function relies on a chain of neural relays that extends from the periphery to the central nervous system and implies neural activity with various timescales. A central question in neuroscience is how information is encoded by the neural activity. In the mammalian olfactory bulb, local neural activity oscillations in the 40–80 Hz range (gamma) may influence the timing of individual neuron activities such that olfactory information may be encoded in this way. In this study, we first characterize in vivo the detailed activity of individual neurons relative to the oscillation and find that, depending on their state, neurons can exhibit periodic activity patterns. We also find, at least qualitatively, a relation between this activity and a particular odor. This is reminiscent of general physical phenomena—the entrainment by an oscillation—and to verify this hypothesis, in a second phase, we build a biologically realistic model mimicking these in vivo conditions. Our model confirms quantitatively this hypothesis and reveals that entrainment is maximal in the gamma range. Taken together, our results suggest that the neuronal activity may be specifically formatted in time during the gamma oscillation in such a way that it could, at this stage, encode the odor.

Adaptive synchronization of activities in a recurrent network

Predictive learning rules, where synaptic changes are driven by the difference between a random input and its reconstruction derived from internal variables, have proven to be very stable and efficient. However, it is not clear how such learning rules could take place in biological synapses. Here we propose an implementation that exploits the synchronization of neural activities within a recurrent network. In this framework, the asymmetric shape of spike-timing-dependent plasticity (STDP) can be interpreted as a self-stabilizing mechanism. Our results suggest a novel hypothesis concerning the computational role of neural synchrony and oscillations.

Type-II phase resetting curve is optimal for stochastic synchrony

The phase-resetting curve (PRC) describes the response of a neural oscillator to small perturbations in membrane potential. Its usefulness for predicting the dynamics of weakly coupled deterministic networks has been well characterized. However, the inputs to real neurons may often be more accurately described as barrages of synaptic noise. Effective connectivity between cells may thus arise in the form of correlations between the noisy input streams. We use constrained optimization and perturbation methods to prove that the PRC shape determines susceptibility to synchrony among otherwise uncoupled noise-driven neural oscillators. PRCs can be placed into two general categories: type-I PRCs are non-negative, while type-II PRCs have a large negative region. Here we show that oscillators with type-II PRCs receiving common noisy input synchronize more readily than those with type-I PRCs.

Predicting synchronous and asynchronous network groupings of hippocampal interneurons coupled with dendritic gap junctions

Direct electrical communication between central nervous system (CNS) neurons including those in the hippocampus is well-established. This form of communication is mediated by gap junctions and it is known that this coupling is important for brain rhythms such as gamma (20-80 Hz) which occur during active behavioural states. It is also known that gap junctions are present at several locations along the dendrites of hippocampal interneurons including parvalbumin-positive basket cell types. Weakly coupled oscillator theory, which uses phase response curves (PRCs), has been used to understand and predict the dynamics of electrically coupled networks. Here we use compartmental models of hippocampal basket cells with different levels of basal and apical spike attenuation together with the theory to show that network output can be broken down into three groupings: synchronous, asynchronous and antiphase-like patterns. Moreover, quantified PRCs can be used as a rule of thumb to determine the occurrence of a particular grouping under weak coupling conditions, which in turn implies that spike delays are critical factors in determining network output. In moving beyond weak coupling to encompass the full physiological regime of coupling strengths with network simulations, we note that it is important to be able to differentiate between these different groupings as it affects how the network responds with modulation. Specifically, an asynchronous grouping provides more dynamic richness as a larger range of phase-locked states can be expressed with strength changes. From a functional viewpoint it may be that modulation of electrically coupled networks are key to controlling cell assemblies that contribute to information coding brain substrates.

Linear control of neuronal spike timing using phase response curves

We propose a simple, robust, linear method to control the spike timing of a periodically firing neuron. The control scheme uses the neuron's phase response curve to identify an area of optimal sensitivity for the chosen stimulation parameters. The spike advance as a function of current pulse amplitude is characterized at the optimal phase and a linear least-squares regression is fit to the data. The inverted regression is used as the control function for this method. The efficacy of this method is demonstrated through numerical simulations of a Hodgkin-Huxley style neuron model as well as in real neurons from rat hippocampal slice preparations. The study shows a proof of concept for the application of a linear control scheme to control neuron spike timing in-vitro. This study was done on an individual cell level, but translation to a tissue or network level is possible. Control schemes of this type could be implemented in a closed loop implantable device to treat neuromotor disorders involving pathologically neuronal activity such as epilepsy or Parkinson's disease.

Integration of broadband conductance input in rat somatosensory cortical inhibitory interneurons: an inhibition-controlled switch between intrinsic and input-driven spiking in fast-spiking cells

Quantitative understanding of the dynamics of particular cell types when responding to complex, natural inputs is an important prerequisite for understanding the operation of the cortical network. Different types of inhibitory neurons are connected by electrical synapses to nearby neurons of the same type, enabling the formation of synchronized assemblies of neurons with distinct dynamical behaviors. Under what conditions is spike timing in such cells determined by their intrinsic dynamics and when is it driven by the timing of external input? In this study, we have addressed this question using a systematic approach to characterizing the input-output relationships of three types of cortical interneurons (fast spiking [FS], low-threshold spiking [LTS], and nonpyramidal regular-spiking [NPRS] cells) in the rat somatosensory cortex, during fluctuating conductance input designed to mimic natural complex activity. We measured the shape of average conductance input trajectories preceding spikes and fitted a two-component linear model of neuronal responses, which included an autoregressive term from its own output, to gain insight into the input-output relationships of neurons. This clearly separated the contributions of stimulus and discharge history, in a cell-type dependent manner. Unlike LTS and NPRS cells, FS cells showed a remarkable switch in dynamics, from intrinsically driven spike timing to input-fluctuation-controlled spike timing, with the addition of even a small amount of inhibitory conductance. Such a switch could play a pivotal role in the function of FS cells in organizing coherent gamma oscillations in the local cortical network. Using both pharmacological perturbations and modeling, we show how this property is a consequence of the particular complement of voltage-dependent conductances in these cells.

Cholinergic neuromodulation changes phase response curve shape and type in cortical pyramidal neurons

Spike generation in cortical neurons depends on the interplay between diverse intrinsic conductances. The phase response curve (PRC) is a measure of the spike time shift caused by perturbations of the membrane potential as a function of the phase of the spike cycle of a neuron. Near the rheobase, purely positive (type I) phase-response curves are associated with an onset of repetitive firing through a saddle-node bifurcation, whereas biphasic (type II) phase-response curves point towards a transition based on a Hopf-Andronov bifurcation. In recordings from layer 2/3 pyramidal neurons in cortical slices, cholinergic action, consistent with down-regulation of slow voltage-dependent potassium currents such as the M-current, switched the PRC from type II to type I. This is the first report showing that cholinergic neuromodulation may cause a qualitative switch in the PRCs type implying a change in the fundamental dynamical mechanism of spike generation.

Detection of subthreshold pulses in neurons with channel noise

Neurons are subject to various kinds of noise. In addition to synaptic noise, the stochastic opening and closing of ion channels represents an intrinsic source of noise that affects the signal-processing properties of the neuron. We study the response of a stochastic Hodgkin-Huxley neuron to transient input subthreshold pulses. It is found that the average response time decreases but variance increases as the amplitude of channel noise increases. In the case of single-pulse detection, we show that channel noise enables one neuron to detect the subthreshold signals and an optimal membrane area (or channel noise intensity) exists for a single neuron to achieve optimal performance. However, the detection ability of a single neuron is limited by large errors. Here, we test a simple neuronal network that can enhance the pulse-detecting abilities of neurons and find that dozens of neurons can perfectly detect subthreshold pulses. The phenomenon of intrinsic stochastic resonance is also found at both the level of single neurons and the level of networks. At the network level, the detection ability of networks can be optimized for the number of neurons comprising the network.

The effects of cholinergic neuromodulation on neuronal phase-response curves of modeled cortical neurons

The response of an oscillator to perturbations is described by its phase-response curve (PRC), which is related to the type of bifurcation leading from rest to tonic spiking. In a recent experimental study, we have shown that the type of PRC in cortical pyramidal neurons can be switched by cholinergic neuromodulation from type II (biphasic) to type I (monophasic). We explored how intrinsic mechanisms affected by acetylcholine influence the PRC using three different types of neuronal models: a theta neuron, single-compartment neurons and a multi-compartment neuron. In all of these models a decrease in the amount of a spike-frequency adaptation current was a necessary and sufficient condition for the shape of the PRC to change from biphasic (type II) to purely positive (type I).

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.

Dendritic excitability and synaptic plasticity

Most synaptic inputs are made onto the dendritic tree. Recent work has shown that dendrites play an active role in transforming synaptic input into neuronal output and in defining the relationships between active synapses. In this review, we discuss how these dendritic properties influence the rules governing the induction of synaptic plasticity. We argue that the location of synapses in the dendritic tree, and the type of dendritic excitability associated with each synapse, play decisive roles in determining the plastic properties of that synapse. Furthermore, since the electrical properties of the dendritic tree are not static, but can be altered by neuromodulators and by synaptic activity itself, we discuss how learning rules may be dynamically shaped by tuning dendritic function. We conclude by describing how this reciprocal relationship between plasticity of dendritic excitability and synaptic plasticity has changed our view of information processing and memory storage in neuronal networks.

Bursting regimes in map-based neuron models coupled through fast threshold modulation

A system consisting of two map-based neurons coupled through reciprocal excitatory or inhibitory chemical synapses is discussed. After a brief explanation of the basic mechanism behind generation and synchronization of bursts, parameter space is explored to determine less obvious but biologically meaningful regimes and effects. Among them, we show how excitatory synapses without any delays may induce antiphase synchronization; that a synapse may change its character from excitatory to inhibitory and vice versa by changing its conductance, without any change in reversal potential; and that small variations in the synaptic threshold may result in drastic changes in the synchronization of spikes within bursts. Finally we show how the synchronization effects found in the two-neuron system carry over to larger networks.

Gamma oscillations and stimulus selection

More coherent excitatory stimuli are known to have a competitive advantage over less coherent ones. We show here that this advantage is amplified greatly when the target includes inhibitory interneurons acting via GABA(A)-receptor-mediated synapses and the coherent input oscillates at gamma frequency. We hypothesize that therein lies, at least in part, the functional significance of the experimentally observed link between attentional biasing of stimulus competition and gamma frequency rhythmicity.

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.

The dynamical response properties of neocortical neurons to temporally modulated noisy inputs in vitro

Cortical neurons are often classified by current-frequency relationship. Such a static description is inadequate to interpret neuronal responses to time-varying stimuli. Theoretical studies suggested that single-cell dynamical response properties are necessary to interpret ensemble responses to fast input transients. Further, it was shown that input-noise linearizes and boosts the response bandwidth, and that the interplay between the barrage of noisy synaptic currents and the spike-initiation mechanisms determine the dynamical properties of the firing rate. To test these model predictions, we estimated the linear response properties of layer 5 pyramidal cells by injecting a superposition of a small-amplitude sinusoidal wave and a background noise. We characterized the evoked firing probability across many stimulation trials and a range of oscillation frequencies (1-1000 Hz), quantifying response amplitude and phase-shift while changing noise statistics. We found that neurons track unexpectedly fast transients, as their response amplitude has no attenuation up to 200 Hz. This cut-off frequency is higher than the limits set by passive membrane properties (approximately 50 Hz) and average firing rate (approximately 20 Hz) and is not affected by the rate of change of the input. Finally, above 200 Hz, the response amplitude decays as a power-law with an exponent that is independent of voltage fluctuations induced by the background noise.

Partial phase synchronization of neural populations due to random Poisson inputs

We show that populations of identical uncoupled neurons exhibit partial phase synchronization when stimulated with independent, random unidirectional current spikes with interspike time intervals drawn from a Poisson distribution. We characterize this partial synchronization using the phase distribution of the population, and consider analytical approximations and numerical simulations of phase-reduced models and the corresponding conductance-based models of typical Type I (Hindmarsh-Rose) and Type II (Hodgkin-Huxley) neurons, showing quantitatively how the extent of the partial phase synchronization depends on the magnitude and mean interspike frequency of the stimulus. Furthermore, we present several simple examples that disprove the notion that phase synchrony must be strongly related to spike synchrony. Instead, the importance of partial phase synchrony is shown to lie in its influence on the response of the population to stimulation, which we illustrate using first spike time histograms.

Relating neural dynamics to neural coding

We demonstrate that two key theoretical objects used widely in computational neuroscience, the phase-resetting curve (PRC) from dynamics and the spike triggered average (STA) from statistical analysis, are closely related when neurons fire in a nearly regular manner and the stimulus is sufficiently small. We prove that the STA due to injected noisy current is proportional to the derivative of the PRC. We compare these analytic results with numerical calculations for the Hodgkin-Huxley neuron and we apply the method to neurons in the olfactory bulb of mice. This observation allows us to relate the stimulus-response properties of a neuron to its dynamics, bridging the gap between dynamical and information theoretic approaches to understanding brain computations and facilitating the interpretation of changes in channels and other cellular properties as influencing the representation of stimuli.

Generation of theta oscillations by weakly coupled neural oscillators in the presence of noise

Neuronal oscillations are a robust phenomenon occurring in a variety of brain regions despite considerable amounts of noise. In this article classical phase-response theory is generalized to the case of noisy weak-coupling regimes by deriving an iterated map for the asynchrony of spikes in an oscillation cycle. Two criteria are introduced to check the validity of our approximations: One criterion tests the assumption that all neurons fire exactly once per cycle, the other criterion tests for linearity. The framework is applied to stellate cells of the medial entorhinal cortex layer II. We find that rhythmogenesis is more robust in the case of excitatory noise as compared to inhibitory noise. It is shown that a network of stellate cells can also act as a generator of theta if the neurons are connected via a fast-oscillating network of inhibitory interneurons.

Layer and frequency dependencies of phase response properties of pyramidal neurons in rat motor cortex

It is postulated that synchronous firing of cortical neurons plays an active role in cognitive functions of the brain. An important issue is whether pyramidal neurons in different cortical layers exhibit similar tendencies to synchronise. To address this issue, we performed intracellular and whole-cell recordings of regular-spiking pyramidal neurons in slice preparations of the rat motor cortex (18-45 days old) and analysed the phase response curves of these pyramidal neurons in layers 2/3 and 5. The phase response curve represents how an external stimulus affects the timing of spikes immediately after the stimulus in repetitively firing neurons. The phase response curve can be classified into two categories, type 1 (the spike is always advanced) and type 2 (the spike is advanced or delayed depending on the stimulus phase), and are important determinants of whether or not rhythmic synchronization of neuron pairs occurs. We found that pyramidal neurons in layer 2/3 tend to display type-2 phase response curves whereas those in layer 5 tend to exhibit type-1 phase response curves. The differences were prominent particularly in the gamma-frequency range (20-45 Hz). Our results imply that the layer-2/3 pyramidal neurons, when coupled mutually through fast excitatory synapses, may exhibit a much stronger tendency for rhythmic synchronization than layer-5 neurons in the gamma-frequency range.

Transition between two excitabilities in mesencephalic V neurons

Neurons can make different responses to identical inputs. According to the emerging frequency of repetitive firing, neurons are classified into two types: type 1 and type 2 excitability. Though in mathematical simulations, minor modifications of parameters describing ionic currents can result in transitions between these two excitabilities, empirical evidence to support these theoretical possibilities is scarce. Here we report a joint theoretical and experimental study to test the hypothesis that changes in parameters describing ionic currents cause predictable transitions between the two excitabilities in mesencephalic V (Mes V) neurons. We developed a simple mathematical model of Mes V neurons. Using bifurcation analysis and model simulation, we then predicted that changes in conductance of two low-threshold currents would result in transitions between type 1 and type 2. Finally, by applying specific channel blockers, we observed the transition between two excitabilities forecast by the mathematical model.

Cortical pyramidal cells as non-linear oscillators: experiment and spike-generation theory

Cortical neurons are capable of generating trains of action potentials in response to current injections. These discharges can take different forms, e.g., repetitive firing that adapts during the period of current injection or bursting behaviors. We have used a combined experimental and computational approach to characterize the dynamics leading to action potential responses in single neurons. Specifically we investigated the origin of complex firing patterns in response to sinusoidal current injections. Using a reduced model, the theta-neuron, alongside recordings from cortical pyramidal cells we show that both real and simulated neurons show phase-locking to sine wave stimuli up to a critical frequency, above which period skipping and 1-to-x phase-locking occurs. The locking behavior follows a complex "devil's staircase" phenomena, where locked modes are interleaved with irregular firing. We further show that the critical frequency depends on the time scale of spike generation and on the level of spike frequency adaptation. These results suggest that phase-locking of neuronal responses to complex input patterns can be explained by basic properties of the spike-generating machinery.

Mechanisms of firing patterns in fast-spiking cortical interneurons

Cortical fast-spiking (FS) interneurons display highly variable electrophysiological properties. Their spike responses to step currents occur almost immediately following the step onset or after a substantial delay, during which subthreshold oscillations are frequently observed. Their firing patterns include high-frequency tonic firing and rhythmic or irregular bursting (stuttering). What is the origin of this variability? In the present paper, we hypothesize that it emerges naturally if one assumes a continuous distribution of properties in a small set of active channels. To test this hypothesis, we construct a minimal, single-compartment conductance-based model of FS cells that includes transient Na⁺, delayed-rectifier K⁺, and slowly inactivating d-type K⁺ conductances. The model is analyzed using nonlinear dynamical system theory. For small Na⁺ window current, the neuron exhibits high-frequency tonic firing. At current threshold, the spike response is almost instantaneous for small d-current conductance, g_d, and it is delayed for larger g_d. As g_d further increases, the neuron stutters. Noise substantially reduces the delay duration and induces subthreshold oscillations. In contrast, when the Na⁺ window current is large, the neuron always fires tonically. Near threshold, the firing rates are low, and the delay to firing is only weakly sensitive to noise; subthreshold oscillations are not observed. We propose that the variability in the response of cortical FS neurons is a consequence of heterogeneities in their g_d and in the strength of their Na⁺ window current. We predict the existence of two types of firing patterns in FS neurons, differing in the sensitivity of the delay duration to noise, in the minimal firing rate of the tonic discharge, and in the existence of subthreshold oscillations. We report experimental results from intracellular recordings supporting this prediction.

About 25% of the neurons in the mammalian neocortex are inhibitory, namely reduce the activity of neurons they contact. These inhibitory neurons exhibit diversity of morphological, chemical, and biophysical properties, and their classification has recently been the focus of much debate. Even neurons belonging to a single class of “fast-spiking” (FS) display a large variety of firing patterns in response to standard square current pulses. Previous works proposed that this class is in fact a discrete set of neuronal subtypes with biophysical properties differing in a discontinuous way. In this work, we propose an alternative theory, according to which the biophysical properties of FS neurons are continuously distributed, but distinct firing patterns emerge due to highly nonlinear dynamics of these neurons. We ascertain this theory by exploring with mathematical techniques a biophysically based model of FS neurons. We demonstrate that variable firing responses of cortical FS neurons can be accounted for if one assumes heterogeneity in the strength of some of the ionic conductances underlying neuronal activity. Our theory predicts the existence of two main firing patterns of FS neurons. This prediction is verified by direct recordings in cortical slices.

Synchrony of neuronal oscillations controlled by GABAergic reversal potentials

GABAergic synapse reversal potential is controlled by the concentration of chloride. This concentration can change significantly during development and as a function of neuronal activity. Thus, GABA inhibition can be hyperpolarizing, shunting, or partially depolarizing. Previous results pinpointed the conditions under which hyperpolarizing inhibition (or depolarizing excitation) can lead to synchrony of neural oscillators. Here we examine the role of the GABAergic reversal potential in generation of synchronous oscillations in circuits of neural oscillators. Using weakly coupled oscillator analysis, we show when shunting and partially depolarizing inhibition can produce synchrony, asynchrony, and coexistence of the two. In particular, we show that this depends critically on such factors as the firing rate, the speed of the synapse, spike frequency adaptation, and, most important, the dynamics of spike generation (type I versus type II). We back up our analysis with simulations of small circuits of conductance-based neurons, as well as large-scale networks of neural oscillators. The simulation results are compatible with the analysis: for example, when bistability is predicted analytically, the large-scale network shows clustered states.

Synchronous and asynchronous bursting states: role of intrinsic neural dynamics

Brain signals such as local field potentials often display gamma-band oscillations (30-70 Hz) in a variety of cognitive tasks. These oscillatory activities possibly reflect synchronization of cell assemblies that are engaged in a cognitive function. A type of pyramidal neurons, i.e., chattering neurons, show fast rhythmic bursting (FRB) in the gamma frequency range, and may play an active role in generating the gamma-band oscillations in the cerebral cortex. Our previous phase response analyses have revealed that the synchronization between the coupled bursting neurons significantly depends on the bursting mode that is defined as the number of spikes in each burst. Namely, a network of neurons bursting through a Ca(2+)-dependent mechanism exhibited sharp transitions between synchronous and asynchronous firing states when the neurons exchanged the bursting mode between singlet, doublet and so on. However, whether a broad class of bursting neuron models commonly show such a network behavior remains unclear. Here, we analyze the mechanism underlying this network behavior using a mathematically tractable neuron model. Then we extend our results to a multi-compartment version of the NaP current-based neuron model and prove a similar tight relationship between the bursting mode changes and the network state changes in this model. Thus, the synchronization behavior couples tightly to the bursting mode in a wide class of networks of bursting neurons.

Synchronization of electrically coupled pairs of inhibitory interneurons in neocortex

We performed a systematic analysis of phase locking in pairs of electrically coupled neocortical fast-spiking (FS) and low-threshold-spiking (LTS) interneurons and in a conductance-based model of a pair of FS cells. Phase-response curves (PRCs) were obtained for real interneurons and the model cells. We used PRCs and the theory of weakly coupled oscillators to make predictions about phase-locking characteristics of cell pairs. Phase locking and the robustness of phase-locked states to differences in intrinsic frequencies of cells were directly examined by driving interneuron pairs through a wide range of firing frequencies. Calculations using PRCs accurately predicted that electrical coupling robustly synchronized the firing of interneurons over all frequencies studied (FS, approximately 25-80 Hz; LTS, approximately 10-30 Hz). The synchronizing ability of electrical coupling and the robustness of the phase-locked states were directly dependent on the strength of coupling but not on firing frequency. The FS cell model also predicted the existence of stable antiphase firing at frequencies below approximately 30 Hz, but no evidence for stable antiphase firing was found using the experimentally determined PRCs or in direct measures of phase locking in pairs of interneurons. Despite significant differences in biophysical properties of FS and LTS cells, their phase-locking behavior was remarkably similar. The wide spikes and shallow action potential afterhyperpolarizations of interneurons, compared with the model, prohibited antiphase behavior. Electrical coupling between cortical interneurons of the same type maintained robust synchronous firing of cell pairs for up to approximately 10% heterogeneity in their intrinsic frequencies.

Phase resetting curves and oscillatory stability in interneurons of rat somatosensory cortex

Synchronous oscillations in neural activity are found over wide areas of the cortex. Specific populations of interneurons are believed to play a significant role in generating these synchronized oscillations through mutual synaptic and gap-junctional interactions. Little is known, though, about the mechanism of how oscillations are maintained stably by particular types of interneurons and by their local networks. To obtain more insight into this, we measured membrane-potential responses to small current-pulse perturbations during regular firing, to construct phase resetting curves (PRCs) for three types of interneurons: nonpyramidal regular-spiking (NPRS), low-threshold spiking (LTS), and fast-spiking (FS) cells. Within each cell type, both monophasic and biphasic PRCs were observed, but the proportions and sensitivities to perturbation amplitude were clearly correlated to cell type. We then analyzed the experimentally measured PRCs to predict oscillation stability, or firing reliability, of cells for a complex stochastic input, as occurs in vivo. To do this, we used a method from random dynamical system theory to estimate Lyapunov exponents of the simplified phase model on the circle. The results indicated that LTS and NPRS cells have greater oscillatory stability (are more reliably entrained) in small noisy inputs than FS cells, which is consistent with their distinct types of threshold dynamics.

Low-dimensional maps encoding dynamics in entorhinal cortex and hippocampus

Cells that produce intrinsic theta oscillations often contain the hyperpolarization-activated current I(h). In this article, we use models and dynamic clamp experiments to investigate the synchronization properties of two such cells (stellate cells of the entorhinal cortex and O-LM cells of the hippocampus) in networks with fast-spiking (FS) interneurons. The model we use for stellate cells and O-LM cells is the same, but the stellate cells are excitatory and the O-LM cells are inhibitory, with inhibitory postsynaptic potential considerably longer than those from FS interneurons. We use spike time response curve methods (STRC), expanding that technique to three-cell networks and giving two different ways in which the analysis of the three-cell network reduces to that of a two-cell network. We show that adding FS cells to a network of stellate cells can desynchronize the stellate cells, while adding them to a network of O-LM cells can synchronize the O-LM cells. These synchronization and desynchronization properties critically depend on I(h). The analysis of the deterministic system allows us to understand some effects of noise on the phase relationships in the stellate networks. The dynamic clamp experiments use biophysical stellate cells and in silico FS cells, with connections that mimic excitation or inhibition, the latter with decay times associated with FS cells or O-LM cells. The results obtained in the dynamic clamp experiments are in a good agreement with the analytical framework.

Estimating receptive fields in the presence of spike-time jitter

Neurons in sensory systems are commonly characterized by their receptive fields. These are experimentally often obtained by reverse-correlation analyses, for example, by calculating the spike-triggered average. The reverse-correlation approach, however, generally assumes a fixed temporal relation between spike-generating stimulus features and measured spikes. Temporal jitter of spikes will therefore distort the estimated receptive fields. Here, a novel extension of widely used reverse-correlation techniques (spike-triggered average as well as spike-triggered covariance) is presented that allows accurate measurements of receptive fields even in the presence of considerable spike-time jitter. It is shown that the method correctly recovers the receptive fields from simulated spike trains. When applied to recordings from auditory receptor cells of locusts, a considerable sharpening of receptive fields as compared to standard spike-triggered averages is observed. In addition, the multiple filters that are obtained from a conventional spike-triggered covariance analysis of these data can be collapsed into a single component if spike jitter is accounted for. Finally, it is shown how further effects on spike timing, such as systematic shifts in spike latency, can be included in the approach.

Inhibitory neurons can facilitate rhythmic activity in a neural network

Synchrony emerges in a population of oscillators interacting through in-phase couplings. We addressed a question of whether inhibitory neurons simply hinder the emergence of the synchronous activity among excitatory neurons, or facilitate it. An analysis of a simple phase model revealed that both cases may take place. Numerical simulations of the more realistic models revealed that inhibitory neurons rather facilitate rhythmic activity.

Conservation rules, their breakdown, and optimality in Caenorhabditis sinusoidal locomotion

Undulatory locomotion is common to nematodes as well as to limbless vertebrates, but its control is not understood in spite of the identification of hundred of genes involved in Caenorhabditis elegans locomotion. To reveal the mechanisms of nematode undulatory locomotion, we quantitatively analysed the movement of C. elegans with genetic perturbations to neurons, muscles, and skeleton (cuticle). We also compared locomotion of different Caenorhabditis species. We constructed a theoretical model that combines mechanics and biophysics, and that is constrained by the observations of propulsion and muscular velocities, as well as wavelength and amplitude of undulations. We find that normalized wavelength is a conserved quantity among wild-type C. elegans individuals, across mutants, and across different species. The velocity of forward propulsion scales linearly with the velocity of the muscular wave and the corresponding slope is also a conserved quantity and almost optimal; the exceptions are in some mutants affecting cuticle structure. In theoretical terms, the optimality of the slope is equivalent to the exact balance between muscular and visco-elastic body reaction bending moments. We find that the amplitude and frequency of undulations are inversely correlated and provide a theoretical explanation for this fact. These experimental results are valid both for young adults and for all larval stages of wild-type C. elegans. In particular, during development, the amplitude scales linearly with the wavelength, consistent with our theory. We also investigated the influence of substrate firmness on motion parameters, and found that it does not affect the above invariants. In general, our biomechanical model can explain the observed robustness of the mechanisms controlling nematode undulatory locomotion.

Factors affecting phase synchronization in integrate-and-fire oscillators

Step changes in input current are known to induce partial phase synchrony in ensembles of leaky integrate-and-fire neurons operating in the oscillatory or "regular firing" regime. An analysis of this phenomenon in the absence of noise is presented based on the probability flux within an ensemble of generalized integrate-and-fire neurons. It is shown that the induction of phase synchrony by a step input can be determined by calculating the ratio of the voltage densities obtained from fully desynchronized ensembles firing at the pre and post-step firing rates. In the limit of low noise and in the absence of phase synchrony, the probability density as a function of voltage is inversely proportional to the time derivative along the voltage trajectory. It follows that the magnitude of phase synchronization depends on the degree to which a change in input leads to a uniform multiplication of the voltage derivative over the range from reset to spike threshold. This analysis is used to investigate several factors affecting phase synchronization including high firing rates, inputs modeled as conductances rather than currents, peri-threshold sodium currents, and spike-triggered potassium currents. Finally, we show that without noise, the equilibrium ensemble density is proportional to the phase response curve commonly used to analyze oscillatory systems.

Bidirectional modulation of neuronal responses by depolarizing GABAergic inputs

The reversal potential of GABAA receptor channels is known to be less negative than the resting membrane potential under some cases. Recent electrophysiological experiments revealed that a GABAergic unitary conductance with such a depolarized reversal potential could not only prevent but also facilitate action potential generation depending on the timing of its application relative to the excitatory unitary conductance. Using a two-dimensional point neuron model, we simulate the experiments regarding the integration of unitary conductances, and execute bifurcation analysis. Then we extend our analysis to the case in which the neuron receives two kinds of periodic input trains-an excitatory one and a GABAergic one. We show that the periodic depolarizing GABAergic input train can modulate the output time-averaged firing rate bidirectionally, namely as an increase or a decrease, in a devil's-staircase-like manner depending on the phase difference with the excitatory input train. Bifurcation analysis reveals the existence of a wide variety of phase-locked solutions underlying such a graded response of the neuron. We examine how the input time-width and the value of the GABAA reversal potential affect the response. Moreover, considering a neuronal population, we show that depolarizing GABAergic inputs bidirectionally modulate the amplitude of the oscillatory population activity.

Matching storage and recall: hippocampal spike timing-dependent plasticity and phase response curves

Hippocampal area CA3 is widely considered to function as an autoassociative memory. However, it is insufficiently understood how it does so. In particular, the extensive experimental evidence for the importance of carefully regulated spiking times poses the question as to how spike timing-based dynamics may support memory functions. Here, we develop a normative theory of autoassociative memory encompassing such network dynamics. Our theory specifies the way that the synaptic plasticity rule of a memory constrains the form of neuronal interactions that will retrieve memories optimally. If memories are stored by spike timing-dependent plasticity, neuronal interactions should be formalized in terms of a phase response curve, indicating the effect of presynaptic spikes on the timing of postsynaptic spikes. We show through simulation that such memories are competent analog autoassociators and demonstrate directly that the attributes of phase response curves of CA3 pyramidal cells recorded in vitro qualitatively conform with the theory.