Dynamical model of olfactory bulb mitral cell (Rubin, Cleland 2006)

 Download zip file   Auto-launch 
Help downloading and running models
Accession:64296
This four-compartment mitral cell exhibits endogenous subthreshold oscillations, phase resetting, and evoked spike phasing properties as described in electrophysiological studies of mitral cells. It is derived from the prior work of Davison et al (2000) and Bhalla and Bower (1993). See readme.txt for details.
Reference:
1 . Rubin DB, Cleland TA (2006) Dynamical mechanisms of odor processing in olfactory bulb mitral cells. J Neurophysiol 96(2):555-568 [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): Olfactory bulb main mitral cell;
Channel(s): I Na,p; I Na,t; I L high threshold; I A; I K; I K,leak; I h; I K,Ca; I Sodium; I Calcium; I Potassium; I A, slow;
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: NEURON;
Model Concept(s): Bursting; Temporal Pattern Generation; Oscillations; Synchronization; Simplified Models; Delay; Olfaction;
Implementer(s): Rubin, Daniel B ;
Search NeuronDB for information about:  Olfactory bulb main mitral cell; I Na,p; I Na,t; I L high threshold; I A; I K; I K,leak; I h; I K,Ca; I Sodium; I Calcium; I Potassium; I A, slow;
/
mitral
readme.txt *
cadecay.mod *
Ih.mod *
INaP.mod
kA.mod
kca3.mod
kfasttab.mod *
kO.mod *
kslowtab.mod *
lcafixed.mod *
nafast.mod
kfast_k.inf *
kfast_k.tau *
kfast_n.inf *
kfast_n.tau *
kslow_k.inf *
kslow_k.tau *
kslow_n.inf *
kslow_n.tau *
mitral.hoc
mosinit.hoc *
tabchannels.hoc *
                            
TITLE Hyperpolarization-activated cation current
: Implemented in Rubin and Cleland (2006) J Neurophysiology
: Adapted from Saraga et al (2003) J. Physiology

UNITS {
        (mA) = (milliamp)
        (mV) = (millivolt)
}
 
NEURON {
        SUFFIX Ih
        USEION h READ eh WRITE ih VALENCE 1
        RANGE gkhbar,ih
        GLOBAL rinf, rexp, tau_r
}
 
INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)}
 
PARAMETER {
        v (mV)
        p = 5 (degC)
        dt (ms)
        gkhbar = 0.00 (mho/cm2)			
        eh = -32.9 (mV)
}
 
STATE {
        r
}
 
ASSIGNED {
        ih (mA/cm2)
	rinf rexp
	tau_r
}
 
BREAKPOINT {
        SOLVE deriv METHOD derivimplicit
        ih = gkhbar*r*(v - eh)
}
 
INITIAL {
	rates(v)
	r = rinf
}

DERIVATIVE deriv { :Computes state variable h at current v and dt.
	rates(v)
	r' = (rinf - r)/tau_r
}

PROCEDURE rates(v) {  :Computes rate and other constants at current v.
                      :Call once from HOC to initialize inf at resting v.
        TABLE rinf, rexp, tau_r DEPEND dt, p FROM -200
TO 100 WITH 300
	rinf = 1/(1 + exp((v+84.1)/10.2))
	rexp = 1 - exp(-dt/(tau_r))
	tau_r = 100 + 1/(exp(-17.9-0.116*v)+exp(-1.84+0.09*v))
}
 
UNITSON


Loading data, please wait...