Spontaneous calcium oscillations in astrocytes (Lavrentovich and Hemkin 2008)

 Download zip file 
Help downloading and running models
Accession:112547
" ... We propose here a mathematical model of how spontaneous Ca2+ oscillations arise in astrocytes. This model uses the calcium-induced calcium release and inositol cross-coupling mechanisms coupled with a receptor-independent method for producing inositol (1,4,5)-trisphosphate as the heart of the model. By computationally mimicking experimental constraints we have found that this model provides results that are qualitatively similar to experiment."
Reference:
1 . Lavrentovich M, Hemkin S (2008) A mathematical model of spontaneous calcium(II) oscillations in astrocytes. J Theor Biol 251:553-60 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Electrogenic pump; Glia;
Brain Region(s)/Organism:
Cell Type(s): Astrocyte;
Channel(s):
Gap Junctions:
Receptor(s): IP3;
Gene(s):
Transmitter(s): Ions;
Simulation Environment: XPP;
Model Concept(s): Oscillations; Simplified Models; Signaling pathways; Calcium dynamics; Calcium waves;
Implementer(s): Wu, Sheng-Nan [snwu at mail.ncku.edu.tw];
Search NeuronDB for information about:  IP3; Ions;
/
Glia-Ca
readme.html
Astro-Ca.pdf
Ca-Oscil-Astrocyte.ode
Glia-Oscill.JPG
                            
% Ca-Oscil-Astrocyte.ode
% Ref: Laventovich and Hemkin, A mathematical model of spontaneous calcium(II)  oscillations 
% in astrocytes.  J Theor Biol 2008; 251(4):553-60.

% Units: microM/sec; microM; sec
% x: cytosolic Ca; y: Ca concencentration in ER; z: IP3 concentration

init x=0.1, y=1.5, z=0.1

parm vm2=15, vm3=40, vin=0.05, vp=0.05, k2=0.1
parm kcaa=0.15, kcat=0.15, kip3=0.1, kp=0.3
parm kdeg=0.08, kout=0.5, kf=0.5, n=2.02, m=2.2

vserca=vm2*(x^2/(x^2+k2^2))
vplc=vp*(x^2/(x^2+kp^2))
vcicr=4*vm3*((kcaa^n)*(x^n)/((x^n+kcaa^n)*(x^n+kcat^n)))*(z^m/(z^m+kip3^m))*(y-x)

dx/dt=vin-kout*x+vcicr-vserca+kf*(y-x)
dy/dt=vserca-vcicr-kf*(y-x)
dz/dt=vplc-kdeg*z

@ meth=Euler, dt=0.01, total=600
@ yp=x, ylo=0, yhi=0.7, xlo=0, xhi=600, bound=10000

done



Loading data, please wait...