Virtual deep brain stimulation: Multiscale co-simulation of a spiking basal ganglia model and a whole-brain mean-field model with The Virtual Brain

Deep brain stimulation (DBS) has been successfully applied in various neurodegenerative diseases as an effective symptomatic treatment. However, its mechanisms of action within the brain network are still poorly understood. Many virtual DBS models analyze a subnetwork around the basal ganglia and its dynamics as a spiking network with their details validated by experimental data. However, connectomic evidence shows widespread effects of DBS affecting many different cortical and subcortical areas. From a clinical perspective, various effects of DBS besides the motoric impact have been demonstrated. The neuroinformatics platform The Virtual Brain (TVB) offers a modeling framework allowing us to virtually perform stimulation, including DBS, and forecast the outcome from a dynamic systems perspective prior to invasive surgery with DBS lead placement. For an accurate prediction of the effects of DBS, we implement a detailed spiking model of the basal ganglia, which we combine with TVB via our previously developed co-simulation environment. This multiscale co-simulation approach builds on the extensive previous literature of spiking models of the basal ganglia while simultaneously offering a whole-brain perspective on widespread effects of the stimulation going beyond the motor circuit. In the first demonstration of our model, we show that virtual DBS can move the firing rates of a Parkinson's disease patient's thalamus - basal ganglia network towards the healthy regime while, at the same time, altering the activity in distributed cortical regions with a pronounced effect in frontal regions. Thus, we provide proof of concept for virtual DBS in a co-simulation environment with TVB. The developed modeling approach has the potential to optimize DBS lead placement and configuration and forecast the success of DBS treatment for individual patients

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

Dopamine-modulated dynamic cell assemblies generated by the GABAergic striatal microcircuit

The striatum, the principal input structure of the basal ganglia, is crucial to both motor control and learning. It receives convergent input from all over the neocortex, hippocampal formation, amygdala and thalamus, and is the primary recipient of dopamine in the brain. Within the striatum is a GABAergic microcircuit that acts upon these inputs, formed by the dominant medium-spiny projection neurons (MSNs) and fast-spiking interneurons (FSIs). There has been little progress in understanding the computations it performs, hampered by the non-laminar structure that prevents identification of a repeating canonical microcircuit. We here begin the identification of potential dynamically-defined computational elements within the striatum. We construct a new three-dimensional model of the striatal microcircuit's connectivity, and instantiate this with our dopamine-modulated neuron models of the MSNs and FSIs. A new model of gap junctions between the FSIs is introduced and tuned to experimental data. We introduce a novel multiple spike-train analysis method, and apply this to the outputs of the model to find groups of synchronised neurons at multiple time-scales. We find that, with realistic in vivo background input, small assemblies of synchronised MSNs spontaneously appear, consistent with experimental observations, and that the number of assemblies and the time-scale of synchronisation is strongly dependent on the simulated concentration of dopamine. We also show that feed-forward inhibition from the FSIs counter-intuitively increases the firing rate of the MSNs. Such small cell assemblies forming spontaneously only in the absence of dopamine may contribute to motor control problems seen in humans and animals following a loss of dopamine cells. (C) 2009 Elsevier Ltd. All rights reserved

Transient and steady-state selection in the striatal microcircuit

Although the basal ganglia have been widely studied and implicated in signal processing and action selection, little information is known about the active role the striatal microcircuit plays in action selection in the basal ganglia-thalamo-cortical loops. To address this knowledge gap we use a large scale three dimensional spiking model of the striatum, combined with a rate coded model of the basal ganglia-thalamo-cortical loop, to asses the computational role the striatum plays in action selection. We identify a robust transient phenomena generated by the striatal microcircuit, which temporarily enhances the difference between two competing cortical inputs. We show that this transient is sufficient to modulate decision making in the basal ganglia-thalamo-cortical circuit. We also find that the transient selection originates from a novel adaptation effect in single striatal projection neurons, which is amenable to experimental testing. Finally, we compared transient selection with models implementing classical steady-state selection. We challenged both forms of model to account for recent reports of paradoxically enhanced response selection in Huntington's disease patients. We found that steady-state selection was uniformly impaired under all simulated Huntington's conditions, but transient selection was enhanced given a sufficient Huntington's-like increase in NMDA receptor sensitivity. Thus our models provide an intriguing hypothesis for the mechanisms underlying the paradoxical cognitive improvements in manifest Huntington's patients

Reconstructing the three-dimensional GABAergic microcircuit of the striatum

A system's wiring constrains its dynamics, yet modelling of neural structures often overlooks the specific networks formed by their neurons. We developed an approach for constructing anatomically realistic networks and reconstructed the GABAergic microcircuit formed by the medium spiny neurons (MSNs) and fast-spiking interneurons (FSIs) of the adult rat striatum. We grew dendrite and axon models for these neurons and extracted probabilities for the presence of these neurites as a function of distance from the soma. From these, we found the probabilities of intersection between the neurites of two neurons given their inter-somatic distance, and used these to construct three-dimensional striatal networks. The MSN dendrite models predicted that half of all dendritic spines are within 100 mu m of the soma. The constructed networks predict distributions of gap junctions between FSI dendrites, synaptic contacts between MSNs, and synaptic inputs from FSIs to MSNs that are consistent with current estimates. The models predict that to achieve this, FSIs should be at most 1% of the striatal population. They also show that the striatum is sparsely connected: FSI-MSN and MSN-MSN contacts respectively form 7% and 1.7% of all possible connections. The models predict two striking network properties: the dominant GABAergic input to a MSN arises from neurons with somas at the edge of its dendritic field; and FSIs are interconnected on two different spatial scales: locally by gap junctions and distally by synapses. We show that both properties influence striatal dynamics: the most potent inhibition of a MSN arises from a region of striatum at the edge of its dendritic field; and the combination of local gap junction and distal synaptic networks between FSIs sets a robust input-output regime for the MSN population. Our models thus intimately link striatal micro-anatomy to its dynamics, providing a biologically grounded platform for further study

Action selection in the striatum: Implications for Huntington's disease

Although the basal ganglia have been widely studied and implicated in signal processing and action selection, little information is known about the active role that the striatal microcircuit plays in action selection in the basal ganglia-cortical-thalamic loops. To address this knowledge gap we use a large scale three dimensional spiking model of the striatum, combined with a rate coded model of the basal ganglia-cortical-thalamic loop, to asses the computational role the striatum plays in action selection. We identify robust transient phenomena generated by the striatal microcircuit, which temporarily enhances the difference between two competing cortical inputs. We show that this transient is sufficient to modulate decision making in the basal ganglia-thalamo-cortical circuit. We also find that the transient selection originates from a novel adaptation effect in single striatal projection neurons, which is amenable to experimental testing. Finally, we compared transient selection with models implementing classical steady-state selection. We challenged both forms of model to account for recent reports of paradoxically enhanced response selection in Huntington's disease patients. We found that steady-state selection was uniformly impaired under all simulated Huntington's conditions, but transient selection was enhanced given a sufficient Huntington's-like increase in NMDA receptor sensitivity. I propose a mechanistic underpinning to a novel neural compensatory mechanism, responsible for improved cognition in severe neuro-degeneration. Thus, our models provide an intriguing hypothesis for the mechanisms underlying the paradoxical cognitive improvements in manifest Huntington's patients, which is consistent with recent behavioural data

Using Phase Response Curves to Optimize Deep Brain Stimulation

University of Minnesota Ph.D. dissertation. April 2016. Major: Neuroscience. Advisor: Theoden Netoff. 1 computer file (PDF); vii, 190 pages.Deep brain stimulation (DBS) is a neuromodulation therapy effective at treating motor symptoms of patients with Parkinson’s disease (PD). Currently, an open-loop approach is used to set stimulus parameters, where stimulation settings are programmed by a clinician using a time intensive trial-and-error process. There is a need for a systematic approach to tuning stimulation parameters based on a patient’s physiology. An effective biomarker in the recorded neural signal is needed for this approach. It is hypothesized that DBS may work by disrupting enhanced oscillatory activity seen in PD. In this thesis I propose and provide evidence for using a simple measure, called a phase response curve, to systematically tune stimulation parameters and develop novel approaches to stimulation to suppress pathological oscillations. In this work I show that PRCs can be used to optimize stimulus frequency, waveform, and stimulus phase to disrupt a pathological oscillation in a computational model of Parkinson’s disease and/or to disrupt entrainment of single neurons in vitro. This approach has the potential to improve efficacy and reduce post-operative programming time

A dynamical model for the basal ganglia-thalamo-cortical oscillatory activity and its implications in Parkinson’s disease

© 2020 The Authors. Published by Springer. This is an open access article available under a Creative Commons licence. The published version can be accessed at the following link on the publisher’s website: https://doi.org/10.1007/s11571-020-09653-y© 2020, The Author(s). We propose to investigate brain electrophysiological alterations associated with Parkinson’s disease through a novel adaptive dynamical model of the network of the basal ganglia, the cortex and the thalamus. The model uniquely unifies the influence of dopamine in the regulation of the activity of all basal ganglia nuclei, the self-organised neuronal interdependent activity of basal ganglia-thalamo-cortical circuits and the generation of subcortical background oscillations. Variations in the amount of dopamine produced in the neurons of the substantia nigra pars compacta are key both in the onset of Parkinson’s disease and in the basal ganglia action selection. We model these dopamine-induced relationships, and Parkinsonian states are interpreted as spontaneous emergent behaviours associated with different rhythms of oscillatory activity patterns of the basal ganglia-thalamo-cortical network. These results are significant because: (1) the neural populations are built upon single-neuron models that have been robustly designed to have eletrophysiologically-realistic responses, and (2) our model distinctively links changes in the oscillatory activity in subcortical structures, dopamine levels in the basal ganglia and pathological synchronisation neuronal patterns compatible with Parkinsonian states, this still remains an open problem and is crucial to better understand the progression of the disease.Published versio

Global neural rhythm control by local neuromodulation

Neural oscillations are a ubiquitous form of neural activity seen across scales and modalities. These neural rhythms correlate with diverse cognitive functions and brain states. One mechanism for changing the oscillatory dynamics of large neuronal populations is through neuromodulator activity. An intriguing phenomenon explored here is when local neuromodulation of a distinct neuron type within a single brain nucleus exerts a powerful influence on global cortical rhythms. One approach to investigate the impact of local circuits on global rhythms is through optogenetic techniques. My first project involves the statistical analysis of electrophysiological recordings of an optogenetically-mediated Parkinsonian phenotype. Empirical studies demonstrate that Parkinsonian motor deficits correlate with the emergence of exaggerated beta frequency (15-30 Hz) oscillations throughout the cortico-basal ganglia-thalamic network. However, the mechanism of these aberrant oscillatory dynamics is not well understood. A previous modeling study predicted that cholinergic neuromodulation of medium spiny neurons in the striatum of the basal ganglia may mediate the pathologic beta rhythm. Here, this hypothesis was tested using selective optogenetic stimulation of striatal cholinergic interneurons in normal mice; stimulation robustly and reversibly amplified beta oscillations and Parkinsonian motor symptoms. The modulation of global rhythms by local networks was further studied using computational modeling in the context of intrathalamic neuromodulation. While intrathalamic vasoactive intestinal peptide (VIP) is known to cause long-lasting excitation in vitro, its in vivo dynamical effects have not been reported. Here, biophysical computational models were used to elucidate the impact of VIP on thalamocortical dynamics during sleep and propofol general anesthesia. The modeling results suggest that VIP can form robust sleep spindle oscillations and control aspects of sleep architecture through a novel homeostatic mechanism. This homeostatic mechanism would be inhibited by general anesthesia, representing a new mechanism contributing to anesthetic-induced loss of consciousness. While the previous two projects differed in their use of empirical versus theoretical methods, a challenge common to both domains is the difficulty in visualizing and analyzing large multi-dimensional datasets. A tool to mitigate these issues is introduced here: GIMBL-Vis is a Graphical Interactive Multi-dimensional extensiBLe Visualization toolbox for Matlab. This toolbox simplifies the process of exploring multi-dimensional data in Matlab by providing a graphical interface for visualization and analysis. Furthermore, it provides an extensible open platform for distributed development by the community