Type I and Type II Neuron Models Are Selectively Driven by Differential Stimulus Features

Reconfigurable Cascaded Thermal Neuristors for Neuromorphic Computing

While the complementary metal-oxide semiconductor (CMOS) technology is the mainstream for the hardware implementation of neural networks, an alternative route is explored based on a new class of spiking oscillators called "thermal neuristors", which operate and interact solely via thermal processes. Utilizing the insulator-to-metal transition (IMT) in vanadium dioxide, a wide variety of reconfigurable electrical dynamics mirroring biological neurons is demonstrated. Notably, inhibitory functionality is achieved just in a single oxide device, and cascaded information flow is realized exclusively through thermal interactions. To elucidate the underlying mechanisms of the neuristors, a detailed theoretical model is developed, which accurately reflects the experimental results. This study establishes the foundation for scalable and energy-efficient thermal neural networks, fostering progress in brain-inspired computing.

Parameter Estimation in the Age of Degeneracy and Unidentifiability

Parameter estimation from observable or experimental data is a crucial stage in any modeling study. Identifiability refers to one’s ability to uniquely estimate the model parameters from the available data. Structural unidentifiability in dynamic models, the opposite of identifiability, is associated with the notion of degeneracy where multiple parameter sets produce the same pattern. Therefore, the inverse function of determining the model parameters from the data is not well defined. Degeneracy is not only a mathematical property of models, but it has also been reported in biological experiments. Classical studies on structural unidentifiability focused on the notion that one can at most identify combinations of unidentifiable model parameters. We have identified a different type of structural degeneracy/unidentifiability present in a family of models, which we refer to as the Lambda-Omega (Λ-Ω) models. These are an extension of the classical lambda-omega (λ-ω) models that have been used to model biological systems, and display a richer dynamic behavior and waveforms that range from sinusoidal to square wave to spike like. We show that the Λ-Ω models feature infinitely many parameter sets that produce identical stable oscillations, except possible for a phase shift (reflecting the initial phase). These degenerate parameters are not identifiable combinations of unidentifiable parameters as is the case in structural degeneracy. In fact, reducing the number of model parameters in the Λ-Ω models is minimal in the sense that each one controls a different aspect of the model dynamics and the dynamic complexity of the system would be reduced by reducing the number of parameters. We argue that the family of Λ-Ω models serves as a framework for the systematic investigation of degeneracy and identifiability in dynamic models and for the investigation of the interplay between structural and other forms of unidentifiability resulting on the lack of information from the experimental/observational data.

Emergence of global synchronization in directed excitatory networks of type I neurons

The collective behaviour of neural networks depends on the cellular and synaptic properties of the neurons. The phase-response curve (PRC) is an experimentally obtainable measure of cellular properties that quantifies the shift in the next spike time of a neuron as a function of the phase at which stimulus is delivered to that neuron. The neuronal PRCs can be classified as having either purely positive values (type I) or distinct positive and negative regions (type II). Networks of type 1 PRCs tend not to synchronize via mutual excitatory synaptic connections. We study the synchronization properties of identical type I and type II neurons, assuming unidirectional synapses. Performing the linear stability analysis and the numerical simulation of the extended Kuramoto model, we show that feedforward loop motifs favour synchronization of type I excitatory and inhibitory neurons, while feedback loop motifs destroy their synchronization tendency. Moreover, large directed networks, either without feedback motifs or with many of them, have been constructed from the same undirected backbones, and a high synchronization level is observed for directed acyclic graphs with type I neurons. It has been shown that, the synchronizability of type I neurons depends on both the directionality of the network connectivity and the topology of its undirected backbone. The abundance of feedforward motifs enhances the synchronizability of the directed acyclic graphs.

The Dynamics of Networks of Identical Theta Neurons

We consider finite and infinite all-to-all coupled networks of identical theta neurons. Two types of synaptic interactions are investigated: instantaneous and delayed (via first-order synaptic processing). Extensive use is made of the Watanabe/Strogatz (WS) ansatz for reducing the dimension of networks of identical sinusoidally-coupled oscillators. As well as the degeneracy associated with the constants of motion of the WS ansatz, we also find continuous families of solutions for instantaneously coupled neurons, resulting from the reversibility of the reduced model and the form of the synaptic input. We also investigate a number of similar related models. We conclude that the dynamics of networks of all-to-all coupled identical neurons can be surprisingly complicated.

Complete Firing-Rate Response of Neurons with Complex Intrinsic Dynamics

The response of a neuronal population over a space of inputs depends on the intrinsic properties of its constituent neurons. Two main modes of single neuron dynamics–integration and resonance–have been distinguished. While resonator cell types exist in a variety of brain areas, few models incorporate this feature and fewer have investigated its effects. To understand better how a resonator’s frequency preference emerges from its intrinsic dynamics and contributes to its local area’s population firing rate dynamics, we analyze the dynamic gain of an analytically solvable two-degree of freedom neuron model. In the Fokker-Planck approach, the dynamic gain is intractable. The alternative Gauss-Rice approach lifts the resetting of the voltage after a spike. This allows us to derive a complete expression for the dynamic gain of a resonator neuron model in terms of a cascade of filters on the input. We find six distinct response types and use them to fully characterize the routes to resonance across all values of the relevant timescales. We find that resonance arises primarily due to slow adaptation with an intrinsic frequency acting to sharpen and adjust the location of the resonant peak. We determine the parameter regions for the existence of an intrinsic frequency and for subthreshold and spiking resonance, finding all possible intersections of the three. The expressions and analysis presented here provide an account of how intrinsic neuron dynamics shape dynamic population response properties and can facilitate the construction of an exact theory of correlations and stability of population activity in networks containing populations of resonator neurons.

Dynamic gain, the amount by which features at specific frequencies in the input to a neuron are amplified or attenuated in its output spiking, is fundamental for the encoding of information by neural populations. Most studies of dynamic gain have focused on neurons without intrinsic degrees of freedom exhibiting integrator-type subthreshold dynamics. Many neuron types in the brain, however, exhibit complex subthreshold dynamics such as resonance, found for instance in cortical interneurons, stellate cells, and mitral cells. A resonator neuron has at least two degrees of freedom for which the classical Fokker-Planck approach to calculating the dynamic gain is largely intractable. Here, we lift the voltage-reset rule after a spike, allowing us to derive a complete expression of the dynamic gain of a resonator neuron model. We find the gain can exhibit only six shapes. The resonant ones have peaks that become large due to intrinsic adaptation and become sharp due to an intrinsic frequency. A resonance can nevertheless result from either property. The analysis presented here helps explain how intrinsic neuron dynamics shape population-level response properties and provides a powerful tool for developing theories of inter-neuron correlations and dynamic responses of neural populations.

Active dendrites mediate stratified gamma-range coincidence detection in hippocampal model neurons

KEY POINTS

Quantitative metrics for the temporal window of integration/coincidence detection, based on the spike-triggered average, were employed to assess the emergence and dependence of gamma-range coincidence detection in hippocampal pyramidal neurons on various ion channel combinations. The presence of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels decreased the coincidence detection window (CDW) of the neuronal compartment to the gamma frequency range. Interaction of HCN channels with T-type calcium channels and persistent sodium channels further reduced the CDW, whereas interaction with A-type potassium channels broadened the CDW. When multiple channel gradients were co-expressed, the high density of resonating conductances in the distal dendrites led to a slow gamma CDW in the proximal dendrites and a fast-gamma CDW in the distal dendrites. The presence of resonating and spike-generating conductances serve as a mechanism underlying the emergence of stratified gamma-range coincidence detection in the dendrites of CA1 pyramidal neurons, enabling them to perform behaviour- and state-dependent gamma frequency multiplexing.

ABSTRACT

Hippocampal pyramidal neurons exhibit gamma-phase preference in their spikes, selectively route inputs through gamma frequency multiplexing and are considered part of gamma-bound cell assemblies. How do these neurons exhibit gamma-frequency coincidence detection capabilities, a feature that is essential for the expression of these physiological observations, despite their slow membrane time constant? In this conductance-based modelling study, we developed quantitative metrics for the temporal window of integration/coincidence detection based on the spike-triggered average (STA) of the neuronal compartment. We employed these metrics in conjunction with quantitative measures for spike initiation dynamics to assess the emergence and dependence of coincidence detection and STA spectral selectivity on various ion channel combinations. We found that the presence of resonating conductances (hyperpolarization-activated cyclic nucleotide-gated or T-type calcium), either independently or synergistically when expressed together, led to the emergence of spectral selectivity in the spike initiation dynamics and a significant reduction in the coincidence detection window (CDW). The presence of A-type potassium channels, along with resonating conductances, reduced the STA characteristic frequency and broadened the CDW, but persistent sodium channels sharpened the CDW by strengthening the spectral selectivity in the STA. Finally, in a morphologically precise model endowed with experimentally constrained channel gradients, we found that somatodendritic compartments expressed functional maps of strong theta-frequency selectivity in spike initiation dynamics and gamma-range CDW. Our results reveal the heavy expression of resonating and spike-generating conductances as the mechanism underlying the robust emergence of stratified gamma-range coincidence detection in the dendrites of hippocampal and cortical pyramidal neurons.

Biophysical Insights into How Spike Threshold Depends on the Rate of Membrane Potential Depolarization in Type I and Type II Neurons

Dynamic spike threshold plays a critical role in neuronal input-output relations. In many neurons, the threshold potential depends on the rate of membrane potential depolarization (dV/dt) preceding a spike. There are two basic classes of neural excitability, i.e., Type I and Type II, according to input-output properties. Although the dynamical and biophysical basis of their spike initiation has been established, the spike threshold dynamic for each cell type has not been well described. Here, we use a biophysical model to investigate how spike threshold depends on dV/dt in two types of neuron. It is observed that Type II spike threshold is more depolarized and more sensitive to dV/dt than Type I. With phase plane analysis, we show that each threshold dynamic arises from the different separatrix and K+ current kinetics. By analyzing subthreshold properties of membrane currents, we find the activation of hyperpolarizing current prior to spike initiation is a major factor that regulates the threshold dynamics. The outward K+ current in Type I neuron does not activate at the perithresholds, which makes its spike threshold insensitive to dV/dt. The Type II K+ current activates prior to spike initiation and there is a large net hyperpolarizing current at the perithresholds, which results in a depolarized threshold as well as a pronounced threshold dynamic. These predictions are further attested in several other functionally equivalent cases of neural excitability. Our study provides a fundamental description about how intrinsic biophysical properties contribute to the threshold dynamics in Type I and Type II neurons, which could decipher their significant functions in neural coding.

Analysis of the role of the low threshold currents IT and Ih in intrinsic delta oscillations of thalamocortical neurons

Thalamocortical neurons are involved in the generation and maintenance of brain rhythms associated with global functional states. The repetitive burst firing of TC neurons at delta frequencies (1–4 Hz) has been linked to the oscillations recorded during deep sleep and during episodes of absence seizures. To get insight into the biophysical properties that are the basis for intrinsic delta oscillations in these neurons, we performed a bifurcation analysis of a minimal conductance-based thalamocortical neuron model including only the I_T channel and the sodium and potassium leak channels. This analysis unveils the dynamics of repetitive burst firing of TC neurons, and describes how the interplay between the amplifying variable m_T and the recovering variable h_T of the calcium channel I_T is sufficient to generate low threshold oscillations in the delta band. We also explored the role of the hyperpolarization activated cationic current I_h in this reduced model and determine that, albeit not required, I_h amplifies and stabilizes the oscillation.

When less is more: non-monotonic spike sequence processing in neurons

Fundamental response properties of neurons centrally underly the computational capabilities of both individual nerve cells and neural networks. Most studies on neuronal input-output relations have focused on continuous-time inputs such as constant or noisy sinusoidal currents. Yet, most neurons communicate via exchanging action potentials (spikes) at discrete times. Here, we systematically analyze the stationary spiking response to regular spiking inputs and reveal that it is generically non-monotonic. Our theoretical analysis shows that the underlying mechanism relies solely on a combination of the discrete nature of the communication by spikes, the capability of locking output to input spikes and limited resources required for spike processing. Numerical simulations of mathematically idealized and biophysically detailed models, as well as neurophysiological experiments confirm and illustrate our theoretical predictions.

Effects of spike-triggered negative feedback on receptive-field properties

Sensory neurons are often described in terms of a receptive field, that is, a linear kernel through which stimuli are filtered before they are further processed. If information transmission is assumed to proceed in a feedforward cascade, the receptive field may be interpreted as the external stimulus' profile maximizing neuronal output. The nervous system, however, contains many feedback loops, and sensory neurons filter more currents than the ones representing the transduced external stimulus. Some of the additional currents are generated by the output activity of the neuron itself, and therefore constitute feedback signals. By means of a time-frequency analysis of the input/output transformation, here we show how feedback modifies the receptive field. The model is applicable to various types of feedback processes, from spike-triggered intrinsic conductances to inhibitory synaptic inputs from nearby neurons. We distinguish between the intrinsic receptive field (filtering all input currents) and the effective receptive field (filtering only external stimuli). Whereas the intrinsic receptive field summarizes the biophysical properties of the neuron associated to subthreshold integration and spike generation, only the effective receptive field can be interpreted as the external stimulus' profile maximizing neuronal output. We demonstrate that spike-triggered feedback shifts low-pass filtering towards band-pass processing, transforming integrator neurons into resonators. For strong feedback, a sharp resonance in the spectral neuronal selectivity may appear. Our results provide a unified framework to interpret a collection of previous experimental studies where specific feedback mechanisms were shown to modify the filtering properties of neurons.

The Competing Benefits of Noise and Heterogeneity in Neural Coding

Noise and heterogeneity are both known to benefit neural coding. Stochastic resonance describes how noise, in the form of random fluctuations in a neuron's membrane voltage, can improve neural representations of an input signal. Neuronal heterogeneity refers to variation in any one of a number of neuron parameters and is also known to increase the information content of a population. We explore the interaction between noise and heterogeneity and find that their benefits to neural coding are not independent. Specifically, a neuronal population better represents an input signal when either noise or heterogeneity is added, but adding both does not always improve representation further. To explain this phenomenon, we propose that noise and heterogeneity operate using two shared mechanisms: (1) temporally desynchronizing the firing of neurons in the population and (2) linearizing the response of a population to a stimulus. We first characterize the effects of noise and heterogeneity on the information content of populations of either leaky integrate-and-fire or FitzHugh-Nagumo neurons. We then examine how the mechanisms of desynchronization and linearization produce these effects, and find that they work to distribute information equally across all neurons in the population in terms of both signal timing (desynchronization) and signal amplitude (linearization). Without noise or heterogeneity, all neurons encode the same aspects of the input signal; adding noise or heterogeneity allows neurons to encode complementary aspects of the input signal, thereby increasing information content. The simulations detailed in this letter highlight the importance of heterogeneity and noise in population coding, demonstrate their complex interactions in terms of the information content of neurons, and explain these effects in terms of underlying mechanisms.

Abrupt and gradual transitions between low and hyperexcited firing frequencies in neuronal models with fast synaptic excitation: a comparative study

Hyperexcitability of neuronal networks is one of the hallmarks of epileptic brain seizure generation, and results from a net imbalance between excitation and inhibition that promotes excessive abnormal firing frequencies. The transition between low and high firing frequencies as the levels of recurrent AMPA excitation change can occur either gradually or abruptly. We used modeling, numerical simulations, and dynamical systems tools to investigate the biophysical and dynamic mechanisms that underlie these two identified modes of transition in recurrently connected neurons via AMPA excitation. We compare our results and demonstrate that these two modes of transition are qualitatively different and can be linked to different intrinsic properties of the participating neurons.

Impact of neuronal properties on network coding: roles of spike initiation dynamics and robust synchrony transfer

Neural networks are more than the sum of their parts, but the properties of those parts are nonetheless important. For instance, neuronal properties affect the degree to which neurons receiving common input will spike synchronously, and whether that synchrony will propagate through the network. Stimulus-evoked synchrony can help or hinder network coding depending on the type of code. In this Perspective, we describe how spike initiation dynamics influence neuronal input-output properties, how those properties affect synchronization, and how synchronization affects network coding. We propose that synchronous and asynchronous spiking can be used to multiplex temporal (synchrony) and rate coding and discuss how pyramidal neurons would be well suited for that task.

Linking dynamical and functional properties of intrinsically bursting neurons

Several studies have shown that bursting neurons can encode information in the number of spikes per burst: As the stimulus varies, so does the length of individual bursts. There presented stimuli, however, vary substantially among different sensory modalities and different neurons.The goal of this paper is to determine which kind of stimulus features can be encoded in burst length, and how those features depend on the mathematical properties of the underlying dynamical system.We show that the initiation and termination of each burst is triggered by specific stimulus features whose temporal characteristsics are determined by the types of bifurcations that initiate and terminate firing in each burst. As only a few bifurcations are possible, only a restricted number of encoded features exists. Here we focus specifically on describing parabolic, square-wave and elliptic bursters. We find that parabolic bursters, whose firing is initiated and terminated by saddle-node bifurcations, behave as prototypical integrators: Firing is triggered by depolarizing stimuli, and lasts for as long as excitation is prolonged. Elliptic bursters, contrastingly, constitute prototypical resonators, since both the initiating and terminating bifurcations possess well-defined oscillation time scales. Firing is therefore triggered by stimulus stretches of matching frequency and terminated by a phase-inversion in the oscillation. The behavior of square-wave bursters is somewhat intermediate, since they are triggered by a fold bifurcation of cycles of well-defined frequency but are terminated by a homoclinic bifurcation lacking an oscillating time scale. These correspondences show that stimulus selectivity is determined by the type of bifurcations. By testing several neuron models, we also demonstrate that additional biological properties that do not modify the bifurcation structure play a minor role in stimulus encoding. Moreover, we show that burst-length variability (and thereby, the capacity to transmit information) depends on a trade-off between the variance of the external signal driving the cell and the strength of the slow internal currents modulating bursts. Thus, our work explicitly links the computational properties of bursting neurons to the mathematical properties of the underlying dynamical systems.

Higher-order spike triggered analysis of neural oscillators

For the purpose of elucidating the neural coding process based on the neural excitability mechanism, researchers have recently investigated the relationship between neural dynamics and the spike triggered stimulus ensemble (STE). Ermentrout et al. analytically derived the relational equation between the phase response curve (PRC) and the spike triggered average (STA). The STA is the first cumulant of the STE. However, in order to understand the neural function as the encoder more explicitly, it is necessary to elucidate the relationship between the PRC and higher-order cumulants of the STE. In this paper, we give a general formulation to relate the PRC and the nth moment of the STE. By using this formulation, we derive a relational equation between the PRC and the spike triggered covariance (STC), which is the covariance of the STE. We show the effectiveness of the relational equation through numerical simulations and use the equation to identify the feature space of the rat hippocampal CA1 pyramidal neurons from their PRCs. Our result suggests that the hippocampal CA1 pyramidal neurons oscillating in the theta frequency range are commonly sensitive to inputs composed of theta and gamma frequency components.

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.

Spike-triggered covariance: geometric proof, symmetry properties, and extension beyond Gaussian stimuli

The space of sensory stimuli is complex and high-dimensional. Yet, single neurons in sensory systems are typically affected by only a small subset of the vast space of all possible stimuli. A proper understanding of the input–output transformation represented by a given cell therefore requires the identification of the subset of stimuli that are relevant in shaping the neuronal response. As an extension to the commonly-used spike-triggered average, the analysis of the spike-triggered covariance matrix provides a systematic methodology to detect relevant stimuli. As originally designed, the consistency of this method is guaranteed only if stimuli are drawn from a Gaussian distribution. Here we present a geometric proof of consistency, which provides insight into the foundations of the method, in particular, into the crucial role played by the geometry of stimulus space and symmetries in the stimulus–response relation. This approach leads to a natural extension of the applicability of the spike-triggered covariance technique to arbitrary spherical or elliptic stimulus distributions. The extension only requires a subtle modification of the original prescription. Furthermore, we present a new resampling method for assessing statistical significance of identified relevant stimuli, applicable to spherical and elliptic stimulus distributions. Finally, we exemplify the modified method and compare it to other prescriptions given in the literature.

A-current and type I/type II transition determine collective spiking from common input

The mechanisms and impact of correlated, or synchronous, firing among pairs and groups of neurons are under intense investigation throughout the nervous system. A ubiquitous circuit feature that can give rise to such correlations consists of overlapping, or common, inputs to pairs and populations of cells, leading to common spike train responses. Here, we use computational tools to study how the transfer of common input currents into common spike outputs is modulated by the physiology of the recipient cells. We focus on a key conductance, g(A), for the A-type potassium current, which drives neurons between "type II" excitability (low g(A)), and "type I" excitability (high g(A)). Regardless of g(A), cells transform common input fluctuations into a tendency to spike nearly simultaneously. However, this process is more pronounced at low g(A) values. Thus, for a given level of common input, type II neurons produce spikes that are relatively more correlated over short time scales. Over long time scales, the trend reverses, with type II neurons producing relatively less correlated spike trains. This is because these cells' increased tendency for simultaneous spiking is balanced by an anticorrelation of spikes at larger time lags. These findings extend and interpret prior findings for phase oscillators to conductance-based neuron models that cover both oscillatory (superthreshold) and subthreshold firing regimes. We demonstrate a novel implication for neural signal processing: downstream cells with long time constants are selectively driven by type I cell populations upstream and those with short time constants by type II cell populations. Our results are established via high-throughput numerical simulations and explained via the cells' filtering properties and nonlinear dynamics.

Membrane properties and the balance between excitation and inhibition control gamma-frequency oscillations arising from feedback inhibition

Computational studies as well as in vivo and in vitro results have shown that many cortical neurons fire in a highly irregular manner and at low average firing rates. These patterns seem to persist even when highly rhythmic signals are recorded by local field potential electrodes or other methods that quantify the summed behavior of a local population. Models of the 30-80 Hz gamma rhythm in which network oscillations arise through 'stochastic synchrony' capture the variability observed in the spike output of single cells while preserving network-level organization. We extend upon these results by constructing model networks constrained by experimental measurements and using them to probe the effect of biophysical parameters on network-level activity. We find in simulations that gamma-frequency oscillations are enabled by a high level of incoherent synaptic conductance input, similar to the barrage of noisy synaptic input that cortical neurons have been shown to receive in vivo. This incoherent synaptic input increases the emergent network frequency by shortening the time scale of the membrane in excitatory neurons and by reducing the temporal separation between excitation and inhibition due to decreased spike latency in inhibitory neurons. These mechanisms are demonstrated in simulations and in vitro current-clamp and dynamic-clamp experiments. Simulation results further indicate that the membrane potential noise amplitude has a large impact on network frequency and that the balance between excitatory and inhibitory currents controls network stability and sensitivity to external inputs.

Establishing Communication between Neuronal Populations through Competitive Entrainment

The role of gamma frequency oscillation in neuronal interaction, and the relationship between oscillation and information transfer between neurons, has been the focus of much recent research. While the biological mechanisms responsible for gamma oscillation and the properties of resulting networks are well studied, the dynamics of changing phase coherence between oscillating neuronal populations are not well understood. To this end we develop a computational model of competitive selection between multiple stimuli, where the selection and transfer of population-encoded information arises from competition between converging stimuli to entrain a target population of neurons. Oscillation is generated by Pyramidal-Interneuronal Network Gamma through the action of recurrent synaptic connections between a locally connected network of excitatory and inhibitory neurons. Competition between stimuli is driven by differences in coherence of oscillation, while transmission of a single selected stimulus is enabled between generating and receiving neurons via Communication-through-Coherence. We explore the effect of varying synaptic parameters on the competitive transmission of stimuli over different neuron models, and identify a continuous region within the parameter space of the recurrent synaptic loop where inhibition-induced oscillation results in entrainment of target neurons. Within this optimal region we find that competition between stimuli of equal coherence results in model output that alternates between representation of the stimuli, in a manner strongly resembling well-known biological phenomena resulting from competitive stimulus selection such as binocular rivalry.

Inferring single neuron properties in conductance based balanced networks

Balanced states in large networks are a usual hypothesis for explaining the variability of neural activity in cortical systems. In this regime the statistics of the inputs is characterized by static and dynamic fluctuations. The dynamic fluctuations have a Gaussian distribution. Such statistics allows to use reverse correlation methods, by recording synaptic inputs and the spike trains of ongoing spontaneous activity without any additional input. By using this method, properties of the single neuron dynamics that are masked by the balanced state can be quantified. To show the feasibility of this approach we apply it to large networks of conductance based neurons. The networks are classified as Type I or Type II according to the bifurcations which neurons of the different populations undergo near the firing onset. We also analyze mixed networks, in which each population has a mixture of different neuronal types. We determine under which conditions the intrinsic noise generated by the network can be used to apply reverse correlation methods. We find that under realistic conditions we can ascertain with low error the types of neurons present in the network. We also find that data from neurons with similar firing rates can be combined to perform covariance analysis. We compare the results of these methods (that do not requite any external input) to the standard procedure (that requires the injection of Gaussian noise into a single neuron). We find a good agreement between the two procedures.

On the role of frustration in excitable systems

We study the role of frustration in excitable systems that allow for oscillations either by construction or in an induced way. We first generalize the notion of frustration to systems whose dynamical equations do not derive from a Hamiltonian. Their couplings can be directed or undirected; they should come in pairs of opposing effects like attractive and repulsive, or activating and repressive, ferromagnetic and antiferromagnetic. As examples we then consider bistable frustrated units as elementary building blocks of our motifs of coupled units. Frustration can be implemented in these systems in various ways: on the level of a single unit via the coupling of a self-loop of positive feedback to a negative feedback loop, on the level of coupled units via the topology or via the type of coupling which may be repressive or activating. In comparison to systems without frustration, we analyze the impact of frustration on the type and number of attractors and observe a considerable enrichment of phase space, ranging from stable fixed-point behavior over different patterns of coexisting options for phase-locked motion to chaotic behavior. In particular we find multistable behavior even for the smallest motifs as long as they are frustrated. Therefore we confirm an enrichment of phase space here for excitable systems with their many applications in biological systems, a phenomenon that is familiar from frustrated spin systems and less known from frustrated phase oscillators. So the enrichment of phase space seems to be a generic effect of frustration in dynamical systems. For a certain range of parameters our systems may be realized in cell tissues. Our results point therefore on a possible generic origin for dynamical behavior that is flexible and functionally stable at the same time, since frustrated systems provide alternative paths for the same set of parameters and at the same "energy costs."

Neurophysiological and computational principles of cortical rhythms in cognition

Synchronous rhythms represent a core mechanism for sculpting temporal coordination of neural activity in the brain-wide network. This review focuses on oscillations in the cerebral cortex that occur during cognition, in alert behaving conditions. Over the last two decades, experimental and modeling work has made great strides in elucidating the detailed cellular and circuit basis of these rhythms, particularly gamma and theta rhythms. The underlying physiological mechanisms are diverse (ranging from resonance and pacemaker properties of single cells to multiple scenarios for population synchronization and wave propagation), but also exhibit unifying principles. A major conceptual advance was the realization that synaptic inhibition plays a fundamental role in rhythmogenesis, either in an interneuronal network or in a reciprocal excitatory-inhibitory loop. Computational functions of synchronous oscillations in cognition are still a matter of debate among systems neuroscientists, in part because the notion of regular oscillation seems to contradict the common observation that spiking discharges of individual neurons in the cortex are highly stochastic and far from being clocklike. However, recent findings have led to a framework that goes beyond the conventional theory of coupled oscillators and reconciles the apparent dichotomy between irregular single neuron activity and field potential oscillations. From this perspective, a plethora of studies will be reviewed on the involvement of long-distance neuronal coherence in cognitive functions such as multisensory integration, working memory, and selective attention. Finally, implications of abnormal neural synchronization are discussed as they relate to mental disorders like schizophrenia and autism.

Switch of encoding characteristics in single neurons by subthreshold and suprathreshold stimuli

Spike-triggered analysis is a statistical method used to elucidate encoding properties in neural systems by estimating the statistical structure of input stimulus preceding spikes. A recent numerical study suggested that the profile of the spike-triggered average (STA) changes depending on whether the mean input stimuli are subthreshold or suprathreshold. Here we analytically verify the difference between subthreshold STA and suprathreshold STA by using the spike response model (SRM). We show by moment expansion that the suprathreshold STA is proportional to the first derivative of the response kernel, and that the subthreshold STA is expressed by a linear combination of the response kernel and its first derivative. We verify whether the analytical results obtained from the SRM can be applied to a multicompartment model with Hodgkin-Huxley type dynamics.

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.

Networks on the edge of chaos: global feedback control of turbulence in oscillator networks

Random networks of coupled phase oscillators with phase shifts in the interaction functions are considered. In such systems, extensive chaos (turbulence) is observed in a wide range of parameters. We show that, by introducing global feedback, the turbulence can be suppressed and a transition to synchronous oscillations can be induced. Our attention is focused on the transition scenario and the properties of patterns, including intermittent turbulence, which are found at the edge of chaos. The emerging coherent patterns represent various self-organized active (sub)networks whose size and behavior can be controlled.