Multistability in the Kuramoto model with synaptic plasticity

Adaptive coupling between neurons widens the entrainment range of the suprachiasmatic nucleus

In many realistic systems, such as neural networks in the brain, the coupling strength between neurons is not fixed, but adaptively adjusts according to their activities. The suprachiasmatic nucleus (SCN), as the main clock in the mammalian brain, has been found to be a plastic neural network, and the coupling strength between neurons is highly dynamical. An important function of the SCN is entrainment, reflecting the ability of the SCN to synchronize with the external light-dark cycle. The entrainment ability is reflected by the entrainment range, which is a period range for the external light-dark cycle to which the SCN can entrain. In this article, we investigated whether the entrainment range of the SCN is affected by the adaptive coupling. We use a modified Kuramoto model with external light-dark cycle. We found that when the light sensitivity is larger than the fixed coupling strength (the coupling strength without adaptive rules), adaptive coupling can widen the entrainment range. Our findings help to understand the impact of the adaptive coupling between oscillatorty neurons on the collective behavior of the SCN, and provides a possible explanation for the plasticity of coupling in the master clock network.

Stimulus-induced dynamical states in an adaptive network with symmetric adaptation

We consider an adaptive network of identical phase oscillators with the symmetric adaptation rule for the evolution of the connection weights under the influence of an external force. We show that the adaptive network exhibits a plethora of self-organizing dynamical states such as the two-cluster state, multiantipodal clusters, splay cluster, splay chimera, forced entrained state, chimera state, bump state, coherent, and incoherent states in the two-parameter phase diagrams. The intriguing structures of the frequency clusters and instantaneous phases of the oscillators characterize the distinct self-organized synchronized and partial synchronized states. The hierarchical organization of the frequency clusters, resulting in strongly coupled subnetworks, is also evident from the dynamics of the coupling weights, where the frequency clusters are either very weakly coupled or even completely decoupled from each other. Additionally, we also deduce the stability condition for the forced entrained state.

Synaptic reorganization of synchronized neuronal networks with synaptic weight and structural plasticity

Abnormally strong neural synchronization may impair brain function, as observed in several brain disorders. We computationally study how neuronal dynamics, synaptic weights, and network structure co-emerge, in particular, during (de)synchronization processes and how they are affected by external perturbation. To investigate the impact of different types of plasticity mechanisms, we combine a network of excitatory integrate-and-fire neurons with different synaptic weight and/or structural plasticity mechanisms: (i) only spike-timing-dependent plasticity (STDP), (ii) only homeostatic structural plasticity (hSP), i.e., without weight-dependent pruning and without STDP, (iii) a combination of STDP and hSP, i.e., without weight-dependent pruning, and (iv) a combination of STDP and structural plasticity (SP) that includes hSP and weight-dependent pruning. To accommodate the diverse time scales of neuronal firing, STDP, and SP, we introduce a simple stochastic SP model, enabling detailed numerical analyses. With tools from network theory, we reveal that structural reorganization may remarkably enhance the network's level of synchrony. When weaker contacts are preferentially eliminated by weight-dependent pruning, synchrony is achieved with significantly sparser connections than in randomly structured networks in the STDP-only model. In particular, the strengthening of contacts from neurons with higher natural firing rates to those with lower rates and the weakening of contacts in the opposite direction, followed by selective removal of weak contacts, allows for strong synchrony with fewer connections. This activity-led network reorganization results in the emergence of degree-frequency, degree-degree correlations, and a mixture of degree assortativity. We compare the stimulation-induced desynchronization of synchronized states in the STDP-only model (i) with the desynchronization of models (iii) and (iv). The latter require stimuli of significantly higher intensity to achieve long-term desynchronization. These findings may inform future pre-clinical and clinical studies with invasive or non-invasive stimulus modalities aiming at inducing long-lasting relief of symptoms, e.g., in Parkinson's disease.

Coordinated reset stimulation of plastic neural networks with spatially dependent synaptic connections

Background

Abnormal neuronal synchrony is associated with several neurological disorders, including Parkinson's disease (PD), essential tremor, dystonia, and epilepsy. Coordinated reset (CR) stimulation was developed computationally to counteract abnormal neuronal synchrony. During CR stimulation, phase-shifted stimuli are delivered to multiple stimulation sites. Computational studies in plastic neural networks reported that CR stimulation drove the networks into an attractor of a stable desynchronized state by down-regulating synaptic connections, which led to long-lasting desynchronization effects that outlasted stimulation. Later, corresponding long-lasting desynchronization and therapeutic effects were found in animal models of PD and PD patients. To date, it is unclear how spatially dependent synaptic connections, as typically observed in the brain, shape CR-induced synaptic downregulation and long-lasting effects.

Methods

We performed numerical simulations of networks of leaky integrate-and-fire neurons with spike-timing-dependent plasticity and spatially dependent synaptic connections to study and further improve acute and long-term responses to CR stimulation.

Results

The characteristic length scale of synaptic connections relative to the distance between stimulation sites plays a key role in CR parameter adjustment. In networks with short synaptic length scales, a substantial synaptic downregulation can be achieved by selecting appropriate stimulus-related parameters, such as the stimulus amplitude and shape, regardless of the employed spatiotemporal pattern of stimulus deliveries. Complex stimulus shapes can induce local connectivity patterns in the vicinity of the stimulation sites. In contrast, in networks with longer synaptic length scales, the spatiotemporal sequence of stimulus deliveries is of major importance for synaptic downregulation. In particular, rapid shuffling of the stimulus sequence is advantageous for synaptic downregulation.

Conclusion

Our results suggest that CR stimulation parameters can be adjusted to synaptic connectivity to further improve the long-lasting effects. Furthermore, shuffling of CR sequences is advantageous for long-lasting desynchronization effects. Our work provides important hypotheses on CR parameter selection for future preclinical and clinical studies.

The influence of synaptic plasticity on critical coupling estimates for neural populations

The presence or absence of synaptic plasticity can dramatically influence the collective behavior of populations of coupled neurons. In this work, we consider spike-timing dependent plasticity (STDP) and its resulting influence on phase cohesion in computational models of heterogeneous populations of conductance-based neurons. STDP allows for the influence of individual synapses to change over time, strengthening or weakening depending on the relative timing of the relevant action potentials. Using phase reduction techniques, we derive an upper bound on the critical coupling strength required to retain phase cohesion for a network of synaptically coupled, heterogeneous neurons with STDP. We find that including STDP can significantly alter phase cohesion as compared to a network with static synaptic connections. Analytical results are validated numerically. Our analysis highlights the importance of the relative ordering of action potentials emitted in a population of tonically firing neurons and demonstrates that order switching can degrade the synchronizing influence of coupling when STDP is considered.

Neuromodulatory effects on synchrony and network reorganization in networks of coupled Kuramoto oscillators

Neuromodulatory processes in the brain can critically change signal processing on a cellular level leading to dramatic changes in network level reorganization. Here, we use coupled non-identical Kuramoto oscillators to investigate how changes in the shape of phase response curves from Type 1 to Type 2, mediated by varying ACh levels, coupled with activity dependent plasticity may alter network reorganization. We first show that when plasticity is absent, the Type 1 networks, as expected, exhibit asynchronous dynamics with oscillators of the highest natural frequency robustly evolving faster in terms of their phase dynamics. At the same time, the Type 2 networks synchronize, with oscillators locked so that the ones with higher natural frequency have a constant phase lead as compared to the ones with lower natural frequency. This relationship establishes a robust mapping between the frequency and oscillators' phases in the network, leading to structure/frequency mapping when plasticity is present. Further we show that while connection plasticity can produce stable synchrony (so called splay states) in Type 1 networks, the structure/frequency reorganization observed in Type 2 networks is not present.

Mean-Field Approximations With Adaptive Coupling for Networks With Spike-Timing-Dependent Plasticity

Understanding the effect of spike-timing-dependent plasticity (STDP) is key to elucidating how neural networks change over long timescales and to design interventions aimed at modulating such networks in neurological disorders. However, progress is restricted by the significant computational cost associated with simulating neural network models with STDP and by the lack of low-dimensional description that could provide analytical insights. Phase-difference-dependent plasticity (PDDP) rules approximate STDP in phase oscillator networks, which prescribe synaptic changes based on phase differences of neuron pairs rather than differences in spike timing. Here we construct mean-field approximations for phase oscillator networks with STDP to describe part of the phase space for this very high-dimensional system. We first show that single-harmonic PDDP rules can approximate a simple form of symmetric STDP, while multiharmonic rules are required to accurately approximate causal STDP. We then derive exact expressions for the evolution of the average PDDP coupling weight in terms of network synchrony. For adaptive networks of Kuramoto oscillators that form clusters, we formulate a family of low-dimensional descriptions based on the mean-field dynamics of each cluster and average coupling weights between and within clusters. Finally, we show that such a two-cluster mean-field model can be fitted to synthetic data to provide a low-dimensional approximation of a full adaptive network with symmetric STDP. Our framework represents a step toward a low-dimensional description of adaptive networks with STDP, and could for example inform the development of new therapies aimed at maximizing the long-lasting effects of brain stimulation.

Synchronization transitions in Kuramoto networks with higher-mode interaction

Synchronization is an omnipresent collective phenomenon in nature and technology, whose understanding is still elusive for real-world systems in particular. We study the synchronization transition in a phase oscillator system with two nonvanishing Fourier-modes in the interaction function, hence going beyond the Kuramoto paradigm. We show that the transition scenarios crucially depend on the interplay of the two coupling modes. We describe the multistability induced by the presence of a second coupling mode. By extending the collective coordinate approach, we describe the emergence of various states observed in the transition from incoherence to coherence. Remarkably, our analysis suggests that, in essence, the two-mode coupling gives rise to states characterized by two independent but interacting groups of oscillators. We believe that these findings will stimulate future research on dynamical systems, including complex interaction functions beyond the Kuramoto-type.

Perspectives on adaptive dynamical systems

Adaptivity is a dynamical feature that is omnipresent in nature, socio-economics, and technology. For example, adaptive couplings appear in various real-world systems, such as the power grid, social, and neural networks, and they form the backbone of closed-loop control strategies and machine learning algorithms. In this article, we provide an interdisciplinary perspective on adaptive systems. We reflect on the notion and terminology of adaptivity in different disciplines and discuss which role adaptivity plays for various fields. We highlight common open challenges and give perspectives on future research directions, looking to inspire interdisciplinary approaches.

Transient Phase Clusters in a Two-Population Network of Kuramoto Oscillators with Heterogeneous Adaptive Interaction

Adaptive interactions are an important property of many real-word network systems. A feature of such networks is the change in their connectivity depending on the current states of the interacting elements. In this work, we study the question of how the heterogeneous character of adaptive couplings influences the emergence of new scenarios in the collective behavior of networks. Within the framework of a two-population network of coupled phase oscillators, we analyze the role of various factors of heterogeneous interaction, such as the rules of coupling adaptation and the rate of their change in the formation of various types of coherent behavior of the network. We show that various schemes of heterogeneous adaptation lead to the formation of transient phase clusters of various types.

Complex dynamics in adaptive phase oscillator networks

Networks of coupled dynamical units give rise to collective dynamics such as the synchronization of oscillators or neurons in the brain. The ability of the network to adapt coupling strengths between units in accordance with their activity arises naturally in a variety of contexts, including neural plasticity in the brain, and adds an additional layer of complexity: the dynamics on the nodes influence the dynamics of the network and vice versa. We study a minimal model of Kuramoto phase oscillators including a general adaptive learning rule with three parameters (strength of adaptivity, adaptivity offset, adaptivity shift), mimicking learning paradigms based on spike-time-dependent plasticity. Importantly, the strength of adaptivity allows to tune the system away from the limit of the classical Kuramoto model, corresponding to stationary coupling strengths and no adaptation and, thus, to systematically study the impact of adaptivity on the collective dynamics. We carry out a detailed bifurcation analysis for the minimal model consisting of N=2 oscillators. The non-adaptive Kuramoto model exhibits very simple dynamic behavior, drift, or frequency-locking; but once the strength of adaptivity exceeds a critical threshold non-trivial bifurcation structures unravel: A symmetric adaptation rule results in multi-stability and bifurcation scenarios, and an asymmetric adaptation rule generates even more intriguing and rich dynamics, including a period-doubling cascade to chaos as well as oscillations displaying features of both librations and rotations simultaneously. Generally, adaptation improves the synchronizability of the oscillators. Finally, we also numerically investigate a larger system consisting of N=50 oscillators and compare the resulting dynamics with the case of N=2 oscillators.

Inhibitory neurons control the consolidation of neural assemblies via adaptation to selective stimuli

Brain circuits display modular architecture at different scales of organization. Such neural assemblies are typically associated to functional specialization but the mechanisms leading to their emergence and consolidation still remain elusive. In this paper we investigate the role of inhibition in structuring new neural assemblies driven by the entrainment to various inputs. In particular, we focus on the role of partially synchronized dynamics for the creation and maintenance of structural modules in neural circuits by considering a network of excitatory and inhibitory [Formula: see text]-neurons with plastic Hebbian synapses. The learning process consists of an entrainment to temporally alternating stimuli that are applied to separate regions of the network. This entrainment leads to the emergence of modular structures. Contrary to common practice in artificial neural networks-where the acquired weights are typically frozen after the learning session-we allow for synaptic adaptation even after the learning phase. We find that the presence of inhibitory neurons in the network is crucial for the emergence and the post-learning consolidation of the modular structures. Indeed networks made of purely excitatory neurons or of neurons not respecting Dale's principle are unable to form or to maintain the modular architecture induced by the stimuli. We also demonstrate that the number of inhibitory neurons in the network is directly related to the maximal number of neural assemblies that can be consolidated, supporting the idea that inhibition has a direct impact on the memory capacity of the neural network.

Heterogeneous Nucleation in Finite-Size Adaptive Dynamical Networks

Phase transitions in equilibrium and nonequilibrium systems play a major role in the natural sciences. In dynamical networks, phase transitions organize qualitative changes in the collective behavior of coupled dynamical units. Adaptive dynamical networks feature a connectivity structure that changes over time, coevolving with the nodes' dynamical state. In this Letter, we show the emergence of two distinct first-order nonequilibrium phase transitions in a finite-size adaptive network of heterogeneous phase oscillators. Depending on the nature of defects in the internal frequency distribution, we observe either an abrupt single-step transition to full synchronization or a more gradual multistep transition. This observation has a striking resemblance to heterogeneous nucleation. We develop a mean-field approach to study the interplay between adaptivity and nodal heterogeneity and describe the dynamics of multicluster states and their role in determining the character of the phase transition. Our work provides a theoretical framework for studying the interplay between adaptivity and nodal heterogeneity.

Interplay of different synchronization modes and synaptic plasticity in a system of class I neurons

We analyze the effect of spike-timing-dependent plasticity (STDP) on a system of pulse-coupled class I neurons. Our research begins with a system of two mutually connected quadratic integrate-and-fire (QIF) neurons, which are canonical representatives of class I neurons. Along with various asymptotic modes previously observed in other neuronal models with plastic synapses, we found a stable synchronous mode characterized by unidirectional link from a slower neuron to a faster neuron. In this frequency-locked mode, the faster neuron emits multiple spikes per cycle of the slower neuron. We analytically obtain the Arnold tongues for this mode without STDP and with STDP. We also consider larger plastic networks of QIF neurons and show that the detected mode can manifest itself in such a way that slow neurons become pacemakers. As a result, slow and fast neurons can form large synchronous clusters that generate low-frequency oscillations. We demonstrate the generality of the results obtained with two connected QIF neurons using Wang-Buzsáki and Morris-Lecar biophysically plausible class I neuron models.

Synaptic reshaping of plastic neuronal networks by periodic multichannel stimulation with single-pulse and burst stimuli

Synaptic dysfunction is associated with several brain disorders, including Alzheimer's disease, Parkinson's disease (PD) and obsessive compulsive disorder (OCD). Utilizing synaptic plasticity, brain stimulation is capable of reshaping synaptic connectivity. This may pave the way for novel therapies that specifically counteract pathological synaptic connectivity. For instance, in PD, novel multichannel coordinated reset stimulation (CRS) was designed to counteract neuronal synchrony and down-regulate pathological synaptic connectivity. CRS was shown to entail long-lasting therapeutic aftereffects in PD patients and related animal models. This is in marked contrast to conventional deep brain stimulation (DBS) therapy, where PD symptoms return shortly after stimulation ceases. In the present paper, we study synaptic reshaping by periodic multichannel stimulation (PMCS) in networks of leaky integrate-and-fire (LIF) neurons with spike-timing-dependent plasticity (STDP). During PMCS, phase-shifted periodic stimulus trains are delivered to segregated neuronal subpopulations. Harnessing STDP, PMCS leads to changes of the synaptic network structure. We found that the PMCS-induced changes of the network structure depend on both the phase lags between stimuli and the shape of individual stimuli. Single-pulse stimuli and burst stimuli with low intraburst frequency down-regulate synapses between neurons receiving stimuli simultaneously. In contrast, burst stimuli with high intraburst frequency up-regulate these synapses. We derive theoretical approximations of the stimulation-induced network structure. This enables us to formulate stimulation strategies for inducing a variety of network structures. Our results provide testable hypotheses for future pre-clinical and clinical studies and suggest that periodic multichannel stimulation may be suitable for reshaping plastic neuronal networks to counteract pathological synaptic connectivity. Furthermore, we provide novel insight on how the stimulus type may affect the long-lasting outcome of conventional DBS. This may strongly impact parameter adjustment procedures for clinical DBS, which, so far, primarily focused on acute effects of stimulation.

Dynamics of phase oscillator networks with synaptic weight and structural plasticity

We study the dynamics of Kuramoto oscillator networks with two distinct adaptation processes, one varying the coupling strengths and the other altering the network structure. Such systems model certain networks of oscillatory neurons where the neuronal dynamics, synaptic weights, and network structure interact with and shape each other. We model synaptic weight adaptation with spike-timing-dependent plasticity (STDP) that runs on a longer time scale than neuronal spiking. Structural changes that include addition and elimination of contacts occur at yet a longer time scale than the weight adaptations. First, we study the steady-state dynamics of Kuramoto networks that are bistable and can settle in synchronized or desynchronized states. To compare the impact of adding structural plasticity, we contrast the network with only STDP to one with a combination of STDP and structural plasticity. We show that the inclusion of structural plasticity optimizes the synchronized state of a network by allowing for synchronization with fewer links than a network with STDP alone. With non-identical units in the network, the addition of structural plasticity leads to the emergence of correlations between the oscillators' natural frequencies and node degrees. In the desynchronized regime, the structural plasticity decreases the number of contacts, leading to a sparse network. In this way, adding structural plasticity strengthens both synchronized and desynchronized states of a network. Second, we use desynchronizing coordinated reset stimulation and synchronizing periodic stimulation to induce desynchronized and synchronized states, respectively. Our findings indicate that a network with a combination of STDP and structural plasticity may require stronger and longer stimulation to switch between the states than a network with STDP only.

Critical Parameters in Dynamic Network Modeling of Sepsis

In this work, we propose a dynamical systems perspective on the modeling of sepsis and its organ-damaging consequences. We develop a functional two-layer network model for sepsis based upon the interaction of parenchymal cells and immune cells via cytokines, and the coevolutionary dynamics of parenchymal, immune cells, and cytokines. By means of the simple paradigmatic model of phase oscillators in a two-layer system, we analyze the emergence of organ threatening interactions between the dysregulated immune system and the parenchyma. We demonstrate that the complex cellular cooperation between parenchyma and stroma (immune layer) either in the physiological or in the pathological case can be related to dynamical patterns of the network. In this way we explain sepsis by the dysregulation of the healthy homeostatic state (frequency synchronized) leading to a pathological state (desynchronized or multifrequency cluster) in the parenchyma. We provide insight into the complex stabilizing and destabilizing interplay of parenchyma and stroma by determining critical interaction parameters. The coupled dynamics of parenchymal cells (metabolism) and nonspecific immune cells (response of the innate immune system) is represented by nodes of a duplex layer. Cytokine interaction is modeled by adaptive coupling weights between nodes representing immune cells (with fast adaptation timescale) and parenchymal cells (slow adaptation timescale), and between pairs of parenchymal and immune cells in the duplex network (fixed bidirectional coupling). The proposed model allows for a functional description of organ dysfunction in sepsis and the recurrence risk in a plausible pathophysiological context.

Long-Lasting Desynchronization of Plastic Neuronal Networks by Double-Random Coordinated Reset Stimulation

Hypersynchrony of neuronal activity is associated with several neurological disorders, including essential tremor and Parkinson's disease (PD). Chronic high-frequency deep brain stimulation (HF DBS) is the standard of care for medically refractory PD. Symptoms may effectively be suppressed by HF DBS, but return shortly after cessation of stimulation. Coordinated reset (CR) stimulation is a theory-based stimulation technique that was designed to specifically counteract neuronal synchrony by desynchronization. During CR, phase-shifted stimuli are delivered to multiple neuronal subpopulations. Computational studies on CR stimulation of plastic neuronal networks revealed long-lasting desynchronization effects obtained by down-regulating abnormal synaptic connectivity. This way, networks are moved into attractors of stable desynchronized states such that stimulation-induced desynchronization persists after cessation of stimulation. Preclinical and clinical studies confirmed corresponding long-lasting therapeutic and desynchronizing effects in PD. As PD symptoms are associated with different pathological synchronous rhythms, stimulation-induced long-lasting desynchronization effects should favorably be robust to variations of the stimulation frequency. Recent computational studies suggested that this robustness can be improved by randomizing the timings of stimulus deliveries. We study the long-lasting effects of CR stimulation with randomized stimulus amplitudes and/or randomized stimulus timing in networks of leaky integrate-and-fire (LIF) neurons with spike-timing-dependent plasticity. Performing computer simulations and analytical calculations, we study long-lasting desynchronization effects of CR with and without randomization of stimulus amplitudes alone, randomization of stimulus times alone as well as the combination of both. Varying the CR stimulation frequency (with respect to the frequency of abnormal target rhythm) and the number of separately stimulated neuronal subpopulations, we reveal parameter regions and related mechanisms where the two qualitatively different randomization mechanisms improve the robustness of long-lasting desynchronization effects of CR. In particular, for clinically relevant parameter ranges double-random CR stimulation, i.e., CR stimulation with the specific combination of stimulus amplitude randomization and stimulus time randomization, may outperform regular CR stimulation with respect to long-lasting desynchronization. In addition, our results provide the first evidence that an effective reduction of the overall stimulation current by stimulus amplitude randomization may improve the frequency robustness of long-lasting therapeutic effects of brain stimulation.

Thalamic bursts modulate cortical synchrony locally to switch between states of global functional connectivity in a cognitive task

Performing a cognitive task requires going through a sequence of functionally diverse stages. Although it is typically assumed that these stages are characterized by distinct states of cortical synchrony that are triggered by sub-cortical events, little reported evidence supports this hypothesis. To test this hypothesis, we first identified cognitive stages in single-trial MEG data of an associative recognition task, showing with a novel method that each stage begins with local modulations of synchrony followed by a state of directed functional connectivity. Second, we developed the first whole-brain model that can simulate cortical synchrony throughout a task. The model suggests that the observed synchrony is caused by thalamocortical bursts at the onset of each stage, targeted at cortical synapses and interacting with the structural anatomical connectivity. These findings confirm that cognitive stages are defined by distinct states of cortical synchrony and explains the network-level mechanisms necessary for reaching stage-dependent synchrony states.

Short- and Long-Range Connections Differentially Modulate the Dynamics and State of Small-World Networks

The human brain contains billions of neurons that flexibly interconnect to support local and global computational spans. As neuronal activity propagates through the neural medium, it approaches a critical state hedged between ordered and disordered system regimes. Recent work demonstrates that this criticality coincides with the small-world topology, a network arrangement that accommodates both local (subcritical) and global (supercritical) system properties. On one hand, operating near criticality is thought to offer several neurocomputational advantages, e.g., high-dynamic range, efficient information capacity, and information transfer fidelity. On the other hand, aberrations from the critical state have been linked to diverse pathologies of the brain, such as post-traumatic epileptiform seizures and disorders of consciousness. Modulation of brain activity, through neuromodulation, presents an attractive mode of treatment to alleviate such neurological disorders, but a tractable neural framework is needed to facilitate clinical progress. Using a variation on the generative small-world model of Watts and Strogatz and Kuramoto's model of coupled oscillators, we show that the topological and dynamical properties of the small-world network are divided into two functional domains based on the range of connectivity, and that these domains play distinct roles in shaping the behavior of the critical state. We demonstrate that short-range network connections shape the dynamics of the system, e.g., its volatility and metastability, whereas long-range connections drive the system state, e.g., a seizure. Together, these findings lend support to combinatorial neuromodulation approaches that synergistically normalize the system dynamic while mobilizing the system state.

Modeling Tumor Disease and Sepsis by Networks of Adaptively Coupled Phase Oscillators

In this study, we provide a dynamical systems perspective to the modelling of pathological states induced by tumors or infection. A unified disease model is established using the innate immune system as the reference point. We propose a two-layer network model for carcinogenesis and sepsis based upon the interaction of parenchymal cells and immune cells via cytokines, and the co-evolutionary dynamics of parenchymal, immune cells, and cytokines. Our aim is to show that the complex cellular cooperation between parenchyma and stroma (immune layer) in the physiological and pathological case can be qualitatively and functionally described by a simple paradigmatic model of phase oscillators. By this, we explain carcinogenesis, tumor progression, and sepsis by destabilization of the healthy homeostatic state (frequency synchronized), and emergence of a pathological state (desynchronized or multifrequency cluster). The coupled dynamics of parenchymal cells (metabolism) and nonspecific immune cells (reaction of innate immune system) are represented by nodes of a duplex layer. The cytokine interaction is modeled by adaptive coupling weights between the nodes representing the immune cells (with fast adaptation time scale) and the parenchymal cells (slow adaptation time scale) and between the pairs of parenchymal and immune cells in the duplex network (fixed bidirectional coupling). Thereby, carcinogenesis, organ dysfunction in sepsis, and recurrence risk can be described in a correct functional context.

Spike-Timing-Dependent Plasticity Mediated by Dopamine and its Role in Parkinson's Disease Pathophysiology

Parkinson's disease (PD) is a multi-systemic neurodegenerative brain disorder. Motor symptoms of PD are linked to the significant dopamine (DA) loss in substantia nigra pars compacta (SNc) followed by basal ganglia (BG) circuit dysfunction. Increasing experimental and computational evidence indicates that (synaptic) plasticity plays a key role in the emergence of PD-related pathological changes following DA loss. Spike-timing-dependent plasticity (STDP) mediated by DA provides a mechanistic model for synaptic plasticity to modify synaptic connections within the BG according to the neuronal activity. To shed light on how DA-mediated STDP can shape neuronal activity and synaptic connectivity in the PD condition, we reviewed experimental and computational findings addressing the modulatory effect of DA on STDP as well as other plasticity mechanisms and discussed their potential role in PD pathophysiology and related network dynamics and connectivity. In particular, reshaping of STDP profiles together with other plasticity-mediated processes following DA loss may abnormally modify synaptic connections in competing pathways of the BG. The cascade of plasticity-induced maladaptive or compensatory changes can impair the excitation-inhibition balance towards the BG output nuclei, leading to the emergence of pathological activity-connectivity patterns in PD. Pre-clinical, clinical as well as computational studies reviewed here provide an understanding of the impact of synaptic plasticity and other plasticity mechanisms on PD pathophysiology, especially PD-related network activity and connectivity, after DA loss. This review may provide further insights into the abnormal structure-function relationship within the BG contributing to the emergence of pathological states in PD. Specifically, this review is intended to provide detailed information for the development of computational network models for PD, serving as testbeds for the development and optimization of invasive and non-invasive brain stimulation techniques. Computationally derived hypotheses may accelerate the development of therapeutic stimulation techniques and potentially reduce the number of related animal experiments.

Long-Lasting Desynchronization Effects of Coordinated Reset Stimulation Improved by Random Jitters

Abnormally strong synchronized activity is related to several neurological disorders, including essential tremor, epilepsy, and Parkinson's disease. Chronic high-frequency deep brain stimulation (HF DBS) is an established treatment for advanced Parkinson's disease. To reduce the delivered integral electrical current, novel theory-based stimulation techniques such as coordinated reset (CR) stimulation directly counteract the abnormal synchronous firing by delivering phase-shifted stimuli through multiple stimulation sites. In computational studies in neuronal networks with spike-timing-dependent plasticity (STDP), it was shown that CR stimulation down-regulates synaptic weights and drives the network into an attractor of a stable desynchronized state. This led to desynchronization effects that outlasted the stimulation. Corresponding long-lasting therapeutic effects were observed in preclinical and clinical studies. Computational studies suggest that long-lasting effects of CR stimulation depend on the adjustment of the stimulation frequency to the dominant synchronous rhythm. This may limit clinical applicability as different pathological rhythms may coexist. To increase the robustness of the long-lasting effects, we study randomized versions of CR stimulation in networks of leaky integrate-and-fire neurons with STDP. Randomization is obtained by adding random jitters to the stimulation times and by shuffling the sequence of stimulation site activations. We study the corresponding long-lasting effects using analytical calculations and computer simulations. We show that random jitters increase the robustness of long-lasting effects with respect to changes of the number of stimulation sites and the stimulation frequency. In contrast, shuffling does not increase parameter robustness of long-lasting effects. Studying the relation between acute, acute after-, and long-lasting effects of stimulation, we find that both acute after- and long-lasting effects are strongly determined by the stimulation-induced synaptic reshaping, whereas acute effects solely depend on the statistics of administered stimuli. We find that the stimulation duration is another important parameter, as effective stimulation only entails long-lasting effects after a sufficient stimulation duration. Our results show that long-lasting therapeutic effects of CR stimulation with random jitters are more robust than those of regular CR stimulation. This might reduce the parameter adjustment time in future clinical trials and make CR with random jitters more suitable for treating brain disorders with abnormal synchronization in multiple frequency bands.

Generalized splay states in phase oscillator networks

Networks of coupled phase oscillators play an important role in the analysis of emergent collective phenomena. In this article, we introduce generalized m-splay states constituting a special subclass of phase-locked states with vanishing mth order parameter. Such states typically manifest incoherent dynamics, and they often create high-dimensional families of solutions (splay manifolds). For a general class of phase oscillator networks, we provide explicit linear stability conditions for splay states and exemplify our results with the well-known Kuramoto-Sakaguchi model. Importantly, our stability conditions are expressed in terms of just a few observables such as the order parameter or the trace of the Jacobian. As a result, these conditions are simple and applicable to networks of arbitrary size. We generalize our findings to phase oscillators with inertia and adaptively coupled phase oscillator models.

NetworkDynamics.jl-Composing and simulating complex networks in Julia

NetworkDynamics.jl is an easy-to-use and computationally efficient package for simulating heterogeneous dynamical systems on complex networks, written in Julia, a high-level, high-performance, dynamic programming language. By combining state-of-the-art solver algorithms from DifferentialEquations.jl with efficient data structures, NetworkDynamics.jl achieves top performance while supporting advanced features such as events, algebraic constraints, time delays, noise terms, and automatic differentiation.

Multistability in a star network of Kuramoto-type oscillators with synaptic plasticity

We analyze multistability in a star-type network of phase oscillators with coupling weights governed by phase-difference-dependent plasticity. It is shown that a network with N leaves can evolve into \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2^N$$\end{document}2N various asymptotic states, characterized by different values of the coupling strength between the hub and the leaves. Starting from the simple case of two coupled oscillators, we develop an analytical approach based on two small parameters \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon$$\end{document}ε and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μ, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon$$\end{document}ε is the ratio of the time scales of the phase variables and synaptic weights, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu$$\end{document}μ defines the sharpness of the plasticity boundary function. The limit \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mu \rightarrow 0$$\end{document}μ→0 corresponds to a hard boundary. The analytical results obtained on the model of two oscillators are generalized for multi-leaf star networks. Multistability with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2^N$$\end{document}2N various asymptotic states is numerically demonstrated for one-, two-, three- and nine-leaf star-type networks.

Long-Lasting Desynchronization of Plastic Neural Networks by Random Reset Stimulation

Excessive neuronal synchrony is a hallmark of neurological disorders such as epilepsy and Parkinson's disease. An established treatment for medically refractory Parkinson's disease is high-frequency (HF) deep brain stimulation (DBS). However, symptoms return shortly after cessation of HF-DBS. Recently developed decoupling stimulation approaches, such as Random Reset (RR) stimulation, specifically target pathological connections to achieve long-lasting desynchronization. During RR stimulation, a temporally and spatially randomized stimulus pattern is administered. However, spatial randomization, as presented so far, may be difficult to realize in a DBS-like setup due to insufficient spatial resolution. Motivated by recently developed segmented DBS electrodes with multiple stimulation sites, we present a RR stimulation protocol that copes with the limited spatial resolution of currently available depth electrodes for DBS. Specifically, spatial randomization is realized by delivering stimuli simultaneously to L randomly selected stimulation sites out of a total of M stimulation sites, which will be called L/M-RR stimulation. We study decoupling by L/M-RR stimulation in networks of excitatory integrate-and-fire neurons with spike-timing dependent plasticity by means of theoretical and computational analysis. We find that L/M-RR stimulation yields parameter-robust decoupling and long-lasting desynchronization. Furthermore, our theory reveals that strong high-frequency stimulation is not suitable for inducing long-lasting desynchronization effects. As a consequence, low and high frequency L/M-RR stimulation affect synaptic weights in qualitatively different ways. Our simulations confirm these predictions and show that qualitative differences between low and high frequency L/M-RR stimulation are present across a wide range of stimulation parameters, rendering stimulation with intermediate frequencies most efficient. Remarkably, we find that L/M-RR stimulation does not rely on a high spatial resolution, characterized by the density of stimulation sites in a target area, corresponding to a large M. In fact, L/M-RR stimulation with low resolution performs even better at low stimulation amplitudes. Our results provide computational evidence that L/M-RR stimulation may present a way to exploit modern segmented lead electrodes for long-lasting therapeutic effects.

Desynchronization Transitions in Adaptive Networks

Adaptive networks change their connectivity with time, depending on their dynamical state. While synchronization in structurally static networks has been studied extensively, this problem is much more challenging for adaptive networks. In this Letter, we develop the master stability approach for a large class of adaptive networks. This approach allows for reducing the synchronization problem for adaptive networks to a low-dimensional system, by decoupling topological and dynamical properties. We show how the interplay between adaptivity and network structure gives rise to the formation of stability islands. Moreover, we report a desynchronization transition and the emergence of complex partial synchronization patterns induced by an increasing overall coupling strength. We illustrate our findings using adaptive networks of coupled phase oscillators and FitzHugh-Nagumo neurons with synaptic plasticity.

Impact of number of stimulation sites on long-lasting desynchronization effects of coordinated reset stimulation

Excessive neuronal synchrony is a hallmark of several neurological disorders, e.g., Parkinson's disease. An established treatment for medically refractory Parkinson's disease is high-frequency deep brain stimulation. However, it provides only acute relief, and symptoms return shortly after cessation of stimulation. A theory-based approach called coordinated reset (CR) has shown great promise in achieving long-lasting effects. During CR stimulation, phase-shifted stimuli are delivered to multiple stimulation sites to counteract neuronal synchrony. Computational studies in plastic neuronal networks reported that synaptic weights reduce during stimulation, which may cause sustained structural changes leading to stabilized desynchronized activity even after stimulation ceases. Corresponding long-lasting effects were found in recent preclinical and clinical studies. We study long-lasting desynchronization by CR stimulation in excitatory recurrent neuronal networks of integrate-and-fire neurons with spike-timing-dependent plasticity (STDP). We focus on the impact of the stimulation frequency and the number of stimulation sites on long-lasting effects. We compare theoretical predictions to simulations of plastic neuronal networks. Our results are important regarding CR calibration for two reasons. We reveal that long-lasting effects become most pronounced when stimulation parameters are adjusted to the characteristics of STDP-rather than to neuronal frequency characteristics. This is in contrast to previous studies where the CR frequency was adjusted to the dominant neuronal rhythm. In addition, we reveal a nonlinear dependence of long-lasting effects on the number of stimulation sites and the CR frequency. Intriguingly, optimal long-lasting desynchronization does not require larger numbers of stimulation sites.

Two paradigmatic scenarios for inverse stochastic resonance

Inverse stochastic resonance comprises a nonlinear response of an oscillatory system to noise where the frequency of noise-perturbed oscillations becomes minimal at an intermediate noise level. We demonstrate two generic scenarios for inverse stochastic resonance by considering a paradigmatic model of two adaptively coupled stochastic active rotators whose local dynamics is close to a bifurcation threshold. In the first scenario, shown for the two rotators in the excitable regime, inverse stochastic resonance emerges due to a biased switching between the oscillatory and the quasi-stationary metastable states derived from the attractors of the noiseless system. In the second scenario, illustrated for the rotators in the oscillatory regime, inverse stochastic resonance arises due to a trapping effect associated with a noise-enhanced stability of an unstable fixed point. The details of the mechanisms behind the resonant effect are explained in terms of slow-fast analysis of the corresponding noiseless systems.

Birth and Stabilization of Phase Clusters by Multiplexing of Adaptive Networks

We propose a concept to generate and stabilize diverse partial synchronization patterns (phase clusters) in adaptive networks which are widespread in neuroscience and social sciences, as well as biology, engineering, and other disciplines. We show by theoretical analysis and computer simulations that multiplexing in a multilayer network with symmetry can induce various stable phase cluster states in a situation where they are not stable or do not even exist in the single layer. Further, we develop a method for the analysis of Laplacian matrices of multiplex networks which allows for insight into the spectral structure of these networks enabling a reduction to the stability problem of single layers. We employ the multiplex decomposition to provide analytic results for the stability of the multilayer patterns. As local dynamics we use the paradigmatic Kuramoto phase oscillator, which is a simple generic model and has been successfully applied in the modeling of synchronization phenomena in a wide range of natural and technological systems.

Frequency cluster formation and slow oscillations in neural populations with plasticity

We report the phenomenon of frequency clustering in a network of Hodgkin-Huxley neurons with spike timing-dependent plasticity. The clustering leads to a splitting of a neural population into a few groups synchronized at different frequencies. In this regime, the amplitude of the mean field undergoes low-frequency modulations, which may contribute to the mechanism of the emergence of slow oscillations of neural activity observed in spectral power of local field potentials or electroencephalographic signals at high frequencies. In addition to numerical simulations of such multi-clusters, we investigate the mechanisms of the observed phenomena using the simplest case of two clusters. In particular, we propose a phenomenological model which describes the dynamics of two clusters taking into account the adaptation of coupling weights. We also determine the set of plasticity functions (update rules), which lead to multi-clustering.

Hierarchical frequency clusters in adaptive networks of phase oscillators

Adaptive dynamical networks appear in various real-word systems. One of the simplest phenomenological models for investigating basic properties of adaptive networks is the system of coupled phase oscillators with adaptive couplings. In this paper, we investigate the dynamics of this system. We extend recent results on the appearance of hierarchical frequency multiclusters by investigating the effect of the time scale separation. We show that the slow adaptation in comparison with the fast phase dynamics is necessary for the emergence of the multiclusters and their stability. Additionally, we study the role of double antipodal clusters, which appear to be unstable for all considered parameter values. We show that such states can be observed for a relatively long time, i.e., they are metastable. A geometrical explanation for such an effect is based on the emergence of a heteroclinic orbit.

Plasticity facilitates pattern selection of networks of chemical oscillations

Rotating wave synchronization patterns are explored with a ring of 20 electrochemical oscillators during nickel electrodissolution in sulfuric acid. With desynchronized initial states, coupling alone yields predominance of nonrotating solutions, i.e., in-phase synchronization. An experimental technique is presented in which, through a combination of temporary alterations in topology, the application of global feedback provides rotational solutions. With phase repulsive global feedback, the in-phase synchronization is destabilized and a rotating wave is obtained. This feedback induced rotating wave can be employed to establish an initial condition for the rotating wave with coupling only. Higher order rotating solutions with 2, 3, and 4 waves corotating around the ring are observed, where the initial conditions are generated by temporary network rewiring to a structure with 2, 3, and 4 loops, respectively, and by global feedback. The experimental observations are supported by numerical simulations with a phase model. The results indicate that while network plasticity is thought to be significant in the operation of neural systems, it can also play a role in pattern selection of chemical systems.

Effects of synaptic and myelin plasticity on learning in a network of Kuramoto phase oscillators

Models of learning typically focus on synaptic plasticity. However, learning is the result of both synaptic and myelin plasticity. Specifically, synaptic changes often co-occur and interact with myelin changes, leading to complex dynamic interactions between these processes. Here, we investigate the implications of these interactions for the coupling behavior of a system of Kuramoto oscillators. To that end, we construct a fully connected, one-dimensional ring network of phase oscillators whose coupling strength (reflecting synaptic strength) as well as conduction velocity (reflecting myelination) are each regulated by a Hebbian learning rule. We evaluate the behavior of the system in terms of structural (pairwise connection strength and conduction velocity) and functional connectivity (local and global synchronization behavior). We find that adaptive myelination is able to both functionally decouple structurally connected oscillators as well as to functionally couple structurally disconnected oscillators. With regard to the latter, we find that for conditions in which a system limited to synaptic plasticity develops two distinct clusters both structurally and functionally, additional adaptive myelination allows for functional communication across these structural clusters. These results confirm that network states following learning may be different when myelin plasticity is considered in addition to synaptic plasticity, pointing toward the relevance of integrating both factors in computational models of learning.

Synchronization in time-varying random networks with vanishing connectivity

A sufficiently connected topology linking the constituent units of a complex system is usually seen as a prerequisite for the emergence of collective phenomena such as synchronization. We present a random network of heterogeneous phase oscillators in which the links mediating the interactions are constantly rearranged with a characteristic timescale and, possibly, an extremely low instantaneous connectivity. We show that with strong coupling and sufficiently fast rewiring the network reaches partial synchronization even in the vanishing connectivity limit. In particular, we provide an approximate analytical argument, based on the comparison between the different characteristic timescales of our system in the low connectivity regime, which is able to predict the transition to synchronization threshold with satisfactory precision beyond the formal fast rewiring limit. We interpret our results as a qualitative mechanism for emergence of consensus in social communities. In particular, our result suggest that groups of individuals are capable of aligning their opinions under extremely sparse exchanges of views, which is reminiscent of fast communications that take place in the modern social media. Our results may also be relevant to characterize the onset of collective behavior in engineered systems of mobile units with limited wireless capabilities.

Propagation delays determine neuronal activity and synaptic connectivity patterns emerging in plastic neuronal networks

In plastic neuronal networks, the synaptic strengths are adapted to the neuronal activity. Specifically, spike-timing-dependent plasticity (STDP) is a fundamental mechanism that modifies the synaptic strengths based on the relative timing of pre- and postsynaptic spikes, taking into account the spikes' temporal order. In many studies, propagation delays were neglected to avoid additional dynamic complexity or computational costs. So far, networks equipped with a classic STDP rule typically rule out bidirectional couplings (i.e., either loops or uncoupled states) and are, hence, not able to reproduce fundamental experimental findings. In this review paper, we consider additional features, e.g., extensions of the classic STDP rule or additional aspects like noise, in order to overcome the contradictions between theory and experiment. In addition, we review in detail recent studies showing that a classic STDP rule combined with realistic propagation patterns is able to capture relevant experimental findings. In two coupled oscillatory neurons with propagation delays, bidirectional synapses can be preserved and potentiated. This result also holds for large networks of type-II phase oscillators. In addition, not only the mean of the initial distribution of synaptic weights, but also its standard deviation crucially determines the emergent structural connectivity, i.e., the mean final synaptic weight, the number of two-neuron loops, and the symmetry of the final connectivity pattern. The latter is affected by the firing rates, where more symmetric synaptic configurations emerge at higher firing rates. Finally, we discuss these findings in the context of the computational neuroscience-based development of desynchronizing brain stimulation techniques.

Beta-Rhythm Oscillations and Synchronization Transition in Network Models of Izhikevich Neurons: Effect of Topology and Synaptic Type

Despite their significant functional roles, beta-band oscillations are least understood. Synchronization in neuronal networks have attracted much attention in recent years with the main focus on transition type. Whether one obtains explosive transition or a continuous transition is an important feature of the neuronal network which can depend on network structure as well as synaptic types. In this study we consider the effect of synaptic interaction (electrical and chemical) as well as structural connectivity on synchronization transition in network models of Izhikevich neurons which spike regularly with beta rhythms. We find a wide range of behavior including continuous transition, explosive transition, as well as lack of global order. The stronger electrical synapses are more conducive to synchronization and can even lead to explosive synchronization. The key network element which determines the order of transition is found to be the clustering coefficient and not the small world effect, or the existence of hubs in a network. These results are in contrast to previous results which use phase oscillator models such as the Kuramoto model. Furthermore, we show that the patterns of synchronization changes when one goes to the gamma band. We attribute such a change to the change in the refractory period of Izhikevich neurons which changes significantly with frequency.

Delay-Induced Multistability and Loop Formation in Neuronal Networks with Spike-Timing-Dependent Plasticity

Spike-timing-dependent plasticity (STDP) adjusts synaptic strengths according to the precise timing of pre- and postsynaptic spike pairs. Theoretical and computational studies have revealed that STDP may contribute to the emergence of a variety of structural and dynamical states in plastic neuronal populations. In this manuscript, we show that by incorporating dendritic and axonal propagation delays in recurrent networks of oscillatory neurons, the asymptotic connectivity displays multistability, where different structures emerge depending on the initial distribution of the synaptic strengths. In particular, we show that the standard deviation of the initial distribution of synaptic weights, besides its mean, determines the main properties of the emergent structural connectivity such as the mean final synaptic weight, the number of two-neuron loops and the symmetry of the final structure. We also show that the firing rates of the neurons affect the evolution of the network, and a more symmetric configuration of the synapses emerges at higher firing rates. We justify the network results based on a two-neuron framework and show how the results translate to large recurrent networks.

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.

Emergence of local synchronization in neuronal networks with adaptive couplings

Local synchronization, both prolonged and transient, of oscillatory neuronal behavior in cortical networks plays a fundamental role in many aspects of perception and cognition. Here we study networks of Hindmarsh-Rose neurons with a new type of adaptive coupling, and show that these networks naturally produce both permanent and transient synchronization of local clusters of neurons. These deterministic systems exhibit complex dynamics with 1/f^η power spectra, which appears to be a consequence of a novel form of self-organized criticality.

The stability of fixed points for a Kuramoto model with Hebbian interactions

We consider a variation of the Kuramoto model with dynamic coupling, where the coupling strengths are allowed to evolve in response to the phase difference between the oscillators, a model first considered by Ha, Noh, and Park. We demonstrate that the fixed points of this model, as well as their stability, can be completely expressed in terms of the fixed points and stability of the analogous classical Kuramoto problem where the coupling strengths are fixed to a constant (the same for all edges). In particular, for the "all-to-all" network, where the underlying graph is the complete graph, the problem reduces to the problem of understanding the fixed points and stability of the all-to-all Kuramoto model with equal edge weights, a problem that is well understood.

Closed-loop deep brain stimulation by pulsatile delayed feedback with increased gap between pulse phases

Computationally it was shown that desynchronizing delayed feedback stimulation methods are effective closed-loop techniques for the control of synchronization in ensembles of interacting oscillators. We here computationally design stimulation signals for electrical stimulation of neuronal tissue that preserve the desynchronizing delayed feedback characteristics and comply with mandatory charge deposit-related safety requirements. For this, the amplitude of the high-frequency (HF) train of biphasic charge-balanced pulses used by the standard HF deep brain stimulation (DBS) is modulated by the smooth feedback signals. In this way we combine the desynchronizing delayed feedback approach with the HF DBS technique. We show that such a pulsatile delayed feedback stimulation can effectively and robustly desynchronize a network of model neurons comprising subthalamic nucleus and globus pallidus external and suggest this approach for desynchronizing closed-loop DBS. Intriguingly, an interphase gap introduced between the recharging phases of the charge-balanced biphasic pulses can significantly improve the stimulation-induced desynchronization and reduce the amount of the administered stimulation. In view of the recent experimental and clinical studies indicating a superiority of the closed-loop DBS to open-loop HF DBS, our results may contribute to a further development of effective stimulation methods for the treatment of neurological disorders characterized by abnormal neuronal synchronization.

Spatiotemporal dynamics of the Kuramoto-Sakaguchi model with time-dependent connectivity

We study the dynamics of the paradigmatic Kuramoto-Sakaguchi model of identical coupled phase oscillators with various kinds of time-dependent connectivity using Eulerian discretization. We explore the parameter spaces for various types of collective states using the phase plots of the two statistical quantities, namely, the strength of incoherence and the discontinuity measure. In the quasistatic limit of the changing of coupling range, we observe how the system relaxes from one state to another and identify a few interesting collective dynamical states along the way. Under a sinusoidal change of the coupling range, the global order parameter characterizing the degree of synchronization in the system is shown to undergo a hysteresis with the coupling range. We also study the low-dimensional spatiotemporal dynamics of the local order parameter in the continuum limit using the recently developed Ott-Antonsen ansatz and justify some of our numerical results. In particular, we identify an intrinsic time scale of the Kuramoto system and show that the simulations exhibit two distinct kinds of qualitative behavior in two cases when the time scale associated with the switching of the coupling radius is very large compared to the intrinsic time scale and when it is comparable to the intrinsic time scale.

Anti-kindling Induced by Two-Stage Coordinated Reset Stimulation with Weak Onset Intensity

Abnormal neuronal synchrony plays an important role in a number of brain diseases. To specifically counteract abnormal neuronal synchrony by desynchronization, Coordinated Reset (CR) stimulation, a spatiotemporally patterned stimulation technique, was designed with computational means. In neuronal networks with spike timing–dependent plasticity CR stimulation causes a decrease of synaptic weights and finally anti-kindling, i.e., unlearning of abnormally strong synaptic connectivity and abnormal neuronal synchrony. Long-lasting desynchronizing aftereffects of CR stimulation have been verified in pre-clinical and clinical proof of concept studies. In general, for different neuromodulation approaches, both invasive and non-invasive, it is desirable to enable effective stimulation at reduced stimulation intensities, thereby avoiding side effects. For the first time, we here present a two-stage CR stimulation protocol, where two qualitatively different types of CR stimulation are delivered one after another, and the first stage comes at a particularly weak stimulation intensity. Numerical simulations show that a two-stage CR stimulation can induce the same degree of anti-kindling as a single-stage CR stimulation with intermediate stimulation intensity. This stimulation approach might be clinically beneficial in patients suffering from brain diseases characterized by abnormal neuronal synchrony where a first treatment stage should be performed at particularly weak stimulation intensities in order to avoid side effects. This might, e.g., be relevant in the context of acoustic CR stimulation in tinnitus patients with hyperacusis or in the case of electrical deep brain CR stimulation with sub-optimally positioned leads or side effects caused by stimulation of the target itself. We discuss how to apply our method in first in man and proof of concept studies.

Noise-enhanced coupling between two oscillators with long-term plasticity

Spike timing-dependent plasticity is a fundamental adaptation mechanism of the nervous system. It induces structural changes of synaptic connectivity by regulation of coupling strengths between individual cells depending on their spiking behavior. As a biophysical process its functioning is constantly subjected to natural fluctuations. We study theoretically the influence of noise on a microscopic level by considering only two coupled neurons. Adopting a phase description for the neurons we derive a two-dimensional system which describes the averaged dynamics of the coupling strengths. We show that a multistability of several coupling configurations is possible, where some configurations are not found in systems without noise. Intriguingly, it is possible that a strong bidirectional coupling, which is not present in the noise-free situation, can be stabilized by the noise. This means that increased noise, which is normally expected to desynchronize the neurons, can be the reason for an antagonistic response of the system, which organizes itself into a state of stronger coupling and counteracts the impact of noise. This mechanism, as well as a high potential for multistability, is also demonstrated numerically for a coupled pair of Hodgkin-Huxley neurons.

Augmented brain function by coordinated reset stimulation with slowly varying sequences

Several brain disorders are characterized by abnormally strong neuronal synchrony. Coordinated Reset (CR) stimulation was developed to selectively counteract abnormal neuronal synchrony by desynchronization. For this, phase resetting stimuli are delivered to different subpopulations in a timely coordinated way. In neural networks with spike timing-dependent plasticity CR stimulation may eventually lead to an anti-kindling, i.e., an unlearning of abnormal synaptic connectivity and abnormal synchrony. The spatiotemporal sequence by which all stimulation sites are stimulated exactly once is called the stimulation site sequence, or briefly sequence. So far, in simulations, pre-clinical and clinical applications CR was applied either with fixed sequences or rapidly varying sequences (RVS). In this computational study we show that appropriate repetition of the sequence with occasional random switching to the next sequence may significantly improve the anti-kindling effect of CR. To this end, a sequence is applied many times before randomly switching to the next sequence. This new method is called SVS CR stimulation, i.e., CR with slowly varying sequences. In a neuronal network with strong short-range excitatory and weak long-range inhibitory dynamic couplings SVS CR stimulation turns out to be superior to CR stimulation with fixed sequences or RVS.

Gradual emergence of spontaneous correlated brain activity during fading of general anesthesia in rats: Evidences from fMRI and local field potentials

Intrinsic brain activity is characterized by the presence of highly structured networks of correlated fluctuations between different regions of the brain. Such networks encompass different functions, whose properties are known to be modulated by the ongoing global brain state and are altered in several neurobiological disorders. In the present study, we induced a deep state of anesthesia in rats by means of a ketamine/medetomidine peritoneal injection, and analyzed the time course of the correlation between the brain activity in different areas while anesthesia spontaneously decreased over time. We compared results separately obtained from fMRI and local field potentials (LFPs) under the same anesthesia protocol, finding that while most profound phases of anesthesia can be described by overall sparse connectivity, stereotypical activity and poor functional integration, during lighter states different frequency-specific functional networks emerge, endowing the gradual restoration of structured large-scale activity seen during rest. Noteworthy, our in vivo results show that those areas belonging to the same functional network (the default-mode) exhibited sustained correlated oscillations around 10 Hz throughout the protocol, suggesting the presence of a specific functional backbone that is preserved even during deeper phases of anesthesia. Finally, the overall pattern of results obtained from both imaging and in vivo-recordings suggests that the progressive emergence from deep anesthesia is reflected by a corresponding gradual increase of organized correlated oscillations across the cortex.

•

Rat brain activity was recorded while anesthesia spontaneously decreased over time.

•

Fading of anesthesia modulates the overall pattern of brain structured co-activations.

•

Areas of the same network show frequency-specific coupling even in deep anesthesia.

•

Light anesthesia is characterized by the gradual emergence of large-scale connectivity.

•

Correlated band-limited oscillations distinguish between states and network properties.

The spacing principle for unlearning abnormal neuronal synchrony

Desynchronizing stimulation techniques were developed to specifically counteract abnormal neuronal synchronization relevant to several neurological and psychiatric disorders. The goal of our approach is to achieve an anti-kindling, where the affected neural networks unlearn abnormal synaptic connectivity and, hence, abnormal neuronal synchrony, by means of desynchronizing stimulation, in particular, Coordinated Reset (CR) stimulation. As known from neuroscience, psychology and education, learning effects can be enhanced by means of the spacing principle, i.e. by delivering repeated stimuli spaced by pauses as opposed to delivering a massed stimulus (in a single long stimulation session). To illustrate that the spacing principle may boost the anti-kindling effect of CR neuromodulation, in this computational study we carry this approach to extremes. To this end, we deliver spaced CR neuromodulation at particularly weak intensities which render permanently delivered CR neuromodulation ineffective. Intriguingly, spaced CR neuromodulation at these particularly weak intensities effectively induces an anti-kindling. In fact, the spacing principle enables the neuronal population to successively hop from one attractor to another one, finally approaching attractors characterized by down-regulated synaptic connectivity and synchrony. Our computational results might open up novel opportunities to effectively induce sustained desynchronization at particularly weak stimulation intensities, thereby avoiding side effects, e.g., in the case of deep brain stimulation.