Page 1 of 1
current between two compartments
Posted: Wed Mar 07, 2007 12:01 pm
by thats_karlo
Hi,
Suppose we have a 2-compartment model with different currents in these two compartments. we inject current and measure membrane potential in the first one.
consider this case: we have injected a constant depolarized current, for a specific duration, in the middle of fist compartment
how can we determine,
1- the value of (total/individual) current that reach to second one and depolarize it?
(after some times the feedback current)
2- the value of (total/individual) current coming from second one to the first compartment?
Thanks in advance,
Karlo
Posted: Wed Mar 07, 2007 12:07 pm
by ted
Posted: Wed Mar 07, 2007 12:15 pm
by thats_karlo
Dear Ted,
Thanks for very fast reply!!
Posted: Wed Mar 07, 2007 12:24 pm
by ted
When you're not sure about the meaning of something in the Programmer's Reference,
do an experiment that removes the ambiguity. In this case, make a toy model. Calculate
the values by hand. Compare with what you get from ri().
Posted: Wed Mar 07, 2007 4:57 pm
by thats_karlo
Dear Ted,
please tell me am i right or no?
I have two compartments model, different currents, with following topology
________
| dend[0] |__________
| dend[0] | dend[1] |
|_______|--------------
and forall nseg=3
dose
Code: Select all
dend[1] { ( v(0.16)- v(0) )/ri(0.16) }
means axial current between to compartments?
Posted: Thu Mar 08, 2007 9:55 am
by ted
thats_karlo wrote:I have two compartments model, different currents, with following topology
________
| dend[0] |__________
| dend[0] | dend[1] |
|_______|--------------
and forall nseg=3
Then it isn't a two compartment model, is it?
Code: Select all
dend[1] { ( v(0.16)-.v(0) )/ri(0.16) }
means axial current between to compartments?
Absolutely correct. In that expression, ri(0.16) is the axial resistance between the first
internal node of dend[1] and the node at the 0 end of dend[1]. v(0.16) and v(0) are
the membrane potentials at those nodes. Therefore v(0.16) - v(0) is the voltage gradient
across the axial resistance, and the ratio v(0.16) - v(0) / ri(0.16) is the axial current that
enters the 0 end of dend[1]. This is equal to the axial current that leaves the 1 end of
dend[0].
But what are the units of this current?
Posted: Thu Mar 08, 2007 10:48 am
by thats_karlo
Thanks Ted,
1- probably, i should say, we have two sections, with nseg=3, which means 6-compartment,right?
2-I think the unite for axial current is nA, right?
Yours,
Karlo
Posted: Thu Mar 08, 2007 12:47 pm
by ted
True and true. Good work, Karlo.
Posted: Tue Mar 13, 2007 10:02 am
by thats_karlo
Hello Ted,
Could you please tell me the following formula is a right method to calculate total current that depolarized cell,(in the soma to reach threshold to firing); I=I1+I2+I3, where
I1=(Ionic current+leak current)*area(0.5)/100 // (like Na+/K+, Ca2+,...) in [nA]
I2=(v(0.5)-v(1))/ri(1) // axial current toward to soma
I3=Iinj // injected current
Karlo
Posted: Tue Mar 13, 2007 4:42 pm
by ted
That would work, as long as the default section has nseg = 1. If nseg is not 1, you need
to do
totalarea = 0
for (x, 0) totalarea += area(x)
and then change the equation for I1 to
I1=(Ionic current+leak current)*totalarea/100 // (like Na+/K+, Ca2+,...) in [nA]