begin name
xu
end name

begin reactions
Jg10 {1.0}R2C2 + {2.0}cAMP = {1.0}R2CcAMP2 + {1.0}Cvar
Jg11 {1.0}R2CcAMP2 + {2.0}cAMP = {1.0}R2cAMP4 + {1.0}Cvar
Jg3 {1.0}$EXTERNAL = {1.0}PKa
Jg4 {1.0}PKa = {1.0}$EXTERNAL
Jg5 {1.0}$EXTERNAL = {1.0}GPa
Jg6 {1.0}GPa = {1.0}$EXTERNAL
Jg7 {1.0}GSa = {1.0}$EXTERNAL
Jg8 {1.0}$EXTERNAL = {1.0}GSa
Jg9 {1.0}PP1 + {1.0}GPa = {1.0}PP1GPa
vF1 {1.0}Fg6p = {1.0}Facyl
vF3 {1.0}Facyl = {1.0}FTG
vF4 {1.0}FTG = {1.0}$Fffa
vF5 {1.0}$Fffa = {1.0}Facyl
vGgluc {1.0}Bglucgn = {1.0}$EXTERNAL
vIgluc {1.0}$EXTERNAL = {1.0}Bins
vL10 {1.0}aK = {1.0}malate
vL11f {1.0}oam = {1.0}malate
vL11r {1.0}malate = {1.0}oam
vL12 {1.0}malate = {1.0}oac
vL13 {1.0}oac = {1.0}pep
vL14 {1.0}citrate = {1.0}oac + {1.0}acetc
vL15 {1.0}acetc = {1.0}malonyl
vL16 {1.0}acetc + {7.0}malonyl = {1.0}palm
vL17 {1.0}palm = {1.0}palmCoA
vL18 {1.0}palmCoA = {8.0}acetm
vL19 {2.0}acetm = {1.0}ket
vL1f {1.0}gluc = {1.0}g6p
vL1r {1.0}g6p = {1.0}gluc
vL20 {1.0}alan = {1.0}pyr
vL21f {1.0}pyr + {1.0}glutamate = {1.0}alan + {1.0}aK
vL21r {1.0}alan + {1.0}aK = {1.0}pyr + {1.0}glutamate
vL22 {1.0}glutamate = {1.0}$EXTERNAL
vL2f {1.0}g6p = {1.0}glycgn
vL2r {1.0}glycgn = {1.0}g6p
vL3f {1.0}g6p = {2.0}pep
vL3r {2.0}pep = {1.0}g6p
vL4 {1.0}pep = {1.0}pyr
vL5f {1.0}pyr = {1.0}lac
vL5r {1.0}lac = {1.0}pyr
vL6 {1.0}pyr = {1.0}oam
vL7 {1.0}pyr = {1.0}acetm
vL8 {1.0}oam + {1.0}acetm = {1.0}citrate
vL9 {1.0}citrate = {1.0}aK
vS1f {1.0}Mg6p = {1.0}Mglycgn
vS1r {1.0}Mglycgn = {1.0}Mg6p
vS2 {1.0}Mg6p = {2.0}Mpyr
vS3f {1.0}Mpyr = {1.0}Mlac
vS3r {1.0}Mlac = {1.0}Mpyr
vS4f {1.0}Mpyr = {1.0}$Malan
vS4r {1.0}$Malan = {1.0}Mpyr
vc1cAMP {1.0}$EXTERNAL = {1.0}cAMP
vc2cAMP {1.0}cAMP = {1.0}$EXTERNAL
vdBalan {1.0}Balan = {1.0}$EXTERNAL
vdBffa {1.0}Bffa = {1.0}$EXTERNAL
vdBgluc {1.0}Bgluc = {1.0}$EXTERNAL
vdBglucgn {1.0}Bglucgn = {1.0}$EXTERNAL
vdBins {1.0}Bins = {1.0}$EXTERNAL
vdBket {1.0}Bket = {1.0}$EXTERNAL
vdBlac {1.0}Blac = {1.0}$EXTERNAL
vfeed {1.0}$EXTERNAL = {1.0}Bgluc
vglucgn {1.0}$EXTERNAL = {1.0}Bglucgn
vins {1.0}$EXTERNAL = {1.0}Bins
vsdket {1.0}Mket = {1.0}$EXTERNAL
vtF1 {1.0}Bgluc = {1.0}Fg6p
vtF3 {1.0}$Fffa = {1.0}Bffa
vtL1 {1.0}Bgluc = {1.0}gluc
vtL2 {1.0}Blac = {1.0}lac
vtL3 {1.0}ket = {1.0}Bket
vtL5 {1.0}Bffa = {1.0}palm
vtL6 {1.0}Balan = {1.0}alan
vtS1 {1.0}Bgluc = {1.0}Mg6p
vtS2 {1.0}Bket = {1.0}Mket
vtS3 {1.0}Mlac = {1.0}Blac
vtS4 {1.0}$Malan = {1.0}Balan
end reactions

begin rate equations
Jg10 = kgc1*cAMP[t]^2*R2C2[t] - kgc1*Cvar[t]*R2CcAMP2[t]
Jg11 = -(kgc2*Cvar[t]*R2cAMP4[t]) + kgc2*cAMP[t]^2*R2CcAMP2[t]
Jg3 = (kg3*Cvar[t]*(kt - PKa[t]))/(kmg3 + kt - PKa[t])
Jg4 = (kg4*PKa[t]*(PP1[t] + PP1GPa[t]))/(kmg4 + PKa[t])
Jg5 = (kg5*(pt - GPa[t])*PKa[t])/(kmg5s + pt - GPa[t])
Jg6 = (kg6*GPa[t]*(PP1[t] + PP1GPa[t]))/(kmg6s + GPa[t])
Jg7 = (kg7*GSa[t]*(Cvar[t] + PKa[t]))/(kmg7s + GSa[t])
Jg8 = (kg8*(st - GSa[t])*PP1[t])/(kmg8s + st - GSa[t])
Jg9 = ka*GPa[t]*PP1[t] - ka*PP1GPa[t]
vF1 = kf1*(1 + Bins[t]^ep14/(kDins^ep14 + Bins[t]^ep14))*Fg6p[t]
vF3 = kf3*(1 + Bins[t]^ep15/(kDins^ep15 + Bins[t]^ep15))*Facyl[t]
vF4 = (kDins^en12*kf4*FTG[t])/(kDins^en12 + Bins[t]^en12)
vF5 = Fffa^3*kf5
vGgluc = (k1glucgn*Bgluc[t]^ng)/(kmGlgn^ng + Bgluc[t]^ng)
vIgluc = (k1ins*Bgluc[t]^ni)/(kmIns^ni + Bgluc[t]^ni)
vL10 = kL10*aK[t]
vL11f = kL11f*(1 + Bglucgn[t]^ep5/(kDglucgn^ep5 + Bglucgn[t]^ep5))*oam[t]
vL11r = kL11r*malate[t]
vL12 = (kDins^en4*kL12*(1 + Bglucgn[t]^ep6/(kDglucgn^ep6 + Bglucgn[t]^ep6))*malate[t])/(kDins^en4 + Bins[t]^en4)
vL13 = (kDins^en5*kL13*(1 + Bglucgn[t]^ep7/(kDglucgn^ep7 + Bglucgn[t]^ep7))*oac[t])/(kDins^en5 + Bins[t]^en5)
vL14 = (kDglucgn^en7*kL14*(1 + Bins[t]^ep10/(kDins^ep10 + Bins[t]^ep10))*citrate[t])/(kDglucgn^en7 + Bglucgn[t]^en7)
vL15 = (ki5*kL15*acetc[t]*(1 + (p1*citrate[t])/(kp1 + citrate[t])))/(ki5 + palmCoA[t])
vL16 = kL16*acetc[t]*malonyl[t]^7
vL17 = (ki1*kL17*palm[t])/(ki1 + malonyl[t])
vL18 = (ki2*kL18*palmCoA[t])/(ki2 + malonyl[t])
vL19 = (kDins^en8*kL19*acetm[t]^2*(1 + cAMP[t]^ep11/(kdcAMP^ep11 + cAMP[t]^ep11)))/(kDins^en8 + Bins[t]^en8)
vL1f = (kL1f*(1 + Bins[t]^ep1/(kDins^ep1 + Bins[t]^ep1))*gluc[t])/(kmL1f + gluc[t])
vL1r = (kL1r*(1 + Bglucgn[t]^ep9/(kDglucgn^ep9 + Bglucgn[t]^ep9))*g6p[t])/(kmL1r + g6p[t])
vL20 = (kDins2^en3*kL20*alan[t])/(kDins2^en3 + Bins[t]^en3)
vL21f = (kDins^en3*kL21f*glutamate[t]*pyr[t])/((kDins^en3 + Bins[t]^en3)*(kmL21g + glutamate[t])*(kmL21p + pyr[t]))
vL21r = (kDins^en3*kL21r*acetm[t]*aK[t]*alan[t])/((kmL21k + aK[t])*(kmL21a + alan[t])*(kDins^en3 + Bins[t]^en3))
vL22 = kL22*glutamate[t]
vL2f = (kL2f*g6p[t]*GSa[t])/(kmL2f + g6p[t])
vL2r = (kL2r*glycgn[t]*GPa[t])/(kmL2r + glycgn[t])
vL3f = (kDglucgn^en1*kL3f*(1 + Bins[t]^ep2/(kDins^ep2 + Bins[t]^ep2))*g6p[t])/((kDglucgn^en1 + Bglucgn[t]^en1)*(kmL3f + g6p[t]))
vL3r = (kDins2^en6*kL3r*(1 + Bglucgn[t]^ep8/(kDglucgn^ep8 + Bglucgn[t]^ep8))*pep[t])/((kDins2^en6 + Bins[t]^en6)*(kmL3r + pep[t]))
vL4 = (kDglucgn^en2*ki13*kL4*(1 + Bins[t]^ep3/(kDins^ep3 + Bins[t]^ep3))*pep[t])/((ki13 + alan[t])*(kDglucgn^en2 + Bglucgn[t]^en2)*(kmL4 + pep[t]))
vL5f = (kL5f*pyr[t])/(kmL5f + pyr[t])
vL5r = (kL5r*lac[t])/(kmL5r + lac[t])
vL6 = (kL6*(1 + (p2*acetm[t])/(kp2 + acetm[t]))*(1 + Bglucgn[t]^ep4/(kDglucgn^ep4 + Bglucgn[t]^ep4))*pyr[t])/(kmL6 + pyr[t])
vL7 = (ki8*kL7*pyr[t])/((ki8 + acetm[t])*(kmL7 + pyr[t]))
vL8 = (ki4*kL8*acetm[t]*oam[t])/(ki4 + palmCoA[t])
vL9 = kL9*citrate[t]
vS1f = ks1f*Mg6p[t]
vS1r = ks1r*Mglycgn[t]
vS2 = ks2*Mg6p[t]
vS3f = ks3f*Mpyr[t]
vS3r = ks3r*Mlac[t]
vS4f = ks4f*Mpyr[t]
vS4r = ks4r*Malan*Mket[t]
vc1cAMP = (kc1*Bglucgn[t]^ng)/(kcm1^ng + Bglucgn[t]^ng)
vc2cAMP = (kc2*Bins[t]^ni*cAMP[t])/(kcm2^ni + Bins[t]^ni)
vdBalan = kdBalan*Balan[t]
vdBffa = kdBffa*Bffa[t]
vdBgluc = kdBgluc*Bgluc[t]
vdBglucgn = kdBglucgn*Bglucgn[t]
vdBins = kdBins*Bins[t]
vdBket = kdBket*Bket[t]
vdBlac = kdBlac*Blac[t]
vfeed = Piecewise[{{kf, t < 20000}, {kresume, t > tresume}}, (1.*kf)/E^(0.0006157521601035995*kf^2*(-20000 + t)^2)]
vglucgn = kglucgn
vins = kins
vsdket = ksdket*Mket[t]
vtF1 = ktF1*Bgluc[t]*(1 + Bins[t]^ep12/(kDins^ep12 + Bins[t]^ep12))
vtF3 = (Fffa^8*kDins2^en9*ktF3)/(kDins2^en9 + Bins[t]^en9)
vtL1 = ktL1*(Bgluc[t] - gluc[t])
vtL2 = ktL2*Blac[t]
vtL3 = ktL3*ket[t]
vtL5 = ktL5*Bffa[t]
vtL6 = (kDins2^en10*ktL6*Balan[t])/(kDins2^en10 + Bins[t]^en10)
vtS1 = ktS1*Bgluc[t]*(1 + Bins[t]^ep13/(kDins^ep13 + Bins[t]^ep13))
vtS2 = ktS2*Bket[t]
vtS3 = ktS3*Mlac[t]
vtS4 = (kDins2^en11*ktS4*Malan)/(kDins2^en11 + Bins[t]^en11)
end rate equations

begin parameters
Gmax = "5e-08"
Gmin = "3e-08"
Imax = "1.3e-06"
Imin = "7e-07"
PP1t = 0.00025
c0 = 5.0
capkt = 0.00025
cmax = 0.0032
cmin = "2e-06"
en1 = 10.0
en10 = 10.0
en11 = 10.0
en12 = 10.0
en2 = 10.0
en3 = 10.0
en4 = 10.0
en5 = 10.0
en6 = 10.0
en7 = 10.0
en8 = 10.0
en9 = 20.0
ep1 = 10.0
ep10 = 10.0
ep11 = 10.0
ep12 = 10.0
ep13 = 10.0
ep14 = 10.0
ep15 = 10.0
ep2 = 10.0
ep3 = 10.0
ep4 = 10.0
ep5 = 10.0
ep6 = 10.0
ep7 = 10.0
ep8 = 10.0
ep9 = 10.0
index = 5.0
k11 = "4.3e-05"
k22 = 0.0007
kDglucgn = "4e-08"
kDins = "1e-06"
kDins2 = "7.5e-07"
kL10 = 0.1
kL11f = 0.6
kL11r = 0.01
kL12 = 0.6
kL13 = 0.5
kL14 = 0.01
kL15 = 0.01
kL16 = 0.01
kL17 = 0.01
kL18 = 0.01
kL19 = 0.01
kL1f = 3.0
kL1r = 4.0
kL20 = 0.2
kL21f = 0.001
kL21r = 0.2
kL22 = 0.01
kL2f = 200.0
kL2r = 20.0
kL3f = 0.1
kL3r = 0.3
kL4 = 2.0
kL5f = 0.0
kL5r = 0.1
kL6 = 1.0
kL7 = 1.0
kL8 = 0.1
kL9 = 0.1
ka = 60.0
kdcAMP = "3.16227e-06"
kf1 = 0.1
kf3 = 0.2
kf4 = 0.1
kf5 = 0.1
kgc1 = 60000.0
kgc2 = 60000.0
ki1 = 0.1
ki13 = 2.0
ki2 = 1.0
ki4 = 3.0
ki5 = 1.0
ki8 = 2.0
kmGlgn = 8.0
kmIns = 8.0
kmL1f = 7.7
kmL1r = 1.3
kmL21a = 21.0
kmL21g = 4.3
kmL21k = 0.22
kmL21p = 0.4
kmL2f = 0.57
kmL2r = 1.4
kmL3f = 0.01
kmL3r = 0.0034
kmL4 = 0.18
kmL5f = 0.03
kmL5r = 0.8
kmL6 = 0.22
kmL7 = 0.0204
kp1 = 0.5
kp2 = 2.0
ks1f = 0.0
ks1r = 0.0
ks2 = 0.02
ks3f = 0.01
ks3r = 0.01
ks4f = 0.07
ks4r = 0.0
ktF1 = 0.01
ktF3 = 0.008
ktL1 = 100.0
ktL2 = 0.1
ktL3 = 0.1
ktL5 = 0.1
ktL6 = 1.0
ktS1 = 0.01
ktS2 = 1.0
ktS3 = 0.01
ktS4 = 3.0
kc1 = 1.0
kcm1 = "4e-08"
kcm2 = "1e-06"
kdBalan = 0.015
kdBffa = 0.015
kdBgluc = 0.015
kdBglucgn = 0.015
kdBins = 0.015
kdBket = 0.015
kdBlac = 0.015
kf = 0.5
kg2 = 0.5
kg3 = 1200.0
kg4 = 300.0
kg5 = 1200.0
kg6 = 300.0
kg7 = 1200.0
kg8 = 300.0
kgi = 10.0
kmg3 = 0.0004
kmg4 = 0.0011
kmg5 = 0.01
kmg6 = 0.005
kmg7 = 0.015
kmg8 = 0.00012
kresume = 0.6
ksdket = 0.01
kt = 0.0025
ng = 10.0
ni = 10.0
p1 = 1.0
p2 = 1.0
pt = 0.07
s1 = 100.0
s2 = 0.001
st = 0.003
tresume = 20250.0
Fffa = 2
Malan = 0.4
PKb = 0.0025
blood = 1.0
fat = 1.0
glycogenregulation = 1.0
liver = 1.0
muscle = 1.0
end parameters

begin initial conditions
Balan[0] = Balani
Bffa[0] = Bffai
Bgluc[0] = Bgluci
Bglucgn[0] = Bglucgni
Bins[0] = Binsi
Bket[0] = Bketi
Blac[0] = Blaci
Cvar[0] = Cvari
FTG[0] = FTGi
Facyl[0] = Facyli
Fg6p[0] = Fg6pi
GPa[0] = GPai
GSa[0] = GSai
Mg6p[0] = Mg6pi
Mglycgn[0] = Mglycgni
Mket[0] = Mketi
Mlac[0] = Mlaci
Mpyr[0] = Mpyri
PKa[0] = PKai
PP1[0] = PP1i
PP1GPa[0] = PP1GPai
R2C2[0] = R2C2i
R2CcAMP2[0] = R2CcAMP2i
R2cAMP4[0] = R2cAMP4i
aK[0] = aKi
acetc[0] = acetci
acetm[0] = acetmi
alan[0] = alani
cAMP[0] = cAMPi
citrate[0] = citratei
g6p[0] = g6pi
gluc[0] = gluci
glutamate[0] = glutamatei
glycgn[0] = glycgni
ket[0] = keti
lac[0] = laci
malate[0] = malatei
malonyl[0] = malonyli
oac[0] = oaci
oam[0] = oami
palm[0] = palmi
palmCoA[0] = palmCoAi
pep[0] = pepi
pyr[0] = pyri
end initial conditions

begin initial values
Balani = 0.1
Bffai = 0.3
Bgluci = 4
Bglucgni = 0
Binsi = 0
Bketi = 0.08
Blaci = 0.2
Cvari = 0
FTGi = 0
Facyli = 0
Fg6pi = 0
GPai = 0.035
GSai = 0.0015
Mg6pi = 8
Mglycgni = 100
Mketi = 0.1
Mlaci = 0.2
Mpyri = 0.2
PKai = 0
PP1i = 0.00025
PP1GPai = 0
R2C2i = 0.00025
R2CcAMP2i = 0
R2cAMP4i = 0
aKi = 1
acetci = 6
acetmi = 7
alani = 0.1
cAMPi = 0
citratei = 0.5
g6pi = 0.2
gluci = 5
glutamatei = 0.1
glycgni = 10
keti = 0.08
laci = 1
malatei = 3
malonyli = 0.3
oaci = 0.0004
oami = 0.0002
palmi = 0.3
palmCoAi = 0.3
pepi = 0.001
pyri = 3
end initial values

begin assignment rules
k1ins := (Imax - Imin)*kdBins
kgc1 := kgc1/k11
kgc2 := kgc2/k22
kmg5s := kmg5*(1 + (s1*g6p[t])/kg2)
kmg6s := kmg6/(1 + (s2*gluc[t])/kgi)
kmg7s := kmg7*(1 + (s1*g6p[t])/kg2)
kmg8s := kmg8/(1 + (s1*g6p[t])/kg2)
kc2 := 316227.7660168379*kc1
ka := ka/kd
kglucgn := Gmax*kdBglucgn
k1glucgn := (Gmax - Gmin)*kdBglucgn
kd := cmin + (c0^index*(cmax - cmin))/(c0^index + glycgn[t]^index)
kins := Imin*kdBins
GPaPlot := 14.285714285714285*GPa[t]
GSaPlot := 333.3333333333333*GSa[t]
end assignment rules

begin function definitions
end function definitions

begin events
unnamed = (t == 20000, {glycgn[t] -> 120}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units