Model Mathematics

Component: environment

Component: volume_i

volume_i=volume_i0

Component: volume_er

volume_er=volume_igamma

Component: volumeCa_i

ddtimevolumeCa_i=Jipt+Jryr+Jer-Jserca+Jin-Jpm

Component: Ca_i

Ca_i=volumeCa_ivolume_i

Component: volumeCa_er

ddtimevolumeCa_er=-Jipt+Jryr+Jer-Jserca

Component: Ca_er

Ca_er=volumeCa_ervolume_er

Component: volume_iIP3

ddtimevolume_iIP3=J_ip3P-J_ip3Dvolume_i

Component: IP3

IP3=volume_iIP3volume_i

Component: PIP2

ddtimePIP2=J_ip3RIPX-J_ip3PPIP2_Total

Component: IPX

IPX=1-IP3PIP2_Total-PIP2

Component: J_ip3R

J_ip3R=J_ip3R0

Component: receptor

V_mech=V_mech0iftime>10time<25MechanicalStimulation>00otherwise V_IP3=C_1ATP_eC_2+ATP_e+C_3ATP_e2C_42+ATP_e2 ddtimeATP_e=-V_ATPATP_eK_ATP+ATP_e

Component: J_ip3P

J_ip3P=V_IP3PIP2

Component: J_ip3D

F_Ca=Ca_iK_rc+Ca_i J_ip3D=Beta_1+Beta_2F_Ca2IP3K_IP3+IP3

Component: Jipt

Pipt=O4 Jipt=Jipt0volume_i0PiptCa_er-Ca_i

Component: Jryr

W_inf=1+KaCa_i4+Ca_iKb31+1Kc+Ca_iKb3+KaCa_i4 Pryr=W_inf1+KaCa_i4+Ca_iKb31+1Kc+Ca_iKb3+KaCa_i4 Jryr=Jryr0volume_i0Ca_er-Ca_i

Component: Jer

Jer=Jer0volume_i0Ca_er-Ca_i

Component: Jserca

Jserca=volume_i0VsercaCa_i2Kserca2+Ca_i21Ca_er

Component: Jin

V_Ca=RTz_CaFlnCa_extCa_i Iin_1=volume_i0P_pm_1V_Ca-Vm Iin_2=volume_i0P_pm_2V_mechV_Ca-Vm Jin_1=Iin_1z_CaF Jin_2=Iin_2z_CaF Jin=Jin_1+Jin_2

Component: Jpm

Jpm=Vpmvolume_i0Ca_i2Kpm2+Ca_i2

Component: O

ddtimeO=k1IP3S-k1_O+k2O S=1-O+I1+I2

Component: I1

ddtimeI1=k2O-k3+k4I1

Component: I2

ddtimeI2=k4I1-k5I2

Component: constants

k1=alpha_1Ca_i3beta_13+Ca_i3 k4=alpha_4IP3beta_4+IP3