Model Mathematics

Component based math viewer is available

Component: environment

Component: cell

i_Stim = \{\begin{matrix} - stim_amplitude if time - (⌊ \frac{time}{stim_period} ⌋ stim_period) \geq stim_start \land time - (⌊ \frac{time}{stim_period} ⌋ stim_period) ≦ stim_start + stim_duration \\ 0 otherwise \end{matrix} \frac{d}{d time} (V) = - (i_Na + i_Nab + i_NaK + i_Kr + i_Ks + i_tos + i_tof + i_K1 + i_Kp + i_NaCa + i_Cl_Ca + i_Clb + i_CaL + i_Cab + i_Cap + i_Stim)

Component: model_parameters

Cm = \frac{\frac{Cm_per_area 2 cell_radius}{10000} π cell_length}{10000} Vol_Cell = \frac{3.141592654 {(\frac{cell_radius}{1000})}^{2} cell_length}{1000^{3}} Vol_SR = 0.035 Vol_Cell Vol_SL = 0.02 Vol_Cell Vol_jct = 0.00051 Vol_Cell Vol_cytosol = 0.65 Vol_Cell Vol_mito = Vol_Cell - (Vol_SR + Vol_SL + Vol_jct + Vol_cytosol)

Component: reversal_potentials

E_Na_jct = \frac{R T}{F} \ln (\frac{Nao}{Na_jct}) E_Na_SL = \frac{R T}{F} \ln (\frac{Nao}{Na_SL}) E_Ca_jct = \frac{R T}{2 F} \ln (\frac{Cao}{Ca_jct}) E_Ca_SL = \frac{R T}{2 F} \ln (\frac{Cao}{Ca_SL}) E_K = \frac{R T}{F} \ln (\frac{Ko}{Ki}) E_Cl = \frac{(- R) T}{F} \ln (\frac{Clo}{Cli})

Component: INa

openProb = m^{3} h j i_Na_jct = Fx_Na_jct G_INa openProb (V - E_Na_jct) i_Na_SL = Fx_Na_SL G_INa openProb (V - E_Na_SL) i_Na = i_Na_jct + i_Na_SL

Component: INa_h_gate

alpha_h = \{\begin{matrix} 0.135 ⅇ^{\frac{80 + V}{- 6.8}} if V < - 40 \\ 0 otherwise \end{matrix} beta_h = \{\begin{matrix} 3.56 ⅇ^{0.079 V} + 3.1e5 ⅇ^{0.35 V} if V < - 40 \\ \frac{1}{0.13 (1 + ⅇ^{\frac{V + 10.66}{- 11.1}})} otherwise \end{matrix} tau_h = \frac{1}{alpha_h + beta_h} h_infinity = \frac{alpha_h}{alpha_h + beta_h} \frac{d}{d time} (h) = \frac{h_infinity - h}{tau_h}

Component: INa_j_gate

alpha_j = \{\begin{matrix} \frac{\frac{((- 1.2714e5) ⅇ^{0.2444 V} - ((3.474e-5) ⅇ^{(- 0.04391) V})) (V + 37.78)}{1}}{1 + ⅇ^{0.311 (V + 79.23)}} if V < - 40 \\ 0 otherwise \end{matrix} beta_j = \{\begin{matrix} \frac{0.1212 ⅇ^{(- 0.01052) V}}{1 + ⅇ^{(- 0.1378) (V + 40.14)}} if V < - 40 \\ \frac{0.3 ⅇ^{(- (2.535e-7)) V}}{1 + ⅇ^{(- 0.1) (V + 32)}} otherwise \end{matrix} tau_j = \frac{1}{alpha_j + beta_j} j_infinity = \frac{alpha_j}{alpha_j + beta_j} \frac{d}{d time} (j) = \frac{j_infinity - j}{tau_j}

Component: INa_m_gate

alpha_m = \frac{\frac{0.32 (V + 47.13)}{1}}{1 - (ⅇ^{(- 0.1) (V + 47.13)})} beta_m = 0.08 ⅇ^{\frac{- V}{11}} tau_m = \frac{1}{alpha_m + beta_m} m_infinity = \frac{alpha_m}{alpha_m + beta_m} \frac{d}{d time} (m) = \frac{m_infinity - m}{tau_m}

Component: INab

i_Nab_jct = Fx_NaBk_jct G_NaBk (V - E_Na_jct) i_Nab_SL = Fx_NaBk_SL G_NaBk (V - E_Na_SL) i_Nab = i_Nab_jct + i_Nab_SL

Component: INaK

Q_NaK = {Q10_NaK}^{\frac{T - 310}{10}} Q_Km_Nai = {Q10_Km_Nai}^{\frac{T - 310}{10}} sigma = \frac{ⅇ^{\frac{Nao}{67.3}} - 1}{7} f_NaK = \frac{1}{1 + 0.1245 ⅇ^{\frac{(- 0.1) V F}{R T}} + 0.0365 sigma ⅇ^{\frac{(- V) F}{R T}}} i_NaK_jct = \frac{\frac{Fx_NaK_jct Q_NaK I_NaK_max f_NaK}{1 + {(\frac{Q_Km_Nai Km_Nai}{Na_jct})}^{H_NaK}} Ko}{Ko + Km_Ko} i_NaK_SL = \frac{\frac{Fx_NaK_SL Q_NaK I_NaK_max f_NaK}{1 + {(\frac{Q_Km_Nai Km_Nai}{Na_SL})}^{H_NaK}} Ko}{Ko + Km_Ko} i_NaK = i_NaK_jct + i_NaK_SL

Component: IKr

G_IKr = 0.03 \sqrt{\frac{Ko}{5.4}} i_Kr = G_IKr Xr Rr (V - E_K)

Component: IKr_Xr_gate

\frac{d}{d time} (Xr) = \frac{Xr_infinity - Xr}{tau_Xr} Xr_infinity = \frac{1}{1 + ⅇ^{\frac{- (50 + V)}{7.5}}} tau_Xr = \frac{1}{\frac{0.00138 (V + 7)}{1 - (ⅇ^{(- 0.123) (V + 7)})} + \frac{0.00061 (V + 10)}{ⅇ^{0.145 (V + 10)} - 1}}

Component: IKr_Rr_gate

Rr = \frac{1}{1 + ⅇ^{\frac{33 + V}{22.4}}}

Component: IKs

pCa_jct = - (\log_{10} (\frac{Ca_jct}{1} + 3)) pCa_SL = - (\log_{10} (\frac{Ca_SL}{1} + 3)) G_Ks_jct = 0.07 (0.057 + \frac{0.19}{1 + ⅇ^{(- 7.2) + \frac{pCa_jct}{0.6}}}) G_Ks_SL = 0.07 (0.057 + \frac{0.19}{1 + ⅇ^{(- 7.2) + \frac{pCa_SL}{0.6}}}) E_Ks_jct = \frac{R T}{F} \ln (\frac{Ko + pKNa Nao}{Ki + pKNa Na_jct}) E_Ks_SL = \frac{R T}{F} \ln (\frac{Ko + pKNa Nao}{Ki + pKNa Na_SL}) E_Ks = \frac{R T}{F} \ln (\frac{Ko + pKNa Nao}{Ki + pKNa Nai}) i_Ks_jct = Fx_Ks_jct G_Ks_jct {Xs}^{2} (V - E_Ks) i_Ks_SL = Fx_Ks_SL G_Ks_SL {Xs}^{2} (V - E_Ks) i_Ks = i_Ks_jct + i_Ks_SL

Component: IKs_Xs_gate

\frac{d}{d time} (Xs) = \frac{Xs_infinity - Xs}{tau_Xs} Xs_infinity = \frac{1}{1 + ⅇ^{\frac{- (V - 1.5)}{16.7}}} tau_Xs = \frac{1}{\frac{(7.19e-5) (V + 30)}{1 - (ⅇ^{(- 0.148) (V + 30)})} + \frac{(1.31e-4) (V + 30)}{(- 1) + ⅇ^{0.0687 (V + 30)}}}

Component: Itos

i_tos = G_tos X_tos (Y_tos + 0.5 R_tos) (V - E_K)

Component: Itos_X_gate

\frac{d}{d time} (X_tos) = \frac{X_tos_infinity - X_tos}{tau_X_tos} X_tos_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 3)}{15}}} tau_X_tos = \frac{9}{1 + ⅇ^{\frac{V + 3}{15}}} + 0.5

Component: Itos_Y_gate

\frac{d}{d time} (Y_tos) = \frac{Y_tos_infinity - Y_tos}{tau_Y_tos} Y_tos_infinity = \frac{1}{1 + ⅇ^{\frac{V + 33.5}{10}}} tau_Y_tos = \frac{3000}{1 + ⅇ^{\frac{V + 60}{10}}} + 30

Component: Itos_R_gate

R_tos_infinity = \frac{1}{1 + ⅇ^{\frac{V + 33.5}{10}}} tau_R_tos = \frac{2.8e3}{1 + ⅇ^{\frac{V + 60}{10}}} + 220 R_tos = R_tos_infinity \frac{d}{d time} (R_tos_other) = \frac{R_tos_infinity - R_tos}{tau_R_tos}

Component: Itof

i_tof = G_tof X_tof Y_tof (V - E_K)

Component: Itof_X_gate

\frac{d}{d time} (X_tof) = \frac{X_tof_infinity - X_tof}{tau_X_tof} X_tof_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 3)}{15}}} tau_X_tof = 3.5 ⅇ^{- ({(\frac{V}{30})}^{2})} + 1.5

Component: Itof_Y_gate

\frac{d}{d time} (Y_tof) = \frac{Y_tof_infinity - Y_tof}{tau_Y_tof} Y_tof_infinity = \frac{1}{1 + ⅇ^{\frac{V + 33.5}{10}}} tau_Y_tof = \frac{20}{1 + ⅇ^{\frac{V + 33.5}{10}}} + 20

Component: IK1

G_K1 = 0.9 \sqrt{\frac{Ko}{5.4}} i_K1 = G_K1 K1_infinity (V - E_K)

Component: IK1_K1_gate

K1_infinity = \frac{alpha_K1}{alpha_K1 + beta_K1} alpha_K1 = \frac{1.02}{1 + ⅇ^{0.2385 (V - (E_K + 59.215))}} beta_K1 = \frac{0.49124 ⅇ^{0.08032 ((V - E_K) + 5.476)} + 1 ⅇ^{0.06175 (V - (E_K + 594.31))}}{1 + ⅇ^{(- 0.5143) ((V - E_K) + 4.753)}}

Component: IKp

i_Kp = \frac{g_Kp (V - E_K)}{1 + ⅇ^{7.488 - (\frac{V}{5.98})}}

Component: ICl_Ca

i_Cl_Ca = G_Cl (V - E_Cl) (\frac{Fx_Cl_jct}{1 + \frac{Kd_ClCa}{Ca_jct}} + \frac{Fx_Cl_SL}{1 + \frac{Kd_ClCa}{Ca_SL}})

Component: IClb

i_Clb = G_ClBk (V - E_Cl)

Component: ICaL

Q_CaL = {Q10_CaL}^{\frac{T - 310}{10}} temp = \frac{d f Q_CaL V F^{2}}{R T} i_CaL_Ca_jct = \frac{temp fCa_jct Fx_ICaL_jct PCa 4 (gamma_Cai Ca_jct ⅇ^{\frac{2 V F}{R T}} - (gamma_Cao Cao))}{ⅇ^{\frac{2 V F}{R T}} - 1} i_CaL_Na_jct = \frac{temp fCa_jct Fx_ICaL_jct PNa (gamma_Nai Na_jct ⅇ^{\frac{V F}{R T}} - (gamma_Nao Nao))}{ⅇ^{\frac{V F}{R T}} - 1} i_CaL_Ca_SL = \frac{temp fCa_SL Fx_ICaL_SL PCa 4 (gamma_Cai Ca_SL ⅇ^{\frac{2 V F}{R T}} - (gamma_Cao Cao))}{ⅇ^{\frac{2 V F}{R T}} - 1} i_CaL_Na_SL = \frac{temp fCa_SL Fx_ICaL_SL PNa (gamma_Nai Na_SL ⅇ^{\frac{V F}{R T}} - (gamma_Nao Nao))}{ⅇ^{\frac{V F}{R T}} - 1} i_CaL_K = \frac{temp (fCa_SL Fx_ICaL_SL + fCa_jct Fx_ICaL_jct) PK (gamma_Ki Ki ⅇ^{\frac{V F}{R T}} - (gamma_Ko Ko))}{ⅇ^{\frac{V F}{R T}} - 1} i_CaL = i_CaL_Ca_SL + i_CaL_Ca_jct + i_CaL_Na_SL + i_CaL_Na_jct + i_CaL_K

Component: ICaL_d_gate

d_infinity = \frac{1}{1 + ⅇ^{\frac{- (V + 14.5)}{6}}} tau_d = \frac{1 d_infinity (1 - (ⅇ^{\frac{- (V + 14.5)}{6}}))}{0.035 (V + 14.5)} \frac{d}{d time} (d) = \frac{d_infinity - d}{tau_d}

Component: ICaL_f_gate

f_infinity = \frac{1}{1 + ⅇ^{\frac{V + 35.06}{3.6}}} + \frac{0.6}{1 + ⅇ^{\frac{50 - V}{20}}} tau_f = \frac{1}{0.0197 ⅇ^{- ({(0.0337 (V + 14.5))}^{2})} + 0.02} \frac{d}{d time} (f) = \frac{f_infinity - f}{tau_f}

Component: ICaL_fCa_gate

fCa_SL = 1 - fCaB_SL fCa_jct = 1 - fCaB_jct \frac{d}{d time} (fCaB_SL) = 1.7 Ca_SL (1 - fCaB_SL) - (11.9 fCaB_SL) \frac{d}{d time} (fCaB_jct) = 1.7 Ca_jct (1 - fCaB_jct) - (11.9 fCaB_jct)

Component: INaCa

temp_jct = \frac{ⅇ^{\frac{eta V F}{R T}} {Na_jct}^{HNa} Cao - (ⅇ^{\frac{(eta - 1) V F}{R T}} {Nao}^{HNa} Ca_jct)}{1 + ksat ⅇ^{\frac{(eta - 1) V F}{R T}}} temp_SL = \frac{ⅇ^{\frac{eta V F}{R T}} {Na_SL}^{HNa} Cao - (ⅇ^{\frac{(eta - 1) V F}{R T}} {Nao}^{HNa} Ca_SL)}{1 + ksat ⅇ^{\frac{(eta - 1) V F}{R T}}} Q_NCX = {Q10_NCX}^{\frac{T - 310}{10}} Ka_SL = \frac{1}{1 + {(\frac{Kd_act}{Ca_SL})}^{3}} Ka_jct = \frac{1}{1 + {(\frac{Kd_act}{Ca_jct})}^{3}} i_NaCa_jct = \frac{Fx_NCX_jct V_max Ka_jct Q_NCX temp_jct}{K_mCai {Nao}^{HNa} (1 + {(\frac{Na_jct}{K_mNai})}^{HNa}) + {K_mNao}^{HNa} Ca_jct (1 + \frac{Ca_jct}{K_mCai}) + K_mCao {Na_jct}^{HNa} + {Na_jct}^{HNa} Cao + {Nao}^{HNa} Ca_jct} i_NaCa_SL = \frac{Fx_NCX_SL V_max Ka_SL Q_NCX temp_SL}{K_mCai {Nao}^{HNa} (1 + {(\frac{Na_SL}{K_mNai})}^{HNa}) + {K_mNao}^{HNa} Ca_SL (1 + \frac{Ca_SL}{K_mCai}) + K_mCao {Na_SL}^{HNa} + {Na_SL}^{HNa} Cao + {Nao}^{HNa} Ca_SL} i_NaCa = i_NaCa_jct + i_NaCa_SL

Component: ICap

V_maxAF_old = \frac{V_max Vol_cytosol 2 F}{Cm} Q_SLCaP = {Q10_SLCaP}^{\frac{T - 310}{10}} i_Cap_jct = \frac{Q_SLCaP V_maxAF Fx_SLCaP_jct}{1 + {(\frac{Km}{Ca_jct})}^{H}} i_Cap_SL = \frac{Q_SLCaP V_maxAF Fx_SLCaP_SL}{1 + {(\frac{Km}{Ca_SL})}^{H}} i_Cap = i_Cap_jct + i_Cap_SL

Component: ICab

i_Cab_jct = G_CaBk Fx_CaBk_jct (V - E_Ca_jct) i_Cab_SL = G_CaBk Fx_CaBk_SL (V - E_Ca_SL) i_Cab = i_Cab_SL + i_Cab_jct

Component: Jrel_SR

kCaSR = Max_SR - (\frac{Max_SR - Min_SR}{1 + {(\frac{EC50_SR}{Ca_SR})}^{HSR}}) koSRCa = \frac{koCa}{kCaSR} kiSRCa = kiCa kCaSR RI = 1 - R - O - I \frac{d}{d time} (R) = kim RI - (kiSRCa Ca_jct R) - (koSRCa {Ca_jct}^{2} R - (kom O)) \frac{d}{d time} (O) = koSRCa {Ca_jct}^{2} R - (kom O) - (kiSRCa Ca_jct O - (kim I)) \frac{d}{d time} (I) = kiSRCa Ca_jct O - (kim I) - (kom I - (koSRCa {Ca_jct}^{2} RI)) j_rel_SR = ks O (Ca_SR - Ca_jct)

Component: Jleak_SR

j_leak_SR = KSRleak (Ca_SR - Ca_jct)

Component: Jpump_SR

Q_SRCaP = {Q10_SRCaP}^{\frac{T - 310}{10}} j_pump_SR = \frac{\frac{Q_SRCaP V_max Vol_cytosol}{Vol_SR} ({(\frac{Cai}{Kmf})}^{H} - ({(\frac{Ca_SR}{Kmr})}^{H}))}{1 + {(\frac{Cai}{Kmf})}^{H} + {(\frac{Ca_SR}{Kmr})}^{H}}

Component: ion_diffusion

j_Na_jct_SL = \frac{D_Na_jct_SL A_jct_SL}{x_jct_SL} j_Na_SL_cytosol = \frac{D_Na_SL_cytosol A_SL_cytosol}{x_SL_cytosol} j_Ca_jct_SL = \frac{D_Ca_jct_SL A_jct_SL}{x_jct_SL} j_Ca_SL_cytosol = \frac{D_Ca_SL_cytosol A_SL_cytosol}{x_SL_cytosol} j_Na_jct_SL_new = 1.8313e-14 j_Na_SL_cytosol_new = 1.6386e-12 j_Ca_jct_SL_new = 8.2413e-13 j_Ca_SL_cytosol_new = 3.7243e-12 J_Na_jct_SL = (Na_jct - Na_SL) j_Na_jct_SL J_Na_SL_cytosol = (Na_SL - Nai) j_Na_SL_cytosol J_Ca_jct_SL = (Ca_jct - Ca_SL) j_Ca_jct_SL J_Ca_SL_cytosol = (Ca_SL - Cai) j_Ca_SL_cytosol

Component: Na_buffer

dNa_jct_buf = kon Na_jct (Bmax_jct - Na_jct_buf) - (koff Na_jct_buf) dNa_SL_buf = kon Na_SL (Bmax_SL - Na_SL_buf) - (koff Na_SL_buf) \frac{d}{d time} (Na_jct_buf) = dNa_jct_buf \frac{d}{d time} (Na_SL_buf) = dNa_SL_buf \frac{d}{d time} (Na_jct) = \frac{(- Cm) (i_Na_jct + 3 i_NaCa_jct + i_Nab_jct + 3 i_NaK_jct + i_CaL_Na_jct)}{Vol_jct F} - (\frac{J_Na_jct_SL}{Vol_jct}) - dNa_jct_buf \frac{d}{d time} (Na_SL) = \frac{(- Cm) (i_Na_SL + 3 i_NaCa_SL + i_Nab_SL + 3 i_NaK_SL + i_CaL_Na_SL)}{Vol_SL F} + \frac{J_Na_jct_SL - J_Na_SL_cytosol}{Vol_SL} - dNa_SL_buf \frac{d}{d time} (Nai) = \frac{J_Na_SL_cytosol}{Vol_cytosol}

Component: Ca_buffer

dCalsequestrin = kon_Calsequestrin Ca_SR (\frac{Bmax_Calsequestrin Vol_cytosol}{Vol_SR} - Ca_Calsequestrin) - (koff_Calsequestrin Ca_Calsequestrin) \frac{d}{d time} (Ca_Calsequestrin) = dCalsequestrin dCa_SLB_SL = kon_SL Ca_SL (\frac{Bmax_SLB_SL Vol_cytosol}{Vol_SL} - Ca_SLB_SL) - (koff_SLB Ca_SLB_SL) dCa_SLB_jct = kon_SL Ca_jct (\frac{Bmax_SLB_jct Vol_cytosol}{Vol_jct} - Ca_SLB_jct) - (koff_SLB Ca_SLB_jct) dCa_SLHigh_SL = kon_SL Ca_SL (\frac{Bmax_SLHigh_SL Vol_cytosol}{Vol_SL} - Ca_SLHigh_SL) - (koff_SLHigh Ca_SLHigh_SL) dCa_SLHigh_jct = kon_SL Ca_jct (\frac{Bmax_SLHigh_jct Vol_cytosol}{Vol_jct} - Ca_SLHigh_jct) - (koff_SLHigh Ca_SLHigh_jct) \frac{d}{d time} (Ca_SLB_SL) = dCa_SLB_SL \frac{d}{d time} (Ca_SLB_jct) = dCa_SLB_jct \frac{d}{d time} (Ca_SLHigh_SL) = dCa_SLHigh_SL \frac{d}{d time} (Ca_SLHigh_jct) = dCa_SLHigh_jct dCa_jct_tot_bound = dCa_SLB_jct + dCa_SLHigh_jct dCa_SL_tot_bound = dCa_SLB_SL + dCa_SLHigh_SL i_Ca_jct_tot = (i_CaL_Ca_jct - (2 i_NaCa_jct)) + i_Cab_jct + i_Cap_jct i_Ca_SL_tot = (i_CaL_Ca_SL - (2 i_NaCa_SL)) + i_Cab_SL + i_Cap_SL \frac{d}{d time} (Ca_SR) = j_pump_SR - (\frac{j_leak_SR Vol_cytosol}{Vol_SR} + j_rel_SR) - dCalsequestrin \frac{d}{d time} (Ca_jct) = (\frac{(- i_Ca_jct_tot) Cm}{Vol_jct 2 F} - (\frac{J_Ca_jct_SL}{Vol_jct})) + \frac{j_rel_SR Vol_SR}{Vol_jct} + \frac{j_leak_SR Vol_cytosol}{Vol_jct} - (1 dCa_jct_tot_bound) \frac{d}{d time} (Ca_SL) = \frac{(- i_Ca_SL_tot) Cm}{Vol_SL 2 F} + \frac{J_Ca_jct_SL - J_Ca_SL_cytosol}{Vol_SL} - (1 dCa_SL_tot_bound) \frac{d}{d time} (Cai) = \frac{(- j_pump_SR) Vol_SR}{Vol_cytosol} + \frac{J_Ca_SL_cytosol}{Vol_cytosol} - (1 dCa_cytosol_tot_bound)

Component: cytosolic_Ca_buffer

dCa_TroponinC = kon_TroponinC Cai (Bmax_TroponinC - Ca_TroponinC) - (koff_TroponinC Ca_TroponinC) dCa_TroponinC_Ca_Mg = kon_TroponinC_Ca_Mg_Ca Cai (Bmax_TroponinC_Ca_Mg_Ca - (Ca_TroponinC_Ca_Mg + Mg_TroponinC_Ca_Mg)) - (koff_TroponinC_Ca_Mg_Ca Ca_TroponinC_Ca_Mg) dMg_TroponinC_Ca_Mg = kon_TroponinC_Ca_Mg_Mg Mgi (Bmax_TroponinC_Ca_Mg_Mg - (Ca_TroponinC_Ca_Mg + Mg_TroponinC_Ca_Mg)) - (koff_TroponinC_Ca_Mg_Mg Mg_TroponinC_Ca_Mg) dCa_Calmodulin = kon_Calmodulin Cai (Bmax_Calmodulin - Ca_Calmodulin) - (koff_Calmodulin Ca_Calmodulin) dCa_Myosin = kon_Myosin_Ca Cai (Bmax_Myosin_Ca - (Ca_Myosin + Mg_Myosin)) - (koff_Myosin_Ca Ca_Myosin) dMg_Myosin = kon_Myosin_Mg Mgi (Bmax_Myosin_Mg - (Ca_Myosin + Mg_Myosin)) - (koff_Myosin_Mg Mg_Myosin) dCa_SRB = kon_SRB Cai (Bmax_SRB - Ca_SRB) - (koff_SRB Ca_SRB) dCa_cytosol_tot_bound = dCa_TroponinC + dCa_TroponinC_Ca_Mg + dMg_TroponinC_Ca_Mg + dCa_Calmodulin + dCa_Myosin + dMg_Myosin + dCa_SRB \frac{d}{d time} (Ca_TroponinC) = dCa_TroponinC \frac{d}{d time} (Ca_TroponinC_Ca_Mg) = dCa_TroponinC_Ca_Mg \frac{d}{d time} (Mg_TroponinC_Ca_Mg) = dMg_TroponinC_Ca_Mg \frac{d}{d time} (Ca_Calmodulin) = dCa_Calmodulin \frac{d}{d time} (Ca_Myosin) = dCa_Myosin \frac{d}{d time} (Mg_Myosin) = dMg_Myosin \frac{d}{d time} (Ca_SRB) = dCa_SRB

Component: indo_fluo_Ca_buffer_not_connected

Indo1Bound = Ca_Indo1_Cytosol + \frac{Ca_Indo1_jct Vol_jct}{Vol_cytosol} + \frac{Ca_Indo1_SL Vol_SL}{Vol_cytosol} Fluo3Bound = Ca_Fluo3_Cytosol + \frac{Ca_Fluo3_jct Vol_jct}{Vol_cytosol} + \frac{Ca_Fluo3_SL Vol_SL}{Vol_cytosol} dCa_Indo1_Cytosol = kon_Indo1 (Indo1 - Indo1Bound) Cai (Bmax_Indo1_Cytosol - Ca_Indo1_Cytosol) - (koff_Indo1 Ca_Indo1_Cytosol) dCa_Indo1_jct = kon_Indo1 (Indo1 - Indo1Bound) Ca_jct (\frac{Bmax_Indo1_jct Vol_cytosol}{Vol_jct} - Ca_Indo1_jct) - (koff_Indo1 Ca_Indo1_jct) dCa_Indo1_SL = kon_Indo1 (Indo1 - Indo1Bound) Ca_SL (\frac{Bmax_Indo1_SL Vol_cytosol}{Vol_SL} - Ca_Indo1_SL) - (koff_Indo1 Ca_Indo1_SL) dCa_Fluo3_Cytosol = kon_Fluo3 (Fluo3 - Fluo3Bound) Cai (Bmax_Fluo3_Cytosol - Ca_Fluo3_Cytosol) - (koff_Fluo3 Ca_Fluo3_Cytosol) dCa_Fluo3_jct = kon_Fluo3 (Fluo3 - Fluo3Bound) Ca_jct (\frac{Bmax_Fluo3_jct Vol_cytosol}{Vol_jct} - Ca_Fluo3_jct) - (koff_Fluo3 Ca_Fluo3_jct) dCa_Fluo3_SL = kon_Fluo3 (Fluo3 - Fluo3Bound) Ca_SL (\frac{Bmax_Fluo3_SL Vol_cytosol}{Vol_SL} - Ca_Fluo3_SL) - (koff_Fluo3 Ca_Fluo3_SL) \frac{d}{d time} (Ca_Indo1_Cytosol) = dCa_Indo1_Cytosol \frac{d}{d time} (Ca_Indo1_SL) = dCa_Indo1_SL \frac{d}{d time} (Ca_Indo1_jct) = dCa_Indo1_jct \frac{d}{d time} (Ca_Fluo3_Cytosol) = dCa_Fluo3_Cytosol \frac{d}{d time} (Ca_Fluo3_SL) = dCa_Fluo3_SL \frac{d}{d time} (Ca_Fluo3_jct) = dCa_Fluo3_jct

Model Mathematics

Component: environment

Component: cell

Component: model_parameters

Component: reversal_potentials

Component: INa

Component: INa_h_gate

Component: INa_j_gate

Component: INa_m_gate

Component: INab

Component: INaK

Component: IKr

Component: IKr_Xr_gate

Component: IKr_Rr_gate

Component: IKs

Component: IKs_Xs_gate

Component: Itos

Component: Itos_X_gate

Component: Itos_Y_gate

Component: Itos_R_gate

Component: Itof

Component: Itof_X_gate

Component: Itof_Y_gate

Component: IK1

Component: IK1_K1_gate

Component: IKp

Component: ICl_Ca

Component: IClb

Component: ICaL

Component: ICaL_d_gate

Component: ICaL_f_gate

Component: ICaL_fCa_gate

Component: INaCa

Component: ICap

Component: ICab

Component: Jrel_SR

Component: Jleak_SR

Component: Jpump_SR

Component: ion_diffusion

Component: Na_buffer

Component: Ca_buffer

Component: cytosolic_Ca_buffer

Component: indo_fluo_Ca_buffer_not_connected

Component: Ions_n_reversal_potentials

Component: transmembrane_currents

Component: intracellular_currents