A Population Model of Deep Brain Stimulation in Movement Disorders From Circuits to Cells

A Model with Dopamine Depletion in Basal Ganglia and Cerebellum Predicts Changes in Thalamocortical Beta Oscillations

Parkinsonism is presented as a motor syndrome characterized by rigidity, tremors, and bradykinesia, with Parkinson's disease (PD) being the predominant cause. The discovery that those motor symptoms result from the death of dopaminergic cells in the substantia nigra led to focus most of parkinsonism research on the basal ganglia (BG). However, recent findings point to an active involvement of the cerebellum in this motor syndrome. Here, we have developed a multiscale computational model of the rodent brain's BG-cerebellar network. Simulations showed that a direct effect of dopamine depletion on the cerebellum must be taken into account to reproduce the alterations of neural activity in parkinsonism, particularly the increased beta oscillations widely reported in PD patients. Moreover, dopamine depletion indirectly impacted spike-time-dependent plasticity at the parallel fiber-Purkinje cell synapses, degrading associative motor learning as observed in parkinsonism. Overall, these results suggest a relevant involvement of cerebellum in parkinsonism associative motor symptoms.

Modeling Instantaneous Firing Rate of Deep Brain Stimulation Target Neuronal Ensembles in the Basal Ganglia and Thalamus

OBJECTIVE

Deep brain stimulation (DBS) is an effective treatment for movement disorders, including Parkinson disease and essential tremor. However, the underlying mechanisms of DBS remain elusive. Despite the capability of existing models in interpreting experimental data qualitatively, there are very few unified computational models that quantitatively capture the dynamics of the neuronal activity of varying stimulated nuclei-including subthalamic nucleus (STN), substantia nigra pars reticulata (SNr), and ventral intermediate nucleus (Vim)-across different DBS frequencies.

MATERIALS AND METHODS

Both synthetic and experimental data were used in the model fitting; the synthetic data were generated by an established spiking neuron model that was reported in our previous work, and the experimental data were provided using single-unit microelectrode recordings (MERs) during DBS (microelectrode stimulation). Based on these data, we developed a novel mathematical model to represent the firing rate of neurons receiving DBS, including neurons in STN, SNr, and Vim-across different DBS frequencies. In our model, the DBS pulses were filtered through a synapse model and a nonlinear transfer function to formulate the firing rate variability. For each DBS-targeted nucleus, we fitted a single set of optimal model parameters consistent across varying DBS frequencies.

RESULTS

Our model accurately reproduced the firing rates observed and calculated from both synthetic and experimental data. The optimal model parameters were consistent across different DBS frequencies.

CONCLUSIONS

The result of our model fitting was in agreement with experimental single-unit MER data during DBS. Reproducing neuronal firing rates of different nuclei of the basal ganglia and thalamus during DBS can be helpful to further understand the mechanisms of DBS and to potentially optimize stimulation parameters based on their actual effects on neuronal activity.

Local field potential (LFP) power and phase-amplitude coupling (PAC) changes in the striatum and motor cortex reflect neural mechanisms associated with bradykinesia and rigidity during D2R suppression in an animal model

Impairments in motor control are the primary feature of Parkinson's disease, which is caused by dopaminergic imbalance in the basal ganglia. Identification of neural biomarkers of dopamine D2 receptor (D2R) suppression would be useful for monitoring the progress of neuropathologies and effects of treatment. Male Swiss albino ICR mice were deeply anesthetized, and electrodes were implanted in the striatum and motor cortex to record local field potential (LFP). Haloperidol (HAL), a D2R antagonist, was administered to induce decreased D2R activity. Following HAL treatment, the mice showed significantly decreased movement velocity in open field test, increased latency to descend in a bar test, and decreased latency to fall in a rotarod test. LFP signals during HAL-induced immobility (open field test) and catalepsy (bar test) were analyzed. Striatal low-gamma (30.3-44.9 Hz) power decreased during immobility periods, but during catalepsy, delta power (1-4 Hz) increased, beta1(13.6-18 Hz) and low-gamma powers decreased, and high-gamma (60.5-95.7 Hz) power increased. Striatal delta-high-gamma phase-amplitude coupling (PAC) was significantly increased during catalepsy but not immobility. In the motor cortex, during HAL-induced immobility, beta1 power significantly increased and low-gamma power decreased, but during HAL-induced catalepsy, low-gamma and beta1 powers decreased and high-gamma power increased. Delta-high-gamma PAC in the motor cortex significantly increased during catalepsy but not during immobility. Altogether, the present study demonstrated changes in delta, beta1 and gamma powers and delta-high-gamma PAC in the striatum and motor cortex in association with D2R suppression. In particular, delta power in the striatum and delta-high-gamma PAC in the striatum and motor cortex appear to represent biomarkers of neural mechanisms associated with bradykinesia and rigidity.

Computational neurosciences and quantitative systems pharmacology: a powerful combination for supporting drug development in neurodegenerative diseases

Successful clinical development of new therapeutic interventions is notoriously difficult, especially in neurodegenerative diseases, where predictive biomarkers are scarce and functional improvement is often based on patient's perception, captured by structured interviews. As a consequence, mechanistic modeling of the processes relevant to therapeutic interventions in CNS disorders has been lagging behind other disease indications, probably because of the perceived complexity of the brain. However in this report, we develop the argument that a combination of Computational Neurosciences and Quantitative Systems Pharmacology (QSP) modeling of molecular pathways is a powerful simulation tool to enhance the probability of successful drug development for neurodegenerative diseases. Computational Neurosciences aims to predict action potential dynamics and neuronal circuit activation that are ultimately linked to behavioral changes and clinically relevant functional outcomes. These processes can not only be affected by the disease state, but also by common genotype variants on neurotransmitter-related proteins and the psycho-active medications often prescribed in these patient populations. Quantitative Systems Pharmacology (QSP) modeling of molecular pathways allows to simulate key pathological drivers of dementia, such as protein aggregation and neuroinflammatory responses. They often impact neurotransmitter homeostasis and voltage-gated ion-channels or lead to mitochondrial dysfunction, ultimately leading to changes in action potential dynamics and clinical readouts. Combining these two modeling approaches can lead to better actionable understanding of the many non-linear pharmacodynamic processes active in the human diseased brain. Practical applications include a rational selection of the optimal doses in combination therapies, identification of subjects more likely to respond to treatment, a more balanced stratification of treatment arms in terms of comedications, disease status and common genotype variants and re-analysis of small clinical trials to uncover a possible clinical signal. Ultimately this will lead to a higher success rate of bringing new therapeutics to the right patient populations.

Subthalamic Nucleus (STN)-Deep Brain Stimulation Reduces the Power of Mu and Beta Rhythms and Enhances Synchrony at the Motor Cortices in Parkinson's Disease: A Report of Two Cases

The motor circuit in Parkinson's disease (PD) involves the basal ganglia, thalamus, motor cortex, and cerebellum. Hence, subthalamic nucleus (STN) or globus pallidus internus deep brain stimulation is commonly used in treating refractory Parkinson's patients. During the procedure, the local field potential (LPF) is commonly made along the trajectory of the STN. Two cases were assessed, where an electroencephalographic recording at the sensorimotor cortices was also performed with and without stimulation at the optimal STN electrode site. The 'on' stimulation state associated with clinical improvement correlated with a marked reduction in the late theta (7.5 Hz), alpha (10.5 Hz) (Mu wave), and beta (20 Hz) wave power. Besides, more synchronized and coherent brainwaves were noted when the stimulation was 'on'.

Dynamical Mechanism Analysis of Three Neuroregulatory Strategies on the Modulation of Seizures

This paper attempts to explore and compare the regulatory mechanisms of optogenetic stimulation (OS), deep brain stimulation (DBS) and electromagnetic induction on epilepsy. Based on the Wilson-Cowan model, we first demonstrate that the external input received by excitatory and inhibitory neural populations can induce rich dynamic bifurcation behaviors such as Hopf bifurcation, and make the system exhibit epileptic and normal states. Then, both OS and DBS are shown to be effective in controlling the epileptic state to a normal low-level state, and the stimulus parameters have a broad effective range. However, electromagnetic induction cannot directly control epilepsy to this desired state, even if it can significantly reduce the oscillation frequency of neural populations. One main difference worth noting is that the high spatiotemporal specificity of OS allows it to target inhibitory neuronal populations, whereas DBS and electromagnetic induction can only stimulate excitatory as well as inhibitory neuronal populations together. Next, the propagation behavior of epilepsy is explored under a typical three-node feedback loop structure. An increase in coupling strength accelerates and exacerbates epileptic activity in other brain regions. Finally, OS and DBS applied to the epileptic focus play similar positive roles in controlling the behavior of the area of seizure propagation, while electromagnetic induction still only achieves unsatisfactory effects. It is hoped that these dynamical results can provide insights into the treatment of epilepsy as well as other neurological disorders.

Hallmarks of neurodegenerative disease: A systems pharmacology perspective

Age-related central neurodegenerative diseases, such as Alzheimer's and Parkinson's disease, are a rising public health concern and have been plagued by repeated drug development failures. The complex nature and poor mechanistic understanding of the etiology of neurodegenerative diseases has hindered the discovery and development of effective disease-modifying therapeutics. Quantitative systems pharmacology models of neurodegeneration diseases may be useful tools to enhance the understanding of pharmacological intervention strategies and to reduce drug attrition rates. Due to the similarities in pathophysiological mechanisms across neurodegenerative diseases, especially at the cellular and molecular levels, we envision the possibility of structural components that are conserved across models of neurodegenerative diseases. Conserved structural submodels can be viewed as building blocks that are pieced together alongside unique disease components to construct quantitative systems pharmacology (QSP) models of neurodegenerative diseases. Model parameterization would likely be different between the different types of neurodegenerative diseases as well as individual patients. Formulating our mechanistic understanding of neurodegenerative pathophysiology as a mathematical model could aid in the identification and prioritization of drug targets and combinatorial treatment strategies, evaluate the role of patient characteristics on disease progression and therapeutic response, and serve as a central repository of knowledge. Here, we provide a background on neurodegenerative diseases, highlight hallmarks of neurodegeneration, and summarize previous QSP models of neurodegenerative diseases.

Portable Neuroimaging-Guided Noninvasive Brain Stimulation of the Cortico-Cerebello-Thalamo-Cortical Loop—Hypothesis and Theory in Cannabis Use Disorder

Background: Maladaptive neuroplasticity-related learned response in substance use disorder (SUD) can be ameliorated using noninvasive brain stimulation (NIBS); however, inter-individual variability needs to be addressed for clinical translation. Objective: Our first objective was to develop a hypothesis for NIBS for learned response in SUD based on a competing neurobehavioral decision systems model. The next objective was to develop the theory by conducting a computational simulation of NIBS of the cortico-cerebello-thalamo-cortical (CCTC) loop in cannabis use disorder (CUD)-related dysfunctional “cue-reactivity”—a construct closely related to “craving”—that is a core symptom. Our third objective was to test the feasibility of a neuroimaging-guided rational NIBS approach in healthy humans. Methods: “Cue-reactivity” can be measured using behavioral paradigms and portable neuroimaging, including functional near-infrared spectroscopy (fNIRS) and electroencephalogram (EEG) metrics of sensorimotor gating. Therefore, we conducted a computational simulation of NIBS, including transcranial direct current stimulation (tDCS) and transcranial alternating current stimulation (tACS) of the cerebellar cortex and deep cerebellar nuclei (DCN) of the CCTC loop for its postulated effects on fNIRS and EEG metrics. We also developed a rational neuroimaging-guided NIBS approach for the cerebellar lobule (VII) and prefrontal cortex based on a healthy human study. Results: Simulation of cerebellar tDCS induced gamma oscillations in the cerebral cortex, while transcranial temporal interference stimulation induced a gamma-to-beta frequency shift. A preliminary healthy human study (N = 10) found that 2 mA cerebellar tDCS evoked similar oxyhemoglobin (HbO) response in the range of 5 × 10⁻⁶ M across the cerebellum and PFC brain regions (α = 0.01); however, infra-slow (0.01–0.10 Hz) prefrontal cortex HbO-driven phase–amplitude-coupled (PAC; 4 Hz, ±2 mA (max)) cerebellar tACS evoked HbO levels in the range of 10⁻⁷ M that were statistically different (α = 0.01) across these brain regions. Conclusion: Our healthy human study showed the feasibility of fNIRS of cerebellum and PFC and closed-loop fNIRS-driven ctACS at 4 Hz, which may facilitate cerebellar cognitive function via the frontoparietal network. Future work needs to combine fNIRS with EEG for multi-modal imaging for closed-loop NIBS during operant conditioning.

Modeling Neurotransmission: Computational Tools to Investigate Neurological Disorders

The investigation of synaptic functions remains one of the most fascinating challenges in the field of neuroscience and a large number of experimental methods have been tuned to dissect the mechanisms taking part in the neurotransmission process. Furthermore, the understanding of the insights of neurological disorders originating from alterations in neurotransmission often requires the development of (i) animal models of pathologies, (ii) invasive tools and (iii) targeted pharmacological approaches. In the last decades, additional tools to explore neurological diseases have been provided to the scientific community. A wide range of computational models in fact have been developed to explore the alterations of the mechanisms involved in neurotransmission following the emergence of neurological pathologies. Here, we review some of the advancements in the development of computational methods employed to investigate neuronal circuits with a particular focus on the application to the most diffuse neurological disorders.