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]
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
                            
// control

// cholinergic


for b=0, 31 {
		// sweep over parameters:
		
		/*
		Iinject.amp = 0.08 + b*0.5/31 
		print "current injection (nA): ", Iinject.amp
		
		gkbar_iM = b*0.0006/31 
		print "g_iM (Scm-2): ", gkbar_iM
		
		e_pas = -62 + b*8/31
		print "passive reversal (mV): ", e_pas 
		*/
		g_pas = 1.6e-5 + b*6e-6/31
		print "passive conductance (Scm-2): ", g_pas
		
		
		// sweep gm, keep f rate constant by adjusting Iinj
		// f curve vs drive
		makeif()
		
		targetf=10
		dif=100	
		for d=0, frequency.size-1 {
			if (abs(frequency.x[d]-targetf)<dif) {
							dif=abs(frequency.x[d]-targetf)	
							Itarget=d*range/sweeps	
								}	
			}
		
		Iinject.amp=Itarget
		
		soma spikesout = new NetCon(&v(.5), nix, 0, 0, 0)// Output spike train
		spikesouttimes = new Vector()
		spikesout.record(spikesouttimes)
		makeprc()
}


Loading data, please wait...