begin name
bachmann2
end name

begin reactions
reaction1 {1.0}EpoRJAK2 = {1.0}EpoRpJAK2
reaction10 {1.0}EpoRJAK2CIS = {1.0}$EXTERNAL
reaction11 {1.0}SHP1 = {1.0}SHP1Act
reaction12 {1.0}SHP1Act = {1.0}SHP1
reaction13 {1.0}STAT5 = {1.0}pSTAT5
reaction14 {1.0}STAT5 = {1.0}pSTAT5
reaction15 {1.0}pSTAT5 = {1.0}npSTAT5
reaction16 {1.0}npSTAT5 = {1.0}STAT5
reaction17 {1.0}$EXTERNAL = {1.0}CISnRNA1
reaction18 {1.0}CISnRNA1 = {1.0}CISnRNA2
reaction19 {1.0}CISnRNA2 = {1.0}CISnRNA3
reaction2 {1.0}EpoRpJAK2 = {1.0}EpoRJAK2
reaction20 {1.0}CISnRNA3 = {1.0}CISnRNA4
reaction21 {1.0}CISnRNA4 = {1.0}CISnRNA5
reaction22 {1.0}CISnRNA5 = {1.0}CISRNA
reaction23 {1.0}CISRNA = {1.0}$EXTERNAL
reaction24 {1.0}$EXTERNAL = {1.0}CIS
reaction25 {1.0}CIS = {1.0}$EXTERNAL
reaction26 {1.0}$EXTERNAL = {1.0}CIS
reaction27 {1.0}$EXTERNAL = {1.0}SOCS3nRNA1
reaction28 {1.0}SOCS3nRNA1 = {1.0}SOCS3nRNA2
reaction29 {1.0}SOCS3nRNA2 = {1.0}SOCS3nRNA3
reaction3 {1.0}EpoRpJAK2 = {1.0}p1EpoRpJAK2
reaction30 {1.0}SOCS3nRNA3 = {1.0}SOCS3nRNA4
reaction31 {1.0}SOCS3nRNA4 = {1.0}SOCS3nRNA5
reaction32 {1.0}SOCS3nRNA5 = {1.0}SOCS3RNA
reaction33 {1.0}SOCS3RNA = {1.0}$EXTERNAL
reaction34 {1.0}$EXTERNAL = {1.0}SOCS3
reaction35 {1.0}SOCS3 = {1.0}$EXTERNAL
reaction36 {1.0}$EXTERNAL = {1.0}SOCS3
reaction4 {1.0}EpoRpJAK2 = {1.0}p2EpoRpJAK2
reaction5 {1.0}p1EpoRpJAK2 = {1.0}p12EpoRpJAK2
reaction6 {1.0}p2EpoRpJAK2 = {1.0}p12EpoRpJAK2
reaction7 {1.0}p1EpoRpJAK2 = {1.0}EpoRJAK2
reaction8 {1.0}p2EpoRpJAK2 = {1.0}EpoRJAK2
reaction9 {1.0}p12EpoRpJAK2 = {1.0}EpoRJAK2
end reactions

begin rate equations
reaction1 = (cyt*JAK2ActEpo*Epo[t]*EpoRJAK2[t])/(1 + (SOCS3Inh*SOCS3[t])/SOCS3Eqc)
reaction10 = (cyt*EpoRCISRemove*EpoRJAK2CIS[t]*(p12EpoRpJAK2[t] + p1EpoRpJAK2[t]))/initEpoRJAK2
reaction11 = (cyt*SHP1ActEpoR*(EpoRpJAK2[t] + p12EpoRpJAK2[t] + p1EpoRpJAK2[t] + p2EpoRpJAK2[t])*SHP1[t])/initEpoRJAK2
reaction12 = cyt*SHP1Dea*SHP1Act[t]
reaction13 = (cyt*STAT5ActJAK2*(EpoRpJAK2[t] + p12EpoRpJAK2[t] + p1EpoRpJAK2[t] + p2EpoRpJAK2[t])*STAT5[t])/(initEpoRJAK2*(1 + (SOCS3Inh*SOCS3[t])/SOCS3Eqc))
reaction14 = (cyt*STAT5ActEpoR*(p12EpoRpJAK2[t] + p1EpoRpJAK2[t])^2*STAT5[t])/(initEpoRJAK2^2*(1 + (CISInh*CIS[t])/CISEqc)*(1 + (SOCS3Inh*SOCS3[t])/SOCS3Eqc))
reaction15 = cyt*STAT5Imp*pSTAT5[t]
reaction16 = nuc*STAT5Exp*npSTAT5[t]
reaction17 = -(((-1 + ActD)*CISRNAEqc*CISRNATurn*nuc*npSTAT5[t])/initSTAT5)
reaction18 = CISRNADelay*nuc*CISnRNA1[t]
reaction19 = CISRNADelay*nuc*CISnRNA2[t]
reaction2 = (cyt*JAK2EpoRDeaSHP1*EpoRpJAK2[t]*SHP1Act[t])/initSHP1
reaction20 = CISRNADelay*nuc*CISnRNA3[t]
reaction21 = CISRNADelay*nuc*CISnRNA4[t]
reaction22 = CISRNADelay*nuc*CISnRNA5[t]
reaction23 = CISRNATurn*cyt*CISRNA[t]
reaction24 = (CISEqc*CISTurn*cyt*CISRNA[t])/CISRNAEqc
reaction25 = CISTurn*cyt*CIS[t]
reaction26 = CISEqc*CISEqcOE*CISoe*CISTurn
reaction27 = -(((-1 + ActD)*nuc*SOCS3RNAEqc*SOCS3RNATurn*npSTAT5[t])/initSTAT5)
reaction28 = nuc*SOCS3RNADelay*SOCS3nRNA1[t]
reaction29 = nuc*SOCS3RNADelay*SOCS3nRNA2[t]
reaction3 = (cyt*EpoRActJAK2*EpoRpJAK2[t])/(1 + (SOCS3Inh*SOCS3[t])/SOCS3Eqc)
reaction30 = nuc*SOCS3RNADelay*SOCS3nRNA3[t]
reaction31 = nuc*SOCS3RNADelay*SOCS3nRNA4[t]
reaction32 = nuc*SOCS3RNADelay*SOCS3nRNA5[t]
reaction33 = cyt*SOCS3RNATurn*SOCS3RNA[t]
reaction34 = (cyt*SOCS3Eqc*SOCS3Turn*SOCS3RNA[t])/SOCS3RNAEqc
reaction35 = cyt*SOCS3Turn*SOCS3[t]
reaction36 = SOCS3Eqc*SOCS3EqcOE*SOCS3oe*SOCS3Turn
reaction4 = (3*cyt*EpoRActJAK2*EpoRpJAK2[t])/((1 + EpoRCISInh*EpoRJAK2CIS[t])*(1 + (SOCS3Inh*SOCS3[t])/SOCS3Eqc))
reaction5 = (3*cyt*EpoRActJAK2*p1EpoRpJAK2[t])/((1 + EpoRCISInh*EpoRJAK2CIS[t])*(1 + (SOCS3Inh*SOCS3[t])/SOCS3Eqc))
reaction6 = (cyt*EpoRActJAK2*p2EpoRpJAK2[t])/(1 + (SOCS3Inh*SOCS3[t])/SOCS3Eqc)
reaction7 = (cyt*JAK2EpoRDeaSHP1*p1EpoRpJAK2[t]*SHP1Act[t])/initSHP1
reaction8 = (cyt*JAK2EpoRDeaSHP1*p2EpoRpJAK2[t]*SHP1Act[t])/initSHP1
reaction9 = (cyt*JAK2EpoRDeaSHP1*p12EpoRpJAK2[t]*SHP1Act[t])/initSHP1
end rate equations

begin parameters
ActD = 0.0
CISEqc = 432.871
CISEqcOE = 0.530261
CISInh = 784653000.0
CISRNADelay = 0.144775
CISRNAEqc = 1.0
CISRNATurn = 1000.0
CISTurn = 0.00839842
CISoe = 0.0
EpoRActJAK2 = 0.267308
EpoRCISInh = 1000000.0
EpoRCISRemove = 5.42932
JAK2ActEpo = 633253.0
JAK2EpoRDeaSHP1 = 142.722
SHP1ActEpoR = 0.001
SHP1Dea = 0.00816391
SOCS3Eqc = 173.653
SOCS3EqcOE = 0.679157
SOCS3Inh = 10.408
SOCS3RNADelay = 1.06465
SOCS3RNAEqc = 1.0
SOCS3RNATurn = 0.00830844
SOCS3Turn = 10000.0
SOCS3oe = 0.0
STAT5ActEpoR = 38.9757
STAT5ActJAK2 = 0.0780965
STAT5Exp = 0.0745155
STAT5Imp = 0.0268889
epolevel = "1.24997e-07"
initEpoRJAK2 = 3.97622
initSHP1 = 26.7251
initSTAT5 = 79.7535
cyt = 0.4
nuc = 0.275
end parameters

begin initial conditions
CIS[0] = CISi
CISRNA[0] = CISRNAi
CISnRNA1[0] = CISnRNA1i
CISnRNA2[0] = CISnRNA2i
CISnRNA3[0] = CISnRNA3i
CISnRNA4[0] = CISnRNA4i
CISnRNA5[0] = CISnRNA5i
Epo[0] = Epoi
EpoRJAK2[0] = EpoRJAK2i
EpoRJAK2CIS[0] = EpoRJAK2CISi
EpoRpJAK2[0] = EpoRpJAK2i
SHP1[0] = SHP1i
SHP1Act[0] = SHP1Acti
SOCS3[0] = SOCS3i
SOCS3RNA[0] = SOCS3RNAi
SOCS3nRNA1[0] = SOCS3nRNA1i
SOCS3nRNA2[0] = SOCS3nRNA2i
SOCS3nRNA3[0] = SOCS3nRNA3i
SOCS3nRNA4[0] = SOCS3nRNA4i
SOCS3nRNA5[0] = SOCS3nRNA5i
STAT5[0] = STAT5i
npSTAT5[0] = npSTAT5i
p12EpoRpJAK2[0] = p12EpoRpJAK2i
p1EpoRpJAK2[0] = p1EpoRpJAK2i
p2EpoRpJAK2[0] = p2EpoRpJAK2i
pSTAT5[0] = pSTAT5i
end initial conditions

begin initial values
CISi = 0.0
CISRNAi = 0.0
CISnRNA1i = 0.0
CISnRNA2i = 0.0
CISnRNA3i = 0.0
CISnRNA4i = 0.0
CISnRNA5i = 0.0
Epoi = 1.24997e-07
EpoRJAK2i = 3.97622
EpoRJAK2CISi = 0.0
EpoRpJAK2i = 0.0
SHP1i = 26.7251
SHP1Acti = 0.0
SOCS3i = 0.0
SOCS3RNAi = 0.0
SOCS3nRNA1i = 0.0
SOCS3nRNA2i = 0.0
SOCS3nRNA3i = 0.0
SOCS3nRNA4i = 0.0
SOCS3nRNA5i = 0.0
STAT5i = 79.7535
npSTAT5i = 0.0
p12EpoRpJAK2i = 0.0
p1EpoRpJAK2i = 0.0
p2EpoRpJAK2i = 0.0
pSTAT5i = 0.0
end initial values

begin assignment rules
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