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

 Download zip file 
Help downloading and running models
Accession:55749
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.
Reference:
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:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Action Potential Initiation; Action Potentials; Pathophysiology; Electrotonus;
Implementer(s): Hennings, Kristian [krist at hst.auc.dk];
Search NeuronDB for information about:  Spinal cord lumbar motor neuron alpha ACh cell; I Na,p; I Na,t; I K;
function [A,S,E] = rcurve(M,TS,Idc)
%ACURVE Accommodation curve
%   [A,S,E] = rcurve(M,TS,R,Idc) this function simulate the accommodation curve [A]
%   for the durations [TS] of ramp. [A] has the same 
%   length as [TS]. The rheobase [R] is used for calculating the accommodation
%   curve [A] and the accommodation slopes [S]. The raw excitation thresholds are
%   returned in [E].

%Create the parameters for the excitation function
Imax = 10e-9; Nmsi = 5; Itol = 0.0001e-9; noAP = 1; 

S = pulse(0,100e-3);
S = setDC(S,Idc);

R = excitation(Imax,Nmsi,Itol,1,[0 101e-3],M,S);

for n = 1:length(TS)
    stim = ramp(0,TS(n));
    stim = setDC(stim,Idc);
    tspan = [0 TS(n)+1e-3];
    E(n) = excitation(Imax,Nmsi,Itol,noAP,tspan,M,stim);
    S(n) = (E(n) / TS(n)) / R;
    %fprintf('.');
end
A = E / R;

Loading data, please wait...