Model Mathematics

Component: environment

Component: Gi

ddtime Gi = Jglut - Jgk

Component: Ge

ddtime Ge = G_infinity - Ge tau_G G_infinity = Gmin + Gmax - Gmin 1.0 +ⅇ I - I_ SG

Component: G6P

ddtime G6P = kappa ⁢ Jgk - JPFK

Component: FBP

ddtime FBP = kappa ⁢ JPFK - 0.5 ⁢ JGPDH

Component: F6P

F6P = 0.3 ⁢ G6P

Component: Jglut

Jglut = Vglut ⁢ Ge - Gi ⁢ Kglut Kglut + Ge ⁢ Kglut - Gi

Component: Jgk

Jgk = Vgk ⁢ Gi ngk Kgk ngk + Gi ngk

Component: JGPDH

JGPDH = 0.2 ⁢ FBP 1.0 0.5

Component: JPFK

JPFK = Vmax ⁢ 1.0 - lambda ⁢ w1110 + lambda ⁢ w0010 + w1010 + w0011 + w1110 + w0111 + w1011 + w1101 + w1111 w0000 + w1000 + w0100 + w0010 + w0001 + w1100 + w1010 + w1001 + w0101 + w0011 + w0110 + w1110 + w0111 + w1011 + w1101 + w1111

Component: w

w0000 = 1.0 w1000 = AMP K1 1.0 w0100 = FBP K2 1.0 w0010 = F6P 2.0 K3 1.0 w0001 = ATP 2.0 K4 1.0 w1100 = AMP K1 1.0 ⁢ FBP K2 1.0 w1010 = 1.0 f13 1.0 ⁢ 1.0 ⁢ AMP K1 1.0 ⁢ F6P 2.0 K3 1.0 w1001 = 1.0 f41 1.0 ⁢ 1.0 ⁢ AMP K1 1.0 ⁢ ATP 2.0 K4 1.0 w0101 = 1.0 f42 1.0 ⁢ 1.0 ⁢ FBP K2 1.0 ⁢ ATP 2.0 K4 1.0 w0011 = 1.0 f43 1.0 ⁢ 1.0 ⁢ F6P 2.0 K3 1.0 ⁢ ATP 2.0 K4 1.0 w0110 = 1.0 f23 1.0 ⁢ 1.0 ⁢ FBP K2 1.0 ⁢ F6P 2.0 K3 1.0 w1110 = 1.0 f13 1.0 ⁢ 1.0 ⁢ f23 1.0 ⁢ 1.0 ⁢ AMP K1 1.0 ⁢ FBP K2 1.0 ⁢ F6P 2.0 K3 1.0 w0111 = 1.0 f23 1.0 ⁢ 1.0 ⁢ f42 1.0 ⁢ 1.0 ⁢ f43 1.0 ⁢ 1.0 ⁢ FBP K2 1.0 ⁢ F6P 2.0 K3 1.0 ⁢ ATP 2.0 K4 1.0 w1011 = 1.0 f13 1.0 ⁢ 1.0 ⁢ f41 1.0 ⁢ 1.0 ⁢ f43 1.0 ⁢ 1.0 ⁢ AMP K1 1.0 ⁢ F6P 2.0 K3 1.0 ⁢ ATP 2.0 K4 1.0 w1101 = 1.0 f41 1.0 ⁢ 1.0 ⁢ f42 1.0 ⁢ 1.0 ⁢ AMP K1 1.0 ⁢ FBP K2 1.0 ⁢ ATP 2.0 K4 1.0 w1111 = 1.0 f13 1.0 ⁢ 1.0 ⁢ f23 1.0 ⁢ 1.0 ⁢ f41 1.0 ⁢ 1.0 ⁢ f42 1.0 ⁢ 1.0 ⁢ f43 1.0 ⁢ 1.0 ⁢ AMP K1 1.0 ⁢ FBP K2 1.0 ⁢ F6P 2.0 K3 1.0 ⁢ ATP 2.0 K4 1.0

Component: ATP

ATP = Atot - ADP + AMP

Component: AMP

AMP = ADP 2.0 ATP

Component: ADP

ddtime ADP = 1.0 tau_a ⁢ ATP - ADP ⁢ⅇ r + gamma ⁢ 1.0 - Ca r1

Component: gamma

gamma = v_gamma ⁢ JGPDH k_gamma + JGPDH

Component: membrane

ddtime v =- I_K + I_Ca + I_K_Ca + I_K_ATP cm

Component: I_K

I_K = gK ⁢ n ⁢ v - vK

Component: n

ddtime n = n_infinity - n tau_n tau_n = 1.0 0.035 ⁢cosh⁡ v + 16.0 22.4 n_infinity = 1.0 1.0 +ⅇ- v + 16.0 5.6

Component: I_Ca

I_Ca = gCa ⁢ m_infinity ⁢ v - vCa

Component: m

m_infinity = 1.0 1.0 +ⅇ- v + 20.0 12.0

Component: I_K_Ca

I_K_Ca = gkCa ⁢ v - vK gkCa = gkCa_ ⁢ Ca 2.0 KD 2.0 + Ca 2.0

Component: I_K_ATP

I_K_ATP = gkATP ⁢ v - vK gkATP = gkATP_ ⁢ o_infinity o_infinity = 0.08 ⁢ 1.0 + 2.0 ⁢ MgADP 17.0 + 0.89 ⁢ MgADP 17.0 2.0 1.0 + 2.0 ⁢ MgADP 17.0 2.0 ⁢ 1.0 + ADP3_ 26.0 + ATP4_ 1.0 MgADP = 0.165 ⁢ ADP ADP3_ = 0.135 ⁢ ADP ATP4_ = 0.005 ⁢ ATP

Component: Ca

ddtime Ca = fcyt ⁢ Jmem + Jer

Component: Caer

ddtime Caer =- fer ⁢ Vcyt_Ver ⁢ Jer

Component: Jmem

Jmem =- alpha ⁢ I_Ca + kPMCA ⁢ Ca

Component: Jer

Jer = Jleak - JSERCA

Component: JSERCA

JSERCA = kSERCA ⁢ Ca

Component: Jleak

Jleak = pleak ⁢ Caer - Ca

Component: I

ddtime I = I_infinity - I tau_I I_infinity = I_slope ⁢ Ca - Canull if Ca ≥ Canull 0.0 if Ca < Canull

Component: model_parameters