Improved conditions for the generation of beta oscillations in the subthalamic nucleus--globus pallidus network

Mathematical derivation and mechanism analysis of beta oscillations in a cortex-pallidum model

In this paper, we develop a new cortex-pallidum model to study the origin mechanism of Parkinson's oscillations in the cortex. In contrast to many previous models, the globus pallidus internal (GPi) and externa (GPe) both exert direct inhibitory feedback to the cortex. Using Hopf bifurcation analysis, two new critical conditions for oscillations, which can include the self-feedback projection of GPe, are obtained. In this paper, we find that the average discharge rate (ADR) is an important marker of oscillations, which can divide Hopf bifurcations into two types that can uniformly be used to explain the oscillation mechanism. Interestingly, the ADR of the cortex first increases and then decreases with increasing coupling weights that are projected to the GPe. Regarding the Hopf bifurcation critical conditions, the quantitative relationship between the inhibitory projection and excitatory projection to the GPe is monotonically increasing; in contrast, the relationship between different coupling weights in the cortex is monotonically decreasing. In general, the oscillation amplitude is the lowest near the bifurcation points and reaches the maximum value with the evolution of oscillations. The GPe is an effective target for deep brain stimulation to alleviate oscillations in the cortex.

Criticality supports cross-frequency cortical-thalamic information transfer during conscious states

Consciousness is thought to be regulated by bidirectional information transfer between the cortex and thalamus, but the nature of this bidirectional communication - and its possible disruption in unconsciousness - remains poorly understood. Here, we present two main findings elucidating mechanisms of corticothalamic information transfer during conscious states. First, we identify a highly preserved spectral channel of cortical-thalamic communication that is present during conscious states, but which is diminished during the loss of consciousness and enhanced during psychedelic states. Specifically, we show that in humans, mice, and rats, information sent from either the cortex or thalamus via δ/θ/α waves (∼1-13 Hz) is consistently encoded by the other brain region by high γ waves (52-104 Hz); moreover, unconsciousness induced by propofol anesthesia or generalized spike-and-wave seizures diminishes this cross-frequency communication, whereas the psychedelic 5-methoxy-N,N-dimethyltryptamine (5-MeO-DMT) enhances this low-to-high frequency interregional communication. Second, we leverage numerical simulations and neural electrophysiology recordings from the thalamus and cortex of human patients, rats, and mice to show that these changes in cross-frequency cortical-thalamic information transfer may be mediated by excursions of low-frequency thalamocortical electrodynamics toward/away from edge-of-chaos criticality, or the phase transition from stability to chaos. Overall, our findings link thalamic-cortical communication to consciousness, and further offer a novel, mathematically well-defined framework to explain the disruption to thalamic-cortical information transfer during unconscious states.

Multivariable closed-loop control of deep brain stimulation for Parkinson's disease

Objective. Closed-loop deep brain stimulation (DBS) methods for Parkinson's disease (PD) to-date modulate either stimulation amplitude or frequency to control a single biomarker. While good performance has been demonstrated for symptoms that are correlated with the chosen biomarker, suboptimal regulation can occur for uncorrelated symptoms or when the relationship between biomarker and symptom varies. Control of stimulation-induced side-effects is typically not considered.Approach.A multivariable control architecture is presented to selectively target suppression of either tremor or subthalamic nucleus beta band oscillations. DBS pulse amplitude and duration are modulated to maintain amplitude below a threshold and avoid stimulation of distal large diameter axons associated with stimulation-induced side effects. A supervisor selects between a bank of controllers which modulate DBS pulse amplitude to control rest tremor or beta activity depending on the level of muscle electromyographic (EMG) activity detected. A secondary controller limits pulse amplitude and modulates pulse duration to target smaller diameter axons lying close to the electrode. The control architecture was investigated in a computational model of the PD motor network which simulated the cortico-basal ganglia network, motoneuron pool, EMG and muscle force signals.Main results.Good control of both rest tremor and beta activity was observed with reduced power delivered when compared with conventional open loop stimulation, The supervisor avoided over- or under-stimulation which occurred when using a single controller tuned to one biomarker. When DBS amplitude was constrained, the secondary controller maintained the efficacy of stimulation by increasing pulse duration to compensate for reduced amplitude. Dual parameter control delivered effective control of the target biomarkers, with additional savings in the power delivered.Significance.Non-linear multivariable control can enable targeted suppression of motor symptoms for PD patients. Moreover, dual parameter control facilitates automatic regulation of the stimulation therapeutic dosage to prevent overstimulation, whilst providing additional power savings.

The Therapeutic Frequency Profile of Subthalamic Nucleus Deep Brain Stimulation in Rats Is Shaped by Antidromic Spike Failure

The therapeutic mechanisms of subthalamic nucleus (STN) deep brain stimulation (DBS) may depend on antidromic activation of cortex via the hyperdirect pathway. However, hyperdirect pathway neurons cannot reliably follow high-stimulation frequencies, and the spike failure rate appears to correlate with symptom relief as a function of stimulation frequency. We hypothesized that antidromic spike failure contributes to the cortical desynchronization caused by DBS. We measured in vivo evoked cortical activity in female Sprague Dawley rats and developed a computational model of cortical activation from STN DBS. We modeled stochastic antidromic spike failure to determine how spike failure affected the desynchronization of pathophysiological oscillatory activity in cortex. We found that high-frequency STN DBS desynchronized pathologic oscillations via the masking of intrinsic spiking through a combination of spike collision, refractoriness, and synaptic depletion. Antidromic spike failure shaped the parabolic relationship between DBS frequency and cortical desynchronization, with maximum desynchronization at ∼130 Hz. These findings reveal that antidromic spike failure plays a critical role in mediating the dependency of symptom relief on stimulation frequency.SIGNIFICANCE STATEMENT Deep brain stimulation (DBS) is a highly effective neuromodulation therapy, yet it remains uncertain why conventionally used stimulation frequencies (e.g., ∼130 Hz) are optimal. In this study, we demonstrate a potential explanation for the stimulation frequency dependency of DBS through a combination of in vivo experimental measurements and computational modeling. We show that high-frequency stimulation can desynchronize pathologic firing patterns in populations of neurons by inducing an informational lesion. However, sporadic spike failure at these high frequencies limits the efficacy of the informational lesion, yielding a parabolic profile with optimal effects at ∼130 Hz. This work provides a potential explanation for the therapeutic mechanism of DBS, and highlights the importance of considering spike failure in mechanistic models of DBS.

A model description of beta oscillations in the external globus pallidus

The external globus pallidus (GPe), a subcortical nucleus located in the indirect pathway of the basal ganglia, is widely considered to have tight associations with abnormal beta oscillations (13-30 Hz) observed in Parkinson's disease (PD). Despite that many mechanisms have been put forward to explain the emergence of these beta oscillations, however, it is still unclear the functional contributions of the GPe, especially, whether the GPe itself can generate beta oscillations. To investigate the role played by the GPe in producing beta oscillations, we employ a well described firing rate model of the GPe neural population. Through extensive simulations, we find that the transmission delay within the GPe-GPe pathway contributes significantly to inducing beta oscillations, and the impacts of the time constant and connection strength of the GPe-GPe pathway on generating beta oscillations are non-negligible. Moreover, the GPe firing patterns can be significantly modulated by the time constant and connection strength of the GPe-GPe pathway, as well as the transmission delay within the GPe-GPe pathway. Interestingly, both increasing and decreasing the transmission delay can push the GPe firing pattern from beta oscillations to other firing patterns, including oscillation and non-oscillation firing patterns. These findings suggest that if the transmission delays within the GPe are at least 9.8 ms, beta oscillations can be produced originally in the GPe neural population, which also may be the origin of PD-related beta oscillations and should be regarded as a promising target for treatments for PD.

The mechanism of Parkinson oscillation in the cortex: Possible evidence in a feedback model projecting from the globus pallidus to the cortex

The origin, location and cause of Parkinson's oscillation are not clear at present. In this paper, we establish a new cortex-basal ganglia model to study the origin mechanism of Parkinson beta oscillation. Unlike many previous models, this model includes two direct inhibitory projections from the globus pallidus external (GPe) segment to the cortex. We first obtain the critical calculation formula of Parkinson's oscillation by using the method of Quasilinear analysis. Different from previous studies, the formula obtained in this paper can include the self-feedback connection of GPe. Then, we use the bifurcation analysis method to systematically explain the influence of some key parameters on the oscillation. We find that the bifurcation principle of different cortical nuclei is different. In general, the increase of the discharge capacity of the nuclei will cause oscillation. In some special cases, the sharp reduction of the discharge rate of the nuclei will also cause oscillation. The direction of bifurcation simulation is consistent with the critical condition curve. Finally, we discuss the characteristics of oscillation amplitude. At the beginning of the oscillation, the amplitude is relatively small; with the evolution of oscillation, the amplitude will gradually strengthen. This is consistent with the experimental phenomenon. In most cases, the amplitude of cortical inhibitory nuclei (CIN) is greater than that of cortical excitatory nuclei (CEX), and the two direct inhibitory projections feedback from GPe can significantly reduce the amplitude gap between them. We calculate the main frequency of the oscillation generated in this model, which basically falls between 13 and 30 Hz, belonging to the typical beta frequency band oscillation. Some new results obtained in this paper can help to better understand the origin mechanism of Parkinson's disease and have guiding significance for the development of experiments.

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

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

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

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

A bidirectional Hopf bifurcation analysis of Parkinson's oscillation in a simplified basal ganglia model

In this paper, we study the parkinson oscillation mechanism in a computational model by bifurcation analysis and numerical simulation. Oscillatory activities can be induced by abnormal coupling weights and delays. The bidirectional Hopf bifurcation phenomena are found in simulations, which can uniformly explain the oscillation mechanism in this model. The Hopf1 represents the transition between the low firing rate stable state (SS) and oscillation state (OS), the Hopf2 represents the transition between the high firing rate stable state (HSS) and the OS, the mechanisms of them are different. The Hopf1 and Hopf2 bifurcations both show that when the state transfers from the stable region to the oscillation region, oscillatory activities always originate from the beta frequency band, and then gradually evolve into the alpha frequency band, the theta frequency band and delta frequency band in this model. We find that the changing trends of the DF and oscillation amplitude (OSAM) are contrary, oscillation activities in lower frequency band are more stable than that in higher frequency band. The effect of the delay in inhibitory pathways is greater than that of in excitatory pathways, and appropriate delays improve the discharge activation level (DAL) of the system. In all, we infer that oscillations can be induced by the follow factors: 1. Improvement of the DAL of the globus pallidus externa (GPe); 2. Reduce the DAL of the GPe from the HSS or the discharge saturation state; 3. The GPe can also resonate with the subthalamic nucleus (STN).

Average beta burst duration profiles provide a signature of dynamical changes between the ON and OFF medication states in Parkinson's disease

Parkinson's disease motor symptoms are associated with an increase in subthalamic nucleus beta band oscillatory power. However, these oscillations are phasic, and there is a growing body of evidence suggesting that beta burst duration may be of critical importance to motor symptoms. This makes insights into the dynamics of beta bursting generation valuable, in particular to refine closed-loop deep brain stimulation in Parkinson's disease. In this study, we ask the question "Can average burst duration reveal how dynamics change between the ON and OFF medication states?". Our analysis of local field potentials from the subthalamic nucleus demonstrates using linear surrogates that the system generating beta oscillations is more likely to act in a non-linear regime OFF medication and that the change in a non-linearity measure is correlated with motor impairment. In addition, we pinpoint the simplest dynamical changes that could be responsible for changes in the temporal patterning of beta oscillations between medication states by fitting to data biologically inspired models, and simpler beta envelope models. Finally, we show that the non-linearity can be directly extracted from average burst duration profiles under the assumption of constant noise in envelope models. This reveals that average burst duration profiles provide a window into burst dynamics, which may underlie the success of burst duration as a biomarker. In summary, we demonstrate a relationship between average burst duration profiles, dynamics of the system generating beta oscillations, and motor impairment, which puts us in a better position to understand the pathology and improve therapies such as deep brain stimulation.

Bifurcation analyses and potential landscapes of a cortex-basal ganglia-thalamus model

The dynamics of cortical neuronal activity plays important roles in controlling body movement and is regulated by connection weights between neurons in a cortex–basal ganglia–thalamus (BGCT) loop. Beta‐band oscillation of cortical activity is closely associated with the movement disorder of Parkinson's disease, which is caused by an imbalance in the connection weights of direct and indirect pathways in the BGCT loop. In this study, the authors investigate how the dynamics of cortical activity are modulated by connection weights of direct and indirect pathways in the BGCT loop under low dopamine levels through bifurcation analyses and potential landscapes. The results reveal that cortical activity displays rich dynamics under different connection weights, including one, two, or three stable steady states, one or two stable limit cycles, and the coexistence of one stable limit cycle with one stable steady state or two stable ones. For a low dopamine level, cortical activity exhibits oscillation for larger connection weights of direct and indirect pathways. The stability of these stable dynamics is explored by the potential landscapes.

Optimizing deep brain stimulation based on isostable amplitude in essential tremor patient models

OBJECTIVE

Deep brain stimulation is a treatment for medically refractory essential tremor. To improve the therapy, closed-loop approaches are designed to deliver stimulation according to the system's state, which is constantly monitored by recording a pathological signal associated with symptoms (e.g. brain signal or limb tremor). Since the space of possible closed-loop stimulation strategies is vast and cannot be fully explored experimentally, how to stimulate according to the state should be informed by modeling. A typical modeling goal is to design a stimulation strategy that aims to maximally reduce the Hilbert amplitude of the pathological signal in order to minimize symptoms. Isostables provide a notion of amplitude related to convergence time to the attractor, which can be beneficial in model-based control problems. However, how isostable and Hilbert amplitudes compare when optimizing the amplitude response to stimulation in models constrained by data is unknown.

APPROACH

We formulate a simple closed-loop stimulation strategy based on models previously fitted to phase-locked deep brain stimulation data from essential tremor patients. We compare the performance of this strategy in suppressing oscillatory power when based on Hilbert amplitude and when based on isostable amplitude. We also compare performance to phase-locked stimulation and open-loop high-frequency stimulation.

MAIN RESULTS

For our closed-loop phase space stimulation strategy, stimulation based on isostable amplitude is significantly more effective than stimulation based on Hilbert amplitude when amplitude field computation time is limited to minutes. Performance is similar when there are no constraints, however constraints on computation time are expected in clinical applications. Even when computation time is limited to minutes, closed-loop phase space stimulation based on isostable amplitude is advantageous compared to phase-locked stimulation, and is more efficient than high-frequency stimulation.

SIGNIFICANCE

Our results suggest a potential benefit to using isostable amplitude more broadly for model-based optimization of stimulation in neurological disorders.

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

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

Emergence of Beta Oscillations of a Resonance Model for Parkinson's Disease

In Parkinson's disease, the excess of beta oscillations in cortical-basal ganglia (BG) circuits has been correlated with normal movement suppression. In this paper, a physiologically based resonance model, generalizing an earlier model of the STN-GPe circuit, is employed to analyze critical dynamics of the occurrence of beta oscillations, which correspond to Hopf bifurcation. With the experimentally measured parameters, conditions for the occurrence of Hopf bifurcation with time delay are deduced by means of linear stability analysis, center manifold theorem, and normal form analysis. It is found that beta oscillations can be induced by increasing synaptic transmission delay. Furthermore, it is revealed that the oscillations originate from interaction among different synaptic connections. Our analytical results are consistent with the previous experimental and simulating findings, thus may provide a more systematic insight into the mechanisms underlying the transient beta bursts.

The theoretical mechanism of Parkinson's oscillation frequency bands: a computational model study

Excessive synchronous oscillation activities appear in the brain is a key pathological feature of Parkinson's disease, the mechanism of which is still unclear. Although some previous studies indicated that β oscillation (13-30 Hz) may directly originate in the network composed of the subthalamic nucleus (STN) and external globus pallidus (GPe) neurons, specific onset mechanisms of which are unclear, especially theoretical evidences in numerical simulation are still little. In this paper, we employ a STN-GPe mean-field model to explore the onset mechanism of Parkinson's oscillation. In addition to β oscillation, we find that some other common oscillation frequency bands can appear in this network, such as the α oscillation band (8-12 Hz), the θ oscillation band (4-7 Hz) and δ oscillation band (1-3 Hz). In addition to the coupling weight between the STN and GPe, the delay is also a critical factor to affect oscillatory activities, which can not be neglected compared to other parameters. Through simulation and analysis, we propose two possible theories may induce the system to transfer from the stable state to the oscillatory state in this model: (1). The oscillation activity can be induced when the firing activation level of the population increases to large enough; (2). In some special cases, the population may stay in the high firing rate stable state and the mean discharge rate of which is too large to induce oscillations, then oscillation activities may be induced as the MD decreases to moderate value. In most situations, the change trends of the MD and oscillation dominant frequency are contrary, which may be an important physiological phenomenon shown in this model. The delays and firing rates were obtained by simulating, which may be verified in the experiment in the future.

Fractional order controllers increase the robustness of closed-loop deep brain stimulation systems

We studied the effects of using fractional order proportional, integral, and derivative (PID) controllers in a closed-loop mathematical model of deep brain stimulation. The objective of the controller was to dampen oscillations from a neural network model of Parkinson's disease. We varied intrinsic parameters, such as the gain of the controller, and extrinsic variables, such as the excitability of the network. We found that in most cases, fractional order components increased the robustness of the model multi-fold to changes in the gains of the controller. Similarly, the controller could be set to a fixed set of gains and remain stable to a much larger range, than for the classical PID case, of changes in synaptic weights that otherwise would cause oscillatory activity. The increase in robustness is a consequence of the properties of fractional order derivatives that provide an intrinsic memory trace of past activity, which works as a negative feedback system. Fractional order PID controllers could provide a platform to develop stand-alone closed-loop deep brain stimulation systems.

Delta oscillations are a robust biomarker of dopamine depletion severity and motor dysfunction in awake mice

Delta oscillations (0.5-4 Hz) are a robust feature of basal ganglia pathophysiology in patients with Parkinson's disease (PD) in relationship to tremor, but their relationship to other parkinsonian symptoms has not been investigated. While delta oscillations have been observed in mouse models of PD, they have only been investigated in anesthetized animals, suggesting that the oscillations may be an anesthesia artifact and limiting the ability to relate them to motor symptoms. Here, we establish a novel approach to detect spike oscillations embedded in noise to provide the first study of delta oscillations in awake, dopamine-depleted mice. We find that approximately half of neurons in the substantia nigra pars reticulata (SNr) exhibit delta oscillations in dopamine depletion and that these oscillations are a strong indicator of dopamine loss and akinesia, outperforming measures such as changes in firing rate, irregularity, bursting, and synchrony. These oscillations are typically weakened, but not ablated, during movement. We further establish that these oscillations are caused by the loss of D2-receptor activation and do not originate from motor cortex, contrary to previous findings in anesthetized animals. Instead, SNr oscillations precede those in M1 at a 100- to 300-ms lag, and these neurons' relationship to M1 oscillations can be used as the basis for a novel classification of SNr into two subpopulations. These results give insight into how dopamine loss leads to motor dysfunction and suggest a reappraisal of delta oscillations as a marker of akinetic symptoms in PD.NEW & NOTEWORTHY This work introduces a novel method to detect spike oscillations amidst neural noise. Using this method, we demonstrate that delta oscillations in the basal ganglia are a defining feature of awake, dopamine-depleted mice and are strongly correlated with dopamine loss and parkinsonian motor symptoms. These oscillations arise from a loss of D2-receptor activation and do not require motor cortex. Similar oscillations in human patients may be an underappreciated marker and target for Parkinson's disease (PD) treatment.

Self-Tuning Deep Brain Stimulation Controller for Suppression of Beta Oscillations: Analytical Derivation and Numerical Validation

Closed-loop control strategies for deep brain stimulation (DBS) in Parkinson's disease offer the potential to provide more effective control of patient symptoms and fewer side effects than continuous stimulation, while reducing battery consumption. Most of the closed-loop methods proposed and tested to-date rely on controller parameters, such as controller gains, that remain constant over time. While the controller may operate effectively close to the operating point for which it is set, providing benefits when compared to conventional open-loop DBS, it may perform sub-optimally if the operating conditions evolve. Such changes may result from, for example, diurnal variation in symptoms, disease progression or changes in the properties of the electrode-tissue interface. In contrast, an adaptive or “self-tuning” control mechanism has the potential to accommodate slowly varying changes in system properties over a period of days, months, or years. Such an adaptive mechanism would automatically adjust the controller parameters to maintain the desired performance while limiting side effects, despite changes in the system operating point. In this paper, two neural modeling approaches are utilized to derive and test an adaptive control scheme for closed-loop DBS, whereby the gain of a feedback controller is continuously adjusted to sustain suppression of pathological beta-band oscillatory activity at a desired target level. First, the controller is derived based on a simplified firing-rate model of the reciprocally connected subthalamic nucleus (STN) and globus pallidus (GPe). Its efficacy is shown both when pathological oscillations are generated endogenously within the STN-GPe network and when they arise in response to exogenous cortical STN inputs. To account for more realistic biological features, the control scheme is then tested in a physiologically detailed model of the cortical basal ganglia network, comprised of individual conductance-based spiking neurons, and simulates the coupled DBS electric field and STN local field potential. Compared to proportional feedback methods without gain adaptation, the proposed adaptive controller was able to suppress beta-band oscillations with less power consumption, even as the properties of the controlled system evolve over time due to alterations in the target for beta suppression, beta fluctuations and variations in the electrode impedance.

Essential chaotic dynamics of chatter in turning processes

Large-amplitude oscillations of turning operations are investigated, which can be modelled by a one degree-of-freedom damped oscillator subjected to the regenerative effect that introduces a relevant time delay in the system. In the case of large oscillations, when the cutting tool loses contact with the surface of the workpiece, the time delay is switched off, leading to a non-smooth delay differential equation. To explore the geometric structure of the global dynamics of the system, the mathematical model is approximated by means of the essential part of the spectrum in the region where the stationary cutting process may lose its stability. The trajectories embedded in the infinite-dimensional phase space are interpreted in a three-dimensional subspace and then analyzed by means of a discrete Lorenz-map. The bifurcation diagrams of the non-smooth system include chaotic windows, which are presented by numerical and semi-analytical tools and compared to the existing results in the literature.

Beta-Band Resonance and Intrinsic Oscillations in a Biophysically Detailed Model of the Subthalamic Nucleus-Globus Pallidus Network

Increased beta-band oscillatory activity in the basal ganglia network is associated with Parkinsonian motor symptoms and is suppressed with medication and deep brain stimulation (DBS). The origins of the beta-band oscillations, however, remains unclear with both intrinsic oscillations arising within the subthalamic nucleus (STN)-external globus pallidus (GPe) network and exogenous beta-activity, originating outside the network, proposed as potential sources of the pathological activity. The aim of this study was to explore the relative contribution of autonomous oscillations and exogenous oscillatory inputs in the generation of pathological oscillatory activity in a biophysically detailed model of the parkinsonian STN-GPe network. The network model accounts for the integration of synaptic currents and their interaction with intrinsic membrane currents in dendritic structures within the STN and GPe. The model was used to investigate the development of beta-band synchrony and bursting within the STN-GPe network by changing the balance of excitation and inhibition in both nuclei, and by adding exogenous oscillatory inputs with varying phase relationships through the hyperdirect cortico-subthalamic and indirect striato-pallidal pathways. The model showed an intrinsic susceptibility to beta-band oscillations that was manifest in weak autonomously generated oscillations within the STN-GPe network and in selective amplification of exogenous beta-band synaptic inputs near the network's endogenous oscillation frequency. The frequency at which this resonance peak occurred was determined by the net level of excitatory drive to the network. Intrinsic or endogenously generated oscillations were too weak to support a pacemaker role for the STN-GPe network, however, they were considerably amplified by sparse cortical beta inputs and were further amplified by striatal beta inputs that promoted anti-phase firing of the cortex and GPe, resulting in maximum transient inhibition of STN neurons. The model elucidates a mechanism of cortical patterning of the STN-GPe network through feedback inhibition whereby intrinsic susceptibility to beta-band oscillations can lead to phase locked spiking under parkinsonian conditions. These results point to resonance of endogenous oscillations with exogenous patterning of the STN-GPe network as a mechanism of pathological synchronization, and a role for the pallido-striatal feedback loop in amplifying beta oscillations.

Regulation and control roles of the basal ganglia in the development of absence epileptiform activities

Absence epileptiform activities are traditionally considered to be primarily induced by abnormal interactions between the cortical and thalamic neurons, which form the thalamocortical circuit in the brain. The basal ganglia, as an organizational unit in the brain, has close input and output relationships with the thalamocortical circuit. Although several studies report that the basal ganglia may participate in controlling and regulating absence epileptiform activities, to date, there have been no studies regarding whether the basal ganglia directly cause absence epileptiform activities. In this paper, we built a basal ganglia-corticothalamic network model to determine the role of basal ganglia in this disease. We determined that absence epileptiform activities might be directly induced by abnormal coupling strengths on certain pivotal pathways in the basal ganglia. These epileptiform activities can be well controlled by the coupling strengths of three major pathways that project from the thalamocortical network to the basal ganglia. The results implied that the substantia nigra pars compacta (SNc) can be considered to be the effective treatment target area for inhibiting epileptiform activities, which supports the observations of previous studies. Particularly, as a major contribution of this paper, we determined that the final presentation position of the epileptic slow spike waves is not limited to the cerebral cortex; these waves may additionally appear in the thalamus, striatal medium spiny neurons, striatal fast spiking interneuron, the SNc, subthalamic nucleus, substantia nigra pars reticulata and globus pallidus pars externa. In addition, consistent with several previous studies, the delay in the network was observed to be a critical factor for inducing transitions between different types of absence epileptiform activities. Our new model not only explains the onset and control mechanism but also provides a unified framework to study similar problems in neuron systems.

Different Dopaminergic Dysfunctions Underlying Parkinsonian Akinesia and Tremor

Although the occurrence of Parkinsonian akinesia and tremor is traditionally associated to dopaminergic degeneration, the multifaceted neural processes that cause these impairments are not fully understood. As a consequence, current dopamine medications cannot be tailored to the specific dysfunctions of patients with the result that generic drug therapies produce different effects on akinesia and tremor. This article proposes a computational model focusing on the role of dopamine impairments in the occurrence of akinesia and resting tremor. The model has three key features, to date never integrated in a single computational system: (a) an architecture constrained on the basis of the relevant known system-level anatomy of the basal ganglia-thalamo-cortical loops; (b) spiking neurons with physiologically-constrained parameters; (c) a detailed simulation of the effects of both phasic and tonic dopamine release. The model exhibits a neural dynamics compatible with that recorded in the brain of primates and humans. Moreover, it suggests that akinesia might involve both tonic and phasic dopamine dysregulations whereas resting tremor might be primarily caused by impairments involving tonic dopamine release and the responsiveness of dopamine receptors. These results could lead to develop new therapies based on a system-level view of the Parkinson's disease and targeting phasic and tonic dopamine in differential ways.

The beta oscillation conditions in a simplified basal ganglia network

Parkinson's disease is a type of motor dysfunction disease that is induced mainly by abnormal interactions between the subthalamic nucleus (STN) and globus pallidus (GP) neurons. Periodic oscillatory activities with frequencies of 13-30 Hz are the main physiological characteristics of Parkinson's disease. In this paper, we built a class of STN-GP networks to explore beta oscillation conditions. A theoretical formula was obtained for generating oscillations without internal GP connections. Based on this formula, we studied the effects of cortex inputs, striatum inputs, coupling weights and delays on oscillation conditions, and the theoretical results are in good agreement with the numerical results. The onset mechanism can be explained by the model, and the internal GP connection has little effect on oscillations. Finally, we compared oscillation conditions with those in previous studies and found that the delays and coupling weights required for generating oscillations may decrease as the number of nuclei increases. We hope that the results obtained will inspire future theoretical and experimental studies.

Insights into Parkinson's disease from computational models of the basal ganglia

Movement disorders arise from the complex interplay of multiple changes to neural circuits. Successful treatments for these disorders could interact with these complex changes in myriad ways, and as a consequence their mechanisms of action and their amelioration of symptoms are incompletely understood. Using Parkinson's disease as a case study, we review here how computational models are a crucial tool for taming this complexity, across causative mechanisms, consequent neural dynamics and treatments. For mechanisms, we review models that capture the effects of losing dopamine on basal ganglia function; for dynamics, we discuss models that have transformed our understanding of how beta-band (15-30 Hz) oscillations arise in the parkinsonian basal ganglia. For treatments, we touch on the breadth of computational modelling work trying to understand the therapeutic actions of deep brain stimulation. Collectively, models from across all levels of description are providing a compelling account of the causes, symptoms and treatments for Parkinson's disease.

Optogenetic Stimulation of GABAergic Neurons in the Globus Pallidus Produces Hyperkinesia

The globus pallidus (GP) is emerging as a critical locus of basal ganglia control of motor activity, but the exact role of GABAergic GP neurons remain to be defined. By targeted expression of channelrhodopsin 2 (ChR2) in GABAergic neurons using the VGAT-ChR2-EYFP transgenic mice, we showed that optogenetic stimulation of GABAergic neurons in the right GP produced hyperkinesia. Optogenetic stimulation of GABAergic GP neurons increased c-Fos-positive cells in GP, M1 cortex, and caudate-putamen (CPu), and decreased c-Fos-positive cells in entopeduncular nucleus (EPN), compared to the contralateral hemisphere. In agreement with the canonical basal ganglia model. Furthermore, we delivered AAV-CaMKIIα-ChR2-mCherry virus to the excitatory neurons of the subthalamic nucleus (STN) and selectively stimulated glutamatergic afferent fibers from the STN onto the GP. This optogenetic stimulation produced abnormal movements, similar to the behaviors that observed in the VGAT-ChR2-EYFP transgenic mice. Meanwhile, we found that the c-Fos expression pattern in the GP, M1, STN, EPN, and CPu produced by optogenetic activation of glutamatergic afferent fibers from the STN in GP was similar to the c-Fos expression pattern in the VGAT-ChR2-EYFP transgenic mice. Taken together, our results suggest that excess GP GABAergic neurons activity could be the neural substrate of abnormal involuntary movements in hyperkinetic movement disorders. The neural circuitry underlying the abnormal involuntary movements is associated with excessive GP, M1, CPu activity, and reduced EPN activity. Inhibition of GP GABAergic neurons represents new treatment targets for hyperkinetic movement disorder.

Generic dynamic causal modelling: An illustrative application to Parkinson's disease

We present a technical development in the dynamic causal modelling of electrophysiological responses that combines qualitatively different neural mass models within a single network. This affords the option to couple various cortical and subcortical nodes that differ in their form and dynamics. Moreover, it enables users to implement new neural mass models in a straightforward and standardized way. This generic framework hence supports flexibility and facilitates the exploration of increasingly plausible models. We illustrate this by coupling a basal ganglia-thalamus model to a (previously validated) cortical model developed specifically for motor cortex. The ensuing DCM is used to infer pathways that contribute to the suppression of beta oscillations induced by dopaminergic medication in patients with Parkinson's disease. Experimental recordings were obtained from deep brain stimulation electrodes (implanted in the subthalamic nucleus) and simultaneous magnetoencephalography. In line with previous studies, our results indicate a reduction of synaptic efficacy within the circuit between the subthalamic nucleus and external pallidum, as well as reduced efficacy in connections of the hyperdirect and indirect pathway leading to this circuit. This work forms the foundation for a range of modelling studies of the synaptic mechanisms (and pathophysiology) underlying event-related potentials and cross-spectral densities.

Model and Analysis for the Onset of Parkinsonian Firing Patterns in a Simplified Basal Ganglia

Parkinson's disease (PD) is a degenerative neurological disease that disrupts the movement cycle in the basal ganglia. As the disease progresses, dopamine depletion leads to changes to how the basal ganglia functions as well as the appearance of abnormal beta oscillations. There is much debate on just exactly how these connection strengths change and just how the oscillations emerge. One leading hypothesis claims that the oscillations develop in the globus pallidus external, subthalamic nucleus, and globus pallidus internal loop. We introduce a mathematical model that calculates the average firing rates of this loop while still accounting for the larger closed loop of the entire basal ganglia system. This model is constructed such that physiologically realistic results can be obtained while not sacrificing the use of analytic methods. Because of this, it is possible to determine how the change in the connection strengths can drive the necessary changes in firing rates seen in recordings and account for the generation of trademark beta oscillations of PD without relying on highly specific time delays, stochastic approaches, or numerical approximations. Additionally, we find that the entire cortico-basal ganglia-thalamo-cortical loop is essential for abnormal oscillations to originate.

Propagation of beta/gamma rhythms in the cortico-basal ganglia circuits of the parkinsonian rat

Much of the motor impairment associated with Parkinson’s disease is thought to arise from pathological activity in the networks formed by the basal ganglia (BG) and motor cortex. To evaluate several hypotheses proposed to explain the emergence of pathological oscillations in parkinsonism, we investigated changes to the directed connectivity in BG networks following dopamine depletion. We recorded local field potentials (LFPs) in the cortex and basal ganglia of rats rendered parkinsonian by injection of 6-hydroxydopamine (6-OHDA) and in dopamine-intact controls. We performed systematic analyses of the networks using a novel tool for estimation of directed interactions (nonparametric directionality, NPD). We used a “conditioned” version of the NPD analysis that reveals the dependence of the correlation between two signals on a third reference signal. We find evidence of the dopamine dependency of both low-beta (14–20 Hz) and high-beta/low-gamma (20–40 Hz) directed network interactions. Notably, 6-OHDA lesions were associated with enhancement of the cortical “hyperdirect” connection to the subthalamic nucleus (STN) and its feedback to the cortex and striatum. We find that pathological beta synchronization resulting from 6-OHDA lesioning is widely distributed across the network and cannot be located to any individual structure. Furthermore, we provide evidence that high-beta/gamma oscillations propagate through the striatum in a pathway that is independent of STN. Rhythms at high beta/gamma show susceptibility to conditioning that indicates a hierarchical organization compared with those at low beta. These results further inform our understanding of the substrates for pathological rhythms in salient brain networks in parkinsonism.

NEW & NOTEWORTHY We present a novel analysis of electrophysiological recordings in the cortico-basal ganglia network with the aim of evaluating several hypotheses concerning the origins of abnormal brain rhythms associated with Parkinson’s disease. We present evidence for changes in the directed connections within the network following chronic dopamine depletion in rodents. These findings speak to the plausibility of a “short-circuiting” of the network that gives rise to the conditions from which pathological synchronization may arise.

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.

A Population of Indirect Pathway Striatal Projection Neurons Is Selectively Entrained to Parkinsonian Beta Oscillations

Classical schemes of basal ganglia organization posit that parkinsonian movement difficulties presenting after striatal dopamine depletion stem from the disproportionate firing rates of spiny projection neurons (SPNs) therein. There remains, however, a pressing need to elucidate striatal SPN firing in the context of the synchronized network oscillations that are abnormally exaggerated in cortical-basal ganglia circuits in parkinsonism. To address this, we recorded unit activities in the dorsal striatum of dopamine-intact and dopamine-depleted rats during two brain states, respectively defined by cortical slow-wave activity (SWA) and activation. Dopamine depletion escalated striatal net output but had contrasting effects on "direct pathway" SPNs (dSPNs) and "indirect pathway" SPNs (iSPNs); their firing rates became imbalanced, and they disparately engaged in network oscillations. Disturbed striatal activity dynamics relating to the slow (∼1 Hz) oscillations prevalent during SWA partly generalized to the exaggerated beta-frequency (15-30 Hz) oscillations arising during cortical activation. In both cases, SPNs exhibited higher incidences of phase-locked firing to ongoing cortical oscillations, and SPN ensembles showed higher levels of rhythmic correlated firing, after dopamine depletion. Importantly, in dopamine-depleted striatum, a widespread population of iSPNs, which often displayed excessive firing rates and aberrant phase-locked firing to cortical beta oscillations, preferentially and excessively synchronized their firing at beta frequencies. Conversely, dSPNs were neither hyperactive nor synchronized to a large extent during cortical activation. These data collectively demonstrate a cell type-selective entrainment of SPN firing to parkinsonian beta oscillations. We conclude that a population of overactive, excessively synchronized iSPNs could orchestrate these pathological rhythms in basal ganglia circuits.SIGNIFICANCE STATEMENT Chronic depletion of dopamine from the striatum, a part of the basal ganglia, causes some symptoms of Parkinson's disease. Here, we elucidate how dopamine depletion alters striatal neuron firing in vivo, with an emphasis on defining whether and how spiny projection neurons (SPNs) engage in the synchronized beta-frequency (15-30 Hz) oscillations that become pathologically exaggerated throughout basal ganglia circuits in parkinsonism. We discovered that a select population of so-called "indirect pathway" SPNs not only fire at abnormally high rates, but are also particularly prone to being recruited to exaggerated beta oscillations. Our results provide an important link between two complementary theories that explain the presentation of disease symptoms on the basis of changes in firing rate or firing synchronization/rhythmicity.

Computational models of basal ganglia dysfunction: the dynamics is in the details

The development, simulation, and analysis of mathematical models offer helpful tools for integrating experimental findings and exploring or suggesting possible explanatory mechanisms. As models relating to basal ganglia dysfunction have proliferated, however, there has not always been consistency among their findings. This work points out several ways in which biological details, relating to ionic currents and synaptic pathways, can influence the dynamics of models of the basal ganglia under parkinsonian conditions and hence may be important for inclusion in models. It also suggests some additional useful directions for future modeling studies relating to basal ganglia dysfunction.

Quantification of motor network dynamics in Parkinson's disease by means of landscape and flux theory

The basal ganglia neural circuit plays an important role in motor control. Despite the significant efforts, the understanding of the principles and underlying mechanisms of this modulatory circuit and the emergence of abnormal synchronized oscillations in movement disorders is still challenging. Dopamine loss has been proved to be responsible for Parkinson's disease. We quantitatively described the dynamics of the basal ganglia-thalamo-cortical circuit in Parkinson's disease in terms of the emergence of both abnormal firing rates and firing patterns in the circuit. We developed a potential landscape and flux framework for exploring the modulatory circuit. The driving force of the circuit can be decomposed into a gradient of the potential, which is associated with the steady-state probability distributions, and the curl probability flux term. We uncovered the underlying potential landscape as a Mexican hat-shape closed ring valley where abnormal oscillations emerge due to dopamine depletion. We quantified the global stability of the network through the topography of the landscape in terms of the barrier height, which is defined as the potential difference between the maximum potential inside the ring and the minimum potential along the ring. Both a higher barrier and a larger flux originated from detailed balance breaking result in more stable oscillations. Meanwhile, more energy is consumed to support the increasing flux. Global sensitivity analysis on the landscape topography and flux indicates how changes in underlying neural network regulatory wirings and external inputs influence the dynamics of the system. We validated two of the main hypotheses(direct inhibition hypothesis and output activation hypothesis) on the therapeutic mechanism of deep brain stimulation (DBS). We found GPe appears to be another effective stimulated target for DBS besides GPi and STN. Our approach provides a general way to quantitatively explore neural networks and may help for uncovering more efficacious therapies for movement disorders.

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.

Ventral Medial Thalamic Nucleus Promotes Synchronization of Increased High Beta Oscillatory Activity in the Basal Ganglia-Thalamocortical Network of the Hemiparkinsonian Rat

Loss of dopamine is associated with increased synchronization and oscillatory activity in the subthalamic nucleus and basal ganglia (BG) output nuclei in both Parkinson's disease (PD) patients and animal models of PD. We have previously observed substantial increases in spectral power in the 25-40 Hz range in LFPs recorded in the substantia nigra pars reticulata (SNpr) and motor cortex (MCx) in the hemiparkinsonian rat during treadmill walking. The current study explores the hypothesis that SNpr output entrains activity in the ventral medial thalamus (VM) in this frequency range after loss of dopamine, which in turn contributes to entrainment of the MCx and BG. Electrode bundles were implanted in MCx, SNpr, and VM of rats with unilateral dopamine cell lesions. Spiking and LFP activity were recorded during epochs of rest and walking on a circular treadmill. After dopamine cell lesion, 30-36 Hz LFP activity in the VM became more robust during treadmill walking and more coherent with LFP activity in the same range in MCx and SNpr. Infusion of the GABAA antagonist picrotoxin into the VM reduced both high beta power in MCx and SNpr and coherence between MCx and SNpr while temporarily restoring walking ability. Infusion of the GABAA agonist muscimol into the VM also reduced MCx-SNpr coherence and beta power but failed to improve walking. These results support the view that synchronized neuronal activity in the VM contributes to the emergence of high beta oscillations throughout the BG-thalamocortical network in the behaving parkinsonian rat.

SIGNIFICANCE STATEMENT

Parkinson's disease symptoms are associated with dramatic increases in synchronized beta range (15-35 Hz) oscillatory local field activity in several brain areas involved in motor control, but the mechanisms promoting this activity and its functional significance remain unresolved. This oscillatory activity can be recorded in awake behaving rats with unilateral dopamine cell lesions using chronically implanted electrodes. Although these rats have motor deficits, they can walk on a circular treadmill in the direction ipsilateral to their lesion. This study establishes a critical role for the ventral medial thalamus in the propagation of this exaggerated beta range oscillatory activity and the sequential entrainment of structures throughout the basal ganglia-thalamocortical loop in the lesioned hemisphere of hemiparkinsonian rats during treadmill walking.

Properties of Neurons in External Globus Pallidus Can Support Optimal Action Selection

The external globus pallidus (GPe) is a key nucleus within basal ganglia circuits that are thought to be involved in action selection. A class of computational models assumes that, during action selection, the basal ganglia compute for all actions available in a given context the probabilities that they should be selected. These models suggest that a network of GPe and subthalamic nucleus (STN) neurons computes the normalization term in Bayes’ equation. In order to perform such computation, the GPe needs to send feedback to the STN equal to a particular function of the activity of STN neurons. However, the complex form of this function makes it unlikely that individual GPe neurons, or even a single GPe cell type, could compute it. Here, we demonstrate how this function could be computed within a network containing two types of GABAergic GPe projection neuron, so-called ‘prototypic’ and ‘arkypallidal’ neurons, that have different response properties in vivo and distinct connections. We compare our model predictions with the experimentally-reported connectivity and input-output functions (f-I curves) of the two populations of GPe neurons. We show that, together, these dichotomous cell types fulfil the requirements necessary to compute the function needed for optimal action selection. We conclude that, by virtue of their distinct response properties and connectivities, a network of arkypallidal and prototypic GPe neurons comprises a neural substrate capable of supporting the computation of the posterior probabilities of actions.

Choosing an appropriate action as quickly and accurately as possible in a given situation is critical for the survival of animals and humans. One of the brain regions involved in action selection is a set of subcortical nuclei known as the basal ganglia. The importance of understanding information processing in the basal ganglia is further emphasised by the fact that their disturbed interactions in Parkinson’s disease results in profound difficulties in movement. Computational models have suggested how the basal ganglia could select actions in the fastest possible way for the required accuracy level. These models further predict that a part of basal ganglia, called the external globus pallidus (GPe), needs to calculate a particular function of its inputs. This paper proposes how this function could be computed in a mathematical model of a network within GPe. Furthermore, it shows that the experimentally observed connectivity and response properties of GPe neurons fulfil the requirements necessary to support optimal action selection. This suggests the GPe neurons have properties that allow them to contribute to optimal action selection in the whole basal ganglia.

The network of causal interactions for beta oscillations in the pedunculopontine nucleus, primary motor cortex, and subthalamic nucleus of walking parkinsonian rats

Oscillatory activity has been well-studied in many structures within cortico-basal ganglia circuits, but it is not well understood within the pedunculopontine nucleus (PPN), which was recently introduced as a potential target for the treatment of gait and postural impairments in advanced stages of Parkinson's disease (PD). To investigate oscillatory activity in the PPN and its relationship with oscillatory activity in cortico-basal ganglia circuits, we simultaneously recorded local field potentials in the PPN, primary motor cortex (M1), and subthalamic nucleus (STN) of 6-hydroxydopamine (6-OHDA)-induced hemiparkinsonian rats during resting and walking. After analysis of power spectral density, coherence, and partial Granger causality, three major findings emerged: 1) after 6-OHDA lesions, beta band oscillations were enhanced in all three regions during walking; 2) the direction of information flow for beta oscillations among the three structures was STN→M1, STN→PPN, and PPN→M1; 3) after the treatment of levodopa, beta activity in the three regions was reduced significantly and the flow of beta band was also abrogated. Our results suggest that beta activity in the PPN is transmitted from the basal ganglia and probably comes from the STN, and the STN plays a dominant role in the network of causal interactions for beta activity. Thus, the STN may be a potential source of aberrant beta band oscillations in PD. Levodopa can inhibit beta activity in the PPN of parkinsonian rats but cannot relieve parkinsonian patients' axial symptoms clinically. Therefore, beta oscillations may not be the major cause of axial symptoms.

Closed-loop firing rate regulation of two interacting excitatory and inhibitory neural populations of the basal ganglia

This paper develops a new closed-loop firing rate regulation strategy for a population of neurons in the subthalamic nucleus, derived using a model-based analysis of the basal ganglia. The system is described using a firing rate model, in order to analyse the generation of beta-band oscillations. On this system, a proportional regulation of the firing rate reduces the gain of the subthalamo-pallidal loop in the parkinsonian case, thus impeding pathological oscillation generation. A filter with a well-chosen frequency is added to this proportional scheme, in order to avoid a potential instability of the feedback loop due to actuation and measurement delays. Our main result is a set of conditions on the parameters of the stimulation strategy that guarantee both its stability and a prescribed delay margin. A discussion on the applicability of the proposed method and a complete set of mathematical proofs is included.

Computational Models Describing Possible Mechanisms for Generation of Excessive Beta Oscillations in Parkinson's Disease

In Parkinson's disease, an increase in beta oscillations within the basal ganglia nuclei has been shown to be associated with difficulty in movement initiation. An important role in the generation of these oscillations is thought to be played by the motor cortex and by a network composed of the subthalamic nucleus (STN) and the external segment of globus pallidus (GPe). Several alternative models have been proposed to describe the mechanisms for generation of the Parkinsonian beta oscillations. However, a recent experimental study of Tachibana and colleagues yielded results which are challenging for all published computational models of beta generation. That study investigated how the presence of beta oscillations in a primate model of Parkinson's disease is affected by blocking different connections of the STN-GPe circuit. Due to a large number of experimental conditions, the study provides strong constraints that any mechanistic model of beta generation should satisfy. In this paper we present two models consistent with the data of Tachibana et al. The first model assumes that Parkinsonian beta oscillation are generated in the cortex and the STN-GPe circuits resonates at this frequency. The second model additionally assumes that the feedback from STN-GPe circuit to cortex is important for maintaining the oscillations in the network. Predictions are made about experimental evidence that is required to differentiate between the two models, both of which are able to reproduce firing rates, oscillation frequency and effects of lesions carried out by Tachibana and colleagues. Furthermore, an analysis of the models reveals how the amplitude and frequency of the generated oscillations depend on parameters.

A Biologically Inspired Computational Model of Basal Ganglia in Action Selection

The basal ganglia (BG) are a subcortical structure implicated in action selection. The aim of this work is to present a new cognitive neuroscience model of the BG, which aspires to represent a parsimonious balance between simplicity and completeness. The model includes the 3 main pathways operating in the BG circuitry, that is, the direct (Go), indirect (NoGo), and hyperdirect pathways. The main original aspects, compared with previous models, are the use of a two-term Hebb rule to train synapses in the striatum, based exclusively on neuronal activity changes caused by dopamine peaks or dips, and the role of the cholinergic interneurons (affected by dopamine themselves) during learning. Some examples are displayed, concerning a few paradigmatic cases: action selection in basal conditions, action selection in the presence of a strong conflict (where the role of the hyperdirect pathway emerges), synapse changes induced by phasic dopamine, and learning new actions based on a previous history of rewards and punishments. Finally, some simulations show model working in conditions of altered dopamine levels, to illustrate pathological cases (dopamine depletion in parkinsonian subjects or dopamine hypermedication). Due to its parsimonious approach, the model may represent a straightforward tool to analyze BG functionality in behavioral experiments.

Role of the indirect pathway of the basal ganglia in perceptual decision making

The basal ganglia (BG) play an important role in motor control, reinforcement learning, and perceptual decision making. Modeling and experimental evidence suggest that, in a speed-accuracy tradeoff, the corticostriatal pathway can adaptively adjust a decision threshold (the amount of information needed to make a choice). In this study, we go beyond the focus of previous works on the direct and hyperdirect pathways to examine the contribution of the indirect pathway of the BG system to decision making in a biophysically based spiking network model. We find that the mechanism of adjusting the decision threshold by plasticity of the corticostriatal connections is effective, provided that the indirect pathway counterbalances the direct pathway in their projections to the output nucleus. Furthermore, in our model, changes within basal ganglia connections similar to those that arise in parkinsonism give rise to strong beta oscillations. Specifically, beta oscillations are produced by an abnormal enhancement of the interactions between the subthalamic nucleus (STN) and the external segment of globus pallidus (GPe) in the indirect pathway, with an oscillation frequency that depends on the excitatory cortical input to the STN and the inhibitory input to the GPe from the striatum. In a parkinsonian state characterized by pronounced beta oscillations, the mean reaction time and range of threshold variation (a measure of behavioral flexibility) are significantly reduced compared with the normal state. Our work thus reveals a specific circuit mechanism for impairments of perceptual decision making associated with Parkinson's disease.

A canonical circuit for generating phase-amplitude coupling

‘Phase amplitude coupling’ (PAC) in oscillatory neural activity describes a phenomenon whereby the amplitude of higher frequency activity is modulated by the phase of lower frequency activity. Such coupled oscillatory activity – also referred to as ‘cross-frequency coupling’ or ‘nested rhythms’ – has been shown to occur in a number of brain regions and at behaviorally relevant time points during cognitive tasks; this suggests functional relevance, but the circuit mechanisms of PAC generation remain unclear. In this paper we present a model of a canonical circuit for generating PAC activity, showing how interconnected excitatory and inhibitory neural populations can be periodically shifted in to and out of oscillatory firing patterns by afferent drive, hence generating higher frequency oscillations phase-locked to a lower frequency, oscillating input signal. Since many brain regions contain mutually connected excitatory-inhibitory populations receiving oscillatory input, the simplicity of the mechanism generating PAC in such networks may explain the ubiquity of PAC across diverse neural systems and behaviors. Analytic treatment of this circuit as a nonlinear dynamical system demonstrates how connection strengths and inputs to the populations can be varied in order to change the extent and nature of PAC activity, importantly which phase of the lower frequency rhythm the higher frequency activity is locked to. Consequently, this model can inform attempts to associate distinct types of PAC with different network topologies and physiologies in real data.

Functional correlates of exaggerated oscillatory activity in basal ganglia output in hemiparkinsonian rats

Exaggerated beta range (13-30Hz) synchronized activity is observed in the basal ganglia of Parkinson's disease (PD) patients during implantation of deep brain stimulation electrodes and is thought to contribute to the motor symptoms of this disorder. To explore the translational potential of similar activity observed in a rat model of PD, local field potentials (LFPs) and spiking activity in basal ganglia output were characterized in rats with unilateral dopamine cell lesion during a range of behaviors. A circular treadmill was used to assess activity during walking; hemiparkinsonian rats could maintain a steady gait when oriented ipsiversive to the lesioned hemisphere, but were less effective at walking when oriented contraversive to lesion. Dramatic increases in substantia nigra pars reticulata (SNpr) LFP oscillatory activity and spike-LFP synchronization were observed within the beta/low gamma range (12-40Hz) in the lesioned hemisphere, relative to the non-lesioned hemisphere, with the dominant frequency of spike-LFP entrainment and LFP power varying with behavioral state. At 3weeks postlesion, the mean dominant entrainment frequency during ipsiversive treadmill walking and grooming was 34Hz. Other behaviors were associated with lower mean entrainment frequencies: 27-28Hz during alert non-walking and REM, 17Hz during rest and 21Hz during urethane anesthesia with sensory stimulation. SNpr spike-LFP entrainment frequency was stable during individual treadmill walking epochs, but increased gradually over weeks postlesion. In contrast, SNpr LFP power in the 25-40Hz range was greatest at the initiation of each walking epoch, and decreased during walking to stabilize by 6min at 49% of initial values. Power was further modulated in conjunction with the 1.5s stepping rhythm. Administration of l-dopa improved contraversive treadmill walking in correlation with a reduction in SNpr 25-40Hz LFP power and spike synchronization in the dopamine cell lesioned hemisphere. These effects were reversed by the serotonergic 1A agonist, 8-OH-DPAT. While the prominent spike-LFP phase locking observed during ongoing motor activity in the hemiparkinsonian rats occurs at frequencies intriguingly higher than in PD patients, the synchronized activity in the SNpr of this animal model has much in common with oscillatory activity recorded from the basal ganglia of the PD patients. Results support the potential of this model for providing insight into relationships between synchronization of basal ganglia output induced by loss of dopamine and motor symptoms in PD.

Closing the loop of deep brain stimulation

High-frequency deep brain stimulation is used to treat a wide range of brain disorders, like Parkinson's disease. The stimulated networks usually share common electrophysiological signatures, including hyperactivity and/or dysrhythmia. From a clinical perspective, HFS is expected to alleviate clinical signs without generating adverse effects. Here, we consider whether the classical open-loop HFS fulfills these criteria and outline current experimental or theoretical research on the different types of closed-loop DBS that could provide better clinical outcomes. In the first part of the review, the two routes followed by HFS-evoked axonal spikes are explored. In one direction, orthodromic spikes functionally de-afferent the stimulated nucleus from its downstream target networks. In the opposite direction, antidromic spikes prevent this nucleus from being influenced by its afferent networks. As a result, the pathological synchronized activity no longer propagates from the cortical networks to the stimulated nucleus. The overall result can be described as a reversible functional de-afferentation of the stimulated nucleus from its upstream and downstream nuclei. In the second part of the review, the latest advances in closed-loop DBS are considered. Some of the proposed approaches are based on mathematical models, which emphasize different aspects of the parkinsonian basal ganglia: excessive synchronization, abnormal firing-rate rhythms, and a deficient thalamo-cortical relay. The stimulation strategies are classified depending on the control-theory techniques on which they are based: adaptive and on-demand stimulation schemes, delayed and multi-site approaches, stimulations based on proportional and/or derivative control actions, optimal control strategies. Some of these strategies have been validated experimentally, but there is still a large reservoir of theoretical work that may point to ways of improving practical treatment.

Effective connectivity of the subthalamic nucleus-globus pallidus network during Parkinsonian oscillations

In Parkinsonism, subthalamic nucleus (STN) neurons and two types of external globus pallidus (GP) neuron inappropriately synchronise their firing in time with slow (∼1 Hz) or beta (13-30 Hz) oscillations in cortex. We recorded the activities of STN, Type-I GP (GP-TI) and Type-A GP (GP-TA) neurons in anaesthetised Parkinsonian rats during such oscillations to constrain a series of computational models that systematically explored the effective connections and physiological parameters underlying neuronal rhythmic firing and phase preferences in vivo. The best candidate model, identified with a genetic algorithm optimising accuracy/complexity measures, faithfully reproduced experimental data and predicted that the effective connections of GP-TI and GP-TA neurons are quantitatively different. Estimated inhibitory connections from striatum were much stronger to GP-TI neurons than to GP-TA neurons, whereas excitatory connections from thalamus were much stronger to GP-TA and STN neurons than to GP-TI neurons. Reciprocal connections between GP-TI and STN neurons were matched in weight, but those between GP-TA and STN neurons were not; only GP-TI neurons sent substantial connections back to STN. Different connection weights between and within the two types of GP neuron were also evident. Adding to connection differences, GP-TA and GP-TI neurons were predicted to have disparate intrinsic physiological properties, reflected in distinct autonomous firing rates. Our results elucidate potential substrates of GP functional dichotomy, and emphasise that rhythmic inputs from striatum, thalamus and cortex are important for setting activity in the STN-GP network during Parkinsonian beta oscillations, suggesting they arise from interactions between most nodes of basal ganglia-thalamocortical circuits.

Delay-induced oscillations in Wilson and Cowan's model: an analysis of the subthalamo-pallidal feedback loop in healthy and parkinsonian subjects

The model proposed by Wilson and Cowan (1972) describes the dynamics of two interacting subpopulations of excitatory and inhibitory neurons. It has been used to model neural structures like the olfactory bulb, whisker barrels, and the subthalamo-pallidal system. It is well-known that this system can exhibit an oscillatory behavior that is amplified by the presence of delays. In the absence of delays, the conditions for stability are well-known. The aim of our paper is to clarify these conditions when delays are included in the model. The first ingredient of our methods is a new necessary and sufficient condition for the existence of multiple equilibria. This condition is related to those for local asymptotic stability. In addition, a sufficient condition for global stability is also proposed. The second and main ingredient is a stability analysis of the system in the frequency-domain, based on the Nyquist criterion, that takes the four independent delays into account. The methods proposed in this paper can be applied to analyse the stability of the subthalamo-pallidal feedback loop, a deep brain structure involved in Parkinson's disease. Our stability conditions are easy to compute and characterize sharply the system's parameters for which spontaneous oscillations appear.

The dynamics of neurological disease: integrating computational, experimental and clinical neuroscience

There is a vast (and rapidly growing) amount of experimental and clinical data of the nervous system at very diverse spatial scales of activity (e.g. from sub-cellular through to whole organ), with many neurological disorders characterized by oscillations in neural activity across these disparate scales. Computer modelling and the development of associated mathematical theories provide us with a unique opportunity to integrate information from across these diverse scales of activity; leading to explanations of the potential mechanisms underlying the time-evolving dynamics and, more importantly, allowing the development of new hypotheses regarding neural function that may be tested experimentally and ultimately translated into the clinic. The purpose of this special issue is to present an overview of current integrative research in the areas of epilepsy, Parkinson's disease and schizophrenia, where multidisciplinary relationships involving theory, experimental and clinical research are becoming increasingly established.