You are here: Home / Exposures / Model of the human Purkinje cell by Sampson et al. / sampson-2010.cellml

Model Mathematics

Component based math viewer is available

Component: caflux

\frac{d}{d time} (HTRPNCa) = khtrpn_plus Cai (1.0 - HTRPNCa) - caflux_HTRPNCa_a1 Jtr = \frac{CaNSR - CaJSR}{tautr} Jtrpn = LTRPNtot \frac{d}{d time} (LTRPNCa) + HTRPNtot \frac{d}{d time} (HTRPNCa) Jxfer = \frac{CaSS - Cai}{tauxfer} \frac{d}{d time} (LTRPNCa) = kltrpn_plus Cai (1.0 - LTRPNCa) - caflux_LTRPNCa_a1 caflux_HTRPNCa_a1 = khtrpn_minus HTRPNCa caflux_LTRPNCa_a1 = kltrpn_minus LTRPNCa

Component: calcium

\frac{d}{d time} (CaJSR) = beta_JSR (Jtr - Jrel) \frac{d}{d time} (CaNSR) = \frac{Jup Vmyo}{VNSR} - (\frac{Jtr VJSR}{VNSR}) \frac{d}{d time} (CaSS) = beta_SS (calcium_CaSS_a3 - (ICa cell_a2) - (icat_I cell_a2)) \frac{d}{d time} (Cai) = beta_i (Jxfer - Jup - Jtrpn - (calcium_Cai_a3 0.5 cell_a1)) beta_JSR = \frac{1.0}{1.0 + calcium_beta_JSR_b1} beta_SS = \frac{1.0}{1.0 + calcium_beta_SS_b1 + calcium_beta_SS_b2} beta_i = \frac{1.0}{1.0 + calcium_beta_i_b1 + calcium_beta_i_b2} calcium_CaSS_a3 = \frac{Jrel VJSR}{VSS} - (\frac{Jxfer Vmyo}{VSS}) calcium_Cai_a3 = (- 2.0) INaCa + IpCa calcium_beta_JSR_b1 = \frac{CSQNtot KmCSQN}{{(CaJSR + KmCSQN)}^{2.0}} calcium_beta_SS_b1 = \frac{CMDNtot KmCMDN}{{(CaSS + KmCMDN)}^{2.0}} calcium_beta_SS_b2 = \frac{EGTAtot KmEGTA}{{(CaSS + KmEGTA)}^{2.0}} calcium_beta_i_b1 = \frac{CMDNtot KmCMDN}{{(Cai + KmCMDN)}^{2.0}} calcium_beta_i_b2 = \frac{EGTAtot KmEGTA}{{(Cai + KmEGTA)}^{2.0}}

Component: cell

a1 = \frac{Acap}{Vmyo F 1000.0} a2 = \frac{Acap}{2.0 VSS F 1000.0}

Component: engine

Component: extra

Component: ical

\frac{d}{d time} (C0) = ical_C0_a2 - ical_C0_a1 C0_to_C1 = 4.0 alpha C0_to_CCa0 = gamma \frac{d}{d time} (C1) = ical_C1_a2 - ical_C1_a1 C1_to_C0 = beta C1_to_C2 = 3.0 alpha C1_to_CCa1 = aL C0_to_CCa0 \frac{d}{d time} (C2) = ical_C2_a2 - ical_C2_a1 C2_to_C1 = 2.0 beta C2_to_C3 = 2.0 alpha C2_to_CCa2 = aL C1_to_CCa1 \frac{d}{d time} (C3) = ical_C3_a2 - ical_C3_a1 C3_to_C2 = 3.0 beta C3_to_C4 = alpha C3_to_CCa3 = aL C2_to_CCa2 \frac{d}{d time} (C4) = ical_C4_a2 - ical_C4_a1 C4_to_C3 = 4.0 beta C4_to_CCa4 = aL C3_to_CCa3 \frac{d}{d time} (CCa0) = ical_CCa0_a2 - ical_CCa0_a1 CCa0_to_C0 = omega CCa0_to_CCa1 = 4.0 alpha_prime \frac{d}{d time} (CCa1) = ical_CCa1_a2 - ical_CCa1_a1 CCa1_to_C1 = \frac{CCa0_to_C0}{bL} CCa1_to_CCa0 = beta_prime CCa1_to_CCa2 = 3.0 alpha_prime \frac{d}{d time} (CCa2) = ical_CCa2_a2 - ical_CCa2_a1 CCa2_to_C2 = \frac{CCa1_to_C1}{bL} CCa2_to_CCa1 = 2.0 beta_prime CCa2_to_CCa3 = 2.0 alpha_prime \frac{d}{d time} (CCa3) = ical_CCa3_a2 - ical_CCa3_a1 CCa3_to_C3 = \frac{CCa2_to_C2}{bL} CCa3_to_CCa2 = 3.0 beta_prime CCa3_to_CCa4 = alpha_prime \frac{d}{d time} (CCa4) = ical_CCa4_a2 - ical_CCa4_a1 CCa4_to_C4 = \frac{CCa3_to_C3}{bL} CCa4_to_CCa3 = 4.0 beta_prime ICa = ICamax yCa Open ICaK = PKprime Open yCa VFFRT \frac{ical_ICaK_a1}{ical_ICaK_a2} ICahalf = - 0.265 ICamax = PCa 4.0 VFFRT \frac{ical_ICamax_a1}{ical_ICamax_a2} \frac{d}{d time} (Open) = fL C4 - (gL Open) PCa = Pscale 0.2469 PK = Pscale 0.0004574 PKprime = \frac{PK}{1.0 + \frac{imax}{ICahalf}} alpha = 4.0 1.2 0.416 ⅇ^{0.012 (V - 35.0)} alpha_prime = aL alpha beta = 4.0 0.45 0.049 ⅇ^{(- 0.065) (V - 22.0)} beta_prime = \frac{beta}{bL} gamma = 0.6 0.09233 CaSS ical_C0_a1 = (C0_to_C1 + C0_to_CCa0) C0 ical_C0_a2 = C1_to_C0 C1 + CCa0_to_C0 CCa0 ical_C1_a1 = (C1_to_C0 + C1_to_C2 + C1_to_CCa1) C1 ical_C1_a2 = C0_to_C1 C0 + C2_to_C1 C2 + CCa1_to_C1 CCa1 ical_C2_a1 = (C2_to_C1 + C2_to_C3 + C2_to_CCa2) C2 ical_C2_a2 = C1_to_C2 C1 + C3_to_C2 C3 + CCa2_to_C2 CCa2 ical_C3_a1 = (C3_to_C2 + C3_to_C4 + C3_to_CCa3) C3 ical_C3_a2 = C2_to_C3 C2 + C4_to_C3 C4 + CCa3_to_C3 CCa3 ical_C4_a1 = (C4_to_C3 + fL + C4_to_CCa4) C4 ical_C4_a2 = C3_to_C4 C3 + gL Open + CCa4_to_C4 CCa4 ical_CCa0_a1 = (CCa0_to_CCa1 + CCa0_to_C0) CCa0 ical_CCa0_a2 = CCa1_to_CCa0 CCa1 + C0_to_CCa0 C0 ical_CCa1_a1 = (CCa1_to_CCa0 + CCa1_to_CCa2 + CCa1_to_C1) CCa1 ical_CCa1_a2 = CCa0_to_CCa1 CCa0 + CCa2_to_CCa1 CCa2 + C1_to_CCa1 C1 ical_CCa2_a1 = (CCa2_to_CCa1 + CCa2_to_CCa3 + CCa2_to_C2) CCa2 ical_CCa2_a2 = CCa1_to_CCa2 CCa1 + CCa3_to_CCa2 CCa3 + C2_to_CCa2 C2 ical_CCa3_a1 = (CCa3_to_CCa2 + CCa3_to_CCa4 + CCa3_to_C3) CCa3 ical_CCa3_a2 = CCa2_to_CCa3 CCa2 + CCa4_to_CCa3 CCa4 + C3_to_CCa3 C3 ical_CCa4_a1 = (CCa4_to_CCa3 + CCa4_to_C4) CCa4 ical_CCa4_a2 = CCa3_to_CCa4 CCa3 + C4_to_CCa4 C4 ical_ICaK_a1 = Ki ⅇ^{VFRT} - Ko ical_ICaK_a2 = ⅇ^{VFRT} - 1.0 ical_ICamax_a1 = 0.001 ⅇ^{2.0 VFRT} - (Cao 0.341) ical_ICamax_a2 = ⅇ^{2.0 VFRT} - 1.0 imax = \{\begin{matrix} 0.0 if ICamax > 0.0 \\ ICamax otherwise \end{matrix} tau_yCa = \frac{1.0}{\frac{3.36336209451999985 e -3}{0.5 + ⅇ^{\frac{V}{- 5.53899874036055007}}} + 7.79046570737000014 e -3 ⅇ^{\frac{V}{- 4.95103963116038628 e 1}}} \frac{d}{d time} (yCa) = \frac{yCa_inf - yCa}{tau_yCa} yCa_inf = \frac{ical_yCa_yCa_inf_a1}{1.0 + ⅇ^{\frac{V + 28.5}{7.8}}} + (1.0 - ical_yCa_yCa_inf_a1)

Component: icat

I = \frac{ICamax}{PCa} Ttypescale n n l icat_l_inf = \frac{1.0}{1.0 + ⅇ^{\frac{V + 75.6}{6.2}}} icat_l_tau = \{\begin{matrix} 500.0 if V < - 60.0 \\ 18.3 + 0.005 ⅇ^{\frac{- V}{6.2}} otherwise \end{matrix} icat_n_inf = \frac{1.0}{1.0 + ⅇ^{\frac{- (V + 48.4)}{5.2}}} icat_n_tau = \{\begin{matrix} 2.44 ⅇ^{\frac{V + 120.0}{40.0}} if V < - 56.0 \\ 1.34 + 0.035 ⅇ^{\frac{- V}{11.8}} otherwise \end{matrix} \frac{d}{d time} (l) = \frac{icat_l_inf - l}{icat_l_tau} \frac{d}{d time} (n) = \frac{icat_n_inf - n}{icat_n_tau}

Component: ihcn

IHCN = IHCNmax (hcn6 + hcn7 + hcn8 + hcn9 + hcn10) (V - (\frac{ENa}{3.0} + \frac{2.0 EK}{3.0})) h_alpha = \frac{0.0001712 ⅇ^{(- 1.465) VFRT}}{1000.0} h_beta = \frac{26.17 ⅇ^{1.465 VFRT}}{1000.0} h_delta = \frac{287.5 ⅇ^{1.242 VFRT}}{1000.0} h_gamma = \frac{0.04025 ⅇ^{(- 1.242) VFRT}}{1000.0} \frac{d}{d time} (hcn1) = (- (4.0 h_gamma + h_alpha)) hcn1 + h_delta hcn2 + h_beta hcn6 \frac{d}{d time} (hcn10) = (- (\frac{h_beta}{{h_f}^{4.0}} + \frac{4.0 h_delta}{h_f})) hcn10 + h_alpha {h_f}^{4.0} hcn5 + h_gamma h_f hcn9 \frac{d}{d time} (hcn2) = (- (h_delta + 3.0 h_gamma + h_alpha h_f)) hcn2 + 4.0 h_gamma hcn1 + 2.0 h_delta hcn3 + \frac{h_beta}{h_f} hcn7 \frac{d}{d time} (hcn3) = (- (2.0 h_delta + 2.0 h_gamma + h_alpha {h_f}^{2.0})) hcn3 + 3.0 h_gamma hcn2 + 3.0 h_delta hcn4 + \frac{h_beta}{{h_f}^{2.0}} hcn8 \frac{d}{d time} (hcn4) = (- (3.0 h_delta + h_gamma + h_alpha {h_f}^{3.0})) hcn4 + 2.0 h_gamma hcn3 + 4.0 h_delta hcn5 + \frac{h_beta}{{h_f}^{3.0}} hcn9 \frac{d}{d time} (hcn5) = (- (4.0 h_delta + h_alpha {h_f}^{4.0})) hcn5 + h_gamma hcn4 + \frac{h_beta}{{h_f}^{4.0}} hcn10 \frac{d}{d time} (hcn6) = (- (h_beta + 4.0 h_gamma h_f)) hcn6 + h_alpha hcn1 + \frac{h_delta}{h_f} hcn7 \frac{d}{d time} (hcn7) = (- (\frac{h_beta}{h_f} + \frac{h_delta}{h_f} + 3.0 h_gamma h_f)) hcn7 + h_alpha h_f hcn2 + 4.0 h_gamma h_f hcn6 + \frac{2.0 h_delta}{h_f} hcn8 \frac{d}{d time} (hcn8) = (- (\frac{h_beta}{{h_f}^{2.0}} + \frac{2.0 h_delta}{h_f} + 2.0 h_gamma h_f)) hcn8 + h_alpha {h_f}^{2.0} hcn3 + 3.0 h_gamma h_f hcn7 + \frac{3.0 h_delta}{h_f} hcn9 \frac{d}{d time} (hcn9) = (- (\frac{h_beta}{{h_f}^{3.0}} + \frac{3.0 h_delta}{h_f} + h_gamma h_f)) hcn9 + h_alpha {h_f}^{3.0} hcn4 + 2.0 h_gamma h_f hcn8 + \frac{4.0 h_delta}{h_f} hcn10

Component: ik1

IK1 = GK1 \sqrt{Ko 1.0} ik1_IK1_inf (V - EK) ik1_IK1_inf = \frac{1.0}{0.9681 + ⅇ^{\frac{V + 82.18623}{15.88642}}}

Component: ikr

B1 = - 4.30605416397999985 e -2 B3 = - 2.69138549839900015 e -2 B5 = - 4.53664295954299968 e -2 \frac{d}{d time} (C1) = C2H_to_C1H C2 - (C1H_to_C2H C1) C1H_to_C2H = T_Const A0 ⅇ^{B0 V} \frac{d}{d time} (C2) = ikr_C2_a1 - ikr_C2_a2 C2H_to_C1H = T_Const A1 ⅇ^{B1 V} C2H_to_C3H = T_Const 2.60836204333700014 e -2 \frac{d}{d time} (C3) = ikr_C3_a1 - ikr_C3_a2 C3H_to_C2H = T_Const 1.48329781321449999 e -1 C3H_to_IH = T_Const A6 ⅇ^{B6 V} C3H_to_OH = T_Const A2 ⅇ^{B2 V} \frac{d}{d time} (I) = ikr_I_a1 - ikr_I_a2 IH_to_C3H = \frac{OH_to_C3H IH_to_OH C3H_to_IH}{C3H_to_OH OH_to_IH} IH_to_OH = T_Const A5 ⅇ^{B5 V} IKr = GKr fKo O (V - EK) \frac{d}{d time} (O) = ikr_O_a1 - ikr_O_a2 OH_to_C3H = T_Const A3 ⅇ^{B3 V} OH_to_IH = T_Const A4 ⅇ^{B4 V} fKo = \sqrt{\frac{Ko}{4.0}} ikr_C2_a1 = C1H_to_C2H C1 + C3H_to_C2H C3 ikr_C2_a2 = (C2H_to_C1H + C2H_to_C3H) C2 ikr_C3_a1 = C2H_to_C3H C2 + OH_to_C3H O + IH_to_C3H I ikr_C3_a2 = (C3H_to_IH + C3H_to_OH + C3H_to_C2H) C3 ikr_I_a1 = C3H_to_IH C3 + OH_to_IH O ikr_I_a2 = (IH_to_C3H + IH_to_OH) I ikr_O_a1 = C3H_to_OH C3 + IH_to_OH I ikr_O_a2 = (OH_to_C3H + OH_to_IH) O

Component: iks

IKs = GKs xs_wt xf_wt (V - EK) iks_xf_wt_alpha = \{\begin{matrix} \frac{0.000146}{0.078} if V = 21.0 \\ \frac{0.000146 (V - 21.0)}{1.0 - (ⅇ^{(- 0.078) (V - 21.0)})} otherwise \end{matrix} iks_xf_wt_beta = 0.00091 ⅇ^{(- 0.028) V} iks_xs_wt_alpha = \{\begin{matrix} \frac{3.3 e -5}{0.13} if V = - 11.0 \\ \frac{3.3 e -5 (V + 11.0)}{1.0 - (ⅇ^{(- 0.13) (V + 11.0)})} otherwise \end{matrix} iks_xs_wt_beta = 0.0001 ⅇ^{(- 0.015) V} \frac{d}{d time} (xf_wt) = iks_xf_wt_alpha (1.0 - xf_wt) - (iks_xf_wt_beta xf_wt) \frac{d}{d time} (xs_wt) = iks_xs_wt_alpha (1.0 - xs_wt) - (iks_xs_wt_beta xs_wt)

Component: inaca

INaCa = \frac{kNaCa a5 (a1 - a2)}{a4 a3} a1 = ⅇ^{eta VFRT} {Nai}^{3.0} Cao a2 = ⅇ^{(eta - 1.0) VFRT} nao3 Cai a3 = 1.0 + ksat ⅇ^{(eta - 1.0) VFRT} a4 = KmCa + Cao a5 = \frac{5000.0}{{KmNa}^{3.0} + nao3} nao3 = {Nao}^{3.0}

Component: inak

INaK = INaKmax fNaK \frac{inak_INaK_a1}{inak_INaK_a2} fNaK = \frac{1.0}{inak_fNaK_a1 + inak_fNaK_a2} inak_INaK_a1 = \frac{Ko}{Ko + KmKo} inak_INaK_a2 = 1.0 + {(\frac{KmNai}{Nai})}^{1.5} inak_fNaK_a1 = 1.0 + 0.1245 ⅇ^{(- 0.1) VFRT} inak_fNaK_a2 = 0.0365 sigma ⅇ^{(- 1.33) VFRT} sigma = \frac{ⅇ^{\frac{Nao}{67.3}} - 1.0}{7.0}

Component: ipca

IpCa = \frac{IpCamax Cai}{KmpCa + Cai}

Component: isus

Isus = Isusmax \frac{1.0}{1.0 + ⅇ^{\frac{- (V - 20.0)}{12.0}}} (V - EK)

Component: ito

\frac{d}{d time} (C0) = ito_C0_a2 - ito_C0_a1 C0_to_C1 = 4.0 alpha_act43 C0_to_CI0 = beta_inact43 \frac{d}{d time} (C1) = ito_C1_a2 - ito_C1_a1 C1_to_C0 = beta_act43 C1_to_C2 = 3.0 alpha_act43 C1_to_CI1 = f1 beta_inact43 \frac{d}{d time} (C2) = ito_C2_a2 - ito_C2_a1 C2_to_C1 = 2.0 beta_act43 C2_to_C3 = 2.0 alpha_act43 C2_to_CI2 = f2 beta_inact43 \frac{d}{d time} (C3) = ito_C3_a2 - ito_C3_a1 C3_to_C2 = 3.0 beta_act43 C3_to_CI3 = f3 beta_inact43 C3_to_O = alpha_act43 \frac{d}{d time} (CI0) = ito_CI0_a2 - ito_CI0_a1 CI0_to_C0 = alpha_inact43 CI0_to_CI1 = 4.0 b1 alpha_act43 \frac{d}{d time} (CI1) = ito_CI1_a2 - ito_CI1_a1 CI1_to_C1 = \frac{alpha_inact43}{b1} CI1_to_CI0 = \frac{beta_act43}{f1} CI1_to_CI2 = \frac{3.0 b2 alpha_act43}{b1} \frac{d}{d time} (CI2) = ito_CI2_a2 - ito_CI2_a1 CI2_to_C2 = \frac{alpha_inact43}{b2} CI2_to_CI1 = \frac{2.0 f1 beta_act43}{f2} CI2_to_CI3 = \frac{2.0 b3 alpha_act43}{b2} \frac{d}{d time} (CI3) = ito_CI3_a2 - ito_CI3_a1 CI3_to_C3 = \frac{alpha_inact43}{b3} CI3_to_CI2 = \frac{3.0 f2 beta_act43}{f3} CI3_to_OI = \frac{b4 alpha_act43}{b3} Ito1 = G O (V - EK) \frac{d}{d time} (O) = ito_O_a2 - ito_O_a1 \frac{d}{d time} (OI) = ito_OI_a2 - ito_OI_a1 OI_to_CI3 = \frac{4.0 f3 beta_act43}{f4} OI_to_O = \frac{alpha_inact43}{b4} O_to_C3 = 4.0 beta_act43 O_to_OI = f4 beta_inact43 alpha_act43 = alphaa0 ⅇ^{aa V} alpha_inact43 = alphai0 ⅇ^{(- ai) V} beta_act43 = betaa0 ⅇ^{(- ba) V} beta_inact43 = betai0 ⅇ^{bi V} ito_C0_a1 = (C0_to_C1 + C0_to_CI0) C0 ito_C0_a2 = C1_to_C0 C1 + CI0_to_C0 CI0 ito_C1_a1 = (C1_to_C2 + C1_to_C0 + C1_to_CI1) C1 ito_C1_a2 = C2_to_C1 C2 + CI1_to_C1 CI1 + C0_to_C1 C0 ito_C2_a1 = (C2_to_C3 + C2_to_C1 + C2_to_CI2) C2 ito_C2_a2 = C3_to_C2 C3 + CI2_to_C2 CI2 + C1_to_C2 C1 ito_C3_a1 = (C3_to_O + C3_to_C2 + C3_to_CI3) C3 ito_C3_a2 = O_to_C3 O + CI3_to_C3 CI3 + C2_to_C3 C2 ito_CI0_a1 = (CI0_to_C0 + CI0_to_CI1) CI0 ito_CI0_a2 = C0_to_CI0 C0 + CI1_to_CI0 CI1 ito_CI1_a1 = (CI1_to_CI2 + CI1_to_C1 + CI1_to_CI0) CI1 ito_CI1_a2 = CI2_to_CI1 CI2 + C1_to_CI1 C1 + CI0_to_CI1 CI0 ito_CI2_a1 = (CI2_to_CI3 + CI2_to_C2 + CI2_to_CI1) CI2 ito_CI2_a2 = CI3_to_CI2 CI3 + C2_to_CI2 C2 + CI1_to_CI2 CI1 ito_CI3_a1 = (CI3_to_OI + CI3_to_C3 + CI3_to_CI2) CI3 ito_CI3_a2 = OI_to_CI3 OI + C3_to_CI3 C3 + CI2_to_CI3 CI2 ito_OI_a1 = (OI_to_O + OI_to_CI3) OI ito_OI_a2 = O_to_OI O + CI3_to_OI CI3 ito_O_a1 = (O_to_C3 + O_to_OI) O ito_O_a2 = C3_to_O C3 + OI_to_O OI

Component: membrane

\frac{d}{d time} (V) = - (i_ion + i_diff + i_stim) a1 = INa + INa1 + ICa + icat_I + ICaK + IKr + IKs + IHCN a2 = IK1 + INaCa + INaK + Ito1 + Isus a3 = IpCa amplitude = - 100.0 i_ion = a1 + a2 + a3 i_stim = \{\begin{matrix} 1.0 if time - offset - (period ⌊ \frac{time - offset}{period} ⌋) < duration \\ 0.0 otherwise \end{matrix} amplitude

Component: nav11

\frac{d}{d time} (BC1) = BO b13 + BC2 a12 + C1 mu1 - (BC1 (a13 + b12 + mu2)) \frac{d}{d time} (BC2) = BC1 b12 + BC3 a11 + C2 mu1 - (BC2 (a12 + b11 + mu2)) \frac{d}{d time} (BC3) = BC2 b11 + C3 mu1 - (BC3 (a11 + mu2)) \frac{d}{d time} (BO) = BC1 a13 + O mu1 - (BO (b13 + mu2)) \frac{d}{d time} (C1) = IF a3 + BC1 mu2 + O b13 + C2 a12 - (C1 (a13 + b3 + mu1 + b12)) \frac{d}{d time} (C2) = IC2 a3 + BC2 mu2 + C1 b12 + C3 a11 - (C2 (a12 + b3 + mu1 + b11)) \frac{d}{d time} (C3) = IC3 a3 + BC3 mu2 + C2 b11 - (C3 (a11 + b3 + mu1)) \frac{d}{d time} (IC2) = IF b12 + C2 b3 + IC3 a11 - (IC2 (a12 + a3 + b11)) \frac{d}{d time} (IC3) = IC2 b11 + C3 b3 - (IC3 (a11 + a3)) \frac{d}{d time} (IF) = IS1 b4 + O a2 + IC2 a12 + C1 b3 - (IF (a4 + b2 + b12 + a3)) INa1 = GNa1 (O + BO) (V - ENa) \frac{d}{d time} (IS1) = IS2 b5 + IF a4 - (IS1 (a5 + b4)) \frac{d}{d time} (IS2) = IS1 a5 - (IS2 b5) \frac{d}{d time} (O) = IF b2 + BO mu2 + C1 a13 - (O (a2 + b13 + mu1)) a11 = \frac{2.802}{0.21 ⅇ^{\frac{- V}{17.0}} + 0.23 ⅇ^{\frac{- V}{150.0}}} a12 = \frac{2.802}{0.23 ⅇ^{\frac{- V}{15.0}} + 0.25 ⅇ^{\frac{- V}{150.0}}} a13 = \frac{2.802}{0.25 ⅇ^{\frac{- V}{12.0}} + 0.27 ⅇ^{\frac{- V}{150.0}}} a2 = \frac{9.178 ⅇ^{\frac{V}{29.68}}}{4.5} a3 = 3.79329999999999996 e -7 ⅇ^{\frac{- V}{7.6}} 3.0 a4 = \frac{a2}{100.0} 1.5 0.285 a5 = \frac{a2}{95000.0} 80.0 b11 = 0.4 ⅇ^{\frac{- V}{20.3}} b12 = 0.4 ⅇ^{\frac{- (V - 5.0)}{20.3}} b13 = \frac{0.4 ⅇ^{\frac{- (V - 10.0)}{20.3}}}{4.5} b2 = \frac{a3 a2 a13}{b3 b13} b3 = 0.0084 + 2 e -5 V b4 = \frac{a3}{5.0} b5 = \frac{\frac{a3}{30.0}}{10.0}

Component: nav15

\frac{d}{d time} (BC1) = BO b13 + BC2 a12 + C1 mu1 - (BC1 (a13 + b12 + mu2)) \frac{d}{d time} (BC2) = BC1 b12 + BC3 a11 + C2 mu1 - (BC2 (a12 + b11 + mu2)) \frac{d}{d time} (BC3) = BC2 b11 + C3 mu1 - (BC3 (a11 + mu2)) \frac{d}{d time} (BO) = BC1 a13 + O mu1 - (BO (b13 + mu2)) \frac{d}{d time} (C1) = IF a3 + BC1 mu2 + O b13 + C2 a12 - (C1 (a13 + b3 + mu1 + b12)) \frac{d}{d time} (C2) = IC2 a3 + BC2 mu2 + C1 b12 + C3 a11 - (C2 (a12 + b3 + mu1 + b11)) \frac{d}{d time} (C3) = IC3 a3 + BC3 mu2 + C2 b11 - (C3 (a11 + b3 + mu1)) \frac{d}{d time} (IC2) = IF b12 + C2 b3 + IC3 a11 - (IC2 (a12 + a3 + b11)) \frac{d}{d time} (IC3) = IC2 b11 + C3 b3 - (IC3 (a11 + a3)) \frac{d}{d time} (IF) = IS1 b4 + O a2 + IC2 a12 + C1 b3 - (IF (a4 + b2 + b12 + a3)) INa = GNa (O + BO) (V - ENa) \frac{d}{d time} (IS1) = IS2 b5 + IF a4 - (IS1 (a5 + b4)) \frac{d}{d time} (IS2) = IS1 a5 - (IS2 b5) \frac{d}{d time} (O) = IF b2 + BO mu2 + C1 a13 - (O (a2 + b13 + mu1)) a11 = \frac{3.802}{0.1027 ⅇ^{\frac{- V}{17.0}} + 0.2 ⅇ^{\frac{- V}{150.0}}} a12 = \frac{3.802}{0.1027 ⅇ^{\frac{- V}{15.0}} + 0.23 ⅇ^{\frac{- V}{150.0}}} a13 = \frac{3.802}{0.1027 ⅇ^{\frac{- V}{12.0}} + 0.25 ⅇ^{\frac{- V}{150.0}}} a2 = 9.178 ⅇ^{\frac{V}{29.68}} a3 = 3.79329999999999996 e -7 ⅇ^{\frac{- V}{7.7}} a4 = \frac{a2}{100.0} a5 = \frac{a2}{95000.0} b11 = 0.1917 ⅇ^{\frac{- V}{20.3}} b12 = 0.2 ⅇ^{\frac{- (V - 5.0)}{20.3}} b13 = 0.22 ⅇ^{\frac{- (V - 10.0)}{20.3}} b2 = \frac{a3 a2 a13}{b3 b13} b3 = 0.0084 + 2 e -5 V b4 = a3 b5 = \frac{a3}{50.0}

Component: nernst

EK = RTF \ln (\frac{Ko}{Ki}) ENa = RTF \ln (\frac{Nao}{Nai})

Component: phys

RTF = \frac{R T}{F} VFFRT = F VFRT VFRT = \frac{V}{RTF}

Component: potassium

IK_tot = (IKr + IKs + IK1 + ICaK + \frac{2.0 IHCN}{3.0} + Isus + i_stim - (2.0 INaK)) + Ito1 \frac{d}{d time} (Ki) = (- cell_a1) IK_tot

Component: ryr

\frac{d}{d time} (C1) = (- kaplus) ryr_C1_a2 C1 + kaminus O1 \frac{d}{d time} (C2) = kcplus O1 - (kcminus C2) Jrel = v1 (O1 + O2) (CaJSR - CaSS) \frac{d}{d time} (O1) = - (\frac{d}{d time} (C1) + \frac{d}{d time} (O2) + \frac{d}{d time} (C2)) \frac{d}{d time} (O2) = kbplus ryr_O2_a1 O1 - (kbminus O2) ryr_C1_a2 = {(CaSS 1000.0)}^{4.0} ryr_O2_a1 = {(CaSS 1000.0)}^{3.0}

Component: serca

Jup = \frac{KSR (vmaxf fb - (vmaxr rb))}{1.0 + fb + rb} fb = {(\frac{Cai}{Kfb})}^{Nfb} rb = {(\frac{CaNSR}{Krb})}^{Nrb}

Component: sodium

INa_tot = INa + INa1 + \frac{IHCN}{3.0} + 3.0 (INaCa + INaK) \frac{d}{d time} (Nai) = (- cell_a1) INa_tot

Source: Derived from workspace Sampson-Iyer-Marks-Kass-2010 at changeset 7e0450c5e3cd.

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

Downloads

Views Available

Tools

License: The terms of use/license for this work is unspecified.

sampson-2010.cellml