Model Mathematics

Component based math viewer is available

Component: Istim_parameters

Component: outputs

Component: g_parameters

Component: Membrane_potential

\frac{d}{d time} (V) = \frac{- (I_stim + I_ion)}{Cm} I_ion = I_CaL + I_CaT + I_Kv + I_BK + I_Na + I_NCX + I_NaK + I_NsNa + I_NsK

Component: Cai

Cai_total_rate = \frac{(- (I_CaL + I_CaT - (2 I_NCX))) Cvt}{2 F V_cell} \frac{d}{d time} (Cai_total) = \frac{(- (I_CaL + I_CaT - (2 I_NCX))) Cvt}{2 F V_cell} \frac{d}{d time} (Cai) = \frac{Cai_total_rate}{1 + \frac{n_CRT CRT_total K_D_CRT {Cai}^{n_CRT - 1}}{{({Cai}^{n_CRT} + K_D_CRT)}^{2}} + \frac{n_CaM CaM_total K_D_CaM {Cai}^{n_CaM - 1}}{{({Cai}^{n_CaM} + K_D_CaM)}^{2}}}

Component: Nai

\frac{d}{d time} (Nai) = \frac{(- (I_Na + 3 I_NaK + 3 I_NCX + I_NsNa)) Cvt}{F V_cell}

Component: Ki

\frac{d}{d time} (Ki) = \frac{(- ((I_Kv + I_BK + I_stim - (2 I_NaK)) + I_NsK)) Cvt}{F V_cell}

Component: Cai

Cai_total_rate = \frac{(- (I_CaL + I_CaT - (2 I_NCX))) Cvt}{2 F V_cell} \frac{d}{d time} (Cai_total) = \frac{(- (I_CaL + I_CaT - (2 I_NCX))) Cvt}{2 F V_cell} \frac{d}{d time} (Cai) = \frac{Cai_total_rate}{1 + \frac{n_CRT CRT_total K_D_CRT {Cai}^{n_CRT - 1}}{{({Cai}^{n_CRT} + K_D_CRT)}^{2}} + \frac{n_CaM CaM_total K_D_CaM {Cai}^{n_CaM - 1}}{{({Cai}^{n_CaM} + K_D_CaM)}^{2}}}

Component: Nai

\frac{d}{d time} (Nai) = \frac{(- (I_Na + 3 I_NaK + 3 I_NCX + I_NsNa)) Cvt}{F V_cell}

Component: Ki

\frac{d}{d time} (Ki) = \frac{(- ((I_Kv + I_BK + I_stim - (2 I_NaK)) + I_NsK)) Cvt}{F V_cell}

Component: Cai

Cai_total_rate = \frac{(- (I_CaL + I_CaT - (2 I_NCX))) Cvt}{2 F V_cell} \frac{d}{d time} (Cai_total) = \frac{(- (I_CaL + I_CaT - (2 I_NCX))) Cvt}{2 F V_cell} \frac{d}{d time} (Cai) = \frac{Cai_total_rate}{1 + \frac{n_CRT CRT_total K_D_CRT {Cai}^{n_CRT - 1}}{{({Cai}^{n_CRT} + K_D_CRT)}^{2}} + \frac{n_CaM CaM_total K_D_CaM {Cai}^{n_CaM - 1}}{{({Cai}^{n_CaM} + K_D_CaM)}^{2}}}

Component: Nai

\frac{d}{d time} (Nai) = \frac{(- (I_Na + 3 I_NaK + 3 I_NCX + I_NsNa)) Cvt}{F V_cell}

Component: Ki

\frac{d}{d time} (Ki) = \frac{(- ((I_Kv + I_BK + I_stim - (2 I_NaK)) + I_NsK)) Cvt}{F V_cell}

Component: ICaL

Component: Ionic_currents

I = g_max PO (V - E)

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: ICaL_channel_states

alpha = phi 0.731 ⅇ^{\frac{V}{30}} beta = phi 0.2149 ⅇ^{\frac{- V}{40}} alpha_0 = 4 alpha alpha_1 = 3 alpha alpha_2 = 2 alpha alpha_3 = 1 alpha beta_0 = 1 beta beta_1 = 2 beta beta_2 = 3 beta beta_3 = 4 beta phi_f = phi 0.4742 ⅇ^{\frac{V}{10}} phi_s = phi 0.005956 ⅇ^{\frac{- V}{40}} xi_f = phi 0.01407 ⅇ^{\frac{- V}{300}} xi_s = phi 0.01213 ⅇ^{\frac{V}{500}} psi_f = phi 0.02197 ⅇ^{\frac{V}{500}} psi_s = phi 0.00232 ⅇ^{\frac{- V}{280}} omega_f = \frac{beta_3 xi_f phi_f}{alpha_3 psi_f} omega_s = \frac{beta_3 xi_s phi_s}{alpha_3 psi_s} omega_sf = \frac{xi_s psi_f}{xi_f} omega_fs = psi_s delta = phi 0.01 theta = \frac{phi 4}{1 + \frac{1}{Cai}} norm = C3 + C2 + C1 + C0 + C3Ca + C2Ca + C1Ca + C0Ca + O_CaL + ICa + IVs + IVf + IVsCa + IVfCa \frac{d}{d time} (C3) = (\frac{alpha_2 C2}{norm} - (\frac{(alpha_3 + beta_2 + phi_s + phi_f + theta) C3}{norm})) + \frac{beta_3 O_CaL}{norm} + \frac{omega_f IVf}{norm} + \frac{omega_s IVs}{norm} + \frac{delta C3Ca}{norm} \frac{d}{d time} (C2) = (\frac{alpha_1 C1}{norm} - (\frac{(alpha_2 + beta_1 + theta) C2}{norm})) + \frac{beta_2 C3}{norm} + \frac{delta C2Ca}{norm} \frac{d}{d time} (C1) = (\frac{alpha_0 C0}{norm} - (\frac{(alpha_1 + beta_0 + theta) C1}{norm})) + \frac{beta_1 C2}{norm} + \frac{delta C1Ca}{norm} \frac{d}{d time} (C0) = \frac{(- (alpha_0 + theta)) C0}{norm} + \frac{beta_0 C1}{norm} + \frac{delta C0Ca}{norm} \frac{d}{d time} (C3Ca) = (\frac{theta C3}{norm} + \frac{alpha_2 C2Ca}{norm} - (\frac{(alpha_3 + beta_2 + phi_f + phi_s + delta) C3Ca}{norm})) + \frac{beta_3 ICa}{norm} + \frac{omega_f IVfCa}{norm} + \frac{omega_s IVsCa}{norm} \frac{d}{d time} (C2Ca) = (\frac{theta C2}{norm} + \frac{alpha_1 C1Ca}{norm} - (\frac{(alpha_2 + beta_1 + delta) C2Ca}{norm})) + \frac{beta_2 C3Ca}{norm} \frac{d}{d time} (C1Ca) = (\frac{theta C1}{norm} + \frac{alpha_0 C0Ca}{norm} - (\frac{(alpha_1 + beta_0 + delta) C1Ca}{norm})) + \frac{beta_1 C2Ca}{norm} \frac{d}{d time} (C0Ca) = (\frac{theta C0}{norm} - (\frac{(alpha_0 + delta) C0Ca}{norm})) + \frac{beta_0 C1Ca}{norm} \frac{d}{d time} (O_CaL) = (\frac{alpha_3 C3}{norm} - (\frac{(beta_3 + psi_f + psi_s + theta) O_CaL}{norm})) + \frac{xi_f IVf}{norm} + \frac{xi_s IVs}{norm} + \frac{delta ICa}{norm} \frac{d}{d time} (ICa) = (\frac{theta O_CaL}{norm} + \frac{alpha_3 C3Ca}{norm} - (\frac{(beta_3 + psi_f + psi_s + delta) ICa}{norm})) + \frac{xi_f IVfCa}{norm} + \frac{xi_s IVsCa}{norm} \frac{d}{d time} (IVs) = (\frac{phi_s C3}{norm} + \frac{psi_s O_CaL}{norm} + \frac{omega_fs IVf}{norm} - (\frac{(omega_s + xi_s + omega_sf + theta) IVs}{norm})) + \frac{delta IVsCa}{norm} \frac{d}{d time} (IVf) = (\frac{phi_f C3}{norm} + \frac{psi_f O_CaL}{norm} - (\frac{(omega_f + xi_f + omega_fs + theta) IVf}{norm})) + \frac{omega_sf IVs}{norm} + \frac{delta IVfCa}{norm} \frac{d}{d time} (IVsCa) = \frac{theta IVs}{norm} + \frac{phi_s C3Ca}{norm} + \frac{psi_s ICa}{norm} + \frac{omega_fs IVfCa}{norm} - (\frac{(omega_s + xi_s + omega_sf + delta) IVsCa}{norm}) \frac{d}{d time} (IVfCa) = (\frac{theta IVf}{norm} + \frac{phi_f C3Ca}{norm} + \frac{psi_f ICa}{norm} - (\frac{(omega_f + xi_f + omega_fs + delta) IVfCa}{norm})) + \frac{omega_sf IVsCa}{norm}

Component: temperature_factor

phi = {Q10}^{\frac{T - T0}{10}}

Component: ICaT

P_CaT = d_CaT f_CaT

Component: Ionic_currents

I = g_max PO (V - E)

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: gating_kinetics

\frac{d}{d time} (X) = \frac{X_inf - X}{tau_X} phi

Component: temperature_factor

phi = {Q10}^{\frac{T - T0}{10}}

Component: gating_kinetics

\frac{d}{d time} (X) = \frac{X_inf - X}{tau_X} phi

Component: temperature_factor

phi = {Q10}^{\frac{T - T0}{10}}

Component: d_inf

X_inf = \frac{1}{1 + ⅇ^{\frac{- (V + 60.5)}{5.3}}}

Component: f_inf

X_inf = \frac{1}{1 + ⅇ^{\frac{V + 75.5}{4.0}}}

Component: tau_f

tau = 0.38117 (8.6 + 14.7 ⅇ^{\frac{- ({(V + 50)}^{2})}{900 1}})

Component: IKv

P_Kv = x_Kv y_Kv

Component: Ionic_currents

I = g_max PO (V - E)

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: gating_kinetics

\frac{d}{d time} (X) = \frac{X_inf - X}{tau_X} phi

Component: temperature_factor

phi = {Q10}^{\frac{T - T0}{10}}

Component: gating_kinetics

\frac{d}{d time} (X) = \frac{X_inf - X}{tau_X} phi

Component: temperature_factor

phi = {Q10}^{\frac{T - T0}{10}}

Component: x_inf

X_inf = \frac{1}{1 + ⅇ^{\frac{- (V + 43.0)}{17.36}}}

Component: y_inf

X_inf = \frac{1}{1 + ⅇ^{\frac{V - 44.9}{12.0096}}}

Component: IBK

Component: Ionic_currents

I = g_max PO (V - E)

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: IBK_channel_states

alpha = 1 ⅇ^{\frac{8.47188 V}{1 T}} beta = 1 ⅇ^{\frac{(- 7.77556) V}{1 T}} k_C0O0 = 0.02162 alpha k_C1O1 = 0.000869 alpha k_C2O2 = 0.0000281 alpha k_C3O3 = 0.000781 alpha k_C4O4 = 0.044324 alpha k_O0C0 = 318.1084 beta k_O1C1 = 144.1736 beta k_O2C2 = 32.6594 beta k_O3C3 = 0.095312 beta k_O4C4 = 0.000106 beta k_C0C1 = 4 k_on Cai k_C1C2 = 3 k_on Cai k_C2C3 = 2 k_on Cai k_C3C4 = 1 k_on Cai k_C4C3 = 4 k_off_C k_C3C2 = 3 k_off_C k_C2C1 = 2 k_off_C k_C1C0 = 1 k_off_C k_O0O1 = 4 k_on Cai k_O1O2 = 3 k_on Cai k_O2O3 = 2 k_on Cai k_O3O4 = 1 k_on Cai k_O4O3 = 4 k_off_O k_O3O2 = 3 k_off_O k_O2O1 = 2 k_off_O k_O1O0 = 1 k_off_O norm = C0 + C1 + C2 + C3 + C4 + O0 + O1 + O2 + O3 + O4 \frac{d}{d time} (C4) = \frac{(- (k_C4C3 + k_C4O4)) C4}{norm} + \frac{k_C3C4 C3}{norm} + \frac{k_O4C4 O4}{norm} \frac{d}{d time} (C3) = \frac{(- (k_C3C2 + k_C3O3 + k_C3C4)) C3}{norm} + \frac{k_C2C3 C2}{norm} + \frac{k_O3C3 O3}{norm} + \frac{k_C4C3 C4}{norm} \frac{d}{d time} (C2) = \frac{(- (k_C2C1 + k_C2O2 + k_C2C3)) C2}{norm} + \frac{k_C1C2 C1}{norm} + \frac{k_O2C2 O2}{norm} + \frac{k_C3C2 C3}{norm} \frac{d}{d time} (C1) = \frac{(- (k_C1C0 + k_C1O1 + k_C1C2)) C1}{norm} + \frac{k_C0C1 C0}{norm} + \frac{k_O1C1 O1}{norm} + \frac{k_C2C1 C2}{norm} \frac{d}{d time} (C0) = \frac{(- (k_C0C1 + k_C0O0)) C0}{norm} + \frac{k_C1C0 C1}{norm} + \frac{k_O0C0 O0}{norm} \frac{d}{d time} (O4) = \frac{(- (k_O4O3 + k_O4C4)) O4}{norm} + \frac{k_O3O4 O3}{norm} + \frac{k_C4O4 C4}{norm} \frac{d}{d time} (O3) = \frac{(- (k_O3O2 + k_O3C3 + k_O3O4)) O3}{norm} + \frac{k_O2O3 O2}{norm} + \frac{k_C3O3 C3}{norm} + \frac{k_O4O3 O4}{norm} \frac{d}{d time} (O2) = \frac{(- (k_O2O1 + k_O2C2 + k_O2O3)) O2}{norm} + \frac{k_O1O2 O1}{norm} + \frac{k_C2O2 C2}{norm} + \frac{k_O3O2 O3}{norm} \frac{d}{d time} (O1) = \frac{(- (k_O1O0 + k_O1C1 + k_O1O2)) O1}{norm} + \frac{k_O0O1 O0}{norm} + \frac{k_C1O1 C1}{norm} + \frac{k_O2O1 O2}{norm} \frac{d}{d time} (O0) = \frac{(- (k_O0O1 + k_O0C0)) O0}{norm} + \frac{k_O1O0 O1}{norm} + \frac{k_C0O0 C0}{norm}

Component: INa

Component: Ionic_currents

I = g_max PO (V - E)

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: INa_channel_states

k_I2I1 = phi_T 0.0039239 ⅇ^{2.6793 + 0.0061468 V} k_I1O = phi_T 0.12052 ⅇ^{(- 9.6028) + 0.083025 V} k_OC1 = phi_T 2.391 ⅇ^{- 13.335 - (0.25289 V)} k_C1C2 = phi_T 3.1566 ⅇ^{0.36352 + 0.077193 V} k_C2C3 = phi_T 0.55432 ⅇ^{(- 0.099074) + 0.036441 V} k_C3C2 = phi_T 0.00052548 ⅇ^{(- 0.069102) + 0.0031945 V} k_C2C1 = phi_T 1.4496 ⅇ^{(- 0.1566) + 0.058353 V} k_C1O = phi_T 1.5329 ⅇ^{0.0093193 + 0.041075 V} k_OI1 = phi_T 1.6164 ⅇ^{0.30763 + 0.0060535 V} k_I1I2 = phi_T 0.027735 ⅇ^{0.05149 - (0.046865 V)} k_I1C1 = phi_T 1.9046 ⅇ^{(- 2.484) + 0.020406 V} k_C1I1 = phi_T 0.00021688 ⅇ^{(- 0.063438) + 0.0046683 V} norm = O_Na + C1 + C2 + C3 + I1 + I2 \frac{d}{d time} (C3) = \frac{(- k_C3C2) C3}{norm} + \frac{k_C2C3 C2}{norm} \frac{d}{d time} (C2) = \frac{(- (k_C2C1 + k_C2C3)) C2}{norm} + \frac{k_C1C2 C1}{norm} + \frac{k_C3C2 C3}{norm} \frac{d}{d time} (C1) = \frac{(- (k_C1C2 + k_C1O + k_C1I1)) C1}{norm} + \frac{k_OC1 O_Na}{norm} + \frac{k_C2C1 C2}{norm} + \frac{k_I1C1 I1}{norm} \frac{d}{d time} (O_Na) = \frac{(- (k_OC1 + k_OI1)) O_Na}{norm} + \frac{k_C1O C1}{norm} + \frac{k_I1O I1}{norm} \frac{d}{d time} (I2) = \frac{(- k_I2I1) I2}{norm} + \frac{k_I1I2 I1}{norm} \frac{d}{d time} (I1) = \frac{(- (k_I1O + k_I1I2 + k_I1C1)) I1}{norm} + \frac{k_I2I1 I2}{norm} + \frac{k_C1I1 C1}{norm} + \frac{k_OI1 O_Na}{norm}

Component: temperature_factor

phi = {Q10}^{\frac{T - T0}{10}}

Component: INCX

I_NCX = \frac{P_NCX (ⅇ^{\frac{gamma V F}{R T}} {Nai}^{3} Cao - (2.5 ⅇ^{\frac{(gamma - 1) V F}{R T}} {Nao}^{3} Cai))}{({K_mNai}^{3} + {Nao}^{3}) (K_mCa + Cao) (1 + k_sat ⅇ^{\frac{(gamma - 1) V F}{R T}})}

Component: INaK

I_NaK = \frac{P_NaK Ko Nai}{(Ko + K_mK) (K_mNa + Nai) (1 + 0.1245 ⅇ^{\frac{(- 0.1) V F}{R T}} + 0.0353 ⅇ^{\frac{(- V) F}{R T}})}

Component: INS_Na

Component: Ionic_currents

I = g_max PO (V - E)

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: INS_K

Component: Ionic_currents

I = g_max PO (V - E)

Component: E_Na

Component: E_Ca

Component: E_K

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: constants

Component: model_parameters

Component: initial_conditions

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Nernst_potential

E = \frac{R T}{z F} \ln (\frac{X_o}{X_i})

Component: Periodic_IStim_protocol

I_stim = Gcouple (V - V_ICC) t = time mod t_period V_ICC = \{\begin{matrix} V_ICCrest + \frac{V_ICCamp t}{f_2} if t < t_ICCpeak \\ V_ICCrest + \frac{V_ICCamp (1 + ⅇ^{\frac{- f_1}{2 t_slope}})}{1 + ⅇ^{\frac{t - f_2 - (0.5 f_1)}{t_slope}}} otherwise \end{matrix}