Continuum model of tubulin-driven neurite elongation (Graham et al 2006)

 Download zip file 
Help downloading and running models
Accession:59581
This model investigates the elongation over time of a single developing neurite (axon or dendrite). Our neurite growth model describes the elongation of a single,unbranched neurite in terms of the rate of extension of the microtubule cytoskeleton. The cytoskeleton is not explicitly modelled, but its construction is assumed to depend on the available free tubulin at the growing neurite tip.
References:
1 . Graham BP, Lauchlan K, McLean DR (2006) Dynamics of outgrowth in a continuum model of neurite elongation. J Comput Neurosci 20:43-60 [PubMed]
2 . Graham BP, van Ooyen A (2006) Mathematical modelling and numerical simulation of the morphological development of neurons. BMC Neurosci 7 Suppl 1:S9 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Axon; Dendrite;
Brain Region(s)/Organism:
Cell Type(s):
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Development; Olfaction;
Implementer(s): Graham, Bruce [B.Graham at cs.stir.ac.uk];
function [simp, modp, calcp] = CMNG_params
% CMNG_parameters Parameter values for Continuum Model 
%   for Autoregulatory Neurite Outgrowth
% simp - simulation parameters
% modp - user-defined (dimensional) parameters
% calcp - calculated (non-dimensional) parameters
% Version 2.0 (BPG & DRM 24-1-03)

% Parameters

% simulation
simp.dt = 0.1;                  % time step
simp.tmax = 1500;                  % simulation time
simp.datat = 10;                 % data collection time step
simp.N = 100;                     % number of spatial points
simp.kmax = 10000;               % maximum corrector steps
simp.mc = 0.0001;              % tolerance on C;
simp.ml = 0.0001;              % tolerance on l;

% user-defined
modp.c0 = 10;                     % concentration scale
modp.c1 = 10;                       % autoregulatory ideal concentration
modp.l0 = 0.01;                   % initial (min) length;
modp.D = 30000;                    % diffusion constant
modp.a = 100;                      % active transport rate
modp.g = 0.002;                    % decay rate
% normally rg=el & sg=zl
modp.rg = 10;                   % growth rate constant
modp.sg = 100;                   % growth rate set point (threshold)
modp.e0 = 2e-5;                  % soma flux-source rate
modp.el = 1e-4;                     % growth tip flux-sink rate
modp.rg = modp.el;                   % growth cone flux-sink rate
modp.zl = modp.sg;                   % growth cone flux-source rate

% calculated
calcp.alpha     = modp.a/(modp.rg*modp.c0)
calcp.beta      = (modp.g*modp.D)/(modp.rg*modp.rg*modp.c0*modp.c0)
calcp.gamma     = modp.sg/(modp.rg*modp.c0)
calcp.phi       = (modp.e0*modp.D)/(modp.rg*modp.c0)
calcp.rho       = modp.el/modp.e0
calcp.theta     = modp.c0/modp.c1
calcp.sigma     = (modp.el*modp.sg)/(modp.e0*modp.rg*modp.c0)
calcp.delaystep = 0                 % no delay of delaystep*dt in x=0 bdy cn

calcp.jmax      = round(simp.tmax/simp.dt)      % final time step index
calcp.dy        = 1/simp.N;                   % nondimensional spatial step

Loading data, please wait...