This in silico study predicts the impact that the single-cell neuronal morphological alterations will have on the striatal microcircuit connectivity. We find that the richness in the topological striatal motifs is significantly reduced in Parkinson's disease (PD), highlighting that just measuring the pairwise connectivity between neurons gives an incomplete description of network connectivity. Moreover, we predict how the resulting electrophysiological changes of striatal projection neuron excitability together with their reduced number of dendritic branches affect their response to the glutamatergic drive from the cortex and thalamus. We find that the effective glutamatergic drive is likely significantly increased in PD, in accordance with the hyperglutamatergic hypothesis.

Biophysically detailed multi-compartment models are powerful tools to explore computational principles of the brain and also serve as a theoretical framework to generate algorithms for artificial intelligence (AI) systems. However, the expensive computational cost severely limits the applications in both the neuroscience and AI fields. The major bottleneck during simulating detailed compartment models is the ability of a simulator to solve large systems of linear equations. Here, we present a novel Dendritic Hierarchical Scheduling (DHS) method to markedly accelerate such a process. We theoretically prove that the DHS implementation is computationally optimal and accurate. This GPU-based method performs with 2-3 orders of magnitude higher speed than that of the classic serial Hines method in the conventional CPU platform. We build a DeepDendrite framework, which integrates the DHS method and the GPU computing engine of the NEURON simulator and demonstrate applications of DeepDendrite in neuroscience tasks. We investigate how spatial patterns of spine inputs affect neuronal excitability in a detailed human pyramidal neuron model with 25,000 spines. Furthermore, we provide a brief discussion on the potential of DeepDendrite for AI, specifically highlighting its ability to enable the efficient training of biophysically detailed models in typical image classification tasks.

The basal ganglia are important for action initiation, selection, and motor learning. The input level, the striatum, receives input preferentially from the cortex and thalamus and is to 95% composed of striatal projection neurons (SPNs) with sparse GABAergic collaterals targeting distal dendrites of neighboring SPNs, in a distance-dependent manner. The remaining 5% are GABAergic and cholinergic interneurons. Our aim here is to investigate the role of surround inhibition for the intrinsic function of the striatum. Large-scale striatal networks of 20 to 40 thousand neurons were simulated with detailed multicompartmental models of different cell types, corresponding to the size of a module of the dorsolateral striatum, like the forelimb area (mouse). The effect of surround inhibition on dendritic computation and network activity was investigated, while groups of SPNs were activated. The SPN-induced surround inhibition in distal dendrites shunted effectively the corticostriatal EPSPs. The size of dendritic plateau-like potentials within the specific dendritic segment was both reduced and enhanced by inhibition, due to the hyperpolarized membrane potential of SPNs and the reversal-potential of GABA. On a population level, the competition between two subpopulations of SPNs was found to depend on the distance between the two units, the size of each unit, the activity level in each subgroup and the dopaminergic modulation of the dSPNs and iSPNs. The SPNs provided the dominating source of inhibition within the striatum, while the fast-spiking interneuron mainly had an initial effect due to short-term synaptic plasticity as shown in with ablation of the synaptic interaction.

Neuromodulation is present throughout the nervous system and serves a critical role for circuit function and dynamics. The computational investigations of neuromodulation in large scale networks require supportive software platforms. Snudda is a software for the creation and simulation of large scale networks of detailed microcircuits consisting of multicompartmental neuron models. We have developed an extension to Snudda to incorporate neuromodulation in large scale simulations. The extended Snudda framework implements neuromodulation at the level of single cells incorporated into large-scale microcircuits. We also developed Neuromodcell, a software for optimizing neuromodulation in detailed multicompartmental neuron models. The software adds parameters within the models modulating the conductances of ion channels and ionotropic receptors. Bath application of neuromodulators is simulated and models which reproduce the experimentally measured effects are selected. In Snudda, we developed an extension to accommodate large scale simulations of neuromodulation. The simulator has two modes of simulation - denoted replay and adaptive. In the replay mode, transient levels of neuromodulators can be defined as a time-varying function which modulates the receptors and ion channels within the network in a cell-type specific manner. In the adaptive mode, spiking neuromodulatory neurons are connected via integrative modulating mechanisms to ion channels and receptors. Both modes of simulating neuromodulation allow for simultaneous modulation by several neuromodulators that can interact dynamically with each other. Here, we used the Neuromodcell software to simulate dopaminergic and muscarinic modulation of neurons from the striatum. We also demonstrate how to simulate different neuromodulatory states with dopamine and acetylcholine using Snudda. All software is freely available on Github, including tutorials on Neuromodcell and Snudda-neuromodulation.

Simulation of large-scale networks of neurons is an important approach to understanding and interpreting experimental data from healthy and diseased brains. Owing to the rapid development of simulation software and the accumulation of quantitative data of different neuronal types, it is possible to predict both computational and dynamical properties of local microcircuits in a ‘bottom-up’ manner. Simulated data from these models can be compared with experiments and ‘top-down’ modelling approaches, successively bridging the scales. Here we describe an open source pipeline, using the software Snudda, for predicting microcircuit connectivity and for setting up simulations using the NEURON simulation environment in a reproducible way. We also illustrate how to further ‘curate’ data on single neuron morphologies acquired from public databases. This model building pipeline was used to set up a first version of a full-scale cellular level model of mouse dorsal striatum. Model components from that work are here used to illustrate the different steps that are needed when modelling subcortical nuclei, such as the basal ganglia.

The basal ganglia play an important role in decision making and selection of action primarily based on input from cortex, thalamus, and the dopamine system. Their main input structure, striatum, is central to this process. It consists of two types of projection neurons, together representing 95% of the neurons, and 5% of interneurons, among which are the cholinergic, fast-spiking, and low threshold-spiking subtypes. The membrane properties, somadendritic shape, and intrastriatal and extrastriatal synaptic interactions of these neurons are quite well described in the mouse, and therefore they can be simulated in sufficient detail to capture their intrinsic properties, as well as the connectivity. We focus on simulation at the striatal cellular/microcircuit level, in which the molecular/subcellular and systems levels meet. We present a nearly full-scale model of the mouse striatum using available data on synaptic connectivity, cellular morphology, and electrophysiological properties to create a microcircuit mimicking the real network. A striatal volume is populated with reconstructed neuronal morphologies with appropriate cell densities, and then we connect neurons together based on appositions between neurites as possible synapses and constrain them further with available connectivity data. Moreover, we simulate a subset of the striatum involving 10,000 neurons, with input from cortex, thalamus, and the dopamine system, as a proof of principle. Simulation at this biological scale should serve as an invaluable tool to understand the mode of operation of this complex structure. This platform will be updated with new data and expanded to simulate the entire striatum.

Generation of neuronal cultures from induced pluripotent stem cells (hiPSCs) serve the studies of human brain disorders. However we lack neuronal networks with balanced excitatory- inhibitory activities, which are suitable for single cell analysis. We generated low-density networks of hPSC-derived GABAergic and glutamatergic cortical neurons. We used two different co-culture models with astrocytes. We show that these cultures have balanced excitatory-inhibitory synaptic identities using confocal microscopy, electrophysiological recordings, calcium imaging and mRNA analysis. These simple and robust protocols offer the opportunity for single-cell to multi-level analysis of patient hiPSC-derived cortical excitatory- inhibitory networks; thereby creating advanced tools to study disease mechanisms underlying neurodevelopmental disorders.

In the striatal microcircuit, fast-spiking (FS) interneurons have an important role in mediating inhibition onto neighboring medium spiny (MS) projection neurons. In this study, we combined computational modeling with in vitro and in vivo electrophysiological measurements to investigate FS cells in terms of their discharge properties and their synaptic efficacies onto MS neurons. In vivo firing of striatal FS interneurons is characterized by a high firing variability. It is not known, however, if this variability results from the input that FS cells receive, or if it is promoted by the stuttering spike behavior of these neurons. Both our model and measurements in vitro show that FS neurons that exhibit random stuttering discharge in response to steady depolarization do not show the typical stuttering behavior when they receive fluctuating input. Importantly, our model predicts that electrically coupled FS cells show substantial spike synchronization only when they are in the stuttering regime. Therefore, together with the lack of synchronized firing of striatal FS interneurons that has been reported in vivo, these results suggest that neighboring FS neurons are not in the stuttering regime simultaneously and that in vivo FS firing variability is more likely determined by the input fluctuations. Furthermore, the variability in FS firing is translated to variability in the postsynaptic amplitudes in MS neurons due to the strong synaptic depression of the FS-to-MS synapse. Our results support the idea that these synapses operate over a wide range from strongly depressed to almost fully recovered. The strong inhibitory effects that FS cells can impose on their postsynaptic targets, and the fact that the FS-to-MS synapse model showed substantial depression over extended periods of time might indicate the importance of cooperative effects of multiple presynaptic FS interneurons and the precise orchestration of their activity.

The intrastriatal microcircuit is a predominantly inhibitory GABAergic network comprised of a majority of projection neurons [medium spiny neurons (MSNs)] and a minority of interneurons. The connectivity within this microcircuit is divided into two main categories: lateral connectivity between MSNs, and inhibition mediated by interneurons, in particular fast spiking (FS) cells. To understand the operation of striatum, it is essential to have a good description of the dynamic properties of these respective pathways and how they affect different types of striatal projection neurons. We recorded from neuronal pairs, triplets, and quadruplets in slices of rat and mouse striatum and analyzed the dynamics of synaptic transmission between MSNs and FS cells. Retrograde fluorescent labeling and transgenic EGFP (enhanced green fluorescent protein) mice were used to distinguish between MSNs of the direct (striatonigral) and indirect (striatopallidal) pathways. Presynaptic neurons were stimulated with trains of action potentials, and activity-dependent depression and facilitation of synaptic efficacy was recorded from postsynaptic neurons. We found that FS cells provide a strong and homogeneously depressing inhibition of both striatonigral and striatopallidal MSN types. Moreover, individual FS cells are connected to MSNs of both types. In contrast, both MSN types receive sparse and variable, depressing and facilitating synaptic transmission from nearby MSNs. The connection probability was higher for pairs with presynaptic striatopallidal MSNs; however, the variability in synaptic dynamics did not depend on the types of interconnected MSNs. The differences between the two inhibitory pathways were clear in both species and at different developmental stages. Our findings show that the two intrastriatal inhibitory pathways have fundamentally different dynamic properties that are, however, similarly applied to both direct and indirect striatal projections.

MUSIC is an API allowing large scale neuron simulators using MPI internally to exchange data during runtime. We provide experiences from the adaptation of two neuronal network simulators of different kinds, NEST and MOOSE, to this API. A multi-simulation of a cortico-striatal network model involving both simulators is performed, demonstrating how MUSIC can promote inter-operability between models written for different simulators and how these can be re-used to build a larger model system. We conclude that MUSIC fulfills the design goals of being portable and simple to adapt to existing simulators. In addition, since the MUSIC API enforces independence between the applications, the multi-simulationcould be built from pluggable component modules without adaptation of the components to each other in terms of simulation time-step or topology of connections between the modules.