Introduction

This package contains code used in
Effects of Ih and TASK-like shunting current on dendritic impedance in layer 5 pyramidal-tract type neurons (Kelley et al. 2020).

This includes generalizable tools for computing neuronal impedance, standalone versions of pyramidal tract neuron models previously implemented in other studies, and scripts for batch processing and analysis.

Basic Use

Make sure you have installed NEURON and NETPYNE.

Clone this repository using 'git'

git clone https://github.com/suny-downstate-medical-center/L5PYR_Resonance.git

Compile .mod files for the cell model(s) you want to use. To run init.py, you need to compile the .mod files from Hay et al. 2011 in './models/Hay'.

cd L5PYR_Resonance/models/Hay/
nrnivmodl mod
cd ../../

To run an example simulation and compute impedance profiles for a single dendritic compartment:

python3 -i init.py

To simulate chirp stimulation in the model neuron from Dura-Bernal et al. 2019, run the following:

cd L5PYR_Resonance/models/DuraBernal
nrnivmodl mod
python3 -i runChirp.py

Code

init.py

Example simulation of chirp current stimulation of the apical trunk and computation of transfer impedance between stimulated dendrite and the soma. The chirp stimulus provides sinusoidal current stimultion with instantaneous frequency linearly increasing from 0.5-10Hz in 10s. Transfer impedance is computed from the resulting somatic membrane potential waveforms.

getCells.py

Contains functions for loading stand-alone versions of pyramidal-tract type neurons used in the manuscript. There are also functions for automatically downloading and instantiating models from the Allen Brain Atlas.

chirpUtils.py

Contains functions for simulating chirp current stimulation and computing neuronal impedance.

./models/

Contains folders with each of the cell models studied in the manuscript. Each of those folders contains .mod files that must be compiled prior to runtime. Also includes simplified model presented in the paper that was built with NEURON's LinearCircuitBuilder: 'cellwithL.ses'. This model can be used interactively by running:

from neuron import h, gui
h.load_file('models/cellwithL.ses')

Note that circuit elements will not be recognized by NEURON unless the 'Simulate' button is selected on the LinearCircuit window.

./simcode/

Contains scripts used to run the batch simulations and analyses used in the study. Many of these are machine specific and weren't intended for general use.