Sufficient Conditions for Fast Switching Synchronization in Time Varying Network Topologies

In previous work, empirical evidence indicated that a time-varying network could propagate sufficient information to allow synchronization of the sometimes coupled oscillators, despite an instantaneously disconnected topology. We prove here that if the network of oscillators synchronizes for the static time-average of the topology, then the network will synchronize with the time-varying topology if the time-average is achieved sufficiently fast. Fast switching, fast on the time-scale of the coupled oscillators, overcomes the descychnronizing decoherence suggested by disconnected instantaneous networks. This result agrees in spirit with that of where empirical evidence suggested that a moving averaged graph Laplacian could be used in the master-stability function analysis. A new fast switching stability criterion here-in gives sufficiency of a fast-switching network leading to synchronization. Although this sufficient condition appears to be very conservative, it provides new insights about the requirements for synchronization when the network topology is time-varying. In particular, it can be shown that networks of oscillators can synchronize even if at every point in time the frozen-time network topology is insufficiently connected to achieve synchronization.Comment: Submitted to SIAD

Response of electrically coupled spiking neurons: a cellular automaton approach

Experimental data suggest that some classes of spiking neurons in the first layers of sensory systems are electrically coupled via gap junctions or ephaptic interactions. When the electrical coupling is removed, the response function (firing rate {\it vs.} stimulus intensity) of the uncoupled neurons typically shows a decrease in dynamic range and sensitivity. In order to assess the effect of electrical coupling in the sensory periphery, we calculate the response to a Poisson stimulus of a chain of excitable neurons modeled by $n$-state Greenberg-Hastings cellular automata in two approximation levels. The single-site mean field approximation is shown to give poor results, failing to predict the absorbing state of the lattice, while the results for the pair approximation are in good agreement with computer simulations in the whole stimulus range. In particular, the dynamic range is substantially enlarged due to the propagation of excitable waves, which suggests a functional role for lateral electrical coupling. For probabilistic spike propagation the Hill exponent of the response function is $\alpha=1$, while for deterministic spike propagation we obtain $\alpha=1/2$, which is close to the experimental values of the psychophysical Stevens exponents for odor and light intensities. Our calculations are in qualitative agreement with experimental response functions of ganglion cells in the mammalian retina.Comment: 11 pages, 8 figures, to appear in the Phys. Rev.

Dynamic clamp with StdpC software

Dynamic clamp is a powerful method that allows the introduction of artificial electrical components into target cells to simulate ionic conductances and synaptic inputs. This method is based on a fast cycle of measuring the membrane potential of a cell, calculating the current of a desired simulated component using an appropriate model and injecting this current into the cell. Here we present a dynamic clamp protocol using free, fully integrated, open-source software (StdpC, for spike timing-dependent plasticity clamp). Use of this protocol does not require specialist hardware, costly commercial software, experience in real-time operating systems or a strong programming background. The software enables the configuration and operation of a wide range of complex and fully automated dynamic clamp experiments through an intuitive and powerful interface with a minimal initial lead time of a few hours. After initial configuration, experimental results can be generated within minutes of establishing cell recording

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

Network inference - with confidence - from multivariate time series

Networks - collections of interacting elements or nodes - abound in the natural and manmade worlds. For many networks, complex spatiotemporal dynamics stem from patterns of physical interactions unknown to us. To infer these interactions, it is common to include edges between those nodes whose time series exhibit sufficient functional connectivity, typically defined as a measure of coupling exceeding a pre-determined threshold. However, when uncertainty exists in the original network measurements, uncertainty in the inferred network is likely, and hence a statistical propagation-of-error is needed. In this manuscript, we describe a principled and systematic procedure for the inference of functional connectivity networks from multivariate time series data. Our procedure yields as output both the inferred network and a quantification of uncertainty of the most fundamental interest: uncertainty in the number of edges. To illustrate this approach, we apply our procedure to simulated data and electrocorticogram data recorded from a human subject during an epileptic seizure. We demonstrate that the procedure is accurate and robust in both the determination of edges and the reporting of uncertainty associated with that determination.Comment: 12 pages, 7 figures (low resolution), submitte

Impact of Adaptation Currents on Synchronization of Coupled Exponential Integrate-and-Fire Neurons

The ability of spiking neurons to synchronize their activity in a network depends on the response behavior of these neurons as quantified by the phase response curve (PRC) and on coupling properties. The PRC characterizes the effects of transient inputs on spike timing and can be measured experimentally. Here we use the adaptive exponential integrate-and-fire (aEIF) neuron model to determine how subthreshold and spike-triggered slow adaptation currents shape the PRC. Based on that, we predict how synchrony and phase locked states of coupled neurons change in presence of synaptic delays and unequal coupling strengths. We find that increased subthreshold adaptation currents cause a transition of the PRC from only phase advances to phase advances and delays in response to excitatory perturbations. Increased spike-triggered adaptation currents on the other hand predominantly skew the PRC to the right. Both adaptation induced changes of the PRC are modulated by spike frequency, being more prominent at lower frequencies. Applying phase reduction theory, we show that subthreshold adaptation stabilizes synchrony for pairs of coupled excitatory neurons, while spike-triggered adaptation causes locking with a small phase difference, as long as synaptic heterogeneities are negligible. For inhibitory pairs synchrony is stable and robust against conduction delays, and adaptation can mediate bistability of in-phase and anti-phase locking. We further demonstrate that stable synchrony and bistable in/anti-phase locking of pairs carry over to synchronization and clustering of larger networks. The effects of adaptation in aEIF neurons on PRCs and network dynamics qualitatively reflect those of biophysical adaptation currents in detailed Hodgkin-Huxley-based neurons, which underscores the utility of the aEIF model for investigating the dynamical behavior of networks. Our results suggest neuronal spike frequency adaptation as a mechanism synchronizing low frequency oscillations in local excitatory networks, but indicate that inhibition rather than excitation generates coherent rhythms at higher frequencies

The response of a classical Hodgkin–Huxley neuron to an inhibitory input pulse

A population of uncoupled neurons can often be brought close to synchrony by a single strong inhibitory input pulse affecting all neurons equally. This mechanism is thought to underlie some brain rhythms, in particular gamma frequency (30–80 Hz) oscillations in the hippocampus and neocortex. Here we show that synchronization by an inhibitory input pulse often fails for populations of classical Hodgkin–Huxley neurons. Our reasoning suggests that in general, synchronization by inhibitory input pulses can fail when the transition of the target neurons from rest to spiking involves a Hopf bifurcation, especially when inhibition is shunting, not hyperpolarizing. Surprisingly, synchronization is more likely to fail when the inhibitory pulse is stronger or longer-lasting. These findings have potential implications for the question which neurons participate in brain rhythms, in particular in gamma oscillations

Synapse Clusters Are Preferentially Formed by Synapses with Large Recycling Pool Sizes

Synapses are distributed heterogeneously in neural networks. The relationship between the spatial arrangement of synapses and an individual synapse's structural and functional features remains to be elucidated. Here, we examined the influence of the number of adjacent synapses on individual synaptic recycling pool sizes. When measuring the discharge of the styryl dye FM1–43 from electrically stimulated synapses in rat hippocampal tissue cultures, a strong positive correlation between the number of neighbouring synapses and recycling vesicle pool sizes was observed. Accordingly, vesicle-rich synapses were found to preferentially reside next to neighbours with large recycling pool sizes. Although these synapses with large recycling pool sizes were rare, they were densely arranged and thus exhibited a high amount of release per volume. To consolidate these findings, functional terminals were marked by live-cell antibody staining with anti-synaptotagmin-1-cypHer or overexpression of synaptopHluorin. Analysis of synapse distributions in these systems confirmed the results obtained with FM 1–43. Our findings support the idea that clustering of synapses with large recycling pool sizes is a distinct developmental feature of newly formed neural networks and may contribute to functional plasticity

The Mechanism of Abrupt Transition between Theta and Hyper-Excitable Spiking Activity in Medial Entorhinal Cortex Layer II Stellate Cells

Recent studies have shown that stellate cells (SCs) of the medial entorhinal cortex become hyper-excitable in animal models of temporal lobe epilepsy. These studies have also demonstrated the existence of recurrent connections among SCs, reduced levels of recurrent inhibition in epileptic networks as compared to control ones, and comparable levels of recurrent excitation among SCs in both network types. In this work, we investigate the biophysical and dynamic mechanism of generation of the fast time scale corresponding to hyper-excitable firing and the transition between theta and fast firing frequency activity in SCs. We show that recurrently connected minimal networks of SCs exhibit abrupt, threshold-like transition between theta and hyper-excitable firing frequencies as the result of small changes in the maximal synaptic (AMPAergic) conductance. The threshold required for this transition is modulated by synaptic inhibition. Similar abrupt transition between firing frequency regimes can be observed in single, self-coupled SCs, which represent a network of recurrently coupled neurons synchronized in phase, but not in synaptically isolated SCs as the result of changes in the levels of the tonic drive. Using dynamical systems tools (phase-space analysis), we explain the dynamic mechanism underlying the genesis of the fast time scale and the abrupt transition between firing frequency regimes, their dependence on the intrinsic SC's currents and synaptic excitation. This abrupt transition is mechanistically different from others observed in similar networks with different cell types. Most notably, there is no bistability involved. ‘In vitro’ experiments using single SCs self-coupled with dynamic clamp show the abrupt transition between firing frequency regimes, and demonstrate that our theoretical predictions are not an artifact of the model. In addition, these experiments show that high-frequency firing is burst-like with a duration modulated by an M-current