Addition of deep brain stimulation signal to a local field potential driven Izhikevich model masks the pathological firing pattern of an STN neuron

β Oscillations of Dorsal STN as a Potential Biomarker in Parkinson's Disease Motor Subtypes: An Exploratory Study

Parkinson's disease (PD) can be divided into postural instability and difficult gait (PIGD) and tremor dominance (TD) subtypes. However, potential neural markers located in the dorsal ventral side of the subthalamic nucleus (STN) for delineating the two subtypes of PIGD and TD have not been demonstrated. Therefore, this study aimed to investigate the spectral characteristics of PD on the dorsal ventral side. The differences in the β oscillation spectrum of the spike signal on the dorsal and ventral sides of the STN during deep brain stimulation (DBS) were investigated in 23 patients with PD, and coherence analysis was performed for both subtypes. Finally, each feature was associated with the Unified Parkinson's Disease Rating Scale (UPDRS). The β power spectral density (PSD) in the dorsal STN was found to be the best predictor of the PD subtype, with 82.6% accuracy. The PSD of dorsal STN β oscillations was greater in the PIGD group than in the TD group (22.17% vs. 18.22%; p < 0.001). Compared with the PIGD group, the TD group showed greater consistency in the β and γ bands. In conclusion, dorsal STN β oscillations could be used as a biomarker to classify PIGD and TD subtypes, guide STN-DBS treatment, and relate to some motor symptoms.

A Multi-Scale Computational Model of Levodopa-Induced Toxicity in Parkinson's Disease

Parkinson's disease (PD) is caused by the progressive loss of dopaminergic cells in substantia nigra pars compacta (SNc). The root cause of this cell loss in PD is still not decisively elucidated. A recent line of thinking has traced the cause of PD neurodegeneration to metabolic deficiency. Levodopa (L-DOPA), a precursor of dopamine, used as a symptom-relieving treatment for PD, leads to positive and negative outcomes. Several researchers inferred that L-DOPA might be harmful to SNc cells due to oxidative stress. The role of L-DOPA in the course of the PD pathogenesis is still debatable. We hypothesize that energy deficiency can lead to L-DOPA-induced toxicity in two ways: by promoting dopamine-induced oxidative stress and by exacerbating excitotoxicity in SNc. We present a systems-level computational model of SNc-striatum, which will help us understand the mechanism behind neurodegeneration postulated above and provide insights into developing disease-modifying therapeutics. It was observed that SNc terminals are more vulnerable to energy deficiency than SNc somas. During L-DOPA therapy, it was observed that higher L-DOPA dosage results in increased loss of terminals in SNc. It was also observed that co-administration of L-DOPA and glutathione (antioxidant) evades L-DOPA-induced toxicity in SNc neurons. Our proposed model of the SNc-striatum system is the first of its kind, where SNc neurons were modeled at a biophysical level, and striatal neurons were modeled at a spiking level. We show that our proposed model was able to capture L-DOPA-induced toxicity in SNc, caused by energy deficiency.

Automatic identification of sweet spots from MERs for electrodes implantation in STN-DBS

PURPOSE

Microelectrode recordings (MERs) are a significant clinical indicator for sweet spots identification of implanted electrodes during deep brain stimulation of the subthalamic nucleus (STN) surgery. As 1D MERs signals have the unboundedness, large-range, large-amount and time-dependent characteristics, the purpose of this study is to propose an automatic and precise identification method of sweet spots from MERs, reducing the time-consuming and labor-intensive human annotations.

METHODS

We propose an automatic identification method of sweet spots from MERs for electrodes implantation in STN-DBS. To better imitate the surgeons' observation and obtain more intuitive contextual information, we first employ the 2D Gramian angular summation field (GASF) images generated from MERs data to perform the sweet spots determination for electrodes implantation. Then, we introduce the convolutional block attention module into convolutional neural network (CNN) to identify the 2D GASF images of sweet spots for electrodes implantation.

RESULTS

Experimental results illustrate that the identification result of our method is consistent with the result of doctor's decision, while our method can achieve the accuracy and precision of 96.72% and 98.97%, respectively, which outperforms state-of-the-art for intraoperative sweet spots determination.

CONCLUSIONS

The proposed method is the first time to automatically and accurately identify sweet spots from MERs for electrodes implantation by the combination an advanced time series-to-image encoding way with CBAM-enhanced networks model. Our method can assist neurosurgeons in automatically detecting the most likely locations of sweet spots for electrodes implantation, which can provide an important indicator for target selection while it reduces the localization error of the target during STN-DBS surgery.

A Multi-Scale Computational Model of Excitotoxic Loss of Dopaminergic Cells in Parkinson's Disease

Parkinson's disease (PD) is a neurodegenerative disorder caused by loss of dopaminergic neurons in substantia nigra pars compacta (SNc). Although the exact cause of cell death is not clear, the hypothesis that metabolic deficiency is a key factor has been gaining attention in recent years. In the present study, we investigated this hypothesis using a multi-scale computational model of the subsystem of the basal ganglia comprising the subthalamic nucleus (STN), globus pallidus externa (GPe), and SNc. The proposed model is a multiscale model in that interaction among the three nuclei are simulated using more abstract Izhikevich neuron models, while the molecular pathways involved in cell death of SNc neurons are simulated in terms of detailed chemical kinetics. Simulation results obtained from the proposed model showed that energy deficiencies occurring at cellular and network levels could precipitate the excitotoxic loss of SNc neurons in PD. At the subcellular level, the models show how calcium elevation leads to apoptosis of SNc neurons. The therapeutic effects of several neuroprotective interventions are also simulated in the model. From neuroprotective studies, it was clear that glutamate inhibition and apoptotic signal blocker therapies were able to halt the progression of SNc cell loss when compared to other therapeutic interventions, which only slowed down the progression of SNc cell loss.

A Computational Model of Loss of Dopaminergic Cells in Parkinson's Disease Due to Glutamate-Induced Excitotoxicity

Parkinson's disease (PD) is a neurodegenerative disease associated with progressive and inexorable loss of dopaminergic cells in Substantia Nigra pars compacta (SNc). Although many mechanisms have been suggested, a decisive root cause of this cell loss is unknown. A couple of the proposed mechanisms, however, show potential for the development of a novel line of PD therapeutics. One of these mechanisms is the peculiar metabolic vulnerability of SNc cells compared to other dopaminergic clusters; the other is the SubThalamic Nucleus (STN)-induced excitotoxicity in SNc. To investigate the latter hypothesis computationally, we developed a spiking neuron network-model of SNc-STN-GPe system. In the model, prolonged stimulation of SNc cells by an overactive STN leads to an increase in ‘stress' variable; when the stress in a SNc neuron exceeds a stress threshold, the neuron dies. The model shows that the interaction between SNc and STN involves a positive-feedback due to which, an initial loss of SNc cells that crosses a threshold causes a runaway-effect, leading to an inexorable loss of SNc cells, strongly resembling the process of neurodegeneration. The model further suggests a link between the two aforementioned mechanisms of SNc cell loss. Our simulation results show that the excitotoxic cause of SNc cell loss might initiate by weak-excitotoxicity mediated by energy deficit, followed by strong-excitotoxicity, mediated by a disinhibited STN. A variety of conventional therapies were simulated to test their efficacy in slowing down SNc cell loss. Among them, glutamate inhibition, dopamine restoration, subthalamotomy and deep brain stimulation showed superior neuroprotective-effects in the proposed model.

A spiking Basal Ganglia model of synchrony, exploration and decision making

To make an optimal decision we need to weigh all the available options, compare them with the current goal, and choose the most rewarding one. Depending on the situation an optimal decision could be to either “explore” or “exploit” or “not to take any action” for which the Basal Ganglia (BG) is considered to be a key neural substrate. In an attempt to expand this classical picture of BG function, we had earlier hypothesized that the Indirect Pathway (IP) of the BG could be the subcortical substrate for exploration. In this study we build a spiking network model to relate exploration to synchrony levels in the BG (which are a neural marker for tremor in Parkinson's disease). Key BG nuclei such as the Sub Thalamic Nucleus (STN), Globus Pallidus externus (GPe) and Globus Pallidus internus (GPi) were modeled as Izhikevich spiking neurons whereas the Striatal output was modeled as Poisson spikes. The model is cast in reinforcement learning framework with the dopamine signal representing reward prediction error. We apply the model to two decision making tasks: a binary action selection task (similar to one used by Humphries et al., 2006) and an n-armed bandit task (Bourdaud et al., 2008). The model shows that exploration levels could be controlled by STN's lateral connection strength which also influenced the synchrony levels in the STN-GPe circuit. An increase in STN's lateral strength led to a decrease in exploration which can be thought as the possible explanation for reduced exploratory levels in Parkinson's patients. Our simulations also show that on complete removal of IP, the model exhibits only Go and No-Go behaviors, thereby demonstrating the crucial role of IP in exploration. Our model provides a unified account for synchronization, action section, and explorative behavior.

Parameter identification for a local field potential driven model of the Parkinsonian subthalamic nucleus spike activity

Several models, with various degrees of complexity have been proposed to model the neuronal activity from different parts of the human brain. We have shown before that various modeling approaches, including a Hammerstein-Wiener (H-W) model, can be used to predict the spike trains from a deep nucleus, the subthalamic nucleus, using the underlying local field potentials. In this article, we present, in depth, the various choices one has to make, and the limitations that they introduce, during the H-W model parameter identification process. From a segment of the recorded data, which contains information about the spike times of a single neuron, we identify and extract the model parameters. We then use those parameters to numerically simulate the spike timing, the rhythm and the inter-spike intervals for the rest of the recording. To assess how well the model fits to the measured data we combine measures of spike train synchrony, namely the Victor-Purpura distance and the Gaussian similarity measure, with time-scale independent train distances. We show that a wise combination of metrics results in models that predict the spikes with temporal accuracy ranging, on average, from 53% to more than 80%, depending on the number of the neurons' spikes recorded. The model's prediction is adequate for estimating accurately the spike rhythm. Quantitative results establish the model's validity as a simple yet biologically plausible model of the spike activity recorded from a deep nucleus inside the human brain.