Dendritic Impedance in Neocortical L5 PT neurons (Kelley et al. 2021)

 Download zip file 
Help downloading and running models
We simulated chirp current stimulation in the apical dendrites of 5 biophysically-detailed multi-compartment models of neocortical pyramidal tract neurons and found that a combination of HCN channels and TASK-like channels produced the best fit to experimental measurements of dendritic impedance. We then explored how HCN and TASK-like channels can shape the dendritic impedance as well as the voltage response to synaptic currents.
1 . Kelley C, Dura-Bernal S, Neymotin SA, Antic SD, Carnevale NT, Migliore M, Lytton WW (2021) Effects of Ih and TASK-like shunting current on dendritic impedance in layer 5 pyramidal-tract neurons. J Neurophysiology 125:1501-1516 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Neocortex L5/6 pyramidal GLU cell; Neocortex M1 L5B pyramidal pyramidal tract GLU cell;
Channel(s): I h; TASK channel;
Gap Junctions:
Simulation Environment: NEURON; Python; NetPyNE;
Model Concept(s): Impedance;
Implementer(s): Kelley, Craig;
Search NeuronDB for information about:  Neocortex L5/6 pyramidal GLU cell; Neocortex M1 L5B pyramidal pyramidal tract GLU cell; I h; TASK channel;
import sys
# sys.path.insert(0,'../my_netpyne')
import os
os.environ["OPENBLAS_NUM_THREADS"] = sys.argv[-1]
N = int(sys.argv[-1])
from os import listdir, mkdir
import numpy as np
import multiprocessing
from math import nan
from import savemat 

def chirpForMulti(invar):
    sec_num, loc, filename = invar
    from getCells import M1Cell
    s = M1Cell()
    seg =[0].secs[sec_num]['hObj'](loc)
    soma_seg =[0].secs['soma']['hObj'](0.5)

    from chirpUtils import applyChirp, getChirp
    amp = 0.025
    f0, f1, t0, Fs, delay = 0.5, 10, 10, 1000, 5 # for looking at bimodal leading phase response in Hay cell
    I, t = getChirp(f0, f1, t0, amp, Fs, delay)
    print('running chirp on ' + str(seg))
    applyChirp(I, t, seg, soma_seg, t0, delay, Fs, f1, out_file_name=filename)
    print(str(seg) + ': done')

#populate data tuple
data = []

from getCells import M1Cell
s = M1Cell()
for sec_key in[0].secLists['apical_maintrunk']:
    nseg =[0].secs[sec_key]['hObj'].nseg 
    if nseg == 1:
        data.append([sec_key, 0.5, '/u/craig/L5PYR_Resonance/M1_PTcell/trunk_data/'+sec_key+'_0.5'])
        for loc in np.linspace(1/(nseg+1), nseg/(nseg+1), nseg):
            data.append([sec_key, loc, '/u/craig/L5PYR_Resonance/M1_PTcell/trunk_data/'+sec_key+'_'+str(loc)])

data = tuple(data)

def mp_handler():
    p = multiprocessing.Pool(N), data)

if __name__ == '__main__':

Loading data, please wait...