Dynamics of spontaneous activity in random networks with multiple neuron subtypes and synaptic noise : Spontaneous activity in networks with synaptic noise

The Roles of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slow-wave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking (RS) cells with frequency adaptation and do not exhibit bursts in current-clamp experiments (in vitro). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings (gsyn) and external current inputs (I). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to the bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold (IL) and low-threshold (IT) ion channels play a role in increasing and decreasing the parameter conditions (gsyn and I) in which bistable dynamics occur, respectively. For high values of IL conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance IT increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

The Role of Potassium and Calcium Currents in the Bistable Firing Transition

Healthy brains display a wide range of firing patterns, from synchronized oscillations during slowwave sleep to desynchronized firing during movement. These physiological activities coexist with periods of pathological hyperactivity in the epileptic brain, where neurons can fire in synchronized bursts. Most cortical neurons are pyramidal regular spiking cells (RS) with frequency adaptation and do not exhibit bursts in current-clamp experiments ( in vitro ). In this work, we investigate the transition mechanism of spike-to-burst patterns due to slow potassium and calcium currents, considering a conductance-based model of a cortical RS cell. The joint influence of potassium and calcium ion channels on high synchronous patterns is investigated for different synaptic couplings ( g syn ) and external current inputs ( I ). Our results suggest that slow potassium currents play an important role in the emergence of high-synchronous activities, as well as in the spike-to-burst firing pattern transitions. This transition is related to bistable dynamics of the neuronal network, where physiological asynchronous states coexist with pathological burst synchronization. The hysteresis curve of the coefficient of variation of the inter-spike interval demonstrates that a burst can be initiated by firing states with neuronal synchronization. Furthermore, we notice that high-threshold ( I L ) and low-threshold ( I T ) ion channels play a role in increasing and decreasing the parameter conditions ( g syn and I ) in which bistable dynamics occur, respectively. For high values of I L conductance, a synchronous burst appears when neurons are weakly coupled and receive more external input. On the other hand, when the conductance I T increases, higher coupling and lower I are necessary to produce burst synchronization. In light of our results, we suggest that channel subtype-specific pharmacological interactions can be useful to induce transitions from pathological high bursting states to healthy states.

Anti-Disturbance of Scale-Free Spiking Neural Network against Impulse Noise

The bio-brain presents robustness function to external stimulus through its self-adaptive regulation and neural information processing. Drawing from the advantages of the bio-brain to investigate the robustness function of a spiking neural network (SNN) is conducive to the advance of brain-like intelligence. However, the current brain-like model is insufficient in biological rationality. In addition, its evaluation method for anti-disturbance performance is inadequate. To explore the self-adaptive regulation performance of a brain-like model with more biological rationality under external noise, a scale-free spiking neural network(SFSNN) is constructed in this study. Then, the anti-disturbance ability of the SFSNN against impulse noise is investigated, and the anti-disturbance mechanism is further discussed. Our simulation results indicate that: (i) our SFSNN has anti-disturbance ability against impulse noise, and the high-clustering SFSNN outperforms the low-clustering SFSNN in terms of anti-disturbance performance. (ii) The neural information processing in the SFSNN under external noise is clarified, which is a dynamic chain effect of the neuron firing, the synaptic weight, and the topological characteristic. (iii) Our discussion hints that an intrinsic factor of the anti-disturbance ability is the synaptic plasticity, and the network topology is a factor that affects the anti-disturbance ability at the level of performance.

Plasticity impairment alters community structure but permits successful pattern separation in a hippocampal network model

Patients who suffer from traumatic brain injury (TBI) often complain of learning and memory problems. Their symptoms are principally mediated by the hippocampus and the ability to adapt to stimulus, also known as neural plasticity. Therefore, one plausible injury mechanism is plasticity impairment, which currently lacks comprehensive investigation across TBI research. For these studies, we used a computational network model of the hippocampus that includes the dentate gyrus, CA3, and CA1 with neuron-scale resolution. We simulated mild injury through weakened spike-timing-dependent plasticity (STDP), which modulates synaptic weights according to causal spike timing. In preliminary work, we found functional deficits consisting of decreased firing rate and broadband power in areas CA3 and CA1 after STDP impairment. To address structural changes with these studies, we applied modularity analysis to evaluate how STDP impairment modifies community structure in the hippocampal network. We also studied the emergent function of network-based learning and found that impaired networks could acquire conditioned responses after training, but the magnitude of the response was significantly lower. Furthermore, we examined pattern separation, a prerequisite of learning, by entraining two overlapping patterns. Contrary to our initial hypothesis, impaired networks did not exhibit deficits in pattern separation with either population- or rate-based coding. Collectively, these results demonstrate how a mechanism of injury that operates at the synapse regulates circuit function.

Oscillations and variability in neuronal systems: interplay of autonomous transient dynamics and fast deterministic fluctuations

Neuronal systems are subject to rapid fluctuations both intrinsically and externally. These fluctuations can be disruptive or constructive. We investigate the dynamic mechanisms underlying the interactions between rapidly fluctuating signals and the intrinsic properties of the target cells to produce variable and/or coherent responses. We use linearized and non-linear conductance-based models and piecewise constant (PWC) inputs with short duration pieces. The amplitude distributions of the constant pieces consist of arbitrary permutations of a baseline PWC function. In each trial within a given protocol we use one of these permutations and each protocol consists of a subset of all possible permutations, which is the only source of uncertainty in the protocol. We show that sustained oscillatory behavior can be generated in response to various forms of PWC inputs independently of whether the stable equilibria of the corresponding unperturbed systems are foci or nodes. The oscillatory voltage responses are amplified by the model nonlinearities and attenuated for conductance-based PWC inputs as compared to current-based PWC inputs, consistent with previous theoretical and experimental work. In addition, the voltage responses to PWC inputs exhibited variability across trials, which is reminiscent of the variability generated by stochastic noise (e.g., Gaussian white noise). Our analysis demonstrates that both oscillations and variability are the result of the interaction between the PWC input and the target cell's autonomous transient dynamics with little to no contribution from the dynamics in vicinities of the steady-state, and do not require input stochasticity.

Heterogeneous Responses to Changes in Inhibitory Synaptic Strength in Networks of Spiking Neurons

How does the dynamics of neurons in a network respond to changes in synaptic weights? Answer to this question would be important for a full understanding of synaptic plasticity. In this article, we report our numerical study of the effects of changes in inhibitory synaptic weights on the spontaneous activity of networks of spiking neurons with conductance-based synapses. Networks with biologically realistic features, which were reconstructed from multi-electrode array recordings taken in a cortical neuronal culture, and their modifications were used in the simulations. The magnitudes of the synaptic weights of all the inhibitory connections are decreased by a uniform amount subjecting to the condition that inhibitory connections would not be turned into excitatory ones. Our simulation results reveal that the responses of the neurons are heterogeneous: while the firing rate of some neurons increases as expected, the firing rate of other neurons decreases or remains unchanged. The same results show that heterogeneous responses also occur for an enhancement of inhibition. This heterogeneity in the responses of neurons to changes in inhibitory synaptic strength suggests that activity-induced modification of synaptic strength does not necessarily generate a positive feedback loop on the dynamics of neurons connected in a network. Our results could be used to understand the effects of bicuculline on spiking and bursting activities of neuronal cultures. Using reconstructed networks with biologically realistic features enables us to identify a long-tailed distribution of average synaptic weights for outgoing links as a crucial feature in giving rise to bursting in neuronal networks and in determining the overall response of the whole network to changes in synaptic strength. For networks whose average synaptic weights for outgoing links have a long-tailed distribution, bursting is observed and the average firing rate of the whole network increases upon inhibition suppression or decreases upon inhibition enhancement. For networks whose average synaptic weights for outgoing links are approximately normally distributed, bursting is not found and the average firing rate of the whole network remains approximately constant upon changes in inhibitory synaptic strength.

Reduction of the collective dynamics of neural populations with realistic forms of heterogeneity

Reduction of collective dynamics of large heterogeneous populations to low-dimensional mean-field models is an important task of modern theoretical neuroscience. Such models can be derived from microscopic equations, for example with the help of Ott-Antonsen theory. An often used assumption of the Lorentzian distribution of the unit parameters makes the reduction especially efficient. However, the Lorentzian distribution is often implausible as having undefined moments, and the collective behavior of populations with other distributions needs to be studied. In the present Letter we propose a method which allows efficient reduction for an arbitrary distribution and show how it performs for the Gaussian distribution. We show that a reduced system for several macroscopic complex variables provides an accurate description of a population of thousands of neurons. Using this reduction technique we demonstrate that the population dynamics depends significantly on the form of its parameter distribution. In particular, the dynamics of populations with Lorentzian and Gaussian distributions with the same center and width differ drastically.

Pharmacological reversal of synaptic and network pathology in human MECP2-KO neurons and cortical organoids

Duplication or deficiency of the X-linked MECP2 gene reliably produces profound neurodevelopmental impairment. MECP2 mutations are almost universally responsible for Rett syndrome (RTT), and particular mutations and cellular mosaicism of MECP2 may underlie the spectrum of RTT symptomatic severity. No clinically approved treatments for RTT are currently available, but human pluripotent stem cell technology offers a platform to identify neuropathology and test candidate therapeutics. Using a strategic series of increasingly complex human stem cell-derived technologies, including human neurons, MECP2-mosaic neurospheres to model RTT female brain mosaicism, and cortical organoids, we identified synaptic dysregulation downstream from knockout of MECP2 and screened select pharmacological compounds for their ability to treat this dysfunction. Two lead compounds, Nefiracetam and PHA 543613, specifically reversed MECP2-knockout cytologic neuropathology. The capacity of these compounds to reverse neuropathologic phenotypes and networks in human models supports clinical studies for neurodevelopmental disorders in which MeCP2 deficiency is the predominant etiology.

Neurodegeneration exposes firing rate dependent effects on oscillation dynamics in computational neural networks

Traumatic brain injury (TBI) can lead to neurodegeneration in the injured circuitry, either through primary structural damage to the neuron or secondary effects that disrupt key cellular processes. Moreover, traumatic injuries can preferentially impact subpopulations of neurons, but the functional network effects of these targeted degeneration profiles remain unclear. Although isolating the consequences of complex injury dynamics and long-term recovery of the circuit can be difficult to control experimentally, computational networks can be a powerful tool to analyze the consequences of injury. Here, we use the Izhikevich spiking neuron model to create networks representative of cortical tissue. After an initial settling period with spike-timing-dependent plasticity (STDP), networks developed rhythmic oscillations similar to those seen in vivo. As neurons were sequentially removed from the network, population activity rate and oscillation dynamics were significantly reduced. In a successive period of network restructuring with STDP, network activity levels returned to baseline for some injury levels and oscillation dynamics significantly improved. We next explored the role that specific neurons have in the creation and termination of oscillation dynamics. We determined that oscillations initiate from activation of low firing rate neurons with limited structural inputs. To terminate oscillations, high activity excitatory neurons with strong input connectivity activate downstream inhibitory circuitry. Finally, we confirm the excitatory neuron population role through targeted neurodegeneration. These results suggest targeted neurodegeneration can play a key role in the oscillation dynamics after injury.

Hot coffee: associative memory with bump attractor cell assemblies of spiking neurons

Networks of spiking neurons can have persistently firing stable bump attractors to represent continuous spaces (like temperature). This can be done with a topology with local excitatory synapses and local surround inhibitory synapses. Activating large ranges in the attractor can lead to multiple bumps, that show repeller and attractor dynamics; however, these bumps can be merged by overcoming the repeller dynamics. A simple associative memory can include these bump attractors, allowing the use of continuous variables in these memories, and these associations can be learned by Hebbian rules. These simulations are related to biological networks, showing that this is a step toward a more complete neural cognitive associative memory.

Modeling Hippocampal CA1 Gabaergic Synapses of Audiogenic Rats

Wistar Audiogenic Rats (WARs) are genetically susceptible to sound-induced seizures that start in the brainstem and, in response to repetitive stimulation, spread to limbic areas, such as hippocampus. Analysis of the distribution of interevent intervals of GABAergic inhibitory postsynaptic currents (IPSCs) in CA1 pyramidal cells showed a monoexponential trend in Wistar rats, suggestive of a homogeneous population of synapses, but a biexponential trend in WARs. Based on this, we hypothesize that there are two populations of GABAergic synaptic release sites in CA1 pyramidal neurons from WARs. To address this hypothesis, we used a well-established neuronal computational model of a CA1 pyramidal neuron previously developed to replicate physiological properties of these cells. Our simulations replicated the biexponential trend only when we decreased the release frequency of synaptic currents by a factor of six in at least 40% of distal synapses. Our results suggest that almost half of the GABAergic synapses of WARs have a drastically reduced spontaneous release frequency. The computational model was able to reproduce the temporal dynamics of GABAergic inhibition that could underlie susceptibility to the spread of seizures.

Optimal Interplay between Synaptic Strengths and Network Structure Enhances Activity Fluctuations and Information Propagation in Hierarchical Modular Networks

In network models of spiking neurons, the joint impact of network structure and synaptic parameters on activity propagation is still an open problem. Here, we use an information-theoretical approach to investigate activity propagation in spiking networks with a hierarchical modular topology. We observe that optimized pairwise information propagation emerges due to the increase of either (i) the global synaptic strength parameter or (ii) the number of modules in the network, while the network size remains constant. At the population level, information propagation of activity among adjacent modules is enhanced as the number of modules increases until a maximum value is reached and then decreases, showing that there is an optimal interplay between synaptic strength and modularity for population information flow. This is in contrast to information propagation evaluated among pairs of neurons, which attains maximum value at the maximum values of these two parameter ranges. By examining the network behavior under the increase of synaptic strength and the number of modules, we find that these increases are associated with two different effects: (i) the increase of autocorrelations among individual neurons and (ii) the increase of cross-correlations among pairs of neurons. The second effect is associated with better information propagation in the network. Our results suggest roles that link topological features and synaptic strength levels to the transmission of information in cortical networks.

Oscillations and collective excitability in a model of stochastic neurons under excitatory and inhibitory coupling

We study a model with excitable neurons modeled as stochastic units with three states, representing quiescence, firing, and refractoriness. The transition rates between quiescence and firing depend exponentially on the number of firing neighbors, whereas all other rates are kept constant. This model class was shown to exhibit collective oscillations (synchronization) if neurons are spiking autonomously, but not if neurons are in the excitable regime. In both cases, neurons were restricted to interact through excitatory coupling. Here we show that a plethora of collective phenomena appear if inhibitory coupling is added. Besides the usual transition between an absorbing and an active phase, the model with excitatory and inhibitory neurons can also undergo reentrant transitions to an oscillatory phase. In the mean-field description, oscillations can emerge through supercritical or subcritical Hopf bifurcations, as well as through infinite period bifurcations. The model has bistability between active and oscillating behavior, as well as collective excitability, a regime where the system can display a peak of global activity when subject to a sufficiently strong perturbation. We employ a variant of the Shinomoto-Kuramoto order parameter to characterize the phase transitions and their system-size dependence.