Page 1 of 1
Sum of three exponential
Posted: Wed Jul 08, 2009 3:59 pm
by Bill Connelly
Hi,
I'm trying to figure out the differential equations that explain the recovery of a process. If the recovery from an inactive state to an active state is explained by the sum of three exponentials (with each time constant being an order of magnitude difference from the last), am I best to think of the process in the following way
Code: Select all
Inactive ---> Active
Changed to:
slow rate
Inactive ---> A
medium rate
A ---> B
fast rate
B ---> Active
I have a feeling I'm completely wrong, so maybe there is a section in the NEURON book I've missed, or some other resource you can point me to?
Re: Sum of three exponential
Posted: Wed Jul 08, 2009 5:02 pm
by ted
That would work regardless of the relative magnitudes of the rate constants. Note that the rate constants and time constants do not have a simple relationship. Rewriting with a slightly different notation, the reactions are
so
Code: Select all
Dt I = -kI I
Dt W = kI I - kW W
Dt X = kW W - kX X
where Dt means d/dt
Re: Sum of three exponential
Posted: Fri Jul 10, 2009 12:46 am
by Bill Connelly
So I have two versions of a mod file that it my mind should produce identical results, but the first one works (the current decays back to zero) and the second one doesn't. Are you able to explain to me why not?
Code: Select all
NEURON {
POINT_PROCESS expi
NONSPECIFIC_CURRENT i
RANGE i, e, g, gbar, kA, kB, kC
}
PARAMETER {
e = -50 (millivolts)
gbar = 0.01 (microsiemens)
kA = 3 (/ms)
kB = 6 (/ms)
kD = 9 (/ms)
}
STATE {
O
A
B
D
C
}
INITIAL {
C = 1
A = 0
B = 0
D = 0
O = 0
}
ASSIGNED {
v (millivolt)
i (nanoamp)
g (microsiemens)
}
BREAKPOINT {
SOLVE state METHOD cnexp
O = 0.2 * B + 0.4 * A + 0.4 * D
g = gbar * O
i = g * (v - e)
}
DERIVATIVE state {
A' = -kA * A
B' = -kB * B
D' = -kD * D
C' = 0.2 * (kB * B) + 0.4 * (kA * A) + 0.4 * (kD * D) :Give that the code below doesn't work, I get the fact that this probably doesn't represent the closed state
}
NET_RECEIVE(weight (microsiemens)) {
O = O + weight
A = O
B = O
D = O
C = 1 - O
}
Code: Select all
...
BREAKPOINT {
SOLVE state METHOD cnexp
g = gbar * O
i = g * (v - e)
}
DERIVATIVE state {
A' = -kA * A
B' = -kB * B
D' = -kD * D
C' = 0.2 * (kB * B) + 0.4 * (kA * A) + 0.4 * (kD * D)
O' = -0.2 * (kB * B) - 0.4 * (kA * A) - 0.4 * (kD * D)
}
...
Re: Sum of three exponential
Posted: Fri Jul 10, 2009 6:59 am
by Bill Connelly
Either way, I came up with an acceptable model of synaptic depression. I'm sure it isn't as good as the other ones that are available in modelDB, but the difference is, I understand it. I'm just posting it up here because I thought someone else might find it useful
Code: Select all
NEURON {
POINT_PROCESS expi
NONSPECIFIC_CURRENT i
RANGE i, e, tau, kA, kB, kC, dec
}
PARAMETER {
e = 0 (millivolts)
tau = 4 (ms)
kA = 0.00025 (/ms)
kB = 0.005 (/ms)
kC = 0.025 (/ms)
dec = 0.5
}
STATE {
O
A
B
C
F
g
}
INITIAL {
A = 0
B = 0
C = 0
O = 1
g = 0
}
ASSIGNED {
v (millivolt)
i (nanoamp)
}
BREAKPOINT {
SOLVE state METHOD cnexp
F = A/3 + B/3 + C/3
O = 1 - F
i = g * (v - e)
}
DERIVATIVE state {
A' = -kA * A
B' = -kB * B
C' = -kC * C
g' = -g/tau
}
NET_RECEIVE(weight (microsiemens)) {
g = (g + weight) * O
O = O * dec
F = 1 - O
A = F
B = F
C = F
}