Abstract:The primary somatosensory cortex (S1) of mammals is critically important in the perception of touch and related sensorimotor behaviors. In 2015, the Blue Brain Project (BBP) developed a groundbreaking rat S1 microcircuit simulation with over 31,000 neurons with 207 morpho-electrical neuron types, and 37 million synapses, incorporating anatomical and physiological information from a wide range of experimental studies. We have implemented this highly detailed and complex S1 model in NetPyNE, using the data available in the Neocortical Microcircuit Collaboration Portal. NetPyNE provides a Python high-level interface to NEURON and allows defining complicated multiscale models using an intuitive declarative standardized language. It also facilitates running parallel simulations, automates the optimization and exploration of parameters using supercomputers, and provides a wide range of built-in analysis functions. This will make the S1 model more accessible and simpler to scale, modify and extend in order to explore research questions or interconnect to other existing models. Despite some implementation differences, the NetPyNE model preserved the original cell morphologies, electrophysiological responses and spatial distribution for all 207 cell types; and the connectivity properties of all 1941 pathways, including synaptic dynamics and short-term plasticity (STP). The NetPyNE S1 simulations produced reasonable physiological firing rates and activity patterns across all populations. When STP was included, the network generated a 1 Hz oscillation comparable to the original model in vitro -like state. By then reducing the extracellular calcium concentration, the model reproduced the original S1 in vivo -like states with asynchronous activity. These results validate the original study using a new modeling tool. Simulated local field potentials (LFPs) exhibited realistic oscillatory patterns and features, including distance- and frequency-dependent attenuation. The model was extended by adding thalamic circuits, including 6 distinct thalamic populations with intrathalamic, thalamocortical (TC) and corticothalamic connectivity derived from experimental data. The thalamic model reproduced single known cell and circuit-level dynamics, including burst and tonic firing modes and oscillatory patterns, providing a more realistic input to cortex and enabling study of TC interactions. Overall, our work provides a widely accessible, data-driven and biophysically-detailed model of the somatosensory TC circuits that can be employed as a community tool for researchers to study neural dynamics, function and disease.

Modeling and Simulation of Neocortical Micro- and Mesocircuitry. Part II: Physiology and Experimentation

Modeling and Simulation of Rat Non-Barrel Somatosensory Cortex. Part I: Modeling Anatomy

Large-scale biophysically detailed model of somatosensory thalamocortical circuits in NetPyNE

Large-scale model of mammalian thalamocortical systems

Large-Scale Mechanistic Models of Brain Circuits with Biophysically and Morphologically Detailed Neurons

A Layered Microcircuit Model of Somatosensory Cortex with Three Interneuron Types and Cell-Type-Specific Short-Term Plasticity

A biophysically-detailed model of inter-areal interactions in cortical sensory processing

Biomimetic model of corticostriatal micro-assemblies discovers new neural code

NeoCOMM: A Neocortical Neuroinspired Computational Model for the Reconstruction and Simulation of Epileptiform Events

Multiscale modeling of neuronal dynamics in hippocampus CA1

Simulation of transcranial magnetic stimulation in head model with morphologically-realistic cortical neurons

Cellular Connectomes As Arbiters of Local Circuit Models in the Cerebral Cortex.

Cortical multi-area model with joint excitatory-inhibitory clusters accounts for spiking statistics, inter-area propagation, and variability dynamics

Assemblies, synapse clustering and network topology interact with plasticity to explain structure-function relationships of the cortical connectome

The dynamic brain: from spiking neurons to neural masses and cortical fields

Multi-scale spiking network model of human cerebral cortex

Cortical, striatal, and thalamic populations self-organize into a functionally connected circuit with long-term memory properties

Multiscale modelling of neuronal dynamics in hippocampus CA1

Stimulation mapping and whole-brain modeling reveal gradients of excitability and recurrence in cortical networks

NeoCoMM: Neocortical Computational Microscale Model