ModelDB is moving. Check out our new site at The corresponding page is

Breakdown of accmmodation in nerve: a possible role for INAp (Hennings et al 2005)

 Download zip file 
Help downloading and running models
The present modeling study suggests that persistent, low-threshold, rapidly activating sodium currents have a key role in breakdown of accommodation, and that breakdown of accommodation can be used as a tool for studying persistent sodium current under normal and pathological conditions. See paper for more and details.
1 . Hennings K, Arendt-Nielsen L, Andersen OK (2005) Breakdown of accommodation in nerve: a possible role for persistent sodium current. Theor Biol Med Model 2:16 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuron or other electrically excitable cell;
Brain Region(s)/Organism:
Cell Type(s): Spinal cord lumbar motor neuron alpha ACh cell; Myelinated neuron;
Channel(s): I Na,p; I Na,t; I K;
Gap Junctions:
Simulation Environment: MATLAB;
Model Concept(s): Action Potential Initiation; Action Potentials; Pathophysiology; Electrotonus;
Implementer(s): Hennings, Kristian [krist at];
Search NeuronDB for information about:  Spinal cord lumbar motor neuron alpha ACh cell; I Na,p; I Na,t; I K;
function [E,E0] = te_m(M,TD,P,Idc)
%TE Threshold Electrotonus
%   [E,E0] = te_m(M,TD,P,Idc)

Imax = 10e-9; Nmsi = 5; Itol = 0.0001e-9; noAP = 1; 

[S] = pulse(0,1e-3);  
S = setDC(S,0);
E0 = excitation(Imax,Nmsi,Itol,noAP,[0 2e-3],M,S);        

for n = 1:length(P)
    for k = 1:length(TD)
        %fprintf('Threshold Electrotonus P = %.2f Tdelay = %.2fms\n',P(n),TD(k)*1e3);
        M = MSAVE;
        S = pulse(-E0*P(n),100e-3);   
        S = setDC(S,Idc);
        if TD(k) ~= 0
            M.X0 = x0patch([0 TD(k)],S,M);
        tspan = [TD(k) TD(k) + 2e-3];
        S = cpulse(0, 1e-3, TD(k), -E0*P(n), 100e-3);
        S = setDC(S,Idc);
        E(n,k) = excitation(Imax,Nmsi,Itol,noAP,tspan,M,S);        

Loading data, please wait...