ModelDB is moving. Check out our new site at The corresponding page is

GP Neuron, somatic and dendritic phase response curves (Schultheiss et al. 2011)

 Download zip file 
Help downloading and running models
Phase response analysis of a GP neuron model showing type I PRCs for somatic inputs and type II PRCs for dendritic excitation. Analysis of intrinsic currents underlying type II dendritic PRCs.
1 . Schultheiss NW, Edgerton JR, Jaeger D (2010) Phase response curve analysis of a full morphological globus pallidus neuron model reveals distinct perisomatic and dendritic modes of synaptic integration. J Neurosci 30:2767-82 [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: Basal ganglia;
Cell Type(s): Globus pallidus neuron;
Gap Junctions:
Gene(s): HCN1; HCN2;
Simulation Environment: GENESIS;
Model Concept(s): Oscillations; Synchronization; Active Dendrites; Synaptic Integration; Parkinson's; Phase Response Curves;
Implementer(s): Edgerton, Jeremy R. [jedgert at]; Hanson, Jesse E.; Schultheiss, Nathan W [nwschultheiss at];
// Script to calculate the time integral for a single synaptic event given
//	tau1, tau2, and gmax.

function calc_syn_integ(gmax, tau1, tau2)
	float gmax, tau1, tau2
	float pktime = ({log {tau1}}-{log {tau2}})*{tau1}*{tau2}/({tau1}-{tau2})
	float pkval = ({gmax} / ({tau1} - {tau2}))	\
		* ({exp {-{pktime} / {tau1}}} - {exp {-{pktime} / {tau2}}})
	float A = gmax/{abs {pkval}}
	// Now have all components of synapse alpha function:
	//	gsyn = ((A * gmax)/(tau1-tau2)) * (exp(-t/tau1) - exp(-t/tau2))

	// Can now integrate alpha function over time. Stop integration when
	//	time reaches 5*tau2, where amplitude < 1% of peak.

	float syninteg = ({A} * {gmax}) / ({tau1} - {tau2})	\
		* (({tau2} * {exp -5}) - ({tau1} * {exp {-5*{tau2}/{tau1}}}) \
		- {tau2} + {tau1})

	return {syninteg}
        echo synInteg is {syninteg}

Loading data, please wait...