Causal information quantification of prominent dynamical features of biological neurons

Neurons tend to fire a spike when they are near a bifurcation from the resting state to spiking activity. It is a delicate balance between noise, dynamic currents and initial condition that determines the phase diagram of neural activity. Many possible ionic mechanisms can be accounted for as the source of spike generation. Moreover, the biophysics and the dynamics behind it can usually be described through a phase diagram that involves membrane voltage versus the activation variable of the ionic channel. In this paper, we present a novel methodology to characterize the dynamics of this system, which takes into account the fine temporal structures of the complex neuronal signals. This allows us to accurately distinguish the most fundamental properties of neurophysiological neurons that were previously described by Izhikevich considering the phase-space trajectory, using a time causal space: statistical complexity versus Fisher information versus Shannon entropy.Instituto de Física de Líquidos y Sistemas Biológico

A simple method for detecting chaos in nature

Chaos, or exponential sensitivity to small perturbations, appears everywhere in nature. Moreover, chaos is predicted to play diverse functional roles in living systems. A method for detecting chaos from empirical measurements should therefore be a key component of the biologist's toolkit. But, classic chaos-detection tools are highly sensitive to measurement noise and break down for common edge cases, making it difficult to detect chaos in domains, like biology, where measurements are noisy. However, newer tools promise to overcome these limitations. Here, we combine several such tools into an automated processing pipeline, and show that our pipeline can detect the presence (or absence) of chaos in noisy recordings, even for difficult edge cases. As a first-pass application of our pipeline, we show that heart rate variability is not chaotic as some have proposed, and instead reflects a stochastic process in both health and disease. Our tool is easy-to-use and freely available

Oscillator-based neuronal modeling for seizure progression investigation and seizure control strategy

The coupled oscillator model has previously been used for the simulation of neuronal activities in in vitro rat hippocampal slice seizure data and the evaluation of seizure suppression algorithms. Each model unit can be described as either an oscillator which can generate action potential spike trains without inputs, or a threshold-based unit. With the change of only one parameter, each unit can either be an oscillator or a threshold-based spiking unit. This would eliminate the need for a new set of equations for each type of unit. Previous analysis has suggested that long kernel duration and imbalance of inhibitory feedback can cause the system to intermittently transition into and out of ictal activities. The state transitions of seizure-like events were investigated here; specifically, how the system excitability may change when the system undergoes transitions in the preictal and postictal processes. Analysis showed that the area of the excitation kernel is positively correlated with the mean firing rate of the ictal activity. The kernel duration is also correlated to the amount of ictal activity. The transition into ictal activity involved the escape from the saddle point foci in the state space trajectory identified by using Newton\u27s method. The ability to accurately anticipate and suppress seizures is an important endeavor that has tremendous impact on improving the quality of lives for epileptic patients. The stimulation studies have suggested that an electrical stimulation strategy that uses the intrinsic high complexity dynamics of the biological system may be more effective in reducing the duration of seizure-like activities in the computer model. In this research, we evaluate this strategy on an in vitro rat hippocampal slice magnesium-free model. Simulated postictal field potential data generated by an oscillator-based hippocampal network model was applied to the CA1 region of the rat hippocampal slices through a multi-electrode array (MEA) system. It was found to suppress and delay the onset of future seizures temporarily. The average inter-seizure time was found to be significantly prolonged after postictal stimulation when compared to the negative control trials and bipolar square wave signals. The result suggests that neural signal-based stimulation related to resetting may be suitable for seizure control in the clinical environment

Mechanisms of Zero-Lag Synchronization in Cortical Motifs

Zero-lag synchronization between distant cortical areas has been observed in a diversity of experimental data sets and between many different regions of the brain. Several computational mechanisms have been proposed to account for such isochronous synchronization in the presence of long conduction delays: Of these, the phenomenon of "dynamical relaying" - a mechanism that relies on a specific network motif - has proven to be the most robust with respect to parameter mismatch and system noise. Surprisingly, despite a contrary belief in the community, the common driving motif is an unreliable means of establishing zero-lag synchrony. Although dynamical relaying has been validated in empirical and computational studies, the deeper dynamical mechanisms and comparison to dynamics on other motifs is lacking. By systematically comparing synchronization on a variety of small motifs, we establish that the presence of a single reciprocally connected pair - a "resonance pair" - plays a crucial role in disambiguating those motifs that foster zero-lag synchrony in the presence of conduction delays (such as dynamical relaying) from those that do not (such as the common driving triad). Remarkably, minor structural changes to the common driving motif that incorporate a reciprocal pair recover robust zero-lag synchrony. The findings are observed in computational models of spiking neurons, populations of spiking neurons and neural mass models, and arise whether the oscillatory systems are periodic, chaotic, noise-free or driven by stochastic inputs. The influence of the resonance pair is also robust to parameter mismatch and asymmetrical time delays amongst the elements of the motif. We call this manner of facilitating zero-lag synchrony resonance-induced synchronization, outline the conditions for its occurrence, and propose that it may be a general mechanism to promote zero-lag synchrony in the brain.Comment: 41 pages, 12 figures, and 11 supplementary figure

Efficiency characterization of a large neuronal network: a causal information approach

When inhibitory neurons constitute about 40% of neurons they could have an important antinociceptive role, as they would easily regulate the level of activity of other neurons. We consider a simple network of cortical spiking neurons with axonal conduction delays and spike timing dependent plasticity, representative of a cortical column or hypercolumn with large proportion of inhibitory neurons. Each neuron fires following a Hodgkin-Huxley like dynamics and it is interconnected randomly to other neurons. The network dynamics is investigated estimating Bandt and Pompe probability distribution function associated to the interspike intervals and taking different degrees of inter-connectivity across neurons. More specifically we take into account the fine temporal ``structures'' of the complex neuronal signals not just by using the probability distributions associated to the inter spike intervals, but instead considering much more subtle measures accounting for their causal information: the Shannon permutation entropy, Fisher permutation information and permutation statistical complexity. This allows us to investigate how the information of the system might saturate to a finite value as the degree of inter-connectivity across neurons grows, inferring the emergent dynamical properties of the system.Comment: 26 pages, 3 Figures; Physica A, in pres

Two-cluster bifurcations in systems of globally pulse-coupled oscillators

For a system of globally pulse-coupled phase-oscillators, we derive conditions for stability of the completely synchronous state and all possible two-cluster states and explain how the different states are naturally connected via bifurcations. The coupling is modeled using the phase-response-curve (PRC), which measures the sensitivity of each oscillator's phase to perturbations. For large systems with a PRC, which turns to zero at the spiking threshold, we are able to find the parameter regions where multiple stable two-cluster states coexist and illustrate this by an example. In addition, we explain how a locally unstable one-cluster state may form an attractor together will its homoclinic connections. This leads to the phenomenon of intermittent, asymptotic synchronization with abating beats away from the perfect synchrony.Comment: 12 pages. 6 figure

Loss of synchrony in an inhibitory network of type-I oscillators

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 (2007) 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. Also, we show that the entire bifurcation structure of the network can be explained by a simple scaling of the STRC (spike- time response curve) amplitude, using a simplified quadratic STRC as an example, and derive the general conditions on the shape of the STRC function that leads to leap-frog firing. Further, 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 (2002) of the order-preserving phase transition map. We show that the extension of STRC to negative values of phase is necessary to predict the response of a model cell to several close non-weak perturbations. This allows us for instance to accurately describe the dynamics of non-weakly coupled network of three model cells. Finally, the phase return map is also extended to the heterogenous network, and is used to analyze both the order-alternating firing and the order-preserving non-zero phase locked state in this case