Predicting the effects of deep brain stimulation using a reduced coupled oscillator model

Deep brain stimulation (DBS) is known to be an effective treatment for a variety of neurological disorders, including Parkinson’s disease and essential tremor (ET). At present, it involves administering a train of pulses with constant frequency via electrodes implanted into the brain. New ‘closed-loop’ approaches involve delivering stimulation according to the ongoing symptoms or brain activity and have the potential to provide improvements in terms of efficiency, efficacy and reduction of side effects. The success of closed-loop DBS depends on being able to devise a stimulation strategy that minimizes oscillations in neural activity associated with symptoms of motor disorders. A useful stepping stone towards this is to construct a mathematical model, which can describe how the brain oscillations should change when stimulation is applied at a particular state of the system. Our work focuses on the use of coupled oscillators to represent neurons in areas generating pathological oscillations. Using a reduced form of the Kuramoto model, we analyse how a patient should respond to stimulation when neural oscillations have a given phase and amplitude. We predict that, provided certain conditions are satisfied, the best stimulation strategy should be phase specific but also that stimulation should have a greater effect if applied when the amplitude of brain oscillations is lower. We compare this surprising prediction with data obtained from ET patients. In light of our predictions, we also propose a new hybrid strategy which effectively combines two of the strategies found in the literature, namely phase-locked and adaptive DBS. Author summary Deep brain stimulation (DBS) involves delivering electrical impulses to target sites within the brain and is a proven therapy for a variety of neurological disorders. Closed loop DBS is a promising new approach where stimulation is applied according to the state of a patient. Crucial to the success of this approach is being able to predict how a patient should respond to stimulation. Our work focusses on DBS as applied to patients with essential tremor (ET). On the basis of a theoretical model, which describes neurons as oscillators that respond to stimulation and have a certain tendency to synchronize, we provide predictions for how a patient should respond when stimulation is applied at a particular phase and amplitude of the ongoing tremor oscillations. Previous experimental studies of closed loop DBS provided stimulation either on the basis of ongoing phase or amplitude of pathological oscillations. Our study suggests how both of these measurements can be used to control stimulation. As part of this work, we also look for evidence for our theories in experimental data and find our predictions to be satisfied in one patient. The insights obtained from this work should lead to a better understanding of how to optimise closed loop DBS strategies.

[1]  G. Fink,et al.  Essential tremor and tremor in Parkinson's disease are associated with distinct ‘tremor clusters’ in the ventral thalamus , 2012, Experimental Neurology.

[2]  P. Brown,et al.  Adaptive Deep Brain Stimulation In Advanced Parkinson Disease , 2013, Annals of neurology.

[3]  B. Doiron,et al.  Axonal and synaptic failure suppress the transfer of firing rate oscillations, synchrony and information during high frequency deep brain stimulation , 2014, Neurobiology of Disease.

[4]  E. Ott,et al.  Low dimensional behavior of large systems of globally coupled oscillators. , 2008, Chaos.

[5]  Rafal Bogacz,et al.  Computational Models Describing Possible Mechanisms for Generation of Excessive Beta Oscillations in Parkinson’s Disease , 2015, PLoS Comput. Biol..

[6]  Marc Timme,et al.  Chaos in symmetric phase oscillator networks. , 2011, Physical review letters.

[7]  Jeff Moehlis,et al.  Phasic Burst Stimulation: A Closed-Loop Approach to Tuning Deep Brain Stimulation Parameters for Parkinson’s Disease , 2016, PLoS Comput. Biol..

[8]  M. Butz,et al.  Synchronized brain network associated with essential tremor as revealed by magnetoencephalography , 2009, Movement disorders : official journal of the Movement Disorder Society.

[9]  Boris Gutkin,et al.  Macroscopic phase-resetting curves for spiking neural networks. , 2017, Physical review. E.

[10]  Yoji Kawamura,et al.  Noise-induced synchronization and clustering in ensembles of uncoupled limit-cycle oscillators. , 2007, Physical review letters.

[11]  C. McIntyre,et al.  Uncovering the mechanism(s) of action of deep brain stimulation: activation, inhibition, or both , 2004, Clinical Neurophysiology.

[12]  Charles J. Wilson,et al.  Chaotic Desynchronization as the Therapeutic Mechanism of Deep Brain Stimulation , 2011, Front. Syst. Neurosci..

[13]  Alexander Münchau,et al.  Activity Parameters of Subthalamic Nucleus Neurons Selectively Predict Motor Symptom Severity in Parkinson's Disease , 2014, The Journal of Neuroscience.

[14]  Daniel B. Forger,et al.  Macroscopic models for networks of coupled biological oscillators , 2017, Science Advances.

[15]  H. Bergman,et al.  Pathological synchronization in Parkinson's disease: networks, models and treatments , 2007, Trends in Neurosciences.

[16]  Eve Marder,et al.  The Functional Consequences of Changes in the Strength and Duration of Synaptic Inputs to Oscillatory Neurons , 2003, The Journal of Neuroscience.

[17]  Carmen C Canavier,et al.  Phase-resetting as a tool of information transmission , 2015, Current Opinion in Neurobiology.

[18]  Ernest Barreto,et al.  Complete Classification of the Macroscopic Behavior of a Heterogeneous Network of Theta Neurons , 2013, Neural Computation.

[19]  R. Bogacz,et al.  Phase dependence of response curves to stimulation and their relationship: from a Wilson-Cowan model to essential tremor patient data , 2019, bioRxiv.

[20]  M. Rosenblum,et al.  Dynamics of weakly inhomogeneous oscillator populations: perturbation theory on top of Watanabe–Strogatz integrability , 2016, 1601.07170.

[21]  S. Haber,et al.  Closed-Loop Deep Brain Stimulation Is Superior in Ameliorating Parkinsonism , 2011, Neuron.

[22]  Jeff Moehlis,et al.  Optimal Chaotic Desynchronization for Neural Populations , 2014, SIAM J. Appl. Dyn. Syst..

[23]  F. Lenz,et al.  Patterns of bursting occurring in thalamic cells during parkinsonian tremor , 1998, Neuroscience.

[24]  S. Rehncrona,et al.  Effects of different electrical parameter settings on the intelligibility of speech in patients with Parkinson's disease treated with subthalamic deep brain stimulation , 2005, Movement disorders : official journal of the Movement Disorder Society.

[25]  Jeff Moehlis,et al.  Clustered Desynchronization from High-Frequency Deep Brain Stimulation , 2015, PLoS Comput. Biol..

[26]  Jonathan E. Rubin,et al.  High Frequency Stimulation of the Subthalamic Nucleus Eliminates Pathological Thalamic Rhythmicity in a Computational Model , 2004, Journal of Computational Neuroscience.

[27]  A. Lang,et al.  High‐frequency unilateral thalamic stimulation in the treatment of essential and parkinsonian tremor , 1997, Annals of neurology.

[28]  H. H. Ku Notes on the Use of Propagation of Error Formulas , 2010 .

[29]  Peter A. Tass,et al.  A model of desynchronizing deep brain stimulation with a demand-controlled coordinated reset of neural subpopulations , 2003, Biological Cybernetics.

[30]  Eduardo Martin Moraud,et al.  Predicting beta bursts from local field potentials to improve closed-loop DBS paradigms in Parkinson’s patients , 2018, 2018 40th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[31]  H. Shill,et al.  Pathology in essential tremor. , 2012, Parkinsonism & related disorders.

[32]  J. Cowan,et al.  Excitatory and inhibitory interactions in localized populations of model neurons. , 1972, Biophysical journal.

[33]  Karl J. Friston,et al.  Stimulating at the right time: phase-specific deep brain stimulation , 2016, Brain : a journal of neurology.

[34]  P. Brown,et al.  Phase dependent modulation of tremor amplitude in essential tremor through thalamic stimulation , 2013, Brain : a journal of neurology.

[35]  Stephen Coombes,et al.  Mathematical Frameworks for Oscillatory Network Dynamics in Neuroscience , 2015, The Journal of Mathematical Neuroscience.

[36]  P. Brown,et al.  Modulation of Long-Range Connectivity Patterns via Frequency-Specific Stimulation of Human Cortex , 2017, Current Biology.

[37]  Eric Shea-Brown,et al.  On the Phase Reduction and Response Dynamics of Neural Oscillator Populations , 2004, Neural Computation.

[38]  E. Foncke,et al.  Pathological gambling after bilateral subthalamic nucleus stimulation in Parkinson disease , 2006, Journal of Neurology, Neurosurgery & Psychiatry.

[39]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1990 .

[40]  Mason A Porter,et al.  Noise-induced synchronization, desynchronization, and clustering in globally coupled nonidentical oscillators. , 2013, Physical review. E, Statistical, nonlinear, and soft matter physics.

[41]  Carlo R. Laing,et al.  Understanding the dynamics of biological and neural oscillator networks through mean-field reductions: a review , 2019, 1902.05307.

[42]  P. Brown,et al.  Reduction in subthalamic 8–35 Hz oscillatory activity correlates with clinical improvement in Parkinson's disease , 2006, The European journal of neuroscience.

[43]  V. Booth,et al.  Collective phase response curves for heterogeneous coupled oscillators. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[44]  T. Sejnowski,et al.  Cholinergic Neuromodulation Changes Phase Response Curve Shape and Type in Cortical Pyramidal Neurons , 2008, PloS one.

[45]  Jan Raethjen,et al.  The oscillating central network of Essential tremor , 2012, Clinical Neurophysiology.

[46]  Yoji Kawamura,et al.  Noise-induced synchronization and clustering in ensembles of uncoupled limit-cycle oscillators (生命リズムと振動子ネットワーク) , 2007 .

[47]  Carlo R Laing,et al.  Derivation of a neural field model from a network of theta neurons. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.