You are here: Home / Exposures / O'Hara-Rudy-CiPA-v1.0 (2017) / O'Hara-Rudy CiPA v1.0 (2017)

Model Mathematics

Component based math viewer is available

Component: environment

Component: extracellular

Component: physical_constants

Component: cell_geometry

vcell = 1000 3.14 rad rad L Ageo = 2 3.14 rad rad + 2 3.14 rad L Acap = 2 Ageo vmyo = 0.68 vcell vnsr = 0.0552 vcell vjsr = 0.0048 vcell vss = 0.02 vcell

Component: membrane

Istim = \{\begin{matrix} i_Stim_Amplitude if time \geq i_Stim_Start \land time ≦ i_Stim_End \land time - i_Stim_Start - (⌊ \frac{time - i_Stim_Start}{i_Stim_Period} ⌋ i_Stim_Period) ≦ i_Stim_PulseDuration \\ 0 otherwise \end{matrix} frt = \frac{F}{R T} ffrt = F frt vfrt = v frt \frac{d}{d time} (v) = - (INa + INaL + Ito + ICaL + ICaNa + ICaK + IKr + IKs + IK1 + INaCa_i + INaCa_ss + INaK + INab + IKb + IpCa + ICab + Istim)

Component: CaMK

CaMKb = \frac{CaMKo (1 - CaMKt)}{1 + \frac{KmCaM}{cass}} CaMKa = CaMKb + CaMKt \frac{d}{d time} (CaMKt) = aCaMK CaMKb (CaMKb + CaMKt) - (bCaMK CaMKt)

Component: intracellular_ions

cmdnmax = \{\begin{matrix} cmdnmax_b 1.3 if celltype = 1 \\ cmdnmax_b otherwise \end{matrix} \frac{d}{d time} (nai) = \frac{(- (INa + INaL + 3 INaCa_i + 3 INaK + INab)) Acap cm}{F vmyo} + \frac{JdiffNa vss}{vmyo} \frac{d}{d time} (nass) = \frac{(- (ICaNa + 3 INaCa_ss)) cm Acap}{F vss} - JdiffNa \frac{d}{d time} (ki) = \frac{(- (Ito + IKr + IKs + IK1 + IKb + Istim - (2 INaK))) cm Acap}{F vmyo} + \frac{JdiffK vss}{vmyo} \frac{d}{d time} (kss) = \frac{(- ICaK) cm Acap}{F vss} - JdiffK Bcai = \frac{1}{1 + \frac{cmdnmax kmcmdn}{{(kmcmdn + cai)}^{2}} + \frac{trpnmax kmtrpn}{{(kmtrpn + cai)}^{2}}} \frac{d}{d time} (cai) = Bcai ((\frac{(- (IpCa + ICab - (2 INaCa_i))) cm Acap}{2 F vmyo} - (\frac{Jup vnsr}{vmyo})) + \frac{Jdiff vss}{vmyo}) Bcass = \frac{1}{1 + \frac{BSRmax KmBSR}{{(KmBSR + cass)}^{2}} + \frac{BSLmax KmBSL}{{(KmBSL + cass)}^{2}}} \frac{d}{d time} (cass) = Bcass (\frac{(- (ICaL - (2 INaCa_ss))) cm Acap}{2 F vss} + \frac{Jrel vjsr}{vss} - Jdiff) \frac{d}{d time} (cansr) = Jup - (\frac{Jtr vjsr}{vnsr}) Bcajsr = \frac{1}{1 + \frac{csqnmax kmcsqn}{{(kmcsqn + cajsr)}^{2}}} \frac{d}{d time} (cajsr) = Bcajsr (Jtr - Jrel)

Component: reversal_potentials

ENa = \frac{R T}{F} \ln (\frac{nao}{nai}) EK = \frac{R T}{F} \ln (\frac{ko}{ki}) EKs = \frac{R T}{F} \ln (\frac{ko + PKNa nao}{ki + PKNa nai})

Component: INa

mss = \frac{1}{1 + ⅇ^{\frac{- (v + mssV1)}{mssV2}}} tm = \frac{1}{mtD1 ⅇ^{\frac{v + mtV1}{mtV2}} + mtD2 ⅇ^{\frac{- (v + mtV3)}{mtV4}}} \frac{d}{d time} (m) = \frac{mss - m}{tm} hss = \frac{1}{1 + ⅇ^{\frac{v + hssV1 - shift_INa_inact}{hssV2}}} thf = \frac{1}{1.432 e -5 ⅇ^{\frac{- (v + 1.196 - shift_INa_inact)}{6.285}} + 6.149 ⅇ^{\frac{v + 0.5096 - shift_INa_inact}{20.27}}} ths = \frac{1}{0.009794 ⅇ^{\frac{- (v + 17.95 - shift_INa_inact)}{28.05}} + 0.3343 ⅇ^{\frac{v + 5.73 - shift_INa_inact}{56.66}}} Ahs = 1 - Ahf \frac{d}{d time} (hf) = \frac{hss - hf}{thf} \frac{d}{d time} (hs) = \frac{hss - hs}{ths} h = Ahf hf + Ahs hs jss = hss tj = 2.038 + \frac{1}{0.02136 ⅇ^{\frac{- (v + 100.6 - shift_INa_inact)}{8.281}} + 0.3052 ⅇ^{\frac{v + 0.9941 - shift_INa_inact}{38.45}}} \frac{d}{d time} (j) = \frac{jss - j}{tj} hssp = \frac{1}{1 + ⅇ^{\frac{v + 89.1 - shift_INa_inact}{6.086}}} thsp = 3 ths \frac{d}{d time} (hsp) = \frac{hssp - hsp}{thsp} hp = Ahf hf + Ahs hsp tjp = 1.46 tj \frac{d}{d time} (jp) = \frac{jss - jp}{tjp} fINap = \frac{1}{1 + \frac{KmCaMK}{CaMKa}} INa = GNa (v - ENa) m^{3} ((1 - fINap) h j + fINap hp jp)

Component: INaL

mLss = \frac{1}{1 + ⅇ^{\frac{- (v + 42.85)}{5.264}}} tmL = tm \frac{d}{d time} (mL) = \frac{mLss - mL}{tmL} hLss = \frac{1}{1 + ⅇ^{\frac{v + 87.61}{7.488}}} \frac{d}{d time} (hL) = \frac{hLss - hL}{thL} hLssp = \frac{1}{1 + ⅇ^{\frac{v + 93.81}{7.488}}} thLp = 3 thL \frac{d}{d time} (hLp) = \frac{hLssp - hLp}{thLp} GNaL = \{\begin{matrix} GNaL_b 0.6 if celltype = 1 \\ GNaL_b otherwise \end{matrix} fINaLp = \frac{1}{1 + \frac{KmCaMK}{CaMKa}} INaL = GNaL (v - ENa) mL ((1 - fINaLp) hL + fINaLp hLp)

Component: Ito

ass = \frac{1}{1 + ⅇ^{\frac{- (v - 14.34)}{14.82}}} ta = \frac{1.0515}{\frac{1}{1.2089 (1 + ⅇ^{\frac{- (v - 18.4099)}{29.3814}})} + \frac{3.5}{1 + ⅇ^{\frac{v + 100}{29.3814}}}} \frac{d}{d time} (a) = \frac{ass - a}{ta} iss = \frac{1}{1 + ⅇ^{\frac{v + 43.94}{5.711}}} delta_epi = \{\begin{matrix} 1 - (\frac{0.95}{1 + ⅇ^{\frac{v + 70}{5}}}) if celltype = 1 \\ 1 otherwise \end{matrix} tiF_b = 4.562 + \frac{1}{0.3933 ⅇ^{\frac{- (v + 100)}{100}} + 0.08004 ⅇ^{\frac{v + 50}{16.59}}} tiS_b = 23.62 + \frac{1}{0.001416 ⅇ^{\frac{- (v + 96.52)}{59.05}} + 1.78 e -8 ⅇ^{\frac{v + 114.1}{8.079}}} tiF = tiF_b delta_epi tiS = tiS_b delta_epi AiF = \frac{1}{1 + ⅇ^{\frac{v - 213.6}{151.2}}} AiS = 1 - AiF \frac{d}{d time} (iF) = \frac{iss - iF}{tiF} \frac{d}{d time} (iS) = \frac{iss - iS}{tiS} i = AiF iF + AiS iS assp = \frac{1}{1 + ⅇ^{\frac{- (v - 24.34)}{14.82}}} \frac{d}{d time} (ap) = \frac{assp - ap}{ta} dti_develop = 1.354 + \frac{1 e -4}{ⅇ^{\frac{v - 167.4}{15.89}} + ⅇ^{\frac{- (v - 12.23)}{0.2154}}} dti_recover = 1 - (\frac{0.5}{1 + ⅇ^{\frac{v + 70}{20}}}) tiFp = dti_develop dti_recover tiF tiSp = dti_develop dti_recover tiS \frac{d}{d time} (iFp) = \frac{iss - iFp}{tiFp} \frac{d}{d time} (iSp) = \frac{iss - iSp}{tiSp} ip = AiF iFp + AiS iSp Gto = \{\begin{matrix} Gto_b 4 if celltype = 1 \\ Gto_b 4 if celltype = 2 \\ Gto_b otherwise \end{matrix} fItop = \frac{1}{1 + \frac{KmCaMK}{CaMKa}} Ito = Gto (v - EK) ((1 - fItop) a i + fItop ap ip)

Component: ICaL

dss = \frac{1}{1 + ⅇ^{\frac{- (v + 3.94)}{4.23}}} td = 0.6 + \frac{1}{ⅇ^{(- 0.05) (v + 6)} + ⅇ^{0.09 (v + 14)}} \frac{d}{d time} (d) = \frac{dss - d}{td} fss = \frac{1}{1 + ⅇ^{\frac{v + 19.58}{3.696}}} tff = 7 + \frac{1}{0.0045 ⅇ^{\frac{- (v + 20)}{10}} + 0.0045 ⅇ^{\frac{v + 20}{10}}} tfs = 1000 + \frac{1}{0.000035 ⅇ^{\frac{- (v + 5)}{4}} + 0.000035 ⅇ^{\frac{v + 5}{6}}} Aff = 0.6 Afs = 1 - Aff \frac{d}{d time} (ff) = \frac{fss - ff}{tff} \frac{d}{d time} (fs) = \frac{fss - fs}{tfs} f = Aff ff + Afs fs fcass = fss tfcaf = 7 + \frac{1}{0.04 ⅇ^{\frac{- (v - 4)}{7}} + 0.04 ⅇ^{\frac{v - 4}{7}}} tfcas = 100 + \frac{1}{0.00012 ⅇ^{\frac{- v}{3}} + 0.00012 ⅇ^{\frac{v}{7}}} Afcaf = 0.3 + \frac{0.6}{1 + ⅇ^{\frac{v - 10}{10}}} Afcas = 1 - Afcaf \frac{d}{d time} (fcaf) = \frac{fcass - fcaf}{tfcaf} \frac{d}{d time} (fcas) = \frac{fcass - fcas}{tfcas} fca = Afcaf fcaf + Afcas fcas tjca = 75 \frac{d}{d time} (jca) = \frac{fcass - jca}{tjca} tffp = 2.5 tff \frac{d}{d time} (ffp) = \frac{fss - ffp}{tffp} fp = Aff ffp + Afs fs tfcafp = 2.5 tfcaf \frac{d}{d time} (fcafp) = \frac{fcass - fcafp}{tfcafp} fcap = Afcaf fcafp + Afcas fcas km2n = jca 1 anca = \frac{1}{\frac{k2n}{km2n} + {(1 + \frac{Kmn}{cass})}^{4}} \frac{d}{d time} (nca) = anca k2n - (nca km2n)

v0 = 0

B_1 = 2 frt A_1 = \frac{4 ffrt (cass ⅇ^{2 vfrt} - (0.341 cao))}{B_1} U_1 = B_1 (v - v0) PhiCaL = \{\begin{matrix} A_1 (1 - (0.5 U_1)) if - 1 e -7 ≦ U_1 \land U_1 ≦ 1 e -7 \\ \frac{A_1 U_1}{ⅇ^{U_1} - 1} otherwise \end{matrix} B_2 = frt A_2 = \frac{0.75 ffrt (nass ⅇ^{vfrt} - nao)}{B_2} U_2 = B_2 (v - v0) PhiCaNa = \{\begin{matrix} A_2 (1 - (0.5 U_2)) if - 1 e -7 ≦ U_2 \land U_2 ≦ 1 e -7 \\ \frac{A_2 U_2}{ⅇ^{U_2} - 1} otherwise \end{matrix} B_3 = frt A_3 = \frac{0.75 ffrt (kss ⅇ^{vfrt} - ko)}{B_3} U_3 = B_3 (v - v0) PhiCaK = \{\begin{matrix} A_3 (1 - (0.5 U_3)) if - 1 e -7 ≦ U_3 \land U_3 ≦ 1 e -7 \\ \frac{A_3 U_3}{ⅇ^{U_3} - 1} otherwise \end{matrix} PCa = \{\begin{matrix} PCa_b 1.2 if celltype = 1 \\ PCa_b 2.5 if celltype = 2 \\ PCa_b otherwise \end{matrix} PCap = 1.1 PCa PCaNa = 0.00125 PCa PCaK = 3.574 e -4 PCa PCaNap = 0.00125 PCap PCaKp = 3.574 e -4 PCap fICaLp = \frac{1}{1 + \frac{KmCaMK}{CaMKa}} ICaL = (1 - fICaLp) PCa PhiCaL d (f (1 - nca) + jca fca nca) + fICaLp PCap PhiCaL d (fp (1 - nca) + jca fcap nca) ICaNa = (1 - fICaLp) PCaNa PhiCaNa d (f (1 - nca) + jca fca nca) + fICaLp PCaNap PhiCaNa d (fp (1 - nca) + jca fcap nca) ICaK = (1 - fICaLp) PCaK PhiCaK d (f (1 - nca) + jca fca nca) + fICaLp PCaKp PhiCaK d (fp (1 - nca) + jca fcap nca)

Component: IKr

\frac{d}{d time} (IC1) = (- (A11 ⅇ^{B11 v} IC1 ⅇ^{\frac{(Temp - 20) \ln q11}{10}} - (A21 ⅇ^{B21 v} IC2 ⅇ^{\frac{(Temp - 20) \ln q21}{10}}))) + A51 ⅇ^{B51 v} C1 ⅇ^{\frac{(Temp - 20) \ln q51}{10}} - (A61 ⅇ^{B61 v} IC1 ⅇ^{\frac{(Temp - 20) \ln q61}{10}}) \frac{d}{d time} (IC2) = (A11 ⅇ^{B11 v} IC1 ⅇ^{\frac{(Temp - 20) \ln q11}{10}} - (A21 ⅇ^{B21 v} IC2 ⅇ^{\frac{(Temp - 20) \ln q21}{10}}) - (A3 ⅇ^{B3 v} IC2 ⅇ^{\frac{(Temp - 20) \ln q3}{10}} - (A4 ⅇ^{B4 v} IO ⅇ^{\frac{(Temp - 20) \ln q4}{10}}))) + A52 ⅇ^{B52 v} C2 ⅇ^{\frac{(Temp - 20) \ln q52}{10}} - (A62 ⅇ^{B62 v} IC2 ⅇ^{\frac{(Temp - 20) \ln q62}{10}}) \frac{d}{d time} (C1) = - (A1 ⅇ^{B1 v} C1 ⅇ^{\frac{(Temp - 20) \ln q1}{10}} - (A2 ⅇ^{B2 v} C2 ⅇ^{\frac{(Temp - 20) \ln q2}{10}})) - (A51 ⅇ^{B51 v} C1 ⅇ^{\frac{(Temp - 20) \ln q51}{10}} - (A61 ⅇ^{B61 v} IC1 ⅇ^{\frac{(Temp - 20) \ln q61}{10}})) \frac{d}{d time} (C2) = A1 ⅇ^{B1 v} C1 ⅇ^{\frac{(Temp - 20) \ln q1}{10}} - (A2 ⅇ^{B2 v} C2 ⅇ^{\frac{(Temp - 20) \ln q2}{10}}) - (A31 ⅇ^{B31 v} C2 ⅇ^{\frac{(Temp - 20) \ln q31}{10}} - (A41 ⅇ^{B41 v} O ⅇ^{\frac{(Temp - 20) \ln q41}{10}})) - (A52 ⅇ^{B52 v} C2 ⅇ^{\frac{(Temp - 20) \ln q52}{10}} - (A62 ⅇ^{B62 v} IC2 ⅇ^{\frac{(Temp - 20) \ln q62}{10}})) \frac{d}{d time} (O) = A31 ⅇ^{B31 v} C2 ⅇ^{\frac{(Temp - 20) \ln q31}{10}} - (A41 ⅇ^{B41 v} O ⅇ^{\frac{(Temp - 20) \ln q41}{10}}) - (A53 ⅇ^{B53 v} O ⅇ^{\frac{(Temp - 20) \ln q53}{10}} - (A63 ⅇ^{B63 v} IO ⅇ^{\frac{(Temp - 20) \ln q63}{10}})) - (\frac{Kmax Ku ⅇ^{n \ln D}}{ⅇ^{n \ln D} + halfmax} O - (Ku Obound)) \frac{d}{d time} (IO) = (A3 ⅇ^{B3 v} IC2 ⅇ^{\frac{(Temp - 20) \ln q3}{10}} - (A4 ⅇ^{B4 v} IO ⅇ^{\frac{(Temp - 20) \ln q4}{10}})) + A53 ⅇ^{B53 v} O ⅇ^{\frac{(Temp - 20) \ln q53}{10}} - (A63 ⅇ^{B63 v} IO ⅇ^{\frac{(Temp - 20) \ln q63}{10}}) - (\frac{Kmax Ku ⅇ^{n \ln D}}{ⅇ^{n \ln D} + halfmax} IO - (\frac{Ku A53 ⅇ^{B53 v} ⅇ^{\frac{(Temp - 20) \ln q53}{10}}}{A63 ⅇ^{B63 v} ⅇ^{\frac{(Temp - 20) \ln q63}{10}}} IObound)) \frac{d}{d time} (IObound) = (\frac{Kmax Ku ⅇ^{n \ln D}}{ⅇ^{n \ln D} + halfmax} IO - (\frac{Ku A53 ⅇ^{B53 v} ⅇ^{\frac{(Temp - 20) \ln q53}{10}}}{A63 ⅇ^{B63 v} ⅇ^{\frac{(Temp - 20) \ln q63}{10}}} IObound)) + \frac{Kt}{1 + ⅇ^{\frac{- (v - Vhalf)}{6.789}}} Cbound - (Kt IObound) \frac{d}{d time} (Obound) = (\frac{Kmax Ku ⅇ^{n \ln D}}{ⅇ^{n \ln D} + halfmax} O - (Ku Obound)) + \frac{Kt}{1 + ⅇ^{\frac{- (v - Vhalf)}{6.789}}} Cbound - (Kt Obound) \frac{d}{d time} (Cbound) = - (\frac{Kt}{1 + ⅇ^{\frac{- (v - Vhalf)}{6.789}}} Cbound - (Kt Obound)) - (\frac{Kt}{1 + ⅇ^{\frac{- (v - Vhalf)}{6.789}}} Cbound - (Kt IObound)) \frac{d}{d time} (D) = 0 GKr = \{\begin{matrix} GKr_b 1.3 if celltype = 1 \\ GKr_b 0.8 if celltype = 2 \\ GKr_b otherwise \end{matrix} IKr = GKr \sqrt{\frac{ko}{5.4}} O (v - EK)

Component: IKs

xs1ss = \frac{1}{1 + ⅇ^{\frac{- (v + 11.6)}{8.932}}} txs1 = txs1_max + \frac{1}{2.326 e -4 ⅇ^{\frac{v + 48.28}{17.8}} + 0.001292 ⅇ^{\frac{- (v + 210)}{230}}} \frac{d}{d time} (xs1) = \frac{xs1ss - xs1}{txs1} xs2ss = xs1ss txs2 = \frac{1}{0.01 ⅇ^{\frac{v - 50}{20}} + 0.0193 ⅇ^{\frac{- (v + 66.54)}{31}}} \frac{d}{d time} (xs2) = \frac{xs2ss - xs2}{txs2} KsCa = 1 + \frac{0.6}{1 + {(\frac{3.8 e -5}{cai})}^{1.4}} GKs = \{\begin{matrix} GKs_b 1.4 if celltype = 1 \\ GKs_b otherwise \end{matrix} IKs = GKs KsCa xs1 xs2 (v - EKs)

Component: IK1

xk1ss = \frac{1}{1 + ⅇ^{\frac{- (v + 2.5538 ko + 144.59)}{1.5692 ko + 3.8115}}} txk1 = \frac{122.2}{ⅇ^{\frac{- (v + 127.2)}{20.36}} + ⅇ^{\frac{v + 236.8}{69.33}}} \frac{d}{d time} (xk1) = \frac{xk1ss - xk1}{txk1} rk1 = \frac{1}{1 + ⅇ^{\frac{v + 105.8 - (2.6 ko)}{9.493}}} GK1 = \{\begin{matrix} GK1_b 1.2 if celltype = 1 \\ GK1_b 1.3 if celltype = 2 \\ GK1_b otherwise \end{matrix} IK1 = GK1 \sqrt{ko} rk1 xk1 (v - EK)

Component: INaCa_i

hca = ⅇ^{\frac{qca v F}{R T}} hna = ⅇ^{\frac{qna v F}{R T}} h1_i = 1 + \frac{nai}{kna3} (1 + hna) h2_i = \frac{nai hna}{kna3 h1_i} h3_i = \frac{1}{h1_i} h4_i = 1 + \frac{nai}{kna1} (1 + \frac{nai}{kna2}) h5_i = \frac{nai nai}{h4_i kna1 kna2} h6_i = \frac{1}{h4_i} h7_i = 1 + \frac{nao}{kna3} (1 + \frac{1}{hna}) h8_i = \frac{nao}{kna3 hna h7_i} h9_i = \frac{1}{h7_i} h10_i = kasymm + 1 + \frac{nao}{kna1} (1 + \frac{nao}{kna2}) h11_i = \frac{nao nao}{h10_i kna1 kna2} h12_i = \frac{1}{h10_i} k1_i = h12_i cao kcaon k2_i = kcaoff k3p_i = h9_i wca k3pp_i = h8_i wnaca k3_i = k3p_i + k3pp_i k4p_i = \frac{h3_i wca}{hca} k4pp_i = h2_i wnaca k4_i = k4p_i + k4pp_i k5_i = kcaoff k6_i = h6_i cai kcaon k7_i = h5_i h2_i wna k8_i = h8_i h11_i wna x1_i = k2_i k4_i (k7_i + k6_i) + k5_i k7_i (k2_i + k3_i) x2_i = k1_i k7_i (k4_i + k5_i) + k4_i k6_i (k1_i + k8_i) x3_i = k1_i k3_i (k7_i + k6_i) + k8_i k6_i (k2_i + k3_i) x4_i = k2_i k8_i (k4_i + k5_i) + k3_i k5_i (k1_i + k8_i) E1_i = \frac{x1_i}{x1_i + x2_i + x3_i + x4_i} E2_i = \frac{x2_i}{x1_i + x2_i + x3_i + x4_i} E3_i = \frac{x3_i}{x1_i + x2_i + x3_i + x4_i} E4_i = \frac{x4_i}{x1_i + x2_i + x3_i + x4_i} allo_i = \frac{1}{1 + {(\frac{KmCaAct}{cai})}^{2}} JncxNa_i = 3 (E4_i k7_i - (E1_i k8_i)) + E3_i k4pp_i - (E2_i k3pp_i) JncxCa_i = E2_i k2_i - (E1_i k1_i) Gncx = \{\begin{matrix} Gncx_b 1.1 if celltype = 1 \\ Gncx_b 1.4 if celltype = 2 \\ Gncx_b otherwise \end{matrix} INaCa_i = 0.8 Gncx allo_i (zna JncxNa_i + zca JncxCa_i) h1_ss = 1 + \frac{nass}{kna3} (1 + hna) h2_ss = \frac{nass hna}{kna3 h1_ss} h3_ss = \frac{1}{h1_ss} h4_ss = 1 + \frac{nass}{kna1} (1 + \frac{nass}{kna2}) h5_ss = \frac{nass nass}{h4_ss kna1 kna2} h6_ss = \frac{1}{h4_ss} h7_ss = 1 + \frac{nao}{kna3} (1 + \frac{1}{hna}) h8_ss = \frac{nao}{kna3 hna h7_ss} h9_ss = \frac{1}{h7_ss} h10_ss = kasymm + 1 + \frac{nao}{kna1} (1 + \frac{nao}{kna2}) h11_ss = \frac{nao nao}{h10_ss kna1 kna2} h12_ss = \frac{1}{h10_ss} k1_ss = h12_ss cao kcaon k2_ss = kcaoff k3p_ss = h9_ss wca k3pp_ss = h8_ss wnaca k3_ss = k3p_ss + k3pp_ss k4p_ss = \frac{h3_ss wca}{hca} k4pp_ss = h2_ss wnaca k4_ss = k4p_ss + k4pp_ss k5_ss = kcaoff k6_ss = h6_ss cass kcaon k7_ss = h5_ss h2_ss wna k8_ss = h8_ss h11_ss wna x1_ss = k2_ss k4_ss (k7_ss + k6_ss) + k5_ss k7_ss (k2_ss + k3_ss) x2_ss = k1_ss k7_ss (k4_ss + k5_ss) + k4_ss k6_ss (k1_ss + k8_ss) x3_ss = k1_ss k3_ss (k7_ss + k6_ss) + k8_ss k6_ss (k2_ss + k3_ss) x4_ss = k2_ss k8_ss (k4_ss + k5_ss) + k3_ss k5_ss (k1_ss + k8_ss) E1_ss = \frac{x1_ss}{x1_ss + x2_ss + x3_ss + x4_ss} E2_ss = \frac{x2_ss}{x1_ss + x2_ss + x3_ss + x4_ss} E3_ss = \frac{x3_ss}{x1_ss + x2_ss + x3_ss + x4_ss} E4_ss = \frac{x4_ss}{x1_ss + x2_ss + x3_ss + x4_ss} allo_ss = \frac{1}{1 + {(\frac{KmCaAct}{cass})}^{2}} JncxNa_ss = 3 (E4_ss k7_ss - (E1_ss k8_ss)) + E3_ss k4pp_ss - (E2_ss k3pp_ss) JncxCa_ss = E2_ss k2_ss - (E1_ss k1_ss) INaCa_ss = 0.2 Gncx allo_ss (zna JncxNa_ss + zca JncxCa_ss)

Component: INaK

Knai = Knai0 ⅇ^{\frac{delta v F}{3 R T}} Knao = Knao0 ⅇ^{\frac{(1 - delta) v F}{3 R T}} P = \frac{eP}{1 + \frac{H}{Khp} + \frac{nai}{Knap} + \frac{ki}{Kxkur}} a1 = \frac{k1p {(\frac{nai}{Knai})}^{3}}{{(1 + \frac{nai}{Knai})}^{3} + {(1 + \frac{ki}{Kki})}^{2} - 1} b1 = k1m MgADP a2 = k2p b2 = \frac{k2m {(\frac{nao}{Knao})}^{3}}{{(1 + \frac{nao}{Knao})}^{3} + {(1 + \frac{ko}{Kko})}^{2} - 1} a3 = \frac{k3p {(\frac{ko}{Kko})}^{2}}{{(1 + \frac{nao}{Knao})}^{3} + {(1 + \frac{ko}{Kko})}^{2} - 1} b3 = \frac{k3m P H}{1 + \frac{MgATP}{Kmgatp}} a4 = \frac{\frac{k4p MgATP}{Kmgatp}}{1 + \frac{MgATP}{Kmgatp}} b4 = \frac{k4m {(\frac{ki}{Kki})}^{2}}{{(1 + \frac{nai}{Knai})}^{3} + {(1 + \frac{ki}{Kki})}^{2} - 1} x1 = a4 a1 a2 + b2 b4 b3 + a2 b4 b3 + b3 a1 a2 x2 = b2 b1 b4 + a1 a2 a3 + a3 b1 b4 + a2 a3 b4 x3 = a2 a3 a4 + b3 b2 b1 + b2 b1 a4 + a3 a4 b1 x4 = b4 b3 b2 + a3 a4 a1 + b2 a4 a1 + b3 b2 a1 E1 = \frac{x1}{x1 + x2 + x3 + x4} E2 = \frac{x2}{x1 + x2 + x3 + x4} E3 = \frac{x3}{x1 + x2 + x3 + x4} E4 = \frac{x4}{x1 + x2 + x3 + x4} JnakNa = 3 (E1 a3 - (E2 b3)) JnakK = 2 (E4 b1 - (E3 a1)) Pnak = \{\begin{matrix} Pnak_b 0.9 if celltype = 1 \\ Pnak_b 0.7 if celltype = 2 \\ Pnak_b otherwise \end{matrix} INaK = Pnak (zna JnakNa + zk JnakK)

Component: IKb

xkb = \frac{1}{1 + ⅇ^{\frac{- (v - 14.48)}{18.34}}} GKb = \{\begin{matrix} GKb_b 0.6 if celltype = 1 \\ GKb_b otherwise \end{matrix} IKb = GKb xkb (v - EK)

Component: INab

B = frt v0 = 0 A = \frac{PNab ffrt (nai ⅇ^{vfrt} - nao)}{B} U = B (v - v0) INab = \{\begin{matrix} A (1 - (0.5 U)) if - 1 e -7 ≦ U \land U ≦ 1 e -7 \\ \frac{A U}{ⅇ^{U} - 1} otherwise \end{matrix}

Component: ICab

B = 2 frt v0 = 0 A = \frac{PCab 4 ffrt (cai ⅇ^{2 vfrt} - (0.341 cao))}{B} U = B (v - v0) ICab = \{\begin{matrix} A (1 - (0.5 U)) if - 1 e -7 ≦ U \land U ≦ 1 e -7 \\ \frac{A U}{ⅇ^{U} - 1} otherwise \end{matrix}

Component: IpCa

IpCa = \frac{GpCa cai}{KmCap + cai}

Component: diff

JdiffNa = \frac{nass - nai}{2} JdiffK = \frac{kss - ki}{2} Jdiff = \frac{cass - cai}{0.2}

Component: ryr

a_rel = 0.5 bt Jrel_inf_temp = \frac{a_rel (- ICaL)}{1 + 1 {(\frac{1.5}{cajsr})}^{8}} Jrel_inf = \{\begin{matrix} Jrel_inf_temp 1.7 if celltype = 2 \\ Jrel_inf_temp otherwise \end{matrix} tau_rel_temp = \frac{bt}{1 + \frac{0.0123}{cajsr}} tau_rel = \{\begin{matrix} 0.001 if tau_rel_temp < 0.001 \\ tau_rel_temp otherwise \end{matrix} \frac{d}{d time} (Jrelnp) = \frac{Jrel_inf - Jrelnp}{tau_rel} btp = 1.25 bt a_relp = 0.5 btp Jrel_temp = \frac{a_relp (- ICaL)}{1 + {(\frac{1.5}{cajsr})}^{8}} Jrel_infp = \{\begin{matrix} Jrel_temp 1.7 if celltype = 2 \\ Jrel_temp otherwise \end{matrix} tau_relp_temp = \frac{btp}{1 + \frac{0.0123}{cajsr}} tau_relp = \{\begin{matrix} 0.001 if tau_relp_temp < 0.001 \\ tau_relp_temp otherwise \end{matrix} \frac{d}{d time} (Jrelp) = \frac{Jrel_infp - Jrelp}{tau_relp} fJrelp = \frac{1}{1 + \frac{KmCaMK}{CaMKa}} Jrel = Jrel_scaling_factor ((1 - fJrelp) Jrelnp + fJrelp Jrelp)

Component: SERCA

upScale = \{\begin{matrix} 1.3 if celltype = 1 \\ 1 otherwise \end{matrix} Jupnp = \frac{upScale 0.004375 cai}{cai + 0.00092} Jupp = \frac{upScale 2.75 0.004375 cai}{cai + 0.00092 - 0.00017} fJupp = \frac{1}{1 + \frac{KmCaMK}{CaMKa}} Jleak = \frac{0.0039375 cansr}{15} Jup = Jup_b ((1 - fJupp) Jupnp + fJupp Jupp - Jleak)

Component: trans_flux

Jtr = \frac{cansr - cajsr}{100}

Source: Derived from workspace O'Hara-Rudy-CiPA-v1.0 (2017) at changeset 2c2629a8e2b6.

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.

O'Hara-Rudy CiPA v1.0 (2017)