A universal description of stochastic oscillators.

Many systems in physics, chemistry, and biology exhibit oscillations with a pronounced random component. Such stochastic oscillations can emerge via different mechanisms, for example, linear dynamics of a stable focus with fluctuations, limit-cycle systems perturbed by noise, or excitable systems in which random inputs lead to a train of pulses. Despite their diverse origins, the phenomenology of random oscillations can be strikingly similar. Here, we introduce a nonlinear transformation of stochastic oscillators to a complex-valued function [Formula: see text](x) that greatly simplifies and unifies the mathematical description of the oscillator's spontaneous activity, its response to an external time-dependent perturbation, and the correlation statistics of different oscillators that are weakly coupled. The function [Formula: see text] (x) is the eigenfunction of the Kolmogorov backward operator with the least negative (but nonvanishing) eigenvalue λ1 = µ1 + iω1. The resulting power spectrum of the complex-valued function is exactly given by a Lorentz spectrum with peak frequency ω1 and half-width µ1; its susceptibility with respect to a weak external forcing is given by a simple one-pole filter, centered around ω1; and the cross-spectrum between two coupled oscillators can be easily expressed by a combination of the spontaneous power spectra of the uncoupled systems and their susceptibilities. Our approach makes qualitatively different stochastic oscillators comparable, provides simple characteristics for the coherence of the random oscillation, and gives a framework for the description of weakly coupled oscillators.

Spikebench: An open benchmark for spike train time-series classification.

Modern well-performing approaches to neural decoding are based on machine learning models such as decision tree ensembles and deep neural networks. The wide range of algorithms that can be utilized to learn from neural spike trains, which are essentially time-series data, results in the need for diverse and challenging benchmarks for neural decoding, similar to the ones in the fields of computer vision and natural language processing. In this work, we propose a spike train classification benchmark, based on open-access neural activity datasets and consisting of several learning tasks such as stimulus type classification, animal's behavioral state prediction, and neuron type identification. We demonstrate that an approach based on hand-crafted time-series feature engineering establishes a strong baseline performing on par with state-of-the-art deep learning-based models for neural decoding. We release the code allowing to reproduce the reported results.

A framework for macroscopic phase-resetting curves for generalised spiking neural networks.

Brain rhythms emerge from synchronization among interconnected spiking neurons. Key properties of such rhythms can be gleaned from the phase-resetting curve (PRC). Inferring the PRC and developing a systematic phase reduction theory for large-scale brain rhythms remains an outstanding challenge. Here we present a theoretical framework and methodology to compute the PRC of generic spiking networks with emergent collective oscillations. We adopt a renewal approach where neurons are described by the time since their last action potential, a description that can reproduce the dynamical feature of many cell types. For a sufficiently large number of neurons, the network dynamics are well captured by a continuity equation known as the refractory density equation. We develop an adjoint method for this equation giving a semi-analytical expression of the infinitesimal PRC. We confirm the validity of our framework for specific examples of neural networks. Our theoretical framework can link key biological properties at the individual neuron scale and the macroscopic oscillatory network properties. Beyond spiking networks, the approach is applicable to a broad class of systems that can be described by renewal processes.

Efficient and robust coding in heterogeneous recurrent networks.

Cortical networks show a large heterogeneity of neuronal properties. However, traditional coding models have focused on homogeneous populations of excitatory and inhibitory neurons. Here, we analytically derive a class of recurrent networks of spiking neurons that close to optimally track a continuously varying input online, based on two assumptions: 1) every spike is decoded linearly and 2) the network aims to reduce the mean-squared error between the input and the estimate. From this we derive a class of predictive coding networks, that unifies encoding and decoding and in which we can investigate the difference between homogeneous networks and heterogeneous networks, in which each neurons represents different features and has different spike-generating properties. We find that in this framework, 'type 1' and 'type 2' neurons arise naturally and networks consisting of a heterogeneous population of different neuron types are both more efficient and more robust against correlated noise. We make two experimental predictions: 1) we predict that integrators show strong correlations with other integrators and resonators are correlated with resonators, whereas the correlations are much weaker between neurons with different coding properties and 2) that 'type 2' neurons are more coherent with the overall network activity than 'type 1' neurons.

Phase response approaches to neural activity models with distributed delay.

In weakly coupled neural oscillator networks describing brain dynamics, the coupling delay is often distributed. We present a theoretical framework to calculate the phase response curve of distributed-delay induced limit cycles with infinite-dimensional phase space. Extending previous works, in which non-delayed or discrete-delay systems were investigated, we develop analytical results for phase response curves of oscillatory systems with distributed delay using Gaussian and log-normal delay distributions. We determine the scalar product and normalization condition for the linearized adjoint of the system required for the calculation of the phase response curve. As a paradigmatic example, we apply our technique to the Wilson-Cowan oscillator model of excitatory and inhibitory neuronal populations under the two delay distributions. We calculate and compare the phase response curves for the Gaussian and log-normal delay distributions. The phase response curves obtained from our adjoint calculations match those compiled by the direct perturbation method, thereby proving that the theory of weakly coupled oscillators can be applied successfully for distributed-delay-induced limit cycles. We further use the obtained phase response curves to derive phase interaction functions and determine the possible phase locked states of multiple inter-coupled populations to illuminate different synchronization scenarios. In numerical simulations, we show that the coupling delay distribution can impact the stability of the synchronization between inter-coupled gamma-oscillatory networks.

Multistability and evolution of chimera states in a network of type II Morris-Lecar neurons with asymmetrical nonlocal inhibitory connections.

Formation of synchronous activity patterns is an essential property of neuronal networks that has been of central interest to synchronization theory. Chimera states, where both synchronous and asynchronous activities of neurons co-exist in a single network, are particularly poignant examples of such patterns, whose dynamics and multistability may underlie brain function, such as cognitive tasks. However, dynamical mechanisms of coherent state formation in spiking neuronal networks as well as ways to control these states remain unclear. In this paper, we take a step in this direction by considering the evolution of chimera states in a network of class II excitable Morris-Lecar neurons with asymmetrical nonlocal inhibitory connections. Using the adaptive coherence measure, we are able to partition the network parameter space into regions of various collective behaviors (antiphase synchronous clusters, traveling waves, different types of chimera states as well as a spiking death regime) and have shown multistability between the various regimes. We track the evolution of the chimera states as a function of changed key network parameters and found transitions between various types of chimera states. We further find that the network can demonstrate long transients leading to quasi-persistence of activity patterns in the border regions hinting at near-criticality behaviors.

Progressive alignment of inhibitory and excitatory delay may drive a rapid developmental switch in cortical network dynamics.

Nervous system maturation occurs on multiple levels-synaptic, circuit, and network-at divergent timescales. For example, many synaptic properties mature gradually, whereas emergent network dynamics can change abruptly. Here we combine experimental and theoretical approaches to investigate a sudden transition in spontaneous and sensory evoked thalamocortical activity necessary for the development of vision. Inspired by in vivo measurements of timescales and amplitudes of synaptic currents, we extend the Wilson and Cowan model to take into account the relative onset timing and amplitudes of inhibitory and excitatory neural population responses. We study this system as these parameters are varied within amplitudes and timescales consistent with developmental observations to identify the bifurcations of the dynamics that might explain the network behaviors in vivo. Our findings indicate that the inhibitory timing is a critical determinant of thalamocortical activity maturation; a gradual decay of the ratio of inhibitory to excitatory onset time drives the system through a bifurcation that leads to a sudden switch of the network spontaneous activity from high-amplitude oscillations to a nonoscillatory active state. This switch also drives a change from a threshold bursting to linear response to transient stimuli, also consistent with in vivo observation. Thus we show that inhibitory timing is likely critical to the development of network dynamics and may underlie rapid changes in activity without similarly rapid changes in the underlying synaptic and cellular parameters.NEW & NOTEWORTHY Relying on a generalization of the Wilson-Cowan model, which allows a solid analytic foundation for the understanding of the link between maturation of inhibition and network dynamics, we propose a potential explanation for the role of developing excitatory/inhibitory synaptic delays in mediating a sudden switch in thalamocortical visual activity preceding vision onset.

Macroscopic phase resetting-curves determine oscillatory coherence and signal transfer in inter-coupled neural circuits.

Macroscopic oscillations of different brain regions show multiple phase relationships that are persistent across time and have been implicated in routing information. While multiple cellular mechanisms influence the network oscillatory dynamics and structure the macroscopic firing motifs, one of the key questions is to identify the biophysical neuronal and synaptic properties that permit such motifs to arise. A second important issue is how the different neural activity coherence states determine the communication between the neural circuits. Here we analyse the emergence of phase-locking within bidirectionally delayed-coupled spiking circuits in which global gamma band oscillations arise from synaptic coupling among largely excitable neurons. We consider both the interneuronal (ING) and the pyramidal-interneuronal (PING) population gamma rhythms and the inter coupling targeting the pyramidal or the inhibitory neurons. Using a mean-field approach together with an exact reduction method, we reduce each spiking network to a low dimensional nonlinear system and derive the macroscopic phase resetting-curves (mPRCs) that determine how the phase of the global oscillation responds to incoming perturbations. This is made possible by the use of the quadratic integrate-and-fire model together with a Lorentzian distribution of the bias current. Depending on the type of gamma (PING vs. ING), we show that incoming excitatory inputs can either speed up the macroscopic oscillation (phase advance; type I PRC) or induce both a phase advance and a delay (type II PRC). From there we determine the structure of macroscopic coherence states (phase-locking) of two weakly synaptically-coupled networks. To do so we derive a phase equation for the coupled system which links the synaptic mechanisms to the coherence states of the system. We show that a synaptic transmission delay is a necessary condition for symmetry breaking, i.e. a non-symmetric phase lag between the macroscopic oscillations. This potentially provides an explanation to the experimentally observed variety of gamma phase-locking modes. Our analysis further shows that symmetry-broken coherence states can lead to a preferred direction of signal transfer between the oscillatory networks where this directionality also depends on the timing of the signal. Hence we suggest a causal theory for oscillatory modulation of functional connectivity between cortical circuits.

Dynamical ventral tegmental area circuit mechanisms of alcohol-dependent dopamine release.

A large body of data has identified numerous molecular targets through which ethanol (EtOH) acts on brain circuits. Yet how these multiple mechanisms interact to result in dysregulated dopamine (DA) release under the influence of alcohol in vivo remains unclear. In this manuscript, we delineate potential circuit-level mechanisms responsible for EtOH-dependent dysregulation of DA release from the ventral tegmental area (VTA) into its projection areas. For this purpose, we constructed a circuit model of the VTA that integrates realistic Glutamatergic (Glu) inputs and reproduces DA release observed experimentally. We modelled the concentration-dependent effects of EtOH on its principal VTA targets. We calibrated the model to reproduce the inverted U-shape dose dependence of DA neuron activity on EtOH concentration. The model suggests a primary role of EtOH-induced boost in the Ih and AMPA currents in the DA firing-rate/bursting increase. This is counteracted by potentiated GABA transmission that decreases DA neuron activity at higher EtOH concentrations. Thus, the model connects well-established in vitro pharmacological EtOH targets with its in vivo influence on neuronal activity. Furthermore, we predict that increases in VTA activity produced by moderate EtOH doses require partial synchrony and relatively low rates of the Glu afferents. We propose that the increased frequency of transient (phasic) DA peaks evoked by EtOH results from synchronous population bursts in VTA DA neurons. Our model predicts that the impact of acute ETOH on dopamine release is critically shaped by the structure of the cortical inputs to the VTA.

Action in auctions: neural and computational mechanisms of bidding behaviour.

Competition for resources is a fundamental characteristic of evolution. Auctions have been widely used to model competition of individuals for resources, and bidding behaviour plays a major role in social competition. Yet, how humans learn to bid efficiently remains an open question. We used model-based neuroimaging to investigate the neural mechanisms of bidding behaviour under different types of competition. Twenty-seven subjects (nine male) played a prototypical bidding game: a double action, with three "market" types, which differed in the number of competitors. We compared different computational learning models of bidding: directional learning models (DL), where the model bid is "nudged" depending on whether it was accepted or rejected, along with standard reinforcement learning models (RL). We found that DL fit the behaviour best and resulted in higher payoffs. We found the binary learning signal associated with DL to be represented by neural activity in the striatum distinctly posterior to a weaker reward prediction error signal. We posited that DL is an efficient heuristic for valuation when the action (bid) space is continuous. Indeed, we found that the posterior parietal cortex represents the continuous action space of the task, and the frontopolar prefrontal cortex distinguishes among conditions of social competition. Based on our findings, we proposed a conceptual model that accounts for a sequence of processes that are required to perform successful and flexible bidding under different types of competition.

Sensory noise predicts divisive reshaping of receptive fields.

In order to respond reliably to specific features of their environment, sensory neurons need to integrate multiple incoming noisy signals. Crucially, they also need to compete for the interpretation of those signals with other neurons representing similar features. The form that this competition should take depends critically on the noise corrupting these signals. In this study we show that for the type of noise commonly observed in sensory systems, whose variance scales with the mean signal, sensory neurons should selectively divide their input signals by their predictions, suppressing ambiguous cues while amplifying others. Any change in the stimulus context alters which inputs are suppressed, leading to a deep dynamic reshaping of neural receptive fields going far beyond simple surround suppression. Paradoxically, these highly variable receptive fields go alongside and are in fact required for an invariant representation of external sensory features. In addition to offering a normative account of context-dependent changes in sensory responses, perceptual inference in the presence of signal-dependent noise accounts for ubiquitous features of sensory neurons such as divisive normalization, gain control and contrast dependent temporal dynamics.

Reduced Efficacy of the KCC2 Cotransporter Promotes Epileptic Oscillations in a Subiculum Network Model.

Pharmacoresistant epilepsy is a chronic neurological condition in which a basal brain hyperexcitability results in paroxysmal hypersynchronous neuronal discharges. Human temporal lobe epilepsy has been associated with dysfunction or loss of the potassium-chloride cotransporter KCC2 in a subset of pyramidal cells in the subiculum, a key structure generating epileptic activities. KCC2 regulates intraneuronal chloride and extracellular potassium levels by extruding both ions. Absence of effective KCC2 may alter the dynamics of chloride and potassium levels during repeated activation of GABAergic synapses due to interneuron activity. In turn, such GABAergic stress may itself affect Cl- regulation. Such changes in ionic homeostasis may switch GABAergic signaling from inhibitory to excitatory in affected pyramidal cells and also increase neuronal excitability. Possibly these changes contribute to periodic bursting in pyramidal cells, an essential component in the onset of ictal epileptic events. We tested this hypothesis with a computational model of a subicular network with realistic connectivity. The pyramidal cell model explicitly incorporated the cotransporter KCC2 and its effects on the internal/external chloride and potassium levels. Our network model suggested the loss of KCC2 in a critical number of pyramidal cells increased external potassium and intracellular chloride concentrations leading to seizure-like field potential oscillations. These oscillations included transient discharges leading to ictal-like field events with frequency spectra as in vitro Restoration of KCC2 function suppressed seizure activity and thus may present a useful therapeutic option. These simulations therefore suggest that reduced KCC2 cotransporter activity alone may underlie the generation of ictal discharges. SIGNIFICANCE STATEMENT: Ion regulation in the brain is a major determinant of neural excitability. Intracellular chloride in neurons, a partial determinant of the resting potential and the inhibitory reversal potentials, is regulated together with extracellular potassium via kation chloride cotransporters. During temporal lobe epilepsy, the homeostatic regulation of intracellular chloride is impaired in pyramidal cells, yet how this dysregulation may lead to seizures has not been explored. Using a realistic neural network model describing ion mechanisms, we show that chloride homeostasis pathology provokes seizure activity analogous to recordings from epileptogenic brain tissue. We show that there is a critical percentage of pathological cells required for seizure initiation. Our model predicts that restoration of the chloride homeostasis in pyramidal cells could be a viable antiepileptic strategy.

Temporal integration and 1/f power scaling in a circuit model of cerebellar interneurons.

Inhibitory interneurons interconnected via electrical and chemical (GABAA receptor) synapses form extensive circuits in several brain regions. They are thought to be involved in timing and synchronization through fast feedforward control of principal neurons. Theoretical studies have shown, however, that whereas self-inhibition does indeed reduce response duration, lateral inhibition, in contrast, may generate slow response components through a process of gradual disinhibition. Here we simulated a circuit of interneurons (stellate and basket cells) of the molecular layer of the cerebellar cortex and observed circuit time constants that could rise, depending on parameter values, to >1 s. The integration time scaled both with the strength of inhibition, vanishing completely when inhibition was blocked, and with the average connection distance, which determined the balance between lateral and self-inhibition. Electrical synapses could further enhance the integration time by limiting heterogeneity among the interneurons and by introducing a slow capacitive current. The model can explain several observations, such as the slow time course of OFF-beam inhibition, the phase lag of interneurons during vestibular rotation, or the phase lead of Purkinje cells. Interestingly, the interneuron spike trains displayed power that scaled approximately as 1/f at low frequencies. In conclusion, stellate and basket cells in cerebellar cortex, and interneuron circuits in general, may not only provide fast inhibition to principal cells but also act as temporal integrators that build a very short-term memory.NEW & NOTEWORTHY The most common function attributed to inhibitory interneurons is feedforward control of principal neurons. In many brain regions, however, the interneurons are densely interconnected via both chemical and electrical synapses but the function of this coupling is largely unknown. Based on large-scale simulations of an interneuron circuit of cerebellar cortex, we propose that this coupling enhances the integration time constant, and hence the memory trace, of the circuit.

Semiclassical Identification of Periodic Orbits in a Quantum Many-Body System.

While a wealth of results has been obtained for chaos in single-particle quantum systems, much less is known about chaos in quantum many-body systems. We contribute to recent efforts to make a semiclassical analysis of such systems feasible, which is nontrivial due to the exponential proliferation of orbits with increasing particle number. Employing a recently discovered duality relation, we focus on the collective, coherent motion that together with the also present incoherent one typically leads to a mixture of regular and chaotic dynamics. We investigate a kicked spin chain as an example of a presently experimentally and theoretically much studied class of systems.

Inverse Stochastic Resonance in Cerebellar Purkinje Cells.

Purkinje neurons play an important role in cerebellar computation since their axons are the only projection from the cerebellar cortex to deeper cerebellar structures. They have complex internal dynamics, which allow them to fire spontaneously, display bistability, and also to be involved in network phenomena such as high frequency oscillations and travelling waves. Purkinje cells exhibit type II excitability, which can be revealed by a discontinuity in their f-I curves. We show that this excitability mechanism allows Purkinje cells to be efficiently inhibited by noise of a particular variance, a phenomenon known as inverse stochastic resonance (ISR). While ISR has been described in theoretical models of single neurons, here we provide the first experimental evidence for this effect. We find that an adaptive exponential integrate-and-fire model fitted to the basic Purkinje cell characteristics using a modified dynamic IV method displays ISR and bistability between the resting state and a repetitive activity limit cycle. ISR allows the Purkinje cell to operate in different functional regimes: the all-or-none toggle or the linear filter mode, depending on the variance of the synaptic input. We propose that synaptic noise allows Purkinje cells to quickly switch between these functional regimes. Using mutual information analysis, we demonstrate that ISR can lead to a locally optimal information transfer between the input and output spike train of the Purkinje cell. These results provide the first experimental evidence for ISR and suggest a functional role for ISR in cerebellar information processing.

Dopamine Neurons Change the Type of Excitability in Response to Stimuli.

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

Contribution of synchronized GABAergic neurons to dopaminergic neuron firing and bursting.

In the ventral tegmental area (VTA), interactions between dopamine (DA) and γ-aminobutyric acid (GABA) neurons are critical for regulating DA neuron activity and thus DA efflux. To provide a mechanistic explanation of how GABA neurons influence DA neuron firing, we developed a circuit model of the VTA. The model is based on feed-forward inhibition and recreates canonical features of the VTA neurons. Simulations revealed that γ-aminobutyric acid (GABA) receptor (GABAR) stimulation can differentially influence the firing pattern of the DA neuron, depending on the level of synchronization among GABA neurons. Asynchronous activity of GABA neurons provides a constant level of inhibition to the DA neuron and, when removed, produces a classical disinhibition burst. In contrast, when GABA neurons are synchronized by common synaptic input, their influence evokes additional spikes in the DA neuron, resulting in increased measures of firing and bursting. Distinct from previous mechanisms, the increases were not based on lowered firing rate of the GABA neurons or weaker hyperpolarization by the GABAR synaptic current. This phenomenon was induced by GABA-mediated hyperpolarization of the DA neuron that leads to decreases in intracellular calcium (Ca2+) concentration, thus reducing the Ca2+-dependent potassium (K+) current. In this way, the GABA-mediated hyperpolarization replaces Ca2+-dependent K+ current; however, this inhibition is pulsatile, which allows the DA neuron to fire during the rhythmic pauses in inhibition. Our results emphasize the importance of inhibition in the VTA, which has been discussed in many studies, and suggest a novel mechanism whereby computations can occur locally.

Implications of cellular models of dopamine neurons for disease.

This review addresses the present state of single-cell models of the firing pattern of midbrain dopamine neurons and the insights that can be gained from these models into the underlying mechanisms for diseases such as Parkinson's, addiction, and schizophrenia. We will explain the analytical technique of separation of time scales and show how it can produce insights into mechanisms using simplified single-compartment models. We also use morphologically realistic multicompartmental models to address spatially heterogeneous aspects of neural signaling and neural metabolism. Separation of time scale analyses are applied to pacemaking, bursting, and depolarization block in dopamine neurons. Differences in subpopulations with respect to metabolic load are addressed using multicompartmental models.

Flexible frequency control of cortical oscillations enables computations required for working memory.

Cognitive effort leads to a seeming cacophony of brain oscillations. For example, during tasks engaging working memory (WM), specific oscillatory frequency bands modulate in space and time. Despite ample data correlating such modulation to task performance, a mechanistic explanation remains elusive. We propose that flexible control of neural oscillations provides a unified mechanism for the rapid and controlled transitions between the computational operations required by WM. We show in a spiking network model that modulating the input oscillation frequency sets the network in different operating modes: rapid memory access and load is enabled by the beta-gamma oscillations, maintaining a memory while ignoring distractors by the theta, rapid memory clearance by the alpha. The various frequency bands determine the dynamic gating regimes enabling the necessary operations for WM, whose succession explains the need for the complex oscillatory brain dynamics during effortful cognition.

A review of methods for identifying stochastic resonance in simulations of single neuron models.

Stochastic resonance (SR) is said to be observed when the presence of noise in a nonlinear system enables an output signal from the system to better represent some feature of an input signal than it does in the absence of noise. The effect has been observed in models of individual neurons, and in experiments performed on real neural systems. Despite the ubiquity of biophysical sources of stochastic noise in the nervous system, however, it has not yet been established whether neuronal computation mechanisms involved in performance of specific functions such as perception or learning might exploit such noise as an integral component, such that removal of the noise would diminish performance of these functions. In this paper we revisit the methods used to demonstrate stochastic resonance in models of single neurons. This includes a previously unreported observation in a multicompartmental model of a CA1-pyramidal cell. We also discuss, as a contrast to these classical studies, a form of 'stochastic facilitation', known as inverse stochastic resonance. We draw on the reviewed examples to argue why new approaches to studying 'stochastic facilitation' in neural systems need to be developed.