ORIGINAL RESEARCH article
Computational Model of Recurrent Subthalamo-Pallidal Circuit for Generation of Parkinsonian Oscillations
- 1Okinawa Institute of Science and Technology Graduate University, Okinawa, Japan
- 2Honda Research Institute Japan Co., Ltd., Saitama, Japan
- 3Division of System Neurophysiology, Department of Physiological Sciences, National Institute for Physiological Sciences, Graduate University for Advanced Studies, Aichi, Japan
Parkinson's disease is a movement disorder caused by dopamine depletion in the basal ganglia. Abnormally synchronized neuronal oscillations between 8 and 15 Hz in the basal ganglia are implicated in motor symptoms of Parkinson's disease. However, how these abnormal oscillations are generated and maintained in the dopamine-depleted state is unknown. Based on neural recordings in a primate model of Parkinson's disease and other experimental and computational evidence, we hypothesized that the recurrent circuit between the subthalamic nucleus (STN) and the external segment of the globus pallidus (GPe) generates and maintains parkinsonian oscillations, and that the cortical excitatory input to the STN amplifies them. To investigate this hypothesis through computer simulations, we developed a spiking neuron model of the STN-GPe circuit by incorporating electrophysiological properties of neurons and synapses. A systematic parameter search by computer simulation identified regions in the space of the intrinsic excitability of GPe neurons and synaptic strength from the GPe to the STN that reproduce normal and parkinsonian states. In the parkinsonian state, reduced firing of GPe neurons and increased GPe-STN inhibition trigger burst activities of STN neurons with strong post-inhibitory rebound excitation, which is usually subject to short-term depression. STN neuronal bursts are shaped into the 8–15 Hz, synchronous oscillations via recurrent interactions of STN and GPe neurons. Furthermore, we show that cortical excitatory input to the STN can amplify or suppress pathological STN oscillations depending on their phase and strength, predicting conditions of cortical inputs to the STN for suppressing oscillations.
Parkinson's disease (PD) is characterized by progressive movement disorders caused by dopamine depletion in the basal ganglia (BG) resulting from degeneration of midbrain dopaminergic neurons. Abnormally synchronized oscillations in local field potentials (LFPs) between 8 and 30 Hz, which are called “the BG beta frequency band” (Hammond et al., 2007), are observed in the subthalamic nucleus (STN) and the internal/external segments of the globus pallidus (GPi/GPe) of PD patients (Miller and DeLong, 1987; Levy et al., 2001; Kühn et al., 2005). These abnormal oscillations are implicated in PD symptoms, such as akinesia, bradykinesia, and rigidity (Boraud et al., 2002; Brown, 2003; Gatev et al., 2006; Rivlin-Etzion et al., 2006; Hammond et al., 2007; but see Leblois et al., 2006b). We constructed a spiking neuron model of the recurrent STN-GPe circuit and investigated mechanisms for generation of parkinsonian beta band oscillations by exploring physiological parameters.
BG beta oscillations in LFPs are divided into upper-beta (>20 Hz) and lower-beta (<20 Hz) oscillations, and only the lower-beta oscillations are sensitive to dopaminergic therapy (Priori et al., 2004). Single-unit recordings from the STN in human PD patients show abnormal neuronal burst oscillations in the 10–25 Hz frequency band, which could be the source of pathological LFP oscillations (Levy et al., 2001, 2002). Neuronal recordings in nonhuman primate PD models show exaggerated oscillatory burst activities in STN, GPe, and GPi that are concentrated in narrow frequency bands of 4–8 Hz and 8–15 Hz (Raz et al., 2000; Heimer et al., 2006; Wichmann and Soares, 2006; Tachibana et al., 2011), the latter of which overlaps with the frequency band recorded from human PD patients. These 8–15 Hz neuronal burst oscillations are reduced by both dopaminergic therapy and deep brain stimulation, which can alleviate PD motor symptoms (Heimer et al., 2006; Rosin et al., 2011; Tachibana et al., 2011). Recently, neuronal recordings from primate PD models under a variety of pharmacological manipulations demonstrated that STN 8–15 Hz neuronal oscillations are the origin of GPe/GPi oscillations and that they are dependent on GABAergic inhibition from GPe and glutamatergic inputs from the cortex/thalamus (Tachibana et al., 2011).
Computational modeling studies have tackled mechanisms of pathological BG oscillations (Gillies et al., 2002; Terman et al., 2002; Rubin and Terman, 2004; Leblois et al., 2006a; Hahn and McIntyre, 2010; Holgado et al., 2010; Kumar et al., 2011). However, their simulated results did not agree well with experimentally-observed neuronal properties of the BG 8–15 Hz burst oscillations (Wichmann and Soares, 2006; Tachibana et al., 2011). Furthermore, the proposed role of strengthened striatal inhibition for the induction of the BG oscillations (Terman et al., 2002; Kumar et al., 2011) is not consistent with the experimental evidence that the blockade of GABAergic synaptic transmissions of GPe neurons from the striatum and other GPe neurons did not affect the 8–15 Hz oscillations (Tachibana et al., 2011). Therefore, the manner in which the pathological 8–15 Hz burst oscillations are generated remains unclear. Based on experimentally-observed neuronal data (Tachibana et al., 2011), we investigated mechanisms for generation and maintenance of pathological 8–15 Hz burst oscillations in recurrent STN-GPe circuits in a dopamine-depleted state. Using computer simulations of a spiking neuron model of the recurrent STN-GPe circuit, we explored the required neuronal, synaptic, and circuit mechanisms for normal firing and pathological burst oscillation. We then investigated how cortical excitatory input to the STN affects BG oscillations generated in the STN-GPe circuit and showed that for two parameter regions, the model reproduces BG neuronal activity in terms of mean firing rates, oscillation frequency, and burst discharges, corresponding to normal and parkinsonian states. The same mechanism for generation of BG pathological oscillations also predicts suppressive effects of non-oscillatory cortical inputs to the BG oscillations that explains suppression of pathological BG oscillations by movements (Amirnovin et al., 2004; Brown and Williams, 2005). We further demonstrated that oscillatory inputs from the cortex to the STN can amplify or diminish pathological BG oscillations in a phase-dependent manner, and can predict necessary conditions for their suppression.
Materials and Methods
Experimental evidence suggests that reciprocal STN-GPe interconnections and glutamatergic inputs to the STN in dopamine-depleted states are both important for generation and amplification of the parkinsonian, oscillatory burst activity of STN neurons that subsequently propagates to the GPi (Tachibana et al., 2011). Post-inhibitory rebound excitation of STN neurons is a candidate driving force of the oscillatory burst activity (Bevan et al., 2002), which could account for the importance of GPe inhibition to the STN for generation of the BG parkinsonian oscillatory burst activity of STN neurons (Tachibana et al., 2011). Short-term plasticity at excitatory synapses of GPe neurons and at GABAergic synapses of STN neurons are implicated in modulating temporal patterns of GPe and STN neurons in the normal state, and as well as in the parkinsonian state (Hanson and Jaeger, 2002; Atherton et al., 2013). Dopamine depletion in the BG is the key phenomenon for the induction of parkinsonian symptoms and BG pathological oscillations (Heimer et al., 2006; Wichmann and Soares, 2006; Tachibana et al., 2011). We considered the following neural phenomena that are related to BG dopamine depletion as candidate causes of the BG pathological oscillations: (1) reduction of GPe autonomous activity (Chan et al., 2011), (2) strengthening of GPe inhibition to the STN due to release from D2-receptor-mediated presynaptic suppression of transmitter release (Shen and Johnson, 2000, 2005), and by an increase in the number of synaptic connections per GPe-STN axon terminal (Fan et al., 2012), (3) increased STN-GPe excitation to compensate for the reduction of GPe autonomous activity (Chan et al., 2011), and (4) loss of active decorrelation ability of GPe neurons (Edgerton and Jaeger, 2011). Although elevated striatal inhibition of GPe neurons by dopamine depletion (Hernandez-Lopez et al., 2000; Liang et al., 2008) is often linked to pathological BG oscillations (Terman et al., 2002; Kumar et al., 2011), we excluded it from candidate neural phenomena related to generation of BG oscillations because experimental results showed that blockade of ionotropic GABAergic transmission either did not affect or even exaggerated GPe 8–15 Hz oscillations (Tachibana et al., 2011). Based on this evidence, we focused on a model of closed STN-GPe circuits incorporating the STN post-inhibitory rebound mechanism and dynamic synapses as mechanisms for generation of parkinsonian BG oscillations in the dopamine-depleted state. First, we explored required conditions of GPe autonomous activity and of connection strengths between the STN and GPe for generation of BG oscillations in the closed STN-GPe circuits. Later, we investigated roles of cortical excitatory inputs in the STN for BG oscillations that were generated in closed STN-GPe circuits.
Anatomically, the STN and GPe are reciprocally connected. STN neurons are a major source of excitatory inputs to the GPe, and GPe neurons are a major source of inhibitory inputs to the STN. These connections are very sparse and their average connectivity is low (Kita and Kitai, 1994; Bevan et al., 1998; Sato et al., 2000; Baufreton et al., 2009). Based on these findings, we developed a random, sparsely-connected architecture of the recurrent circuit model between the STN an GPe. The architecture of our network model is shown in Figure 1A. The model consisted of 64 STN neurons and 192 GPe neurons unless otherwise noted. Each STN neuron received GABAergic projections from 12 randomly-selected GPe neurons. GPe-STN connectivity is of the same order of magnitude as experimental results (Baufreton et al., 2009). Each GPe neuron received glutamatergic projections from 6 randomly selected STN neurons. We did not include GABAergic connections between GPe neurons in our model because they were controversial: Anatomical studies showed inhibitory axonal collaterals of GPe neurons (Bevan et al., 1998; Sato et al., 2000), while functional interactions between GPe neurons were very weak or nonexistent. Each glutamatergic projection onto a GPe neuron gave a single synaptic contact, and exhibited short-term plasticity (Hanson and Jaeger, 2002; Figure 1B). Each GABAergic projection from the GPe onto an STN neuron gave N synaptic contacts that were randomly drawn from a log-normal distribution of mean = 15.50 and standard deviation = 14.53 (Figure 1C), which fit experimental data well (Baufreton et al., 2009). Each GPe neuron had intrinsic noisy conductance (GPe Int.), as will hereinafter be described in detail. A population of cortical excitatory inputs to the STN was modeled as a stochastic spike train that obeyed a Poisson process, the mean rate of which was constant at 100 Hz, unless otherwise noted. GPe-STN GABAergic synapses exhibited short-term depression (Atherton et al., 2013). Axonal conduction delays of STN-GPe, and GPe-STN were 5 ms (Fujimoto and Kita, 1993; Kita et al., 2005).
Figure 1. (A) A random, sparsely-connected architecture of the recurrent circuit model between the STN and GPe. (B) Properties of three different types of short-term plasticity at STN-GPe synapses. A type of short-term plasticity is randomly chosen for a STN-GPe synapse with equal probability from these three types regardless of a post-synaptic neuron. (C) A single GPe axon gives multiple synaptic contacts on a STN neuron. Transmission probability of a unitary GPe projection is dynamically modulated by short-term depression. The number of synaptic contacts that a single GPe axon established on an STN neuron obeys the log-normal distribution.
Neuron and Synapse Models
An STN neuron was simulated by using a single-compartment, conductance-based model (Otsuka et al., 2004) that contains several sets of currents:
where Cm is membrane capacitance and takes 1 μF/cm2. V is the membrane potential. INa is a voltage-dependent sodium current. IK is a voltage-dependent, delayed rectifier potassium current. IA is an A-type potassium current. ILCa is an L-like, long-lasting calcium current. IT is a low threshold T-type calcium current. ICa−k is a calcium-activated potassium current. a, b, c, d1, d2, h, m, n, p, q, and r are activation or inactivation variables of these ionic currents. VNa, VK, and VCa are the equilibrium potentials of sodium, potassium, and calcium, respectively. Vleak is the resting potential of the leak current. Vexc and Vinh are excitatory and inhibitory synaptic reversal potentials. gNa, gK, gA, gLCa, gT, and gleak are the maximal conductances of INa, IK, IA, ILCa, IT and Ileak, respectively. Parameter values of the STN model cell were identical to the original ones except for excitatory and inhibitory synaptic reversal potentials (Vexc = 0.0 mV and Vinh = −84.0 mV), and the surface area of a cell, Am. Although the original model is for a unit-surface area of membrane, we set Am = 0.2 μm2 so as to match the experimentally observed relationship between firing rate and applied current intensity (Otsuka et al., 2004). Excitatory synaptic conductance mediated by AMPA-type glutamatergic receptors, gAMPA, was modeled as the alpha function:
We chose a synaptic time constant and the maximal synaptic conductance of an excitatory synapse as τAMPA = 1.0 ms, ḡAMPA−Ctx−STN = 1.0 nS. Although a physiological study shows that an STN neuron receives a substantial number of excitatory inputs via NMDA-type glutamatergic receptors in addition to AMPA-type-receptor-mediated inputs (Nambu et al., 2000), we did not simulate NMDA-type conductance in this study because preliminary simulation results did not show any significant differences with or without it. The GABAergic synaptic conductance, gGABA, was modeled as the beta function:
We chose τdecay = 7.7 ms and τrise = 0.380 ms for GABAergic synapses to fit experimental data (Baufreton et al., 2009). A GPe neuron was implemented as a single-compartment conductance-based leaky integrate-and-fire model:
Parameter values of GPe model neurons were as follows: Cm = 0.25 pF, gleak = 16.66 nS, Vleak = −70.0 mV, Vexc = 0.0 mV. The absolute refractory period was set to 2.0 ms. Threshold and reset potentials were −55.0 and −60.0 mV, respectively. Synaptic conductances of GPe neurons were modeled as the beta function, after Baufreton et al. (2009):
We chose τdecay = 12.4 ms and τrise = 5.0 ms (Hanson and Jaeger, 2002). ḡAMPA−STN−GPe and ḡAMPA−GPe−Int represent the maximal synaptic conductances of STN-to-GPe synapses and the GPe intrinsic noisy conductance. The GPe intrinsic noisy conductance (GPe Int. for short) was introduced to model autonomous activity of GPe neurons, and was implemented as a pseudo excitatory synaptic conductance that responds to a pseudo spike train obeying a Poisson process in order to reproduce both facilitating and desynchronizing effects on GPe neuronal activity. The mean rate of a pseudo spike train was set as 100 Hz, unless otherwise noted. The dopamine depletion-related reduction of GPe autonomous activity was modeled as the reduction of the maximal conductance of the GPe intrinsic noisy conductance. Although there were other ways to implement autonomous activity, such as using single-compartment or multi-compartment conductance-based models incorporating hyperpolarization-activated, cyclic nucleotide-gated (HCN) channels, we chose this simple implementation because our focus was generation of STN parkinsonian oscillations. The means by which significant variability in intrinsic electrophysiological properties of GPe neurons (Günay et al., 2008) affects generation of parkinsonian oscillations was beyond the scope of this paper.
Short-Term Plasticity Models
Short-term depression at a GPe-to-STN GABAergic synapse was simulated using a phenomenological model:
where P was the transmission probability of a unitary GPe projection. At every arrival of a presynaptic spike, P was decreased by multiplication with a decrement factor as follows:
where the parameter fD < 1 and the depression of P was limited by a lower bound, Pbound. This implementation smoothly increased the decrement factor from fD to 1 as the lower bound for P was reached, and smaller values of the parameter, nbound, retarded the increase of the decrement factor. At each synaptic contact of a unitary GPe projection, a spike was successfully transmitted with the probability,
assuming that the transmission probability of a unitary GPe projection, P, obeyed the binomial distribution of transmissions at Ncontact synaptic contacts made by a unitary GPe projection. Parameters of these models were manually tuned to fit experimental data (Atherton et al., 2013; Figure 2).
Figure 2. Properties of the short-term depression model at a GPe-STN synapse. Plots show simulated transmission probability of a single GPe axon stimulated at 1, 10, 20, 33, and 100 Hz for 10 s. (A) Depression during stimulation. (B) Recovery from the depression. Crosses represent simulation results, and circles with error bars represent means and standard deviations of experimental data (Atherton et al., 2013).
We used the phenomenological model by Hanson and Jaeger (2002) to simulate short-term plasticity at STN-to-GPe glutamatergic synapses. In this model, response amplitude, A, was a product of a facilitation variable, F, and a depression variable, D (F > 1 and D < 1) that decayed exponentially:
At every arrival of a presynaptic spike, the facilitation variable, F, was increased through
multiplication by an incremental factor as follows:
where the parameter fF > 1 and the increase of F was limited by an upper bound, Fbound. At every arrival of a presynaptic spike, the depression variable, D, was decreased by multiplication by the parameter, fD (fD < 1). According to Hanson and Jaeger (2002), each synapse was randomly assigned to one of three types of plasticity, facilitation-dominant (type 1, f F = 1.4, f D = 0.9, τF = 241.0 ms and τD = 491.0 ms), intermediate (type 2, fF = 1.34, fD = 0.86, τF = 345.0 ms and τD = 700.0 ms), and depression-dominant (type 3, f F = 1.64, f D = 0.55, τF = 148.0 ms, and τD = 764.0 ms) with the same probability regardless of a post-synaptic neuron.
Data Analysis and Computer Simulations
Oscillatory activity of neurons, and oscillatory correlations of pairs of neurons were estimated by the power and cross-spectral analyses of spike trains with the shuffling method (Rivlin-Etzion et al., 2006). The compensated power spectral density (PSD) was calculated by dividing the spectrum of the original spike train by the mean spectrum of locally (T = 175–225 ms) shuffled (n = 50) spike trains (Tachibana et al., 2011). A confidence level for compensated PSD was calculated based on the mean and standard deviation of the spectrum in the range of 270–300 Hz (Rivlin-Etzion et al., 2006). A cell was classified as “oscillatory” when its compensated PSD had values above the confidence level (p = 0.01) in at least two adjacent frequency bins of 1,000 Hz. The compensated cross-spectral density (CSD) was calculated by dividing the CSD of the original spike train by the mean CSD of globally shuffled (n = 20) spike trains. A confidence level for compensated CSD was constructed based on the range of 270–300 Hz. A pair of neurons was classified as having “correlated oscillations” if the compensated CSD contained at least two adjacent bins within a frequency range, the values of which were above the confidence level (p = 0.01). Bursts were detected by the “Poisson surprise” method (Wichmann and Soares, 2006). A Poisson surprise value of 3 for a given sequence of at least two interspike intervals (three spikes) was used to separate burst from non-burst spike trains.
All neuronal network simulations were performed in NEST (http://www.nest-initiative.org), and analyses of simulation results were performed in R and Python using Numpy, Scipy, Matplotlib, rpy2, and statsmodels libraries.
We constructed a spiking neuron model of the STN-GPe network by incorporating known physiological and anatomical properties of their neurons and synapses (Figures 1, 2, see Methods for details). Specifically, we incorporated short-term facilitation and depression of STN-to-GPe synapses (Figure 1B; Hanson and Jaeger, 2002) and short-term depression of GPe-to-STN synapses (Figures 1C, 2; Atherton et al., 2013). We considered changes in strengths of STN-GPe and GPe-STN synapses and intrinsic excitability of GPe neurons (GPe Int.) as parameters possibly affected by dopamine depletion, and thereby potentially responsible for the 8–15 Hz oscillatory bursts.
We systematically searched the parameter space for regions that reproduce neuronal behaviors reported in Tachibana et al. (2011), namely, the mean firing rate, the peak power of 8–15 Hz oscillation, and the number of bursts per second (burst frequency) in normal and dopamine-depleted monkeys (Table 1). We found that there are distinct parameter regions where simulated neuronal activities of STN and GPe exhibited properties similar to those experimentally observed in normal and parkinsonian states (Figure 3). In the normal state (ḡAMPA−GPe−Int = 1.125 nS, ḡGABA−GPe−STN = 3.65 nS, ḡAMPA−STN−GPe = 1.05 nS; open squares in Figure 3), electrophysiological properties of STN and GPe neurons of the model were similar to those in normal animals (Table 1), although the proportion of GPe neurons showing 8–15 Hz oscillation was >10x larger than those observed in animal experiments. The parkinsonian state was reproduced with a lower GPe intrinsic excitability and a higher GPe-STN synaptic strength (ḡAMPA−GPe−Int = 0.0 nS, ḡGABA−GPe−STN = 7.3 nS, ḡAMPA−STN−GPe = 1.35 nS; open circles in Figure 3), where STN and GPe neuronal activities showed higher burst frequency and stronger 8–15 Hz oscillations than in the normal state. In this parameter region, mean firing rates of simulated neurons in both structures (STN and GPe) were well matched to mean firing rates of neurons experimentally recorded from dopamine-depleted, parkinsonian-state monkeys (Tachibana et al., 2011, Table 1; Figure 3). The peak frequency of 8–15 Hz oscillations with maximum power ranged from 14 to 15 Hz in both structures (Figure 4B, Table 1). This range of peak frequencies is very similar to pathological oscillations observed in dopamine-depleted monkeys (Raz et al., 2000; Heimer et al., 2006; Rosin et al., 2011; Tachibana et al., 2011). No or very few neurons exhibited significant 3–8 Hz oscillations under the conditions tested.
Table 1. Electrophysiological properties of STN and GPe neurons from animal experiments (A), and model simulations (B).
Figure 3. Comparison of conditions for properties of neuronal activity with simulations of a network model for different parameter values. A column of two panels corresponds to the feature of neuronal activity noted above. In each panel, each data point represents the mean of 8 distinct simulations with values of maximal conductances of GPe-STN synapses and GPe intrinsic noisy conductances (GPe Int.) that control levels of GPe autonomous activity, shown on the axes. In the left column, a contour plot in each panel represents the mean firing rates of neurons in simulation results, and light red and light blue patches behind the contour plot indicate parameter areas for which simulation results are within the means ± standard deviation of firing rates of neurons recorded from the BG of parkinsonian and normal monkeys (Tachibana et al., 2011), respectively. An overlap between light red and light blue patches appears slate gray. Grey circles and squares indicate parameters that correspond to the parkinsonian and normal states, respectively.
Figure 4. Typical network and neuronal activities in simulations in the normal (A) and parkinsonian states (B) (A,B), Left, four panels show raster plots of STN and GPe neurons. Although a network model contains 192 neurons for the GPe, activities of 64 neurons are shown for the sake of clarity. In the middle column, each panel shows the power spectral density functions of two neurons in the STN or GPe. Heavy and light blue traces indicate the power spectral density functions of two neurons, respectively. Heavy and light red dotted lines indicate the confidence levels of the power spectral density functions of two neurons, respectively. In the right column, a top panel shows the absolute cross-spectral density functions of 196 randomly-selected pairs of STN and GPe neurons in a simulation, and a bottom panel shows the instantaneous phase values for the correlated pairs of STN and GPe neurons.
In the normal state, synchronous oscillatory burst discharges were not observed, and no neurons exhibited 8–15 Hz oscillations (Figure 4A). Burst discharges were also detected in the normal state, but they occurred in a rather irregular manner and less frequently, and their durations were longer than those in the parkinsonian state. In contrast, in the parkinsonian state, neurons showed synchronous oscillatory burst discharges, with different timing in different nuclei (Figure 4B). Each neuron exhibited 14–15 Hz burst oscillation and frequent burst firings. Cross-spectral density functions of pairs of STN and GPe neurons peaked at 14–15 Hz, as in the power spectral density functions of STN and GPe neurons (Figure 4B). Instantaneous phase differences for STN vs. GPe pairs clustered around 90°, indicating that the maximal firing of GPe neurons lagged behind that of STN neurons by ~17 ms (Figure 4B, right panels), which is comparable to experimental observations in rodents (Mallet et al., 2008).
Mechanisms for Generation of 8–15 Hz Oscillatory Burst Firings
To understand neuronal and synaptic mechanisms underlying 14–15 Hz oscillatory burst firings in the parkinsonian state, we focused on a local sub-circuit in the STN-GPe network, an STN neuron (STN #46), an afferent (presynaptic) GPe neuron (GPe #41), and an efferent (postsynaptic) GPe neuron (GPe #31). Under normal conditions, STN, and GPe neurons did not exhibit clear synchronous oscillations (Figure 5A). GPe neurons showed sustained, high-frequency discharges with occasional pauses, and the STN neuron emitted isolated spikes in an irregular manner. In spite of high firing rates of presynaptic GPe neurons, their GABAergic influences on postsynaptic STN neurons were weak (Figure 5A) due to short-term depression of GPe-STN synapses (Figure 5A).
Figure 5. (A,B) Interactions between an STN neuron and its afferent and efferent GPe neurons in the normal (A) and parkinsonian states (B). The top panel shows a raster plot of an STN neuron (red) and afferent and efferent GPe neurons (blue). The middle three panels show membrane potential, GABA-receptor-mediated synaptic conductance and an inactivation variable of T-current as a function of time for the STN neuron. The bottom panel shows the transmission probability of a unitary GPe projection to the STN neuron. A gray dotted trace, and light gray area indicate the mean and mean ± standard deviation of transmission probabilities of all GPe projections to the STN neuron. A blue trace indicates a transmission probability of the projection from the afferent GPe neuron as a function of time.
In the parkinsonian state (Figure 5B), STN neurons received strong inhibition due to burst discharges of their presynaptic GPe neurons, and were rapidly hyperpolarized so that T-currents of STN neurons were de-inactivated, as shown by an increase of the T-current inactivation variable (Figure 5B). This induced strong post-inhibitory rebound excitation of STN neurons. Consequent burst firing of STN neurons then strongly excited their postsynaptic GPe neurons to discharge in bursts. Timing of burst discharges of efferent GPe neurons overlapped with subsequent bursts of afferent GPe neurons (Figures 4B, 5B). Although an afferent GPe neuron seldom received a reciprocal excitatory connection from its postsynaptic STN neurons under the sparse random connection architecture, these synchronous bursts of GPe neurons shut off rebound burst firing of STN neurons and then hyperpolarized them again. By repeating this cycle of events, STN-GPe circuits maintained rhythmic burst discharges.
The timing of shut-off of STN burst discharges by inhibition from GPe neurons is critical for the 14–15 Hz oscillations because post-inhibitory rebound burst discharges could last longer without active shut-offs. The slow time scale of the post-inhibitory rebound mechanism alone could lead to much slower oscillations. Although almost none of the pairs of STN and GPe neurons in the network are reciprocally connected under the connection probabilities between STN and GPe neurons in the network model, network synchronization makes presynaptic and postsynaptic GPe neurons to an STN neuron discharge in bursts in a synchronous manner.
In order to assess the role of this sparsely closed circuit of STN and GPe neurons for the generation of oscillatory burst discharges, we ran simulations of network models with different proportions of reciprocally-connected pairs of STN and GPe neurons, while the connection probability was kept constant (Figure 6A). The power of 8–15 Hz oscillations increased with an increasing proportion of reciprocally connected neurons. However, oscillations remained even when reciprocal connections were absent.
Figure 6. Effects of probability of closed-loop connections between STN and GPe neurons (A) sparseness in connections between STN and GPe neurons (B) properties of short-term depression at GABAergic synapses of STN neurons (C) and stochasticity in GPe intrinsic noisy conductances (D) on the 8–15 Hz oscillatory activity in the STN. Each panel represents dependencies of STN parkinsonian oscillations over maximal conductances of GPe-STN synapses and the GPe Int. in the same manner as in Figure 3 with a modified parameter value noted above. Grey circles and squares indicate parameters that correspond to the parkinsonian and normal states, respectively. (A) Probability of closed-loop connections between the STN and GPe neurons increased from near-zero probability of closed-loop connections as a result of the default connection rule, random-convergent connections, denoted as “rand.” (B) Sparseness in connections between the STN and GPe neurons was enhanced by scaling network size 2, 4, 8, and 16 times under the same connection rule. (C), Strength of short-term depression at GABAergic synapses of STN neurons was enhanced by reducing values of the short-term depression parameter, nbound. (D), Increasing the mean rate of a Poisson spike train noted above the panels with decreasing maximal conductance so as to keep the mean conductance constant over time, reduced the variation of GPe intrinsic noisy conductances. See details in the main text.
We also tested whether synchronized bursts are due to highly convergent connections between STN and GPe neurons in the small network model we used. In simulations of network models with reduced connection probabilities between STN and GPe neurons by increasing the total numbers of STN and GPe neurons while maintaining the numbers of target neurons for each STN and GPe neuron constant, we observed stronger oscillatory burst activity (Figure 6B). This indicates that generation of oscillatory burst firings is not the result of the dense connectivity between STN and GPe due to the small network size.
STN-GPe network models with higher levels of GPe intrinsic excitability, including those in the normal state, exhibited much weaker oscillations in the 8–15 Hz range. Greater GPe intrinsic excitability can reduce GABAergic conductance of STN neurons due to the strong short-term depression of GPe-STN synapses (Figure 5B). The time constant of short-term depression was long enough to keep the synaptic transmission probability at STN GABAergic synapses almost fully depressed, indicating that short-term depression can act as a switch that turns on and off the 8–15 Hz burst oscillations in response to the level of presynaptic GPe activity. We tested the sensitivity of 8–15 Hz oscillatory burst discharges to the short-term depression parameter, nbound (Figure 6C). Smaller nbound, which causes stronger depression, reduced the 8–15 Hz oscillations with lower levels of GPe intrinsic excitability, and larger nbound enhanced the oscillations.
Another possible effect of higher GPe intrinsic excitability is enhanced stochasticity in GPe firing, because we implemented that using noisy conductance that obeyed a Poisson process (see Method). To address this issue, we investigated the effect of varying GPe intrinsic noisy conductance by increasing the mean rate of the Poisson process and reducing the maximal conductance inversely so as to keep the mean value constant. With decreased variance of GPe noisy intrinsic conductance, a slight enhancement of 8–15 Hz oscillations occurred at intermediate GPe intrinsic conductances (0.75 nS), but very little enhancement was observed at higher levels (1.5 nS), even when the conductance did not fluctuate stochastically, but remained constant at the mean value (Figure 6D). This result suggests that stochasticity in GPe intrinsic conductance is not a major factor in preventing the occurrence of 8–15 Hz oscillations.
Roles of Cortical Excitatory Inputs on STN Parkinsonian Oscillations
Next, we investigated how cortical excitatory inputs to the STN affect parkinsonian oscillations. Anatomically, layer-V pyramidal-tract neurons in the motor cortex are a major excitatory source for the STN (Parent and Parent, 2006; Kita and Kita, 2012). Layer-V pyramidal-tract neurons exhibit spontaneous firing and phasically increase their activities just before and during movements (Turner and DeLong, 2000; Pasquereau and Turner, 2011). First, we investigated the effects of cortical tonic background excitatory inputs to STN neurons on STN parkinsonian oscillations. We used homogeneous Poisson spike trains of fixed mean rate (100 Hz) as the background excitatory inputs to STN neurons, which were not correlated with BG parkinsonian oscillations. In contrast to experimental evidence (Tachibana et al., 2011), reducing synaptic strengths of background inputs to the STN enhanced 8–15 Hz neuronal oscillations in the STN (Figure 7A) and GPe (data not shown). This is because tonic excitatory inputs to STN neurons antagonized GPe inhibition, which depolarized STN neurons and diminished post-inhibitory rebound excitation. The suppressive effect of the background excitatory inputs on the parkinsonian oscillations suggests that stronger cortical excitatory inputs can reduce BG parkinsonian oscillations. This is consistent with experimental reports showing that pathological oscillations in LFPs recorded from the STN of PD patients were reduced by voluntary movements (Amirnovin et al., 2004; Brown and Williams, 2005) because just before and during voluntary movement, cortical inputs to the STN are phasically enhanced (Turner and DeLong, 2000).
Figure 7. Effects of cortical excitatory inputs on the STN 8-15 Hz oscillatory activity. (A) Reduced background excitatory inputs to STN neurons in the model increase the peak power of STN 8–15 Hz neuronal oscillations. Each panel represents dependencies of the STN parkinsonian oscillations over the maximal conductances of GPe-STN synapses and GPe Int. in the same manner as in Figure 3. The maximal conductance of an AMPA-mediated synapse of an STN neuron that receives background excitatory inputs from the cortex, ḡAMPA−Ctx−STN, noted above each panel. (B,C) Phasically increased excitatory inputs to STN neurons suppressed 8–15 Hz oscillations of STN and GPe neurons. (B) Typical spectrogram of an STN neuron in the model around 1,000 Hz excitatory Poisson spike inputs (PT-stimulation) mimicking peri-movement activations of cortical, layer-V pyramidal tract neurons. (C) Changes of mean power of 8–15 Hz oscillations of the STN and GPe neurons of the model by 1,000 Hz PT-stimulations. Blue, red, and green columns represent normalized power of 8–15 Hz oscillations of 1-s period before, just after, and 1-s after onset of PT-stimulation of 1,000 Hz, respectively. * indicates statistically significant changes (Friedman rank sum test, p < 0.001; pair-wise Wilcoxon rank sum test, p < 0.001, n = 256 neurons). (D) Correlation of changes of T-current conductances of and normalized power of oscillations of STN neurons by PT-stimulation. Each circle represents mean values of a neuron at a 250-ms bin during 1-s after onset of PT-stimulation. (E) Peak power of 8-15 Hz oscillations in the STN neurons depended on different GPi-thalamo-cortical delays of cortico-STN feedback oscillatory inputs. Each statistical value was calculated from pooled spike trains of STN neurons in 8 independent simulations (total n = 512 per condition). Red hatched and blue bars correspond to network parameters for the parkinsonian and normal conditions shown in Figure 3, respectively, except as otherwise noted. “delay” indicates the GPi-thalamus-cortex delays of AMPA-mediated recurrent feedback oscillatory inputs to the STN. Presence and absence of tonic background (bkg) and recurrent (rec) cortical excitatory inputs are shown as o and -, respectively. The leftmost bar represents data from the parkinsonian state without cortical recurrent inputs to the STN, shown in Figure 3, and the second bar from the left represents data of the parkinsonian state without either the tonic background or recurrent cortical inputs to the STN. Asterisks, hash marks, and plus marks represent significant differences (pair-wise t-test with Bonferroni adjustment, p < 0.001, n = 1,024) from the parkinsonian state without AMPA-mediated, recurrent oscillatory inputs to the STN, from the parkinsonian state without either the AMPA-mediated, background or recurrent oscillatory inputs to the STN, and from the parkinsonian state with AMPA-mediated, background inputs to the STN, and recurrent oscillatory inputs to the STN with a delay of 70 ms, respectively.
Simulation results indeed showed that parkinsonian neuronal oscillations dropped at about the time of phasic increases (up to 1000 Hz) of excitatory inputs to STN neurons (Figures 7B,C). Reduction of STN parkinsonian neuronal oscillations correlated with reduction of STN T-current conductances that underlie post-inhibitory rebound excitation (r2 = 0.26, p < 0.001; Figure 7D). Smaller phasic increases (up to 250 Hz or 500 Hz) of excitatory inputs to STN neurons also showed drops of STN oscillations (Friedman test, p < 0.001; pairwise comparisons using Wilcoxon rank sum test, p < 0.001, data not shown).
We also investigated effects of oscillatory cortical inputs to STN neurons on STN parkinsonian oscillations. Because the manner in which cortical oscillations relate to STN parkinsonian oscillations is still unclear (Brittain and Brown, 2014), we simply assumed that STN neuronal oscillations propagate through the STN-GPi-thalamus-cortex loop and go back to the STN, while maintaining their oscillation frequency. We investigated the effect of the time delay in the STN-GPi-thalamus-cortex feedback loop on STN 8–15 Hz oscillations by introducing delay lines via which STN neurons were recurrently connected. More precisely, each STN neuron was randomly connected with 16 afferent STN neurons (average) and spikes were stochastically transmitted (p = 0.25) via delay lines. Synaptic strength of delayed feedback inputs was the same as that of background inputs to the STN. Simulation results showed that STN 8–15 Hz neuronal oscillations were affected by both the delay and the strength of the recurrent excitatory input to STN neurons (Welch's one-way ANOVA F = 4741.3, df = 17, p < 0.001). When the delay was within 60–80 ms, which was close to the cycle of STN parkinsonian oscillations, the STN oscillation in the parkinsonian state was stronger with feedback recurrent inputs than without (Figure 7E, pair-wise t-test with Bonferroni adjustment, p < 0.001, n = 512). In contrast, when the delay was shorter than 60 ms or longer than 80 ms, feedback recurrent inputs diminished STN oscillations to less than those without feedback. The power of STN oscillations in the parkinsonian state with feedback recurrent inputs, and with a delay between 60 and 80 ms was even greater than with a blockade of AMPA transmission both at background and at recurrent inputs to the STN (Figure 7E, pair-wise t-test with Bonferroni adjustment, p < 0.001, n = 512). Feedback excitatory inputs did not affect STN 8–15 Hz oscillations in the normal state (Figure 7E, pair-wise t-test with Bonferroni adjustment, p = 1.0, n = 512). These results indicate that cortical excitatory inputs affect STN parkinsonian oscillations in a complex manner: cortical non-oscillatory inputs are suppressive to STN parkinsonian oscillations, but feedback oscillatory inputs can amplify or diminish STN parkinsonian oscillations in a delay-dependent manner. The complex effects of cortical excitatory inputs on STN parkinsonian oscillations explain two lines of experimental evidence suggesting suppressive or facilitative effects of cortical-STN excitatory inputs, respectively (Amirnovin et al., 2004; Brown and Williams, 2005; Tachibana et al., 2011).
In this study, we explored mechanisms of generation of parkinsonian 8–15 Hz neuronal burst oscillations in the STN-GPe circuit by building a spiking neural network model, incorporating anatomical and physiological constraints. Through a systematic search of the parameter space, we identified parameter regions that reproduce normal and parkinsonian states as reported by Tachibana et al. (2011). The parkinsonian state was characterized by reduced intrinsic excitability of GPe neurons and increased conductance of GPe-STN synapses. Reduction of autonomous activity of GPe neurons by dopamine depletion paradoxically elevated inhibition of STN neurons via strong short-term depression and consequently induced strong post-inhibitory rebound burst firing of STN neurons, which were then shaped into synchronous oscillations in the 14–15 Hz frequency range via interactions between STN and GPe neurons.
Characteristics of neuronal activities observed in our simulations agreed well with experimental findings. First, generation of neuronal burst oscillations is critically dependent on GABAergic inhibition of STN neurons (Tachibana et al., 2011) and is sensitive to a level of GABAergic inhibition of STN neurons (Figure 3). Second, the frequency range of these neuronal burst oscillations is in the 8–15 Hz frequency band (Figure 4B), in which pathological oscillations were observed in the BG of dopamine-depleted monkeys (Heimer et al., 2006; Wichmann and Soares, 2006; Tachibana et al., 2011). Third, mean firing rates of neurons in the model exhibited in simulations under normal or dopamine-depleted states (Figure 3) accord well with those experimentally observed in the BG of monkeys (Tachibana et al., 2011). Fourth, inactivation of the STN diminishes 8–15 Hz neuronal oscillations in the GPe and GPi (Tachibana et al., 2011), which is consistent with our simulations with smaller synaptic strengths of STN projections (Figure 6). Fifth, acute administration of the dopamine precursor, L-DOPA, reduces the power of 8–15 Hz oscillations, the strength of burst discharges in the STN and GPi, and firing rates in the STN (Tachibana et al., 2011). These experimental results are consistent with our simulation results (Figures 3, 4) if L-DOPA reduces GABAergic inhibition of the STN by restoring D2-receptor-mediated presynaptic suppression (Shen and Johnson, 2000, 2005), but does not acutely elevate GPe autonomous activity through dopamine-dependent expression of HCN channels, which is supposed to be a slower process (Chan et al., 2011). Sixth, burst frequency in the STN and GPe (but not in the GPi) in the model were higher in the pathological state than that in the normal state (Figures 3, 4), which is consistent with experiments (Wichmann and Soares, 2006; Tachibana et al., 2011). Seventh, our simulation results suggest that contributions of cortical excitatory inputs to the STN for PD burst oscillations are complex and depend on their temporal patterns and firing rates (Figure 7). Enhancing non-oscillatory cortical inputs suppressed STN PD oscillations, which is consistent with experimental evidence from PD patients (Amirnovin et al., 2004; Brown and Williams, 2005), although there are also experimental results showing that the blockade of ionotropic glutamatergic transmission in the STN diminished STN PD oscillations (Tachibana et al., 2011). Additional oscillatory inputs, the frequency of which was similar to STN PD oscillations, could amplify or diminish STN PD oscillations in a phase-dependent manner. These amplifying effects of oscillatory cortical inputs are consistent with experimental results (Tachibana et al., 2011).
We verified the critical role of strong post-inhibitory rebound excitation of an STN neuron for generation and maintenance of PD 8–15 Hz neuronal burst oscillations (Figure 5). The same mechanism was already explored in the pioneering study of a spiking neuron model for PD oscillations (Terman et al., 2002), which explained the slower oscillations, but not the 8–15 Hz oscillations. The slowly-changing nature of the rebound mechanism is ascribed to the slower oscillations (Holgado et al., 2010). In our model, in response to STN rebound firing, GPe neurons actively terminate the rebound excitation of the STN despite sparsely closed loops between STN and GPe neurons, where each population was synchronously activated. This network mechanism enabled faster neuronal oscillations in spite of the slow rebound mechanism. The suppressive role of excitatory inputs on the STN for the 8–15 Hz oscillations we observed can be explained as the damping of the post-inhibitory rebound excitation by preventing required hyperpolarization. It is consistent with peri-movement, phasic activation of layer-V pyramidal tract neurons in the motor cortex projecting to both the spinal cord and the STN (Turner and DeLong, 2000). It is also consistent with the therapeutic effects of deep brain stimulation to the STN that is able to suppress 8–15 Hz oscillations by virtue of its high-frequency, direct stimulation of cortical afferents to the STN (Miocinovic et al., 2006).
The decline of GPe autonomous activity due to down regulation of the HCN channel caused by dopamine depletion (Chan et al., 2011) is the key phenomenon in our model for converting normal neuronal activity into a pathological state. Normal autonomous activity of GPe neurons could have two possible effects: maintenance of a high level of neuronal firing, and decorrelation of GPe neuronal activities by collateral inhibition (Edgerton and Jaeger, 2011). Our simulations showed that GPe neuronal firing maintained at a high level reduces GABAergic inhibition of the STN via strong short-term depression (Atherton et al., 2013) and prevents the generation of pathological oscillations. Once GPe neuronal firing is reduced due to reduced autonomous activity caused by dopamine depletion, GABAergic inhibition becomes less depressed and induces the rebound excitation of STN neurons, and consequently pathological oscillations occur. Decorrelated GPe neuron activities may prevent the generation or maintenance of 8–15 Hz burst oscillations, but it showed limited effects in a PD state in our simulations. Decorrelation of GPe neurons might be more implicated in normal BG functioning.
Although elevated striatal inhibition of GPe neurons by dopamine depletion (Liang et al., 2008) has been considered the key cause of PD pathologies (Albin et al., 1989; DeLong and Wichmann, 2007), we did not explore its role relative to pathological oscillations because of negative experimental evidence showing that blockade of ionotropic GABAergic transmission in the GPe does not change or elevate the power of the pathological oscillations (Tachibana et al., 2011). However, elevated striatal inhibition could have a similar effect, since the decline of GPe autonomous activity in our model reduced GPe firing to cause pathological oscillations.
Our simulation results showed that in-phase oscillatory excitatory inputs from the cortex amplify pathological oscillations in the STN, as observed experimentally (Tachibana et al., 2011). Indeed, cortical and STN oscillatory activities are mainly coherent in the low (12–20 Hz) and high (20–35 Hz) beta range (Hirschmann et al., 2011). Whether STN PD lower-beta oscillations proliferate throughout the cortex-BG-thalamus loops as an entrainment loop is currently a matter of debate (Brittain and Brown, 2014). However, enhancement of cortical lower-beta oscillations in PD patients (Pollok et al., 2012) and neuronal lower-beta (and upper-beta) oscillations of pyramidal tract neurons in the motor cortex of dopamine-depleted primates (Pasquereau and Turner, 2011) are consistent with the entrainment loop hypothesis. The required delay from the STN to the cortex might be too large considering just conduction and synaptic delays. However, because GPi projections to the thalamus are inhibitory, oscillations could be delayed by more than the sum of conduction and synaptic delays, as observed in closed STN-GPe circuits. The reduction of STN PD lower-beta oscillations by out-of-phase cortical oscillations in our simulations is a result of damping of STN post-inhibitory rebound excitation by counteraction of phasically-elevated excitatory inputs from the motor cortex to inhibition from the GPe. It is similar to the suppression by dynamically increased cortical inputs during movements (Turner and DeLong, 2000), and account for mechanisms of closed-loop deep brain stimulation (Rosin et al., 2011).
Our model has several limitations in reproducing experimentally observed phenomena. First, our model did not show neuronal oscillations in the 3–8 Hz frequency band in the BG of dopamine-depleted monkeys (Heimer et al., 2006; Wichmann and Soares, 2006; Tachibana et al., 2011), although the relationship of such oscillations to dopamine depletion is controversial. Second, burst frequency of the STN in the parkinsonian state in our simulations was lower than in experimental observations (Wichmann and Soares, 2006; Tachibana et al., 2011).
In summary, our model provides plausible mechanisms for generation and maintenance of PD oscillations in the BG. Our model also showed that cortical inputs to the STN can amplify PD oscillations and it predicted required conditions for amplification or reduction of pathological BG oscillations.
OS and KD planned the study. YT and AN provided physiological data. OS constructed the model and performed numerical simulation. All authors discussed the results and OS and KD wrote the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by HPCI Strategic Program for Supercomputational Life Science and Post-K Application Development for Exploratory Challenges from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT) of Japan, and Okinawa Institute of Science and Technology Graduate University research support to KD, and Grant-in-Aid for Scientific Research on Innovative Areas, “Non-linear Neuro-oscillology” (15H05873) and Grant-in-Aid for Scientific Research (A) (26250009) to AN from MEXT. We thank J. Igarashi, J. Moren, J. Yoshimoto and B. Girard for helpful comments and discussions.
Amirnovin, R., Williams, Z., Cosgrove, R., and Eskandar, E. (2004). Visually guided movements suppress subthalamic oscillations in parkinson's disease patients. J. Neurosci. 24, 11302–11306. doi: 10.1523/JNEUROSCI.3242-04.2004
Atherton, J., Menard, A., Urbain, N., and Bevan, M. (2013). Short-term depression of external globus pallidus-subthalamic nucleus synaptic transmission and implications for patterning subthalamic activity. J. Neurosci. 33, 7130–7144. doi: 10.1523/JNEUROSCI.3576-12.2013
Baufreton, J., Kirkham, E., Atherton, J., Menard, A., Magill, P., Bolam, P., et al. (2009). Sparse but selective and potent synaptic transmission from the globus pallidus to the subthalamic nucleus. J. Neurophysiol. 102, 532–545. doi: 10.1152/jn.00305.2009
Bevan, M. D., Booth, P. A., Eaton, S. A., and Bolam, J. P. (1998). Selective innervation of neostriatal interneurons by a subclass of neuron in the globus pallidus of the rat. J. Neurosci. 18, 9438–9452.
Bevan, M. D., Magill, P. J., Terman, D., Bolam, J. P., and Wilson, C. J. (2002). Move to the rhythm: oscillations in the subthalamic nucleus-external globus pallidus network. Trends Neurosci. 25, 525–531. doi: 10.1016/S0166-2236(02)02235-X
Boraud, T., Bezard, E., Bioulac, B., Gross, C. E. (2002). From single extracellular unit recording in experimental and human parkinsonism to the development of a functional concept of the role played by the basal ganglia in motor control. Prog. Neurobiol. 66, 265–283.
Brown, P., and Williams, D. (2005). Basal ganglia local field potential activity: character and functional significance in the human. Clin. Neurophysiol. 116, 2510–2519. doi: 10.1016/j.clinph.2005.05.009
Chan, S., Glajch, K., Gertler, T., Guzman, J., Mercer, J., Lewis, A., et al. (2011). HCN channelopathy in external globus pallidus neurons in models of Parkinson's disease. Nat. Neurosci. 14, 85–92. doi: 10.1038/nn.2692
Edgerton, J., and Jaeger, D. (2011). Dendritic sodium channels promote active decorrelation and reduce phase locking to parkinsonian input oscillations in model globus pallidus neurons. J. Neurosci. 31, 10919–10936. doi: 10.1523/JNEUROSCI.6062-10.2011
Fan, K., Baufreton, J., Surmeier, J., Chan, S., and Bevan, M. (2012). Proliferation of external globus pallidus-subthalamic nucleus synapses following degeneration of midbrain dopamine neurons. J. Neurosci. 32, 13718–13728. doi: 10.1523/JNEUROSCI.5750-11.2012
Gillies, A., Willshaw, D., and Li, Z. (2002). Subthalamic-pallidal interactions are critical indetermining normal and abnormal functioning of the basal ganglia. Proc. Biol. Sci. 269, 545–551. doi: 10.1098/rspb.2001.1817
Günay, C., Edgerton, J., and Jaeger, D. (2008). Channel density distributions explain spiking variability in the globus pallidus: a combined physiology and computer simulation database approach. J. Neurosci. 28, 7476–7491. doi: 10.1523/JNEUROSCI.4198-07.2008
Hahn, P. J., and McIntyre, C. C. (2010). Modeling shifts in the rate and pattern of subthalamopallidal network activity during deep brain stimulation. J. Comput. Neurosci. 28, 425–441. doi: 10.1007/s10827-010-0225-8
Heimer, G., Rivlin-Etzion, M., Bar-Gad, I., Goldberg, J., Haber, S., and Bergman, H. (2006). Dopamine replacement therapy does not restore the full spectrum of normal pallidal activity in the 1-Methyl-4-Phenyl-1,2,3,6-Tetra-Hydropyridine primate model of parkinsonism. J. Neurosci. 26, 8101–8114. doi: 10.1523/JNEUROSCI.5140-05.2006
Hernandez-Lopez, S., Tkatch, T., Perez-Garci, E., Galarraga, E., Bargas, J., Hamm, H., et al. (2000). D2 dopamine receptors in striatal medium spiny neurons reduce L-type Ca2+ currents and excitability via a novel PLCβ1-IP3–calcineurin-signaling cascade. J. Neurosci. 20, 8987–8995.
Hirschmann, J., Özkurt, T. E., Butz, M., and Homburger, M. (2011). Distinct oscillatory STN-cortical loops revealed by simultaneous MEG and local field potential recordings in patients with Parkinson's disease. Neuroimage 55, 1159–1168. doi: 10.1016/j.neuroimage.2010.11.063
Holgado, A., Terry, J., and Bogacz, R. (2010). Conditions for the generation of beta oscillations in the subthalamic nucleus-globus pallidus network. J. Neurosci. 30, 12340–12352. doi: 10.1523/JNEUROSCI.0817-10.2010
Kita, H., Tachibana, Y., Nambu, A., and Chiken, S. (2005). Balance of monosynaptic excitatory and disynaptic inhibitory responses of the globus pallidus induced after stimulation of the subthalamic nucleus in the monkey. J. Neurosci. 25, 8611–8619. doi: 10.1523/JNEUROSCI.1719-05.2005
Kita, T., and Kita, H. (2012). The subthalamic nucleus is one of multiple innervation sites for long-range corticofugal axons: a single-axon tracing study in the rat. J. Neurosci. 32, 5990–5999. doi: 10.1523/JNEUROSCI.5717-11.2012
Kühn, A., Trottenberg, T., Kivi, A., Kupsch, A., Schneider, G.-H., and Brown, P. (2005). The relationship between local field potential and neuronal discharge in the subthalamic nucleus of patients with Parkinson's disease. Exp. Neurol. 194, 212–220. doi: 10.1016/j.expneurol.2005.02.010
Kumar, A., Cardanobile, S., Rotter, S., and Aertsen, A. (2011). The role of inhibition in generating and controlling parkinson's disease oscillations in the basal ganglia. Front. Syst. Neurosci. 5:86. doi: 10.3389/fnsys.2011.00086
Leblois, A., Boraud, T., Meissner, W., Bergman, H., and Hansel, D. (2006a). Competition between feedback loops underlies normal and pathological dynamics in the basal ganglia. J. Neurosci. 26, 3567–3583. doi: 10.1523/JNEUROSCI.5050-05.2006
Leblois, A., Meissner, W., Bezard, E., Bioulac, B., Gross, C., and Boraud, T. (2006b). Temporal and spatial alterations in GPi neuronal encoding might contribute to slow down movement in Parkinsonian monkeys. Eur. J. Neurosci. 24, 1201–1208. doi: 10.1111/j.1460-9568.2006.04984.x
Levy, R., Ashby, P., Hutchison, W., Lang, A., Lozano, A., and Dostrovsky, J. (2002). Dependence of subthalamic nucleus oscillations on movement and dopamine in Parkinson's disease. Brain 125, 1196–1209. doi: 10.1093/brain/awf128
Liang, L., DeLong, M. R., and Papa, S. M. (2008). Inversion of dopamine responses in striatal medium spiny neurons and involuntary movements. J. Neurosci. 28, 7537–7547. doi: 10.1523/JNEUROSCI.1176-08.2008
Mallet, N., Pogosyan, A., Márton, L., Bolam, P., Brown, P., and Magill, P. (2008). Parkinsonian beta oscillations in the external globus pallidus and their relationship with subthalamic nucleus activity. J. Neurosci. 28, 14245–14258. doi: 10.1523/JNEUROSCI.4199-08.2008
Miller, W. C., and DeLong, M. R. (1987). “Altered tonic activity of neurons in the globus pallidus and subthalamic nucleus in the primate MPTP model of parkinsonism,” in The Basal Ganglia II, eds M. B. Carpenter and A. Jayraman (New York, NY: Plenum), 415–427.
Miocinovic, S., Parent, M., Butson, C., Hahn, P., Russo, G., Vitek, J., et al. (2006). Computational analysis of subthalamic nucleus and lenticular fasciculus activation during therapeutic deep brain stimulation. J. Neurophysiol. 96, 1569–1580. doi: 10.1152/jn.00305.2006
Nambu, A., Tokuno, H., Hamada, I., Kita, H., Imanishi, M., Akazawa, T., et al. (2000). Excitatory cortical inputs to pallidal neurons via the subthalamic nucleus in the monkey. J. Neurophysiol. 84, 289–300.
Otsuka, T., Abe, T., Tsukagawa, T., and Song, W.-J. (2004). Conductance-based model of the voltage-dependent generation of a plateau potential in subthalamic neurons. J. Neurophysiol. 92, 255–264. doi: 10.1152/jn.00508.2003
Pasquereau, B., and Turner, R. (2011). Primary motor cortex of the parkinsonian monkey: differential effects on the spontaneous activity of pyramidal tract-type neurons. Cereb. Cortex 21, 1362–1378. doi: 10.1093/cercor/bhq217
Pollok, B., Krause, V., Martsch, W., Wach, C., Schnitzler, A., and Südmeyer, M. (2012). Motor-cortical oscillations in early stages of Parkinson's disease. J. Physiol. 590, 3203–3212. doi: 10.1113/jphysiol.2012.231316
Priori, A., Foffani, G., Pesenti, A., Tamma, F., Bianchi, A. M., Pellegrini, M., et al. (2004). Rhythm-specific pharmacological modulation of subthalamic activity in Parkinson's disease. Exp. Neurol. 189, 369–379. doi: 10.1016/j.expneurol.2004.06.001
Raz, A., Vaadia, E., and Bergman, H. (2000). Firing patterns and correlations of spontaneous discharge of pallidal neurons in the normal and the tremulous 1-methyl-4-phenyl-1, 2, 3, 6-tetrahydropyridine vervet model of Parkinsonism. J. Neurosci. 20, 8559–8571.
Rivlin-Etzion, M., Ritov, Y., Heimer, G., Bergman, H., and Bar-Gad, I. (2006). Local shuffling of spike trains boosts the accuracy of spike train spectral analysis. J. Neurophysiol. 95, 3245–3256. doi: 10.1152/jn.00055.2005
Rosin, B., Slovik, M., Mitelman, R., Rivlin-Etzion, M., Haber, S., Israel, Z., et al. (2011). Closed-loop deep brain stimulation is superior in ameliorating Parkinsonism. Neuron 72, 370–384. doi: 10.1016/j.neuron.2011.08.023
Rubin, J. E., Terman, D. (2004). High frequency stimulation of the subthalamic nucleus eliminates pathological thalamic rhythmicity in a computational model. J. Comput. Neurosci. 16, 211–235. doi: 10.1023/B:JCNS.0000025686.47117.67
Sato, F., Lavallee, P., Levesque, M., and Parent, A. (2000). Single-axon tracing study of neurons of the external segment of the globus pallidus in primate. J. Comp. Neurol. 417, 17–31. doi: 10.1002/(SICI)1096-9861(20000131)417:1<17::AID-CNE2>3.0.CO;2-I
Shen, K., and Johnson, S. (2000). Presynaptic dopamine D2 and muscarine M3 receptors inhibit excitatory and inhibitory transmission to rat subthalamic neurones in vitro. J. Physiol. 525, 331–341. doi: 10.1111/j.1469-7793.2000.00331.x
Tachibana, Y., Iwamuro, H., Kita, H., Takada, M., and Nambu, A. (2011). Subthalamo-pallidal interactions underlying parkinsonian neuronal oscillations in the primate basal ganglia. Eur. J. Neurosci. 34, 1470–1484. doi: 10.1111/j.1460-9568.2011.07865.x
Wichmann, T., and Soares, J. (2006). Neuronal firing before and after burst discharges in the monkey basal ganglia is predictably patterned in the normal state and altered in Parkinsonism. J. Neurophysiol. 95, 2120–2133. doi: 10.1152/jn.01013.2005
Keywords: basal ganglia, parkinsonian oscillation, subthalamic nucleus, globus pallidus, beta oscillation, rebound excitaton, short-term plasticity
Citation: Shouno O, Tachibana Y, Nambu A and Doya K (2017) Computational Model of Recurrent Subthalamo-Pallidal Circuit for Generation of Parkinsonian Oscillations. Front. Neuroanat. 11:21. doi: 10.3389/fnana.2017.00021
Received: 03 October 2016; Accepted: 06 March 2017;
Published: 21 March 2017.
Edited by:Fumino Fujiyama, Doshisha University, Japan
Reviewed by:Michela Ferrucci, University of Pisa, Italy
Javier Blesa, Centro Integral de Neurociencias, Spain
Yasumasa Ueda, Kansai Medical University, Japan
Copyright © 2017 Shouno, Tachibana, Nambu and Doya. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kenji Doya, firstname.lastname@example.org
†Present Address: Yoshihisa Tachibana, Division of System Neuroscience, Kobe University Graduate School of Medicine, Kobe, Japan