This is the model used in:

Chambers, J. D., J. C. Bornstein, et al. (2014). "A detailed,
conductance based computer model of intrinsic sensory neurons of the
gastrointestinal tract." Am J Physiol Gastrointest Liver Physiol. pii:
ajpgi.00228.2013. [Epub ahead of print]

These simulations require the NEURON simulator available for free from

To run the model, first compile all the files in the ./mod directory
(for example, type "nrnivmodl ./mod" at the command line after
changing into in the directory containing the model). Then run one of
the relevent (fig1.hoc, fig2.hoc, etc.) hoc files by typing "nrngui
fig1.hoc" at the command line or by opening neuron and using the
load_file function.

Each of the fig#.hoc files will produce a graph of the membrane
potential (or membrane current for Fig5.hoc) and produce 2 output
files. The output file fig#mem.txt contains time in the first column
and membrane potential in the second column. The output file
fig#currents.txt contains time in the first first column and currents
or concentrations in the subsequent columns (see individual hoc file
for details).

Fig1.hoc will also produce a second graph of the currents in the
model (1A, 1B):
screenshot 1A
screenshot 1B

Fig2.hoc will produce Panel C2:
screenshot 2C2

Comment out the stimuli in the hoc to produce Panels A2 and
B2. Uncomment conductance changes highlighted in the hoc file to
produce Panel D2.

Fig3.hoc can take command line arguments. For example, to change the
strength of the depolarising current pulse use "nrngui -c stimamp=150
fig3.hoc" for a 150 pA depoalrising current or "nrngui -c stimamp=350
fig3.hoc" for a 350 pA depolarising current. The total conductance of
all currents (and parameters of Ih varied in Figure 8) can also be
altered using command line arguments. For example, "nrngui -c
gnav13=500 fig3.hoc" will change the total conductance of Nav1.3 (for
a full list of variables see Define default variables section of the
fig3.hoc). The value of the variable is divided by 1000 (or 100 for
properties of Ih another than total conductance, see conductance
changes section of fig3.hoc for details) and then multiplied with the
baseline value for that current. So, passing a value of 500 for gnav13
decreases the total conductance to half, whereas passing a value of
2000 will double the total conductance. Without modification Fig3.hoc
produces a graph like Fig3C2:

screenshot 3C2