begin name
venturelli2
end name

begin reactions
v1 {1.0}$x1 = {1.0}xBH
v10 {1.0}$x10 = {1.0}xCH
v11 {1.0}$x11 = {1.0}xDP
v12 {1.0}$x12 = {1.0}xER
v2 {1.0}$x2 = {1.0}xCA
v3 {1.0}$x3 = {1.0}xBU
v4 {1.0}$x4 = {1.0}xPC
v5 {1.0}$x5 = {1.0}xBO
v6 {1.0}$x6 = {1.0}xBV
v7 {1.0}$x7 = {1.0}xBT
v8 {1.0}$x8 = {1.0}xEL
v9 {1.0}$x9 = {1.0}xFP
end reactions

begin rate equations
v1 = xBH[t]*(uBH + aBHBH*xBH[t] + aBHBO*xBO[t] + aBHBT*xBT[t] + aBHBU*xBU[t] + aBHBV*xBV[t] + aBHCA*xCA[t] + aBHCH*xCH[t] + aBHDP*xDP[t] + aBHEL*xEL[t] + aBHER*xER[t] + aBHFP*xFP[t] + aBHPC*xPC[t])
v10 = xCH[t]*(uCH + aCHBH*xBH[t] + aCHBO*xBO[t] + aCHBT*xBT[t] + aCHBU*xBU[t] + aCHBV*xBV[t] + aCHCA*xCA[t] + aCHCH*xCH[t] + aCHDP*xDP[t] + aCHEL*xEL[t] + aCHER*xER[t] + aCHFP*xFP[t] + aCHPC*xPC[t])
v11 = xDP[t]*(uDP + aDPBH*xBH[t] + aDPBO*xBO[t] + aDPBT*xBT[t] + aDPBU*xBU[t] + aDPBV*xBV[t] + aDPCA*xCA[t] + aDPCH*xCH[t] + aDPDP*xDP[t] + aDPEL*xEL[t] + aDPER*xER[t] + aDPFP*xFP[t] + aDPPC*xPC[t])
v12 = xER[t]*(uER + aERBH*xBH[t] + aERBO*xBO[t] + aERBT*xBT[t] + aERBU*xBU[t] + aERBV*xBV[t] + aERCA*xCA[t] + aERCH*xCH[t] + aERDP*xDP[t] + aEREL*xEL[t] + aERER*xER[t] + aERFP*xFP[t] + aERPC*xPC[t])
v2 = xCA[t]*(uCA + aCABH*xBH[t] + aCABO*xBO[t] + aCABT*xBT[t] + aCABU*xBU[t] + aCABV*xBV[t] + aCACA*xCA[t] + aCACH*xCH[t] + aCADP*xDP[t] + aCAEL*xEL[t] + aCAER*xER[t] + aCAFP*xFP[t] + aCAPC*xPC[t])
v3 = xBU[t]*(uBU + aBUBH*xBH[t] + aBUBO*xBO[t] + aBUBT*xBT[t] + aBUBU*xBU[t] + aBUBV*xBV[t] + aBUCA*xCA[t] + aBUCH*xCH[t] + aBUDP*xDP[t] + aBUEL*xEL[t] + aBUER*xER[t] + aBUFP*xFP[t] + aBUPC*xPC[t])
v4 = xPC[t]*(uPC + aPCBH*xBH[t] + aPCBO*xBO[t] + aPCBT*xBT[t] + aPCBU*xBU[t] + aPCBV*xBV[t] + aPCCA*xCA[t] + aPCCH*xCH[t] + aPCDP*xDP[t] + aPCEL*xEL[t] + aPCER*xER[t] + aPCFP*xFP[t] + aPCPC*xPC[t])
v5 = xBO[t]*(uBO + aBOBH*xBH[t] + aBOBO*xBO[t] + aBOBT*xBT[t] + aBOBU*xBU[t] + aBOBV*xBV[t] + aBOCA*xCA[t] + aBOCH*xCH[t] + aBODP*xDP[t] + aBOEL*xEL[t] + aBOER*xER[t] + aBOFP*xFP[t] + aBOPC*xPC[t])
v6 = xBV[t]*(uBV + aBVBH*xBH[t] + aBVBO*xBO[t] + aBVBT*xBT[t] + aBVBU*xBU[t] + aBVBV*xBV[t] + aBVCA*xCA[t] + aBVCH*xCH[t] + aBVDP*xDP[t] + aBVEL*xEL[t] + aBVER*xER[t] + aBVFP*xFP[t] + aBVPC*xPC[t])
v7 = xBT[t]*(uBT + aBTBH*xBH[t] + aBTBO*xBO[t] + aBTBT*xBT[t] + aBTBU*xBU[t] + aBTBV*xBV[t] + aBTCA*xCA[t] + aBTCH*xCH[t] + aBTDP*xDP[t] + aBTEL*xEL[t] + aBTER*xER[t] + aBTFP*xFP[t] + aBTPC*xPC[t])
v8 = xEL[t]*(uEL + aELBH*xBH[t] + aELBO*xBO[t] + aELBT*xBT[t] + aELBU*xBU[t] + aELBV*xBV[t] + aELCA*xCA[t] + aELCH*xCH[t] + aELDP*xDP[t] + aELEL*xEL[t] + aELER*xER[t] + aELFP*xFP[t] + aELPC*xPC[t])
v9 = xFP[t]*(uFP + aFPBH*xBH[t] + aFPBO*xBO[t] + aFPBT*xBT[t] + aFPBU*xBU[t] + aFPBV*xBV[t] + aFPCA*xCA[t] + aFPCH*xCH[t] + aFPDP*xDP[t] + aFPEL*xEL[t] + aFPER*xER[t] + aFPFP*xFP[t] + aFPPC*xPC[t])
end rate equations

begin parameters
aBHBH = -0.802
aBHBO = -0.1407
aBHBT = -0.1829
aBHBU = -0.1492
aBHBV = -0.0542
aBHCA = -0.1616
aBHCH = -0.0674
aBHDP = -0.5281
aBHEL = 1.1317
aBHER = 0.0
aBHFP = -0.0714
aBHPC = -0.0016
aBOBH = 0.0002
aBOBO = -0.8101
aBOBT = -0.6393
aBOBU = -0.5679
aBOBV = -0.5467
aBOCA = 0.0
aBOCH = -0.5602
aBODP = 0.0
aBOEL = 1.9238
aBOER = 0.0
aBOFP = -0.172
aBOPC = -0.0022
aBTBH = 0.0
aBTBO = -0.894
aBTBT = -1.0301
aBTBU = -1.0671
aBTBV = -0.7607
aBTCA = 0.0
aBTCH = 0.0
aBTDP = 0.0
aBTEL = 1.7918
aBTER = 0.0
aBTFP = -0.5519
aBTPC = 0.0
aBUBH = 0.0
aBUBO = -0.7768
aBUBT = -0.9798
aBUBU = -0.9673
aBUBV = -0.5769
aBUCA = 0.0
aBUCH = 0.0509
aBUDP = 0.0
aBUEL = 2.9487
aBUER = 0.0
aBUFP = -0.8667
aBUPC = 0.0
aBVBH = 0.2644
aBVBO = -0.6707
aBVBT = -0.7173
aBVBU = -0.6842
aBVBV = -0.7864
aBVCA = -0.0638
aBVCH = -0.0365
aBVDP = -0.0024
aBVEL = 1.5582
aBVER = 0.0
aBVFP = -0.7502
aBVPC = 0.0
aCABH = 0.3093
aCABO = -0.3117
aCABT = -0.2976
aCABU = -0.2897
aCABV = -0.1381
aCACA = -0.8541
aCACH = 0.2757
aCADP = 0.0
aCAEL = -0.3389
aCAER = 0.0
aCAFP = -0.9018
aCAPC = -0.552
aCHBH = 0.0
aCHBO = 0.2353
aCHBT = -0.7164
aCHBU = -0.2685
aCHBV = 0.0
aCHCA = -0.3833
aCHCH = -1.7395
aCHDP = 0.0
aCHEL = 0.0
aCHER = -0.5151
aCHFP = -0.8546
aCHPC = -0.498
aDPBH = 0.0
aDPBO = -0.1243
aDPBT = -0.132
aDPBU = -0.1096
aDPBV = -0.0218
aDPCA = 0.0902
aDPCH = -0.124
aDPDP = -0.8923
aDPEL = 0.2009
aDPER = 0.0
aDPFP = 0.0
aDPPC = 0.0
aELBH = 0.7903
aELBO = -0.0508
aELBT = -0.0554
aELBU = -0.0282
aELBV = 0.0
aELCA = -0.0021
aELCH = 0.0
aELDP = -0.0002
aELEL = -3.4228
aELER = 0.0
aELFP = 0.0
aELPC = 0.0
aERBH = 0.0003
aERBO = 0.0008
aERBT = 0.0
aERBU = 0.0
aERBV = 0.0009
aERCA = 0.0
aERCH = 0.8619
aERDP = 0.0
aEREL = 0.0
aERER = -0.7175
aERFP = 0.0
aERPC = 0.0
aFPBH = 0.8821
aFPBO = -0.1144
aFPBT = -0.0707
aFPBU = 0.1825
aFPBV = 0.6556
aFPCA = -0.0001
aFPCH = 0.3854
aFPDP = 0.9454
aFPEL = -0.7612
aFPER = 0.0
aFPFP = -0.9455
aFPPC = -0.504
aPCBH = 0.1861
aPCBO = -0.1822
aPCBT = -0.2322
aPCBU = -0.2527
aPCBV = -0.119
aPCCA = 0.4026
aPCCH = 0.3598
aPCDP = 0.0
aPCEL = -0.1694
aPCER = 0.0
aPCFP = 0.0
aPCPC = -0.8721
uBH = 0.1995
uBO = 0.5128
uBT = 0.6794
uBU = 0.6371
uBV = 0.5279
uCA = 0.2639
uCH = 0.6287
uDP = 0.1919
uEL = 0.4773
uER = 0.1066
uFP = 0.2144
uPC = 0.2144
x1 = 1.0
x10 = 1.0
x11 = 1.0
x12 = 1.0
x2 = 1.0
x3 = 1.0
x4 = 1.0
x5 = 1.0
x6 = 1.0
x7 = 1.0
x8 = 1.0
x9 = 1.0
default = 1.0
end parameters

begin initial conditions
xBH[0] = xBHi
xBO[0] = xBOi
xBT[0] = xBTi
xBU[0] = xBUi
xBV[0] = xBVi
xCA[0] = xCAi
xCH[0] = xCHi
xDP[0] = xDPi
xEL[0] = xELi
xER[0] = xERi
xFP[0] = xFPi
xPC[0] = xPCi
end initial conditions

begin initial values
xBHi = 0.01
xBOi = 0.01
xBTi = 0.01
xBUi = 0.01
xBVi = 0.01
xCAi = 0.01
xCHi = 0.01
xDPi = 0.01
xELi = 0.01
xERi = 0.01
xFPi = 0.01
xPCi = 0.01
end initial values

begin assignment rules
fDP := xDP[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fEL := xEL[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBO := xBO[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fPC := xPC[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBH := xBH[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBU := xBU[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBV := xBV[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fER := xER[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fCH := xCH[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBT := xBT[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fFP := xFP[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fCA := xCA[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
end assignment rules

begin function definitions
end function definitions

begin events
event24 = (t > 24, {xCA[t] -> 0.05*xCA[t], xBH[t] -> 0.05*xBH[t], xPC[t] -> 0.05*xPC[t], xBO[t] -> 0.05*xBO[t], xER[t] -> 0.05*xER[t], xCH[t] -> 0.05*xCH[t], xFP[t] -> 0.05*xFP[t], xEL[t] -> 0.05*xEL[t], xBV[t] -> 0.05*xBV[t], xBU[t] -> 0.05*xBU[t], xBT[t] -> 0.05*xBT[t], xDP[t] -> 0.05*xDP[t]}, 0)
event48 = (t > 48, {xDP[t] -> 0.05*xDP[t], xCH[t] -> 0.05*xCH[t], xFP[t] -> 0.05*xFP[t], xBT[t] -> 0.05*xBT[t], xBV[t] -> 0.05*xBV[t], xBO[t] -> 0.05*xBO[t], xBU[t] -> 0.05*xBU[t], xCA[t] -> 0.05*xCA[t], xBH[t] -> 0.05*xBH[t], xPC[t] -> 0.05*xPC[t], xEL[t] -> 0.05*xEL[t], xER[t] -> 0.05*xER[t]}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units