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,NOI] = excitation(Imax,Nmsi,Itol,noAP,tspan,model,stim)
%EXCITATION Estimate the excitation threshold
%   [E,NOI,y1] = excitation(Imax,Nmsi,Itol,noAP,tspan,model,stim) this function
%   estimates the excitation threshold for a nerve fiber [model] to a stimulus
%   [stim]. The timespand for the simulations are [tspan] = [tstart tend]. 
%   The threshold is estimated by a binary search algorithm, where the Itop 
%   first are determined. The intial guess for Itop are [Imax], and this 
%   value are used unless it fails to evoke an action potential (AP). If it fails
%   to evoke an AP a linear downwards seach are used, with [Nmsi] steps. If none
%   of these evokes an AP the function returns NaN. The excitation thresholds 
%   are estimated with an precision of [Itol], and the number of simulations
%   are returned in NOI (No of iterations).

%Find the value which produce an action potential.
Etest = Imax:-Imax/Nmsi:Imax/Nmsi;
Et = NaN;
NOI = 0;
n = Nmsi;

%fprintf('Excitation determination: ');
while isnan(Et) & n > 0
    NOI = NOI + 1;
    if isAP(-Etest(n),noAP,tspan,model,stim)
        Et = Etest(n);
    n = n-1;

Eb = 0;

if ~isnan(Et)
    while (Et - Eb) > Itol
        NOI = NOI + 1;
        Etest = Eb+(Et-Eb)/2;
        if isAP(-Etest,noAP,tspan,model,stim)
           Et = Etest;
           Eb = Etest;
    E = Et;
    %fprintf(' = %.2fnA\n',E*1e9);
    E = NaN;
    %fprintf(' = failed\n');

Loading data, please wait...