Motor cortex microcircuit simulation based on brain activity mapping

Deep Brain Stimulation restores information processing in parkinsonian cortical networks

Parkinson's disease (PD) is a neurodegenerative disorder associated with alterations of neural activity and information processing primarily in the basal ganglia and cerebral cortex. Deep brain stimulation (DBS) of the subthalamic nucleus (STN-DBS) is the most effective therapy when patients experience levodopa-induced motor complications. A growing body of evidence points towards a cortical effect of STN-DBS, restoring key electrophysiological markers, such as excessive beta band oscillations, commonly observed in PD. However, the mechanisms of STN-DBS remain elusive. Here, we aim to better characterize the cortical substrates underlying STN-DBS-induced improvement in motor symptoms. We recorded electroencephalograms (EEG) from PD patients and found that, although apparent EEG features were not different with or without therapy, EEG signals could more accurately predict limb movements under STN-DBS. To understand the origins of this enhanced information transmission under STN-DBS in the human EEG data, we investigated the information capacity and dynamics of a variety of computational models of cortical networks. The extent of improvement in decoding accuracy of complex naturalistic inputs under STN-DBS depended on the synaptic parameters of the network as well as its excitability and synchronization levels. Additionally, decoding accuracy could be optimized by adjusting STN-DBS parameters. Altogether, this work draws a comprehensive link between known alterations in cortical activity and the degradation of information processing capacity, as well as its restoration under DBS. These results also offer new perspectives for optimizing STN-DBS parameters based on clinically accessible measures of cortical information processing capacity.

Multiscale model of primary motor cortex circuits predicts in vivo cell-type-specific, behavioral state-dependent dynamics

Understanding cortical function requires studying multiple scales: molecular, cellular, circuit, and behavioral. We develop a multiscale, biophysically detailed model of mouse primary motor cortex (M1) with over 10,000 neurons and 30 million synapses. Neuron types, densities, spatial distributions, morphologies, biophysics, connectivity, and dendritic synapse locations are constrained by experimental data. The model includes long-range inputs from seven thalamic and cortical regions and noradrenergic inputs. Connectivity depends on cell class and cortical depth at sublaminar resolution. The model accurately predicts in vivo layer- and cell-type-specific responses (firing rates and LFP) associated with behavioral states (quiet wakefulness and movement) and experimental manipulations (noradrenaline receptor blockade and thalamus inactivation). We generate mechanistic hypotheses underlying the observed activity and analyzed low-dimensional population latent dynamics. This quantitative theoretical framework can be used to integrate and interpret M1 experimental data and sheds light on the cell-type-specific multiscale dynamics associated with several experimental conditions and behaviors.

Ultrafast simulation of large-scale neocortical microcircuitry with biophysically realistic neurons

Understanding the activity of the mammalian brain requires an integrative knowledge of circuits at distinct scales, ranging from ion channel gating to circuit connectomics. Computational models are regularly employed to understand how multiple parameters contribute synergistically to circuit behavior. However, traditional models of anatomically and biophysically realistic neurons are computationally demanding, especially when scaled to model local circuits. To overcome this limitation, we trained several artificial neural network (ANN) architectures to model the activity of realistic multicompartmental cortical neurons. We identified an ANN architecture that accurately predicted subthreshold activity and action potential firing. The ANN could correctly generalize to previously unobserved synaptic input, including in models containing nonlinear dendritic properties. When scaled, processing times were orders of magnitude faster compared with traditional approaches, allowing for rapid parameter-space mapping in a circuit model of Rett syndrome. Thus, we present a novel ANN approach allowing for rapid, detailed network experiments using inexpensive and commonly available computational resources.

Deep brain stimulation-guided optogenetic rescue of parkinsonian symptoms

Deep brain stimulation (DBS) of the subthalamic nucleus is a symptomatic treatment of Parkinson’s disease but benefits only to a minority of patients due to stringent eligibility criteria. To investigate new targets for less invasive therapies, we aimed at elucidating key mechanisms supporting deep brain stimulation efficiency. Here, using in vivo electrophysiology, optogenetics, behavioral tasks and mathematical modeling, we found that subthalamic stimulation normalizes pathological hyperactivity of motor cortex pyramidal cells, while concurrently activating somatostatin and inhibiting parvalbumin interneurons. In vivo opto-activation of cortical somatostatin interneurons alleviates motor symptoms in a parkinsonian mouse model. A computational model highlights that a decrease in pyramidal neuron activity induced by DBS or by a stimulation of cortical somatostatin interneurons can restore information processing capabilities. Overall, these results demonstrate that activation of cortical somatostatin interneurons may constitute a less invasive alternative than subthalamic stimulation.

Deep brain stimulation (DBS) is a symptomatic treatment of Parkinson’s disease (PD) that benefits only a minority of patients. Here, the authors show that activation of cortical somatostatin interneurons alleviates motor symptoms in a mouse model of PD and may constitute a less invasive alternative than DBS.

Relating network connectivity to dynamics: opportunities and challenges for theoretical neuroscience

We review recent work relating network connectivity to the dynamics of neural activity. While concepts stemming from network science provide a valuable starting point, the interpretation of graph-theoretic structures and measures can be highly dependent on the dynamics associated to the network. Properties that are quite meaningful for linear dynamics, such as random walk and network flow models, may be of limited relevance in the neuroscience setting. Theoretical and computational neuroscience are playing a vital role in understanding the relationship between network connectivity and the nonlinear dynamics associated to neural networks.

Evolutionary algorithm optimization of biological learning parameters in a biomimetic neuroprosthesis

Biomimetic simulation permits neuroscientists to better understand the complex neuronal dynamics of the brain. Embedding a biomimetic simulation in a closed-loop neuroprosthesis, which can read and write signals from the brain, will permit applications for amelioration of motor, psychiatric, and memory-related brain disorders. Biomimetic neuroprostheses require real-time adaptation to changes in the external environment, thus constituting an example of a dynamic data-driven application system. As model fidelity increases, so does the number of parameters and the complexity of finding appropriate parameter configurations. Instead of adapting synaptic weights via machine learning, we employed major biological learning methods: spike-timing dependent plasticity and reinforcement learning. We optimized the learning metaparameters using evolutionary algorithms, which were implemented in parallel and which used an island model approach to obtain sufficient speed. We employed these methods to train a cortical spiking model to utilize macaque brain activity, indicating a selected target, to drive a virtual musculoskeletal arm with realistic anatomical and biomechanical properties to reach to that target. The optimized system was able to reproduce macaque data from a comparable experimental motor task. These techniques can be used to efficiently tune the parameters of multiscale systems, linking realistic neuronal dynamics to behavior, and thus providing a useful tool for neuroscience and neuroprosthetics.

Multitarget Multiscale Simulation for Pharmacological Treatment of Dystonia in Motor Cortex

A large number of physiomic pathologies can produce hyperexcitability in cortex. Depending on severity, cortical hyperexcitability may manifest clinically as a hyperkinetic movement disorder or as epilpesy. We focus here on dystonia, a movement disorder that produces involuntary muscle contractions and involves pathology in multiple brain areas including basal ganglia, thalamus, cerebellum, and sensory and motor cortices. Most research in dystonia has focused on basal ganglia, while much pharmacological treatment is provided directly at muscles to prevent contraction. Motor cortex is another potential target for therapy that exhibits pathological dynamics in dystonia, including heightened activity and altered beta oscillations. We developed a multiscale model of primary motor cortex, ranging from molecular, up to cellular, and network levels, containing 1715 compartmental model neurons with multiple ion channels and intracellular molecular dynamics. We wired the model based on electrophysiological data obtained from mouse motor cortex circuit mapping experiments. We used the model to reproduce patterns of heightened activity seen in dystonia by applying independent random variations in parameters to identify pathological parameter sets. These models demonstrated degeneracy, meaning that there were many ways of obtaining the pathological syndrome. There was no single parameter alteration which would consistently distinguish pathological from physiological dynamics. At higher dimensions in parameter space, we were able to use support vector machines to distinguish the two patterns in different regions of space and thereby trace multitarget routes from dystonic to physiological dynamics. These results suggest the use of in silico models for discovery of multitarget drug cocktails.

Calcium regulation of HCN channels supports persistent activity in a multiscale model of neocortex

Neuronal persistent activity has been primarily assessed in terms of electrical mechanisms, without attention to the complex array of molecular events that also control cell excitability. We developed a multiscale neocortical model proceeding from the molecular to the network level to assess the contributions of calcium (Ca(2+)) regulation of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels in providing additional and complementary support of continuing activation in the network. The network contained 776 compartmental neurons arranged in the cortical layers, connected using synapses containing AMPA/NMDA/GABAA/GABAB receptors. Metabotropic glutamate receptors (mGluR) produced inositol triphosphate (IP3) which caused the release of Ca(2+) from endoplasmic reticulum (ER) stores, with reuptake by sarco/ER Ca(2+)-ATP-ase pumps (SERCA), and influence on HCN channels. Stimulus-induced depolarization led to Ca(2+) influx via NMDA and voltage-gated Ca(2+) channels (VGCCs). After a delay, mGluR activation led to ER Ca(2+) release via IP3 receptors. These factors increased HCN channel conductance and produced firing lasting for ∼1min. The model displayed inter-scale synergies among synaptic weights, excitation/inhibition balance, firing rates, membrane depolarization, Ca(2+) levels, regulation of HCN channels, and induction of persistent activity. The interaction between inhibition and Ca(2+) at the HCN channel nexus determined a limited range of inhibition strengths for which intracellular Ca(2+) could prepare population-specific persistent activity. Interactions between metabotropic and ionotropic inputs to the neuron demonstrated how multiple pathways could contribute in a complementary manner to persistent activity. Such redundancy and complementarity via multiple pathways is a critical feature of biological systems. Mediation of activation at different time scales, and through different pathways, would be expected to protect against disruption, in this case providing stability for persistent activity.

Restoring Behavior via Inverse Neurocontroller in a Lesioned Cortical Spiking Model Driving a Virtual Arm

Neural stimulation can be used as a tool to elicit natural sensations or behaviors by modulating neural activity. This can be potentially used to mitigate the damage of brain lesions or neural disorders. However, in order to obtain the optimal stimulation sequences, it is necessary to develop neural control methods, for example by constructing an inverse model of the target system. For real brains, this can be very challenging, and often unfeasible, as it requires repeatedly stimulating the neural system to obtain enough probing data, and depends on an unwarranted assumption of stationarity. By contrast, detailed brain simulations may provide an alternative testbed for understanding the interactions between ongoing neural activity and external stimulation. Unlike real brains, the artificial system can be probed extensively and precisely, and detailed output information is readily available. Here we employed a spiking network model of sensorimotor cortex trained to drive a realistic virtual musculoskeletal arm to reach a target. The network was then perturbed, in order to simulate a lesion, by either silencing neurons or removing synaptic connections. All lesions led to significant behvaioral impairments during the reaching task. The remaining cells were then systematically probed with a set of single and multiple-cell stimulations, and results were used to build an inverse model of the neural system. The inverse model was constructed using a kernel adaptive filtering method, and was used to predict the neural stimulation pattern required to recover the pre-lesion neural activity. Applying the derived neurostimulation to the lesioned network improved the reaching behavior performance. This work proposes a novel neurocontrol method, and provides theoretical groundwork on the use biomimetic brain models to develop and evaluate neurocontrollers that restore the function of damaged brain regions and the corresponding motor behaviors.

Towards real-time communication between in vivo neurophysiological data sources and simulator-based brain biomimetic models

Development of more sophisticated implantable brain-machine interface (BMI) will require both interpretation of the neurophysiological data being measured and subsequent determination of signals to be delivered back to the brain. Computational models are the heart of the machine of BMI and therefore an essential tool in both of these processes. One approach is to utilize brain biomimetic models (BMMs) to develop and instantiate these algorithms. These then must be connected as hybrid systems in order to interface the BMM with in vivo data acquisition devices and prosthetic devices. The combined system then provides a test bed for neuroprosthetic rehabilitative solutions and medical devices for the repair and enhancement of damaged brain. We propose here a computer network-based design for this purpose, detailing its internal modules and data flows. We describe a prototype implementation of the design, enabling interaction between the Plexon Multichannel Acquisition Processor (MAP) server, a commercial tool to collect signals from microelectrodes implanted in a live subject and a BMM, a NEURON-based model of sensorimotor cortex capable of controlling a virtual arm. The prototype implementation supports an online mode for real-time simulations, as well as an offline mode for data analysis and simulations without real-time constraints, and provides binning operations to discretize continuous input to the BMM and filtering operations for dealing with noise. Evaluation demonstrated that the implementation successfully delivered monkey spiking activity to the BMM through LAN environments, respecting real-time constraints.

Major remaining gaps in models of sensorimotor systems

Experimental descriptions of the anatomy and physiology of individual components of sensorimotor systems have revealed substantial complexity, making it difficult to intuit how complete systems might work. This has led to increasing efforts to develop and employ mathematical models to study the emergent properties of such systems. Conversely, the development of such models tends to reveal shortcomings in the experimental database upon which models must be constructed and validated. In both cases models are most useful when they point up discrepancies between what we think we know and possibilities that we may have overlooked. This overview considers those components of complete sensorimotor systems that currently appear to be potentially important but poorly understood. These are generally omitted completely from modeled systems or buried in implicit assumptions that underlie the design of the model.