NetPyNE, a tool for data-driven multiscale modeling of brain circuits
Figures
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig1-v2.tif/full/617,/0/default.jpg)
Overview of NetPyNE components and workflow.
Users start by specifying the network parameters and simulation configuration using a high-level JSON-like format. Existing NEURON and NeuroML models can be imported. Next, a NEURON network model is instantiated based on these specifications. This model can be simulated in parallel using NEURON as the underlying simulation engine. Simulation results are gathered in the master node. Finally, the user can analyze the network and simulation results using a variety of plots; save to multiple formats or export to NeuroML. The Batch Simulation module enables automating this process to run multiple simulations on HPCs and explore a range of parameter values.
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig2-v2.tif/full/617,/0/default.jpg)
High-level specification of network parameters.
(A) Programmatic parameter specification using standardized declarative JSON-like format. i,ii: specification of two populations iii,iv: cell parameters; v: reaction-diffusion parameters; vi,vii,viii: synapse parameters and connectivity rules. (B) GUI-based parameter specification, showing the definition of populations equivalent to those in panel A. (C) Schematic of network model resulting from the specifications in A.
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig3-v2.tif/full/617,/0/default.jpg)
Specification of dendritic distribution of synapses.
(A) Optogenetic data provides synapse density across the 2D grid shown at left (Suter and Shepherd, 2015). (B) Data are imported directly into NetPyNE which automatically calculates synapse location in simplified or full multicompartmental representations of a pyramidal cell. (C) Corresponding synaptic density plot generated by NetPyNE.
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig4-v2.tif/full/617,/0/default.jpg)
Instantiated network hierarchical data model.
The instantiated network is represented using a standardized hierarchically organized Python structure generated from NetPyNE ’s high-level specifications. This data structure provides direct access to all elements, state variables and parameters to be simulated. Defined NEURON simulator objects (represented as boxes with red borders) are included within the Python data structure.
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig5-v2.tif/full/617,/0/default.jpg)
NetPyNE visualization and analysis plots for a simple three-layer network example.
(A) Connectivity matrix, (B) stacked bar graph, (C) 2D graph representation of cells and connections, (D) voltage traces of three cells, (E) spike raster plot, (F) population firing rate statistics (boxplot).
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig6-v2.tif/full/617,/0/default.jpg)
LFP recording and analysis.
(A) LFP signals (left) from 10 extracellular recording electrodes located around a morphologically detailed cell (right) producing a single action potential (top-right). (B) LFP signals, PSDs and spectrograms (left and center) from four extracellular recording electrodes located at different depths of a network of 120 five-compartment neurons (right) producing oscillatory activity (top-left).
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig7-v2.tif/full/617,/0/default.jpg)
NetPyNE graphical user interface (GUI) showing a multiscale model.
Background shows 3D representation of example network with 6 populations of multi-channel, multi-compartment neurons; results panels from left to right: singe-neuron traces (voltage, intracellular and extracellular calcium concentration, and potassium current); spike raster plot; extracellular potassium concentration; LFP signals recorded from three electrodes; and 3D location of the LFP electrodes within network.
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig8-v2.tif/full/617,/0/default.jpg)
Model of M1 microcircuits developed using NetPyNE (scaled down version).
NetPyNE GUI showing 3D representation of M1 network (background), spike raster plot and population firing rate statistics (top left), voltage traces (bottom left) and firing rate power spectral density (top right).
![](https://iiif.elifesciences.org/lax/44494%2Felife-44494-fig9-v2.tif/full/617,/0/default.jpg)
Published models converted to NetPyNE.
All figures were generated using the NetPyNE version of the models. (A) Spike raster plot and boxplot statistics of the Potjans and Diesmann thalamocortical network originally implemented in NEST (Potjans and Diesmann, 2014; Romaro et al., 2018). (B) Spike raster plot and voltage traces of a thalamocortical network exhibiting epileptic activity originally implemented in NEURON/hoc (Knox et al., 2018). (C) 3D representation of the cell types and network topology, and spike raster plots of a dentate gyrus model originally implemented in NEURON/hoc (Rodriguez, 2018; Tejada et al., 2014). (D) Connectivity rules (top) and voltage traces of 2 cell types (bottom) in a hippocampal CA1 model originally implemented in NEURON/hoc (Cutsuridis et al., 2010; Tepper et al., 2018).
Tables
Number of lines of code in the original models and the NetPyNE reimplementations.
https://doi.org/10.7554/eLife.44494.012Model description (reference) | Original language | Original num lines | NetPyNE num lines |
---|---|---|---|
Dentate gyrus (Tejada et al., 2014) | NEURON/hoc | 1029 | 261 |
CA1 microcircuits (Cutsuridis et al., 2010) | NEURON/hoc | 642 | 306 |
Epilepsy in thalamo (Knox et al., 2018) | NEURON/hoc | 556 | 201 |
EEG and MEG in cortex/HNN model (Jones et al., 2009) | NEURON/Python | 2288 | 924 |
Motor cortex with RL (Dura-Bernal et al., 2017) | NEURON/Python | 1171 | 362 |
Cortical microcircuits (Potjans and Diesmann, 2014) | PyNEST | 689 | 198 |
Additional files
-
Transparent reporting form
- https://doi.org/10.7554/eLife.44494.013