TITLE Basic sodium current COMMENT from "Gamma Oscillation by Synaptic Inhibition in a Hippocampal Interneuronal Network Model" (Wang and Buzsaki 1996) Used in Role of a Striatal Slowly Inactivating Potassion Current in Short-term Facilitation of Corticostriatal Inputs" A computer Simulation Study" (Mahon et al. 2000) Implemented by Kevin M. Biddell kevin.biddell@gmail.com 7/11/06 NOTE: 1S=1mho Neuron wants the units in mhos not millisiemens, please note the conversion! Phi =5 and no q10 or temp adjustment according to Bruno Delord 11/13/06 ENDCOMMENT UNITS { (mA) = (milliamp) (mV) = (millivolt) } NEURON { SUFFIX Nam USEION na WRITE ina RANGE gnabar, gna, minf, hinf, ah, Bh } INDEPENDENT {t FROM 0 TO 1 WITH 1 (ms)} PARAMETER { ena = 55 (mV) gnabar = 0.035 (mho/cm2) : 35mS phi = 5 < 0, 1e9 > : from delord 11/13/06 Vam = -28 :NOT the original value from wang and Buzsaki Kam = 1 Vbm = -53 :NOT the original value from wang and Buzsaki Kbm = 18 Vah = -51 :NOT the original value from wang and Buzsaki Kah = 20 Vbh = -21 :NOT the original value from wang and Buzsaki Kbh = 1 } STATE { m h } ASSIGNED { v (mV) ina (mA/cm2) celsius (degC) minf hinf ah Bh gna } BREAKPOINT { SOLVE states METHOD cnexp gna = gnabar*m^3*h ina = gna*(v - ena) } UNITSOFF INITIAL { rates(v) m = minf h= hinf } DERIVATIVE states { :Computes states variable m and h rates(v) : at the current v and dt. h'=phi*(ah*(1-h)-Bh*h) } PROCEDURE rates(v) { :Computes rate and other constants at current v. :Call once from HOC to initialize inf at resting v. LOCAL am, Bm am = (-0.1*(v-Vam)/Kam/(exp(-0.1*(v-Vam)/Kam)-1)) Bm = 4*exp(-(v-Vbm)/Kbm) ah = 0.07*exp(-(v-Vah)/Kah) Bh= 1/(1+exp(-0.1*(v-Vbh)/Kbh)) minf = am/(am+Bm) hinf = ah/(ah+Bh) m=minf } UNITSON