Cortical pyramidal neuron, phase response curve (Stiefel et al 2009)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:144372
Three models of increasing complexity all showing a switch from type II (biphasic) to type I (monophasic) phase response curves with a cholinergic down-modulation of K+ conductances.
Reference:
1 . Stiefel KM, Gutkin BS, Sejnowski TJ (2009) The effects of cholinergic neuromodulation on neuronal phase-response curves of modeled cortical neurons. J Comput Neurosci 26:289-301 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type:
Brain Region(s)/Organism:
Cell Type(s): Neocortex L2/3 pyramidal GLU cell;
Channel(s): I Na,p; I Na,t; I M;
Gap Junctions:
Receptor(s): Muscarinic;
Gene(s):
Transmitter(s): Acetylcholine;
Simulation Environment: NEURON;
Model Concept(s): Action Potentials;
Implementer(s): Stiefel, Klaus [stiefel at salk.edu];
Search NeuronDB for information about:  Neocortex L2/3 pyramidal GLU cell; Muscarinic; I Na,p; I Na,t; I M; Acetylcholine;
/
StiefelEtAl2009
README.txt
ca.mod *
cacum.mod
cad.mod *
H.mod
iahp2.mod *
il.mod *
im.mod *
KA.mod
kca.mod *
Kdr.mod
km.mod *
Ks.mod
kv.mod *
Na.mod *
NaP.mod
cell.ses
displayshape.hoc
fig4A.hoc
fig4A_new.hoc
fig5A.hoc
fig5B.hoc
fig5C.hoc
gui.hoc
j8.hoc *
ksprc.ses
makeIF.hoc
multi.hoc
PRC.hoc
PRCsweep.hoc
PY-golomb_original.hoc
PY-golomb_plus.hoc
PY-golomb_simple.hoc
PyMainen.hoc
single.hoc
single_plus.hoc
single1.ses
surface.hoc
synproxy_cch.hoc
synproxy_sweeps.hoc
                            
load_file("nrngui.hoc")
objectvar save_window_, rvp_
objectvar scene_vector_[5]
objectvar ocbox_, ocbox_list_, scene_, scene_list_
{ocbox_list_ = new List()  scene_list_ = new List()}
{pwman_place(0,0,0)}
{
xpanel("NEURON Main Menu", 0)
xmenu("File", 0)
xbutton("load session","load(\"*.ses\", 1)")
xbutton("load hoc","load(\"*.hoc\", 0)")
xbutton("load dll","load(\"*.dll\", 2)")
xbutton("save session","save_ses_nrnmainmenu()")
xbutton("working dir","change_working_dir()")
xmenu("recent dir", 0)
xbutton("Print","print_session(0)")
xbutton("Quit","quit()")
xmenu()
xmenu("Edit", 1)
xmenu()
xmenu("Build", 1)
xbutton("single compartment","load_file(\"single.hoc\") makesinglecompartment()")
xbutton("Cell Builder","load_file(\"celbild.hoc\") makecellbuilder()")
xmenu("NetWork Cell", 0)
xbutton("From Cell Builder","load_file(\"netbild.hoc\") makenetreadycell()")
xbutton("Artificial Cell","load_file(\"netbild.hoc\") makeartcelltype()")
xmenu()
xbutton("NetWork Builder","load_file(\"netbild.hoc\") makenetbuilder()")
xbutton("Linear Circuit","load_file(\"lincir.hoc\") makeCircuit()")
xbutton("Channel Builder","load_file(\"chanbild.hoc\", \"ChannelBuild\") chanbild()")
xmenu()
xmenu("Tools", 1)
xbutton("RunControl","nrncontrolmenu()")
xbutton("RunButton","xpanel(\"Run\") xcheckbox(\"Init & Run\", &running_,\"runbutton()\") xpanel()")
xbutton("VariableStepControl","numericalmethodpanel.map()")
xmenu("Point Processes", 0)
xmenu("Managers", 0)
xbutton("Point Manager","load_file(\"pointman.hoc\") makeppm()")
xbutton("Point Group","load_file(\"pointgrp.hoc\") makeppgm()")
xbutton("Electrode","load_file(\"electrod.hoc\") makeelectrode()")
xmenu()
xmenu("Viewers", 0)
xmenu()
xmenu("Distributed Mechanisms", 0)
xmenu("Managers", 0)
xbutton("Inserter","load_file(\"inserter.hoc\", \"Inserter\") makeinserter()")
xbutton("Homogeneous Spec","load_file(\"showmech.hoc\") makeshowmechanism()")
xmenu()
xmenu("Viewers", 0)
xbutton("celsius","celsius_panel()")
xbutton("globalRa","globalra_panel()")
xmenu()
xmenu("Fitting", 0)
xbutton("Parameterized Function","load_file(\"funfit.hoc\") makefitter()")
xbutton("Run Fitter","load_file(\"runfit.hoc\") makerunfitter()")
xbutton("Multiple Run Fitter","load_file(\"mulfit.hoc\") makemulrunfitter()")
xmenu()
xmenu("Impedance", 0)
xbutton("Frequency","load_file(\"impratio.hoc\") makeImpRatio()")
xbutton("Path","load_file(\"impedanx.hoc\") makeImpx()")
xbutton("log(A) vs x","load_file(\"logax.hoc\") makelogax()")
xbutton("Shape","load_file(\"attshape.hoc\") makeImpShape()")
xmenu()
xbutton("Model View","load_file(\"mview.hoc\", \"ModelView\") execute(\"mview()\")")
xmenu("Miscellaneous", 0)
xmenu()
xmenu("Graph", 1)
xbutton("Voltage axis","newPlotV()")
xbutton("Current axis","newPlotI()")
xbutton("State axis","newPlotS()")
xbutton("Shape plot","newshapeplot()")
xbutton("Vector movie","newvectorplot()")
xbutton("Phase Plane","newphaseplane()")
xbutton("Grapher","load_file(\"grapher.hoc\") makegrapher(1)")
xmenu()
xmenu("Vector", 1)
xbutton("Save to File","load_file(\"stdlib.hoc\",\"clipboard_save\") clipboard_save()")
xbutton("Retrieve from File","load_file(\"stdlib.hoc\",\"clipboard_retrieve\") clipboard_retrieve()")
xbutton("Gather Values","load_file(\"gatherv.hoc\") makeGatherVec()")
xbutton("Play","load_file(\"vplay.hoc\") makeVectorPlay()")
xbutton("Display","load_file(\"vecwrap.hoc\") makeVecWrap()")
xbutton("Draw","load_file(\"mkcurve.hoc\") makeMakeCurve()")
xmenu()
xmenu("Window", 1)
xpanel(0,0)
}
{
xpanel("RunControl", 0)
v_init = -70
xvalue("Init","v_init", 1,"stdinit()", 1, 1 )
xbutton("Init & Run","run()")
xbutton("Stop","stoprun=1")
runStopAt = 5
xvalue("Continue til","runStopAt", 1,"{continuerun(runStopAt) stoprun=1}", 1, 1 )
runStopIn = 1
xvalue("Continue for","runStopIn", 1,"{continuerun(t + runStopIn) stoprun=1}", 1, 1 )
xbutton("Single Step","steprun()")
t = 698.929
xvalue("t","t", 2 )
tstop = 1000
xvalue("Tstop","tstop", 1,"tstop_changed()", 0, 1 )
dt = 0.00218452
xvalue("dt","dt", 1,"setdt()", 0, 1 )
steps_per_ms = 40
xvalue("Points plotted/ms","steps_per_ms", 1,"setdt()", 0, 1 )
xcheckbox("Quiet",&stdrun_quiet,"")
realtime = 53
xvalue("Real Time","realtime", 0,"", 0, 1 )
xpanel(226,225)
}
{
save_window_ = new Graph(0)
save_window_.size(0,1000,-80,40)
scene_vector_[2] = save_window_
{save_window_.view(0, -80, 1000, 120, 900, 213, 299.52, 199.38)}
graphList[0].append(save_window_)
save_window_.save_name("graphList[0].")
save_window_.addexpr("v(.5)", 1, 1, 0.8, 0.9, 2)
}
{
xpanel("iseg(0 - 1) (Parameters)", 0)
xlabel("iseg(0 - 1) (Parameters)")
xlabel("nseg = 5")
iseg.L = 15
xvalue("L","iseg.L", 1,"define_shape()", 0, 0 )
iseg.cm = 0.75
xvalue("cm","iseg.cm", 1,"", 0, 0 )
iseg.diam = 0.99279
xvalue("diam","iseg.diam", 1,"", 0, 0 )
iseg.g_pas = 3.3333e-05
xvalue("g_pas","iseg.g_pas", 1,"", 0, 0 )
iseg.e_pas = -70
xvalue("e_pas","iseg.e_pas", 1,"", 0, 0 )
iseg.gbar_na = 30000
xvalue("gbar_na","iseg.gbar_na", 1,"", 0, 0 )
iseg.ena = 60
xvalue("ena","iseg.ena", 1,"", 0, 0 )
iseg.gbar_kv = 2000
xvalue("gbar_kv","iseg.gbar_kv", 1,"", 0, 0 )
iseg.gbar_km = 10
xvalue("gbar_km","iseg.gbar_km", 1,"", 0, 0 )
iseg.ek = -90
xvalue("ek","iseg.ek", 1,"", 0, 0 )
xpanel(620,684)
}

//Begin VariableTimeStep
{
ocbox_ = NumericalMethodPanel[0]
}
{object_push(ocbox_)}
{
atol_ = 0.001  CVode[0].atol(atol_)
restore(301, 1)
}
{object_pop()}
{
ocbox_.map("VariableTimeStep", 565, -22, 234, 121.68)
}
objref ocbox_
//End VariableTimeStep

objectvar scene_vector_[1]
{doNotify()}