e_pas for every compartment set to 112.55mV in the model!!!
I just ran the H_current.hoc file in the model given at modeldb: https://senselab.med.yale.edu/modeldb/s ... 212#tabs1
When I printed e_pas using the following command in hoc: soma[0].e_pas(0.5) , It showed 112.55mV. I know this was done by the currentbalance.hoc file to initialize the model to v_init, but this seems unphysiological.
Furthermore, if you calculate the input capacitance (current calmp of 10pA at soma, then calculate input resistance  ~65Mohms and tau  ~15ms. Then calculated Cm), it's coming out to be 420pF!! (Neuroelectro and my inlab data says it should be around 65pF)
Can anybody help me with this?
e_pas in Poirazi2003 CA1 pyr cell model
Moderator: tom_morse

 Site Admin
 Posts: 5632
 Joined: Wed May 18, 2005 4:50 pm
 Location: Yale University School of Medicine
 Contact:
Re: e_pas in Poirazi2003 CA1 pyr cell model
It was done to force the resting potential of the model to be v_init everywhere. The basic idea is expressed in this algorithm:e_pas for every compartment set to 112.55mV in the model!!! . . . I know this was done by the currentbalance.hoc file to initialize the model to v_init, but this seems unphysiological.
Code: Select all
initialize the model to v_init
for each compartment in the model
calculate the local net ionic membrane current
add a local constant membrane current of equal magnitude but opposite sign
One way to accomplish this result would be to insert into all sections a density mechanism that delivers a constant current, and adjust the value of this current in each compartment to force local charge balance.
The authors of the model you cite accomplished the same result by taking advantage of the fact that every compartment in their model had a mechanism with a linear constant conductance (the pas mechanism, in this particular case). Then they probably did something equivalent to
Code: Select all
initialize the model to v_init
for each compartment in the model
calculate the sum of all local ionic membrane currents EXCEPT i_pas
set local e_pas to a value that makes local i_pas equal to 1 times that sum
"Can you prove that this second approach, i.e. changing e_pas in each compartment, is equivalent to the first approach, i.e. adding a constant current to each compartment?"
Yes. Consider the model with unperturbed e_pas. The equivalent circuit for each compartment includes a branch with voltage source e_pas in series with a conductance g_pas
Code: Select all
og_pase_paso
which can be analyzed like any other single port linear circuit. In particular, it has an open circuit voltage, a closed circuit current, and a source resistance, the values of which are e_pas, g_pas*e_pas, and 1/g_pas, resepectively.
Now consider the model with perturbed e_pas. The equivalent circuit for each compartment has a branch with voltage source e_pas+delta in series with conductance g_pas
Code: Select all
og_pase_pas+deltao
But from linear circuit theory we know that this is equivalent to the parallel combination of the original
Code: Select all
og_pase_paso
Good question. Experimental measurements are prone to all kinds of vagaries, and so is computational modeling. Did Poirazi et al. say where they got the anatomical and biophysical data on which they based their model? What species were used, what were the ages of the animals? What about for the data you found in NeuroElectro? Were the experiments done in slices, and if so, were the dendritic trees intact?if you calculate the input capacitance (current calmp of 10pA at soma, then calculate input resistance  ~65Mohms and tau  ~15ms. Then calculated Cm), it's coming out to be 420pF!! (Neuroelectro and my inlab data says it should be around 65pF)