How does NEURON calculate the membrane potential?
Posted: Tue Apr 25, 2006 5:30 am
How does NEURON calculate the membrane potential? Is there any diffrences between ion current (such like 'ica') and nonspecific current i?
We have build a two compartement model containing soma and dend. In the dend we have inserted r-Type voltage gated channels 'car.mod' by the commend "insert car", and AMPA receptors. To trigger the AMPA receptors we used 'NetCon()' function.
By plotting the membrane potential and the currents we found although the r-Type VGCC (ica) is big, it doesn't seem to contribute as much as it should to the membrane potential. In the plot the ion current from r-Type VGCC (ica) is as big as the the nonspecific AMPA current (i) but they do not contrubute the same. Also by blocking r-Type VGCC (ica) there is nearly no difference in the membrane potential. Both the ion current and the non specific current are in current density (mA/cm^2).
Could the problem be the way NEURON handles the two kinds of current (nonspecific and ion current)?
the main code for the two compartment model
--------------------------------------------------------------
...
access soma
access dend
celldef()
forsec "dend" {
insert car
insert canmda
insert cad
}
dend[0].gcabar_car=5
......
objref synAmpa
dend[0] synAmpa = new ampa(0.5)
objref nc1, nil
nc1 = new NetCon(nil, synAmpa)
nc1.weight = 10
nc1.delay = 0
objref syntimes
syntimes = new Vector(1)
syntimes.x[0] = 10
objref fih
fih = new FInitializeHandler("loadqueue()")
proc loadqueue() { local ii
for ii=0,syntimes.size()-1 {
nc1.event(syntimes.x[ii])
}
}
------------------------------------------------
We have build a two compartement model containing soma and dend. In the dend we have inserted r-Type voltage gated channels 'car.mod' by the commend "insert car", and AMPA receptors. To trigger the AMPA receptors we used 'NetCon()' function.
By plotting the membrane potential and the currents we found although the r-Type VGCC (ica) is big, it doesn't seem to contribute as much as it should to the membrane potential. In the plot the ion current from r-Type VGCC (ica) is as big as the the nonspecific AMPA current (i) but they do not contrubute the same. Also by blocking r-Type VGCC (ica) there is nearly no difference in the membrane potential. Both the ion current and the non specific current are in current density (mA/cm^2).
Could the problem be the way NEURON handles the two kinds of current (nonspecific and ion current)?
the main code for the two compartment model
--------------------------------------------------------------
...
access soma
access dend
celldef()
forsec "dend" {
insert car
insert canmda
insert cad
}
dend[0].gcabar_car=5
......
objref synAmpa
dend[0] synAmpa = new ampa(0.5)
objref nc1, nil
nc1 = new NetCon(nil, synAmpa)
nc1.weight = 10
nc1.delay = 0
objref syntimes
syntimes = new Vector(1)
syntimes.x[0] = 10
objref fih
fih = new FInitializeHandler("loadqueue()")
proc loadqueue() { local ii
for ii=0,syntimes.size()-1 {
nc1.event(syntimes.x[ii])
}
}
------------------------------------------------