Activity dependent changes in dendritic spine density and spine structure (Crook et al. 2007)

 Download zip file 
Help downloading and running models
Accession:114342
"... In this work, we extend previous modeling studies [27] by combining a model for activity-dependent spine density with one for calcium-mediated spine stem restructuring. ... Additional equations characterize the change in spine density along the dendrite, the current balance equation for an individual spine head, the change in calcium concentration in the spine head, and the dynamics of spine stem resistance. We use computational studies to investigate the changes in spine density and structure for differing synaptic inputs and demonstrate the effects of these changes on the input-output properties of the dendritic branch. ... "
Reference:
1 . Crook SM, Dur-E-Ahmad M, Baer SM (2007) A model of activity-dependent changes in dendritic spine density and spine structure. Math Biosci Eng 4:617-31 [PubMed]
Citations  Citation Browser
Model Information (Click on a link to find other models with that property)
Model Type: Dendrite;
Brain Region(s)/Organism:
Cell Type(s):
Channel(s): I Na,t; I K;
Gap Junctions:
Receptor(s): AMPA;
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Synaptic Plasticity; Synaptic Integration; Calcium dynamics;
Implementer(s):
Search NeuronDB for information about:  AMPA; I Na,t; I K;
/
mbe
readme.html
dmc.m
fig8a.jpg
fig8c.jpg
gsyn.m
MBEv4n4.pdf
run_sp.m
sbf_sp.m
                            
% This program generates the frequancy stimulation  used in the second
% equation for V_sh in the term I_syn
 
function y=gsyn(x,t)

% Set the period

if(t<=500)
   period = 100;         % 10 Hz
  end
if (t>500 & t<1500) 
   period = 8;           % 125 Hz
end
if (t>=1500)
   period = 100;         % 10 Hz
end

% Other parameters

Dx=0.2;
tp = 0.2;               % time to peak for g_syn (ms)
gp = 0.074*10^(-9);     % Peak synaptic conductance  (S);
tloc=mod(t,period);

% Set the conductance for the region of synaptic input

if (x>=0 & x <= Dx)
   y=gp*(tloc/tp)*exp(1-tloc/tp);
else
   y=0;
end