Motoneuron simulations for counting motor units (Major and Jones 2005)

 Download zip file 
Help downloading and running models
Accession:53425
Simulations of clinical methods to count the number of motoneurons/motor units in human patients. Models include stimulation of motor axons or voluntary activation and responses are measured as muscle tension or EMG.
Reference:
1 . Major LA, Jones KE (2005) Simulations of motor unit number estimation techniques. J Neural Eng 2:17-34 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neuromuscular Junction;
Brain Region(s)/Organism:
Cell Type(s): Spinal cord lumbar motor neuron alpha ACh cell;
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Simplified Models; Methods;
Implementer(s): Major, Lora [lora.major at ualberta.ca];
Search NeuronDB for information about:  Spinal cord lumbar motor neuron alpha ACh cell;
function [ltw,emg,nrec] = IM_pool(excite,p,ttime,etime)
%%--------------------------------------------------------
%% this function simulates synchronous stimulation of a motor neuron pool 
%% inputs are:
%%       excite = absolute level of excitation
%%       p = params of pool
%%       ttime = time vector for force output
%%       etime = time vector for emg output
%% outputs are:
%%       ltw = muscle force obtained from stimulus
%%       emg = emg output
%%       nrec = number of units recruited for each level of excitation

twlin = zeros(length(ttime),p.n);
emglin = zeros(length(etime),p.n);
nrec = 0;

%% sim for each MN
for i=1:p.n
    
    if excite >= p.rte(i)
        %% do twitch function
        twlin(:,i) = twitch(p.del,ttime,p.twtforce(i),p.tc(i));
        emglin(:,i) = muap2(p.del,etime,p.twtforce(i));
        nrec = nrec + 1;
    end
    
end

ltw = sum(twlin,2)';
emg = sum(emglin,2)';

return

Loading data, please wait...