A model for interaural time difference sensitivity in the medial superior olive (Zhou et al 2005)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:53869
This model simulates responses of neurons to interaural time difference (ITD) in the medial superior olive (MSO) of the mammalian brainstem. The model has a bipolar cell structure and incorporates two anatomic observations in the MSO: (1) the axon arises from the dendrite that receives ipsilateral inputs and (2) inhibitory synapses are located primarily on the soma in adult animals. Fine adjustment of the best ITD is achieved by the interplay of somatic sodium currents and synaptic inhibitory currents. The model suggests a mechanism for dynamically "fine-tuning" the ITD sensitivity of MSO cells by the opponency between depolarizing sodium currents and hyperpolarizing inhibitory currents.
Reference:
1 . Zhou Y, Carney LH, Colburn HS (2005) A model for interaural time difference sensitivity in the medial superior olive: interaction of excitatory and inhibitory synaptic inputs, channel dynamics, and cellular morphology. J Neurosci 25:3046-58 [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): Medial Superior Olive (MSO) cell;
Channel(s): I h; I Sodium; I Potassium;
Gap Junctions:
Receptor(s): AMPA; Glycine;
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Coincidence Detection; Axonal Action Potentials;
Implementer(s): Zhou, Yi [yizhou at bu.edu];
Search NeuronDB for information about:  AMPA; Glycine; I h; I Sodium; I Potassium;
TITLE Ih channels in VCN auditory neurons
: ih=ghmax*r*(v-eh) 
: based on Rothman and Manis (2003c)
: Modifications by Yi Zhou for an MSO model

INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}

NEURON {
	SUFFIX ih_VCN2003
	USEION h READ eh WRITE ih VALENCE 1	
	RANGE ghbar 
	RANGE r_inf,tau_r,r_exp
	RANGE ih,gh
	
}


UNITS {
	(mA) = (milliamp)
	(mV) = (millivolt)
}

PARAMETER {
	ghbar	= 0.002	(mho/cm2)  
	eh=-43		(mV) 
	celsius =22		(degC)
	dt              (ms)
	v               (mV)
	
}

STATE {
	r
}

ASSIGNED {
	gh	(mho/cm2)
	ih	(mA/cm2)
	r_inf
	tau_r
	r_exp
	tadj
}


BREAKPOINT {
	SOLVE states
	gh=ghbar *r
	ih  = gh*(v-eh)
}


PROCEDURE states() {	: this discretized form is more stable
	evaluate_fct(v)
	r = r + r_exp * (r_inf - r)
	VERBATIM
	return 0;
	ENDVERBATIM
}

UNITSOFF
INITIAL {
:
:  Q10 was assumed to be 3 for both currents
:
	tadj = 3.0 ^ ((celsius-22)/ 10 )
	evaluate_fct(v)
	r= r_inf
	}

PROCEDURE evaluate_fct(v(mV)) {
	
	tau_r = (100000/(237*exp((v+60)/12)+17*exp(-(v+60)/14))+25)/ tadj
	r_inf = 1/(1+exp((v+76)/7))
   
	r_exp = 1 - exp(-dt/tau_r)
	
}

UNITSON

Loading data, please wait...