Page 1 of 1
Specifying a finite volume of Extracellular fluid
Posted: Mon Jan 23, 2012 9:30 am
by physrob
Hi all
I am very new to Neuron, so apologies if I am asking a stupid question.
Does anyone know of a way to specify a finite volume of extracellular fluid?
Thanks
Rob
Re: Specifying a finite volume of Extracellular fluid
Posted: Tue Jan 24, 2012 10:29 am
by ted
You want a model of ion accumulation in extracellular space. Doable with NMODL, for example see the implementation of Frankenhaeuser-Hodgkin space in chapter 9 of The NEURON Book. For your convenience, I reproduce kext.mod here
Code: Select all
: Extracellular potassium ion accumulation
: fhspace and kxfer are from Frankenhaeuser & Hodgkin
NEURON {
SUFFIX kext
USEION k READ ik WRITE ko
GLOBAL kbath
RANGE fhspace, txfer
}
UNITS {
(mV) = (millivolt)
(mA) = (milliamp)
FARADAY = (faraday) (coulombs)
(molar) = (1/liter)
(mM) = (millimolar)
}
PARAMETER {
kbath = 2.5 (mM)
: Frankenhaeuser-Hodgkin space parameters
fhspace = 300 (angstrom) : thickness
txfer = 50 (ms) : fhspace <-> bath transfer time
}
ASSIGNED {
ik (mA/cm2)
}
STATE {
ko (mM)
}
BREAKPOINT {
SOLVE state METHOD derivimplicit
}
DERIVATIVE state {
ko' = (1e8)*ik/(fhspace*FARADAY) + (kbath - ko)/txfer
}
An aside: astute readers may notice that the METHOD is now derivimplicit rather than cnexp. derivimplicit is appropriate for ion accumulation mechanisms specified by ODEs and will work regardless of the equilibration time constant.
Re: Specifying a finite volume of Extracellular fluid
Posted: Wed Jan 25, 2012 3:33 am
by physrob
Thanks Ted
I will have a good look at the relevant sections in the NEURON book.
Rob