Modelling enteric neuron populations and muscle fed-state motor patterns (Chambers et al. 2011)

 Download zip file 
Help downloading and running models
Accession:136024
"After a meal, the gastrointestinal tract exhibits a set of behaviours known as the fed state. ... Segmentation manifests as rhythmic local constrictions that do not propagate along the intestine. ... We investigated the enteric circuits that regulate segmentation focusing on a central feature of the ENS: a recurrent excitatory network of intrinsic sensory neurons (ISNs) which are characterized by prolonged after-hyperpolarizing potentials (AHPs) following their action potentials. ..."
Reference:
1 . Chambers JD, Bornstein JC, Thomas EA (2011) Multiple Neural Oscillators and Muscle Feedback Are Required for the Intestinal Fed State Motor Program PLoS ONE 6(5):e119597-.
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Realistic Network; Neural mass;
Brain Region(s)/Organism:
Cell Type(s):
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB; Simulink;
Model Concept(s): Activity Patterns;
Implementer(s): Chambers, Jordan [jordandchambers at gmail.com];
/
ChamberEtAl2010
README.txt
Chambers et al 2010.mdl
muscleresponse.m
                            
function [A Q] = muscleresponse(l)
    frs = 0;
    frf = 0;
    fds = 0;
    fdf = 0;
    fneg = 0;
    for i = 1:1:length(l)
        if l(i,1) >= 200
            if fneg == 0
                if l(i,2) < 0
                    fneg = 1;
                end
            elseif fneg == 1
                if frs == 0
                    if l(i,2) > 0
                        frs = l(i,1);
                    end
                else
                    if l(i,2) > 0
                        frf = l(i,1);
                    else
                        fds = l(i,1);
                        fneg = 2;
                    end
                end
            elseif fneg == 2
                if l(i,2) < 0
                    fdf = l(i,1);
                else
                    fneg = 3;
                end
            end
        end
    end
    sprintf('Activity = %g (from %g to %g)\nQuiescence = %g (form %g to %g)', (frf - frs), frs, frf, (fdf - fds), fds, fdf)
	if nargout>0
		A = frf - frs;
		Q = fdf - fds;
	end
end