A population level computational model of the basal ganglia that generates parkinsonian Local Field Potential activity

Estimation of the transmission delays in the basal ganglia of the macaque monkey and subsequent predictions about oscillatory activity under dopamine depletion

The timescales of the dynamics of a system depend on the combination of the timescales of its components and of its transmission delays between components. Here, we combine experimental stimulation data from 10 studies in macaque monkeys that reveal the timing of excitatory and inhibitory events in the basal ganglia circuit, to estimate its set of transmission delays. In doing so, we reveal possible inconsistencies in the existing data, calling for replications, and we propose two possible sets of transmission delays. We then integrate these delays in a model of the primate basal ganglia that does not rely on direct and indirect pathways' segregation and show that extrastriatal dopaminergic depletion in the external part of the globus pallidus and in the subthalamic nucleus is sufficient to generate β-band oscillations (in the high part, 20-35 Hz, of the band). More specifically, we show that D2 and D5 dopamine receptors in these nuclei play opposing roles in the emergence of these oscillations, thereby explaining how completely deactivating D5 receptors in the subthalamic nucleus can, paradoxically, cancel oscillations.

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.

A biologically constrained spiking neural network model of the primate basal ganglia with overlapping pathways exhibits action selection

Action selection has been hypothesized to be a key function of the basal ganglia, yet the nuclei involved, their interactions and the importance of the direct/indirect pathway segregation in such process remain debated. Here, we design a spiking computational model of the monkey basal ganglia derived from a previously published population model, initially parameterized to reproduce electrophysiological activity at rest and to embody as much quantitative anatomical data as possible. As a particular feature, both models exhibit the strong overlap between the direct and indirect pathways that has been documented in non-human primates. Here, we first show how the translation from a population to an individual neuron model was achieved, with the addition of a minimal number of parameters. We then show that our model performs action selection, even though it was built without any assumption on the activity carried out during behaviour. We investigate the mechanisms of this selection through circuit disruptions and found an instrumental role of the off-centre/on-surround structure of the MSN-STN-GPi circuit, as well as of the MSN-MSN and FSI-MSN projections. This validates their potency in enabling selection. We finally study the pervasive centromedian and parafascicular thalamic inputs that reach all basal ganglia nuclei and whose influence is therefore difficult to anticipate. Our model predicts that these inputs modulate the responsiveness of action selection, making them a candidate for the regulation of the speed-accuracy trade-off during decision-making.

Neural mass models describing possible origin of the excessive beta oscillations correlated with Parkinsonian state

In Parkinson's disease, the enhanced beta rhythm is closely associated with akinesia/bradykinesia and rigidity. An increase in beta oscillations (12-35 Hz) within the basal ganglia (BG) nuclei does not proliferate throughout the cortico-basal ganglia loop in uniform fashion; rather it can be subdivided into two distinct frequency bands, i.e. the lower beta (12-20 Hz) and upper beta (21-35 Hz). A computational model of the excitatory and inhibitory neural network that focuses on the population properties is proposed to explore the mechanism underlying the pathological beta oscillations. Simulation results show several findings. The upper beta frequency in the BG originates from a high frequency cortical beta, while the emergence of exaggerated lower beta frequency in the BG depends greatly on the enhanced excitation of a reciprocal network consisting of the globus pallidus externus (GPe) and the subthalamic nucleus (STN). There is also a transition mechanism between the upper and lower beta oscillatory activities, and we explore the impact of self-inhibition within the GPe on the relationship between the upper beta and lower beta oscillations. It is shown that increased self-inhibition within the GPe contributes to increased upper beta oscillations driven by the cortical rhythm, while decrease in the self-inhibition within the GPe facilitates an enhancement of the lower beta oscillations induced by the increased excitability of the BG. This work provides an analysis for understanding the mechanism underlying pathological synchronization in neurological diseases.

A neural mass model of basal ganglia nuclei simulates pathological beta rhythm in Parkinson's disease

An increase in beta oscillations within the basal ganglia nuclei has been shown to be associated with movement disorder, such as Parkinson's disease. The motor cortex and an excitatory-inhibitory neuronal network composed of the subthalamic nucleus (STN) and the external globus pallidus (GPe) are thought to play an important role in the generation of these oscillations. In this paper, we propose a neuron mass model of the basal ganglia on the population level that reproduces the Parkinsonian oscillations in a reciprocal excitatory-inhibitory network. Moreover, it is shown that the generation and frequency of these pathological beta oscillations are varied by the coupling strength and the intrinsic characteristics of the basal ganglia. Simulation results reveal that increase of the coupling strength induces the generation of the beta oscillation, as well as enhances the oscillation frequency. However, for the intrinsic properties of each nucleus in the excitatory-inhibitory network, the STN primarily influences the generation of the beta oscillation while the GPe mainly determines its frequency. Interestingly, describing function analysis applied on this model theoretically explains the mechanism of pathological beta oscillations.

Analysis of Oscillatory Neural Activity in Series Network Models of Parkinson's Disease During Deep Brain Stimulation

Parkinson's disease is a progressive, neurodegenerative disorder, characterized by hallmark motor symptoms. It is associated with pathological, oscillatory neural activity in the basal ganglia. Deep brain stimulation (DBS) is often successfully used to treat medically refractive Parkinson's disease. However, the selection of stimulation parameters is based on qualitative assessment of the patient, which can result in a lengthy tuning period and a suboptimal choice of parameters. This study explores fourth-order, control theory-based models of oscillatory activity in the basal ganglia. Describing function analysis is applied to examine possible mechanisms for the generation of oscillations in interacting nuclei and to investigate the suppression of oscillations with high-frequency stimulation. The theoretical results for the suppression of the oscillatory activity obtained using both the fourth-order model, and a previously described second-order model, are optimized to fit clinically recorded local field potential data obtained from Parkinsonian patients with implanted DBS. Close agreement between the power of oscillations recorded for a range of stimulation amplitudes is observed ( R(2)=0.69-0.99 ). The results suggest that the behavior of the system and the suppression of pathological neural oscillations with DBS is well described by the macroscopic models presented. The results also demonstrate that in this instance, a second-order model is sufficient to model the clinical data, without the need for added complexity. Describing the system behavior with computationally efficient models could aid in the identification of optimal stimulation parameters for patients in a clinical environment.

A biologically constrained model of the whole basal ganglia addressing the paradoxes of connections and selection

The basal ganglia nuclei form a complex network of nuclei often assumed to perform selection, yet their individual roles and how they influence each other is still largely unclear. In particular, the ties between the external and internal parts of the globus pallidus are paradoxical, as anatomical data suggest a potent inhibitory projection between them while electrophysiological recordings indicate that they have similar activities. Here we introduce a theoretical study that reconciles both views on the intra-pallidal projection, by providing a plausible characterization of the relationship between the external and internal globus pallidus. Specifically, we developed a mean-field model of the whole basal ganglia, whose parameterization is optimized to respect best a collection of numerous anatomical and electrophysiological data. We first obtained models respecting all our constraints, hence anatomical and electrophysiological data on the intrapallidal projection are globally consistent. This model furthermore predicts that both aforementioned views about the intra-pallidal projection may be reconciled when this projection is weakly inhibitory, thus making it possible to support similar neural activity in both nuclei and for the entire basal ganglia to select between actions. Second, we predicts that afferent projections are substantially unbalanced towards the external segment, as it receives the strongest excitation from STN and the weakest inhibition from the striatum. Finally, our study strongly suggests that the intrapallidal connection pattern is not focused but diffuse, as this latter pattern is more efficient for the overall selection performed in the basal ganglia.

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.

System identification of local field potentials under deep brain stimulation in a healthy primate

High frequency (HF) Deep Brain Stimulation (DBS) in the Sub-Thalamic Nucleus (STN) is a clinically recognized therapy for the treatment of motor disorders in Parkinson Disease (PD). The underlying mechanisms of DBS and how it impacts neighboring nuclei, however, are not yet completely understood. Electrophysiological data has been collected in PD patients and primates to better understand the impact of DBS on STN and the entire Basal Ganglia (BG) motor circuit. We use single unit recordings from Globus Pallidus, both pars interna and externa segments (GPi and GPe) in the BG, in a normal primate before and after DBS to reconstruct Local Field Potentials (LFPs) in the region. We then use system identification techniques to understand how GPe LFP activity and the DBS signal applied to STN influence GPi LFP activity. Our models suggest that when no stimulation is applied, the GPe LFPs have an inhibitory effect on GPi LFPs with a 2-3 ms delay, as is the case for single unit neuronal activity. On the other hand, when DBS is ON the models suggest that stimulation has a dominant effect on GPi LFPs which mask the inhibitory effects of GPe.