begin name
maiwald
end name

begin reactions
v1 {0.5}ReceptorIFNAR1 + {0.5}TYK = {0.5}ReceptorTYKcomplex
v10 {0.5}ISGF3c = {1.0}ISGF3n
v11 {0.5}STAT1cstarSTAT2cstar = {1.0}STAT1nstarSTAT2nstar
v12 {1.0}STAT1nstarSTAT2nstar + {1.0}IRF9n = {1.0}ISGF3n
v13 {1.0}OpenISGF3nbindingsites + {1.0}ISGF3n = {1.0}OccupiedISGF3nbindingsites
v14 {1.0}$extVariable = {1.0}mRNAn
v15 {1.0}mRNAn = {0.5}mRNAc
v16 {1.0}$extVariable = {0.5}IRF9c
v17 {1.0}$extVariable = {0.5}SOCS1
v18 {0.5}ActivatedReceptorComplex = {0.5}ReceptorIFNAR1 + {0.5}species6 + {0.5}JAK + {0.5}TYK
v19 {0.5}ActivatedReceptorComplex = {0.5}IFNARdimer
v2 {0.5}JAK + {0.5}species6 = {0.5}ReceptorJAKcomplex
v20 {1.0}IRF9n = {1.0}$extVariable
v21 {0.5}STAT2cIRF9 = {0.5}STAT2c
v22 {1.0}STAT2nIRF9 = {1.0}STAT2n
v23 {0.5}ISGF3c + {0.5}CP = {0.5}ISGF3cCP
v24 {0.5}ISGF3cCP = {0.5}STAT1c + {0.5}STAT2c + {0.5}CP + {0.5}IRF9c
v25 {0.5}STAT1cstarSTAT2cstar + {0.5}CP = {0.5}STAT1cstarSTAT2cstarCP
v26 {0.5}STAT1cstarSTAT2cstarCP = {0.5}STAT1c + {0.5}STAT2c + {0.5}CP
v27 {1.0}NP + {1.0}STAT1nstarSTAT2nstar = {1.0}STAT1nstarSTAT2nstarNP
v28 {1.0}STAT1nstarSTAT2nstarNP = {1.0}STAT1n + {1.0}STAT2n + {1.0}NP
v29 {1.0}ISGF3n + {1.0}NP = {1.0}ISGF3nNP
v3 {0.5}ReceptorJAKcomplex + {0.5}ReceptorTYKcomplex + {0.5}IFN = {0.5}IFNARdimer
v30 {1.0}ISGF3nNP = {1.0}STAT1n + {1.0}STAT2n + {1.0}NP + {1.0}IRF9n
v31 {1.0}NP + {1.0}OccupiedISGF3nbindingsites = {1.0}OccupiedISGF3nbindingsitesNP
v32 {1.0}OccupiedISGF3nbindingsitesNP = {1.0}STAT1n + {1.0}STAT2n + {1.0}IRF9n + {1.0}OpenISGF3nbindingsites + {1.0}NP
v33 {1.0}PIAS + {1.0}ISGF3n = {1.0}PIASISGF3n
v34 {0.5}mRNAc = {1.0}$extVariable
v35 {1.0}STAT1n = {0.5}STAT1c
v36 {1.0}STAT2n = {0.5}STAT2c
v37 {0.5}STAT2c + {0.5}IRF9c = {0.5}STAT2cIRF9
v38 {1.0}STAT2n + {1.0}IRF9n = {1.0}STAT2nIRF9
v39 {0.5}STAT2cIRF9 = {1.0}STAT2nIRF9
v4 {0.5}IFNARdimer = {0.5}ActivatedReceptorComplex
v40 {0.5}IRF9c = {1.0}IRF9n
v41 {0.5}preIFN = {0.5}IFN
v5 {0.5}STAT2cIRF9 + {0.5}ActivatedReceptorComplex = {0.5}ActivatedReceptorComplexSTAT2c + {0.5}IRF9c
v6 {0.5}STAT2c + {0.5}ActivatedReceptorComplex = {0.5}ActivatedReceptorComplexSTAT2c
v7 {0.5}STAT1c + {0.5}ActivatedReceptorComplexSTAT2c = {0.5}ActivatedReceptorComplexSTAT2cSTAT1c
v8 {0.5}ActivatedReceptorComplexSTAT2cSTAT1c = {0.5}ActivatedReceptorComplex + {0.5}STAT1cstarSTAT2cstar
v9 {0.5}IRF9c + {0.5}STAT1cstarSTAT2cstar = {0.5}ISGF3c
end reactions

begin rate equations
v1 = cytoplasm*(-(k2*ReceptorTYKcomplex[t]) + k1*ReceptorIFNAR1[t]*TYK[t])
v10 = k19*ISGF3c[t] - k29*ISGF3n[t]
v11 = PhosphoTransport*STAT1cstarSTAT2cstar[t] - k210*STAT1nstarSTAT2nstar[t]
v12 = nucleus*(-(k211*ISGF3n[t]) + k111*IRF9n[t]*STAT1nstarSTAT2nstar[t])
v13 = nucleus*(-(k212*OccupiedISGF3nbindingsites[t]) + k112*ISGF3n[t]*OpenISGF3nbindingsites[t])
v14 = function4Transcription*nucleus
v15 = -(k214*mRNAc[t]) + k114*mRNAn[t]
v16 = cytoplasm*function4
v17 = cytoplasm*function6
v18 = cytoplasm*function5
v19 = cytoplasm*k118*ActivatedReceptorComplex[t]
v2 = cytoplasm*(-(k21*ReceptorJAKcomplex[t]) + k11*JAK[t]*species6[t])
v20 = IRF9deg*nucleus*IRF9n[t]
v21 = cytoplasm*IRF9deg*STAT2cIRF9[t]
v22 = IRF9deg*nucleus*STAT2nIRF9[t]
v23 = cytoplasm*(k122*CP[t]*ISGF3c[t] - k222*ISGF3cCP[t])
v24 = cytoplasm*k123*ISGF3cCP[t]
v25 = cytoplasm*(k124*CP[t]*STAT1cstarSTAT2cstar[t] - k224*STAT1cstarSTAT2cstarCP[t])
v26 = cytoplasm*k125*STAT1cstarSTAT2cstarCP[t]
v27 = nucleus*(NPfreebinding*NP[t]*STAT1nstarSTAT2nstar[t] - NPfreeunbinding*STAT1nstarSTAT2nstarNP[t])
v28 = k127*nucleus*STAT1nstarSTAT2nstarNP[t]
v29 = nucleus*(-(NPfreeunbinding*ISGF3nNP[t]) + NPfreebinding*ISGF3n[t]*NP[t])
v3 = cytoplasm*(-(k22*IFNARdimer[t]) + k12*IFN[t]*ReceptorJAKcomplex[t]*ReceptorTYKcomplex[t])
v30 = k129*nucleus*ISGF3nNP[t]
v31 = nucleus*(k130*NP[t]*OccupiedISGF3nbindingsites[t] - k230*OccupiedISGF3nbindingsitesNP[t])
v32 = k131*nucleus*OccupiedISGF3nbindingsitesNP[t]
v33 = nucleus*(k132*ISGF3n[t]*PIAS[t] - k232*PIASISGF3n[t])
v34 = cytoplasm*k133*mRNAc[t]
v35 = -(k234*STAT1c[t]) + k134*STAT1n[t]
v36 = -(k235*STAT2c[t]) + k135*STAT2n[t]
v37 = cytoplasm*(k136*IRF9c[t]*STAT2c[t] - k236*STAT2cIRF9[t])
v38 = nucleus*(k137*IRF9n[t]*STAT2n[t] - k237*STAT2nIRF9[t])
v39 = k138*STAT2cIRF9[t] - k238*STAT2nIRF9[t]
v4 = cytoplasm*k13*IFNARdimer[t]
v40 = k139*IRF9c[t] - k239*IRF9n[t]
v41 = cytoplasm*k140*preIFN[t]
v5 = cytoplasm*statrecbinding*ActivatedReceptorComplex[t]*STAT2cIRF9[t]
v6 = cytoplasm*(-(statrecunbinding*ActivatedReceptorComplexSTAT2c[t]) + statrecbinding*ActivatedReceptorComplex[t]*STAT2c[t])
v7 = cytoplasm*(-(statrecunbinding*ActivatedReceptorComplexSTAT2cSTAT1c[t]) + statrecbinding*ActivatedReceptorComplexSTAT2c[t]*STAT1c[t])
v8 = cytoplasm*k17*ActivatedReceptorComplexSTAT2cSTAT1c[t]
v9 = cytoplasm*(-(k28*ISGF3c[t]) + k18*IRF9c[t]*STAT1cstarSTAT2cstar[t])
end rate equations

begin parameters
IRF9deg = 0.0001
NPfreebinding = 0.01
NPfreeunbinding = 0.1
PhosphoTransport = 0.015
ReceptorSOCSdeg = "1.2e-05"
Receptordegradation = 0.0003
SOCSbinding = 0.02
SOCSunbinding = 0.1
STAT1Start = 1500.0
STAT2Start = 250.0
StartingJak = 26.0
StartingReceptor = 1000.0
StartingTyk = 13.0
activatedsynth = 0.045
constdeg = 0.0005
cytoplasm = 2.0
k1 = 0.1
k11 = 0.1
k111 = 0.01
k112 = 0.1
k113 = "2.5e-05"
k114 = 0.001
k118 = 0.01
k12 = 0.01
k122 = 0.001
k123 = 0.003
k124 = 0.001
k125 = 0.003
k127 = 0.1
k129 = 0.002
k13 = 0.005
k130 = 0.0001
k131 = 0.1
k132 = 0.1
k133 = 0.0005
k134 = 0.01
k135 = 0.0014
k136 = 0.01
k137 = 0.01
k138 = 0.00125
k139 = 0.02
k140 = 0.0009
k17 = 8.0
k18 = 0.1
k19 = 0.015
k2 = 0.05
k21 = 0.05
k210 = 0.0
k211 = 0.01
k212 = 0.1
k214 = 0.0
k22 = 0.01
k222 = 0.2
k224 = 0.2
k230 = 0.1
k232 = 0.1
k234 = 0.00125
k235 = "8.17e-05"
k236 = 0.01
k237 = 0.01
k238 = 0.0014
k239 = 0.005
k28 = 0.1
k29 = 0.0
nucleus = 1.0
startkin = 13.0
statrecbinding = 0.002
statrecunbinding = 4.0
vactivated = 0.045
vconst = 0.005
extVariable = 0.0
defaultcompartment = 1.0
end parameters

begin initial conditions
ActivatedReceptorComplex[0] = ActivatedReceptorComplexi
ActivatedReceptorComplexSTAT2c[0] = ActivatedReceptorComplexSTAT2ci
ActivatedReceptorComplexSTAT2cSTAT1c[0] = ActivatedReceptorComplexSTAT2cSTAT1ci
CP[0] = CPi
IFN[0] = IFNi
IFNARdimer[0] = IFNARdimeri
IRF9c[0] = IRF9ci
IRF9n[0] = IRF9ni
ISGF3c[0] = ISGF3ci
ISGF3cCP[0] = ISGF3cCPi
ISGF3n[0] = ISGF3ni
ISGF3nNP[0] = ISGF3nNPi
JAK[0] = JAKi
NP[0] = NPi
OccupiedISGF3nbindingsites[0] = OccupiedISGF3nbindingsitesi
OccupiedISGF3nbindingsitesNP[0] = OccupiedISGF3nbindingsitesNPi
OpenISGF3nbindingsites[0] = OpenISGF3nbindingsitesi
PIAS[0] = PIASi
PIASISGF3n[0] = PIASISGF3ni
ReceptorIFNAR1[0] = ReceptorIFNAR1i
ReceptorJAKcomplex[0] = ReceptorJAKcomplexi
ReceptorTYKcomplex[0] = ReceptorTYKcomplexi
SOCS1[0] = SOCS1i
STAT1c[0] = STAT1ci
STAT1cstarSTAT2cstar[0] = STAT1cstarSTAT2cstari
STAT1cstarSTAT2cstarCP[0] = STAT1cstarSTAT2cstarCPi
STAT1n[0] = STAT1ni
STAT1nstarSTAT2nstar[0] = STAT1nstarSTAT2nstari
STAT1nstarSTAT2nstarNP[0] = STAT1nstarSTAT2nstarNPi
STAT2c[0] = STAT2ci
STAT2cIRF9[0] = STAT2cIRF9i
STAT2n[0] = STAT2ni
STAT2nIRF9[0] = STAT2nIRF9i
TYK[0] = TYKi
mRNAc[0] = mRNAci
mRNAn[0] = mRNAni
preIFN[0] = preIFNi
species6[0] = species6i
end initial conditions

begin initial values
ActivatedReceptorComplexi = 0.0
ActivatedReceptorComplexSTAT2ci = 0.0
ActivatedReceptorComplexSTAT2cSTAT1ci = 0.0
CPi = 20.0
IFNi = 0.0
IFNARdimeri = 0.0
IRF9ci = 45.0
IRF9ni = 0.0
ISGF3ci = 0.0
ISGF3cCPi = 0.0
ISGF3ni = 0.0
ISGF3nNPi = 0.0
JAKi = 26.0
NPi = 40.0
OccupiedISGF3nbindingsitesi = 0.0
OccupiedISGF3nbindingsitesNPi = 0.0
OpenISGF3nbindingsitesi = 500.0
PIASi = 10.0
PIASISGF3ni = 0.0
ReceptorIFNAR1i = 1000.0
ReceptorJAKcomplexi = 0.0
ReceptorTYKcomplexi = 0.0
SOCS1i = 0.0
STAT1ci = 1500.0
STAT1cstarSTAT2cstari = 0.0
STAT1cstarSTAT2cstarCPi = 0.0
STAT1ni = 0.0
STAT1nstarSTAT2nstari = 0.0
STAT1nstarSTAT2nstarNPi = 0.0
STAT2ci = 250.0
STAT2cIRF9i = 0.0
STAT2ni = 250.0
STAT2nIRF9i = 0.0
TYKi = 13.0
mRNAci = 0.0
mRNAni = 0.0
preIFNi = 6.5
species6i = 1000.0
end initial values

begin assignment rules
DNAbinding := 2*OccupiedISGF3nbindingsites[t]
function6 := activatedsynth*mRNAc[t] - constdeg*SOCS1[t]
function4Transcription := k1*OccupiedISGF3nbindingsites[t]
function5 := Receptordegradation*ActivatedReceptorComplex[t] + ReceptorSOCSdeg*ActivatedReceptorComplex[t]*SOCS1[t]
function4 := vconst - IRF9deg*IRF9c[t] + vactivated*mRNAc[t]
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