Well, first of all I'd like to say that this is my first approach to NMODL and I'm currently in the learning process...
Now, the issue is that I'm trying to insert some conductances in a model but most of the available examples on the internet are wrote in terms of the alpha and beta functions. My problem is that the conductances I'm trying to insert are expressed in terms of the Hudgkin & Huxley model for ionic currents, according to Belluzzi et al., 1991:
Ina=Ina'[1-exp(-t/taum)]^3[exp(-t/tauh)]
My doubt is about how my code file should be structured in order to work properly? Specifically, if in this case I should include an INITIAL and DERIVATIVE blocks or go directly to the PROCEDURE block?
Basically, I'm wondering how my code file should look like and what to include if I'm working with this kind data.
At this point this is what I have, but before I move on I'd like to ask for your expertise advice.
Code: Select all
TITLE Sodium conductance for sympathetic neuron
UNITS {
(ms) = (millisecond)
(mV) = (millivolt)
(mA) = (milliamp)
(mS) = (millisiemen)
}
NEURON {
SUFFIX Ina
USEION na READ Ina2 tm th WRITE Ina
}
PARAMETER {
Vh = -34.45 (mV)
k = 6.99 (mV)
gna = 0.00928 (mS)
}
ASSIGNED {
tm (ms)
th (ms)
m (mV)
Ina2 (mA)
Ina (mA)
}
STATE {
Ina2 tm th
}
BREAKPOINT {
SOLVE states METHOD cnexp
Ina = Ina2*((1-exp(-t/tm))*(1-exp(-t/tm))*(1-exp(-t/tm)))*(exp(-t/th))
}
: I don't know where to put this
tm = 0.06+(1/((42.98*exp(0.08915*V))+(0.9230*exp(-0.03351*V)))
th = ((-0.0046*V+0.26)+50.85)/((1+exp((-59.46-V)/7.91))+(1+exp((V+40.94)/1.556)))
m = (1/(1+exp((Vh-V)/k)))
:
Thanks guys.