Relative contributions of local cell and passing fiber activation and silencing to changes in thalamic fidelity during deep brain stimulation and lesioning: a computational modeling study

Computational study of associations between the synaptic conductance of STN and GPe and the development of Parkinson's disease

There is evidence that the subthalamic nucleus (STN) and globus pallidus pars externa (GPe) involve in the development of Parkinson's disease, a neurodegenerative disorder characterized by motor and non-motor symptoms and loss of dopaminergic neurons in which the error index (EI) in firing patterns is widely used to address the related issues. Whether and how this interaction mechanism of STN and GPe affects EI in Parkinson's disease is uncertain. To account for this, we propose a kind of basal ganglia-thalamic network model associated with Parkinson's disease coupled with neurons, and investigate the effect of synaptic conductance of STN and GPe on EI in this network, as well as their internal relationship under EI as an index. The results show a relationship like a piecewise function between the error index and the slope of the state transition function of synaptic conductance from STN to GPe ( g snge ) and from GPe to STN ( g gesn ). And there is an approximate negative correlation between EI and g gesn . Increasing g snge and decreasing g gesn can improve the fidelity of thalamus information transmission and alleviate Parkinson's disease effectively. These obtained results can give some theoretical evidence that the abnormal synaptic releases of STN and GPe may be the symptoms of the development of Parkinson's disease, and further enrich the understanding of the pathogenesis and treatment mechanism of Parkinson's disease.

Computational models advance deep brain stimulation for Parkinson's disease

Deep brain stimulation(DBS) has become an effective intervention for advanced Parkinson's disease(PD), but the exact mechanism of DBS is still unclear. In this review, we discuss the history of DBS, the anatomy and internal architecture of the basal ganglia (BG), the abnormal pathological changes of the BG in PD, and how computational models can help understand and advance DBS. We also describe two types of models: mathematical theoretical models and clinical predictive models. Mathematical theoretical models simulate neurons or neural networks of BG to shed light on the mechanistic principle underlying DBS, while clinical predictive models focus more on patients' outcomes, helping to adapt treatment plans for each patient and advance novel electrode designs. Finally, we provide insights and an outlook on future technologies.

Insights from a model based study on optimizing non invasive brain electrical stimulation for Parkinson's disease

Parkinson's Disease (PD) is a disorder in the central nervous system which includes symptoms such as tremor, rigidity, and Bradykinesia. Deep brain stimulation (DBS) is the most effective method to treat PD motor symptoms especially when the patient is not responsive to other treatments. However, its invasiveness and high risk, involving electrode implantation in the Basal Ganglia (BG), prompt recent research to emphasize non-invasive Transcranial Electrical Stimulation (TES). TES proves to be effective in treating some PD symptoms with inherent safety and no associated risks. This study explores the potential of using TES, to modify the firing pattern of cells in BG that are responsible for motor symptoms in PD. The research employs a mathematical model of the BG to examine the impact of applying TES to the brain. This is conducted using a realistic head model incorporating the Finite Element Method (FEM). According to our findings, the firing pattern associated with Parkinson's disease shifted towards a healthier firing pattern through the use of tACS. Employing an adaptive algorithm that continually monitored the behavior of BG cells (specifically, Globus Pallidus Pars externa (GPe)), we determined the optimal electrode number and placement to concentrate the current within the intended region. This resulted in a peak induced electric field of 1.9 v/m at the BG area. Our mathematical modeling together with precise finite element simulation of the brain and BG suggests that proposed method effectively mitigates Parkinsonian behavior in the BG cells. Furthermore, this approach ensures an improvement in the condition while adhering to all safety constraints associated with the current injection into the brain.

Towards an optimised deep brain stimulation using a large-scale computational network and realistic volume conductor model

Objective. Constructing a theoretical framework to improve deep brain stimulation (DBS) based on the neuronal spatiotemporal patterns of the stimulation-affected areas constitutes a primary target.Approach. We develop a large-scale biophysical network, paired with a realistic volume conductor model, to estimate theoretically efficacious stimulation protocols. Based on previously published anatomically defined structural connectivity, a biophysical basal ganglia-thalamo-cortical neuronal network is constructed using Hodgkin-Huxley dynamics. We define a new biomarker describing the thalamic spatiotemporal activity as a ratio of spiking vs. burst firing. The per cent activation of the different pathways is adapted in the simulation to minimise the differences of the biomarker with respect to its value under healthy conditions.Main results.This neuronal network reproduces spatiotemporal patterns that emerge in Parkinson's disease. Simulations of the fibre per cent activation for the defined biomarker propose desensitisation of pallido-thalamic synaptic efficacy, induced by high-frequency signals, as one possible crucial mechanism for DBS action. Based on this activation, we define both an optimal electrode position and stimulation protocol using pathway activation modelling.Significance. A key advantage of this research is that it combines different approaches, i.e. the spatiotemporal pattern with the electric field and axonal response modelling, to compute the optimal DBS protocol. By correlating the inherent network dynamics with the activation of white matter fibres, we obtain new insights into the DBS therapeutic action.

Closed-loop modulation of model parkinsonian beta oscillations based on CAR-fuzzy control algorithm

Closed-loop deep brain stimulation (DBS) can apply on-demand stimulation based on the feedback signal (e.g. beta band oscillation), which is deemed to lower side effects of clinically used open-loop DBS. To facilitate the application of model-based closed-loop DBS in clinical, studies must consider state variations, e.g., variation of desired signal with different movement conditions and variation of model parameters with time. This paper proposes to use the controlled autoregressive (CAR)-fuzzy control algorithm to modulate the pathological beta band (13-35 Hz) oscillation of a basal ganglia-cortex-thalamus model. The CAR model is used to identify the relationship between DBS frequency parameter and beta oscillation power. Then the error between the one-step-ahead predicted beta power of CAR model and the desired value is innovatively used as the input of fuzzy controller to calculate the next step stimulation frequency. Compared with 130 Hz open-loop DBS, the proposed closed-loop DBS method could lower the mean stimulation frequency to 74.04 Hz with similar beta oscillation suppression performance. The Mamdani fuzzy controller is selected because which could establish fuzzy controller rules according to human operation experience. Adding prediction module to closed-loop control improves the accuracy of fuzzy control, compared with proportional-integral control and fuzzy control, the proposed CAR-fuzzy control algorithm has higher tracking reliability, response speed and robustness.

Offline Learning of Closed-Loop Deep Brain Stimulation Controllers for Parkinson Disease Treatment

Deep brain stimulation (DBS) has shown great promise toward treating motor symptoms caused by Parkinson's disease (PD), by delivering electrical pulses to the Basal Ganglia (BG) region of the brain. However, DBS devices approved by the U.S. Food and Drug Administration (FDA) can only deliver continuous DBS (cDBS) stimuli at a fixed amplitude; this energy inefficient operation reduces battery lifetime of the device, cannot adapt treatment dynamically for activity, and may cause significant side-effects (e.g., gait impairment). In this work, we introduce an offline reinforcement learning (RL) framework, allowing the use of past clinical data to train an RL policy to adjust the stimulation amplitude in real time, with the goal of reducing energy use while maintaining the same level of treatment (i.e., control) efficacy as cDBS. Moreover, clinical protocols require the safety and performance of such RL controllers to be demonstrated ahead of deployments in patients. Thus, we also introduce an offline policy evaluation (OPE) method to estimate the performance of RL policies using historical data, before deploying them on patients. We evaluated our framework on four PD patients equipped with the RC+S DBS system, employing the RL controllers during monthly clinical visits, with the overall control efficacy evaluated by severity of symptoms (i.e., bradykinesia and tremor), changes in PD biomakers (i.e., local field potentials), and patient ratings. The results from clinical experiments show that our RL-based controller maintains the same level of control efficacy as cDBS, but with significantly reduced stimulation energy. Further, the OPE method is shown effective in accurately estimating and ranking the expected returns of RL controllers.

Basal ganglia network dynamics and function: Role of direct, indirect and hyper-direct pathways in action selection

Basal ganglia (BG) are a widely recognized neural basis for action selection, but its decision-making mechanism is still a difficult problem for researchers. Therefore, we constructed a spiking neural network inspired by the BG anatomical data. Simulation experiments were based on the principle of dis-inhibition and our functional hypothesis within the BG: the direct pathway, the indirect pathway, and the hyper-direct pathway of the BG jointly implement the initiation execution and termination of motor programs. Firstly, we studied the dynamic process of action selection with the network, which contained intra-group competition and inter-group competition. Secondly, we focused on the effects of the stimulus intensity and the proportion of excitation and inhibition on the GPi/SNr. The results suggested that inhibition and excitation shape action selection. They also explained why the firing rate of GPi/SNr did not continue to increase in the action-selection experiment. Finally, we discussed the experimental results with the functional hypothesis. Uniquely, this paper summarized the decision-making neural mechanism of action selection based on the direct pathway, the indirect pathway, and the hyper-direct pathway within BG.

Efficient suppression of parkinsonian beta oscillations in a closed-loop model of deep brain stimulation with amplitude modulation

Introduction

Parkinson's disease (PD) is a movement disorder characterized by the pathological beta band (15-30 Hz) neural oscillations within the basal ganglia (BG). It is shown that the suppression of abnormal beta oscillations is correlated with the improvement of PD motor symptoms, which is a goal of standard therapies including deep brain stimulation (DBS). To overcome the stimulation-induced side effects and inefficiencies of conventional DBS (cDBS) and to reduce the administered stimulation current, closed-loop adaptive DBS (aDBS) techniques were developed. In this method, the frequency and/or amplitude of stimulation are modulated based on various disease biomarkers.

Methods

Here, by computational modeling of a cortico-BG-thalamic network in normal and PD conditions, we show that closed-loop aDBS of the subthalamic nucleus (STN) with amplitude modulation leads to a more effective suppression of pathological beta oscillations within the parkinsonian BG.

Results

Our results show that beta band neural oscillations are restored to their normal range and the reliability of the response of the thalamic neurons to motor cortex commands is retained due to aDBS with amplitude modulation. Furthermore, notably less stimulation current is administered during aDBS compared with cDBS due to a closed-loop control of stimulation amplitude based on the STN local field potential (LFP) beta activity.

Discussion

Efficient models of closed-loop stimulation may contribute to the clinical development of optimized aDBS techniques designed to reduce potential stimulation-induced side effects of cDBS in PD patients while leading to a better therapeutic outcome.

Remembrance of things perceived: Adding thalamocortical function to artificial neural networks

Recent research has illuminated the complexity and importance of the thalamocortical system but it has been difficult to identify what computational functions it performs. Meanwhile, deep-learning artificial neural networks (ANNs) based on bio-inspired models of purely cortical circuits have achieved surprising success solving sophisticated cognitive problems associated historically with human intelligence. Nevertheless, the limitations and shortcomings of artificial intelligence (AI) based on such ANNs are becoming increasingly clear. This review considers how the addition of thalamocortical connectivity and its putative functions related to cortical attention might address some of those shortcomings. Such bio-inspired models are now providing both testable theories of biological cognition and improved AI technology, much of which is happening outside the usual academic venues.

Intensity-Varied Closed-Loop Noise Stimulation for Oscillation Suppression in the Parkinsonian State

This work explores the effectiveness of the intensity-varied closed-loop noise stimulation on the oscillation suppression in the Parkinsonian state. Deep brain stimulation (DBS) is the standard therapy for Parkinson's disease (PD), but its effects need to be improved. The noise stimulation has compelling results in alleviating the PD state. However, in the open-loop control scheme, the noise stimulation parameters cannot be self-adjusted to adapt to the amplitude of the synchronized neuronal activities in real time. Thus, based on the delayed-feedback control algorithm, an intensity-varied closed-loop noise stimulation strategy is proposed. Based on a computational model of the basal ganglia (BG) that can present the intrinsic properties of the BG neurons and their interactions with the thalamic neurons, the proposed stimulation strategy is tested. Simulation results show that the noise stimulation suppresses the pathological beta (12-35 Hz) oscillations without any new rhythms in other bands compared with traditional high-frequency DBS. The intensity-varied closed-loop noise stimulation has a more profound role in removing the pathological beta oscillations and improving the thalamic reliability than open-loop noise stimulation, especially for different PD states. And the closed-loop noise stimulation enlarges the parameter space of the delayed-feedback control algorithm due to the randomness of noise signals. We also provide a theoretical analysis of the effective parameter domain of the delayed-feedback control algorithm by simplifying the BG model to an oscillator model. This exploration may guide a new approach to treating PD by optimizing the noise-induced improvement of the BG dysfunction.

Model-based quantitative optimization of deep brain stimulation and prediction of parkinson's states

Although the exact etiology of Parkinson's disease (PD) is still unknown, there are a variety of treatments available to alleviate its symptoms according to the development stage of PD. Deep brain stimulation (DBS), the most common surgical treatment for advanced PD, accurately locates and implants stimulating electrodes at specific targets in the brain to deliver high-frequency electrical stimulation that alters the excitability of the corresponding nuclei. However, for different patients and stages of PD development, there exists a choice of the optimal DBS protocol. In this paper, we propose a quantitative method (multi-dimensional feature indexes) to determine the stimulation pattern, stimulation parameters, and target of DBS from the perspective of the network model. On the other hand, based on this method, the development of PD can be predicted so that timely treatment can be given to patients. Simulation results show that, first, different network states can be distinguished by extracting features of the firing activity of neuronal populations within the basal ganglia network system. Secondly, the optimal DBS treatment can be selected by comparing the feature indexes vectors of the pre- and post-state of the network after the action of different modes of DBS. Lastly, the evolution of the network state from normal to pathological is simulated. The critical point of network state transitions is determined. These results provide a quantitative and qualitative method for determining the optimal regimen for DBS for PD, which is helpful for clinical practice.

Excitatory deep brain stimulation quenches beta oscillations arising in a computational model of the subthalamo-pallidal loop

Parkinson’s disease (PD) is associated with abnormal \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β band oscillations (13–30 Hz) in the cortico-basal ganglia circuits. Abnormally increased striato-pallidal inhibition and strengthening the synaptic coupling between subthalamic nucleus (STN) and globus pallidus externa (GPe), due to the loss of dopamine, are considered as the potential sources of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β oscillations in the basal ganglia. Deep brain stimulation (DBS) of the basal ganglia subregions is known as a way to reduce the pathological \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β oscillations and motor deficits related to PD. Despite the success of the DBS, its underlying mechanism is poorly understood and, there is controversy about the inhibitory or excitatory role of the DBS in the literature. Here, we utilized a computational network model of basal ganglia which consists of STN, GPe, globus pallidus interna, and thalamic neuronal population. This model can reproduce healthy and pathological \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β oscillations similar to what has been observed in experimental studies. Using this model, we investigated the effect of DBS to understand whether its effect is excitatory or inhibitory. Our results show that the excitatory DBS is able to quench the pathological synchrony and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta$$\end{document}β oscillations, while, applying inhibitory DBS failed to quench the PD signs. In light of simulation results, we conclude that the effect of the DBS on its target is excitatory.

Deep brain stimulation for movement disorder treatment: exploring frequency-dependent efficacy in a computational network model

A large-scale computational model of the basal ganglia network and thalamus is proposed to describe movement disorders and treatment effects of deep brain stimulation (DBS). The model of this complex network considers three areas of the basal ganglia region: the subthalamic nucleus (STN) as target area of DBS, the globus pallidus, both pars externa and pars interna (GPe-GPi), and the thalamus. Parkinsonian conditions are simulated by assuming reduced dopaminergic input and corresponding pronounced inhibitory or disinhibited projections to GPe and GPi. Macroscopic quantities are derived which correlate closely to thalamic responses and hence motor programme fidelity. It can be demonstrated that depending on different levels of striatal projections to the GPe and GPi, the dynamics of these macroscopic quantities (synchronisation index, mean synaptic activity and response efficacy) switch from normal to Parkinsonian conditions. Simulating DBS of the STN affects the dynamics of the entire network, increasing the thalamic activity to levels close to normal, while differing from both normal and Parkinsonian dynamics. Using the mentioned macroscopic quantities, the model proposes optimal DBS frequency ranges above 130 Hz.

Improved alpha-beta power reduction via combined electrical and ultrasonic stimulation in a parkinsonian cortex-basal ganglia-thalamus computational model

Objective. To investigate computationally the interaction of combined electrical and ultrasonic modulation of isolated neurons and of the parkinsonian cortex-basal ganglia-thalamus loop.Approach. Continuous-wave or pulsed electrical and ultrasonic neuromodulation is applied to isolated Otsuka plateau-potential generating subthalamic nucleus (STN) and Pospischil regular, fast and low-threshold spiking cortical cells in a temporally alternating or simultaneous manner. Similar combinations of electrical/ultrasonic waveforms are applied to a parkinsonian biophysical cortex-basal ganglia-thalamus neuronal network. Ultrasound-neuron interaction is modelled respectively for isolated neurons and the neuronal network with the NICE and SONIC implementations of the bilayer sonophore underlying mechanism. Reduction inα-βspectral energy is used as a proxy to express improvement in Parkinson's disease by insonication and electrostimulation.Main results. Simultaneous electro-acoustic stimulation achieves a given level of neuronal activity at lower intensities compared to the separate stimulation modalities. Conversely, temporally alternating stimulation with50 Hzelectrical and ultrasound pulses is capable of eliciting100 HzSTN firing rates. Furthermore, combination of ultrasound with hyperpolarizing currents can alter cortical cell relative spiking regimes. In the parkinsonian neuronal network, continuous-wave and pulsed ultrasound reduce pathological oscillations by different mechanisms. High-frequency pulsed separated electrical and ultrasonic deep brain stimulation (DBS) reduce pathologicalα-βpower by entraining STN-neurons. In contrast, continuous-wave ultrasound reduces pathological oscillations by silencing the STN. Compared to the separated stimulation modalities, temporally simultaneous or alternating electro-acoustic stimulation can achieve higher reductions inα-βpower for the same safety contraints on electrical/ultrasonic intensity.Significance. Focused ultrasound has the potential of becoming a non-invasive alternative of conventional DBS for the treatment of Parkinson's disease. Here, we elaborate on proposed benefits of combined electro-acoustic stimulation in terms of improved dynamic range, efficiency, spatial resolution, and neuronal selectivity.

Robust multi-input multi-output adaptive fuzzy terminal sliding mode control of deep brain stimulation in Parkinson's disease: a simulation study

Deep brain stimulation (DBS) has become an effective therapeutic solution for Parkinson’s disease (PD). Adaptive closed-loop DBS can be used to minimize stimulation-induced side effects by automatically determining the stimulation parameters based on the PD dynamics. In this paper, by modeling the interaction between the neurons in populations of the thalamic, the network-level modulation of thalamic is represented in a standard canonical form as a multi-input multi-output (MIMO) nonlinear first-order system with uncertainty and external disturbances. A class of fast and robust MIMO adaptive fuzzy terminal sliding mode control (AFTSMC) has been presented for control of membrane potential of thalamic neuron populations through continuous adaptive DBS current applied to the thalamus. A fuzzy logic system (FLS) is used to estimate the unknown nonlinear dynamics of the model, and the weights of FLS are adjusted online to guarantee the convergence of FLS parameters to optimal values. The simulation results show that the proposed AFTSMC not only significantly produces lower tracking errors in comparison with the classical adaptive fuzzy sliding mode control (AFSMC), but also makes more robust and reliable outputs. The results suggest that the proposed AFTSMC provides a more robust and smooth control input which is highly desirable for hardware design and implementation.

The Influence of Transcranial Magnetoacoustic Stimulation Parameters on the Basal Ganglia-Thalamus Neural Network in Parkinson's Disease

Objective: Parkinson's disease (PD) is a degenerative disease of the nervous system that frequently occurs in the aged. Transcranial magnetoacoustic stimulation (TMAS) is a neuronal adjustment method that combines sound fields and magnetic fields. It has the characteristics of high spatial resolution and noninvasive deep brain focusing. Methods: This paper constructed a simulation model of TMAS based on volunteer's skull computer tomography, phased controlled transducer and permanent magnet. It simulates a transcranial focused sound pressure field with the Westervelt equation and builds a basal ganglia and thalamus neural network model in the PD state based on the Hodgkin-Huxley model. Results: A biased sinusoidal pulsed ultrasonic TMAS induced current with 0.3 T static magnetic field induction and 0.2 W⋅cm-2 sound intensity can effectively modulate PD states with RI ≥ 0.633. The magnitude of magnetic induction strength was changed to 0.2 and 0.4 T. The induced current was the same when the sound intensity was 0.4 and 0.1 W⋅cm-2. And the sound pressure level is in the range of -1 dB (the induced current difference is less than or equal to 0.019 μA⋅cm-2). TMAS with a duty cycle of approximately 50% can effectively modulates the error firings in the PD neural network with a relay reliability not less than 0.633. Conclusion: TMAS can modulates the state of PD.

Subspace-based predictive control of Parkinson's disease: A model-based study

Deep brain stimulation (DBS) of the Basal Ganglia (BG) is an effective treatment to suppress the symptoms of Parkinson's disease (PD). Using a closed-loop scheme in DBS can not only improve its therapeutic effects but it can also reduce its energy consumption and possible side effects. In this paper, a predictive closed loop control strategy is employed to suppress the PD in real-time. A linear multi-input multi-output (MIMO) state-delayed system is considered as a simplified model of the BG neuronal network relating the stimulation signals as inputs to the beta power of local field potentials as PD biomarkers. The effect of time delay in different areas of the BG is incorporated into this model and a real-time subspace-based identification is implemented to continuously model the state of the BG neuronal network and drive the predictive control strategy. Simulation results show that the proposed MIMO subspace based predictive controller can suppress PD symptoms more effectively and with less power consumption compared to the conventional open-loop DBS and a recently proposed single-input single-output closed loop controller.

An Embedded Multi-Core Real-Time Simulation Platform of Basal Ganglia for Deep Brain Stimulation

Closed-loop deep brain stimulation (DBS) paradigm is gaining tremendous favor due to its potential capability of further and more efficient improvements in neurological diseases. Preclinical validation of closed-loop controller is quite necessary in order to minimize injury risks of clinical trials to patients, which can greatly benefit from real-time computational models and thus potentially reduce research and development costs and time. Here we developed an embedded multi-core real-time simulation platform (EMC-RTP) for a biological-faithful computational network model of basal ganglia (BG). The single neuron model is implemented in a highly real-time manner using a reasonable simplification. A modular mapping architecture with hierarchical routing organization was constructed to mimic the pathological neural activities of BG observed in parkinsonian conditions. A closed-loop simulation testbed for DBS validation was then set up using a host computer as the DBS controller. The availability of EMC-RTP and the testbed system was validated by comparing the performance of open-loop and proportional-integral (PI) controllers. Our experimental results showed that the proposed EMC-RTP reproduces abnormal beta bursts of BG in parkinsonian conditions while meets requirements of both real-time and computational accuracy as well. Closed-loop DBS experiments using the EMC-RTP suggested that the platform could perform reasonable output under different kinds of DBS strategies, indicating the usability of the platform.

Interaction of Indirect and Hyperdirect Pathways on Synchrony and Tremor-Related Oscillation in the Basal Ganglia

Low-frequency oscillatory activity (3-9 Hz) and increased synchrony in the basal ganglia (BG) are recognized to be crucial for Parkinsonian tremor. However, the dynamical mechanism underlying the tremor-related oscillations still remains unknown. In this paper, the roles of the indirect and hyperdirect pathways on synchronization and tremor-related oscillations are considered based on a modified Hodgkin-Huxley model. Firstly, the effects of indirect and hyperdirect pathways are analysed individually, which show that increased striatal activity to the globus pallidus external (GPe) or strong cortical gamma input to the subthalamic nucleus (STN) is sufficient to promote synchrony and tremor-related oscillations in the BG network. Then, the mutual effects of both pathways are analysed by adjusting the related currents simultaneously. Our results suggest that synchrony and tremor-related oscillations would be strengthened if the current of these two paths are in relative imbalance. And the network tends to be less synchronized and less tremulous when the frequency of cortical input is in the theta band. These findings may provide novel treatments in the cortex and striatum to alleviate symptoms of tremor in Parkinson's disease.

Model-Based Closed-Loop Suppression of Parkinsonian Beta Band Oscillations Through Origin Analysis

Excessive beta band (13-30 Hz) oscillations have been observed in the basal ganglia (BG) of patients with Parkinson's disease (PD). Understanding the origin and transmission of beta band oscillations are important to improve treatments of PD, such as closed-loop deep brain stimulation (DBS). This paper proposed a model-based closed-loop GPi stimulation system to suppress pathological beta band oscillations of BG. The feedback nucleus was selected through the analysis of GPi oscillations variation when different synaptic currents were blocked, mainly projections from globus pallidus external (GPe), the subthalamic nucleus (STN) and striatum. Since simulation results proved the important role of synaptic current from GPe in shaping the excessive GPi beta band oscillations, the local field potential (LFP) of GPe was chosen as the feedback signal. That is to say, the feedback nucleus was selected based on the origin analysis of the pathological GPi beta band oscillation. The closed-loop algorithm was the multiplication of linear delayed feedback of the filtered GPe-LFP and modeled synaptic dynamics from GPe to GPi. Thus, the formed stimulation waveform was synaptic current like shape, which was proved to be more energy efficient than open-loop continuous DBS in suppressing GPi beta band oscillation. With the development of DBS devices, the efficiency of this closed-loop stimulation could be testified in animal model and clinical.

A review of computational modeling and deep brain stimulation: applications to Parkinson's disease

Biophysical computational models are complementary to experiments and theories, providing powerful tools for the study of neurological diseases. The focus of this review is the dynamic modeling and control strategies of Parkinson's disease (PD). In previous studies, the development of parkinsonian network dynamics modeling has made great progress. Modeling mainly focuses on the cortex-thalamus-basal ganglia (CTBG) circuit and its sub-circuits, which helps to explore the dynamic behavior of the parkinsonian network, such as synchronization. Deep brain stimulation (DBS) is an effective strategy for the treatment of PD. At present, many studies are based on the side effects of the DBS. However, the translation from modeling results to clinical disease mitigation therapy still faces huge challenges. Here, we introduce the progress of DBS improvement. Its specific purpose is to develop novel DBS treatment methods, optimize the treatment effect of DBS for each patient, and focus on the study in closed-loop DBS. Our goal is to review the inspiration and insights gained by combining the system theory with these computational models to analyze neurodynamics and optimize DBS treatment.

Application of Reinforcement Learning to Deep Brain Stimulation in a Computational Model of Parkinson's Disease

Deep brain stimulation (DBS) has been proven to be an effective treatment to deal with the symptoms of Parkinson's disease (PD). Currently, the DBS is in an open-loop pattern with which the stimulation parameters remain constant regardless of fluctuations in the disease state, and adjustments of parameters rely mostly on trial and error of experienced clinicians. This could bring adverse effects to patients due to possible overstimulation. Thus closed-loop DBS of which stimulation parameters are automatically adjusted based on variations in the ongoing neurophysiological signals is desired. In this paper, we present a closed-loop DBS method based on reinforcement learning (RL) to regulate stimulation parameters based on a computational model. The network model consists of interconnected biophysically-based spiking neurons, and the PD state is described as distorted relay reliability of thalamus (TH). Results show that the RL-based closed-loop control strategy can effectively restore the distorted relay reliability of the TH but with less DBS energy expenditure.

Model-based optimized phase-deviation deep brain stimulation for Parkinson 's disease

High-frequency deep brain stimulation (HF-DBS) of the subthalamic nucleus (STN), globus pallidus interna (GPi) and globus pallidus externa (GPe) are often considered as effective methods for the treatment of Parkinson's disease (PD). However, the stimulation of a single nucleus by HF-DBS can cause specific physical damage, produce side effects and usually consume more electrical energy. Therefore, we use a biophysically-based model of basal ganglia-thalamic circuits to explore more effective stimulation patterns to reduce adverse effects and save energy. In this paper, we computationally investigate the combined DBS of two nuclei with the phase deviation between two stimulation waveforms (CDBS). Three different stimulation combination strategies are proposed, i.e., STN and GPe CDBS (SED), STN and GPi CDBS (SID), as well as GPi and GPe CDBS (GGD). Resultantly, it is found that anti-phase CDBS is more effective in improving parkinsonian dynamical properties, including desynchronization of neurons and the recovery of the thalamus relay ability. Detailed simulation investigation shows that anti-phase SED and GGD are superior to SID. Besides, the energy consumption can be largely reduced by SED and GGD (72.5% and 65.5%), compared to HF-DBS. These results provide new insights into the optimal stimulation parameter and target choice of PD, which may be helpful for the clinical practice.

Noise-Induced Improvement of the Parkinsonian State: A Computational Study

The benefit of noise in improving the basal ganglia (BG) dysfunctions, especially Parkinsonian state, is explored in this paper. High frequency (≥ 100 Hz) deep brain stimulation (DBS), as a clinical effective stimulation method, has compelling and fantastic results in alleviating the motor symptoms of Parkinson's disease (PD). However, the mechanism of DBS is still unclear. And the selection of the DBS waveform parameters faces great challenges to further optimize the stimulation effects and to reduce its energy expenditure. Considering that the desynchronization of the BG neuronal activities is benefited from the forced high frequency regular spikes driven by standard high frequency DBS, we expect to explore a novel stimulation method that has capability of restoring the BG physiological firing patterns without introducing artificial high-frequency fires. In this paper, a colored noise stimulation is used as a neuromodulation method to disrupt the firing patterns of the pathological neuronal activities. A computational model of the BG that exhibits the intrinsic properties of the BG neurons and their interactions with the thalamic (Th) cells is employed. Based on the model, we investigate the effects of noise stimulation and explore the impacts of the noise stimulation parameters on both relay reliability of the Th neurons and energy expenditure of the stimulation. By comparison, it can be found that noise stimulation does not entrain the network to an artificial high-frequency firing state, but induces the pathological increased synchronous activities back to a normal physiological level. Moreover, besides the capability of restoring the neuronal state, the benefits of the noise also include its balanced waveform to avert potential tissue or electrode damage and its ability to reduce the energy expenditure to 50% less than that of the standard DBS, when the noise stimulation has low frequency (≤ 100 Hz) and appropriate intensity. Thus, the exploration of the optimal noise-induced improvement of the BG dysfunction is of great significance in treating symptoms of neurological disorders such as PD.

Real-Time Neuromorphic System for Large-Scale Conductance-Based Spiking Neural Networks

The investigation of the human intelligence, cognitive systems and functional complexity of human brain is significantly facilitated by high-performance computational platforms. In this paper, we present a real-time digital neuromorphic system for the simulation of large-scale conductance-based spiking neural networks (LaCSNN), which has the advantages of both high biological realism and large network scale. Using this system, a detailed large-scale cortico-basal ganglia-thalamocortical loop is simulated using a scalable 3-D network-on-chip (NoC) topology with six Altera Stratix III field-programmable gate arrays simulate 1 million neurons. Novel router architecture is presented to deal with the communication of multiple data flows in the multinuclei neural network, which has not been solved in previous NoC studies. At the single neuron level, cost-efficient conductance-based neuron models are proposed, resulting in the average utilization of 95% less memory resources and 100% less DSP resources for multiplier-less realization, which is the foundation of the large-scale realization. An analysis of the modified models is conducted, including investigation of bifurcation behaviors and ionic dynamics, demonstrating the required range of dynamics with a more reduced resource cost. The proposed LaCSNN system is shown to outperform the alternative state-of-the-art approaches previously used to implement the large-scale spiking neural network, and enables a broad range of potential applications due to its real-time computational power.

Robust desynchronization of Parkinson's disease pathological oscillations by frequency modulation of delayed feedback deep brain stimulation

The hyperkinetic symptoms of Parkinson's Disease (PD) are associated with the ensembles of interacting oscillators that cause excess or abnormal synchronous behavior within the Basal Ganglia (BG) circuitry. Delayed feedback stimulation is a closed loop technique shown to suppress this synchronous oscillatory activity. Deep Brain Stimulation (DBS) via delayed feedback is known to destabilize the complex intermittent synchronous states. Computational models of the BG network are often introduced to investigate the effect of delayed feedback high frequency stimulation on partially synchronized dynamics. In this study, we develop a reduced order model of four interacting nuclei of the BG as well as considering the Thalamo-Cortical local effects on the oscillatory dynamics. This model is able to capture the emergence of 34 Hz beta band oscillations seen in the Local Field Potential (LFP) recordings of the PD state. Train of high frequency pulses in a delayed feedback stimulation has shown deficiencies such as strengthening the synchronization in case of highly fluctuating neuronal activities, increasing the energy consumed as well as the incapability of activating all neurons in a large-scale network. To overcome these drawbacks, we propose a new feedback control variable based on the filtered and linearly delayed LFP recordings. The proposed control variable is then used to modulate the frequency of the stimulation signal rather than its amplitude. In strongly coupled networks, oscillations reoccur as soon as the amplitude of the stimulus signal declines. Therefore, we show that maintaining a fixed amplitude and modulating the frequency might ameliorate the desynchronization process, increase the battery lifespan and activate substantial regions of the administered DBS electrode. The charge balanced stimulus pulse itself is embedded with a delay period between its charges to grant robust desynchronization with lower amplitudes needed. The efficiency of the proposed Frequency Adjustment Stimulation (FAS) protocol in a delayed feedback method might contribute to further investigation of DBS modulations aspired to address a wide range of abnormal oscillatory behavior observed in neurological disorders.

Computational Modeling of Ultrasonic Subthalamic Nucleus Stimulation

OBJECTIVE

To explore the potential of ultrasonic modulation of plateau-potential generating subthalamic nucleus neurons (STN), by modeling their interaction with continuous and pulsed ultrasonic waves.

METHODS

A computational model for ultrasonic stimulation of the STN is created by combining the Otsuka-model with the bilayer sonophore model. The neuronal response to continuous and pulsed ultrasonic waves is computed in parallel for a range of frequencies, duty cycles, pulse repetition frequencies, and intensities.

RESULTS

Ultrasonic intensity in continuous-wave stimulation determines the firing pattern of the STN. Three observed spiking modes in order of increasing intensity are low frequency spiking, high frequency spiking with significant spike-frequency and spike-amplitude adaptation, and a silenced mode. Continuous-wave stimulation has little capability to manipulate the saturated spiking rate in the high frequency spiking mode. In contrast, STN firing rates induced by pulsed ultrasound insonication will saturate to the pulse repetition frequency with short latencies, for sufficiently large intensity and repetition frequency.

CONCLUSION

Computational results show that the activity of plateau-potential generating STN can be modulated by selection of the stimulus parameters. Low intensities result in repetitive firing, while higher intensities silence the STN. Pulsed ultrasonic stimulation results in a shorter saturation latency and is able to modulate spiking rates.

SIGNIFICANCE

Stimulation or suppresion of the STN is important in the treatment of Parkinson's disease, e.g., in deep brain stimulation. This explorative study on ultrasonic modulation of the STN, could be a step in the direction of minimally invasive alternatives to conventional deep brain stimulation.

Systems approaches to optimizing deep brain stimulation therapies in Parkinson's disease

Over the last 30 years, deep brain stimulation (DBS) has been used to treat chronic neurological diseases like dystonia, obsessive-compulsive disorders, essential tremor, Parkinson's disease, and more recently, dementias, depression, cognitive disorders, and epilepsy. Despite its wide use, DBS presents numerous challenges for both clinicians and engineers. One challenge is the design of novel, more efficient DBS therapies, which are hampered by the lack of complete understanding about the cellular mechanisms of therapeutic DBS. Another challenge is the existence of redundancy in clinical outcomes, that is, different DBS programs can result in similar clinical benefits but very little information (e.g., predictive models, longitudinal data, metrics, etc.) is available to select one program over another. Finally, there is high variability in patients' responses to DBS, which forces clinicians to carefully adjust the stimulation settings to each patient via lengthy programming sessions. Researchers in neural engineering and systems biology have been tackling these challenges over the past few years with the specific goal of developing novel DBS therapies, design methodologies, and computational tools that optimize the therapeutic effects of DBS in each patient. Furthermore, efforts are being made to automatically adapt the DBS treatment to the fluctuations of disease symptoms. A review of the quantitative approaches currently available for the treatment of Parkinson's disease is presented here with an emphasis on the contributions that systems theoretical approaches have provided to understand the global dynamics of complex neuronal circuits in the brain under DBS. This article is categorized under: Translational, Genomic, and Systems Medicine > Therapeutic Methods Analytical and Computational Methods > Computational Methods Analytical and Computational Methods > Dynamical Methods Physiology > Mammalian Physiology in Health and Disease.

Patterned low-frequency deep brain stimulation induces motor deficits and modulates cortex-basal ganglia neural activity in healthy rats

Deep brain stimulation (DBS) is an effective therapy for movement disorders, including Parkinson's disease (PD), although the mechanisms of action remain unclear. Abnormal oscillatory neural activity is correlated with motor symptoms, and pharmacological or DBS treatment that alleviates motor symptoms appears to suppress abnormal oscillations. However, whether such oscillatory activity is causal of motor deficits such as tremor remains unclear. Our goal was to generate abnormal oscillatory activity in the cortex-basal ganglia loop using patterned subthalamic nucleus DBS and to quantify motor behavior in awake healthy rats. Stimulation patterns were designed via model-based optimization to increase power in the low-frequency (7-11 Hz) band because these oscillations are associated with the emergence of motor symptoms in the 6-hydroxydopamine lesioned rat model of parkinsonism. We measured motor activity using a head-mounted accelerometer, as well as quantified neural activity in cortex and globus pallidus (GP), in response to 5 stimulation patterns that generated a range of 7- to 11-Hz spectral power. Stimulation patterns induced oscillatory activity in the low-frequency band in the cortex and GP and caused tremor, whereas control patterns and regular 50-Hz DBS did not generate any such effects. Neural and motor-evoked responses observed during stimulation were synchronous and time-locked to stimulation bursts within the patterns. These results identified elements of irregular patterns of stimulation that were correlated with tremor and tremor-related neural activity in the cortex and basal ganglia and may lead to the identification of the oscillatory activity and structures associated with the generation of tremor activity. NEW & NOTEWORTHY Subthalamic nucleus deep brain stimulation is a promising therapy for movement disorders such as Parkinson's disease. Several groups reported correlation between suppression of abnormal oscillatory activity in the cortex-basal ganglia and motor symptoms, but it remains unclear whether such oscillations play a causal role in the emergence of motor symptoms. We demonstrate generation of tremor and pathological oscillatory activity in otherwise healthy rats by stimulation with patterns that produced increases in low-frequency oscillatory activity.

Closed-loop low-frequency DBS restores thalamocortical relay fidelity in a computational model of the motor loop

Closed-loop modulation of deep brain stimulation (DBS) of the subthalamic nucleus (STN) in Parkinson's disease (PD) is investigated to automatically adjust the stimulation to the patients' conditions, optimize the clinical outcomes, and reduce the energy requirements. This study proposes a closed-loop control system for real-time adaptation of the STN-DBS amplitude based on the neural activity in the motor thalamus. Population-averaged post-stimulus time histograms are used to measure the average effects of STN-DBS on the thalamocortical neurons and a L₂-norm minimization problem is solved to design the control algorithm, while the frequency of stimulation is kept constant. Applied on a large-scale, biophysically-based, anatomically-compliant model of the cortico-basal ganglia-thalamo-cortical motor loop under PD conditions, our adaptive DBS significantly (P-value P<;0.05) improved the relay fidelity of the thalamocortical neurons and restored the average power of the thalamocortical spike trains in the band [3, 100] Hz, two indicators of restored thalamocortical activity. Furthermore, adaptive-DBS significantly decreased the energy requirements when compared with non-adaptive-DBS at the same frequency. Finally, 30- and 60-Hz-adaptive-DBS determined the maximal restoration of thalamocortical activity and outperformed high-frequency, non-adaptive-DBS. Overall, results suggest that a feedback-controlled, low-frequency DBS pattern may result in significant restoration of the thalamocortical encoding while lowering the energy requirements.

Modeling and Analysis of Beta Oscillations in the Basal Ganglia

Enhanced beta (12-30 Hz) oscillatory activity in the basal ganglia (BG) is a prominent feature of the Parkinsonian state in animal models and in patients with Parkinson's disease. Increased beta oscillations are associated with severe dopaminergic striatal depletion. However, the mechanisms underlying these pathological beta oscillations remain elusive. Inspired by the experimental observation that only subsets of neurons within each nucleus in the BG exhibit oscillatory activities, a computational model of the BG-thalamus neuronal network is proposed, which is characterized by subdivided nuclei within the BG. Using different currents externally applied to the neurons within a given nucleus, neurons behave according to one of the two subgroups, named "-N" and "-P," where "-N" and "-P" denote the normal and the Parkinsonian states, respectively. The ratio of "-P" to "-N" neurons indicates the degree of the Parkinsonian state. Simulation results show that if "-P" neurons have a high degree of connectivity in the subthalamic nucleus (STN), they will have a significant downstream effect on the generation of beta oscillations in the globus pallidus. Interestingly, however, the generation of beta oscillations in the STN is independent of the selection of the "-P" neurons in the external segment of the globus pallidus (GPe), despite the reciprocal structure between STN and GPe. This computational model may pave the way to revealing the mechanism of such pathological behaviors in a realistic way that can replicate experimental observations. The simulation results suggest that the STN is more suitable than GPe as a deep brain stimulation target.

Model-based deconstruction of cortical evoked potentials generated by subthalamic nucleus deep brain stimulation

Parkinson's disease is associated with altered neural activity in the motor cortex. Chronic high-frequency deep brain stimulation (DBS) of the subthalamic nucleus (STN) is effective in suppressing parkinsonian motor symptoms and modulates cortical activity. However, the anatomical pathways responsible for STN DBS-mediated cortical modulation remain unclear. Cortical evoked potentials (cEP) generated by STN DBS reflect the response of cortex to subcortical stimulation, and the goal of this study was to determine the neural origin of STN DBS-generated cEP using a two-step approach. First, we recorded cEP over ipsilateral primary motor cortex during different frequencies of STN DBS in awake healthy and unilateral 6-OHDA-lesioned parkinsonian rats. Second, we used a detailed, biophysically based model of the thalamocortical network to deconstruct the neural origin of the recorded cEP. The in vivo cEP included short (R1)-, intermediate (R2)-, and long-latency (R3) responses. Model-based cortical responses to simulated STN DBS matched remarkably well the in vivo responses. The short-latency response was generated by antidromic activation of layer 5 pyramidal neurons, whereas recurrent activation of layer 5 pyramidal neurons via excitatory axon collaterals reproduced the intermediate-latency response. The long-latency response was generated by polysynaptic activation of layer 2/3 pyramidal neurons via the cortico-thalamic-cortical pathway. Antidromic activation of the hyperdirect pathway and subsequent intracortical and cortico-thalamo-cortical synaptic interactions were sufficient to generate cortical potential evoked by STN DBS, and orthodromic activation through basal ganglia-thalamus-cortex pathways was not required. These results demonstrate the utility of cEP to determine the neural elements activated by STN DBS that might modulate cortical activity and contribute to the suppression of parkinsonian symptoms. NEW & NOTEWORTHY Subthalamic nucleus (STN) deep brain stimulation (DBS) is increasingly used to treat Parkinson's disease (PD). Cortical potentials evoked by STN DBS in patients with PD exhibit consistent short-latency (1-3 ms), intermediate-latency (5-15 ms), and long-latency (18-25 ms) responses. The short-latency response occurs as a result of antidromic activation of the hyperdirect pathway comprising corticosubthalamic axons. However, the neural origins of intermediate- and long-latency responses remain elusive, and the dominant view is that these are produced through the orthodromic pathway (basal ganglia-thalamus-cortex). By combining in vivo electrophysiology with computational modeling, we demonstrate that antidromic activation of the cortico-thalamic-cortical pathway is sufficient to generate the intermediate- and long-latency cortical responses to STN DBS.

Nonlinear predictive control for adaptive adjustments of deep brain stimulation parameters in basal ganglia-thalamic network

The efficacy of deep brain stimulation (DBS) for Parkinson's disease (PD) depends in part on the post-operative programming of stimulation parameters. Closed-loop stimulation is one method to realize the frequent adjustment of stimulation parameters. This paper introduced the nonlinear predictive control method into the online adjustment of DBS amplitude and frequency. This approach was tested in a computational model of basal ganglia-thalamic network. The autoregressive Volterra model was used to identify the process model based on physiological data. Simulation results illustrated the efficiency of closed-loop stimulation methods (amplitude adjustment and frequency adjustment) in improving the relay reliability of thalamic neurons compared with the PD state. Besides, compared with the 130Hz constant DBS the closed-loop stimulation methods can significantly reduce the energy consumption. Through the analysis of inter-spike-intervals (ISIs) distribution of basal ganglia neurons, the evoked network activity by the closed-loop frequency adjustment stimulation was closer to the normal state.

The Effects of Medium Spiny Neuron Morphologcial Changes on Basal Ganglia Network under External Electric Field: A Computational Modeling Study

The damage of dopaminergic neurons that innervate the striatum has been considered to be the proximate cause of Parkinson's disease (PD). In the dopamine-denervated state, the loss of dendritic spines and the decrease of dendritic length may prevent medium spiny neuron (MSN) from receiving too much excitatory stimuli from the cortex, thereby reducing the symptom of Parkinson's disease. However, the reduction in dendritic spine density obtained by different experiments is significantly different. We developed a biological-based network computational model to quantify the effect of dendritic spine loss and dendrites tree degeneration on basal ganglia (BG) signal regulation. Through the introduction of error index (EI), which was used to measure the attenuation of the signal, we explored the amount of dendritic spine loss and dendritic trees degradation required to restore the normal regulatory function of the network, and found that there were two ranges of dendritic spine loss that could reduce EI to normal levels in the case of dopamine at a certain level, this was also true for dendritic trees. However, although these effects were the same, the mechanisms of these two cases were significant difference. Using the method of phase diagram analysis, we gained insight into the mechanism of signal degradation. Furthermore, we explored the role of cortex in MSN morphology changes dopamine depletion-induced and found that proper adjustments to cortical activity do stop the loss in dendritic spines induced by dopamine depleted. These results suggested that modifying cortical drive onto MSN might provide a new idea on clinical therapeutic strategies for Parkinson's disease.

Computational Stimulation of the Basal Ganglia Neurons with Cost Effective Delayed Gaussian Waveforms

Deep brain stimulation (DBS) has compelling results in the desynchronization of the basal ganglia neuronal activities and thus, is used in treating the motor symptoms of Parkinson's disease (PD). Accurate definition of DBS waveform parameters could avert tissue or electrode damage, increase the neuronal activity and reduce energy cost which will prolong the battery life, hence avoiding device replacement surgeries. This study considers the use of a charge balanced Gaussian waveform pattern as a method to disrupt the firing patterns of neuronal cell activity. A computational model was created to simulate ganglia cells and their interactions with thalamic neurons. From the model, we investigated the effects of modified DBS pulse shapes and proposed a delay period between the cathodic and anodic parts of the charge balanced Gaussian waveform to desynchronize the firing patterns of the GPe and GPi cells. The results of the proposed Gaussian waveform with delay outperformed that of rectangular DBS waveforms used in in-vivo experiments. The Gaussian Delay Gaussian (GDG) waveforms achieved lower number of misses in eliciting action potential while having a lower amplitude and shorter length of delay compared to numerous different pulse shapes. The amount of energy consumed in the basal ganglia network due to GDG waveforms was dropped by 22% in comparison with charge balanced Gaussian waveforms without any delay between the cathodic and anodic parts and was also 60% lower than a rectangular charged balanced pulse with a delay between the cathodic and anodic parts of the waveform. Furthermore, by defining a Synchronization Level metric, we observed that the GDG waveform was able to reduce the synchronization of GPi neurons more effectively than any other waveform. The promising results of GDG waveforms in terms of eliciting action potential, desynchronization of the basal ganglia neurons and reduction of energy consumption can potentially enhance the performance of DBS devices.

The role of the basal ganglia in the control of seizure

Epilepsy is a network disorder and each type of seizure involves distinct cortical and subcortical network, differently implicated in the control and propagation of the ictal activity. The role of the basal ganglia has been revealed in several cases of focal and generalized seizures. Here, we review the data that show the implication of the basal ganglia in absence, temporal lobe, and neocortical seizures in animal models (rodent, cat, and non-human primate) and in human. Based on these results and the advancement of deep brain stimulation for Parkinson's disease, basal ganglia neuromodulation has been tested with some success that can be equally seen as promising or disappointing. The effect of deep brain stimulation can be considered promising with a 76% in seizure reduction in temporal lobe epilepsy patients, but also disappointing, since only few patients have become seizure free and the antiepileptic effects have been highly variable among patients. This variability could probably be explained by the heterogeneity among the patients included in these clinical studies. To illustrate the importance of specific network identification, electrophysiological activity of the putamen and caudate nucleus has been recorded during penicillin-induced pre-frontal and motor seizures in one monkey. While an increase of the firing rate was found in putamen and caudate nucleus during pre-frontal seizures, only the activity of the putamen cells was increased during motor seizures. These preliminary results demonstrate the implication of the basal ganglia in two types of neocortical seizures and the necessity of studying the network to identify the important nodes implicated in the propagation and control of each type of seizure.

Eliminating Absence Seizures through the Deep Brain Stimulation to Thalamus Reticular Nucleus

Deep brain stimulation (DBS) can play a crucial role in the modulation of absence seizures, yet relevant biophysical mechanisms are not completely established. In this paper, on the basis of a biophysical mean-field model, we investigate a typical absence epilepsy activity by introducing slow kinetics of GABA_B receptors on thalamus reticular nucleus (TRN). We find that the region of spike and slow-wave discharges (SWDs) can be reduced greatly when we add the DBS to TRN. Furthermore, we systematically explore how the corresponding stimulation parameters including frequency, amplitude and positive input duration suppress the SWDs under certain conditions. It is shown that the SWDs can be controlled as key stimulation parameters are suitably chosen. The results in this paper can be helpful for researchers to understand the thalamus stimulation in treating epilepsy patients, and provide theoretical basis for future experimental and clinical studies.

Closed-Loop Modulation of the Pathological Disorders of the Basal Ganglia Network

A generalized predictive closed-loop control strategy to improve the basal ganglia activity patterns in Parkinson's disease (PD) is explored in this paper. Based on system identification, an input-output model is established to reveal the relationship between external stimulation and neuronal responses. The model contributes to the implementation of the generalized predictive control (GPC) algorithm that generates the optimal stimulation waveform to modulate the activities of neuronal nuclei. By analyzing the roles of two critical control parameters within the GPC law, optimal closed-loop control that has the capability of restoring the normal relay reliability of the thalamus with the least stimulation energy expenditure can be achieved. In comparison with open-loop deep brain stimulation and traditional static control schemes, the generalized predictive closed-loop control strategy can optimize the stimulation waveform without requiring any particular knowledge of the physiological properties of the system. This type of closed-loop control strategy generates an adaptive stimulation waveform with low energy expenditure with the potential to improve the treatments for PD.

Efficient implementation of a real-time estimation system for thalamocortical hidden Parkinsonian properties

Real-time estimation of dynamical characteristics of thalamocortical cells, such as dynamics of ion channels and membrane potentials, is useful and essential in the study of the thalamus in Parkinsonian state. However, measuring the dynamical properties of ion channels is extremely challenging experimentally and even impossible in clinical applications. This paper presents and evaluates a real-time estimation system for thalamocortical hidden properties. For the sake of efficiency, we use a field programmable gate array for strictly hardware-based computation and algorithm optimization. In the proposed system, the FPGA-based unscented Kalman filter is implemented into a conductance-based TC neuron model. Since the complexity of TC neuron model restrains its hardware implementation in parallel structure, a cost efficient model is proposed to reduce the resource cost while retaining the relevant ionic dynamics. Experimental results demonstrate the real-time capability to estimate thalamocortical hidden properties with high precision under both normal and Parkinsonian states. While it is applied to estimate the hidden properties of the thalamus and explore the mechanism of the Parkinsonian state, the proposed method can be useful in the dynamic clamp technique of the electrophysiological experiments, the neural control engineering and brain-machine interface studies.

A Network Model of Local Field Potential Activity in Essential Tremor and the Impact of Deep Brain Stimulation

Essential tremor (ET), a movement disorder characterised by an uncontrollable shaking of the affected body part, is often professed to be the most common movement disorder, affecting up to one percent of adults over 40 years of age. The precise cause of ET is unknown, however pathological oscillations of a network of a number of brain regions are implicated in leading to the disorder. Deep brain stimulation (DBS) is a clinical therapy used to alleviate the symptoms of a number of movement disorders. DBS involves the surgical implantation of electrodes into specific nuclei in the brain. For ET the targeted region is the ventralis intermedius (Vim) nucleus of the thalamus. Though DBS is effective for treating ET, the mechanism through which the therapeutic effect is obtained is not understood. To elucidate the mechanism underlying the pathological network activity and the effect of DBS on such activity, we take a computational modelling approach combined with electrophysiological data. The pathological brain activity was recorded intra-operatively via implanted DBS electrodes, whilst simultaneously recording muscle activity of the affected limbs. We modelled the network hypothesised to underlie ET using the Wilson-Cowan approach. The modelled network exhibited oscillatory behaviour within the tremor frequency range, as did our electrophysiological data. By applying a DBS-like input we suppressed these oscillations. This study shows that the dynamics of the ET network support oscillations at the tremor frequency and the application of a DBS-like input disrupts this activity, which could be one mechanism underlying the therapeutic benefit.

Essential tremor (ET) is acknowledged to be the most common movement disorder affecting 1% of the population. Although the underlying mechanisms remain elusive, the thalamus, cortex and cerebellum are implicated in the underlying pathology. More recently, it has been shown that ET can be successfully treated by deep brain stimulation (DBS). This clinical treatment involves the surgical implantation of electrodes into the brain, through which current is applied. However, the mechanisms of how DBS achieves clinical benefit continue to be debated. A key question is whether ET can be modeled as a pathological network behavior as has been suggested previously. If so, we can then ask how DBS would modulate this brain activity. Our study combines: (i) simultaneous electrophysiological recordings from the brain and muscle; (ii) computational modelling; (iii) mathematical analysis. We found that the network supports oscillations in the tremor range, and the application of high frequency DBS switches this to low amplitude, high-frequency activity. We propose that our model can be used to predict DBS parameter settings that suppress pathological network activity and consequently tremor. In summary, we provide the first population level model of essential tremor including the effect of DBS on network behaviour.

Optimized temporal pattern of brain stimulation designed by computational evolution

Brain stimulation is a promising therapy for several neurological disorders, including Parkinson's disease. Stimulation parameters are selected empirically and are limited to the frequency and intensity of stimulation. We varied the temporal pattern of deep brain stimulation to ameliorate symptoms in a parkinsonian animal model and in humans with Parkinson's disease. We used model-based computational evolution to optimize the stimulation pattern. The optimized pattern produced symptom relief comparable to that from standard high-frequency stimulation (a constant rate of 130 or 185 Hz) and outperformed frequency-matched standard stimulation in a parkinsonian rat model and in patients. Both optimized and standard high-frequency stimulation suppressed abnormal oscillatory activity in the basal ganglia of rats and humans. The results illustrate the utility of model-based computational evolution of temporal patterns to increase the efficiency of brain stimulation in treating Parkinson's disease and thereby reduce the energy required for successful treatment below that of current brain stimulation paradigms.

A biophysical model of the cortex-basal ganglia-thalamus network in the 6-OHDA lesioned rat model of Parkinson's disease

Electrical stimulation of sub-cortical brain regions (the basal ganglia), known as deep brain stimulation (DBS), is an effective treatment for Parkinson's disease (PD). Chronic high frequency (HF) DBS in the subthalamic nucleus (STN) or globus pallidus interna (GPi) reduces motor symptoms including bradykinesia and tremor in patients with PD, but the therapeutic mechanisms of DBS are not fully understood. We developed a biophysical network model comprising of the closed loop cortical-basal ganglia-thalamus circuit representing the healthy and parkinsonian rat brain. The network properties of the model were validated by comparing responses evoked in basal ganglia (BG) nuclei by cortical (CTX) stimulation to published experimental results. A key emergent property of the model was generation of low-frequency network oscillations. Consistent with their putative pathological role, low-frequency oscillations in model BG neurons were exaggerated in the parkinsonian state compared to the healthy condition. We used the model to quantify the effectiveness of STN DBS at different frequencies in suppressing low-frequency oscillatory activity in GPi. Frequencies less than 40 Hz were ineffective, low-frequency oscillatory power decreased gradually for frequencies between 50 Hz and 130 Hz, and saturated at frequencies higher than 150 Hz. HF STN DBS suppressed pathological oscillations in GPe/GPi both by exciting and inhibiting the firing in GPe/GPi neurons, and the number of GPe/GPi neurons influenced was greater for HF stimulation than low-frequency stimulation. Similar to the frequency dependent suppression of pathological oscillations, STN DBS also normalized the abnormal GPi spiking activity evoked by CTX stimulation in a frequency dependent fashion with HF being the most effective. Therefore, therapeutic HF STN DBS effectively suppresses pathological activity by influencing the activity of a greater proportion of neurons in the output nucleus of the BG.

Failure to suppress low-frequency neuronal oscillatory activity underlies the reduced effectiveness of random patterns of deep brain stimulation

Subthalamic nucleus (STN) deep brain stimulation (DBS) is an established treatment for the motor symptoms of Parkinson's disease (PD). However, the mechanisms of action of DBS are unknown. Random temporal patterns of DBS are less effective than regular DBS, but the neuronal basis for this dependence on temporal pattern of stimulation is unclear. Using a rat model of PD, we quantified the changes in behavior and single-unit activity in globus pallidus externa and substantia nigra pars reticulata during high-frequency STN DBS with different degrees of irregularity. Although all stimulus trains had the same average rate, 130-Hz regular DBS more effectively reversed motor symptoms, including circling and akinesia, than 130-Hz irregular DBS. A mixture of excitatory and inhibitory neuronal responses was present during all stimulation patterns, and mean firing rate did not change during DBS. Low-frequency (7-10 Hz) oscillations of single-unit firing times present in hemiparkinsonian rats were suppressed by regular DBS, and neuronal firing patterns were entrained to 130 Hz. Irregular patterns of DBS less effectively suppressed 7- to 10-Hz oscillations and did not regularize firing patterns. Random DBS resulted in a larger proportion of neuron pairs with increased coherence at 7-10 Hz compared with regular 130-Hz DBS, which suggested that long pauses (interpulse interval >50 ms) during random DBS facilitated abnormal low-frequency oscillations in the basal ganglia. These results suggest that the efficacy of high-frequency DBS stems from its ability to regularize patterns of neuronal firing and thereby suppress abnormal oscillatory neural activity within the basal ganglia.