Interplay between excitability type and distributions of neuronal connectivity determines neuronal network synchronization

Modeling cortical synaptic effects of anesthesia and their cholinergic reversal

General anesthetics work through a variety of molecular mechanisms while resulting in the common end point of sedation and loss of consciousness. Generally, the administration of common anesthetics induces reduction in synaptic excitation while promoting synaptic inhibition. Exogenous modulation of the anesthetics' synaptic effects can help determine the neuronal pathways involved in anesthesia. For example, both animal and human studies have shown that exogenously induced increases in acetylcholine in the brain can elicit wakeful-like behavior despite the continued presence of the anesthetic. However, the underlying mechanisms of anesthesia reversal at the cellular level have not been investigated. Here we apply a computational model of a network of excitatory and inhibitory neurons to simulate the network-wide effects of anesthesia, due to changes in synaptic inhibition and excitation, and their reversal by cholinergic activation through muscarinic receptors. We use a differential evolution algorithm to fit model parameters to match measures of spiking activity, neuronal connectivity, and network dynamics recorded in the visual cortex of rodents during anesthesia with desflurane in vivo. We find that facilitating muscarinic receptor effects of acetylcholine on top of anesthetic-induced synaptic changes predicts the reversal of anesthetic suppression of neurons' spiking activity, functional connectivity, as well as pairwise and population interactions. Thus, our model predicts a specific neuronal mechanism for the cholinergic reversal of anesthesia consistent with experimental behavioral observations.

Electrocorticography reveals thalamic control of cortical dynamics following traumatic brain injury

The return of consciousness after traumatic brain injury (TBI) is associated with restoring complex cortical dynamics; however, it is unclear what interactions govern these complex dynamics. Here, we set out to uncover the mechanism underlying the return of consciousness by measuring local field potentials (LFP) using invasive electrophysiological recordings in patients recovering from TBI. We found that injury to the thalamus, and its efferent projections, on MRI were associated with repetitive and low complexity LFP signals from a highly structured phase space, resembling a low-dimensional ring attractor. But why do thalamic injuries in TBI patients result in a cortical attractor? We built a simplified thalamocortical model, which connotes that thalamic input facilitates the formation of cortical ensembles required for the return of cognitive function and the content of consciousness. These observations collectively support the view that thalamic input to the cortex enables rich cortical dynamics associated with consciousness.

Dynamics of a neuronal pacemaker in the weakly electric fish Apteronotus

The precise timing of neuronal activity is critical for normal brain function. In weakly electric fish, the medullary pacemaker network (PN) sets the timing for an oscillating electric organ discharge (EOD) used for electric sensing. This network is the most precise biological oscillator known, with sub-microsecond variation in oscillator period. The PN consists of two principle sets of neurons, pacemaker and relay cells, that are connected by gap junctions and normally fire in synchrony, one-to-one with each EOD cycle. However, the degree of gap junctional connectivity between these cells appears insufficient to provide the population averaging required for the observed temporal precision of the EOD. This has led to the hypothesis that individual cells themselves fire with high precision, but little is known about the oscillatory dynamics of these pacemaker cells. As a first step towards testing this hypothesis, we have developed a biophysical model of a pacemaker neuron action potential based on experimental recordings. We validated the model by comparing the changes in oscillatory dynamics produced by different experimental manipulations. Our results suggest that this relatively simple model can capture a large range of channel dynamics exhibited by pacemaker cells, and will thus provide a basis for future work on network synchrony and precision.

Effects of network topologies on stochastic resonance in feedforward neural network

The effects of network topologies on signal propagation are studied in noisy feedforward neural network in detail, where the network topologies are modulated by changing both the in-degree and out-degree distributions of FFNs as identical, uniform and exponential respectively. Stochastic resonance appeared in three FFNs when the same external stimuli and noise are applied to the three different network topologies. It is found that optimal noise intensity decreases with the increase of network's layer index. Meanwhile, the Q index of FFN with identical distribution is higher than that of the other two FFNs, indicating that the synchronization between the neuronal firing activities and the external stimuli is more obvious in FFN with identical distribution. The optimal parameter regions for the time cycle of external stimuli and the noise intensity are found for three FFNs, in which the resonance is more easily induced when the parameters of stimuli are set in this region. Furthermore, the relationship between the in-degree, the average membrane potential and the resonance performance is studied at the neuronal level, where it is found that both the average membrane potentials and the Q indexes of neurons in FFN with identical degree distribution is more consistent with each other than that of the other two FFNs due to their network topologies. In summary, the simulations here indicate that the network topologies play essential roles in affecting the signal propagation of FFNs.

Emergence of global synchronization in directed excitatory networks of type I neurons

The collective behaviour of neural networks depends on the cellular and synaptic properties of the neurons. The phase-response curve (PRC) is an experimentally obtainable measure of cellular properties that quantifies the shift in the next spike time of a neuron as a function of the phase at which stimulus is delivered to that neuron. The neuronal PRCs can be classified as having either purely positive values (type I) or distinct positive and negative regions (type II). Networks of type 1 PRCs tend not to synchronize via mutual excitatory synaptic connections. We study the synchronization properties of identical type I and type II neurons, assuming unidirectional synapses. Performing the linear stability analysis and the numerical simulation of the extended Kuramoto model, we show that feedforward loop motifs favour synchronization of type I excitatory and inhibitory neurons, while feedback loop motifs destroy their synchronization tendency. Moreover, large directed networks, either without feedback motifs or with many of them, have been constructed from the same undirected backbones, and a high synchronization level is observed for directed acyclic graphs with type I neurons. It has been shown that, the synchronizability of type I neurons depends on both the directionality of the network connectivity and the topology of its undirected backbone. The abundance of feedforward motifs enhances the synchronizability of the directed acyclic graphs.

High-order couplings in geometric complex networks of neurons

We explore the consequences of introducing higher-order interactions in a geometric complex network of Morris-Lecar neurons. We focus on the regime where traveling synchronization waves are observed from a first-neighbors-based coupling to evaluate the changes induced when higher-order dynamical interactions are included. We observe that the traveling-wave phenomenon gets enhanced by these interactions, allowing the activity to travel further in the system without generating pathological full synchronization states. This scheme could be a step toward a simple phenomenological modelization of neuroglial networks.

Acetylcholine Mediates Dynamic Switching Between Information Coding Schemes in Neuronal Networks

Rate coding and phase coding are the two major coding modes seen in the brain. For these two modes, network dynamics must either have a wide distribution of frequencies for rate coding, or a narrow one to achieve stability in phase dynamics for phase coding. Acetylcholine (ACh) is a potent regulator of neural excitability. Acting through the muscarinic receptor, ACh reduces the magnitude of the potassium M-current, a hyperpolarizing current that builds up as neurons fire. The M-current contributes to several excitability features of neurons, becoming a major player in facilitating the transition between Type 1 (integrator) and Type 2 (resonator) excitability. In this paper we argue that this transition enables a dynamic switch between rate coding and phase coding as levels of ACh release change. When a network is in a high ACh state variations in synaptic inputs will lead to a wider distribution of firing rates across the network and this distribution will reflect the network structure or pattern of external input to the network. When ACh is low, network frequencies become narrowly distributed and the structure of a network or pattern of external inputs will be represented through phase relationships between firing neurons. This work provides insights into how modulation of neuronal features influences network dynamics and information processing across brain states.

Effects of topologies on signal propagation in feedforward networks

We systematically investigate the effects of topologies on signal propagation in feedforward networks (FFNs) based on the FitzHugh-Nagumo neuron model. FFNs with different topological structures are constructed with same number of both in-degrees and out-degrees in each layer and given the same input signal. The propagation of firing patterns and firing rates are found to be affected by the distribution of neuron connections in the FFNs. Synchronous firing patterns emerge in the later layers of FFNs with identical, uniform, and exponential degree distributions, but the number of synchronous spike trains in the output layers of the three topologies obviously differs from one another. The firing rates in the output layers of the three FFNs can be ordered from high to low according to their topological structures as exponential, uniform, and identical distributions, respectively. Interestingly, the sequence of spiking regularity in the output layers of the three FFNs is consistent with the firing rates, but their firing synchronization is in the opposite order. In summary, the node degree is an important factor that can dramatically influence the neuronal network activity.

Development of structural correlations and synchronization from adaptive rewiring in networks of Kuramoto oscillators

Synchronization of non-identical oscillators coupled through complex networks is an important example of collective behavior, and it is interesting to ask how the structural organization of network interactions influences this process. Several studies have explored and uncovered optimal topologies for synchronization by making purposeful alterations to a network. On the other hand, the connectivity patterns of many natural systems are often not static, but are rather modulated over time according to their dynamics. However, this co-evolution and the extent to which the dynamics of the individual units can shape the organization of the network itself are less well understood. Here, we study initially randomly connected but locally adaptive networks of Kuramoto oscillators. In particular, the system employs a co-evolutionary rewiring strategy that depends only on the instantaneous, pairwise phase differences of neighboring oscillators, and that conserves the total number of edges, allowing the effects of local reorganization to be isolated. We find that a simple rule-which preserves connections between more out-of-phase oscillators while rewiring connections between more in-phase oscillators-can cause initially disordered networks to organize into more structured topologies that support enhanced synchronization dynamics. We examine how this process unfolds over time, finding a dependence on the intrinsic frequencies of the oscillators, the global coupling, and the network density, in terms of how the adaptive mechanism reorganizes the network and influences the dynamics. Importantly, for large enough coupling and after sufficient adaptation, the resulting networks exhibit interesting characteristics, including degree-frequency and frequency-neighbor frequency correlations. These properties have previously been associated with optimal synchronization or explosive transitions in which the networks were constructed using global information. On the contrary, by considering a time-dependent interplay between structure and dynamics, this work offers a mechanism through which emergent phenomena and organization can arise in complex systems utilizing local rules.