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

TMS-induced phase resets depend on TMS intensity and EEG phase

Objective. The phase of the electroencephalographic (EEG) signal predicts performance in motor, somatosensory, and cognitive functions. Studies suggest that brain phase resets align neural oscillations with external stimuli, or couple oscillations across frequency bands and brain regions. Transcranial Magnetic Stimulation (TMS) can cause phase resets noninvasively in the cortex, thus providing the potential to control phase-sensitive cognitive functions. However, the relationship between TMS parameters and phase resetting is not fully understood. This is especially true of TMS intensity, which may be crucial to enabling precise control over the amount of phase resetting that is induced. Additionally, TMS phase resetting may interact with the instantaneous phase of the brain. Understanding these relationships is crucial to the development of more powerful and controllable stimulation protocols.Approach.To test these relationships, we conducted a TMS-EEG study. We applied single-pulse TMS at varying degrees of stimulation intensity to the motor area in an open loop. Offline, we used an autoregressive algorithm to estimate the phase of the intrinsicµ-Alpha rhythm of the motor cortex at the moment each TMS pulse was delivered.Main results. We identified post-stimulation epochs whereµ-Alpha phase resetting and N100 amplitude depend parametrically on TMS intensity and are significantversusperipheral auditory sham stimulation. We observedµ-Alpha phase inversion after stimulations near peaks but not troughs in the endogenousµ-Alpha rhythm.Significance. These data suggest that low-intensity TMS primarily resets existing oscillations, while at higher intensities TMS may activate previously silent neurons, but only when endogenous oscillations are near the peak phase. These data can guide future studies that seek to induce phase resetting, and point to a way to manipulate the phase resetting effect of TMS by varying only the timing of the pulse with respect to ongoing brain activity.

Effects of Neuromodulation on Excitatory-Inhibitory Neural Network Dynamics Depend on Network Connectivity Structure

Acetylcholine (ACh), one of the brain's most potent neuromodulators, can affect intrinsic neuron properties through blockade of an M-type potassium current. The effect of ACh on excitatory and inhibitory cells with this potassium channel modulates their membrane excitability, which in turn affects their tendency to synchronize in networks. Here, we study the resulting changes in dynamics in networks with inter-connected excitatory and inhibitory populations (E-I networks), which are ubiquitous in the brain. Utilizing biophysical models of E-I networks, we analyze how the network connectivity structure in terms of synaptic connectivity alters the influence of ACh on the generation of synchronous excitatory bursting. We investigate networks containing all combinations of excitatory and inhibitory cells with high (Type I properties) or low (Type II properties) modulatory tone. To vary network connectivity structure, we focus on the effects of the strengths of inter-connections between excitatory and inhibitory cells (E-I synapses and I-E synapses), and the strengths of intra-connections among excitatory cells (E-E synapses) and among inhibitory cells (I-I synapses). We show that the presence of ACh may or may not affect the generation of network synchrony depending on the network connectivity. Specifically, strong network inter-connectivity induces synchronous excitatory bursting regardless of the cellular propensity for synchronization, which aligns with predictions of the PING model. However, when a network's intra-connectivity dominates its inter-connectivity, the propensity for synchrony of either inhibitory or excitatory cells can determine the generation of network-wide bursting.

Perampanel reduces paroxysmal depolarizing shift and inhibitory synaptic input in excitatory neurons to inhibit epileptic network oscillations

BACKGROUND AND PURPOSE

Perampanel is a newly approved anticonvulsant uniquely targeting AMPA receptors, which mediate the most abundant form of excitatory synaptic transmission in the brain. However, the network mechanism underlying the anti-epileptic effect of the AMPAergic inhibition remains to be explored.

EXPERIMENTAL APPROACH

The mechanism of perampanel action was studied with the basolateral amygdala network containing pyramidal-inhibitory neuronal resonators in seizure models of 4-aminopyridine (4-AP) and electrical kindling.

KEY RESULTS

Application of either 4-AP or electrical kindling to the basolateral amygdala readily induces AMPAergic transmission-dependent reverberating activities between pyramidal-inhibitory neuronal resonators, which are chiefly characterized by burst discharges in inhibitory neurons and corresponding recurrent inhibitory postsynaptic potentials in pyramidal neurons. Perampanel reduces post-kindling "paroxysmal depolarizing shift" especially in pyramidal neurons and, counterintuitively, eliminates burst activities in inhibitory neurons and inhibitory synaptic inputs onto excitatory pyramidal neurons to result in prevention of epileptiform discharges and seizure behaviours. Intriguingly, similar effects can be obtained with not only the AMPA receptor antagonist CNQX but also the GABA_A receptor antagonist bicuculline, which is usually considered as a proconvulsant.

CONCLUSION AND IMPLICATIONS

Ictogenesis depends on the AMPA receptor-dependent recruitment of pyramidal-inhibitory neuronal network oscillations tuned by dynamic glutamatergic and GABAergic transmission. The anticonvulsant effect of perampanel then stems from disruption of the coordinated network activities rather than simply decreased neuronal excitability or excitatory transmission. Positive or negative modulation of epileptic network reverberations may be pro-ictogenic or anti-ictogenic, respectively, constituting a more applicable rationale for the therapy against seizures.

Bridging the Molecular-Cellular Gap in Understanding Ion Channel Clustering

The clustering of many voltage-dependent ion channel molecules at unique neuronal membrane sites such as axon initial segments, nodes of Ranvier, or the post-synaptic density, is an active process mediated by the interaction of ion channels with scaffold proteins and is of immense importance for electrical signaling. Growing evidence indicates that the density of ion channels at such membrane sites may affect action potential conduction properties and synaptic transmission. However, despite the emerging importance of ion channel density for electrical signaling, how ion channel-scaffold protein molecular interactions lead to cellular ion channel clustering, and how this process is regulated are largely unknown. In this review, we emphasize that voltage-dependent ion channel density at native clustering sites not only affects the density of ionic current fluxes but may also affect the conduction properties of the channel and/or the physical properties of the membrane at such locations, all changes that are expected to affect action potential conduction properties. Using the concrete example of the prototypical Shaker voltage-activated potassium channel (Kv) protein, we demonstrate how insight into the regulation of cellular ion channel clustering can be obtained when the molecular mechanism of ion channel-scaffold protein interaction is known. Our review emphasizes that such mechanistic knowledge is essential, and when combined with super-resolution imaging microscopy, can serve to bridge the molecular-cellular gap in understanding the regulation of ion channel clustering. Pressing questions, challenges and future directions in addressing ion channel clustering and its regulation are discussed.

Analyzing the competition of gamma rhythms with delayed pulse-coupled oscillators in phase representation

Networks of neurons can generate oscillatory activity as result of various types of coupling that lead to synchronization. A prominent type of oscillatory activity is gamma (30-80 Hz) rhythms, which may play an important role in neuronal information processing. Two mechanisms have mainly been proposed for their generation: (1) interneuron network gamma (ING) and (2) pyramidal-interneuron network gamma (PING). In vitro and in vivo experiments have shown that both mechanisms can exist in the same cortical circuits. This raises the questions: How do ING and PING interact when both can in principle occur? Are the network dynamics a superposition, or do ING and PING interact in a nonlinear way and if so, how? In this article, we first generalize the phase representation for nonlinear one-dimensional pulse coupled oscillators as introduced by Mirollo and Strogatz to type II oscillators whose phase response curve (PRC) has zero crossings. We then give a full theoretical analysis for the regular gamma-like oscillations of simple networks consisting of two neural oscillators, an "E neuron" mimicking a synchronized group of pyramidal cells, and an "I neuron" representing such a group of interneurons. Motivated by experimental findings, we choose the E neuron to have a type I PRC [leaky integrate-and-fire (LIF) neuron], while the I neuron has either a type I or type II PRC (LIF or "sine" neuron). The phase representation allows us to define in a simple manner scenarios of interaction between the two neurons, which are independent of the types and the details of the neuron models. The presence of delay in the couplings leads to an increased number of scenarios relevant for gamma-like oscillatory patterns. We analytically derive the set of such scenarios and describe their occurrence in terms of parameter values such as synaptic connectivity and drive to the E and I neurons. The networks can be tuned to oscillate in an ING or PING mode. We focus particularly on the transition region where both rhythms compete to govern the network dynamics and compare with oscillations in reduced networks, which can only generate either ING or PING. Our analytically derived oscillation frequency diagrams indicate that except for small coexistence regions, the networks generate ING if the oscillation frequency of the reduced ING network exceeds that of the reduced PING network, and vice versa. For networks with the LIF I neuron, the network oscillation frequency slightly exceeds the frequencies of corresponding reduced networks, while it lies between them for networks with the sine I neuron. In networks oscillating in ING (PING) mode, the oscillation frequency responds faster to changes in the drive to the I (E) neuron than to changes in the drive to the E (I) neuron. This finding suggests a method to analyze which mechanism governs an observed network oscillation. Notably, also when the network operates in ING mode, the E neuron can spike before the I neuron such that relative spike times of the pyramidal cells and the interneurons alone are not conclusive for distinguishing ING and PING.

Firing rate equations require a spike synchrony mechanism to correctly describe fast oscillations in inhibitory networks

Recurrently coupled networks of inhibitory neurons robustly generate oscillations in the gamma band. Nonetheless, the corresponding Wilson-Cowan type firing rate equation for such an inhibitory population does not generate such oscillations without an explicit time delay. We show that this discrepancy is due to a voltage-dependent spike-synchronization mechanism inherent in networks of spiking neurons which is not captured by standard firing rate equations. Here we investigate an exact low-dimensional description for a network of heterogeneous canonical Class 1 inhibitory neurons which includes the sub-threshold dynamics crucial for generating synchronous states. In the limit of slow synaptic kinetics the spike-synchrony mechanism is suppressed and the standard Wilson-Cowan equations are formally recovered as long as external inputs are also slow. However, even in this limit synchronous spiking can be elicited by inputs which fluctuate on a time-scale of the membrane time-constant of the neurons. Our meanfield equations therefore represent an extension of the standard Wilson-Cowan equations in which spike synchrony is also correctly described.

Population models describing the average activity of large neuronal ensembles are a powerful mathematical tool to investigate the principles underlying cooperative function of large neuronal systems. However, these models do not properly describe the phenomenon of spike synchrony in networks of neurons. In particular, they fail to capture the onset of synchronous oscillations in networks of inhibitory neurons. We show that this limitation is due to a voltage-dependent synchronization mechanism which is naturally present in spiking neuron models but not captured by traditional firing rate equations. Here we investigate a novel set of macroscopic equations which incorporate both firing rate and membrane potential dynamics, and that correctly generate fast inhibition-based synchronous oscillations. In the limit of slow-synaptic processing oscillations are suppressed, and the model reduces to an equation formally equivalent to the Wilson-Cowan model.

Dichotomous Dynamics in E-I Networks with Strongly and Weakly Intra-connected Inhibitory Neurons

The interconnectivity between excitatory and inhibitory neural networks informs mechanisms by which rhythmic bursts of excitatory activity can be produced in the brain. One such mechanism, Pyramidal Interneuron Network Gamma (PING), relies primarily upon reciprocal connectivity between the excitatory and inhibitory networks, while also including intra-connectivity of inhibitory cells. The causal relationship between excitatory activity and the subsequent burst of inhibitory activity is of paramount importance to the mechanism and has been well studied. However, the role of the intra-connectivity of the inhibitory network, while important for PING, has not been studied in detail, as most analyses of PING simply assume that inhibitory intra-connectivity is strong enough to suppress subsequent firing following the initial inhibitory burst. In this paper we investigate the role that the strength of inhibitory intra-connectivity plays in determining the dynamics of PING-style networks. We show that networks with weak inhibitory intra-connectivity exhibit variations in burst dynamics of both the excitatory and inhibitory cells that are not obtained with strong inhibitory intra-connectivity. Networks with weak inhibitory intra-connectivity exhibit excitatory rhythmic bursts with weak excitatory-to-inhibitory synapses for which classical PING networks would show no rhythmic activity. Additionally, variations in dynamics of these networks as the excitatory-to-inhibitory synaptic weight increases illustrates the important role that consistent pattern formation in the inhibitory cells serves in maintaining organized and periodic excitatory bursts. Finally, motivated by these results and the known diversity of interneurons, we show that a PING-style network with two inhibitory subnetworks, one strongly intra-connected and one weakly intra-connected, exhibits organized and periodic excitatory activity over a larger parameter regime than networks with a homogeneous inhibitory population. Taken together, these results serve to better articulate the role of inhibitory intra-connectivity in generating PING-like rhythms, while also revealing how heterogeneity amongst inhibitory synapses might make such rhythms more robust to a variety of network parameters.

Salicylate-Induced Suppression of Electrically Driven Activity in Brain Slices from the Auditory Cortex of Aging Mice

The prevalence of tinnitus is known to increase with age. The age-dependent mechanisms of tinnitus may have important implications for the development of new therapeutic treatments. High doses of salicylate can be used experimentally to induce transient tinnitus and hearing loss. Although accumulating evidence indicates that salicylate induces tinnitus by directly targeting neurons in the peripheral and central auditory systems, the precise effect of salicylate on neural networks in the auditory cortex (AC) is unknown. Here, we examined salicylate-induced changes in stimulus-driven laminar responses of AC slices with salicylate superfusion in young and aged senescence-accelerated-prone (SAMP) and -resistant (SAMR) mice. Of the two strains, SAMP1 is known to be a more suitable model of presbycusis. We recorded stimulus-driven laminar local field potential (LFP) responses at multi sites in AC slice preparations. We found that for all AC slices in the two strains, salicylate always reduced stimulus-driven LFP responses in all layers. However, for the amplitudes of the LFP responses, the two senescence-accelerated mice (SAM) strains showed different laminar properties between the pre- and post-salicylate conditions, reflecting strain-related differences in local circuits. As for the relationships between auditory brainstem response (ABR) thresholds and the LFP amplitude ratios in the pre- vs. post-salicylate condition, we found negative correlations in layers 2/3 and 4 for both older strains, and in layer 5 (L5) in older SAMR1. In contrast, the GABAergic agonist muscimol (MSC) led to positive correlations between ABR thresholds and LFP amplitude ratios in the pre- vs. post-MSC condition in younger SAM mice from both strains. Further, in younger mice, salicylate decreased the firing rate in AC L4 pyramidal neurons. Thus, salicylate can directly reduce neural excitability of L4 pyramidal neurons and thereby influence AC neural circuit activity. That we observed age-dependent effects of salicylate and varied GABAergic sensitivity in the AC among mouse strains with hearing loss implies that potential therapeutic mechanisms for tinnitus may operate differently in young vs. aged subjects. Therefore, scientists developing new therapeutic modalities for tinnitus treatment should consider using both aged and young animals.

Synchronization scenarios in the Winfree model of coupled oscillators

Fifty years ago Arthur Winfree proposed a deeply influential mean-field model for the collective synchronization of large populations of phase oscillators. Here we provide a detailed analysis of the model for some special, analytically tractable cases. Adopting the thermodynamic limit, we derive an ordinary differential equation that exactly describes the temporal evolution of the macroscopic variables in the Ott-Antonsen invariant manifold. The low-dimensional model is then thoroughly investigated for a variety of pulse types and sinusoidal phase response curves (PRCs). Two structurally different synchronization scenarios are found, which are linked via the mutation of a Bogdanov-Takens point. From our results, we infer a general rule of thumb relating pulse shape and PRC offset with each scenario. Finally, we compare the exact synchronization threshold with the prediction of the averaging approximation given by the Kuramoto-Sakaguchi model. At the leading order, the discrepancy appears to behave as an odd function of the PRC offset.

Effect of Phase Response Curve Shape and Synaptic Driving Force on Synchronization of Coupled Neuronal Oscillators

The role of the phase response curve (PRC) shape on the synchrony of synaptically coupled oscillating neurons is examined. If the PRC is independent of the phase, because of the synaptic form of the coupling, synchrony is found to be stable for both excitatory and inhibitory coupling at all rates, whereas the antisynchrony becomes stable at low rates. A faster synaptic rise helps extend the stability of antisynchrony to higher rates. If the PRC is not constant but has a profile like that of a leaky integrate-and-fire model, then, in contrast to the earlier reports that did not include the voltage effects, mutual excitation could lead to stable synchrony provided the synaptic reversal potential is below the voltage level the neuron would have reached in the absence of the interaction and threshold reset. This level is controlled by the applied current and the leakage parameters. Such synchrony is contingent on significant phase response (that would result, for example, by a sharp PRC jump) occurring during the synaptic rising phase. The rising phase, however, does not contribute significantly if it occurs before the voltage spike reaches its peak. Then a stable near-synchronous state can still exist between type 1 PRC neurons if the PRC shows a left skewness in its shape. These results are examined comprehensively using perfect integrate-and-fire, leaky integrate-and-fire, and skewed PRC shapes under the assumption of the weakly coupled oscillator theory applied to synaptically coupled neuron models.

A Consistent Definition of Phase Resetting Using Hilbert Transform

A phase resetting curve (PRC) measures the transient change in the phase of a neural oscillator subject to an external perturbation. The PRC encapsulates the dynamical response of a neural oscillator and, as a result, it is often used for predicting phase-locked modes in neural networks. While phase is a fundamental concept, it has multiple definitions that may lead to contradictory results. We used the Hilbert Transform (HT) to define the phase of the membrane potential oscillations and HT amplitude to estimate the PRC of a single neural oscillator. We found that HT's amplitude and its corresponding instantaneous frequency are very sensitive to membrane potential perturbations. We also found that the phase shift of HT amplitude between the pre- and poststimulus cycles gives an accurate estimate of the PRC. Moreover, HT phase does not suffer from the shortcomings of voltage threshold or isochrone methods and, as a result, gives accurate and reliable estimations of phase resetting.

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.

Mathematical Frameworks for Oscillatory Network Dynamics in Neuroscience

The tools of weakly coupled phase oscillator theory have had a profound impact on the neuroscience community, providing insight into a variety of network behaviours ranging from central pattern generation to synchronisation, as well as predicting novel network states such as chimeras. However, there are many instances where this theory is expected to break down, say in the presence of strong coupling, or must be carefully interpreted, as in the presence of stochastic forcing. There are also surprises in the dynamical complexity of the attractors that can robustly appear-for example, heteroclinic network attractors. In this review we present a set of mathematical tools that are suitable for addressing the dynamics of oscillatory neural networks, broadening from a standard phase oscillator perspective to provide a practical framework for further successful applications of mathematics to understanding network dynamics in neuroscience.

Intrinsic Cellular Properties and Connectivity Density Determine Variable Clustering Patterns in Randomly Connected Inhibitory Neural Networks

The plethora of inhibitory interneurons in the hippocampus and cortex play a pivotal role in generating rhythmic activity by clustering and synchronizing cell firing. Results of our simulations demonstrate that both the intrinsic cellular properties of neurons and the degree of network connectivity affect the characteristics of clustered dynamics exhibited in randomly connected, heterogeneous inhibitory networks. We quantify intrinsic cellular properties by the neuron's current-frequency relation (IF curve) and Phase Response Curve (PRC), a measure of how perturbations given at various phases of a neurons firing cycle affect subsequent spike timing. We analyze network bursting properties of networks of neurons with Type I or Type II properties in both excitability and PRC profile; Type I PRCs strictly show phase advances and IF curves that exhibit frequencies arbitrarily close to zero at firing threshold while Type II PRCs display both phase advances and delays and IF curves that have a non-zero frequency at threshold. Type II neurons whose properties arise with or without an M-type adaptation current are considered. We analyze network dynamics under different levels of cellular heterogeneity and as intrinsic cellular firing frequency and the time scale of decay of synaptic inhibition are varied. Many of the dynamics exhibited by these networks diverge from the predictions of the interneuron network gamma (ING) mechanism, as well as from results in all-to-all connected networks. Our results show that randomly connected networks of Type I neurons synchronize into a single cluster of active neurons while networks of Type II neurons organize into two mutually exclusive clusters segregated by the cells' intrinsic firing frequencies. Networks of Type II neurons containing the adaptation current behave similarly to networks of either Type I or Type II neurons depending on network parameters; however, the adaptation current creates differences in the cluster dynamics compared to those in networks of Type I or Type II neurons. To understand these results, we compute neuronal PRCs calculated with a perturbation matching the profile of the synaptic current in our networks. Differences in profiles of these PRCs across the different neuron types reveal mechanisms underlying the divergent network dynamics.

Neurons as oscillators

Regularly spiking neurons can be described as oscillators. In this article we review some of the insights gained from this conceptualization and their relevance for systems neuroscience. First, we explain how a regularly spiking neuron can be viewed as an oscillator and how the phase-response curve (PRC) describes the response of the neuron's spike times to small perturbations. We then discuss the meaning of the PRC for a single neuron's spiking behavior and review the PRCs measured from a variety of neurons in a range of spiking regimes. Next, we show how the PRC can be related to a number of common measures used to quantify neuronal firing, such as the spike-triggered average and the peristimulus histogram. We further show that the response of a neuron to correlated inputs depends on the shape of the PRC. We then explain how the PRC of single neurons can be used to predict neural network behavior. Given the PRC, conduction delays, and the waveform and time course of the synaptic potentials, it is possible to predict neural population behavior such as synchronization. The PRC also allows us to quantify the robustness of the synchronization to heterogeneity and noise. We finally ask how to combine the measured PRCs and the predictions based on PRC to further the understanding of systems neuroscience. As an example, we discuss how the change of the PRC by the neuromodulator acetylcholine could lead to a destabilization of cortical network dynamics. Although all of these studies are grounded in mathematical abstractions that do not strictly hold in biology, they provide good estimates for the emergence of the brain's network activity from the properties of individual neurons. The study of neurons as oscillators can provide testable hypotheses and mechanistic explanations for systems neuroscience.

Phase response curves for models of earthquake fault dynamics

We systematically study effects of external perturbations on models describing earthquake fault dynamics. The latter are based on the framework of the Burridge-Knopoff spring-block system, including the cases of a simple mono-block fault, as well as the paradigmatic complex faults made up of two identical or distinct blocks. The blocks exhibit relaxation oscillations, which are representative for the stick-slip behavior typical for earthquake dynamics. Our analysis is carried out by determining the phase response curves of first and second order. For a mono-block fault, we consider the impact of a single and two successive pulse perturbations, further demonstrating how the profile of phase response curves depends on the fault parameters. For a homogeneous two-block fault, our focus is on the scenario where each of the blocks is influenced by a single pulse, whereas for heterogeneous faults, we analyze how the response of the system depends on whether the stimulus is applied to the block having a shorter or a longer oscillation period.

Cooperation and competition of gamma oscillation mechanisms

Oscillations of neuronal activity in different frequency ranges are thought to reflect important aspects of cortical network dynamics. Here we investigate how various mechanisms that contribute to oscillations in neuronal networks may interact. We focus on networks with inhibitory, excitatory, and electrical synapses, where the subnetwork of inhibitory interneurons alone can generate interneuron gamma (ING) oscillations and the interactions between interneurons and pyramidal cells allow for pyramidal-interneuron gamma (PING) oscillations. What type of oscillation will such a network generate? We find that ING and PING oscillations compete: The mechanism generating the higher oscillation frequency "wins"; it determines the frequency of the network oscillation and suppresses the other mechanism. For type I interneurons, the network oscillation frequency is equal to or slightly above the higher of the ING and PING frequencies in corresponding reduced networks that can generate only either of them; if the interneurons belong to the type II class, it is in between. In contrast to ING and PING, oscillations mediated by gap junctions and oscillations mediated by inhibitory synapses may cooperate or compete, depending on the type (I or II) of interneurons and the strengths of the electrical and chemical synapses. We support our computer simulations by a theoretical model that allows a full theoretical analysis of the main results. Our study suggests experimental approaches to deciding to what extent oscillatory activity in networks of interacting excitatory and inhibitory neurons is dominated by ING or PING oscillations and of which class the participating interneurons are.

Modeling the Influence of Ion Channels on Neuron Dynamics in Drosophila

Voltage gated ion channels play a major role in determining a neuron's firing behavior, resulting in the specific processing of synaptic input patterns. Drosophila and other invertebrates provide valuable model systems for investigating ion channel kinetics and their impact on firing properties. Despite the increasing importance of Drosophila as a model system, few computational models of its ion channel kinetics have been developed. In this study, experimentally observed biophysical properties of voltage gated ion channels from the fruitfly Drosophila melanogaster are used to develop a minimal, conductance based neuron model. We investigate the impact of the densities of these channels on the excitability of the model neuron. Changing the channel densities reproduces different in situ observed firing patterns and induces a switch from integrator to resonator properties. Further, we analyze the preference to input frequency and how it depends on the channel densities and the resulting bifurcation type the system undergoes. An extension to a three dimensional model demonstrates that the inactivation kinetics of the sodium channels play an important role, allowing for firing patterns with a delayed first spike and subsequent high frequency firing as often observed in invertebrates, without altering the kinetics of the delayed rectifier current.

Dynamics of on-off neural firing patterns and stochastic effects near a sub-critical Hopf bifurcation

On-off firing patterns, in which repetition of clusters of spikes are interspersed with epochs of subthreshold oscillations or quiescent states, have been observed in various nervous systems, but the dynamics of this event remain unclear. Here, we report that on-off firing patterns observed in three experimental models (rat sciatic nerve subject to chronic constrictive injury, rat CA1 pyramidal neuron, and rabbit blood pressure baroreceptor) appeared as an alternation between quiescent state and burst containing multiple period-1 spikes over time. Burst and quiescent state had various durations. The interspike interval (ISI) series of on-off firing pattern was suggested as stochastic using nonlinear prediction and autocorrelation function. The resting state was changed to a period-1 firing pattern via on-off firing pattern as the potassium concentration, static pressure, or depolarization current was changed. During the changing process, the burst duration of on-off firing pattern increased and the duration of the quiescent state decreased. Bistability of a limit cycle corresponding to period-1 firing and a focus corresponding to resting state was simulated near a sub-critical Hopf bifurcation point in the deterministic Morris-Lecar (ML) model. In the stochastic ML model, noise-induced transitions between the coexisting regimes formed an on-off firing pattern, which closely matched that observed in the experiment. In addition, noise-induced exponential change in the escape rate from the focus, and noise-induced coherence resonance were identified. The distinctions between the on-off firing pattern and stochastic firing patterns generated near three other types of bifurcations of equilibrium points, as well as other viewpoints on the dynamics of on-off firing pattern, are discussed. The results not only identify the on-off firing pattern as noise-induced stochastic firing pattern near a sub-critical Hopf bifurcation point, but also offer practical indicators to discriminate bifurcation types and neural excitability types.

On the firing rate dependency of the phase response curve of rat Purkinje neurons in vitro

Synchronous spiking during cerebellar tasks has been observed across Purkinje cells: however, little is known about the intrinsic cellular mechanisms responsible for its initiation, cessation and stability. The Phase Response Curve (PRC), a simple input-output characterization of single cells, can provide insights into individual and collective properties of neurons and networks, by quantifying the impact of an infinitesimal depolarizing current pulse on the time of occurrence of subsequent action potentials, while a neuron is firing tonically. Recently, the PRC theory applied to cerebellar Purkinje cells revealed that these behave as phase-independent integrators at low firing rates, and switch to a phase-dependent mode at high rates. Given the implications for computation and information processing in the cerebellum and the possible role of synchrony in the communication with its post-synaptic targets, we further explored the firing rate dependency of the PRC in Purkinje cells. We isolated key factors for the experimental estimation of the PRC and developed a closed-loop approach to reliably compute the PRC across diverse firing rates in the same cell. Our results show unambiguously that the PRC of individual Purkinje cells is firing rate dependent and that it smoothly transitions from phase independent integrator to a phase dependent mode. Using computational models we show that neither channel noise nor a realistic cell morphology are responsible for the rate dependent shift in the phase response curve.

Critical slowing down governs the transition to neuron spiking

Many complex systems have been found to exhibit critical transitions, or so-called tipping points, which are sudden changes to a qualitatively different system state. These changes can profoundly impact the functioning of a system ranging from controlled state switching to a catastrophic break-down; signals that predict critical transitions are therefore highly desirable. To this end, research efforts have focused on utilizing qualitative changes in markers related to a system’s tendency to recover more slowly from a perturbation the closer it gets to the transition—a phenomenon called critical slowing down. The recently studied scaling of critical slowing down offers a refined path to understand critical transitions: to identify the transition mechanism and improve transition prediction using scaling laws.

Here, we outline and apply this strategy for the first time in a real-world system by studying the transition to spiking in neurons of the mammalian cortex. The dynamical system approach has identified two robust mechanisms for the transition from subthreshold activity to spiking, saddle-node and Hopf bifurcation. Although theory provides precise predictions on signatures of critical slowing down near the bifurcation to spiking, quantitative experimental evidence has been lacking. Using whole-cell patch-clamp recordings from pyramidal neurons and fast-spiking interneurons, we show that 1) the transition to spiking dynamically corresponds to a critical transition exhibiting slowing down, 2) the scaling laws suggest a saddle-node bifurcation governing slowing down, and 3) these precise scaling laws can be used to predict the bifurcation point from a limited window of observation. To our knowledge this is the first report of scaling laws of critical slowing down in an experiment. They present a missing link for a broad class of neuroscience modeling and suggest improved estimation of tipping points by incorporating scaling laws of critical slowing down as a strategy applicable to other complex systems.

Neurons efficiently convey information by being able to switch rapidly between two different states: quiescence and spiking. Such sudden shifts to a qualitatively different state are observed in many complex systems; the often dramatic consequences of these tipping points for diverse fields such as economics, ecology, and the brain have spurred interest to better understand their transition mechanisms and predict their sudden occurrences. By studying the transition from neuronal quiescence to spiking, we show that the quantitative scaling laws for critical slowing down, i.e., a system’s tendency to recover more slowly from perturbations upon approaching its transition point, inform about the underlying bifurcation mechanism and can be used to improve the prediction of a system’s tipping point.

Maladaptive neural synchrony in tinnitus: origin and restoration

Tinnitus is the conscious perception of sound heard in the absence of physical sound sources external or internal to the body, reflected in aberrant neural synchrony of spontaneous or resting-state brain activity. Neural synchrony is generated by the nearly simultaneous firing of individual neurons, of the synchronization of membrane-potential changes in local neural groups as reflected in the local field potentials, resulting in the presence of oscillatory brain waves in the EEG. Noise-induced hearing loss, often resulting in tinnitus, causes a reorganization of the tonotopic map in auditory cortex and increased spontaneous firing rates and neural synchrony. Spontaneous brain rhythms rely on neural synchrony. Abnormal neural synchrony in tinnitus appears to be confined to specific frequency bands of brain rhythms. Increases in delta-band activity are generated by deafferented/deprived neuronal networks resulting from hearing loss. Coordinated reset (CR) stimulation was developed in order to specifically counteract such abnormal neuronal synchrony by desynchronization. The goal of acoustic CR neuromodulation is to desynchronize tinnitus-related abnormal delta-band oscillations. CR neuromodulation does not require permanent stimulus delivery in order to achieve long-lasting desynchronization or even a full-blown anti-kindling but may have cumulative effects, i.e., the effect of different CR epochs separated by pauses may accumulate. Unlike other approaches, acoustic CR neuromodulation does not intend to reduce tinnitus-related neuronal activity by employing lateral inhibition. The potential efficacy of acoustic CR modulation was shown in a clinical proof of concept trial, where effects achieved in 12 weeks of treatment delivered 4–6 h/day persisted through a preplanned 4-week therapy pause and showed sustained long-term effects after 10 months of therapy, leading to 75% responders.

Control of abnormal synchronization in neurological disorders

In the nervous system, synchronization processes play an important role, e.g., in the context of information processing and motor control. However, pathological, excessive synchronization may strongly impair brain function and is a hallmark of several neurological disorders. This focused review addresses the question of how an abnormal neuronal synchronization can specifically be counteracted by invasive and non-invasive brain stimulation as, for instance, by deep brain stimulation for the treatment of Parkinson’s disease, or by acoustic stimulation for the treatment of tinnitus. On the example of coordinated reset (CR) neuromodulation, we illustrate how insights into the dynamics of complex systems contribute to successful model-based approaches, which use methods from synergetics, non-linear dynamics, and statistical physics, for the development of novel therapies for normalization of brain function and synaptic connectivity. Based on the intrinsic multistability of the neuronal populations induced by spike timing-dependent plasticity (STDP), CR neuromodulation utilizes the mutual interdependence between synaptic connectivity and dynamics of the neuronal networks in order to restore more physiological patterns of connectivity via desynchronization of neuronal activity. The very goal is to shift the neuronal population by stimulation from an abnormally coupled and synchronized state to a desynchronized regime with normalized synaptic connectivity, which significantly outlasts the stimulation cessation, so that long-lasting therapeutic effects can be achieved.

Predicting the responses of repetitively firing neurons to current noise

We used phase resetting methods to predict firing patterns of rat subthalamic nucleus (STN) neurons when their rhythmic firing was densely perturbed by noise. We applied sequences of contiguous brief (0.5–2 ms) current pulses with amplitudes drawn from a Gaussian distribution (10–100 pA standard deviation) to autonomously firing STN neurons in slices. Current noise sequences increased the variability of spike times with little or no effect on the average firing rate. We measured the infinitesimal phase resetting curve (PRC) for each neuron using a noise-based method. A phase model consisting of only a firing rate and PRC was very accurate at predicting spike timing, accounting for more than 80% of spike time variance and reliably reproducing the spike-to-spike pattern of irregular firing. An approximation for the evolution of phase was used to predict the effect of firing rate and noise parameters on spike timing variability. It quantitatively predicted changes in variability of interspike intervals with variation in noise amplitude, pulse duration and firing rate over the normal range of STN spontaneous rates. When constant current was used to drive the cells to higher rates, the PRC was altered in size and shape and accurate predictions of the effects of noise relied on incorporating these changes into the prediction. Application of rate-neutral changes in conductance showed that changes in PRC shape arise from conductance changes known to accompany rate increases in STN neurons, rather than the rate increases themselves. Our results show that firing patterns of densely perturbed oscillators cannot readily be distinguished from those of neurons randomly excited to fire from the rest state. The spike timing of repetitively firing neurons may be quantitatively predicted from the input and their PRCs, even when they are so densely perturbed that they no longer fire rhythmically.

Most neurons receive thousands of synaptic inputs per second. Each of these may be individually weak but collectively they shape the temporal pattern of firing by the postsynaptic neuron. If the postsynaptic neuron fires repetitively, its synaptic inputs need not directly trigger action potentials, but may instead control the timing of action potentials that would occur anyway. The phase resetting curve encapsulates the influence of an input on the timing of the next action potential, depending on its time of arrival. We measured the phase resetting curves of neurons in the subthalamic nucleus and used them to accurately predict the timing of action potentials in a phase model subjected to complex input patterns. A simple approximation to the phase model accurately predicted the changes in firing pattern evoked by dense patterns of noise pulses varying in amplitude and pulse duration, and by changes in firing rate. We also showed that the phase resetting curve changes systematically with changes in total neuron conductance, and doing so predicts corresponding changes in firing pattern. Our results indicate that the phase model may accurately represent the temporal integration of complex patterns of input to repetitively firing neurons.

Interplay of intrinsic and synaptic conductances in the generation of high-frequency oscillations in interneuronal networks with irregular spiking

High-frequency oscillations (above 30 Hz) have been observed in sensory and higher-order brain areas, and are believed to constitute a general hallmark of functional neuronal activation. Fast inhibition in interneuronal networks has been suggested as a general mechanism for the generation of high-frequency oscillations. Certain classes of interneurons exhibit subthreshold oscillations, but the effect of this intrinsic neuronal property on the population rhythm is not completely understood. We study the influence of intrinsic damped subthreshold oscillations in the emergence of collective high-frequency oscillations, and elucidate the dynamical mechanisms that underlie this phenomenon. We simulate neuronal networks composed of either Integrate-and-Fire (IF) or Generalized Integrate-and-Fire (GIF) neurons. The IF model displays purely passive subthreshold dynamics, while the GIF model exhibits subthreshold damped oscillations. Individual neurons receive inhibitory synaptic currents mediated by spiking activity in their neighbors as well as noisy synaptic bombardment, and fire irregularly at a lower rate than population frequency. We identify three factors that affect the influence of single-neuron properties on synchronization mediated by inhibition: i) the firing rate response to the noisy background input, ii) the membrane potential distribution, and iii) the shape of Inhibitory Post-Synaptic Potentials (IPSPs). For hyperpolarizing inhibition, the GIF IPSP profile (factor iii)) exhibits post-inhibitory rebound, which induces a coherent spike-mediated depolarization across cells that greatly facilitates synchronous oscillations. This effect dominates the network dynamics, hence GIF networks display stronger oscillations than IF networks. However, the restorative current in the GIF neuron lowers firing rates and narrows the membrane potential distribution (factors i) and ii), respectively), which tend to decrease synchrony. If inhibition is shunting instead of hyperpolarizing, post-inhibitory rebound is not elicited and factors i) and ii) dominate, yielding lower synchrony in GIF networks than in IF networks.

Neurons in the brain engage in collective oscillations at different frequencies. Gamma and high-gamma oscillations (30–100 Hz and higher) have been associated with cognitive functions, and are altered in psychiatric disorders such as schizophrenia and autism. Our understanding of how high-frequency oscillations are orchestrated in the brain is still limited, but it is necessary for the development of effective clinical approaches to the treatment of these disorders. Some neuron types exhibit dynamical properties that can favour synchronization. The theory of weakly coupled oscillators showed how the phase response of individual neurons can predict the patterns of phase relationships that are observed at the network level. However, neurons in vivo do not behave like regular oscillators, but fire irregularly in a regime dominated by fluctuations. Hence, which intrinsic dynamical properties matter for synchronization, and in which regime, is still an open question. Here, we show how single-cell damped subthreshold oscillations enhance synchrony in interneuronal networks by introducing a depolarizing component, mediated by post-inhibitory rebound, that is correlated among neurons due to common inhibitory input.

Global dynamics of a stochastic neuronal oscillator

Nonlinear oscillators have been used to model neurons that fire periodically in the absence of input. These oscillators, which are called neuronal oscillators, share some common response structures with other biological oscillations such as cardiac cells. In this study, we analyze the dependence of the global dynamics of an impulse-driven stochastic neuronal oscillator on the relaxation rate to the limit cycle, the strength of the intrinsic noise, and the impulsive input parameters. To do this, we use a Markov operator that both reflects the density evolution of the oscillator and is an extension of the phase transition curve, which describes the phase shift due to a single isolated impulse. Previously, we derived the Markov operator for the finite relaxation rate that describes the dynamics of the entire phase plane. Here, we construct a Markov operator for the infinite relaxation rate that describes the stochastic dynamics restricted to the limit cycle. In both cases, the response of the stochastic neuronal oscillator to time-varying impulses is described by a product of Markov operators. Furthermore, we calculate the number of spikes between two consecutive impulses to relate the dynamics of the oscillator to the number of spikes per unit time and the interspike interval density. Specifically, we analyze the dynamics of the number of spikes per unit time based on the properties of the Markov operators. Each Markov operator can be decomposed into stationary and transient components based on the properties of the eigenvalues and eigenfunctions. This allows us to evaluate the difference in the number of spikes per unit time between the stationary and transient responses of the oscillator, which we show to be based on the dependence of the oscillator on past activity. Our analysis shows how the duration of the past neuronal activity depends on the relaxation rate, the noise strength, and the impulsive input parameters.

Interaction function of oscillating coupled neurons

Large scale simulations of electrically coupled neuronal oscillators often employ the phase coupled oscillator paradigm to understand and predict network behavior. We study the nature of the interaction between such coupled oscillators using weakly coupled oscillator theory. By employing piecewise linear approximations for phase response curves and voltage time courses and parametrizing their shapes, we compute the interaction function for all such possible shapes and express it in terms of discrete Fourier modes. We find that reasonably good approximation is achieved with four Fourier modes that comprise of both sine and cosine terms.

Optimal number of stimulation contacts for coordinated reset neuromodulation

In this computational study we investigate coordinated reset (CR) neuromodulation designed for an effective control of synchronization by multi-site stimulation of neuronal target populations. This method was suggested to effectively counteract pathological neuronal synchrony characteristic for several neurological disorders. We study how many stimulation sites are required for optimal CR-induced desynchronization. We found that a moderate increase of the number of stimulation sites may significantly prolong the post-stimulation desynchronized transient after the stimulation is completely switched off. This can, in turn, reduce the amount of the administered stimulation current for the intermittent ON-OFF CR stimulation protocol, where time intervals with stimulation ON are recurrently followed by time intervals with stimulation OFF. In addition, we found that the optimal number of stimulation sites essentially depends on how strongly the administered current decays within the neuronal tissue with increasing distance from the stimulation site. In particular, for a broad spatial stimulation profile, i.e., for a weak spatial decay rate of the stimulation current, CR stimulation can optimally be delivered via a small number of stimulation sites. Our findings may contribute to an optimization of therapeutic applications of CR neuromodulation.

Effect of phase response curve skewness on synchronization of electrically coupled neuronal oscillators

We investigate why electrically coupled neuronal oscillators synchronize or fail to synchronize using the theory of weakly coupled oscillators. Stability of synchrony and antisynchrony is predicted analytically and verified using numerical bifurcation diagrams. The shape of the phase response curve (PRC), the shape of the voltage time course, and the frequency of spiking are freely varied to map out regions of parameter spaces that hold stable solutions. We find that type 1 and type 2 PRCs can hold both synchronous and antisynchronous solutions, but the shape of the PRC and the voltage determine the extent of their stability. This is achieved by introducing a five-piecewise linear model to the PRC and a three-piecewise linear model to the voltage time course, and then analyzing the resultant eigenvalue equations that determine the stability of the phase-locked solutions. A single time parameter defines the skewness of the PRC, and another single time parameter defines the spike width and frequency. Our approach gives a comprehensive picture of the relation of the PRC shape, voltage time course, and stability of the resultant synchronous and antisynchronous solutions.

Neural dynamics and information representation in microcircuits of motor cortex

The brain has to analyze and respond to external events that can change rapidly from time to time, suggesting that information processing by the brain may be essentially dynamic rather than static. The dynamical features of neural computation are of significant importance in motor cortex that governs the process of movement generation and learning. In this paper, we discuss these features based primarily on our recent findings on neural dynamics and information coding in the microcircuit of rat motor cortex. In fact, cortical neurons show a variety of dynamical behavior from rhythmic activity in various frequency bands to highly irregular spike firing. Of particular interest are the similarity and dissimilarity of the neuronal response properties in different layers of motor cortex. By conducting electrophysiological recordings in slice preparation, we report the phase response curves (PRCs) of neurons in different cortical layers to demonstrate their layer-dependent synchronization properties. We then study how motor cortex recruits task-related neurons in different layers for voluntary arm movements by simultaneous juxtacellular and multiunit recordings from behaving rats. The results suggest an interesting difference in the spectrum of functional activity between the superficial and deep layers. Furthermore, the task-related activities recorded from various layers exhibited power law distributions of inter-spike intervals (ISIs), in contrast to a general belief that ISIs obey Poisson or Gamma distributions in cortical neurons. We present a theoretical argument that this power law of in vivo neurons may represent the maximization of the entropy of firing rate with limited energy consumption of spike generation. Though further studies are required to fully clarify the functional implications of this coding principle, it may shed new light on information representations by neurons and circuits in motor cortex.

Origin of intrinsic irregular firing in cortical interneurons

Cortical spike trains are highly irregular both during ongoing, spontaneous activity and when driven at high firing rates. There is uncertainty about the source of this irregularity, ranging from intrinsic noise sources in neurons to collective effects in large-scale cortical networks. Cortical interneurons display highly irregular spike times (coefficient of variation of the interspike intervals >1) in response to dc-current injection in vitro. This is in marked contrast to cortical pyramidal cells, which spike highly irregularly in vivo, but regularly in vitro. We show with in vitro recordings and computational models that this is due to the fast activation kinetics of interneuronal K(+) currents. This explanation holds over a wide parameter range and with Gaussian white, power-law, and Ornstein-Uhlenbeck noise. The intrinsically irregular spiking of interneurons could contribute to the irregularity of the cortical network.

Toggling between gamma-frequency activity and suppression of cell assemblies

Gamma (30–80 Hz) rhythms in hippocampus and neocortex resulting from the interaction of excitatory and inhibitory cells (E- and I-cells), called Pyramidal-Interneuronal Network Gamma (PING), require that the I-cells respond to the E-cells, but don't fire on their own. In idealized models, there is a sharp boundary between a parameter regime where the I-cells have weak-enough drive for PING, and one where they have so much drive that they fire without being prompted by the E-cells. In the latter regime, they often de-synchronize and suppress the E-cells; the boundary was therefore called the “suppression boundary” by Börgers and Kopell (2005). The model I-cells used in the earlier work by Börgers and Kopell have a “type 1” phase response, i.e., excitatory input always advances them. However, fast-spiking inhibitory basket cells often have a “type 2” phase response: Excitatory input arriving soon after they fire delays them. We study the effect of the phase response type on the suppression transition, under the additional assumption that the I-cells are kept synchronous by gap junctions. When many E-cells participate on a given cycle, the resulting excitation advances the I-cells on the next cycle if their phase response is of type 1, and this can result in suppression of more E-cells on the next cycle. Therefore, strong E-cell spike volleys tend to be followed by weaker ones, and vice versa. This often results in erratic fluctuations in the strengths of the E-cell spike volleys. When the phase response of the I-cells is of type 2, the opposite happens: strong E-cell spike volleys delay the inhibition on the next cycle, therefore tend to be followed by yet stronger ones. The strengths of the E-cell spike volleys don't oscillate, and there is a nearly abrupt transition from PING to ING (a rhythm involving I-cells only).

Effect of sharp jumps at the edges of phase response curves on synchronization of electrically coupled neuronal oscillators

We study synchronization phenomenon of coupled neuronal oscillators using the theory of weakly coupled oscillators. The role of sudden jumps in the phase response curve profiles found in some experimental recordings and models on the ability of coupled neurons to exhibit synchronous and antisynchronous behavior is investigated, when the coupling between the neurons is electrical. The level of jumps in the phase response curve at either end, spike width and frequency of voltage time course of the coupled neurons are parameterized using piecewise linear functional forms, and the conditions for stable synchrony and stable antisynchrony in terms of those parameters are computed analytically. The role of the peak position of the phase response curve on phase-locking is also investigated.

Hippocampal CA1 pyramidal neurons exhibit type 1 phase-response curves and type 1 excitability

Phase-resetting properties of neurons determine their functionality as integrators (type 1) vs. resonators (type 2), as well as their synchronization tendencies. We introduce a novel, bias-correction method to estimate the infinitesimal phase-resetting curve (iPRC) and confirm type 1 excitability in hippocampal pyramidal CA1 neurons in vitro by two independent methods. First, PRCs evoked using depolarizing pulses consisted only of advances, consistent with type 1. Second, the frequency/current (f/I) plots showed no minimum frequency, again consistent with type 1. Type 1 excitability was also confirmed by the absence of a resonant peak in the interspike interval histograms derived from the f/I data. The PRC bias correction assumed that the distribution of noisy phase resetting is truncated, because an input cannot advance a spike to a point in time before the input (the causal limit) and successfully removed the statistical bias for delays in the null PRC in response to zero-magnitude input by computing the phase resetting as the mean of the untruncated distribution. The PRC for depolarization peaked at late phases and decreased to zero by the end of the cycle, whereas delays observed in response to hyperpolarization increased monotonically. The bias correction did not affect this difference in shape, which was due instead to the causal limit obscuring the iPRC for depolarization but not hyperpolarization. Our results suggest that weak periodic hyperpolarizing drive can theoretically entrain CA1 pyramidal neurons at any phase but that strong excitation will preferentially phase-lock them with zero time lag.

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.

Probabilistic convergence guarantees for type-II pulse-coupled oscillators

We show that a large class of pulse-coupled oscillators converge with high probability from random initial conditions on a large class of graphs with time delays. Our analysis combines previous local convergence results, probabilistic network analysis, and a classification scheme for type-II phase response curves to produce rigorous lower bounds for convergence probabilities based on network density. These results suggest methods for the analysis of pulse-coupled oscillators, and provide insights into the balance of excitation and inhibition in the operation of biological type-II phase response curves and also the design of decentralized and minimal clock synchronization schemes in sensor nets.

Phase response curves of subthalamic neurons measured with synaptic input and current injection

Infinitesimal phase response curves (iPRCs) provide a simple description of the response of repetitively firing neurons and may be used to predict responses to any pattern of synaptic input. Their simplicity makes them useful for understanding the dynamics of neurons when certain conditions are met. For example, the sizes of evoked phase shifts should scale linearly with stimulus strength, and the form of the iPRC should remain relatively constant as firing rate varies. We measured the PRCs of rat subthalamic neurons in brain slices using corticosubthalamic excitatory postsynaptic potentials (EPSPs; mediated by both AMPA- and NMDA-type receptors) and injected current pulses and used them to calculate the iPRC. These were relatively insensitive to both the size of the stimulus and the cell's firing rate, suggesting that the iPRC can predict the response of subthalamic nucleus cells to extrinsic inputs. However, the iPRC calculated using EPSPs differed from that obtained using current pulses. EPSPs (normalized for charge) were much more effective at altering the phase of subthalamic neurons than current pulses. The difference was not attributable to the extended time course of NMDA receptor-mediated currents, being unaffected by blockade of NMDA receptors. The iPRC provides a good description of subthalamic neurons' response to input, but iPRCs are best estimated using synaptic inputs rather than somatic current injection.

Efficient estimation of phase-response curves via compressive sensing

The phase-response curve (PRC), relating the phase shift of an oscillator to external perturbation, is an important tool to study neurons and their population behavior. It can be experimentally estimated by measuring the phase changes caused by probe stimuli. These stimuli, usually short pulses or continuous noise, have a much wider frequency spectrum than that of neuronal dynamics. This makes the experimental data high dimensional while the number of data samples tends to be small. Current PRC estimation methods have not been optimized for efficiently discovering the relevant degrees of freedom from such data. We propose a systematic and efficient approach based on a recently developed signal processing theory called compressive sensing (CS). CS is a framework for recovering sparsely constructed signals from undersampled data and is suitable for extracting information about the PRC from finite but high-dimensional experimental measurements. We illustrate how the CS algorithm can be translated into an estimation scheme and demonstrate that our CS method can produce good estimates of the PRCs with simulated and experimental data, especially when the data size is so small that simple approaches such as naive averaging fail. The tradeoffs between degrees of freedom vs. goodness-of-fit were systematically analyzed, which help us to understand better what part of the data has the most predictive power. Our results illustrate that finite sizes of neuroscientific data in general compounded by large dimensionality can hamper studies of the neural code and suggest that CS is a good tool for overcoming this challenge.

Identifying type I excitability using dynamics of stochastic neural firing patterns

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

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.

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

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

Desynchronizing electrical and sensory coordinated reset neuromodulation

Coordinated reset (CR) stimulation is a desynchronizing stimulation technique based on timely coordinated phase resets of sub-populations of a synchronized neuronal ensemble. It has initially been computationally developed for electrical deep brain stimulation (DBS), to enable an effective desynchronization and unlearning of pathological synchrony and connectivity (anti-kindling). Here we computationally show for ensembles of spiking and bursting model neurons interacting via excitatory and inhibitory adaptive synapses that a phase reset of neuronal populations as well as a desynchronization and an anti-kindling can robustly be achieved by direct electrical stimulation or indirect (synaptically-mediated) excitatory and inhibitory stimulation. Our findings are relevant for DBS as well as for sensory stimulation in neurological disorders characterized by pathological neuronal synchrony. Based on the obtained results, we may expect that the local effects in the vicinity of a depth electrode (realized by direct stimulation of the neurons' somata or stimulation of axon terminals) and the non-local CR effects (realized by stimulation of excitatory or inhibitory efferent fibers) of deep brain CR neuromodulation may be similar or even identical. Furthermore, our results indicate that an effective desynchronization and anti-kindling can even be achieved by non-invasive, sensory CR neuromodulation. We discuss the concept of sensory CR neuromodulation in the context of neurological disorders.

Engineering the synchronization of neuron action potentials using global time-delayed feedback stimulation

We experimentally demonstrate the use of continuous, time-delayed, feedback stimulation for controlling the synchronization of neuron action potentials. Phase-based models were experimentally constructed from a single synaptically isolated cultured hippocampal neuron. These models were used to determine the stimulation parameters necessary to produce the desired synchronization behavior in the action potentials of a pair of neurons coupled through a global time-delayed interaction. Measurements made using a dynamic clamp system confirm the generation of the synchronized states predicted by the experimentally constructed phase model. This model was then utilized to extrapolate the feedback stimulation parameters necessary to disrupt the action potential synchronization of a large population of globally interacting neurons.

Measurement of infinitesimal phase response curves from noisy real neurons

We sought to measure infinitesimal phase response curves (iPRCs) from rat hippocampal CA1 pyramidal neurons. It is difficult to measure iPRCs from noisy neurons because of the dilemma that either the linearity or the signal-to-noise ratio of responses to external perturbations must be sacrificed. To overcome this difficulty, we used an iPRC measurement model formulated as the Langevin phase equation (LPE) to extract iPRCs in the Bayesian scheme. We then simultaneously verified the effectiveness of the measurement model and the reliability of the estimated iPRCs by demonstrating that LPEs with the estimated iPRCs could predict the stochastic behaviors of the same neurons, whose iPRCs had been measured, when they were perturbed by periodic stimulus currents. Our results suggest that the LPE is an effective model for real oscillating neurons and that many theoretical frameworks based on it may be applicable to real nerve systems.

Sensory feedback, error correction, and remapping in a multiple oscillator model of place-cell activity

Mammals navigate by integrating self-motion signals ("path integration") and occasionally fixing on familiar environmental landmarks. The rat hippocampus is a model system of spatial representation in which place cells are thought to integrate both sensory and spatial information from entorhinal cortex. The localized firing fields of hippocampal place cells and entorhinal grid-cells demonstrate a phase relationship with the local theta (6-10 Hz) rhythm that may be a temporal signature of path integration. However, encoding self-motion in the phase of theta oscillations requires high temporal precision and is susceptible to idiothetic noise, neuronal variability, and a changing environment. We present a model based on oscillatory interference theory, previously studied in the context of grid cells, in which transient temporal synchronization among a pool of path-integrating theta oscillators produces hippocampal-like place fields. We hypothesize that a spatiotemporally extended sensory interaction with external cues modulates feedback to the theta oscillators. We implement a form of this cue-driven feedback and show that it can retrieve fixed points in the phase code of position. A single cue can smoothly reset oscillator phases to correct for both systematic errors and continuous noise in path integration. Further, simulations in which local and global cues are rotated against each other reveal a phase-code mechanism in which conflicting cue arrangements can reproduce experimentally observed distributions of "partial remapping" responses. This abstract model demonstrates that phase-code feedback can provide stability to the temporal coding of position during navigation and may contribute to the context-dependence of hippocampal spatial representations. While the anatomical substrates of these processes have not been fully characterized, our findings suggest several signatures that can be evaluated in future experiments.

Optimal phase response curves for stochastic synchronization of limit-cycle oscillators by common Poisson noise

We consider optimization of phase response curves for stochastic synchronization of noninteracting limit-cycle oscillators by common Poisson impulsive signals. The optimal functional shape for sufficiently weak signals is sinusoidal, but can differ for stronger signals. By solving the Euler-Lagrange equation associated with the minimization of the Lyapunov exponent characterizing synchronization efficiency, the optimal phase response curve is obtained. We show that the optimal shape mutates from a sinusoid to a sawtooth as the constraint on its squared amplitude is varied.

Cellularly-driven differences in network synchronization propensity are differentially modulated by firing frequency

Spatiotemporal pattern formation in neuronal networks depends on the interplay between cellular and network synchronization properties. The neuronal phase response curve (PRC) is an experimentally obtainable measure that characterizes the cellular response to small perturbations, and can serve as an indicator of cellular propensity for synchronization. Two broad classes of PRCs have been identified for neurons: Type I, in which small excitatory perturbations induce only advances in firing, and Type II, in which small excitatory perturbations can induce both advances and delays in firing. Interestingly, neuronal PRCs are usually attenuated with increased spiking frequency, and Type II PRCs typically exhibit a greater attenuation of the phase delay region than of the phase advance region. We found that this phenomenon arises from an interplay between the time constants of active ionic currents and the interspike interval. As a result, excitatory networks consisting of neurons with Type I PRCs responded very differently to frequency modulation compared to excitatory networks composed of neurons with Type II PRCs. Specifically, increased frequency induced a sharp decrease in synchrony of networks of Type II neurons, while frequency increases only minimally affected synchrony in networks of Type I neurons. These results are demonstrated in networks in which both types of neurons were modeled generically with the Morris-Lecar model, as well as in networks consisting of Hodgkin-Huxley-based model cortical pyramidal cells in which simulated effects of acetylcholine changed PRC type. These results are robust to different network structures, synaptic strengths and modes of driving neuronal activity, and they indicate that Type I and Type II excitatory networks may display two distinct modes of processing information.

Synchronization of the firing of neurons in the brain is related to many cognitive functions, such as recognizing faces, discriminating odors, and coordinating movement. It is therefore important to understand what properties of neuronal networks promote synchrony of neural firing. One measure that is often used to determine the contribution of individual neurons to network synchrony is called the phase response curve (PRC). PRCs describe how the timing of neuronal firing changes depending on when input, such as a synaptic signal, is received by the neuron. A characteristic of PRCs that has previously not been well understood is that they change dramatically as the neuron's firing frequency is modulated. This effect carries potential significance, since cognitive functions are often associated with specific frequencies of network activity in the brain. We showed computationally that the frequency dependence of PRCs can be explained by the relative timing of ionic membrane currents with respect to the time between spike firings. Our simulations also showed that the frequency dependence of neuronal PRCs leads to frequency-dependent changes in network synchronization that can be different for different neuron types. These results further our understanding of how synchronization is generated in the brain to support various cognitive functions.

Robust convergence in pulse-coupled oscillators with delays

We show that for pulse-coupled oscillators a class of phase response curves with both excitation and inhibition exhibit robust convergence to synchrony on arbitrary aperiodic connected graphs with delays. We describe the basins of convergence and give explicit bounds on the convergence times. These results provide new and more robust methods for synchronization of sensor nets and also have biological implications.