Dynamics of globally coupled inhibitory neurons with heterogeneity

Synchrony in Networks of Type 2 Interneurons Is More Robust to Noise with Hyperpolarizing Inhibition Compared to Shunting Inhibition in Both the Stochastic Population Oscillator and the Coupled Oscillator Regimes

Synchronization in the gamma band (25-150 Hz) is mediated by PV+ inhibitory interneurons, and evidence is accumulating for the essential role of gamma oscillations in cognition. Oscillations can arise in inhibitory networks via synaptic interactions between individual oscillatory neurons (mean-driven) or via strong recurrent inhibition that destabilizes the stationary background firing rate in the fluctuation-driven balanced state, causing an oscillation in the population firing rate. Previous theoretical work focused on model neurons with Hodgkin's Type 1 excitability (integrators) connected by current-based synapses. Here we show that networks comprised of simple Type 2 oscillators (resonators) exhibit a supercritical Hopf bifurcation between synchrony and asynchrony and a gradual transition via cycle skipping from coupled oscillators to stochastic population oscillator (SPO), as previously shown for Type 1. We extended our analysis to homogeneous networks with conductance rather than current based synapses and found that networks with hyperpolarizing inhibitory synapses were more robust to noise than those with shunting synapses, both in the coupled oscillator and SPO regime. Assuming that reversal potentials are uniformly distributed between shunting and hyperpolarized values, as observed in one experimental study, converting synapses to purely hyperpolarizing favored synchrony in all cases, whereas conversion to purely shunting synapses made synchrony less robust except at very high conductance strengths. In mature neurons the synaptic reversal potential is controlled by chloride cotransporters that control the intracellular concentrations of chloride and bicarbonate ions, suggesting these transporters as a potential therapeutic target to enhance gamma synchrony and cognition.

Neural heterogeneity controls computations in spiking neural networks

The brain is composed of complex networks of interacting neurons that express considerable heterogeneity in their physiology and spiking characteristics. How does this neural heterogeneity influence macroscopic neural dynamics, and how might it contribute to neural computation? In this work, we use a mean-field model to investigate computation in heterogeneous neural networks, by studying how the heterogeneity of cell spiking thresholds affects three key computational functions of a neural population: the gating, encoding, and decoding of neural signals. Our results suggest that heterogeneity serves different computational functions in different cell types. In inhibitory interneurons, varying the degree of spike threshold heterogeneity allows them to gate the propagation of neural signals in a reciprocally coupled excitatory population. Whereas homogeneous interneurons impose synchronized dynamics that narrow the dynamic repertoire of the excitatory neurons, heterogeneous interneurons act as an inhibitory offset while preserving excitatory neuron function. Spike threshold heterogeneity also controls the entrainment properties of neural networks to periodic input, thus affecting the temporal gating of synaptic inputs. Among excitatory neurons, heterogeneity increases the dimensionality of neural dynamics, improving the network's capacity to perform decoding tasks. Conversely, homogeneous networks suffer in their capacity for function generation, but excel at encoding signals via multistable dynamic regimes. Drawing from these findings, we propose intra-cell-type heterogeneity as a mechanism for sculpting the computational properties of local circuits of excitatory and inhibitory spiking neurons, permitting the same canonical microcircuit to be tuned for diverse computational tasks.

Synchrony in Networks of Type 2 Interneurons is More Robust to Noise with Hyperpolarizing Inhibition Compared to Shunting Inhibition in Both the Stochastic Population Oscillator and the Coupled Oscillator Regimes

Synchronization in the gamma band (30-80 Hz) is mediated by PV+ inhibitory interneurons, and evidence is accumulating for the essential role of gamma oscillations in cognition. Oscillations can arise in inhibitory networks via synaptic interactions between individual oscillatory neurons (mean-driven) or via strong recurrent inhibition that destabilizes the stationary background firing rate in the fluctuation-driven balanced state, causing an oscillation in the population firing rate. Previous theoretical work focused on model neurons with Hodgkin's type 1 excitability (integrators) connected by current-based synapses. Here we show that networks comprised of simple type 2 oscillators (resonators) exhibit a supercritical Hopf bifurcation between synchrony and asynchrony and a gradual transition via cycle skipping from coupled oscillators to stochastic population oscillator, as previously shown for type 1. We extended our analysis to homogeneous networks with conductance rather than current based synapses and found that networks with hyperpolarizing inhibitory synapses were more robust to noise than those with shunting synapses, both in the coupled oscillator and stochastic population oscillator regime. Assuming that reversal potentials are uniformly distributed between shunting and hyperpolarized values, as observed in one experimental study, converting synapses to purely hyperpolarizing favored synchrony in all cases, whereas conversion to purely shunting synapses made synchrony less robust except at very high conductance strengths. In mature neurons the synaptic reversal potential is controlled by chloride cotransporters that control the intracellular concentrations of chloride and bicarbonate ions, suggesting these transporters as a potential therapeutic target to enhance gamma synchrony and cognition.

Significance Statement

Brain rhythms in the gamma frequency band (30-80 Hz) depend on the activity of inhibitory interneurons and evidence for a causal role for gamma oscillations in cognitive functions is accumulating. Here we extend previous studies on synchronization mechanisms to interneurons that have an abrupt threshold frequency below which they cannot sustain firing. In addition to current based synapses, we examined inhibitory networks with conductance based synapses. We found that if the reversal potential for inhibition was below the average membrane potential (hyperpolarizing), synchrony was more robust to noise than if the reversal potential was very close to the average potential (shunting). These results have implications for therapies to ameliorate cognitive deficits.

Recognizing intertwined patterns using a network of spiking pattern recognition platforms

Artificial intelligence computing adapted from biology is a suitable platform for the development of intelligent machines by imitating the functional mechanisms of the nervous system in creating high-level activities such as learning, decision making and cognition in today's systems. Here, the concentration is on improvement the cognitive potential of artificial intelligence network with a bio-inspired structure. In this regard, four spiking pattern recognition platforms for recognizing digits and letters of EMNIST, patterns of YALE, and ORL datasets are proposed. All networks are developed based on a similar structure in the input image coding, model of neurons (pyramidal neurons and interneurons) and synapses (excitatory AMPA and inhibitory GABA currents), and learning procedure. Networks 1-4 are trained on Digits, Letters, faces of YALE and ORL, respectively, with the proposed un-supervised, spatial-temporal, and sparse spike-based learning mechanism based on the biological observation of the brain learning. When the networks have reached the highest recognition accuracy in the relevant patterns, the main goal of the article, which is to achieve high-performance pattern recognition system with higher cognitive ability, is followed. The pattern recognition network that is able to detect the combination of multiple patterns which called intertwined patterns has not been discussed yet. Therefore, by integrating four trained spiking pattern recognition platforms in one system configuration, we are able to recognize intertwined patterns. These results are presented for the first time and could be the pioneer of a new generation of pattern recognition networks with a significant ability in smart machines.

Loss of neuronal heterogeneity in epileptogenic human tissue impairs network resilience to sudden changes in synchrony

A myriad of pathological changes associated with epilepsy can be recast as decreases in cell and circuit heterogeneity. We thus propose recontextualizing epileptogenesis as a process where reduction in cellular heterogeneity, in part, renders neural circuits less resilient to seizure. By comparing patch clamp recordings from human layer 5 (L5) cortical pyramidal neurons from epileptogenic and non-epileptogenic tissue, we demonstrate significantly decreased biophysical heterogeneity in seizure-generating areas. Implemented computationally, this renders model neural circuits prone to sudden transitions into synchronous states with increased firing activity, paralleling ictogenesis. This computational work also explains the surprising finding of significantly decreased excitability in the population-activation functions of neurons from epileptogenic tissue. Finally, mathematical analyses reveal a bifurcation structure arising only with low heterogeneity and associated with seizure-like dynamics. Taken together, this work provides experimental, computational, and mathematical support for the theory that ictogenic dynamics accompany a reduction in biophysical heterogeneity.

Dynamical Characteristics of Recurrent Neuronal Networks Are Robust Against Low Synaptic Weight Resolution

The representation of the natural-density, heterogeneous connectivity of neuronal network models at relevant spatial scales remains a challenge for Computational Neuroscience and Neuromorphic Computing. In particular, the memory demands imposed by the vast number of synapses in brain-scale network simulations constitute a major obstacle. Limiting the number resolution of synaptic weights appears to be a natural strategy to reduce memory and compute load. In this study, we investigate the effects of a limited synaptic-weight resolution on the dynamics of recurrent spiking neuronal networks resembling local cortical circuits and develop strategies for minimizing deviations from the dynamics of networks with high-resolution synaptic weights. We mimic the effect of a limited synaptic weight resolution by replacing normally distributed synaptic weights with weights drawn from a discrete distribution, and compare the resulting statistics characterizing firing rates, spike-train irregularity, and correlation coefficients with the reference solution. We show that a naive discretization of synaptic weights generally leads to a distortion of the spike-train statistics. If the weights are discretized such that the mean and the variance of the total synaptic input currents are preserved, the firing statistics remain unaffected for the types of networks considered in this study. For networks with sufficiently heterogeneous in-degrees, the firing statistics can be preserved even if all synaptic weights are replaced by the mean of the weight distribution. We conclude that even for simple networks with non-plastic neurons and synapses, a discretization of synaptic weights can lead to substantial deviations in the firing statistics unless the discretization is performed with care and guided by a rigorous validation process. For the network model used in this study, the synaptic weights can be replaced by low-resolution weights without affecting its macroscopic dynamical characteristics, thereby saving substantial amounts of memory.

Optimal responsiveness and information flow in networks of heterogeneous neurons

Cerebral cortex is characterized by a strong neuron-to-neuron heterogeneity, but it is unclear what consequences this may have for cortical computations, while most computational models consider networks of identical units. Here, we study network models of spiking neurons endowed with heterogeneity, that we treat independently for excitatory and inhibitory neurons. We find that heterogeneous networks are generally more responsive, with an optimal responsiveness occurring for levels of heterogeneity found experimentally in different published datasets, for both excitatory and inhibitory neurons. To investigate the underlying mechanisms, we introduce a mean-field model of heterogeneous networks. This mean-field model captures optimal responsiveness and suggests that it is related to the stability of the spontaneous asynchronous state. The mean-field model also predicts that new dynamical states can emerge from heterogeneity, a prediction which is confirmed by network simulations. Finally we show that heterogeneous networks maximise the information flow in large-scale networks, through recurrent connections. We conclude that neuronal heterogeneity confers different responsiveness to neural networks, which should be taken into account to investigate their information processing capabilities.

Neuronal heterogeneity modulates phase synchronization between unidirectionally coupled populations with excitation-inhibition balance

Several experiments and models have highlighted the importance of neuronal heterogeneity in brain dynamics and function. However, how such a cell-to-cell diversity can affect cortical computation, synchronization, and neuronal communication is still under debate. Previous studies have focused on the effect of neuronal heterogeneity in one neuronal population. Here we are specifically interested in the effect of neuronal variability on the phase relations between two populations, which can be related to different cortical communication hypotheses. It has been recently shown that two spiking neuron populations unidirectionally connected in a sender-receiver configuration can exhibit anticipated synchronization (AS), which is characterized by a negative phase lag. This phenomenon has been reported in electrophysiological data of nonhuman primates and human EEG during a visual discrimination cognitive task. In experiments, the unidirectional coupling could be accessed by Granger causality and can be accompanied by either positive or negative phase difference between cortical areas. Here we propose a model of two coupled populations in which the neuronal heterogeneity can determine the dynamical relation between the sender and the receiver and can reproduce phase relations reported in experiments. Depending on the distribution of parameters characterizing the neuronal firing patterns, the system can exhibit both AS and the usual delayed synchronization regime (DS, with positive phase) as well as a zero-lag synchronization regime and phase bistability between AS and DS. Furthermore, we show that our network can present diversity in their phase relations maintaining the excitation-inhibition balance.

Interaction of Indirect and Hyperdirect Pathways on Synchrony and Tremor-Related Oscillation in the Basal Ganglia

Low-frequency oscillatory activity (3-9 Hz) and increased synchrony in the basal ganglia (BG) are recognized to be crucial for Parkinsonian tremor. However, the dynamical mechanism underlying the tremor-related oscillations still remains unknown. In this paper, the roles of the indirect and hyperdirect pathways on synchronization and tremor-related oscillations are considered based on a modified Hodgkin-Huxley model. Firstly, the effects of indirect and hyperdirect pathways are analysed individually, which show that increased striatal activity to the globus pallidus external (GPe) or strong cortical gamma input to the subthalamic nucleus (STN) is sufficient to promote synchrony and tremor-related oscillations in the BG network. Then, the mutual effects of both pathways are analysed by adjusting the related currents simultaneously. Our results suggest that synchrony and tremor-related oscillations would be strengthened if the current of these two paths are in relative imbalance. And the network tends to be less synchronized and less tremulous when the frequency of cortical input is in the theta band. These findings may provide novel treatments in the cortex and striatum to alleviate symptoms of tremor in Parkinson's disease.

Stability of motor cortex network states during learning-associated neural reorganizations

A substantial reorganization of neural activity and neuron-to-movement relationship in motor cortical circuits accompanies the emergence of reproducible movement patterns during motor learning. Little is known about how this tempest of neural activity restructuring impacts the stability of network states in recurrent cortical circuits. To investigate this issue, we reanalyzed data in which we recorded for 14 days via population calcium imaging the activity of the same neural populations of pyramidal neurons in layer 2/3 and layer 5 of forelimb motor and premotor cortex in mice during the daily learning of a lever-press task. We found that motor cortex network states remained stable with respect to the critical network state during the extensive reorganization of both neural population activity and its relation to lever movement throughout learning. Specifically, layer 2/3 cortical circuits unceasingly displayed robust evidence for operating at the critical network state, a regime that maximizes information capacity and transmission and provides a balance between network robustness and flexibility. In contrast, layer 5 circuits operated away from the critical network state for all 14 days of recording and learning. In conclusion, this result indicates that the wide-ranging malleability of synapses, neurons, and neural connectivity during learning operates within the constraint of a stable and layer-specific network state regarding dynamic criticality, and suggests that different cortical layers operate under distinct constraints because of their specialized goals.NEW & NOTEWORTHY The neural activity reorganizes throughout motor learning, but how this reorganization impacts the stability of network states is unclear. We used two-photon calcium imaging to investigate how the network states in layer 2/3 and layer 5 of forelimb motor and premotor cortex are modulated by motor learning. We show that motor cortex network states are layer-specific and constant regarding criticality during neural activity reorganization, and suggests that layer-specific constraints could be motivated by different functions.

Phasic cholinergic signaling promotes emergence of local gamma rhythms in excitatory-inhibitory networks

Recent experimental results have shown that the detection of cues in behavioral attention tasks relies on transient increases of acetylcholine (ACh) release in frontal cortex and cholinergically driven oscillatory activity in the gamma frequency band (Howe et al. Journal of Neuroscience, 2017, 37, 3215). The cue-induced gamma rhythmic activity requires stimulation of M1 muscarinic receptors. Using biophysical computational modeling, we show that a network of excitatory (E) and inhibitory (I) neurons that initially displays asynchronous firing can generate transient gamma oscillatory activity in response to simulated brief pulses of ACh. ACh effects are simulated as transient modulation of the conductance of an M-type K⁺ current which is blocked by activation of muscarinic receptors and has significant effects on neuronal excitability. The ACh-induced effects on the M current conductance, g_Ks , change network dynamics to promote the emergence of network gamma rhythmicity through a Pyramidal-Interneuronal Network Gamma mechanism. Depending on connectivity strengths between and among E and I cells, gamma activity decays with the simulated g_Ks transient modulation or is sustained in the network after the g_Ks transient has completely dissipated. We investigated the sensitivity of the emergent gamma activity to synaptic strengths, external noise and simulated levels of g_Ks modulation. To address recent experimental findings that cholinergic signaling is likely spatially focused and dynamic, we show that localized g_Ks modulation can induce transient changes of cellular excitability in local subnetworks, subsequently causing population-specific gamma oscillations. These results highlight dynamical mechanisms underlying localization of ACh-driven responses and suggest that spatially localized, cholinergically induced gamma may contribute to selectivity in the processing of competing external stimuli, as occurs in attentional tasks.

Information Transmitted From Bioinspired Neuron-Astrocyte Network Improves Cortical Spiking Network's Pattern Recognition Performance

We trained two spiking neural networks (SNNs), the cortical spiking network (CSN) and the cortical neuron-astrocyte network (CNAN), using a spike-based unsupervised method, on the MNIST and alpha-digit data sets and achieve an accuracy of 96.1% and 77.35%, respectively. We then connected CNAN to CSN by preserving maximum synchronization between them thanks to the concept of prolate spheroidal wave functions (PSWF). As a result, CSN receives additional information from CNAN without retraining. The important outcome is that CSN reaches 70.57% correct classification rate on capital letters without being trained on them. The overall contribution of transfer is 87.47%. We observed that for CSN the classifying neurons that relate to digits 0-9 of the alpha-digit data set are completely supported by the ones that relate to digits 0-9 of the MNIST data set. This means that CSN recognizes the similarity between the digits of the MNIST and alpha-digit data sets and classifies each digit of both data sets in the same class.

Inhibitory Network Bistability Explains Increased Interneuronal Activity Prior to Seizure Onset

Recent experimental literature has revealed that GABAergic interneurons exhibit increased activity prior to seizure onset, alongside additional evidence that such activity is synchronous and may arise abruptly. These findings have led some to hypothesize that this interneuronal activity may serve a causal role in driving the sudden change in brain activity that heralds seizure onset. However, the mechanisms predisposing an inhibitory network toward increased activity, specifically prior to ictogenesis, without a permanent change to inputs to the system remain unknown. We address this question by comparing simulated inhibitory networks containing control interneurons and networks containing hyperexcitable interneurons modeled to mimic treatment with 4-Aminopyridine (4-AP), an agent commonly used to model seizures in vivo and in vitro. Our in silico study demonstrates that model inhibitory networks with 4-AP interneurons are more prone than their control counterparts to exist in a bistable state in which asynchronously firing networks can abruptly transition into synchrony driven by a brief perturbation. This transition into synchrony brings about a corresponding increase in overall firing rate. We further show that perturbations driving this transition could arise in vivo from background excitatory synaptic activity in the cortex. Thus, we propose that bistability explains the increase in interneuron activity observed experimentally prior to seizure via a transition from incoherent to coherent dynamics. Moreover, bistability explains why inhibitory networks containing hyperexcitable interneurons are more vulnerable to this change in dynamics, and how such networks can undergo a transition without a permanent change in the drive. We note that while our comparisons are between networks of control and ictogenic neurons, the conclusions drawn specifically relate to the unusual dynamics that arise prior to seizure, and not seizure onset itself. However, providing a mechanistic explanation for this phenomenon specifically in a pro-ictogenic setting generates experimentally testable hypotheses regarding the role of inhibitory neurons in pre-ictal neural dynamics, and motivates further computational research into mechanisms underlying a newly hypothesized multi-step pathway to seizure initiated by inhibition.

Realistic spiking neural network: Non-synaptic mechanisms improve convergence in cell assembly

Learning in neural networks inspired by brain tissue has been studied for machine learning applications. However, existing works primarily focused on the concept of synaptic weight modulation, and other aspects of neuronal interactions, such as non-synaptic mechanisms, have been neglected. Non-synaptic interaction mechanisms have been shown to play significant roles in the brain, and four classes of these mechanisms can be highlighted: (i) electrotonic coupling; (ii) ephaptic interactions; (iii) electric field effects; and iv) extracellular ionic fluctuations. In this work, we proposed simple rules for learning inspired by recent findings in machine learning adapted to a realistic spiking neural network. We show that the inclusion of non-synaptic interaction mechanisms improves cell assembly convergence. By including extracellular ionic fluctuation represented by the extracellular electrodiffusion in the network, we showed the importance of these mechanisms to improve cell assembly convergence. Additionally, we observed a variety of electrophysiological patterns of neuronal activity, particularly bursting and synchronism when the convergence is improved.

Enabling Large-Scale Simulations With the GENESIS Neuronal Simulator

In this paper, we evaluate the computational performance of the GEneral NEural SImulation System (GENESIS) for large scale simulations of neural networks. While many benchmark studies have been performed for large scale simulations with leaky integrate-and-fire neurons or neuronal models with only a few compartments, this work focuses on higher fidelity neuronal models represented by 50–74 compartments per neuron. After making some modifications to the source code for GENESIS and its parallel implementation, PGENESIS, particularly to improve memory usage, we find that PGENESIS is able to efficiently scale on supercomputing resources to network sizes as large as 9 × 10⁶ neurons with 18 × 10⁹ synapses and 2.2 × 10⁶ neurons with 45 × 10⁹ synapses. The modifications to GENESIS that enabled these large scale simulations have been incorporated into the May 2019 Official Release of PGENESIS 2.4 available for download from the GENESIS web site (genesis-sim.org).

Synaptic Failure Differentially Affects Pattern Formation in Heterogenous Networks

The communication of neurons is primarily maintained by synapses, which play a crucial role in the functioning of the nervous system. Therefore, synaptic failure may critically impair information processing in the brain and may underlie many neurodegenerative diseases. A number of studies have suggested that synaptic failure may preferentially target neurons with high connectivity (i.e., network hubs). As a result, the activity of these highly connected neurons can be significantly affected. It has been speculated that anesthetics regulate conscious state by affecting synaptic transmission at these network hubs and subsequently reducing overall coherence in the network activity. In addition, hubs in cortical networks are shown to be more vulnerable to amyloid deposition because of their higher activity within the network, causing decrease in coherence patterns and eventually Alzheimer’s disease (AD). Here, we investigate how synaptic failure can affect spatio-temporal dynamics of scale free networks, having a power law scaling of number of connections per neuron – a relatively few neurons (hubs) with a lot of emanating or incoming connections and many cells with low connectivity. We studied two types of synaptic failure: activity-independent and targeted, activity-dependent synaptic failure. We defined scale-free network structures based on the dominating direction of the connections at the hub neurons: incoming and outgoing. We found that the two structures have significantly different dynamical properties. We show that synaptic failure may not only lead to the loss of coherence but unintuitively also can facilitate its emergence. We show that this is because activity-dependent synaptic failure homogenizes the activity levels in the network creating a dynamical substrate for the observed coherence increase. Obtained results may lead to better understanding of changes in large-scale pattern formation during progression of neuro-degenerative diseases targeting synaptic transmission.

How stimulation frequency and intensity impact on the long-lasting effects of coordinated reset stimulation

Several brain diseases are characterized by abnormally strong neuronal synchrony. Coordinated Reset (CR) stimulation was computationally designed to specifically counteract abnormal neuronal synchronization processes by desynchronization. In the presence of spike-timing-dependent plasticity (STDP) this may lead to a decrease of synaptic excitatory weights and ultimately to an anti-kindling, i.e. unlearning of abnormal synaptic connectivity and abnormal neuronal synchrony. The long-lasting desynchronizing impact of CR stimulation has been verified in pre-clinical and clinical proof of concept studies. However, as yet it is unclear how to optimally choose the CR stimulation frequency, i.e. the repetition rate at which the CR stimuli are delivered. This work presents the first computational study on the dependence of the acute and long-term outcome on the CR stimulation frequency in neuronal networks with STDP. For this purpose, CR stimulation was applied with Rapidly Varying Sequences (RVS) as well as with Slowly Varying Sequences (SVS) in a wide range of stimulation frequencies and intensities. Our findings demonstrate that acute desynchronization, achieved during stimulation, does not necessarily lead to long-term desynchronization after cessation of stimulation. By comparing the long-term effects of the two different CR protocols, the RVS CR stimulation turned out to be more robust against variations of the stimulation frequency. However, SVS CR stimulation can obtain stronger anti-kindling effects. We revealed specific parameter ranges that are favorable for long-term desynchronization. For instance, RVS CR stimulation at weak intensities and with stimulation frequencies in the range of the neuronal firing rates turned out to be effective and robust, in particular, if no closed loop adaptation of stimulation parameters is (technically) available. From a clinical standpoint, this may be relevant in the context of both invasive as well as non-invasive CR stimulation.

Short-Term Dosage Regimen for Stimulation-Induced Long-Lasting Desynchronization

In this paper, we computationally generate hypotheses for dose-finding studies in the context of desynchronizing neuromodulation techniques. Abnormally strong neuronal synchronization is a hallmark of several brain disorders. Coordinated Reset (CR) stimulation is a spatio-temporally patterned stimulation technique that specifically aims at disrupting abnormal neuronal synchrony. In networks with spike-timing-dependent plasticity CR stimulation may ultimately cause an anti-kindling, i.e., an unlearning of abnormal synaptic connectivity and neuronal synchrony. This long-lasting desynchronization was theoretically predicted and verified in several pre-clinical and clinical studies. We have shown that CR stimulation with rapidly varying sequences (RVS) robustly induces an anti-kindling at low intensities e.g., if the CR stimulation frequency (i.e., stimulus pattern repetition rate) is in the range of the frequency of the neuronal oscillation. In contrast, CR stimulation with slowly varying sequences (SVS) turned out to induce an anti-kindling more strongly, but less robustly with respect to variations of the CR stimulation frequency. Motivated by clinical constraints and inspired by the spacing principle of learning theory, in this computational study we propose a short-term dosage regimen that enables a robust anti-kindling effect of both RVS and SVS CR stimulation, also for those parameter values where RVS and SVS CR stimulation previously turned out to be ineffective. Intriguingly, for the vast majority of parameter values tested, spaced multishot CR stimulation with demand-controlled variation of stimulation frequency and intensity caused a robust and pronounced anti-kindling. In contrast, spaced CR stimulation with fixed stimulation parameters as well as singleshot CR stimulation of equal integral duration failed to improve the stimulation outcome. In the model network under consideration, our short-term dosage regimen enables to robustly induce long-term desynchronization at comparably short stimulation duration and low integral stimulation duration. Currently, clinical proof of concept is available for deep brain CR stimulation for Parkinson's therapy and acoustic CR stimulation for tinnitus therapy. Promising first in human data is available for vibrotactile CR stimulation for Parkinson's treatment. For the clinical development of these treatments it is mandatory to perform dose-finding studies to reveal optimal stimulation parameters and dosage regimens. Our findings can straightforwardly be tested in human dose-finding studies.

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.

Coding of time-dependent stimuli in homogeneous and heterogeneous neural populations

We compare the information transmission of a time-dependent signal by two types of uncoupled neuron populations that differ in their sources of variability: i) a homogeneous population whose units receive independent noise and ii) a deterministic heterogeneous population, where each unit exhibits a different baseline firing rate ('disorder'). Our criterion for making both sources of variability quantitatively comparable is that the interspike-interval distributions are identical for both systems. Numerical simulations using leaky integrate-and-fire neurons unveil that a non-zero amount of both noise or disorder maximizes the encoding efficiency of the homogeneous and heterogeneous system, respectively, as a particular case of suprathreshold stochastic resonance. Our findings thus illustrate that heterogeneity can render similarly profitable effects for neuronal populations as dynamic noise. The optimal noise/disorder depends on the system size and the properties of the stimulus such as its intensity or cutoff frequency. We find that weak stimuli are better encoded by a noiseless heterogeneous population, whereas for strong stimuli a homogeneous population outperforms an equivalent heterogeneous system up to a moderate noise level. Furthermore, we derive analytical expressions of the coherence function for the cases of very strong noise and of vanishing intrinsic noise or heterogeneity, which predict the existence of an optimal noise intensity. Our results show that, depending on the type of signal, noise as well as heterogeneity can enhance the encoding performance of neuronal populations.

Synchrony and pattern formation of coupled genetic oscillators on a chip of artificial cells

Understanding how biochemical networks lead to large-scale nonequilibrium self-organization and pattern formation in life is a major challenge, with important implications for the design of programmable synthetic systems. Here, we assembled cell-free genetic oscillators in a spatially distributed system of on-chip DNA compartments as artificial cells, and measured reaction-diffusion dynamics at the single-cell level up to the multicell scale. Using a cell-free gene network we programmed molecular interactions that control the frequency of oscillations, population variability, and dynamical stability. We observed frequency entrainment, synchronized oscillatory reactions and pattern formation in space, as manifestation of collective behavior. The transition to synchrony occurs as the local coupling between compartments strengthens. Spatiotemporal oscillations are induced either by a concentration gradient of a diffusible signal, or by spontaneous symmetry breaking close to a transition from oscillatory to nonoscillatory dynamics. This work offers design principles for programmable biochemical reactions with potential applications to autonomous sensing, distributed computing, and biomedical diagnostics.

CA1 pyramidal cells have diverse biophysical properties, affected by development, experience, and aging

Neuron types (e.g., pyramidal cells) within one area of the brain are often considered homogeneous, despite variability in their biophysical properties. Here we review literature demonstrating variability in the electrical activity of CA1 hippocampal pyramidal cells (PCs), including responses to somatic current injection, synaptic stimulation, and spontaneous network-related activity. In addition, we describe how responses of CA1 PCs vary with development, experience, and aging, and some of the underlying ionic currents responsible. Finally, we suggest directions that may be the most impactful in expanding this knowledge, including the use of text and data mining to systematically study cellular heterogeneity in more depth; dynamical systems theory to understand and potentially classify neuron firing patterns; and mathematical modeling to study the interaction between cellular properties and network output. Our goals are to provide a synthesis of the literature for experimentalists studying CA1 PCs, to give theorists an idea of the rich diversity of behaviors models may need to reproduce to accurately represent these cells, and to provide suggestions for future research.

Boundary in the dynamic phase of globally coupled oscillatory and excitable units

There is a crucial boundary between dynamic phase 1 and dynamic phase 2 of globally coupled oscillatory and excitable units, where the mean field is constant and oscillates in the former and the latter, respectively. This boundary is theoretically derived here for a large population of dynamical units, each having only a phase variable, where it is assumed that both the coupling strength and the distribution width of bifurcation parameters are equally small. This theory, which is applicable only if all or most of the units are intrinsically oscillatory, is confirmed to agree with simulation results for two different distribution densities.

Inverse stochastic resonance in networks of spiking neurons

Inverse Stochastic Resonance (ISR) is a phenomenon in which the average spiking rate of a neuron exhibits a minimum with respect to noise. ISR has been studied in individual neurons, but here, we investigate ISR in scale-free networks, where the average spiking rate is calculated over the neuronal population. We use Hodgkin-Huxley model neurons with channel noise (i.e., stochastic gating variable dynamics), and the network connectivity is implemented via electrical or chemical connections (i.e., gap junctions or excitatory/inhibitory synapses). We find that the emergence of ISR depends on the interplay between each neuron's intrinsic dynamical structure, channel noise, and network inputs, where the latter in turn depend on network structure parameters. We observe that with weak gap junction or excitatory synaptic coupling, network heterogeneity and sparseness tend to favor the emergence of ISR. With inhibitory coupling, ISR is quite robust. We also identify dynamical mechanisms that underlie various features of this ISR behavior. Our results suggest possible ways of experimentally observing ISR in actual neuronal systems.

Pallidostriatal Projections Promote β Oscillations in a Dopamine-Depleted Biophysical Network Model

In the basal ganglia, focused rhythmicity is an important feature of network activity at certain stages of motor processing. In disease, however, the basal ganglia develop amplified rhythmicity. Here, we demonstrate how the cellular architecture and network dynamics of an inhibitory loop in the basal ganglia yield exaggerated synchrony and locking to β oscillations, specifically in the dopamine-depleted state. A key component of this loop is the pallidostriatal pathway, a well-characterized anatomical projection whose function has long remained obscure. We present a synaptic characterization of this pathway in mice and incorporate these data into a computational model that we use to investigate its influence over striatal activity under simulated healthy and dopamine-depleted conditions. Our model predicts that the pallidostriatal pathway influences striatal output preferentially during periods of synchronized activity within GPe. We show that, under dopamine-depleted conditions, this effect becomes a key component of a positive feedback loop between the GPe and striatum that promotes synchronization and rhythmicity. Our results generate novel predictions about the role of the pallidostriatal pathway in shaping basal ganglia activity in health and disease.

SIGNIFICANCE STATEMENT

This work demonstrates that functional connections from the globus pallidus externa (GPe) to striatum are substantially stronger onto fast-spiking interneurons (FSIs) than onto medium spiny neurons. Our circuit model suggests that when GPe spikes are synchronous, this pallidostriatal pathway causes synchronous FSI activity pauses, which allow a transient window of disinhibition for medium spiny neurons. In simulated dopamine-depletion, this GPe-FSI activity is necessary for the emergence of strong synchronization and the amplification and propagation of β oscillations, which are a hallmark of parkinsonian circuit dysfunction. These results suggest that GPe may play a central role in propagating abnormal circuit activity to striatum, which in turn projects to downstream basal ganglia structures. These findings warrant further exploration of GPe as a target for interventions for Parkinson's disease.

Thalamic amplification of cortical connectivity sustains attentional control

While interactions between the thalamus and cortex are critical for cognitive function^1–3, the exact contribution of the thalamus to these interactions is often unclear. Recent studies have shown diverse connectivity patterns across the thalamus ^4,5, but whether this diversity translates to thalamic functions beyond relaying information to or between cortical regions⁶ is unknown. Here, by investigating prefrontal cortical (PFC) representation of two rules used to guide attention, we find that the mediodorsal thalamus (MD) sustains these representations without relaying categorical information. Specifically, MD input amplifies local PFC connectivity, enabling rule-specific neural sequences to emerge and thereby maintain rule representations. Consistent with this notion, broadly enhancing PFC excitability diminishes rule specificity and behavioral performance, while enhancing MD excitability improves both. Overall, our results define a previously unknown principle in neuroscience; thalamic control of functional cortical connectivity. This function indicates that the thalamus plays much more central roles in cognition than previously thought.

Spike timing-dependent plasticity induces non-trivial topology in the brain

We study the capacity of Hodgkin-Huxley neuron in a network to change temporarily or permanently their connections and behavior, the so called spike timing-dependent plasticity (STDP), as a function of their synchronous behavior. We consider STDP of excitatory and inhibitory synapses driven by Hebbian rules. We show that the final state of networks evolved by a STDP depend on the initial network configuration. Specifically, an initial all-to-all topology evolves to a complex topology. Moreover, external perturbations can induce co-existence of clusters, those whose neurons are synchronous and those whose neurons are desynchronous. This work reveals that STDP based on Hebbian rules leads to a change in the direction of the synapses between high and low frequency neurons, and therefore, Hebbian learning can be explained in terms of preferential attachment between these two diverse communities of neurons, those with low-frequency spiking neurons, and those with higher-frequency spiking neurons.

Understanding spatial and temporal patterning of astrocyte calcium transients via interactions between network transport and extracellular diffusion

Astrocytes form interconnected networks in the brain and communicate via calcium signaling. We investigate how modes of coupling between astrocytes influence the spatio-temporal patterns of calcium signaling within astrocyte networks and specifically how these network interactions promote coordination within this group of cells. To investigate these complex phenomena, we study reduced cultured networks of astrocytes and neurons. We image the spatial temporal patterns of astrocyte calcium activity and quantify how perturbing the coupling between astrocytes influences astrocyte activity patterns. To gain insight into the pattern formation observed in these cultured networks, we compare the experimentally observed calcium activity patterns to the patterns produced by a reduced computational model, where we represent astrocytes as simple units that integrate input through two mechanisms: gap junction coupling (network transport) and chemical release (extracellular diffusion). We examine the activity patterns in the simulated astrocyte network and their dependence upon these two coupling mechanisms. We find that gap junctions and extracellular chemical release interact in astrocyte networks to modulate the spatiotemporal patterns of their calcium dynamics. We show agreement between the computational and experimental findings, which suggests that the complex global patterns can be understood as a result of simple local coupling mechanisms.

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.

Bifurcation and scaling at the aging transition boundary in globally coupled excitable and oscillatory units

Following a previous paper [Phys. Rev. E 88, 052907 (2013)PLEEE81539-375510.1103/PhysRevE.88.052907], we study in detail the mechanism of aging transition in globally and diffusively coupled excitable and oscillatory units. Here two of the three models taken up in the earlier work are used, each composed of a large number of units with their bifurcation parameters forming a uniform distribution. The control parameters are the coupling strength and the average of the bifurcation parameters. The present work is mostly devoted to a region of the phase diagram near the aging transition boundary with the coupling strength greater than its critical value for the onset of bistability and hysteresis. The bifurcation structure of each system at the aging transition boundary is investigated theoretically as well as numerically. Moreover, we show that critical scaling laws of order parameters S and M used in the previous paper are different depending on which region of the coupling strength to be chosen.

Coarse-Grained Clustering Dynamics of Heterogeneously Coupled Neurons

The formation of oscillating phase clusters in a network of identical Hodgkin-Huxley neurons is studied, along with their dynamic behavior. The neurons are synaptically coupled in an all-to-all manner, yet the synaptic coupling characteristic time is heterogeneous across the connections. In a network of N neurons where this heterogeneity is characterized by a prescribed random variable, the oscillatory single-cluster state can transition-through [Formula: see text] (possibly perturbed) period-doubling and subsequent bifurcations-to a variety of multiple-cluster states. The clustering dynamic behavior is computationally studied both at the detailed and the coarse-grained levels, and a numerical approach that can enable studying the coarse-grained dynamics in a network of arbitrarily large size is suggested. Among a number of cluster states formed, double clusters, composed of nearly equal sub-network sizes are seen to be stable; interestingly, the heterogeneity parameter in each of the double-cluster components tends to be consistent with the random variable over the entire network: Given a double-cluster state, permuting the dynamical variables of the neurons can lead to a combinatorially large number of different, yet similar "fine" states that appear practically identical at the coarse-grained level. For weak heterogeneity we find that correlations rapidly develop, within each cluster, between the neuron's "identity" (its own value of the heterogeneity parameter) and its dynamical state. For single- and double-cluster states we demonstrate an effective coarse-graining approach that uses the Polynomial Chaos expansion to succinctly describe the dynamics by these quickly established "identity-state" correlations. This coarse-graining approach is utilized, within the equation-free framework, to perform efficient computations of the neuron ensemble dynamics.

Emergent gamma synchrony in all-to-all interneuronal networks

We investigate the emergence of in-phase synchronization in a heterogeneous network of coupled inhibitory interneurons in the presence of spike timing dependent plasticity (STDP). Using a simple network of two mutually coupled interneurons (2-MCI), we first study the effects of STDP on in-phase synchronization. We demonstrate that, with STDP, the 2-MCI network can evolve to either a state of stable 1:1 in-phase synchronization or exhibit multiple regimes of higher order synchronization states. We show that the emergence of synchronization induces a structural asymmetry in the 2-MCI network such that the synapses onto the high frequency firing neurons are potentiated, while those onto the low frequency firing neurons are de-potentiated, resulting in the directed flow of information from low frequency firing neurons to high frequency firing neurons. Finally, we demonstrate that the principal findings from our analysis of the 2-MCI network contribute to the emergence of robust synchronization in the Wang-Buzsaki network (Wang and Buzsáki, 1996) of all-to-all coupled inhibitory interneurons (100-MCI) for a significantly larger range of heterogeneity in the intrinsic firing rate of the neurons in the network. We conclude that STDP of inhibitory synapses provide a viable mechanism for robust neural synchronization.

A positive feedback at the cellular level promotes robustness and modulation at the circuit level

This article highlights the role of a positive feedback gating mechanism at the cellular level in the robustness and modulation properties of rhythmic activities at the circuit level. The results are presented in the context of half-center oscillators, which are simple rhythmic circuits composed of two reciprocally connected inhibitory neuronal populations. Specifically, we focus on rhythms that rely on a particular excitability property, the postinhibitory rebound, an intrinsic cellular property that elicits transient membrane depolarization when released from hyperpolarization. Two distinct ionic currents can evoke this transient depolarization: a hyperpolarization-activated cation current and a low-threshold T-type calcium current. The presence of a slow activation is specific to the T-type calcium current and provides a slow positive feedback at the cellular level that is absent in the cation current. We show that this slow positive feedback is required to endow the network rhythm with physiological modulation and robustness properties. This study thereby identifies an essential cellular property to be retained at the network level in modeling network robustness and modulation.

Synergistic effect of repulsive inhibition in synchronization of excitatory networks

We show that the addition of pairwise repulsive inhibition to excitatory networks of bursting neurons induces synchrony, in contrast to one's expectations. Through stability analysis, we reveal the mechanism underlying this purely synergistic phenomenon and demonstrate that it originates from the transition between different types of bursting, caused by excitatory-inhibitory synaptic coupling. This effect is generic and observed in different models of bursting neurons and fast synaptic interactions. We also find a universal scaling law for the synchronization stability condition for large networks in terms of the number of excitatory and inhibitory inputs each neuron receives, regardless of the network size and topology. This general law is in sharp contrast with linearly coupled networks with positive (attractive) and negative (repulsive) coupling where the placement and structure of negative connections heavily affect synchronization.

Measuring predictability of autonomous network transitions into bursting dynamics

Understanding spontaneous transitions between dynamical modes in a network is of significant importance. These transitions may separate pathological and normal functions of the brain. In this paper, we develop a set of measures that, based on spatio-temporal features of network activity, predict autonomous network transitions from asynchronous to synchronous dynamics under various conditions. These metrics quantify spike-timing distributions within a narrow time window as a function of the relative location of the active neurons. We applied these metrics to investigate the properties of these transitions in excitatory-only and excitatory-and-inhibitory networks and elucidate how network topology, noise level, and cellular heterogeneity affect both the reliability and the timeliness of the predictions. The developed measures can be calculated in real time and therefore potentially applied in clinical situations.

Synchronization properties of heterogeneous neuronal networks with mixed excitability type

We study the synchronization of neuronal networks with dynamical heterogeneity, showing that network structures with the same propensity for synchronization (as quantified by master stability function analysis) may develop dramatically different synchronization properties when heterogeneity is introduced with respect to neuronal excitability type. Specifically, we investigate networks composed of neurons with different types of phase response curves (PRCs), which characterize how oscillating neurons respond to excitatory perturbations. Neurons exhibiting type 1 PRC respond exclusively with phase advances, while neurons exhibiting type 2 PRC respond with either phase delays or phase advances, depending on when the perturbation occurs. We find that Watts-Strogatz small world networks transition to synchronization gradually as the proportion of type 2 neurons increases, whereas scale-free networks may transition gradually or rapidly, depending upon local correlations between node degree and excitability type. Random placement of type 2 neurons results in gradual transition to synchronization, whereas placement of type 2 neurons as hubs leads to a much more rapid transition, showing that type 2 hub cells easily "hijack" neuronal networks to synchronization. These results underscore the fact that the degree of synchronization observed in neuronal networks is determined by a complex interplay between network structure and the dynamical properties of individual neurons, indicating that efforts to recover structural connectivity from dynamical correlations must in general take both factors into account.

Dynamics and Reliability of Bistable Neurons Driven with Time-Dependent Stimuli

The reliability of a spiking neuron depends on the frequency content of the driving input signal. Previous studies have shown that well above threshold, regularly firing neurons generate reliable responses when the input signal resonates with the firing frequency of the cell. Instead, well below threshold, reliable responses are obtained when the input frequency resonates with the subthreshold oscillations of the neuron. Previous theories, however, provide no clear prediction for the input frequency giving rise to maximally reliable spiking at threshold, which is probably the most relevant firing regime in mammalian cortex under physiological conditions. In particular, when the firing onset is governed by a subcritical Hopf bifurcation, the frequency of subthreshold oscillations often differs from the firing rate at threshold. The predictions of previous studies, hence, cannot be smoothly merged at threshold. Here we explore the behavior of reliability in bistable neurons near threshold using three types of driving stimuli: constant, periodic, and stochastic. We find that the two natural frequencies of the system, associated with the two coexisting attractors, provide a rich variety of possible locking modes with the external signal. Reliability is determined by the sensitivity to noise of each locking mode and by the transition probabilities between modes. Noise increases the amount of spike time jitter, and minimal jitter is obtained for input frequencies coinciding with the suprathreshold firing rate of the cell. In addition, noise may either enhance or inhibit transitions between the two attractors, depending on the input frequency. The dual role played by noise in bistable systems implies that reliability is determined by a delicate balance between spike time jitter and the rate of transitions between attractors.

Differential effects of excitatory and inhibitory heterogeneity on the gain and asynchronous state of sparse cortical networks

Recent experimental and theoretical studies have highlighted the importance of cell-to-cell differences in the dynamics and functions of neural networks, such as in different types of neural coding or synchronization. It is still not known, however, how neural heterogeneity can affect cortical computations, or impact the dynamics of typical cortical circuits constituted of sparse excitatory and inhibitory networks. In this work, we analytically and numerically study the dynamics of a typical cortical circuit with a certain level of neural heterogeneity. Our circuit includes realistic features found in real cortical populations, such as network sparseness, excitatory, and inhibitory subpopulations of neurons, and different cell-to-cell heterogeneities for each type of population in the system. We find highly differentiated roles for heterogeneity, depending on the subpopulation in which it is found. In particular, while heterogeneity among excitatory neurons non-linearly increases the mean firing rate and linearizes the f-I curves, heterogeneity among inhibitory neurons may decrease the network activity level and induces divisive gain effects in the f-I curves of the excitatory cells, providing an effective gain control mechanism to influence information flow. In addition, we compute the conditions for stability of the network activity, finding that the synchronization onset is robust to inhibitory heterogeneity, but it shifts to lower input levels for higher excitatory heterogeneity. Finally, we provide an extension of recently reported heterogeneity-induced mechanisms for signal detection under rate coding, and we explore the validity of our findings when multiple sources of heterogeneity are present. These results allow for a detailed characterization of the role of neural heterogeneity in asynchronous cortical networks.

Dynamical robustness of coupled heterogeneous oscillators

We study tolerance of dynamic behavior in networks of coupled heterogeneous oscillators to deterioration of the individual oscillator components. As the deterioration proceeds with reduction in dynamic behavior of the oscillators, an order parameter evaluating the level of global oscillation decreases and then vanishes at a certain critical point. We present a method to analytically derive a general formula for this critical point and an approximate formula for the order parameter in the vicinity of the critical point in networks of coupled Stuart-Landau oscillators. Using the critical point as a measure for dynamical robustness of oscillator networks, we show that the more heterogeneous the oscillator components are, the more robust the oscillatory behavior of the network is to the component deterioration. This property is confirmed also in networks of Morris-Lecar neuron models coupled through electrical synapses. Our approach could provide a useful framework for theoretically understanding the role of population heterogeneity in robustness of biological networks.

SYNCHRONIZATION IN A RING OF UNIDIRECTIONALLY COUPLED FITZHUGH–NAGUMO NEURONS

Synchronization behavior of an ensemble of unidirectionally coupled neurons with a constant input is investigated. Chemical synapses are considered for coupling. Each neuron is also considered to be exposed to a self-delayed feedback. The synchronization phenomenon is analyzed by the error dynamics of the response trajectories of the system. The effect of various model parameters e.g. coupling strength, feedback gain and time delay, on synchronization is also investigated and a measure of synchrony is computed in each cases. It is shown that the synchronization is not only achieved by increasing the coupling strength, the system also required to have a suitable feedback gain and time delay for synchrony. Robustness of the parameters for synchrony is verified for larger systems.

Reliable neuronal systems: the importance of heterogeneity

For every engineer it goes without saying: in order to build a reliable system we need components that consistently behave precisely as they should. It is also well known that neurons, the building blocks of brains, do not satisfy this constraint. Even neurons of the same type come with huge variances in their properties and these properties also vary over time. Synapses, the connections between neurons, are highly unreliable in forwarding signals. In this paper we argue that both these fact add variance to neuronal processes, and that this variance is not a handicap of neural systems, but that instead predictable and reliable functional behavior of neural systems depends crucially on this variability. In particular, we show that higher variance allows a recurrently connected neural population to react more sensitively to incoming signals, and processes them faster and more energy efficient. This, for example, challenges the general assumption that the intrinsic variability of neurons in the brain is a defect that has to be overcome by synaptic plasticity in the process of learning.

High-frequency stimulation of excitable cells and networks

High-frequency (HF) stimulation has been shown to block conduction in excitable cells including neurons and cardiac myocytes. However, the precise mechanisms underlying conduction block are unclear. Using a multi-scale method, the influence of HF stimulation is investigated in the simplified FitzhHugh-Nagumo and biophysically-detailed Hodgkin-Huxley models. In both models, HF stimulation alters the amplitude and frequency of repetitive firing in response to a constant applied current and increases the threshold to evoke a single action potential in response to a brief applied current pulse. Further, the excitable cells cannot evoke a single action potential or fire repetitively above critical values for the HF stimulation amplitude. Analytical expressions for the critical values and thresholds are determined in the FitzHugh-Nagumo model. In the Hodgkin-Huxley model, it is shown that HF stimulation alters the dynamics of ionic current gating, shifting the steady-state activation, inactivation, and time constant curves, suggesting several possible mechanisms for conduction block. Finally, we demonstrate that HF stimulation of a network of neurons reduces the electrical activity firing rate, increases network synchronization, and for a sufficiently large HF stimulation, leads to complete electrical quiescence. In this study, we demonstrate a novel approach to investigate HF stimulation in biophysically-detailed ionic models of excitable cells, demonstrate possible mechanisms for HF stimulation conduction block in neurons, and provide insight into the influence of HF stimulation on neural networks.

Onset of dynamic activity in globally coupled excitable and oscillatory units

The onset of dynamic activity is studied in heterogeneous populations of globally coupled units such that each unit is either excitable or oscillatory depending on its parameter. By varying the mean of the parameter distribution as well as coupling strength, we show that all or part of the populations studied have common features: similar phase diagrams, hysteresis near or at the onset of dynamic activity, and a fractional power law obeyed by an order parameter. A simplified model is proposed to explain these results.

Inferring network properties of cortical neurons with synaptic coupling and parameter dispersion

Computational models at different space-time scales allow us to understand the fundamental mechanisms that govern neural processes and relate uniquely these processes to neuroscience data. In this work, we propose a novel neurocomputational unit (a mesoscopic model which tell us about the interaction between local cortical nodes in a large scale neural mass model) of bursters that qualitatively captures the complex dynamics exhibited by a full network of parabolic bursting neurons. We observe that the temporal dynamics and fluctuation of mean synaptic action term exhibits a high degree of correlation with the spike/burst activity of our population. With heterogeneity in the applied drive and mean synaptic coupling derived from fast excitatory synapse approximations we observe long term behavior in our population dynamics such as partial oscillations, incoherence, and synchrony. In order to understand the origin of multistability at the population level as a function of mean synaptic coupling and heterogeneity in the firing rate threshold we employ a simple generative model for parabolic bursting recently proposed by Ghosh et al. (2009). Further, we use here a mean coupling formulated for fast spiking neurons for our analysis of generic model. Stability analysis of this mean field network allow us to identify all the relevant network states found in the detailed biophysical model. We derive here analytically several boundary solutions, a result which holds for any number of spikes per burst. These findings illustrate the role of oscillations occurring at slow time scales (bursts) on the global behavior of the network.

Intrinsic heterogeneity in oscillatory dynamics limits correlation-induced neural synchronization

Synchronous neural oscillations are found throughout the brain and are thought to contribute to neural coding and the propagation of activity. Several proposed mechanisms of synchronization have gained support through combined theoretical and experimental investigation, including mechanisms based on coupling and correlated input. Here, we ask how correlation-induced synchrony is affected by physiological heterogeneity across neurons. To address this question, we examined cell-to-cell differences in phase-response curves (PRCs), which characterize the response of periodically firing neurons to weak perturbations. Using acute slice electrophysiology, we measured PRCs across a single class of principal neurons capable of sensory-evoked oscillations in vivo: the olfactory bulb mitral cells (MCs). Periodically firing MCs displayed a broad range of PRCs, each of which was well fit by a simple three-parameter model. MCs also displayed differences in firing rate-current relationships and in preferred firing rate ranges. Both the observed PRC heterogeneity and moderate firing rate differences (∼10 Hz) separately reduced the maximum correlation-induced synchrony between MCs by up to 25-30%. Simulations further demonstrated that these components of heterogeneity alone were sufficient to account for the difference in synchronization among heterogeneous vs. homogeneous populations in vitro. Within this simulation framework, independent modulation of specific PRC features additionally revealed which aspects of PRC heterogeneity most strongly impact correlation-induced synchronization. Finally, we demonstrated good agreement of novel mathematical theory with our experimental and simulation results, providing a theoretical basis for the influence of heterogeneity on correlation-induced neural synchronization.

Optimal heterogeneity for coding in spiking neural networks

The effect of cellular heterogeneity on the coding properties of neural populations is studied analytically and numerically. We find that heterogeneity decreases the threshold for synchronization, and its strength is nonlinearly related to the network mean firing rate. In addition, conditions are shown under which heterogeneity optimizes network information transmission for either temporal or rate coding, with high input frequencies leading to different effects for each coding strategy. The results are shown to be robust for more realistic conditions.