1. Computational and Systems Biology
  2. Neuroscience

NetPyNE, a tool for data-driven multiscale modeling of brain circuits

  1. Salvador Dura-Bernal  Is a corresponding author
  2. Benjamin A Suter
  3. Padraig Gleeson
  4. Matteo Cantarelli
  5. Adrian Quintana
  6. Facundo Rodriguez
  7. David J Kedziora
  8. George L Chadderdon
  9. Cliff C Kerr
  10. Samuel A Neymotin
  11. Robert A McDougal
  12. Michael Hines
  13. Gordon M G Shepherd
  14. William W Lytton
  1. State University of New York Downstate Medical Center, United States
  2. Northwestern University, United States
  3. University College London, United Kingdom
  4. Metacell LLC, United States
  5. EyeSeeTea Ltd, United Kingdom
  6. University of Sydney, Australia
  7. Yale University, United States
https://doi.org/10.7554/eLife.44494
Share this article
Cite this article
  1. Salvador Dura-Bernal
  2. Benjamin A Suter
  3. Padraig Gleeson
  4. Matteo Cantarelli
  5. Adrian Quintana
  6. Facundo Rodriguez
  7. David J Kedziora
  8. George L Chadderdon
  9. Cliff C Kerr
  10. Samuel A Neymotin
  11. Robert A McDougal
  12. Michael Hines
  13. Gordon M G Shepherd
  14. William W Lytton
(2019)
NetPyNE, a tool for data-driven multiscale modeling of brain circuits
eLife 8:e44494.
https://doi.org/10.7554/eLife.44494
  2. Download BibTeX
  3. Download .RIS

Abstract

Biophysical modeling of neuronal networks helps to integrate and interpret rapidly growing and disparate experimental datasets at multiple scales. The NetPyNE tool (www.netpyne.org) provides both programmatic and graphical interfaces to develop data-driven multiscale network models in NEURON. NetPyNE clearly separates model parameters from implementation code. Users provide specifications at a high level via a standardized declarative language, e.g. connectivity rules, to create millions of cell-to-cell connections. NetPyNE then enables users to generate the NEURON network, run efficiently parallelized simulations, optimize and explore network parameters through automated batch runs, and use built-in functions for visualization and analysis - connectivity matrices, voltage traces, spike raster plots, local field potentials, and information theoretic measures. NetPyNE also facilitates model sharing by exporting and importing standardized formats (NeuroML and SONATA). NetPyNE is already being used to teach computational neuroscience students and by modelers to investigate brain regions and phenomena.

Data availability

All data and models used in this work are publicly available from the following GitHub and ModelDB links:- Fig 3: https://github.com/Neurosim-lab/netpyne/tree/paper/examples/paper/fig3- Fig 5: https://github.com/Neurosim-lab/netpyne/tree/paper/examples/paper/fig5- Fig 6: https://github.com/Neurosim-lab/netpyne/tree/paper/examples/paper/fig6- Fig 7: https://github.com/Neurosim-lab/netpyne/tree/paper/examples/paper/fig7- Fig 8: https://github.com/Neurosim-lab/netpyne/tree/paper/examples/paper/fig8- Fig 9A: https://github.com/ceciliaromaro/PD_in_NetPyNE- Fig 9B: https://github.com/rtekin/myKnoxRepo- Fig 9C: https://github.com/rodriguez-facundo/LASCON-project- Fig 9D: https://github.com/angietep/CA1-NetPyNE-modelTable 1:Dentate gyrus:- Original: https://modeldb.yale.edu/155568- NetPyNE: https://github.com/rodriguez-facundo/LASCON-projectCA1 microcircuits:- Original: https://modeldb.yale.edu/123815- NetPyNE: https://github.com/angietep/CA1-NetPyNE-modelEpilepsy in thalamocortex:- Original: https://modeldb.yale.edu/234233- NetPyNE: https://github.com/rodriguez-facundo/LASCON-projectEEG and MEG in cortex / HNN model:- Original: https://github.com/jonescompneurolab/hnn- NetPyNE: https://github.com/jonescompneurolab/hnn/tree/netpyneMotor cortex with RL:
- Original: https://modeldb.yale.edu/183014- NetPyNE: https://github.com/Neurosim-lab/netpyne/tree/development/examples/RL_armCortical microcircuits:- Original: https://github.com/OpenSourceBrain/PotjansDiesmann2014/tree/master/PyNEST- NetPyNE: https://github.com/ceciliaromaro/PD_in_NetPyNE

Article and author information

Author details

  1. Salvador Dura-Bernal

    Department of Physiology and Pharmacology, State University of New York Downstate Medical Center, Brooklyn, United States
    For correspondence
    salvadordura@gmail.com
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8561-5324

  2. Benjamin A Suter

    Department of Physiology, Northwestern University, Chicago, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9885-6936

  3. Padraig Gleeson

    Department of Neuroscience, Physiology and Pharmacology, University College London, London, United Kingdom
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5963-8576

  4. Matteo Cantarelli

    Metacell LLC, Boston, United States
    Competing interests
    Matteo Cantarelli, is affiliated with Metacell LLC. The author has no other competing interests to declare..

  5. Adrian Quintana

    EyeSeeTea Ltd, Cheltenham, United Kingdom
    Competing interests
    Adrian Quintana, is affiliated with EyeSeeTea Ltd. The author has no other competing interests to declare..

  6. Facundo Rodriguez

    Department of Physiology and Pharmacology, State University of New York Downstate Medical Center, Brooklyn, United States
    Competing interests
    No competing interests declared.

  7. David J Kedziora

    Complex Systems Group, School of Physics, University of Sydney, Sydney, Australia
    Competing interests
    No competing interests declared.

  8. George L Chadderdon

    Department of Physiology and Pharmacology, State University of New York Downstate Medical Center, Brooklyn, United States
    Competing interests
    No competing interests declared.

  9. Cliff C Kerr

    Complex Systems Group, School of Physics, University of Sydney, Sydney, Australia
    Competing interests
    No competing interests declared.

  10. Samuel A Neymotin

    Department of Physiology and Pharmacology, State University of New York Downstate Medical Center, Brooklyn, United States
    Competing interests
    No competing interests declared.

  11. Robert A McDougal

    Department of Neuroscience, Yale University, New Haven, United States
    Competing interests
    No competing interests declared.

  12. Michael Hines

    Department of Neuroscience, Yale University, New Haven, United States
    Competing interests
    No competing interests declared.

  13. Gordon M G Shepherd

    Department of Physiology, Northwestern University, Chicago, United States
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1455-8262

  14. William W Lytton

    Department of Physiology and Pharmacology, State University of New York Downstate Medical Center, Brooklyn, United States
    Competing interests
    No competing interests declared.

Funding

National Institute of Biomedical Imaging and Bioengineering (U01EB017695)

  • Salvador Dura-Bernal
  • Benjamin A Suter
  • Matteo Cantarelli
  • Adrian Quintana
  • Facundo Rodriguez
  • Samuel A Neymotin
  • Michael Hines
  • Gordon M G Shepherd
  • William W Lytton

New York State Department of Health (DOH01-C32250GG-3450000)

  • Salvador Dura-Bernal
  • Facundo Rodriguez
  • William W Lytton

Wellcome (101445)

  • Padraig Gleeson

National Institute of Biomedical Imaging and Bioengineering (2R01DC012947-06A1)

  • Samuel A Neymotin

National Institute of Biomedical Imaging and Bioengineering (R01EB022903)

  • Salvador Dura-Bernal
  • Michael Hines
  • William W Lytton

National Institute of Biomedical Imaging and Bioengineering (R01MH086638)

  • Robert A McDougal
  • Michael Hines
  • William W Lytton

Wellcome (212941)

  • Padraig Gleeson

National Institute of Biomedical Imaging and Bioengineering (3R01EB022889)

  • Salvador Dura-Bernal
  • Matteo Cantarelli
  • Adrian Quintana
  • Facundo Rodriguez
  • Samuel A Neymotin
  • Michael Hines

Australian Research Council (DE140101375)

  • David J Kedziora
  • Cliff C Kerr

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. Upinder Singh Bhalla, Tata Institute of Fundamental Research, India

Publication history

  1. Received: December 19, 2018
  2. Accepted: April 25, 2019
  3. Accepted Manuscript published: April 26, 2019 (version 1)
  4. Version of Record published: May 24, 2019 (version 2)
  5. Version of Record updated: May 31, 2019 (version 3)

Copyright

© 2019, Dura-Bernal et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 5,634
    Page views
  • 646
    Downloads
  • 54
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, Scopus, PubMed Central.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Salvador Dura-Bernal
  2. Benjamin A Suter
  3. Padraig Gleeson
  4. Matteo Cantarelli
  5. Adrian Quintana
  6. Facundo Rodriguez
  7. David J Kedziora
  8. George L Chadderdon
  9. Cliff C Kerr
  10. Samuel A Neymotin
  11. Robert A McDougal
  12. Michael Hines
  13. Gordon M G Shepherd
  14. William W Lytton
(2019)
NetPyNE, a tool for data-driven multiscale modeling of brain circuits
eLife 8:e44494.
https://doi.org/10.7554/eLife.44494

Categories and tags

Further reading

    1. Cancer Biology
    2. Computational and Systems Biology

    Predictive nonlinear modeling of malignant myelopoiesis and tyrosine kinase inhibitor therapy

    Jonathan Rodriguez, Abdon Iniguez ... Richard A Van Etten
    Research Article Updated

    Chronic myeloid leukemia (CML) is a blood cancer characterized by dysregulated production of maturing myeloid cells driven by the product of the Philadelphia chromosome, the BCR-ABL1 tyrosine kinase. Tyrosine kinase inhibitors (TKIs) have proved effective in treating CML, but there is still a cohort of patients who do not respond to TKI therapy even in the absence of mutations in the BCR-ABL1 kinase domain that mediate drug resistance. To discover novel strategies to improve TKI therapy in CML, we developed a nonlinear mathematical model of CML hematopoiesis that incorporates feedback control and lineage branching. Cell–cell interactions were constrained using an automated model selection method together with previous observations and new in vivo data from a chimeric BCR-ABL1 transgenic mouse model of CML. The resulting quantitative model captures the dynamics of normal and CML cells at various stages of the disease and exhibits variable responses to TKI treatment, consistent with those of CML patients. The model predicts that an increase in the proportion of CML stem cells in the bone marrow would decrease the tendency of the disease to respond to TKI therapy, in concordance with clinical data and confirmed experimentally in mice. The model further suggests that, under our assumed similarities between normal and leukemic cells, a key predictor of refractory response to TKI treatment is an increased maximum probability of self-renewal of normal hematopoietic stem cells. We use these insights to develop a clinical prognostic criterion to predict the efficacy of TKI treatment and design strategies to improve treatment response. The model predicts that stimulating the differentiation of leukemic stem cells while applying TKI therapy can significantly improve treatment outcomes.

    1. Computational and Systems Biology

    Sex, strain, and lateral differences in brain cytoarchitecture across a large mouse population

    David Elkind, Hannah Hochgerner ... Amit Zeisel
    Research Article Updated

    The mouse brain is by far the most intensively studied among mammalian brains, yet basic measures of its cytoarchitecture remain obscure. For example, quantifying cell numbers, and the interplay of sex, strain, and individual variability in cell density and volume is out of reach for many regions. The Allen Mouse Brain Connectivity project produces high-resolution full brain images of hundreds of brains. Although these were created for a different purpose, they reveal details of neuroanatomy and cytoarchitecture. Here, we used this population to systematically characterize cell density and volume for each anatomical unit in the mouse brain. We developed a DNN-based segmentation pipeline that uses the autofluorescence intensities of images to segment cell nuclei even within the densest regions, such as the dentate gyrus. We applied our pipeline to 507 brains of males and females from C57BL/6J and FVB.CD1 strains. Globally, we found that increased overall brain volume does not result in uniform expansion across all regions. Moreover, region-specific density changes are often negatively correlated with the volume of the region; therefore, cell count does not scale linearly with volume. Many regions, including layer 2/3 across several cortical areas, showed distinct lateral bias. We identified strain-specific or sex-specific differences. For example, males tended to have more cells in extended amygdala and hypothalamic regions (MEA, BST, BLA, BMA, and LPO, AHN) while females had more cells in the orbital cortex (ORB). Yet, inter-individual variability was always greater than the effect size of a single qualifier. We provide the results of this analysis as an accessible resource for the community.

    1. Computational and Systems Biology
    2. Immunology and Inflammation

    Statistical inference reveals the role of length, GC content, and local sequence in V(D)J nucleotide trimming

    Magdalena L Russell, Noah Simon ... Frederick A Matsen IV
    Research Article

    To appropriately defend against a wide array of pathogens, humans somatically generate highly diverse repertoires of B cell and T cell receptors (BCRs and TCRs) through a random process called V(D)J recombination. Receptor diversity is achieved during this process through both the combinatorial assembly of V(D)J-genes and the junctional deletion and insertion of nucleotides. While the Artemis protein is often regarded as the main nuclease involved in V(D)J recombination, the exact mechanism of nucleotide trimming is not understood. Using a previously published TCRβ repertoire sequencing data set, we have designed a flexible probabilistic model of nucleotide trimming that allows us to explore various mechanistically interpretable sequence-level features. We show that local sequence context, length, and GC nucleotide content in both directions of the wider sequence, together, can most accurately predict the trimming probabilities of a given V-gene sequence. Because GC nucleotide content is predictive of sequence-breathing, this model provides quantitative statistical evidence regarding the extent to which double-stranded DNA may need to be able to breathe for trimming to occur. We also see evidence of a sequence motif that appears to get preferentially trimmed, independent of GC-content-related effects. Further, we find that the inferred coefficients from this model provide accurate prediction for V- and J-gene sequences from other adaptive immune receptor loci. These results refine our understanding of how the Artemis nuclease may function to trim nucleotides during V(D)J recombination and provide another step toward understanding how V(D)J recombination generates diverse receptors and supports a powerful, unique immune response in healthy humans.