The influence of limit cycle topology on the phase resetting curve

Dopamine receptor antagonists effects on low-dimensional attractors of local field potentials in optogenetic mice

The goal of this study was to investigate the effects of acute cocaine injection or dopamine (DA) receptor antagonists on the medial prefrontal cortex (mPFC) gamma oscillations and their relationship to short term neuroadaptation that may mediate addiction. For this purpose, optogenetically evoked local field potentials (LFPs) in response to a brief 10 ms laser light pulse were recorded from 17 mice. D1-like receptor antagonist SCH 23390 or D2-like receptor antagonist sulpiride, or both, were administered either before or after cocaine. A Euclidian distance-based dendrogram classifier separated the 100 trials for each animal in disjoint clusters. When baseline and DA receptor antagonists trials were combined in a single trial, a minimum of 20% overlap occurred in some dendrogram clusters, which suggests a possible common, invariant, dynamic mechanism shared by both baseline and DA receptor antagonists data. The delay-embedding method of neural activity reconstruction was performed using the correlation time and mutual information to determine the lag/correlation time of LFPs and false nearest neighbors to determine the embedding dimension. We found that DA receptor antagonists applied before cocaine cancels out the effect of cocaine and leaves the lag time distributions at baseline values. On the other hand, cocaine applied after DA receptor antagonists shifts the lag time distributions to longer durations, i.e. increase the correlation time of LFPs. Fourier analysis showed that a reasonable accurate decomposition of the LFP data can be obtained with a relatively small (less than ten) Fourier coefficients.

Cocaine-Induced Changes in Low-Dimensional Attractors of Local Field Potentials in Optogenetic Mice

Optogenetically evoked local field potential (LFP) recorded from the medial prefrontal cortex (mPFC) of mice during basal conditions and following a systemic cocaine administration were analyzed. Blue light stimuli were delivered to mPFC through a fiber optic every 2 s and each trial was repeated 100 times. As in the previous study, we used a surrogate data method to check that nonlinearity was present in the experimental LFPs and only used the last 1.5 s of steady activity to measure the LFPs phase resetting induced by the brief 10 ms light stimulus. We found that the steady dynamics of the mPFC in response to light stimuli could be reconstructed in a three-dimensional phase space with topologically similar "8"-shaped attractors across different animals. Therefore, cocaine did not change the complexity of the recorded nonlinear data compared to the control case. The phase space of the reconstructed attractor is determined by the LFP time series and its temporally shifted versions by a multiple of some lag time. We also compared the change in the attractor shape between cocaine-injected and control using (1) dendrogram clustering and (2) Frechet distance. We found about 20% overlap between control and cocaine trials when classified using dendrogram method, which suggest that it may be possible to describe mathematically both data sets with the same model and slightly different model parameters. We also found that the lag times are about three times shorter for cocaine trials compared to control. As a result, although the phase space trajectories for control and cocaine may look similar, their dynamics is significantly different.

Predicting the Existence and Stability of Phase-Locked Mode in Neural Networks Using Generalized Phase-Resetting Curve

We used the phase-resetting method to study a biologically relevant three-neuron network in which one neuron receives multiple inputs per cycle. For this purpose, we first generalized the concept of phase resetting to accommodate multiple inputs per cycle. We explicitly showed how analytical conditions for the existence and the stability of phase-locked modes are derived. In particular, we solved newly derived recursive maps using as an example a biologically relevant driving-driven neural network with a dynamic feedback loop. We applied the generalized phase-resetting definition to predict the relative-phase and the stability of a phase-locked mode in open loop setup. We also compared the predicted phase-locked mode against numerical simulations of the fully connected network.

Globally attracting synchrony in a network of oscillators with all-to-all inhibitory pulse coupling

The synchronization tendencies of networks of oscillators have been studied intensely. We assume a network of all-to-all pulse-coupled oscillators in which the effect of a pulse is independent of the number of oscillators that simultaneously emit a pulse and the normalized delay (the phase resetting) is a monotonically increasing function of oscillator phase with the slope everywhere less than 1 and a value greater than 2φ-1, where φ is the normalized phase. Order switching cannot occur; the only possible solutions are globally attracting synchrony and cluster solutions with a fixed firing order. For small conduction delays, we prove the former stable and all other possible attractors nonexistent due to the destabilizing discontinuity of the phase resetting at a phase of 0.

Low-dimensional attractor for neural activity from local field potentials in optogenetic mice

We used optogenetic mice to investigate possible nonlinear responses of the medial prefrontal cortex (mPFC) local network to light stimuli delivered by a 473 nm laser through a fiber optics. Every 2 s, a brief 10 ms light pulse was applied and the local field potentials (LFPs) were recorded with a 10 kHz sampling rate. The experiment was repeated 100 times and we only retained and analyzed data from six animals that showed stable and repeatable response to optical stimulations. The presence of nonlinearity in our data was checked using the null hypothesis that the data were linearly correlated in the temporal domain, but were random otherwise. For each trail, 100 surrogate data sets were generated and both time reversal asymmetry and false nearest neighbor (FNN) were used as discriminating statistics for the null hypothesis. We found that nonlinearity is present in all LFP data. The first 0.5 s of each 2 s LFP recording were dominated by the transient response of the networks. For each trial, we used the last 1.5 s of steady activity to measure the phase resetting induced by the brief 10 ms light stimulus. After correcting the LFPs for the effect of phase resetting, additional preprocessing was carried out using dendrograms to identify “similar” groups among LFP trials. We found that the steady dynamics of mPFC in response to light stimuli could be reconstructed in a three-dimensional phase space with topologically similar “8”-shaped attractors across different animals. Our results also open the possibility of designing a low-dimensional model for optical stimulation of the mPFC local network.

Local linear approximation of the Jacobian matrix better captures phase resetting of neural limit cycle oscillators

One effect of any external perturbations, such as presynaptic inputs, received by limit cycle oscillators when they are part of larger neural networks is a transient change in their firing rate, or phase resetting. A brief external perturbation moves the figurative point outside the limit cycle, a geometric perturbation that we mapped into a transient change in the firing rate, or a temporal phase resetting. In order to gain a better qualitative understanding of the link between the geometry of the limit cycle and the phase resetting curve (PRC), we used a moving reference frame with one axis tangent and the others normal to the limit cycle. We found that the stability coefficients associated with the unperturbed limit cycle provided good quantitative predictions of both the tangent and the normal geometric displacements induced by external perturbations. A geometric-to-temporal mapping allowed us to correctly predict the PRC while preserving the intuitive nature of this geometric approach.

Biophysical basis of the phase response curve of subthalamic neurons with generalization to other cell types

Experimental evidence indicates that the response of subthalamic neurons to excitatory postsynaptic potentials (EPSPs) is well described by their infinitesimal phase response curves (iPRC). However, the factors controlling the shape of that iPRC, and hence controlling the way subthalamic neurons respond to synaptic input, are unclear. We developed a biophysical model of subthalamic neurons to aid in the understanding of their iPRCs; this model exhibited an iPRC type common to many subthalamic cells. We devised a method for deriving its iPRC from its biophysical properties that clarifies how these different properties interact to shape the iPRC. This method revealed why the response of subthalamic neurons is well approximated by their iPRCs and how that approximation becomes less accurate under strong fluctuating input currents. It also connected iPRC structure to aspects of cellular physiology that could be estimated in simple current-clamp experiments. This allowed us to directly compare the iPRC predicted by our theory with the iPRC estimated from the response to EPSPs or current pulses in individual cells. We found that theoretically predicted iPRCs agreed well with estimates derived from synaptic stimuli, but not with those estimated from the response to somatic current injection. The difference between synaptic currents and those applied experimentally at the soma may arise from differences in the dynamics of charge redistribution on the dendrites and axon. Ultimately, our approach allowed us to identify novel ways in which voltage-dependent conductances interact with AHP conductances to influence synaptic integration that will apply to a wide range of cell types.

Inhibitory feedback promotes stability in an oscillatory network

Reliability and variability of neuronal activity are both thought to be important for the proper function of neuronal networks. The crustacean pyloric rhythm (∼1 Hz) is driven by a group of pacemaker neurons (AB/PD) that inhibit and burst out of phase with all follower pyloric neurons. The only known chemical synaptic feedback to the pacemakers is an inhibitory synapse from the follower lateral pyloric (LP) neuron. Although this synapse has been studied extensively, its role in the generation and coordination of the pyloric rhythm is unknown. We examine the hypothesis that this synapse acts to stabilize the oscillation by reducing the variability in cycle period on a cycle-by-cycle basis. Our experimental data show that functionally removing the LP-pyloric dilator (PD) synapse by hyperpolarizing the LP neuron increases the pyloric period variability. The increase in pyloric rhythm stability in the presence of the LP-PD synapse is demonstrated by a decrease in the amplitude of the phase response curve of the PD neuron. These experimental results are explained by a reduced mathematical model. Phase plane analysis of this model demonstrates that the effect of the periodic inhibition is to produce asymptotic stability in the oscillation phase, which leads to a reduction in variability of the oscillation cycle period.

Stability of two cluster solutions in pulse coupled networks of neural oscillators

Phase response curves (PRCs) have been widely used to study synchronization in neural circuits comprised of pacemaking neurons. They describe how the timing of the next spike in a given spontaneously firing neuron is affected by the phase at which an input from another neuron is received. Here we study two reciprocally coupled clusters of pulse coupled oscillatory neurons. The neurons within each cluster are presumed to be identical and identically pulse coupled, but not necessarily identical to those in the other cluster. We investigate a two cluster solution in which all oscillators are synchronized within each cluster, but in which the two clusters are phase locked at nonzero phase with each other. Intuitively, one might expect this solution to be stable only when synchrony within each isolated cluster is stable, but this is not the case. We prove rigorously the stability of the two cluster solution and show how reciprocal coupling can stabilize synchrony within clusters that cannot synchronize in isolation. These stability results for the two cluster solution suggest a mechanism by which reciprocal coupling between brain regions can induce local synchronization via the network feedback loop.

Synaptic and intrinsic determinants of the phase resetting curve for weak coupling

A phase resetting curve (PRC) keeps track of the extent to which a perturbation at a given phase advances or delays the next spike, and can be used to predict phase locking in networks of oscillators. The PRC can be estimated by convolving the waveform of the perturbation with the infinitesimal PRC (iPRC) under the assumption of weak coupling. The iPRC is often defined with respect to an infinitesimal current as z(i)(ϕ), where ϕ is phase, but can also be defined with respect to an infinitesimal conductance change as z(g)(ϕ). In this paper, we first show that the two approaches are equivalent. Coupling waveforms corresponding to synapses with different time courses sample z(g)(ϕ) in predictably different ways. We show that for oscillators with Type I excitability, an anomalous region in z(g)(ϕ) with opposite sign to that seen otherwise is often observed during an action potential. If the duration of the synaptic perturbation is such that it effectively samples this region, PRCs with both advances and delays can be observed despite Type I excitability. We also show that changing the duration of a perturbation so that it preferentially samples regions of stable or unstable slopes in z(g)(ϕ) can stabilize or destabilize synchrony in a network with the corresponding dynamics.

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

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

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

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

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

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

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

Rate constants rather than biochemical mechanism determine behaviour of genetic clocks

Many biological systems contain both positive and negative feedbacks. These are often classified as resonators or integrators. Resonators respond preferentially to oscillating signals of a particular frequency. Integrators, on the other hand, accumulate a response to signals. Computational neuroscientists often refer to neurons showing integrator properties as type I neurons and those showing resonator properties as type II neurons. Guantes & Poyatos have shown that type I or type II behaviour can be seen in genetic clocks. They argue that when negative feedback occurs through transcription regulation and post-translationally, genetic clocks act as integrators and resonators, respectively. Here we show that either behaviour can be seen with either design and in a wide range of genetic clocks. This highlights the importance of parameters rather than biochemical mechanism in determining the system behaviour.

Multifunctional pattern-generating circuits

The ability of distinct anatomical circuits to generate multiple behavioral patterns is widespread among vertebrate and invertebrate species. These multifunctional neuronal circuits are the result of multistable neural dynamics and modular organization. The evidence suggests multifunctional circuits can be classified by distinct architectures, yet the activity patterns of individual neurons involved in more than one behavior can vary dramatically. Several mechanisms, including sensory input, the parallel activity of projection neurons, neuromodulation, and biomechanics, are responsible for the switching between patterns. Recent advances in both analytical and experimental tools have aided the study of these complex circuits.

Loss of phase-locking in non-weakly coupled inhibitory networks of type-I model neurons

Synchronization of excitable cells coupled by reciprocal inhibition is a topic of significant interest due to the important role that inhibitory synaptic interaction plays in the generation and regulation of coherent rhythmic activity in a variety of neural systems. While recent work revealed the synchronizing influence of inhibitory coupling on the dynamics of many networks, it is known that strong coupling can destabilize phase-locked firing. Here we examine the loss of synchrony caused by an increase in inhibitory coupling in networks of type-I Morris-Lecar model oscillators, which is characterized by a period-doubling cascade and leads to mode-locked states with alternation in the firing order of the two cells, as reported recently by Maran and Canavier (J Comput Nerosci, 2008) for a network of Wang-Buzsáki model neurons. Although alternating-order firing has been previously reported as a near-synchronous state, we show that the stable phase difference between the spikes of the two Morris-Lecar cells can constitute as much as 70% of the unperturbed oscillation period. Further, we examine the generality of this phenomenon for a class of type-I oscillators that are close to their excitation thresholds, and provide an intuitive geometric description of such "leap-frog" dynamics. In the Morris-Lecar model network, the alternation in the firing order arises under the condition of fast closing of K( + ) channels at hyperpolarized potentials, which leads to slow dynamics of membrane potential upon synaptic inhibition, allowing the presynaptic cell to advance past the postsynaptic cell in each cycle of the oscillation. Further, we show that non-zero synaptic decay time is crucial for the existence of leap-frog firing in networks of phase oscillators. However, we demonstrate that leap-frog spiking can also be obtained in pulse-coupled inhibitory networks of one-dimensional oscillators with a multi-branched phase domain, for instance in a network of quadratic integrate-and-fire model cells. Finally, for the case of a homogeneous network, we establish quantitative conditions on the phase resetting properties of each cell necessary for stable alternating-order spiking, complementing the analysis of Goel and Ermentrout (Physica D 163:191-216, 2002) of the order-preserving phase transition map.

Phase response curves in the characterization of epileptiform activity

Coordinated cellular activity is a major characteristic of nervous system function. Coupled oscillator theory offers unique avenues to address cellular coordination phenomena. In this study, we focus on the characterization of the dynamics of epileptiform activity, based on some seizures that manifest themselves with very periodic rhythmic activity, termed absence seizures. Our approach consists in obtaining experimentally the phase response curves (PRCs) in the neocortex and thalamus, and incorporating these PRCs into a model of coupled oscillators. Phase preferences of the stationary states and their stability are determined, and these results from the model are compared with the experimental recordings, and interpreted in physiological terms.

Using phase resetting to predict 1:1 and 2:2 locking in two neuron networks in which firing order is not always preserved

Our goal is to understand how nearly synchronous modes arise in heterogenous networks of neurons. In heterogenous networks, instead of exact synchrony, nearly synchronous modes arise, which include both 1:1 and 2:2 phase-locked modes. Existence and stability criteria for 2:2 phase-locked modes in reciprocally coupled two neuron circuits were derived based on the open loop phase resetting curve (PRC) without the assumption of weak coupling. The PRC for each component neuron was generated using the change in synaptic conductance produced by a presynaptic action potential as the perturbation. Separate derivations were required for modes in which the firing order is preserved and for those in which it alternates. Networks composed of two model neurons coupled by reciprocal inhibition were examined to test the predictions. The parameter regimes in which both types of nearly synchronous modes are exhibited were accurately predicted both qualitatively and quantitatively provided that the synaptic time constant is short with respect to the period and that the effect of second order resetting is considered. In contrast, PRC methods based on weak coupling could not predict 2:2 modes and did not predict the 1:1 modes with the level of accuracy achieved by the strong coupling methods. The strong coupling prediction methods provide insight into what manipulations promote near-synchrony in a two neuron network and may also have predictive value for larger networks, which can also manifest changes in firing order. We also identify a novel route by which synchrony is lost in mildly heterogenous networks.

Beyond Two-Cell Networks: Experimental Measurement of Neuronal Responses to Multiple Synaptic Inputs

Oscillations of large populations of neurons are thought to be important in the normal functioning of the brain. We have used phase response curve (PRC) methods to characterize the dynamics of single neurons and predict population dynamics. Our past experimental work was limited to special circumstances (e.g., 2-cell networks of periodically firing neurons). Here, we explore the feasibility of extending our methods to predict the synchronization properties of stellate cells (SCs) in the rat entorhinal cortex under broader conditions. In particular, we test the hypothesis that PRCs in SCs scale linearly with changes in synaptic amplitude, and measure how well responses to Poisson process-driven inputs can be predicted in terms of PRCs. Although we see nonlinear responses to excitatory and inhibitory inputs, we find that models based on weak coupling account for scaling and Poisson process-driven inputs reasonably accurately.

The combined effects of inhibitory and electrical synapses in synchrony

Recent experimental results have shown that GABAergic interneurons in the central nervous system are frequently connected via electrical synapses. Hence, depending on the area or the subpopulation, interneurons interact via inhibitory synapses or electrical synapses alone or via both types of interactions. The theoretical work presented here addresses the significance of these different modes of interactions for the interneuron networks dynamics. We consider the simplest system in which this issue can be investigated in models or in experiments: a pair of neurons, interacting via electrical synapses, inhibitory synapses, or both, and activated by the injection of a noisy external current. Assuming that the couplings and the noise are weak, we derive an analytical expression relating the cross-correlation (CC) of the activity of the two neurons to the phase response function of the neurons. When electrical and inhibitory interactions are not too strong, they combine their effect in a linear manner. In this regime, the effect of electrical and inhibitory interactions when combined can be deduced knowing the effects of each of the interactions separately. As a consequence, depending on intrinsic neuronal properties, electrical and inhibitory synapses may cooperate, both promoting synchrony, or may compete, with one promoting synchrony while the other impedes it. In contrast, for sufficiently strong couplings, the two types of synapses combine in a nonlinear fashion. Remarkably, we find that in this regime, combining electrical synapses with inhibition amplifies synchrony, whereas electrical synapses alone would desynchronize the activity of the neurons. We apply our theory to predict how the shape of the CC of two neurons changes as a function of ionic channel conductances, focusing on the effect of persistent sodium conductance, of the firing rate of the neurons and the nature and the strength of their interactions. These predictions may be tested using dynamic clamp techniques.

Efficient estimation of phase-resetting curves in real neurons and its significance for neural-network modeling

The phase-resetting curve (PRC) of a neural oscillator describes the effect of a perturbation on its periodic motion and is therefore useful to study how the neuron responds to stimuli and whether it phase locks to other neurons in a network. Combining theory, computer simulations and electrophysiological experiments we present a simple method for estimating the PRC of real neurons. This allows us to simplify the complex dynamics of a single neuron to a phase model. We also illustrate how to infer the existence of coherent network activity from the estimated PRC.

Scaling of prediction error does not confirm chaotic dynamics underlying irregular firing using interspike intervals from midbrain dopamine neurons

Dopamine neurons in the substantia nigra pars compacta often fire in an irregular, single spike mode in vivo, and a similar firing pattern can be observed in vitro when small conductance calcium-activated potassium channel blockers are applied to the bath. It is not clear whether the irregular firing is due to stochastic processes or nonlinear deterministic processes. A previous study [Neuroscience 104 (2001) 829] used nonlinear forecasting methods applied to a continuous function derived from the interspike interval (ISI) data from irregularly firing dopamine neurons to show that the predictability scaled exponentially with forecast horizon and was consistent with nonlinear deterministic chaos. However, we show here that the observed exponential scaling is also consistent with a stochastic process, because it did not differ significantly from that of shuffled surrogate data. On the other hand, nonlinear forecasting directly from the ISI data using the package TISEAN provided some evidence for nonlinear deterministic structure in four of five records obtained in vitro and in two of nine records obtained in vivo. Although we cannot rule out a role for nonlinear chaotic dynamics in structuring the firing pattern, we suggest an alternate hypothesis that includes a mechanism by which the firing pattern can become more variable in the presence of a constant level of background noise.

Annihilation of single cell neural oscillations by feedforward and feedback control

Annihilation of neural oscillation by localized electrical stimulation has been shown to be a promising treatment modality in a number of neural diseases like Parkinson disease or epilepsy. The contributions presented in this manuscript comprise newly developed stimulation schemes to achieve annihilation of action potential generation and action potential propagation. The ability to achieve oscillation annihilation is demonstrated in computer simulations with a single compartment nerve cell model (annihilation of action potential generation), and with a multi-compartment nerve fiber model (annihilation of action potential propagation). Additionally, we show superiority of the new feedback based schemes over the feedforward schemes in terms of effectiveness, phase robustness, and reduced sensitivity to disturbances. At the end we propose a conditioned switched feedback control regime to be applied in actual applications, where oscillation annihilation is needed.

Synchronization in hybrid neuronal networks of the hippocampal formation

Understanding the mechanistic bases of neuronal synchronization is a current challenge in quantitative neuroscience. We studied this problem in two putative cellular pacemakers of the mammalian hippocampal theta rhythm: glutamatergic stellate cells (SCs) of the medial entorhinal cortex and GABAergic oriens-lacunosum-molecular (O-LM) interneurons of hippocampal region CA1. We used two experimental methods. First, we measured changes in spike timing induced by artificial synaptic inputs applied to individual neurons. We then measured responses of free-running hybrid neuronal networks, consisting of biological neurons coupled (via dynamic clamp) to biological or virtual counterparts. Results from the single-cell experiments predicted network behaviors well and are compatible with previous model-based predictions of how specific membrane mechanisms give rise to empirically measured synchronization behavior. Both cell types phase lock stably when connected via homogeneous excitatory-excitatory (E-E) or inhibitory-inhibitory (I-I) connections. Phase-locked firing is consistently synchronous for either cell type with E-E connections and nearly anti-synchronous with I-I connections. With heterogeneous connections (e.g., excitatory-inhibitory, as might be expected if members of a given population had heterogeneous connections involving intermediate interneurons), networks often settled into phase locking that was either stable or unstable, depending on the order of firing of the two cells in the hybrid network. Our results imply that excitatory SCs, but not inhibitory O-LM interneurons, are capable of synchronizing in phase via monosynaptic mutual connections of the biologically appropriate polarity. Results are largely independent of synaptic strength and synaptic kinetics, implying that our conclusions are robust and largely unaffected by synaptic plasticity.

Dynamics from a time series: can we extract the phase resetting curve from a time series?

Recordings of the membrane potential from a bursting neuron were used to reconstruct the phase curve for that neuron for a limited set of perturbations. These perturbations were inhibitory synaptic conductance pulses able to shift the membrane potential below the most hyperpolarized level attained in the free running mode. The extraction of the phase resetting curve from such a one-dimensional time series requires reconstruction of the periodic activity in the form of a limit cycle attractor. Resetting was found to have two components. In the first component, if the pulse was applied during a burst, the burst was truncated, and the time until the next burst was shortened in a manner predicted by movement normal to the limit cycle. By movement normal to the limit cycle, we mean a switch between two well-defined solution branches of a relaxation-like oscillator in a hysteretic manner enabled by the existence of a singular dominant slow process (variable). In the second component, the onset of the burst was delayed until the end of the hyperpolarizing pulse. Thus, for the pulse amplitudes we studied, resetting was independent of amplitude but increased linearly with pulse duration. The predicted and the experimental phase resetting curves for a pyloric dilator neuron show satisfactory agreement. The method was applied to only one pulse per cycle, but our results suggest it could easily be generalized to accommodate multiple inputs.