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

Non-conventional deep brain stimulation in a network model of movement disorders

Conventional deep brain stimulation (DBS) for movement disorders is a well-established clinical treatment. Over the last few decades, over 200,000 people have been treated by DBS worldwide for several neurological conditions, including Parkinson's disease and Essential Tremor. DBS involves implanting electrodes into disorder-specific targets in the brain and applying an electric current. Although the hardware has developed in recent years, the clinically used stimulation pattern has remained as a regular frequency square pulse. Recent studies have suggested that phase-locking, coordinated reset or irregular patterns may be as or more effective at desynchronising the pathological neural activity. Such studies have shown efficacy using detailed neuron models or highly simplified networks and considered one frequency band. We previously described a population level model which generates oscillatory activity in both the beta band (20 Hz) and the tremor band (4 Hz). Here we use this model to look at the impact of applying regular, irregular and phase dependent bursts of stimulation, and show how this influences both tremor- and beta-band activity. We found that bursts are as or more effective at suppressing the pathological oscillations compared to continuous DBS. Importantly however, at higher amplitudes we found that the stimulus drove the network activity, as seen previously. Strikingly, this suppression was most apparent for the tremor band oscillations, with beta band pathological activity being more resistant to the burst stimulation compared to continuous, conventional DBS. Furthermore, our simulations showed that phase-locked bursts of stimulation did not convey much improvement on regular bursts of oscillation. Using a genetic algorithm optimisation approach to find the best stimulation parameters for regular, irregular and phase-locked bursts, we confirmed that tremor band oscillations could be more readily suppressed. Our results allow exploration of stimulation mechanisms at the network level to formulate testable predictions regarding parameter settings in DBS.

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.

Synaptic network structure shapes cortically evoked spatio-temporal responses of STN and GPe neurons in a computational model

Introduction

The basal ganglia (BG) are involved in motor control and play an essential role in movement disorders such as hemiballismus, dystonia, and Parkinson's disease. Neurons in the motor part of the BG respond to passive movement or stimulation of different body parts and to stimulation of corresponding cortical regions. Experimental evidence suggests that the BG are organized somatotopically, i.e., specific areas of the body are associated with specific regions in the BG nuclei. Signals related to the same body part that propagate along different pathways converge onto the same BG neurons, leading to characteristic shapes of cortically evoked responses. This suggests the existence of functional channels that allow for the processing of different motor commands or information related to different body parts in parallel. Neurological disorders such as Parkinson's disease are associated with pathological activity in the BG and impaired synaptic connectivity, together with reorganization of somatotopic maps. One hypothesis is that motor symptoms are, at least partly, caused by an impairment of network structure perturbing the organization of functional channels.

Methods

We developed a computational model of the STN-GPe circuit, a central part of the BG. By removing individual synaptic connections, we analyzed the contribution of signals propagating along different pathways to cortically evoked responses. We studied how evoked responses are affected by systematic changes in the network structure. To quantify the BG's organization in the form of functional channels, we suggested a two-site stimulation protocol.

Results

Our model reproduced the cortically evoked responses of STN and GPe neurons and the contributions of different pathways suggested by experimental studies. Cortical stimulation evokes spatio-temporal response patterns that are linked to the underlying synaptic network structure. Our two-site stimulation protocol yielded an approximate functional channel width.

Discussion/conclusion

The presented results provide insight into the organization of BG synaptic connectivity, which is important for the development of computational models. The synaptic network structure strongly affects the processing of cortical signals and may impact the generation of pathological rhythms. Our work may motivate further experiments to analyze the network structure of BG nuclei and their organization in functional channels.

How to entrain a selected neuronal rhythm but not others: open-loop dithered brain stimulation for selective entrainment

Objective.While brain stimulation therapies such as deep brain stimulation for Parkinson's disease (PD) can be effective, they have yet to reach their full potential across neurological disorders. Entraining neuronal rhythms using rhythmic brain stimulation has been suggested as a new therapeutic mechanism to restore neurotypical behaviour in conditions such as chronic pain, depression, and Alzheimer's disease. However, theoretical and experimental evidence indicate that brain stimulation can also entrain neuronal rhythms at sub- and super-harmonics, far from the stimulation frequency. Crucially, these counterintuitive effects could be harmful to patients, for example by triggering debilitating involuntary movements in PD. We therefore seek a principled approach to selectively promote rhythms close to the stimulation frequency, while avoiding potential harmful effects by preventing entrainment at sub- and super-harmonics.Approach.Our open-loop approach to selective entrainment, dithered stimulation, consists in adding white noise to the stimulation period.Main results.We theoretically establish the ability of dithered stimulation to selectively entrain a given brain rhythm, and verify its efficacy in simulations of uncoupled neural oscillators, and networks of coupled neural oscillators. Furthermore, we show that dithered stimulation can be implemented in neurostimulators with limited capabilities by toggling within a finite set of stimulation frequencies.Significance.Likely implementable across a variety of existing brain stimulation devices, dithering-based selective entrainment has potential to enable new brain stimulation therapies, as well as new neuroscientific research exploiting its ability to modulate higher-order entrainment.

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.

No robust online effects of transcranial direct current stimulation on corticospinal excitability

Transcranial direct current stimulation (tDCS) has been used for over twenty years to modulate cortical (particularly motor corticospinal) excitability both during (online) and outlasting (offline) the stimulation, with the former effects associated to the latter. However, tDCS effects are highly variable, partially because stimulation intensity is commonly not adjusted individually (in contrast to transcranial magnetic stimulation, TMS). In Experiment 1, we therefore explored an empirical approach of personalizing tDCS intensity for the primary motor cortex (M1) based on dose-response curves (DRCs), individually relating tDCS Intensity (in steps from 0.3 to 2.0 mA) and Polarity (anodal, cathodal) to the online modulation of concurrent TMS motor evoked potentials (MEP), assessing DRC reliability across two separate days. No robust DRCs could be observed, neither at the individual nor at the group level, with the only robust effect being a (paradoxical) MEP facilitation during cathodal tDCS at 2.0 mA, but no modulation at traditional intensities of or near 1 mA. In Experiment 2, we therefore attempted to replicate the classical bidirectional online MEP modulation during 1 mA tDCS that had been reported by several of the early seminal tDCS papers. We either closely recreated stimulation parameters and temporal protocol of these original studies (Experiment 2A) or slightly modernized them according to current standards (Experiment 2B). In neither experiment did we observed any significant online MEP modulation. We conclude that an empirical titration of individually effective tDCS intensities may not be feasible as online tDCS effects do not appear to be sufficiently robust.

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.

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.

Treatment Tone Spacing and Acute Effects of Acoustic Coordinated Reset Stimulation in Tinnitus Patients

Acoustic coordinated reset (aCR) therapy for tinnitus aims to desynchronize neuronal populations in the auditory cortex that exhibit pathologically increased coincident firing. The original therapeutic paradigm involves fixed spacing of four low-intensity tones centered around the frequency of a tone matching the tinnitus pitch, f T , but it is unknown whether these tones are optimally spaced for induction of desynchronization. Computational and animal studies suggest that stimulus amplitude, and relatedly, spatial stimulation profiles, of coordinated reset pulses can have a major impact on the degree of desynchronization achievable. In this study, we transform the tone spacing of aCR into a scale that takes into account the frequency selectivity of the auditory system at each therapeutic tone's center frequency via a measure called the gap index. Higher gap indices are indicative of more loosely spaced aCR tones. The gap index was found to be a significant predictor of symptomatic improvement, with larger gap indices, i.e., more loosely spaced aCR tones, resulting in reduction of tinnitus loudness and annoyance scores in the acute stimulation setting. A notable limitation of this study is the intimate relationship of hearing impairment with the gap index. Particularly, the shape of the audiogram in the vicinity of the tinnitus frequency can have a major impact on tone spacing. However, based on our findings we suggest hypotheses-based experimental protocols that may help to disentangle the impact of hearing loss and tone spacing on clinical outcome, to assess the electrophysiologic correlates of clinical improvement, and to elucidate the effects following chronic rather than acute stimulation.