MRG model & bipolar stimulation - asymmetric response?

Discussions of particular models.

Moderator: tom_morse

Shamsie
Posts: 1
Joined: Wed Sep 20, 2017 7:25 am

MRG model & bipolar stimulation - asymmetric response?

Hello,

I have implemented the MRG model in Python and was able to reproduce some simple results. I use it mainly by applying a vector of extracellular potential along the axon (using e_extracellular).

However, I observed an unexpected outcome with bipolar stimulation and I cannot pinpoint its origin, whether it is an implementation problem, a misunderstanding or something else:

Let's say I apply an extracellular potential (see below) induced by a bipolar stimulation of of a given amplitude (1.3 mA), and the axon (200 mm long, 11.5 um diameter) does not fire.
https://ibb.co/gOYKik
However if I apply the exact same vector of extracellular potential induced by the same stimulation, but I inverse the polarity (see below), then the axon gets activated (an action potential propagates along the axon).
https://ibb.co/fgdKik

I expected a symmetric behavior, i.e. the same threshold amplitude for activation with both extracellular potential vectors. It turns out the difference in the amplitude required to activate the axon is about 0.2 mA, which I do not believe can be attributed to interpolation or other approximations.
Here are the graphs of the membrane potential with the initial and inverted bipolar stimulation:
https://ibb.co/mz2bw5
https://ibb.co/ekO5pQ

and their corresponding GIF with membrane potential over time:
https://ibb.co/dqSn3k
https://ibb.co/czKjik

Here are more details about what I tried to highlight potential cues of where this difference comes from, but in vain:
I made sure the extracellular potential vectors have exactly symmetric peaks (I just duplicated half of the vector and inverted it before concatenating it with the first half), and that the axon is long enough (200 mm) compared to the internode length (1.3 mm) and the extracellular potential decay on the tails. I also tried to shift this bipolar shape of the extracellular potential along the axon, swept a range of axon node offsets, none of which influenced the difference.

I would be so grateful if anyone could give me any lead, clue, idea, anything that could help me understand why I see this difference, or how I could at least identify it as an error of implementation or as an intrinsic characteristic of the model, because I run out of ideas.

I thank you in advance,
Shamsie
ted