You are here: Home / Exposures / Katsnelson, Nikitina, Chemla, Solovyova, Coirault, Lecarpentier, Markhasin, 2004 / Time is expressed in milliseconds, calcium variables and parameters are normalised to total TnC concentration.

Model Mathematics

Component based math viewer is available

Component: environment

Component: parameters

Component: isotonic

flag = \{\begin{matrix} 0 if l \geq l_0 \land time < 150 \\ 1 otherwise \end{matrix} isotonic_mode = \{\begin{matrix} 0 if isotonic = 0 \\ 1 if isotonic = 1 \land F_muscle \geq F_afterload \\ 0 if isotonic = 1 \land l \geq l_0 \land flag = 1 \end{matrix}

Component: parameters_izakov_et_al_1991

q_v = \{\begin{matrix} q_1 - (\frac{q_2 v}{v_max}) if v ≦ 0 \\ \frac{(q_4 - q_3) v}{v_star} + q_3 if v ≦ v_star \land 0 < v \\ \frac{q_4}{{(1 + \frac{beta_Q (v - v_star)}{v_max})}^{alpha_Q}} otherwise \end{matrix} v_1 = \frac{v_max}{10}

Component: force

F_CE = lambda p_v N F_muscle = F_XSE F_SE = beta_1 (ⅇ^{alpha_1 (l_2 - l_1)} - 1) F_PE = beta_2 (ⅇ^{alpha_2 l_2} - 1) F_XSE = beta_3 (ⅇ^{alpha_3 l_3} - 1)

Component: crossbridge_kinetics

M_A = \frac{A^{mu}}{A^{mu} + {A_half}^{mu}} n_1 = 0.6 l_1 + 0.5 L_oz = \frac{l_1 + S_0}{0.46 + S_0} k_p_v = chi chi_0 q_v m_0 G_star k_m_v = \{\begin{matrix} chi_0 q_v (1 - (chi m_0 G_star)) if v ≦ v_star \\ chi_0 (q_4 (1 - (chi m_0 G_star)) + \frac{q_star (v - v_star)}{v_0 - v_star}) otherwise \end{matrix} K_chi = k_p_v M_A n_1 L_oz (1 - N) - (k_m_v N) \frac{d}{d time} (N) = K_chi

Component: length

l = l_2 + l_3 dl_1_dt = v \frac{d}{d time} (l_1) = dl_1_dt dl_2_dt = \{\begin{matrix} phi_chi_2 if k_S_vis = 0 \\ w otherwise \end{matrix} \frac{d}{d time} (l_2) = dl_2_dt dl_3_dt = \{\begin{matrix} 0 if isotonic_mode = 1 \\ - phi_chi_2 if isotonic_mode = 0 \land k_S_vis = 0 \\ - w otherwise \end{matrix} \frac{d}{d time} (l_3) = dl_3_dt

Component: CE_velocity

alp_p = \{\begin{matrix} alpha_P_lengthening if v ≦ 0 \\ alpha_P_shortening otherwise \end{matrix} k_P_vis = \{\begin{matrix} beta_P_lengthening ⅇ^{alpha_P_lengthening l_1} if v ≦ 0 \\ beta_P_shortening ⅇ^{alpha_P_shortening l_1} otherwise \end{matrix} phi_chi = \{\begin{matrix} \frac{- (lambda K_chi p_v + alp_p k_P_vis v^{2} + alpha_2 beta_2 ⅇ^{alpha_2 l_2} w)}{lambda N p_prime_v + k_P_vis} if isotonic_mode = 1 \\ \frac{- (lambda K_chi p_v + alp_p k_P_vis v^{2} + (alpha_2 beta_2 ⅇ^{alpha_2 l_2} + alpha_3 beta_3 ⅇ^{alpha_3 l_3}) w)}{lambda N p_prime_v + k_P_vis} otherwise \end{matrix} phi_chi_2 = \{\begin{matrix} \frac{alpha_1 beta_1 ⅇ^{alpha_1 (l_2 - l_1)} v}{alpha_1 beta_1 ⅇ^{alpha_1 (l_2 - l_1)} + alpha_2 beta_2 ⅇ^{alpha_2 l_2}} if isotonic_mode = 1 \\ \frac{alpha_1 beta_1 ⅇ^{alpha_1 (l_2 - l_1)} v}{alpha_1 beta_1 ⅇ^{alpha_1 (l_2 - l_1)} + alpha_2 beta_2 ⅇ^{alpha_2 l_2} + alpha_3 beta_3 ⅇ^{alpha_3 l_3}} otherwise \end{matrix} \frac{d}{d time} (v) = \{\begin{matrix} \frac{alpha_1 beta_1 ⅇ^{alpha_1 (l_2 - l_1)} (phi_chi_2 - v) - (lambda K_chi p_v + alp_p k_P_vis v^{2})}{lambda N p_prime_v + k_P_vis} if k_S_vis = 0 \\ phi_chi otherwise \end{matrix}

Component: PE_velocity

alp_s = \{\begin{matrix} alpha_S_lengthening if w ≦ v \\ alpha_S_shortening otherwise \end{matrix} k_S_vis = \{\begin{matrix} beta_S_lengthening ⅇ^{alpha_S_lengthening (l_2 - l_1)} if w ≦ v \\ beta_S_shortening ⅇ^{alpha_S_shortening (l_2 - l_1)} otherwise \end{matrix} \frac{d}{d time} (w) = \{\begin{matrix} \frac{k_S_vis (phi_chi - (alp_s {(w - v)}^{2})) - (alpha_1 beta_1 ⅇ^{alpha_1 (l_2 - l_1)} (w - v)) - (alpha_2 beta_2 ⅇ^{alpha_2 l_2} w)}{k_S_vis} if isotonic_mode = 1 \\ phi_chi - (alp_s {(w - v)}^{2}) - (\frac{alpha_1 beta_1 ⅇ^{alpha_1 (l_2 - l_1)} (w - v) + (alpha_2 beta_2 ⅇ^{alpha_2 l_2} + alpha_3 beta_3 ⅇ^{alpha_3 l_3}) w}{k_S_vis}) otherwise \end{matrix}

Component: average_crossbridge_force

gamma = \frac{a d_h {(\frac{v_1}{v_max})}^{2}}{3 a d_h - (\frac{(a + 1) v_1}{v_max})} P_star = \{\begin{matrix} \frac{a (1 + \frac{v}{v_max})}{a - (\frac{v}{v_max})} if v ≦ 0 \\ 1 + d_h - (\frac{{d_h}^{2} a}{\frac{a d_h}{gamma} {(\frac{v}{v_max})}^{2} + \frac{(a + 1) v}{v_max} + a d_h}) otherwise \end{matrix} G_star = \{\begin{matrix} 1 + \frac{0.6 v}{v_max} if - v_max ≦ v \land v ≦ 0 \\ \frac{P_star}{\frac{(0.4 a + 1) v}{a v_max} + 1} if 0 < v \land v ≦ v_1 \\ \frac{P_star ⅇ^{(- alpha_G) {(\frac{v - v_1}{v_max})}^{alpha_P}}}{\frac{(0.4 a + 1) v}{a v_max} + 1} otherwise \end{matrix} case_1 = \frac{a (0.4 + 0.4 a)}{v_max {((a + 1) 0.4)}^{2}} case_2 = \frac{a 1 (1 + 0.4 a + \frac{1.2 v}{v_max} + 0.6 {(\frac{v}{v_max})}^{2})}{v_max {((a - (\frac{v}{v_max})) (1 + \frac{0.6 v}{v_max}))}^{2}} case_3 = \frac{0.4 a + 1}{a v_max} case_4 = \frac{1}{v_max} ⅇ^{(- alpha_G) {(\frac{v}{v_max} - (\frac{v_1}{v_max}))}^{alpha_P}} (\frac{0.4 a + 1}{a} + alpha_G alpha_P (1 + \frac{(0.4 a + 1) v}{a v_max}) {(\frac{v}{v_max} - (\frac{v_1}{v_max}))}^{alpha_P - 1}) p_v = \frac{P_star}{G_star} p_prime_v = \{\begin{matrix} case_1 if v ≦ - v_max \\ case_2 if - v_max < v \land v ≦ 0 \\ case_3 if 0 < v \land v ≦ v_1 \\ case_4 otherwise \end{matrix}

Component: calcium_handling

N_A = \frac{N}{L_oz A} pi_N_A = \{\begin{matrix} 1 if N_A \geq 1 \\ {0.02}^{N_A} otherwise \end{matrix} dA_dt = a_on (A_tot - A) Ca_C - (a_off ⅇ^{(- k_A) A} pi_N_A A) \frac{d}{d time} (A) = dA_dt dB_dt = b_on (B_tot - B) Ca_C - (b_off B) \frac{d}{d time} (B) = dB_dt \frac{d}{d time} (Ca_C) = \{\begin{matrix} 4 a_c Ca_m time (1 - (ⅇ^{(- a_c) {time}^{2}})) ⅇ^{(- a_c) {time}^{2}} if time < t_d \\ - dA_dt - dB_dt - (r_Ca ⅇ^{(- q_Ca) Ca_C} Ca_C) otherwise \end{matrix}

Source: Derived from workspace Katsnelson, Nikitina, Chemla, Solovyova, Coirault, Lecarpentier, Markhasin, 2004 at changeset ff37a293c456.

Collaboration: To begin collaborating on this work, please use your git client and issue this command:

Downloads

Views Available

Tools

License: This work is licensed under a Creative Commons Attribution 3.0 Unported License.