Chapter 21 Mechanisms of synchrony of neural activity in large networks

Electroacupuncture stimulation to modulate neural oscillations in promoting neurological rehabilitation

Electroacupuncture (EA) stimulation is a modern neuromodulation technique that integrates traditional Chinese acupuncture therapy with contemporary electrical stimulation. It involves the application of electrical currents to specific acupoints on the body following acupuncture. EA has been widely used in the treatment of various neurological disorders, including epilepsy, stroke, Parkinson's disease, and Alzheimer's disease. Recent research suggests that EA stimulation may modulate neural oscillations, correcting abnormal brain electrical activity, therefore promoting brain function and aiding in neurological rehabilitation. This paper conducted a comprehensive search in databases such as PubMed, Web of Science, and CNKI using keywords like "electroacupuncture," "neural oscillations," and "neurorehabilitation", covering the period from year 1980 to 2023. We provide a detailed overview of how electroacupuncture stimulation modulates neural oscillations, including maintaining neural activity homeostasis, influencing neurotransmitter release, improving cerebral hemodynamics, and enhancing specific neural functional networks. The paper also discusses the current state of research, limitations of electroacupuncture-induced neural oscillation techniques, and explores prospects for their combined application, aiming to offer broader insights for both basic and clinical research.

Transition from inhomogeneous limit cycles to oscillation death in nonlinear oscillators with similarity-dependent coupling

We consider a system of coupled nonlinear oscillators in which the interaction is modulated by a measure of the similarity between the oscillators. Such a coupling is common in treating spatially mobile dynamical systems where the interaction is distance dependent or in resonance-enhanced interactions, for instance. For a system of Stuart-Landau oscillators coupled in this manner, we observe a novel route to oscillation death via a Hopf bifurcation. The individual oscillators are confined to inhomogeneous limit cycles initially and are damped to different fixed points after the bifurcation. Analytical and numerical results are presented for this case, while numerical results are presented for coupled Rössler and Sprott oscillators.

Multistability and evolution of chimera states in a network of type II Morris-Lecar neurons with asymmetrical nonlocal inhibitory connections

Formation of synchronous activity patterns is an essential property of neuronal networks that has been of central interest to synchronization theory. Chimera states, where both synchronous and asynchronous activities of neurons co-exist in a single network, are particularly poignant examples of such patterns, whose dynamics and multistability may underlie brain function, such as cognitive tasks. However, dynamical mechanisms of coherent state formation in spiking neuronal networks as well as ways to control these states remain unclear. In this paper, we take a step in this direction by considering the evolution of chimera states in a network of class II excitable Morris-Lecar neurons with asymmetrical nonlocal inhibitory connections. Using the adaptive coherence measure, we are able to partition the network parameter space into regions of various collective behaviors (antiphase synchronous clusters, traveling waves, different types of chimera states as well as a spiking death regime) and have shown multistability between the various regimes. We track the evolution of the chimera states as a function of changed key network parameters and found transitions between various types of chimera states. We further find that the network can demonstrate long transients leading to quasi-persistence of activity patterns in the border regions hinting at near-criticality behaviors.

Collective dynamics in the presence of finite-width pulses

The idealization of neuronal pulses as δ-spikes is a convenient approach in neuroscience but can sometimes lead to erroneous conclusions. We investigate the effect of a finite pulse width on the dynamics of balanced neuronal networks. In particular, we study two populations of identical excitatory and inhibitory neurons in a random network of phase oscillators coupled through exponential pulses with different widths. We consider three coupling functions inspired by leaky integrate-and-fire neurons with delay and type I phase-response curves. By exploring the role of the pulse widths for different coupling strengths, we find a robust collective irregular dynamics, which collapses onto a fully synchronous regime if the inhibitory pulses are sufficiently wider than the excitatory ones. The transition to synchrony is accompanied by hysteretic phenomena (i.e., the co-existence of collective irregular and synchronous dynamics). Our numerical results are supported by a detailed scaling and stability analysis of the fully synchronous solution. A conjectured first-order phase transition emerging for δ-spikes is smoothed out for finite-width pulses.

Does layer 4 in the barrel cortex function as a balanced circuit when responding to whisker movements?

Neurons in one barrel in layer 4 (L4) in the mouse vibrissa somatosensory cortex are innervated mostly by neurons from the VPM nucleus and by other neurons within the same barrel. During quiet wakefulness or whisking in air, thalamic inputs vary slowly in time, and excitatory neurons rarely fire. A barrel in L4 contains a modest amount of neurons; the synaptic conductances are not very strong and connections are not sparse. Are the dynamical properties of the L4 circuit similar to those expected from fluctuation-dominated, balanced networks observed for large, strongly coupled and sparse cortical circuits? To resolve this question, we analyze a network of 150 inhibitory parvalbumin-expressing fast-spiking inhibitory interneurons innervated by the VPM thalamus with random connectivity, without or with 1600 low-firing excitatory neurons. Above threshold, the population-average firing rate of inhibitory cortical neurons increases linearly with the thalamic firing rate. The coefficient of variation CV is somewhat less than 1. Moderate levels of synchrony are induced by converging VPM inputs and by inhibitory interaction among neurons. The strengths of excitatory and inhibitory currents during whisking are about three times larger than threshold. We identify values of numbers of presynaptic neurons, synaptic delays between inhibitory neurons, and electrical coupling within the experimentally plausible ranges for which spike synchrony levels are low. Heterogeneity in in-degrees increases the width of the firing rate distribution to the experimentally observed value. We conclude that an L4 circuit in the low-synchrony regime exhibits qualitative dynamical properties similar to those of balanced networks.

Complex dynamics of an oscillator ensemble with uniformly distributed natural frequencies and global nonlinear coupling

We consider large populations of phase oscillators with global nonlinear coupling. For identical oscillators such populations are known to demonstrate a transition from completely synchronized state to the state of self-organized quasiperiodicity. In this state phases of all units differ, yet the population is not completely incoherent but produces a nonzero mean field; the frequency of the latter differs from the frequency of individual units. Here we analyze the dynamics of such populations in case of uniformly distributed natural frequencies. We demonstrate numerically and describe theoretically (i) states of complete synchrony, (ii) regimes with coexistence of a synchronous cluster and a drifting subpopulation, and (iii) self-organized quasiperiodic states with nonzero mean field and all oscillators drifting with respect to it. We analyze transitions between different states with the increase of the coupling strength; in particular we show that the mean field arises via a discontinuous transition. For a further illustration we compare the results for the nonlinear model with those for the Kuramoto-Sakaguchi model.

Persistent Synchronized Bursting Activity in Cortical Tissues With Low Magnesium Concentration: A Modeling Study

We explore the mechanism of synchronized bursting activity with frequency of ∼10 Hz that appears in cortical tissues at low extracellular magnesium concentration [Mg2+]o. We hypothesize that this activity is persistent, namely coexists with the quiescent state and depends on slow N-methyl-d-aspartate (NMDA) conductances. To explore this hypothesis, we construct and investigate a conductance-based model of excitatory cortical networks. Population bursting activity can persist for physiological values of the NMDA decay time constant (∼100 ms). Neurons are synchronized at the time scale of bursts but not of single spikes. A reduced model of a cell coupled to itself can encompass most of this highly synchronized network behavior and is analyzed using the fast-slow method. Synchronized bursts appear for intermediate values of the NMDA conductance gNMDA if NMDA conductances are not too fast. Regular spiking activity appears for larger gNMDA. If the single cell is a conditional burster, persistent synchronized bursts become more robust. Weakly synchronized states appear for zero AMPA conductance gAMPA. Enhancing gAMPA increases both synchrony and the number of spikes within bursts and decreases the bursting frequency. Too strong gAMPA, however, prevents the activity because it enhances neuronal intrinsic adaptation. When [Mg2+]o is increased, higher gNMDA values are needed to maintain bursting activity. Bursting frequency decreases with [Mg2+]o, and the network is silent with physiological [Mg2+]o. Inhibition weakly decreases the bursting frequency if inhibitory cells receive enough NMDA-mediated excitation. This study explains the importance of conditional bursters in layer V in supporting epileptiform activity at low [Mg2+]o.

Nonlinear multivariate analysis of neurophysiological signals

Multivariate time series analysis is extensively used in neurophysiology with the aim of studying the relationship between simultaneously recorded signals. Recently, advances on information theory and nonlinear dynamical systems theory have allowed the study of various types of synchronization from time series. In this work, we first describe the multivariate linear methods most commonly used in neurophysiology and show that they can be extended to assess the existence of nonlinear interdependence between signals. We then review the concepts of entropy and mutual information followed by a detailed description of nonlinear methods based on the concepts of phase synchronization, generalized synchronization and event synchronization. In all cases, we show how to apply these methods to study different kinds of neurophysiological data. Finally, we illustrate the use of multivariate surrogate data test for the assessment of the strength (strong or weak) and the type (linear or nonlinear) of interdependence between neurophysiological signals.

Analytical solution of the cable equation with synaptic reversal potentials for passive neurons with tip-to-tip dendrodendritic coupling

A passive cable model is presented for a pair of electrotonically coupled neurons in order to investigate the effects of tip-to-tip dendrodendritic gap junctions on the interaction between excitation and either pre or postsynaptic inhibition. The model represents each dendritic tree by a tapered equivalent cylinder attached to an isopotential soma. Analytical solution of the cable equation with synaptic reversal potentials is considered for each neuron to yield a system of Volterra integral equations for the voltage. The solution to the system of linear integral equations (expressed as a Neumann series) is used to determine the current spread within the two coupled neurons, and to re-examine the sensitivity of the soma potentials (in particular) to the coupling resistance for various loci of synaptic inputs. The model is actually posed generally, so that active as well as passive properties could be considered. In the active case, a system of non-linear integral equations is derived for the voltage.