Spiny neuron model with dopamine-induced bistability (Gruber et al 2003)

 Download zip file 
Help downloading and running models
Accession:39949
These files implement a model of dopaminergic modulation of voltage-gated currents (called kir2 and caL in the original paper). See spinycell.html for details of usage and implementation. For questions about this implementation, contact Ted Carnevale (ted.carnevale@yale.edu)
Reference:
1 . Gruber AJ, Solla SA, Surmeier DJ, Houk JC (2003) Modulation of striatal single units by expected reward: a spiny neuron model displaying dopamine-induced bistability. J Neurophysiol 90:1095-114 [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): Neostriatum medium spiny direct pathway GABA cell;
Channel(s): I K,leak; I h; I A, slow;
Gap Junctions:
Receptor(s): D1; Dopaminergic Receptor;
Gene(s): Cav1.3 CACNA1D; D1 DRD1A;
Transmitter(s): Dopamine;
Simulation Environment: NEURON;
Model Concept(s): Intrinsic plasticity;
Implementer(s): Carnevale, Ted [Ted.Carnevale at Yale.edu];
Search NeuronDB for information about:  Neostriatum medium spiny direct pathway GABA cell; D1; Dopaminergic Receptor; I K,leak; I h; I A, slow; Dopamine;
objectvar save_window_, rvp_
objectvar scene_vector_[9]
objectvar ocbox_, ocbox_list_, scene_, scene_list_
{ocbox_list_ = new List()  scene_list_ = new List()}

//Begin PointProcessManager
{
load_file("pointman.hoc")
}
{
C_Cell[0].soma ocbox_ = new PointProcessManager(0)
}
{object_push(ocbox_)}
{
mt.select("GSpike") i = mt.selected()
ms[i] = new MechanismStandard("GSpike")
ms[i].set("vt", -58, 0)
ms[i].set("vh", -55, 0)
ms[i].set("vc", 2.5, 0)
ms[i].set("minisi", 20, 0)
mt.select("GSpike") i = mt.selected() maction(i)
hoc_ac_ = 0
sec.sec move() d1.flip_to(0)
}
{object_pop() doNotify()}
{
ocbox_ = ocbox_.v1
ocbox_.map("PointProcessManager", 504, 414, 208.32, 326.4)
}
objref ocbox_
//End PointProcessManager

{
xpanel("RunControl", 0)
v_init = -90
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 = 0
xvalue("t","t", 2 )
tstop = 2000
xvalue("Tstop","tstop", 1,"tstop_changed()", 0, 1 )
dt = 0.025
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 = 0
xvalue("Real Time","realtime", 0,"", 0, 1 )
xpanel(439,22)
}
{
save_window_ = new Graph(0)
save_window_.size(0,2000,-90,-40)
scene_vector_[6] = save_window_
{save_window_.view(0, -90, 2000, 50, 733, 161, 300.48, 200.32)}
graphList[0].append(save_window_)
save_window_.save_name("graphList[0].")
save_window_.addvar("C_Cell[0].soma.v( 0.5 )", 1, 1, 0.454952, 0.895208, 2)
}
{
save_window_ = new Graph(0)
save_window_.size(0,2000,-0.0159167,1.01491)
scene_vector_[7] = save_window_
{save_window_.view(0, -0.0159167, 2000, 1.03083, 733, 426, 300.48, 200.32)}
graphList[2].append(save_window_)
save_window_.save_name("graphList[2].")
save_window_.addvar("GSpike[0].spike", 1, 1, 0.538019, 0.933546, 2)
}
{
save_window_ = new Graph(0)
save_window_.size(0,2000,1,1.4)
scene_vector_[8] = save_window_
{save_window_.view(0, 1, 2000, 0.4, 733, 691, 300.48, 200.32)}
graphList[2].append(save_window_)
save_window_.save_name("graphList[2].")
save_window_.addvar("DAsyn[0].msg", 2, 3, 0.595527, 0.9, 2)
}
objectvar scene_vector_[1]
{doNotify()}

Loading data, please wait...