begin name
fribourg2
end name

begin reactions
J1 {1.0}$w = {1.0}IFNbmRNA
J10 {1.0}$w = {1.0}TNFenv
J11 {1.0}$w = {1.0}STATm
J12 {1.0}$w = {1.0}STAT
J2 {1.0}$w = {1.0}IFNbenv
J3 {1.0}$w = {1.0}STATP2n
J4 {1.0}$w = {1.0}SOCS1m
J5 {1.0}$w = {1.0}IRF7m
J6 {1.0}$w = {1.0}IRF7Pn
J7 {1.0}$w = {1.0}IFNamRNA
J8 {1.0}$w = {1.0}IFNaenv
J9 {1.0}$w = {1.0}TNFam
end reactions

begin rate equations
J1 = IC2*(IC1*r0 + k15*IRF7Pn[t]) - (IFNbmRNA[t]*Log[2])/tao1
J10 = (1000000000*C*vmax19*TNFam[t])/(NA*(K19 + TNFam[t]))
J11 = -((Log[2]*STATm[t])/tao12) + IC2*(IC1*r4 + k26*STATP2n[t])
J12 = -((Log[2]*STAT[t])/tao13) + k28*STATm[t]
J2 = (1000000000*C*vmax2*IFNbmRNA[t])/(NA*(K2 + IFNbmRNA[t]))
J3 = (K5*TJ*STAT[t])/(2*(K5 + STAT[t])) - (Log[2]*STATP2n[t])/tao3
J4 = -((Log[2]*SOCS1m[t])/tao4) + IC2*(IC1*r3 + k8*STATP2n[t])
J5 = -((IRF7m[t]*Log[2])/tao6) + IC2*(k14*IRF7Pn[t] + k11*STATP2n[t])
J6 = IC1*k12*IRF7m[t]
J7 = IC2ifa*k16*IRF7Pn[t] - (IFNamRNA[t]*Log[2])/tao8
J8 = (1000000000*C*vmax17*IFNamRNA[t])/(NA*(K17 + IFNamRNA[t]))
J9 = -((Log[2]*TNFam[t])/tao9) + IC2*(IC1*r1 + (rmax20*TNFenv[t])/(K20 + TNFenv[t]))
end rate equations

begin parameters
C = 500000.0
K17 = 0.002
K19 = 0.004
K2 = 72000.0
K20 = 0.0006
K3 = 0.0043
K5 = 0.01
K9 = 780.0
NA = "6.023e+23"
TJtot = 0.0001
bm = 4.5
delta = 0.0003
delta1 = 0.1
delta2 = 0.4
k11 = 0.00036
k12 = 3600.0
k14 = "3.204e-07"
k15 = "3.6e-05"
k16 = 3600.0
k26 = 0.360085
k28 = 360.0
k8 = 0.0036
n1 = 5.0
n2 = 5.0
n3 = 8.0
r0 = 0.003
r1 = 0.00025
r3 = "1e-07"
r4 = "1e-06"
r5 = 1.0
rmax20 = 0.001
sp = 0.1
sv = 0.1
tao1 = 2.5
tao12 = 1.0
tao13 = 15.0
tao3 = 0.56
tao4 = 0.46
tao6 = 1.0
tao8 = 2.0
tao9 = 2.0
vmax17 = 72000.0
vmax19 = 154800.0
vmax2 = 72000.0
w = 0.0
c2 = 1.0
compartment = 1.0
default = 1.0
end parameters

begin initial conditions
IFNaenv[0] = IFNaenvi
IFNamRNA[0] = IFNamRNAi
IFNbenv[0] = IFNbenvi
IFNbmRNA[0] = IFNbmRNAi
IRF7Pn[0] = IRF7Pni
IRF7m[0] = IRF7mi
SOCS1m[0] = SOCS1mi
STAT[0] = STATi
STATP2n[0] = STATP2ni
STATm[0] = STATmi
TNFam[0] = TNFami
TNFenv[0] = TNFenvi
end initial conditions

begin initial values
IFNaenvi = 0.0
IFNamRNAi = 0.0
IFNbenvi = 0.0
IFNbmRNAi = 0.0
IRF7Pni = 0.0
IRF7mi = 0.0
SOCS1mi = 0.0
STATi = 0.1
STATP2ni = 0.0
STATmi = 0.0
TNFami = 0.0
TNFenvi = 0.0
end initial values

begin assignment rules
NS := (r5*t^n3)/(bm^n3 + t^n3)
IC1 := (1 + (NS/delta1)^n1*sp)/(1 + (NS/delta1)^n1)
IC2 := (1 + (NS/delta2)^n2*sv)/(1 + (NS/delta2)^n2)
TJ := (TJtot*(IFNaenv[t] + IFNbenv[t]))/((K3 + IFNaenv[t] + IFNbenv[t])*(1 + (K9*SOCS1m[t])/delta))
IC2ifa := (1 + 3*(NS/delta2)^n2*sv)/(1 + (NS/delta2)^n2)
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units