#fribourg1

FIX: w

#reactions
J1:
 {1.0}w = {1.0}IFNb_mRNA
  (r0 * IC1 + k15 * IRF7Pn) * IC2 - IFNb_mRNA * log(2) / tao1
J10:
 {1.0}w = {1.0}TNFenv
  1000000000 * C * vmax19 / NA * TNFam / (K19 + TNFam)
J11:
 {1.0}w = {1.0}STATm
  (r4 * IC1 + k26 * STATP2n) * IC2 - STATm * log(2) / tao12
J12:
 {1.0}w = {1.0}STAT
  k28 * STATm - STAT * log(2) / tao13
J2:
 {1.0}w = {1.0}IFNb_env
  1000000000 * C * vmax2 / NA * IFNb_mRNA / (K2 + IFNb_mRNA)
J3:
 {1.0}w = {1.0}STATP2n
  K5 * TJ * STAT / 2 / (K5 + STAT) - STATP2n * log(2) / tao3
J4:
 {1.0}w = {1.0}SOCS1m
  (r3 * IC1 + k8 * STATP2n) * IC2 - SOCS1m * log(2) / tao4
J5:
 {1.0}w = {1.0}IRF7m
  (k11 * STATP2n + k14 * IRF7Pn) * IC2 - IRF7m * log(2) / tao6
J6:
 {1.0}w = {1.0}IRF7Pn
  k12 * IC1 * IRF7m
J7:
 {1.0}w = {1.0}IFNa_mRNA
  k16 * IRF7Pn * IC2ifa - IFNa_mRNA * log(2) / tao8
J8:
 {1.0}w = {1.0}IFNa_env
  1000000000 * C * vmax17 / NA * IFNa_mRNA / (K17 + IFNa_mRNA)
J9:
 {1.0}w = {1.0}TNFam
  (r1 * IC1 + rmax20 * TNFenv / (K20 + TNFenv)) * IC2 - TNFam * log(2) / tao9


#parameters
C = 500000.0
K17 = 0.002
K19 = 0.004
K2 = 0.002
K20 = 0.0006
K3 = 0.0043
K5 = 0.01
K9 = 780.0
NA = 6.023e+23
TJtot = 0.0001
bm = 6.0
delta = 0.0003
delta1 = 0.1
delta2 = 0.4
k11 = 0.00036
k12 = 360.0
k14 = 3.204e-07
k15 = 3.6e-08
k16 = 0.36
k26 = 0.360085
k28 = 360.0
k8 = 0.0036
n1 = 1.0
n2 = 5.0
n3 = 5.0
r0 = 0.001
r1 = 0.0001
r3 = 1e-07
r4 = 1e-05
r5 = 1.0
rmax20 = 0.001
sp = 0.3
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


#initial values
IFNa_env = 0.0
IFNa_mRNA = 0.0
IFNb_env = 0.0
IFNb_mRNA = 0.0
IRF7Pn = 0.0
IRF7m = 0.0
SOCS1m = 0.0
STAT = 0.1
STATP2n = 0.0
STATm = 0.0
TNFam = 0.0
TNFenv = 0.0


#assignment rules
!F IC1 = (1.0 + sp * pow(NS / delta1, n1)) / (1.0 + pow(NS / delta1, n1))
!F NS = r5 * pow(time, n3) / (pow(bm, n3) + pow(time, n3))
!F TJ = TJtot * (IFNb_env + IFNa_env) / (K3 + IFNb_env + IFNa_env) * 1.0 / (1.0 + K9 * SOCS1m / delta)
!F IC2ifa = 1.0
!F IC2 = 1.0