Neural dynamics in parkinsonian brain: the boundary between synchronized and nonsynchronized dynamics

Intermittent generalized synchronization and modified system approach: Discrete maps

The present work deals with the intermittent generalized synchronization regime observed near the boundary of generalized synchronization. The intermittent behavior is considered in the context of two observable phenomena, namely, (i) the birth of the asynchronous stages of motion from the complete synchronous state and (ii) the multistability in detection of the synchronous and asynchronous states. The mechanisms governing these phenomena are revealed and described in this paper with the help of the modified system approach for unidirectionally coupled model oscillators with discrete time.

Synaptic Changes in Pallidostriatal Circuits Observed in the Parkinsonian Model Triggers Abnormal Beta Synchrony with Accurate Spatio-temporal Properties across the Basal Ganglia

Excessive oscillatory activity across basal ganglia (BG) nuclei in the β frequencies (12-30 Hz) is a hallmark of Parkinson's disease (PD). While the link between oscillations and symptoms remains debated, exaggerated β oscillations constitute an important biomarker for therapeutic effectiveness in PD. The neuronal mechanisms of β-oscillation generation however remain unknown. Many existing models rely on a central role of the subthalamic nucleus (STN) or cortical inputs to BG. Contrarily, neural recordings and optogenetic manipulations in normal and parkinsonian rats recently highlighted the central role of the external pallidum (GPe) in abnormal β oscillations, while showing that the integrity of STN or motor cortex is not required. Here, we evaluate the mechanisms for the generation of abnormal β oscillations in a BG network model where neuronal and synaptic time constants, connectivity, and firing rate distributions are strongly constrained by experimental data. Guided by a mean-field approach, we show in a spiking neural network that several BG sub-circuits can drive oscillations. Strong recurrent STN-GPe connections or collateral intra-GPe connections drive γ oscillations (>40 Hz), whereas strong pallidostriatal loops drive low-β (10-15 Hz) oscillations. We show that pathophysiological strengthening of striatal and pallidal synapses following dopamine depletion leads to the emergence of synchronized oscillatory activity in the mid-β range with spike-phase relationships between BG neuronal populations in-line with experiments. Furthermore, inhibition of GPe, contrary to STN, abolishes oscillations. Our modeling study uncovers the neural mechanisms underlying PD β oscillations and may thereby guide the future development of therapeutic strategies.

Biophysical and Architectural Mechanisms of Subthalamic Theta under Response Conflict

The cortico-basal ganglia circuit is needed to suppress prepotent actions and to facilitate controlled behavior. Under conditions of response conflict, the frontal cortex and subthalamic nucleus (STN) exhibit increased spiking and theta band power, which are linked to adaptive regulation of behavioral output. The electrophysiological mechanisms underlying these neural signatures of impulse control remain poorly understood. To address this lacuna, we constructed a novel large-scale, biophysically principled model of the subthalamopallidal (STN-globus pallidus externus) network and examined the mechanisms that modulate theta power and spiking in response to cortical input. Simulations confirmed that theta power does not emerge from intrinsic network dynamics but is robustly elicited in response to cortical input as burst events representing action selection dynamics. Rhythmic burst events of multiple cortical populations, representing a state of conflict where cortical motor plans vacillate in the theta range, led to prolonged STN theta and increased spiking, consistent with empirical literature. Notably, theta band signaling required NMDA, but not AMPA, currents, which were in turn related to a triphasic STN response characterized by spiking, silence, and bursting periods. Finally, theta band resonance was also strongly modulated by architectural connectivity, with maximal theta arising when multiple cortical populations project to individual STN "conflict detector" units because of an NMDA-dependent supralinear response. Our results provide insights into the biophysical principles and architectural constraints that give rise to STN dynamics during response conflict, and how their disruption can lead to impulsivity and compulsivity.SIGNIFICANCE STATEMENT The subthalamic nucleus exhibits theta band power modulation related to cognitive control over motor actions during conditions of response conflict. However, the mechanisms of such dynamics are not understood. Here we developed a novel biophysically detailed and data-constrained large-scale model of the subthalamopallidal network, and examined the impacts of cellular and network architectural properties that give rise to theta dynamics. Our investigations implicate an important role for NMDA receptors and cortico-subthalamic nucleus topographical connectivities in theta power modulation.

Subspace-based predictive control of Parkinson's disease: A model-based study

Deep brain stimulation (DBS) of the Basal Ganglia (BG) is an effective treatment to suppress the symptoms of Parkinson's disease (PD). Using a closed-loop scheme in DBS can not only improve its therapeutic effects but it can also reduce its energy consumption and possible side effects. In this paper, a predictive closed loop control strategy is employed to suppress the PD in real-time. A linear multi-input multi-output (MIMO) state-delayed system is considered as a simplified model of the BG neuronal network relating the stimulation signals as inputs to the beta power of local field potentials as PD biomarkers. The effect of time delay in different areas of the BG is incorporated into this model and a real-time subspace-based identification is implemented to continuously model the state of the BG neuronal network and drive the predictive control strategy. Simulation results show that the proposed MIMO subspace based predictive controller can suppress PD symptoms more effectively and with less power consumption compared to the conventional open-loop DBS and a recently proposed single-input single-output closed loop controller.

Flexible brain dynamics underpins complex behaviours as observed in Parkinson's disease

Rapid reconfigurations of brain activity support efficient neuronal communication and flexible behaviour. Suboptimal brain dynamics is associated to impaired adaptability, possibly leading to functional deficiencies. We hypothesize that impaired flexibility in brain activity can lead to motor and cognitive symptoms of Parkinson’s disease (PD). To test this hypothesis, we studied the ‘functional repertoire’—the number of distinct configurations of neural activity—using source-reconstructed magnetoencephalography in PD patients and controls. We found stereotyped brain dynamics and reduced flexibility in PD. The intensity of this reduction was proportional to symptoms severity, which can be explained by beta-band hyper-synchronization. Moreover, the basal ganglia were prominently involved in the abnormal patterns of brain activity. Our findings support the hypotheses that: symptoms in PD relate to impaired brain flexibility, this impairment preferentially involves the basal ganglia, and beta-band hypersynchronization is associated with reduced brain flexibility. These findings highlight the importance of extensive functional repertoires for correct behaviour.

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.

Uncoupling the roles of firing rates and spike bursts in shaping the STN-GPe beta band oscillations

The excess of 15-30 Hz (β-band) oscillations in the basal ganglia is one of the key signatures of Parkinson's disease (PD). The STN-GPe network is integral to generation and modulation of β band oscillations in basal ganglia. However, the role of changes in the firing rates and spike bursting of STN and GPe neurons in shaping these oscillations has remained unclear. In order to uncouple their effects, we studied the dynamics of STN-GPe network using numerical simulations. In particular, we used a neuron model, in which firing rates and spike bursting can be independently controlled. Using this model, we found that while STN firing rate is predictive of oscillations, GPe firing rate is not. The effect of spike bursting in STN and GPe neurons was state-dependent. That is, only when the network was operating in a state close to the border of oscillatory and non-oscillatory regimes, spike bursting had a qualitative effect on the β band oscillations. In these network states, an increase in GPe bursting enhanced the oscillations whereas an equivalent proportion of spike bursting in STN suppressed the oscillations. These results provide new insights into the mechanisms underlying the transient β bursts and how duration and power of β band oscillations may be controlled by an interplay of GPe and STN firing rates and spike bursts.

Basal ganglia role in learning rewarded actions and executing previously learned choices: Healthy and diseased states

The basal ganglia (BG) is a collection of nuclei located deep beneath the cerebral cortex that is involved in learning and selection of rewarded actions. Here, we analyzed BG mechanisms that enable these functions. We implemented a rate model of a BG-thalamo-cortical loop and simulated its performance in a standard action selection task. We have shown that potentiation of corticostriatal synapses enables learning of a rewarded option. However, these synapses became redundant later as direct connections between prefrontal and premotor cortices (PFC-PMC) were potentiated by Hebbian learning. After we switched the reward to the previously unrewarded option (reversal), the BG was again responsible for switching to the new option. Due to the potentiated direct cortical connections, the system was biased to the previously rewarded choice, and establishing the new choice required a greater number of trials. Guided by physiological research, we then modified our model to reproduce pathological states of mild Parkinson's and Huntington's diseases. We found that in the Parkinsonian state PMC activity levels become extremely variable, which is caused by oscillations arising in the BG-thalamo-cortical loop. The model reproduced severe impairment of learning and predicted that this is caused by these oscillations as well as a reduced reward prediction signal. In the Huntington state, the potentiation of the PFC-PMC connections produced better learning, but altered BG output disrupted expression of the rewarded choices. This resulted in random switching between rewarded and unrewarded choices resembling an exploratory phase that never ended. Along with other computational studies, our results further reconcile the apparent contradiction between the critical involvement of the BG in execution of previously learned actions and yet no impairment of these actions after BG output is ablated by lesions or deep brain stimulation. We predict that the cortico-BG-thalamo-cortical loop conforms to previously learned choice in healthy conditions, but impedes those choices in disease states.

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.

Adaptive delivery of continuous and delayed feedback deep brain stimulation - a computational study

Adaptive deep brain stimulation (aDBS) is a closed-loop method, where high-frequency DBS is turned on and off according to a feedback signal, whereas conventional high-frequency DBS (cDBS) is delivered permanently. Using a computational model of subthalamic nucleus and external globus pallidus, we extend the concept of adaptive stimulation by adaptively controlling not only continuous, but also demand-controlled stimulation. Apart from aDBS and cDBS, we consider continuous pulsatile linear delayed feedback stimulation (cpLDF), specifically designed to induce desynchronization. Additionally, we combine adaptive on-off delivery with continuous delayed feedback modulation by introducing adaptive pulsatile linear delayed feedback stimulation (apLDF), where cpLDF is turned on and off using pre-defined amplitude thresholds. By varying the stimulation parameters of cDBS, aDBS, cpLDF, and apLDF we obtain optimal parameter ranges. We reveal a simple relation between the thresholds of the local field potential (LFP) for aDBS and apLDF, the extent of the stimulation-induced desynchronization, and the integral stimulation time required. We find that aDBS and apLDF can be more efficient in suppressing abnormal synchronization than continuous simulation. However, apLDF still remains more efficient and also causes a stronger reduction of the LFP beta burst length. Hence, adaptive on-off delivery may further improve the intrinsically demand-controlled pLDF.

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.

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

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

Neuron dynamics variability and anomalous phase synchronization of neural networks

Anomalous phase synchronization describes a synchronization phenomenon occurring even for the weakly coupled network and characterized by a non-monotonous dependence of the synchronization strength on the coupling strength. Its existence may support a theoretical framework to some neurological diseases, such as Parkinson's and some episodes of seizure behavior generated by epilepsy. Despite the success of controlling or suppressing the anomalous phase synchronization in neural networks applying external perturbations or inducing ambient changes, the origin of the anomalous phase synchronization as well as the mechanisms behind the suppression is not completely known. Here, we consider networks composed of N = 2000 coupled neurons in a small-world topology for two well known neuron models, namely, the Hodgkin-Huxley-like and the Hindmarsh-Rose models, both displaying the anomalous phase synchronization regime. We show that the anomalous phase synchronization may be related to the individual behavior of the coupled neurons; particularly, we identify a strong correlation between the behavior of the inter-bursting-intervals of the neurons, what we call neuron variability, to the ability of the network to depict anomalous phase synchronization. We corroborate the ideas showing that external perturbations or ambient parameter changes that eliminate anomalous phase synchronization and at the same time promote small changes in the individual dynamics of the neurons, such that an increasing individual variability of neurons implies a decrease of anomalous phase synchronization. Finally, we demonstrate that this effect can be quantified using a well known recurrence quantifier, the "determinism." Moreover, the results obtained by the determinism are based on only the mean field potential of the network, turning these measures more suitable to be used in experimental situations.

Multisite Delayed Feedback for Electrical Brain Stimulation

Demand-controlled deep brain stimulation (DBS) appears to be a promising approach for the treatment of Parkinson's disease (PD) as revealed by computational, pre-clinical and clinical studies. Stimulation delivery is adapted to brain activity, for example, to the amount of neuronal activity considered to be abnormal. Such a closed-loop stimulation setup might help to reduce the amount of stimulation current, thereby maintaining therapeutic efficacy. In the context of the development of stimulation techniques that aim to restore desynchronized neuronal activity on a long-term basis, specific closed-loop stimulation protocols were designed computationally. These feedback techniques, e.g., pulsatile linear delayed feedback (LDF) or pulsatile nonlinear delayed feedback (NDF), were computationally developed to counteract abnormal neuronal synchronization characteristic for PD and other neurological disorders. By design, these techniques are intrinsically demand-controlled methods, where the amplitude of the stimulation signal is reduced when the desired desynchronized regime is reached. We here introduce a novel demand-controlled stimulation method, pulsatile multisite linear delayed feedback (MLDF), by employing MLDF to modulate the pulse amplitude of high-frequency (HF) DBS, in this way aiming at a specific, MLDF-related desynchronizing impact, while maintaining safety requirements with the charge-balanced HF DBS. Previously, MLDF was computationally developed for the control of spatio-temporal synchronized patterns and cluster states in neuronal populations. Here, in a physiologically motivated model network comprising neurons from subthalamic nucleus (STN) and external globus pallidus (GPe), we compare pulsatile MLDF to pulsatile LDF for the case where the smooth feedback signals are used to modulate the amplitude of charge-balanced HF DBS and suggest a modification of pulsatile MLDF which enables a pronounced desynchronizing impact. Our results may contribute to further clinical development of closed-loop DBS techniques.

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.

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.

Potential Mechanisms and Functions of Intermittent Neural Synchronization

Neural synchronization is believed to play an important role in different brain functions. Synchrony in cortical and subcortical circuits is frequently variable in time and not perfect. Few long intervals of desynchronized dynamics may be functionally different from many short desynchronized intervals although the average synchrony may be the same. Recent analysis of imperfect synchrony in different neural systems reported one common feature: neural oscillations may go out of synchrony frequently, but primarily for a short time interval. This study explores potential mechanisms and functional advantages of this short desynchronizations dynamics using computational neuroscience techniques. We show that short desynchronizations are exhibited in coupled neurons if their delayed rectifier potassium current has relatively large values of the voltage-dependent activation time-constant. The delayed activation of potassium current is associated with generation of quickly-rising action potential. This “spikiness” is a very general property of neurons. This may explain why very different neural systems exhibit short desynchronization dynamics. We also show how the distribution of desynchronization durations may be independent of the synchronization strength. Finally, we show that short desynchronization dynamics requires weaker synaptic input to reach a pre-set synchrony level. Thus, this dynamics allows for efficient regulation of synchrony and may promote efficient formation of synchronous neural assemblies.

Closed-loop deep brain stimulation by pulsatile delayed feedback with increased gap between pulse phases

Computationally it was shown that desynchronizing delayed feedback stimulation methods are effective closed-loop techniques for the control of synchronization in ensembles of interacting oscillators. We here computationally design stimulation signals for electrical stimulation of neuronal tissue that preserve the desynchronizing delayed feedback characteristics and comply with mandatory charge deposit-related safety requirements. For this, the amplitude of the high-frequency (HF) train of biphasic charge-balanced pulses used by the standard HF deep brain stimulation (DBS) is modulated by the smooth feedback signals. In this way we combine the desynchronizing delayed feedback approach with the HF DBS technique. We show that such a pulsatile delayed feedback stimulation can effectively and robustly desynchronize a network of model neurons comprising subthalamic nucleus and globus pallidus external and suggest this approach for desynchronizing closed-loop DBS. Intriguingly, an interphase gap introduced between the recharging phases of the charge-balanced biphasic pulses can significantly improve the stimulation-induced desynchronization and reduce the amount of the administered stimulation. In view of the recent experimental and clinical studies indicating a superiority of the closed-loop DBS to open-loop HF DBS, our results may contribute to a further development of effective stimulation methods for the treatment of neurological disorders characterized by abnormal neuronal synchronization.

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.

Pulsatile desynchronizing delayed feedback for closed-loop deep brain stimulation

High-frequency (HF) deep brain stimulation (DBS) is the gold standard for the treatment of medically refractory movement disorders like Parkinson’s disease, essential tremor, and dystonia, with a significant potential for application to other neurological diseases. The standard setup of HF DBS utilizes an open-loop stimulation protocol, where a permanent HF electrical pulse train is administered to the brain target areas irrespectively of the ongoing neuronal dynamics. Recent experimental and clinical studies demonstrate that a closed-loop, adaptive DBS might be superior to the open-loop setup. We here combine the notion of the adaptive high-frequency stimulation approach, that aims at delivering stimulation adapted to the extent of appropriately detected biomarkers, with specifically desynchronizing stimulation protocols. To this end, we extend the delayed feedback stimulation methods, which are intrinsically closed-loop techniques and specifically designed to desynchronize abnormal neuronal synchronization, to pulsatile electrical brain stimulation. We show that permanent pulsatile high-frequency stimulation subjected to an amplitude modulation by linear or nonlinear delayed feedback methods can effectively and robustly desynchronize a STN-GPe network of model neurons and suggest this approach for desynchronizing closed-loop DBS.

Synchronized Beta-Band Oscillations in a Model of the Globus Pallidus-Subthalamic Nucleus Network under External Input

Hypokinetic symptoms of Parkinson's disease are usually associated with excessively strong oscillations and synchrony in the beta frequency band. The origin of this synchronized oscillatory dynamics is being debated. Cortical circuits may be a critical source of excessive beta in Parkinson's disease. However, subthalamo-pallidal circuits were also suggested to be a substantial component in generation and/or maintenance of Parkinsonian beta activity. Here we study how the subthalamo-pallidal circuits interact with input signals in the beta frequency band, representing cortical input. We use conductance-based models of the subthalamo-pallidal network and two types of input signals: artificially-generated inputs and input signals obtained from recordings in Parkinsonian patients. The resulting model network dynamics is compared with the dynamics of the experimental recordings from patient's basal ganglia. Our results indicate that the subthalamo-pallidal model network exhibits multiple resonances in response to inputs in the beta band. For a relatively broad range of network parameters, there is always a certain input strength, which will induce patterns of synchrony similar to the experimentally observed ones. This ability of the subthalamo-pallidal network to exhibit realistic patterns of synchronous oscillatory activity under broad conditions may indicate that these basal ganglia circuits are directly involved in the expression of Parkinsonian synchronized beta oscillations. Thus, Parkinsonian synchronized beta oscillations may be promoted by the simultaneous action of both cortical (or some other) and subthalamo-pallidal network mechanisms. Hence, these mechanisms are not necessarily mutually exclusive.

Probing the Role of Medication, DBS Electrode Position, and Antidromic Activation on Impulsivity Using a Computational Model of Basal Ganglia

Everyday, we encounter situations where available choices are nearly equally rewarding (high conflict) calling for some tough decision making. Experimental recordings showed that the activity of Sub Thalamic Nucleus (STN) increases during such situations providing the extra time needed to make the right decision, teasing apart the most rewarding choice from the runner up closely trailing behind. This prolonged deliberation necessary for decision making under high conflict was absent in Parkinson's disease (PD) patients who underwent Deep Brain Stimulation (DBS) surgery of STN. In an attempt to understand the underlying cause of such adverse response, we built a 2D spiking network model (50 × 50 lattice) of Basal ganglia incorporating the key nuclei. Using the model we studied the Probabilistic learning task (PLT) in untreated, treated (L-Dopa and Dopamine Agonist) and STN-DBS PD conditions. Based on the experimental observation that dopaminergic activity is analogous to temporal difference (TD) and induces cortico-striatal plasticity, we introduced learning in the cortico-striatal weights. The results show that healthy and untreated conditions of PD model were able to more or less equally select (avoid) the rewarding (punitive) choice, a behavior that was absent in treated PD condition. The time taken to select a choice in high conflict trials was high in normal condition, which is in agreement with experimental results. The treated PD (Dopamine Agonist) patients made impulsive decisions (small reaction time) which in turn led to poor performance. The underlying cause of the observed impulsivity in DBS patients was studied in the model by (1) varying the electrode position within STN, (2) causing antidromic activation of GPe neurons. The effect of electrode position on reaction time was analyzed by studying the activity of STN neurons where, a decrease in STN neural activity was observed for certain electrode positions. We also observed that a higher antidromic activation of GPe neurons does not impact the learning ability but decreases reaction time as reported in DBS patients. These results suggest a probable role of electrode and antidromic activation in modulating the STN activity and eventually affecting the patient's performance on PLT.

Anti-kindling Induced by Two-Stage Coordinated Reset Stimulation with Weak Onset Intensity

Abnormal neuronal synchrony plays an important role in a number of brain diseases. To specifically counteract abnormal neuronal synchrony by desynchronization, Coordinated Reset (CR) stimulation, a spatiotemporally patterned stimulation technique, was designed with computational means. In neuronal networks with spike timing–dependent plasticity CR stimulation causes a decrease of synaptic weights and finally anti-kindling, i.e., unlearning of abnormally strong synaptic connectivity and abnormal neuronal synchrony. Long-lasting desynchronizing aftereffects of CR stimulation have been verified in pre-clinical and clinical proof of concept studies. In general, for different neuromodulation approaches, both invasive and non-invasive, it is desirable to enable effective stimulation at reduced stimulation intensities, thereby avoiding side effects. For the first time, we here present a two-stage CR stimulation protocol, where two qualitatively different types of CR stimulation are delivered one after another, and the first stage comes at a particularly weak stimulation intensity. Numerical simulations show that a two-stage CR stimulation can induce the same degree of anti-kindling as a single-stage CR stimulation with intermediate stimulation intensity. This stimulation approach might be clinically beneficial in patients suffering from brain diseases characterized by abnormal neuronal synchrony where a first treatment stage should be performed at particularly weak stimulation intensities in order to avoid side effects. This might, e.g., be relevant in the context of acoustic CR stimulation in tinnitus patients with hyperacusis or in the case of electrical deep brain CR stimulation with sub-optimally positioned leads or side effects caused by stimulation of the target itself. We discuss how to apply our method in first in man and proof of concept studies.

Separation of coexisting dynamical regimes in multistate intermittency based on wavelet spectrum energies in an erbium-doped fiber laser

We propose a method for the detection and localization of different types of coexisting oscillatory regimes that alternate with each other leading to multistate intermittency. Our approach is based on consideration of wavelet spectrum energies. The proposed technique is tested in an erbium-doped fiber laser with four coexisting periodic orbits, where external noise induces intermittent switches between the coexisting states. Statistical characteristics of multistate intermittency, such as the mean duration of the phases for every oscillation type, are examined with the help of the developed method. We demonstrate strong advantages of the proposed technique over previously used amplitude methods.

Coexistence of intermittencies in the neuronal network of the epileptic brain

Intermittent behavior occurs widely in nature. At present, several types of intermittencies are known and well-studied. However, consideration of intermittency has usually been limited to the analysis of cases when only one certain type of intermittency takes place. In this paper, we report on the temporal behavior of the complex neuronal network in the epileptic brain, when two types of intermittent behavior coexist and alternate with each other. We prove the presence of this phenomenon in physiological experiments with WAG/Rij rats being the model living system of absence epilepsy. In our paper, the deduced theoretical law for distributions of the lengths of laminar phases prescribing the power law with a degree of -2 agrees well with the experimental neurophysiological data.

Closed-loop stimulation of a delayed neural fields model of parkinsonian STN-GPe network: a theoretical and computational study

Several disorders are related to pathological brain oscillations. In the case of Parkinson's disease, sustained low-frequency oscillations (especially in the β-band, 13-30 Hz) correlate with motor symptoms. It is still under debate whether these oscillations are the cause of parkinsonian motor symptoms. The development of techniques enabling selective disruption of these β-oscillations could contribute to the understanding of the underlying mechanisms, and could be exploited for treatments. A particularly appealing technique is Deep Brain Stimulation (DBS). With clinical electrical DBS, electrical currents are delivered at high frequency to a region made of potentially heterogeneous neurons (the subthalamic nucleus (STN) in the case of Parkinson's disease). Even more appealing is DBS with optogenetics, which is until now a preclinical method using both gene transfer and deep brain light delivery and enabling neuromodulation at the scale of one given neural network. In this work, we rely on delayed neural fields models of STN and the external Globus Pallidus (GPe) to develop, theoretically validate and test in silico a closed-loop stimulation strategy to disrupt these sustained oscillations with optogenetics. First, we rely on tools from control theory to provide theoretical conditions under which sustained oscillations can be attenuated by a closed-loop stimulation proportional to the measured activity of STN. Second, based on this theoretical framework, we show numerically that the proposed closed-loop stimulation efficiently attenuates sustained oscillations, even in the case when the photosensitization effectively affects only 50% of STN neurons. We also show through simulations that oscillations disruption can be achieved when the same light source is used for the whole STN population. We finally test the robustness of the proposed strategy to possible acquisition and processing delays, as well as parameters uncertainty.

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

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

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.

Designing a deep brain stimulator to suppress pathological neuronal synchrony

Some of neuropathologies are believed to be related to abnormal synchronization of neurons. In the line of therapy, designing effective deep brain stimulators to suppress the pathological synchrony among neuronal ensembles is a challenge of high clinical relevance. The stimulation should be able to disrupt the synchrony in the presence of latencies due to imperfect knowledge about parameters of a neuronal ensemble and stimulation impacts on the ensemble. We propose an adaptive desynchronizing deep brain stimulator capable of dealing with these uncertainties. We analyze the collective behavior of the stimulated neuronal ensemble and show that, using the designed stimulator, the resulting asynchronous state is stable. Simulation results reveal the efficiency of the proposed technique.

Effects of reciprocal inhibitory coupling in model neurons

Central pattern generators are neuron networks that produce vital rhythmic motor outputs such as those observed in mastication, walking and breathing. Their activity patterns depend on the tuning of their intrinsic ionic conductances, their synaptic interconnectivity and entrainment by extrinsic neurons. The influence of two commonly found synaptic connectivities--reciprocal inhibition and electrical coupling--are investigated here using a neuron model with subthreshold oscillation capability, in different firing and entrainment regimes. We study the dynamics displayed by a network of a pair of neurons with various firing regimes, coupled by either (i) only reciprocal inhibition or by (ii) electrical coupling first and then reciprocal inhibition. In both scenarios a range of coupling strengths for the reciprocal inhibition is tested, and in general the neuron with the lower firing rate stops spiking for strong enough inhibitory coupling, while the faster neuron remains active. However, in scenario (ii) the originally slower neuron stops spiking at weaker inhibitory coupling strength, suggesting that the electrical coupling introduces an element of instability to the two-neuron network.

Parameter estimation and control for a neural mass model based on the unscented Kalman filter

Recent progress in Kalman filters to estimate states and parameters in nonlinear systems has provided the possibility of applying such approaches to neural systems. We here apply the nonlinear method of unscented Kalman filters (UKFs) to observe states and estimate parameters in a neural mass model that can simulate distinct rhythms in electroencephalography (EEG) including dynamical evolution during epilepsy seizures. We demonstrate the efficiency of the UKF in estimating states and parameters. We also develop an UKF-based control strategy to modulate the dynamics of the neural mass model. In this strategy the UKF plays the role of observing states, and the control law is constructed via the estimated states. We demonstrate the feasibility of using such a strategy to suppress epileptiform spikes in the neural mass model.

Cortical information flow in Parkinson's disease: a composite network/field model

The basal ganglia play a crucial role in the execution of movements, as demonstrated by the severe motor deficits that accompany Parkinson's disease (PD). Since motor commands originate in the cortex, an important question is how the basal ganglia influence cortical information flow, and how this influence becomes pathological in PD. To explore this, we developed a composite neuronal network/neural field model. The network model consisted of 4950 spiking neurons, divided into 15 excitatory and inhibitory cell populations in the thalamus and cortex. The field model consisted of the cortex, thalamus, striatum, subthalamic nucleus, and globus pallidus. Both models have been separately validated in previous work. Three field models were used: one with basal ganglia parameters based on data from healthy individuals, one based on data from individuals with PD, and one purely thalamocortical model. Spikes generated by these field models were then used to drive the network model. Compared to the network driven by the healthy model, the PD-driven network had lower firing rates, a shift in spectral power toward lower frequencies, and higher probability of bursting; each of these findings is consistent with empirical data on PD. In the healthy model, we found strong Granger causality between cortical layers in the beta and low gamma frequency bands, but this causality was largely absent in the PD model. In particular, the reduction in Granger causality from the main “input” layer of the cortex (layer 4) to the main “output” layer (layer 5) was pronounced. This may account for symptoms of PD that seem to reflect deficits in information flow, such as bradykinesia. In general, these results demonstrate that the brain's large-scale oscillatory environment, represented here by the field model, strongly influences the information processing that occurs within its subnetworks. Hence, it may be preferable to drive spiking network models with physiologically realistic inputs rather than pure white noise.

Short desynchronization episodes prevail in synchronous dynamics of human brain rhythms

Neural synchronization is believed to be critical for many brain functions. It frequently exhibits temporal variability, but it is not known if this variability has a specific temporal patterning. This study explores these synchronization/desynchronization patterns. We employ recently developed techniques to analyze the fine temporal structure of phase-locking to study the temporal patterning of synchrony of the human brain rhythms. We study neural oscillations recorded by electroencephalograms in α and β frequency bands in healthy human subjects at rest and during the execution of a task. While the phase-locking strength depends on many factors, dynamics of synchrony has a very specific temporal pattern: synchronous states are interrupted by frequent, but short desynchronization episodes. The probability for a desynchronization episode to occur decreased with its duration. The transition matrix between synchronized and desynchronized states has eigenvalues close to 0 and 1 where eigenvalue 1 has multiplicity 1, and therefore if the stationary distribution between these states is perturbed, the system converges back to the stationary distribution very fast. The qualitative similarity of this patterning across different subjects, brain states and electrode locations suggests that this may be a general type of dynamics for the brain. Earlier studies indicate that not all oscillatory networks have this kind of patterning of synchronization/desynchronization dynamics. Thus, the observed prevalence of short (but potentially frequent) desynchronization events (length of one cycle of oscillations) may have important functional implications for the brain. Numerous short desynchronizations (as opposed to infrequent, but long desynchronizations) may allow for a quick and efficient formation and break-up of functionally significant neuronal assemblies.

Failure of delayed feedback deep brain stimulation for intermittent pathological synchronization in Parkinson's disease

Suppression of excessively synchronous beta-band oscillatory activity in the brain is believed to suppress hypokinetic motor symptoms of Parkinson's disease. Recently, a lot of interest has been devoted to desynchronizing delayed feedback deep brain stimulation (DBS). This type of synchrony control was shown to destabilize the synchronized state in networks of simple model oscillators as well as in networks of coupled model neurons. However, the dynamics of the neural activity in Parkinson's disease exhibits complex intermittent synchronous patterns, far from the idealized synchronous dynamics used to study the delayed feedback stimulation. This study explores the action of delayed feedback stimulation on partially synchronized oscillatory dynamics, similar to what one observes experimentally in parkinsonian patients. We employ a computational model of the basal ganglia networks which reproduces experimentally observed fine temporal structure of the synchronous dynamics. When the parameters of our model are such that the synchrony is unphysiologically strong, the feedback exerts a desynchronizing action. However, when the network is tuned to reproduce the highly variable temporal patterns observed experimentally, the same kind of delayed feedback may actually increase the synchrony. As network parameters are changed from the range which produces complete synchrony to those favoring less synchronous dynamics, desynchronizing delayed feedback may gradually turn into synchronizing stimulation. This suggests that delayed feedback DBS in Parkinson's disease may boost rather than suppress synchronization and is unlikely to be clinically successful. The study also indicates that delayed feedback stimulation may not necessarily exhibit a desynchronization effect when acting on a physiologically realistic partially synchronous dynamics, and provides an example of how to estimate the stimulation effect.

Potential mechanisms for imperfect synchronization in parkinsonian basal ganglia

Neural activity in the brain of parkinsonian patients is characterized by the intermittently synchronized oscillatory dynamics. This imperfect synchronization, observed in the beta frequency band, is believed to be related to the hypokinetic motor symptoms of the disorder. Our study explores potential mechanisms behind this intermittent synchrony. We study the response of a bursting pallidal neuron to different patterns of synaptic input from subthalamic nucleus (STN) neuron. We show how external globus pallidus (GPe) neuron is sensitive to the phase of the input from the STN cell and can exhibit intermittent phase-locking with the input in the beta band. The temporal properties of this intermittent phase-locking show similarities to the intermittent synchronization observed in experiments. We also study the synchronization of GPe cells to synaptic input from the STN cell with dependence on the dopamine-modulated parameters. Earlier studies showed how the strengthening of dopamine-modulated coupling may lead to transitions from non-synchronized to partially synchronized dynamics, typical in Parkinson's disease. However, dopamine also affects the cellular properties of neurons. We show how the changes in firing patterns of STN neuron due to the lack of dopamine may lead to transition from a lower to a higher coherent state, roughly matching the synchrony levels observed in basal ganglia in normal and parkinsonian states. The intermittent nature of the neural beta band synchrony in Parkinson's disease is achieved in the model due to the interplay of the timing of STN input to pallidum and pallidal neuronal dynamics, resulting in sensitivity of pallidal output to the phase of the arriving STN input. Thus the mechanism considered here (the change in firing pattern of subthalamic neurons through the dopamine-induced change of membrane properties) may be one of the potential mechanisms responsible for the generation of the intermittent synchronization observed in Parkinson's disease.

Modeling effect of GABAergic current in a basal ganglia computational model

Electrical high frequency stimulation (HFS) of deep brain regions is a method shown to be clinically effective in different types of movement and neurological disorders. In order to shed light on its mode of action a computational model of the basal ganglia network coupled the HFS as injection current into the cells of the subthalamic nucleus (STN). Its overall increased activity rendered a faithful transmission of sensorimotor input through thalamo-cortical relay cells possible. Our contribution uses this model by Rubin and Terman (J Comput Neurosci, 16, 211-223, 2004) as a starting point and integrates recent findings on the importance of the extracellular concentrations of the inhibiting neurotransmitter GABA. We are able to show in this computational study that besides electrical stimulation a high concentration of GABA and its resulting conductivity in STN cells is able to re-establish faithful thalamocortical relaying, which otherwise broke down in the simulated parkinsonian state.

On the origin of tremor in Parkinson's disease

The exact origin of tremor in Parkinson’s disease remains unknown. We explain why the existing data converge on the basal ganglia-thalamo-cortical loop as a tremor generator and consider a conductance-based model of subthalamo-pallidal circuits embedded into a simplified representation of the basal ganglia-thalamo-cortical circuit to investigate the dynamics of this loop. We show how variation of the strength of dopamine-modulated connections in the basal ganglia-thalamo-cortical loop (representing the decreasing dopamine level in Parkinson’s disease) leads to the occurrence of tremor-like burst firing. These tremor-like oscillations are suppressed when the connections are modulated back to represent a higher dopamine level (as it would be the case in dopaminergic therapy), as well as when the basal ganglia-thalamo-cortical loop is broken (as would be the case for ablative anti-parkinsonian surgeries). Thus, the proposed model provides an explanation for the basal ganglia-thalamo-cortical loop mechanism of tremor generation. The strengthening of the loop leads to tremor oscillations, while the weakening or disconnection of the loop suppresses them. The loop origin of parkinsonian tremor also suggests that new tremor-suppression therapies may have anatomical targets in different cortical and subcortical areas as long as they are within the basal ganglia-thalamo-cortical loop.

Mathematical insights into the effects of levodopa

Parkinson’s disease has been traditionally thought of as a dopaminergic disease in which cells of the substantia nigra pars compacta (SNc) die. However, accumulating evidence implies an important role for the serotonergic system in Parkinson’s disease in general and in physiological responses to levodopa therapy, the first line of treatment. We use a mathematical model to investigate the consequences of levodopa therapy on the serotonergic system and on the pulsatile release of dopamine (DA) from dopaminergic and serotonergic terminals in the striatum. Levodopa competes with tyrosine and tryptophan at the blood-brain barrier and is taken up by serotonin neurons in which it competes for aromatic amino acid decarboxylase. The DA produced competes with serotonin (5HT) for packaging into vesicles. We predict the time courses of LD, cytosolic DA, and vesicular DA in 5HT neurons during an LD dose. We predict the time courses of DA and 5HT release from 5HT cell bodies and 5HT terminals as well as the changes in 5HT firing rate due to lower 5HT release. We compute the time course of DA release in the striatum from both 5HT and DA neurons and show how the time course changes as more and more SNc cells die. This enables us to explain the shortening of the therapeutic time window for the efficacy of levodopa as Parkinson’s disease progresses. Finally, we study the effects 5HT1a and 5HT1b autoreceptor agonists and explain why they have a synergistic effect and why they lengthen the therapeutic time window for LD therapy. Our results are consistent with and help explain results in the experimental literature and provide new predictions that can be tested experimentally.

Intermittent neural synchronization in Parkinson's disease

Motor symptoms of Parkinson's disease are related to the excessive synchronized oscillatory activity in the beta frequency band (around 20Hz) in the basal ganglia and other parts of the brain. This review explores the dynamics and potential mechanisms of these oscillations employing ideas and methods from nonlinear dynamics. We present extensive experimental documentation of the relevance of synchronized oscillations to motor behavior in Parkinson's disease, and we discuss the intermittent character of this synchronization. The reader is introduced to novel time-series analysis techniques aimed at the detection of the fine temporal structure of intermittent phase locking observed in the brains of parkinsonian patients. Modeling studies of brain networks are reviewed, which may describe the observed intermittent synchrony, and we discuss what these studies reveal about brain dynamics in Parkinson's disease. The parkinsonian brain appears to exist on the boundary between phase-locked and nonsynchronous dynamics. Such a situation may be beneficial in the healthy state, as it may allow for easy formation and dissociation of transient patterns of synchronous activity which are required for normal motor behavior. Dopaminergic degeneration in Parkinson's disease may shift the brain networks closer to this boundary, which would still permit some motor behavior while accounting for the associated motor deficits. Understanding the mechanisms of the intermittent synchrony in Parkinson's disease is also important for biomedical engineering since efficient control strategies for suppression of pathological synchrony through deep brain stimulation require knowledge of the dynamics of the processes subjected to control.

β band stability over time correlates with Parkinsonian rigidity and bradykinesia

Abnormal oscillatory activity in the basal ganglia is increasingly implicated in the pathophysiology of Parkinson's disease. Such activity is recorded in patients in the form of oscillations in the local field potential (LFP) picked up in the subthalamic nucleus. Previous studies have focused on correlations between features of the time averaged power or amplitude spectrum of the LFP and the clinical state, either off medication or in response to levodopa. However, average spectral densities do not take account of time variant spectral properties and we hypothesised that these dynamic properties of the spectrum of the LFP would contain additional information about clinical state. Here we assess the variability in LFP amplitude over time using the coefficient of variation (CV), evaluating this with regard to clinical state off medication and in response to levodopa in two datasets. The CV of activity in the high beta frequency band was found to be correlated with clinical state off levodopa (rho=-0.59, p<0.001) and this was shown to be complementary, rather than redundant, to spectral amplitude in a multiple regression analysis, selective for rigidity-bradykinesia and highly focal. Similarly, a strong correlation was found between change in clinical scores and change in high beta CV following levodopa (rho=-0.66, p=0.004). This too was selective for rigidity-bradykinesia and non-redundant to spectral power in a multiple regression model. Our results indicate that temporal stability in the beta band is correlated with rigidity-bradykinesia. It is suggested that loss of beta reactivity is deleterious to basal ganglia function over and above any concomitant change in absolute level of beta synchrony. The CV of LFP beta band amplitude may potentially provide an additional index of clinical state suitable for feedback control in closed loop stimulation therapy.

Intermittent synchronization in a network of bursting neurons

Synchronized oscillations in networks of inhibitory and excitatory coupled bursting neurons are common in a variety of neural systems from central pattern generators to human brain circuits. One example of the latter is the subcortical network of the basal ganglia, formed by excitatory and inhibitory bursters of the subthalamic nucleus and globus pallidus, involved in motor control and affected in Parkinson's disease. Recent experiments have demonstrated the intermittent nature of the phase-locking of neural activity in this network. Here, we explore one potential mechanism to explain the intermittent phase-locking in a network. We simplify the network to obtain a model of two inhibitory coupled elements and explore its dynamics. We used geometric analysis and singular perturbation methods for dynamical systems to reduce the full model to a simpler set of equations. Mathematical analysis was completed using three slow variables with two different time scales. Intermittently, synchronous oscillations are generated by overlapped spiking which crucially depends on the geometry of the slow phase plane and the interplay between slow variables as well as the strength of synapses. Two slow variables are responsible for the generation of activity patterns with overlapped spiking, and the other slower variable enhances the robustness of an irregular and intermittent activity pattern. While the analyzed network and the explored mechanism of intermittent synchrony appear to be quite generic, the results of this analysis can be used to trace particular values of biophysical parameters (synaptic strength and parameters of calcium dynamics), which are known to be impacted in Parkinson's disease.

Detecting the temporal structure of intermittent phase locking

This study explores a method to characterize the temporal structure of intermittent phase locking in oscillatory systems. When an oscillatory system is in a weakly synchronized regime away from a synchronization threshold, it spends most of the time in parts of its phase space away from the synchronization state. Therefore characteristics of dynamics near this state (such as its stability properties and Lyapunov exponents or distributions of the durations of synchronized episodes) do not describe the system's dynamics for most of the time. We consider an approach to characterize the system dynamics in this case by exploring the relationship between the phases on each cycle of oscillations. If some overall level of phase locking is present, one can quantify when and for how long phase locking is lost, and how the system returns back to the phase-locked state. We consider several examples to illustrate this approach: coupled skewed tent maps, the stability of which can be evaluated analytically; coupled Rössler and Lorenz oscillators, undergoing through different intermittency types on the way to phase synchronization; and a more complex example of coupled neurons. We show that the obtained measures can describe the differences in the dynamics and temporal structure of synchronization and desynchronization events for the systems with a similar overall level of phase locking and similar stability of the synchronized state.