Distinct current modules shape cellular dynamics in model neurons (Alturki et al 2016)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:223649
" ... We hypothesized that currents are grouped into distinct modules that shape specific neuronal characteristics or signatures, such as resting potential, sub-threshold oscillations, and spiking waveforms, for several classes of neurons. For such a grouping to occur, the currents within one module should have minimal functional interference with currents belonging to other modules. This condition is satisfied if the gating functions of currents in the same module are grouped together on the voltage axis; in contrast, such functions are segregated along the voltage axis for currents belonging to different modules. We tested this hypothesis using four published example case models and found it to be valid for these classes of neurons. ..."
Reference:
1 . Alturki A, Feng F, Nair A, Guntu V, Nair SS (2016) Distinct current modules shape cellular dynamics in model neurons. Neuroscience 334:309-331 [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: Hippocampus; Amygdala;
Cell Type(s): Abstract single compartment conductance based cell;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Simplified Models; Activity Patterns; Oscillations; Methods; Olfaction;
Implementer(s):
/
AlturkiEtAl2016
4_LA
Segregated
cadyn.mod *
capool.mod *
currentclamp.mod *
h.mod *
im.mod *
kaprox.mod *
kdrca1.mod *
kdrca1DA.mod *
leak.mod *
leakDA.mod *
na3.mod *
na3DA.mod *
nap.mod *
nat.mod *
sahp.mod *
sahpNE.mod *
graphics_lib.hoc *
main.hoc
main_HTO.hoc
main_LTO.hoc
onecompartment_template_with_osc.hoc
                            
{load_file("nrngui.hoc")}
{load_file("onecompartment_template_with_osc.hoc")}
{load_file("graphics_lib.hoc")}

v_init = -70//-70//-69
stim_amp = 0.19/4.6632 // 0.22/4.6632             //0.26//0.19 
objref Cell
Cell = new Cell_C()

Cell.soma.gbar_na3=0
Cell.soma.gbar_nap=0    //The sodium channels (Nap and Nat) were blocked to reveal HTO as in biology

objref stim, stim2, stim3, stim4, stim5, stim6, stim7, stim8, stim9, stim10
Cell.soma stim = new IClamp(0.5) //creating current clamp in the middle of the soma
stim.del = 300 // in ms
stim.dur = 2500//1200 // 600 // in ms
stim.amp = stim_amp // 0.123 // in nA

Cell.soma stim2 = new IClamp(0.5)
stim2.del = 400//350 // in ms
stim2.dur = 2400//1150 // 550 // in ms
stim2.amp = stim_amp // 0.123 // in nA

Cell.soma stim3 = new IClamp(0.5)
stim3.del = 500//400 // in ms
stim3.dur = 2300//1100 // 500 // in ms
stim3.amp = stim_amp // 0.123 // in nA

Cell.soma stim4 = new IClamp(0.5)
stim4.del = 600//450 //in ms
stim4.dur = 2200//1050 // 450 // in ms
stim4.amp = stim_amp // 0.123 // in nA

Cell.soma stim5 = new IClamp(0.5)
stim5.del = 700//500 // in ms
stim5.dur = 2100//1000 // 550 // in ms
stim5.amp = stim_amp // 0.123 // in nA

Cell.soma stim6 = new IClamp(0.5)
stim6.del = 800//550 // in ms
stim6.dur = 2000//950 // 550 // in ms
stim6.amp = stim_amp // 0.123 // in nA

Cell.soma stim7 = new IClamp(0.5)
stim7.del = 900//600 // in ms
stim7.dur = 1900//900 // 550 // in ms
stim7.amp = stim_amp // 0.123 // in nA

Cell.soma stim8 = new IClamp(0.5)
stim8.del = 1000//650 // in ms
stim8.dur = 1800//850 // 550 // in ms
stim8.amp = stim_amp // 0.123 // in nA

Cell.soma stim9 = new IClamp(0.5)
stim9.del = 1100//700 // in ms
stim9.dur = 1700//800 // 550 // in ms
stim9.amp = stim_amp // 0.123 // in nA

Cell.soma stim10 = new IClamp(0.5)
stim10.del = 1200//750 // in ms
stim10.dur = 1600//750 // 550 // in ms
stim10.amp = stim_amp // 0.123 // in nA

tstop = 3000//2500 // 1000 // in ms

SubVBoxNum = 3
objref MainHBoxObj,SubVBoxObj[SubVBoxNum]

proc MainBox() { local i
		MainHBoxObj = new HBox()
		for i=0,SubVBoxNum-1 SubVBoxObj[i] = new VBox()
		SubVBoxObj[0].intercept(1)
			newPlotVoltage()
			newPlotCurrent() 
			newPlotCurrent2()
		SubVBoxObj[0].intercept(0)
		SubVBoxObj[1].intercept(1)			
			access Cell.soma
			nrnsecmenu(.5,1)
		SubVBoxObj[1].intercept(0)
		SubVBoxObj[2].intercept(1)
			nrncontrolmenu()
			xpanel(" ")
				xbutton("QUIT","quit()")
			xpanel(0)
			nrnpointmenu(stim)	
		SubVBoxObj[2].intercept(0)
		MainHBoxObj.intercept(1)
			for i=0,SubVBoxNum-1 SubVBoxObj[i].map()
		MainHBoxObj.intercept(0)
		MainHBoxObj.map("Re-Tune",10,25,800,600)
}

proc newPlotVoltage() {
	newPlot(0,tstop,-100,50)
	graphItem.save_name("graphList[0].")
	graphList[0].append(graphItem)
	graphItem.addexpr("Cell.soma.v(.5)")
	graphItem.label(.08,.925,"mV")
}
proc newPlotCurrent() {
	newPlot(0,tstop,-0.2,0.2)
	graphItem.save_name("graphList[2].")
	graphList[1].append(graphItem)
	graphItem.addexpr("Cell.soma.i_kdr",1,2)
	//graphItem.addexpr("Cell.soma.i_ICa",2,2)
	//graphItem.addexpr("Cell.soma.i_iCat2",2,2)
	graphItem.label(.08,.925,tstr1)
}

proc newPlotCurrent2() {
	newPlot(0,tstop,-0.2,0.2)
	graphItem.save_name("graphList[2].")
	graphList[1].append(graphItem)
	graphItem.addexpr("Cell.soma.i_im",1,2)
	//graphItem.addexpr("Cell.soma.i_nap",2,2)
	//graphItem.addexpr("Cell.soma.i_kap",3,2)
	graphItem.addexpr("Cell.soma.i_kdr",4,2)
	//graphItem.addexpr("Cell.soma.i_sAHP",5,2)
	//graphItem.addexpr("Cell.soma.i_hd",6,2)
	//graphItem.addexpr("Cell.soma.i_na3",7,2)
	graphItem.addexpr("Cell.soma.i_cadyn",3,2)
	//graphItem.addexpr("Cell.soma.i_ICa",4,2)
	//graphItem.addexpr("Cell.soma.i_iCat2",4,2)
	graphItem.addexpr("Cell.soma.il_leak",9,2)
	graphItem.label(.08,.925,tstr1)
}

MainBox()
run()

// nrnpointmenu(stim)
// nrncontrolmenu()
// nrnsecmenu(.5,1)



// Plot1V("Cell")

// MainBox()
// init()
// run()

Loading data, please wait...