Signal enhancement in the output stage of the basal ganglia by synaptic short-term plasticity in the direct, indirect, and hyperdirect pathways

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

Introduction

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

Methods

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

Results

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

Discussion/conclusion

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

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.

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.

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.

Transient Response of Basal Ganglia Network in Healthy and Low-Dopamine State

The basal ganglia (BG) are crucial for a variety of motor and cognitive functions. Changes induced by persistent low-dopamine (e.g., in Parkinson’s disease; PD) result in aberrant changes in steady-state population activity (β band oscillations) and the transient response of the BG. Typically, a brief cortical stimulation results in a triphasic response in the substantia nigra pars reticulata (SNr; an output of the BG). The properties of the triphasic responses are shaped by dopamine levels. While mechanisms underlying aberrant steady state activity are well studied, it is still unclear which BG interactions are crucial for the aberrant transient responses in the BG. Moreover, it is also unclear whether mechanisms underlying the aberrant changes in steady-state activity and transient response are the same. Here, we used numerical simulations of a network model of BG to identify the key factors that determine the shape of the transient responses. We show that an aberrant transient response of the SNr in the low-dopamine state involves changes in the direct pathway and the recurrent interactions within the globus pallidus externa (GPe) and between GPe and subthalamic nucleus (STN). However, the connections from D2-type spiny projection neurons (D2-SPN) to GPe are most crucial in shaping the transient response and by restoring them to their healthy level, we could restore the shape of transient response even in low-dopamine state. Finally, we show that the changes in BG that result in aberrant transient response are also sufficient to generate pathologic oscillatory activity in the steady state.

Connectivity and Functionality of the Globus Pallidus Externa Under Normal Conditions and Parkinson's Disease

The globus pallidus externa (GPe) functions as a central hub in the basal ganglia for processing motor and non-motor information through the creation of complex connections with the other basal ganglia nuclei and brain regions. Recently, with the adoption of sophisticated genetic tools, substantial advances have been made in understanding the distinct molecular, anatomical, electrophysiological, and functional properties of GPe neurons and non-neuronal cells. Impairments in dopamine transmission in the basal ganglia contribute to Parkinson's disease (PD), the most common movement disorder that severely affects the patients' life quality. Altered GPe neuron activity and synaptic connections have also been found in both PD patients and pre-clinical models. In this review, we will summarize the main findings on the composition, connectivity and functionality of different GPe cell populations and the potential GPe-related mechanisms of PD symptoms to better understand the cell type and circuit-specific roles of GPe in both normal and PD conditions.

A Basal Ganglia Computational Model to Explain the Paradoxical Sensorial Improvement in the Presence of Huntington's Disease

The basal ganglia (BG) represent a critical center of the nervous system for sensorial discrimination. Although it is known that Huntington's disease (HD) affects this brain area, it still remains unclear how HD patients achieve paradoxical improvement in sensorial discrimination tasks. This paper presents a computational model of the BG including the main nuclei and the typical firing properties of their neurons. The BG model has been embedded within an auditory signal detection task. We have emulated the effect that the altered levels of dopamine and the degree of HD affectation have in information processing at different layers of the BG, and how these aspects shape transient and steady states differently throughout the selection task. By extracting the independent components of the BG activity at different populations, it is evidenced that early and medium stages of HD affectation may enhance transient activity in the striatum and the substantia nigra pars reticulata. These results represent a possible explanation for the paradoxical improvement that HD patients present in discrimination task performance. Thus, this paper provides a novel understanding on how the fast dynamics of the BG network at different layers interact and enable transient states to emerge throughout the successive neuron populations.

Significance and Translational Value of High-Frequency Cortico-Basal Ganglia Oscillations in Parkinson's Disease

The mechanisms and significance of basal ganglia oscillations is a fundamental research question engaging both clinical and basic investigators. In Parkinson’s disease (PD), neural activity in basal ganglia nuclei is characterized by oscillatory patterns that are believed to disrupt the dynamic processing of movement-related information and thus generate motor symptoms. Beta-band oscillations associated with hypokinetic states have been reviewed in several excellent previous articles. Here we focus on faster oscillatory phenomena that have been reported in association with a diverse range of motor states. We review the occurrence of different types of fast oscillations and the evidence supporting their pathophysiological role. We also provide a general discussion on the definition, possible mechanisms, and translational value of synchronized oscillations of different frequencies in cortico-basal ganglia structures. Revealing how oscillatory phenomena are caused and spread in cortico-basal ganglia-thalamocortical networks will offer a key to unlock the neural codes of both motor and non-motor symptoms in PD. In preclinical therapeutic research, recording of oscillatory neural activities holds the promise to unravel mechanisms of action of current and future treatments.

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.

Roles for globus pallidus externa revealed in a computational model of action selection in the basal ganglia

The basal ganglia are considered vital to action selection - a hypothesis supported by several biologically plausible computational models. Of the several subnuclei of the basal ganglia, the globus pallidus externa (GPe) has been thought of largely as a relay nucleus, and its intrinsic connectivity has not been incorporated in significant detail, in any model thus far. Here, we incorporate newly revealed subgroups of neurons within the GPe into an existing computational model of the basal ganglia, and investigate their role in action selection. Three main results ensued. First, using previously used metrics for selection, the new extended connectivity improved the action selection performance of the model. Second, low frequency theta oscillations were observed in the subpopulation of the GPe (the TA or 'arkypallidal' neurons) which project exclusively to the striatum. These oscillations were suppressed by increased dopamine activity - revealing a possible link with symptoms of Parkinson's disease. Third, a new phenomenon was observed in which the usual monotonic relationship between input to the basal ganglia and its output within an action 'channel' was, under some circumstances, reversed. Thus, at high levels of input, further increase of this input to the channel could cause an increase of the corresponding output rather than the more usually observed decrease. Moreover, this phenomenon was associated with the prevention of multiple channel selection, thereby assisting in optimal action selection. Examination of the mechanistic origin of our results showed the so-called 'prototypical' GPe neurons to be the principal subpopulation influencing action selection. They control the striatum via the arkypallidal neurons and are also able to regulate the output nuclei directly. Taken together, our results highlight the role of the GPe as a major control hub of the basal ganglia, and provide a mechanistic account for its control function.

The role of cortical oscillations in a spiking neural network model of the basal ganglia

Although brain oscillations involving the basal ganglia (BG) have been the target of extensive research, the main focus lies disproportionally on oscillations generated within the BG circuit rather than other sources, such as cortical areas. We remedy this here by investigating the influence of various cortical frequency bands on the intrinsic effective connectivity of the BG, as well as the role of the latter in regulating cortical behaviour. To do this, we construct a detailed neural model of the complete BG circuit based on fine-tuned spiking neurons, with both electrical and chemical synapses as well as short-term plasticity between structures. As a measure of effective connectivity, we estimate information transfer between nuclei by means of transfer entropy. Our model successfully reproduces firing and oscillatory behaviour found in both the healthy and Parkinsonian BG. We found that, indeed, effective connectivity changes dramatically for different cortical frequency bands and phase offsets, which are able to modulate (or even block) information flow in the three major BG pathways. In particular, alpha (8–12Hz) and beta (13–30Hz) oscillations activate the direct BG pathway, and favour the modulation of the indirect and hyper-direct pathways via the subthalamic nucleus—globus pallidus loop. In contrast, gamma (30–90Hz) frequencies block the information flow from the cortex completely through activation of the indirect pathway. Finally, below alpha, all pathways decay gradually and the system gives rise to spontaneous activity generated in the globus pallidus. Our results indicate the existence of a multimodal gating mechanism at the level of the BG that can be entirely controlled by cortical oscillations, and provide evidence for the hypothesis of cortically-entrained but locally-generated subthalamic beta activity. These two findings suggest new insights into the pathophysiology of specific BG disorders.

Homologous Basal Ganglia Network Models in Physiological and Parkinsonian Conditions

The classical model of basal ganglia has been refined in recent years with discoveries of subpopulations within a nucleus and previously unknown projections. One such discovery is the presence of subpopulations of arkypallidal and prototypical neurons in external globus pallidus, which was previously considered to be a primarily homogeneous nucleus. Developing a computational model of these multiple interconnected nuclei is challenging, because the strengths of the connections are largely unknown. We therefore use a genetic algorithm to search for the unknown connectivity parameters in a firing rate model. We apply a binary cost function derived from empirical firing rate and phase relationship data for the physiological and Parkinsonian conditions. Our approach generates ensembles of over 1,000 configurations, or homologies, for each condition, with broad distributions for many of the parameter values and overlap between the two conditions. However, the resulting effective weights of connections from or to prototypical and arkypallidal neurons are consistent with the experimental data. We investigate the significance of the weight variability by manipulating the parameters individually and cumulatively, and conclude that the correlation observed between the parameters is necessary for generating the dynamics of the two conditions. We then investigate the response of the networks to a transient cortical stimulus, and demonstrate that networks classified as physiological effectively suppress activity in the internal globus pallidus, and are not susceptible to oscillations, whereas parkinsonian networks show the opposite tendency. Thus, we conclude that the rates and phase relationships observed in the globus pallidus are predictive of experimentally observed higher level dynamical features of the physiological and parkinsonian basal ganglia, and that the multiplicity of solutions generated by our method may well be indicative of a natural diversity in basal ganglia networks. We propose that our approach of generating and analyzing an ensemble of multiple solutions to an underdetermined network model provides greater confidence in its predictions than those derived from a unique solution, and that projecting such homologous networks on a lower dimensional space of sensibly chosen dynamical features gives a better chance than a purely structural analysis at understanding complex pathologies such as Parkinson's disease.

Untangling Basal Ganglia Network Dynamics and Function: Role of Dopamine Depletion and Inhibition Investigated in a Spiking Network Model

The basal ganglia are a crucial brain system for behavioral selection, and their function is disturbed in Parkinson's disease (PD), where neurons exhibit inappropriate synchronization and oscillations. We present a spiking neural model of basal ganglia including plausible details on synaptic dynamics, connectivity patterns, neuron behavior, and dopamine effects. Recordings of neuronal activity in the subthalamic nucleus and Type A (TA; arkypallidal) and Type I (TI; prototypical) neurons in globus pallidus externa were used to validate the model. Simulation experiments predict that both local inhibition in striatum and the existence of an indirect pathway are important for basal ganglia to function properly over a large range of cortical drives. The dopamine depletion-induced increase of AMPA efficacy in corticostriatal synapses to medium spiny neurons (MSNs) with dopamine receptor D2 synapses (CTX-MSN D2) and the reduction of MSN lateral connectivity (MSN-MSN) were found to contribute significantly to the enhanced synchrony and oscillations seen in PD. Additionally, reversing the dopamine depletion-induced changes to CTX-MSN D1, CTX-MSN D2, TA-MSN, and MSN-MSN couplings could improve or restore basal ganglia action selection ability. In summary, we found multiple changes of parameters for synaptic efficacy and neural excitability that could improve action selection ability and at the same time reduce oscillations. Identification of such targets could potentially generate ideas for treatments of PD and increase our understanding of the relation between network dynamics and network function.

Functional Relevance of Different Basal Ganglia Pathways Investigated in a Spiking Model with Reward Dependent Plasticity

The brain enables animals to behaviorally adapt in order to survive in a complex and dynamic environment, but how reward-oriented behaviors are achieved and computed by its underlying neural circuitry is an open question. To address this concern, we have developed a spiking model of the basal ganglia (BG) that learns to dis-inhibit the action leading to a reward despite ongoing changes in the reward schedule. The architecture of the network features the two pathways commonly described in BG, the direct (denoted D1) and the indirect (denoted D2) pathway, as well as a loop involving striatum and the dopaminergic system. The activity of these dopaminergic neurons conveys the reward prediction error (RPE), which determines the magnitude of synaptic plasticity within the different pathways. All plastic connections implement a versatile four-factor learning rule derived from Bayesian inference that depends upon pre- and post-synaptic activity, receptor type, and dopamine level. Synaptic weight updates occur in the D1 or D2 pathways depending on the sign of the RPE, and an efference copy informs upstream nuclei about the action selected. We demonstrate successful performance of the system in a multiple-choice learning task with a transiently changing reward schedule. We simulate lesioning of the various pathways and show that a condition without the D2 pathway fares worse than one without D1. Additionally, we simulate the degeneration observed in Parkinson's disease (PD) by decreasing the number of dopaminergic neurons during learning. The results suggest that the D1 pathway impairment in PD might have been overlooked. Furthermore, an analysis of the alterations in the synaptic weights shows that using the absolute reward value instead of the RPE leads to a larger change in D1.

Control effects of stimulus paradigms on characteristic firings of parkinsonism

Experimental studies have shown that neuron population located in the basal ganglia of parkinsonian primates can exhibit characteristic firings with certain firing rates differing from normal brain activities. Motivated by recent experimental findings, we investigate the effects of various stimulation paradigms on the firing rates of parkinsonism based on the proposed dynamical models. Our results show that the closed-loop deep brain stimulation is superior in ameliorating the firing behaviors of the parkinsonism, and other control strategies have similar effects according to the observation of electrophysiological experiments. In addition, in conformity to physiological experiments, we found that there exists optimal delay of input in the closed-loop GPtrain|M1 paradigm, where more normal behaviors can be obtained. More interestingly, we observed that W-shaped curves of the firing rates always appear as stimulus delay varies. We furthermore verify the robustness of the obtained results by studying three pallidal discharge rates of the parkinsonism based on the conductance-based model, as well as the integrate-and-fire-or-burst model. Finally, we show that short-term plasticity can improve the firing rates and optimize the control effects on parkinsonism. Our conclusions may give more theoretical insight into Parkinson's disease studies.