(* Generated by JWS Online *) (* This is an experimental feature of JWS Online. Please report any mistakes.*) (* Note that the following notable SBML entities or features are not supported in notebook outputyet: *) (* Events *) (* Constraints *) (* Units and UnitDefinitions *) (* AlgebraicRules *) (* conversionFactors *) variables = { IFNR[t], IFNRJ[t], IFNRJ2[t], IFNRJ2\[LetterSpace]star[t], IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2[t], IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1[t], IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t], IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]STAT1c[t], IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1[t], IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t], IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c[t], IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[t], JAK[t], PPN[t], PPX[t], R[t], RJ[t], SHP2[t], SOCS1[t], STAT1c[t], STAT1c\[LetterSpace]STAT1c\[LetterSpace]star[t], STAT1c\[LetterSpace]star[t], STAT1c\[LetterSpace]star\[LetterSpace]PPX[t], STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[t], STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star\[LetterSpace]PPX[t], STAT1n[t], STAT1n\[LetterSpace]STAT1n\[LetterSpace]star[t], STAT1n\[LetterSpace]star[t], STAT1n\[LetterSpace]star\[LetterSpace]PPN[t], STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star[t], STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star\[LetterSpace]PPN[t], mRNAc[t], mRNAn[t] }; initialValues = { IFNR[0] == 0.0, IFNRJ[0] == 0.0, IFNRJ2[0] == 0.0, IFNRJ2\[LetterSpace]star[0] == 0.0, IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2[0] == 0.0, IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1[0] == 0.0, IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1\[LetterSpace]STAT1c[0] == 0.0, IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]STAT1c[0] == 0.0, IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1[0] == 0.0, IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1\[LetterSpace]STAT1c[0] == 0.0, IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c[0] == 0.0, IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[0] == 0.0, JAK[0] == 10.0, PPN[0] == 60.0, PPX[0] == 50.0, R[0] == 10.0, RJ[0] == 0.0, SHP2[0] == 100.0, SOCS1[0] == 0.0, STAT1c[0] == 1000.0, STAT1c\[LetterSpace]STAT1c\[LetterSpace]star[0] == 0.0, STAT1c\[LetterSpace]star[0] == 0.0, STAT1c\[LetterSpace]star\[LetterSpace]PPX[0] == 0.0, STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[0] == 0.0, STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star\[LetterSpace]PPX[0] == 0.0, STAT1n[0] == 0.0, STAT1n\[LetterSpace]STAT1n\[LetterSpace]star[0] == 0.0, STAT1n\[LetterSpace]star[0] == 0.0, STAT1n\[LetterSpace]star\[LetterSpace]PPN[0] == 0.0, STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star[0] == 0.0, STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star\[LetterSpace]PPN[0] == 0.0, mRNAc[0] == 0.0, mRNAn[0] == 0.0 }; rates = { v1, v10, v11, v12, v13, v14, v15, v16, v17, v18, v19, v2, v20, v21, v22, v23, v25, v26, v27, v28, v29, v3, v30, v31, v32, v33, v34, v35, v36, v37, v38, v39, v4, v40, v41, v42, v43, v44, v45, v46, v5, v6, v7, v8, v9 }; rateEquations = { v1 -> cytoplasm*(v1\[LetterSpace]kf*JAK[t]*R[t] - v1\[LetterSpace]kb*RJ[t]), v10 -> cytoplasm*v10\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2[t], v11 -> cytoplasm*(v11\[LetterSpace]kf*PPX[t]*STAT1c\[LetterSpace]star[t] - v11\[LetterSpace]kb*STAT1c\[LetterSpace]star\[LetterSpace]PPX[t]), v12 -> cytoplasm*v12\[LetterSpace]kf*STAT1c\[LetterSpace]star\[LetterSpace]PPX[t], v13 -> cytoplasm*(v13\[LetterSpace]kf*PPX[t]*STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[t] - v13\[LetterSpace]kb*STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star\[LetterSpace]PPX[t]), v14 -> cytoplasm*v14\[LetterSpace]kf*STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star\[LetterSpace]PPX[t], v15 -> cytoplasm*(v15\[LetterSpace]kf*STAT1c[t]*STAT1c\[LetterSpace]star[t] - v15\[LetterSpace]kb*STAT1c\[LetterSpace]STAT1c\[LetterSpace]star[t]), v16 -> cytoplasm*v16\[LetterSpace]kf*STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[t], v17 -> nucleus*(v17\[LetterSpace]kf*STAT1n\[LetterSpace]star[t]^2 - v17\[LetterSpace]kb*STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star[t]), v18 -> nucleus*(v18\[LetterSpace]kf*PPN[t]*STAT1n\[LetterSpace]star[t] - v18\[LetterSpace]kb*STAT1n\[LetterSpace]star\[LetterSpace]PPN[t]), v19 -> nucleus*v19\[LetterSpace]kf*STAT1n\[LetterSpace]star\[LetterSpace]PPN[t], v2 -> cytoplasm*(-(v2\[LetterSpace]kb*IFNRJ[t]) + IFN*v2\[LetterSpace]kf*RJ[t]), v20 -> nucleus*(v20\[LetterSpace]kf*PPN[t]*STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star[t] - v20\[LetterSpace]kb*STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star\[LetterSpace]PPN[t]), v21 -> nucleus*v21\[LetterSpace]kf*STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star\[LetterSpace]PPN[t], v22 -> nucleus*(v22\[LetterSpace]kf*STAT1n[t]*STAT1n\[LetterSpace]star[t] - v22\[LetterSpace]kb*STAT1n\[LetterSpace]STAT1n\[LetterSpace]star[t]), v23 -> nucleus*v23\[LetterSpace]kf*STAT1n[t], v25 -> nucleus*v25\[LetterSpace]kf*mRNAn[t], v26 -> cytoplasm*v26\[LetterSpace]kf*mRNAc[t], v27 -> cytoplasm*v27\[LetterSpace]kf*mRNAc[t], v28 -> cytoplasm*v28\[LetterSpace]kf*SOCS1[t], v29 -> cytoplasm*(-(v29\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1[t]) + v29\[LetterSpace]kf*IFNRJ2\[LetterSpace]star[t]*SOCS1[t]), v3 -> cytoplasm*(v3\[LetterSpace]kf*IFNRJ[t]^2 - v3\[LetterSpace]kb*IFNRJ2[t]), v30 -> cytoplasm*(-(v30\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t]) + v30\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1[t]*STAT1c[t]), v31 -> cytoplasm*(-(v31\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t]) + v31\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t]*SHP2[t]), v32 -> nucleus*v32\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t], v33 -> cytoplasm*v33\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t], v34 -> cytoplasm*(-(v34\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1[t]) + v34\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1[t]*SHP2[t]), v35 -> cytoplasm*(-(v35\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t]) + v35\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1[t]*STAT1c[t]), v36 -> cytoplasm*(-(v36\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]STAT1c[t]) + v36\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c[t]*SHP2[t]), v37 -> cytoplasm*v37\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]STAT1c[t], v38 -> cytoplasm*v38\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t], v39 -> cytoplasm*v39\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1[t], v4 -> cytoplasm*v4\[LetterSpace]kf*IFNRJ2[t], v40 -> cytoplasm*v40\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1[t], v41 -> cytoplasm*v41\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1[t], v42 -> cytoplasm*(-(v42\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t]) + v42\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c[t]*SOCS1[t]), v43 -> cytoplasm*(-(v43\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1[t]) + v43\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2[t]*SOCS1[t]), v44 -> cytoplasm*(-(v44\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1\[LetterSpace]STAT1c[t]) + v44\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]STAT1c[t]*SOCS1[t]), v45 -> cytoplasm*(-(v45\[LetterSpace]kb*IFNR[t]) + IFN*v45\[LetterSpace]kf*R[t]), v46 -> cytoplasm*(-(v46\[LetterSpace]kb*IFNRJ[t]) + v46\[LetterSpace]kf*IFNR[t]*JAK[t]), v5 -> cytoplasm*(-(v5\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c[t]) + v5\[LetterSpace]kf*IFNRJ2\[LetterSpace]star[t]*STAT1c[t]), v6 -> cytoplasm*v6\[LetterSpace]kf*IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c[t], v7 -> cytoplasm*(-(v7\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[t]) + v7\[LetterSpace]kf*IFNRJ2\[LetterSpace]star[t]*STAT1c\[LetterSpace]star[t]), v8 -> cytoplasm*(v8\[LetterSpace]kf*STAT1c\[LetterSpace]star[t]^2 - v8\[LetterSpace]kb*STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[t]), v9 -> cytoplasm*(-(v9\[LetterSpace]kb*IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2[t]) + v9\[LetterSpace]kf*IFNRJ2\[LetterSpace]star[t]*SHP2[t]) }; parameters = { IFN -> 10.0, v1\[LetterSpace]kf -> 0.1, v1\[LetterSpace]kb -> 0.05, v2\[LetterSpace]kf -> 0.02, v2\[LetterSpace]kb -> 0.02, v3\[LetterSpace]kf -> 0.04, v3\[LetterSpace]kb -> 0.2, v4\[LetterSpace]kf -> 0.005, v5\[LetterSpace]kf -> 0.008, v5\[LetterSpace]kb -> 0.8, v6\[LetterSpace]kf -> 0.4, v7\[LetterSpace]kf -> 0.005, v7\[LetterSpace]kb -> 0.5, v8\[LetterSpace]kf -> 0.02, v8\[LetterSpace]kb -> 0.1, v9\[LetterSpace]kf -> 0.001, v9\[LetterSpace]kb -> 0.2, v10\[LetterSpace]kf -> 0.003, v11\[LetterSpace]kf -> 0.001, v11\[LetterSpace]kb -> 0.2, v12\[LetterSpace]kf -> 0.003, v13\[LetterSpace]kf -> 0.001, v13\[LetterSpace]kb -> 0.2, v14\[LetterSpace]kf -> 0.003, v15\[LetterSpace]kf -> 2*^-07, v15\[LetterSpace]kb -> 0.2, v16\[LetterSpace]kf -> 0.005, v17\[LetterSpace]kf -> 0.02, v17\[LetterSpace]kb -> 0.1, v18\[LetterSpace]kf -> 0.001, v18\[LetterSpace]kb -> 0.2, v19\[LetterSpace]kf -> 0.005, v20\[LetterSpace]kf -> 0.001, v20\[LetterSpace]kb -> 0.2, v21\[LetterSpace]kf -> 0.005, v22\[LetterSpace]kf -> 2*^-07, v22\[LetterSpace]kb -> 0.2, v23\[LetterSpace]kf -> 0.05, v25\[LetterSpace]kf -> 0.001, v26\[LetterSpace]kf -> 0.01, v27\[LetterSpace]kf -> 0.0005, v28\[LetterSpace]kf -> 0.0005, v29\[LetterSpace]kf -> 0.02, v29\[LetterSpace]kb -> 0.1, v30\[LetterSpace]kf -> 0.008, v30\[LetterSpace]kb -> 0.8, v31\[LetterSpace]kf -> 0.001, v31\[LetterSpace]kb -> 0.2, v32\[LetterSpace]kf -> 0.003, v33\[LetterSpace]kf -> 0.0005, v34\[LetterSpace]kf -> 0.001, v34\[LetterSpace]kb -> 0.2, v35\[LetterSpace]kf -> 0.008, v35\[LetterSpace]kb -> 0.8, v36\[LetterSpace]kf -> 0.001, v36\[LetterSpace]kb -> 0.2, v37\[LetterSpace]kf -> 0.003, v38\[LetterSpace]kf -> 0.0005, v39\[LetterSpace]kf -> 0.0005, v40\[LetterSpace]kf -> 0.003, v41\[LetterSpace]kf -> 0.0005, v42\[LetterSpace]kf -> 0.02, v42\[LetterSpace]kb -> 0.1, v43\[LetterSpace]kf -> 0.02, v43\[LetterSpace]kb -> 0.1, v44\[LetterSpace]kf -> 0.02, v44\[LetterSpace]kb -> 0.1, v45\[LetterSpace]kf -> 0.02, v45\[LetterSpace]kb -> 0.02, v46\[LetterSpace]kf -> 0.1, v46\[LetterSpace]kb -> 0.05, cytoplasm -> 1.0, nucleus -> 1.0 }; assignments = { }; events = { }; speciesAnnotations = { IFN[t]->"http://identifiers.org/uniprot/P01572", IFNRJ[t]->"http://identifiers.org/uniprot/P52332", IFNRJ[t]->"http://identifiers.org/uniprot/P15261", IFNRJ[t]->"http://identifiers.org/uniprot/P01572", JAK[t]->"http://identifiers.org/uniprot/P52332", JAK[t]->"http://identifiers.org/interpro/IPR009127", R[t]->"http://identifiers.org/uniprot/P15261", R[t]->"http://identifiers.org/interpro/IPR008355", SHP2[t]->"http://identifiers.org/uniprot/P35235", SOCS1[t]->"http://identifiers.org/uniprot/O35716", SOCS1[t]->"http://identifiers.org/interpro/IPR001496", STAT1c[t]->"http://identifiers.org/uniprot/P42225", STAT1c[t]->"http://identifiers.org/interpro/IPR001217", STAT1c\[LetterSpace]STAT1c\[LetterSpace]star[t]->"http://identifiers.org/uniprot/P42225", STAT1c\[LetterSpace]star[t]->"http://identifiers.org/uniprot/P42225", STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star[t]->"http://identifiers.org/uniprot/P42225", STAT1n[t]->"http://identifiers.org/uniprot/P42225", STAT1n\[LetterSpace]STAT1n\[LetterSpace]star[t]->"http://identifiers.org/uniprot/P42225", STAT1n\[LetterSpace]star[t]->"http://identifiers.org/uniprot/P42225", STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star[t]->"http://identifiers.org/uniprot/P42225", mRNAc[t]->"http://identifiers.org/chebi/CHEBI:33699", mRNAc[t]->"http://identifiers.org/kegg.compound/C00046", mRNAn[t]->"http://identifiers.org/chebi/CHEBI:33699", mRNAn[t]->"http://identifiers.org/kegg.compound/C00046" }; reactionAnnotations = { v1->"http://identifiers.org/go/GO:0005515", v1->"http://identifiers.org/go/GO:0005133", v10->"http://identifiers.org/ec-code/3.1.3.48", v10->"http://identifiers.org/kegg.reaction/R02585", v10->"http://identifiers.org/go/GO:0005001", v11->"http://identifiers.org/go/GO:0019902", v13->"http://identifiers.org/go/GO:0019902", v14->"http://identifiers.org/go/GO:0042512", v15->"http://identifiers.org/go/GO:0005515", v16->"http://identifiers.org/go/GO:0007262", v17->"http://identifiers.org/go/GO:0007261", v18->"http://identifiers.org/go/GO:0005515", v18->"http://identifiers.org/go/GO:0019902", v19->"http://identifiers.org/go/GO:0042512", v2->"http://identifiers.org/go/GO:0004906", v2->"http://identifiers.org/go/GO:0019964", v20->"http://identifiers.org/go/GO:0019902", v21->"http://identifiers.org/go/GO:0042512", v22->"http://identifiers.org/go/GO:0007261", v23->"http://identifiers.org/go/GO:0065002", v25->"http://identifiers.org/go/GO:0006406", v26->"http://identifiers.org/go/GO:0006412", v27->"http://identifiers.org/go/GO:0006402", v28->"http://identifiers.org/go/GO:0044257", v29->"http://identifiers.org/go/GO:0005515", v3->"http://identifiers.org/go/GO:0007171", v30->"http://identifiers.org/go/GO:0005515", v31->"http://identifiers.org/go/GO:0005515", v31->"http://identifiers.org/go/GO:0019902", v32->"http://identifiers.org/go/GO:0005001", v34->"http://identifiers.org/go/GO:0005515", v34->"http://identifiers.org/go/GO:0019902", v35->"http://identifiers.org/go/GO:0005515", v36->"http://identifiers.org/go/GO:0005515", v36->"http://identifiers.org/go/GO:0019902", v37->"http://identifiers.org/go/GO:0043624", v37->"http://identifiers.org/go/GO:0005001", v4->"http://identifiers.org/ec-code/2.7.10.2", v4->"http://identifiers.org/kegg.reaction/R02584", v4->"http://identifiers.org/go/GO:0004713", v40->"http://identifiers.org/go/GO:0005001", v42->"http://identifiers.org/go/GO:0005515", v43->"http://identifiers.org/go/GO:0005515", v44->"http://identifiers.org/go/GO:0005515", v45->"http://identifiers.org/go/GO:0005133", v45->"http://identifiers.org/go/GO:0019964", v46->"http://identifiers.org/go/GO:0019964", v5->"http://identifiers.org/go/GO:0005515", v6->"http://identifiers.org/go/GO:0007260", v7->"http://identifiers.org/go/GO:0005515", v8->"http://identifiers.org/go/GO:0007261", v9->"http://identifiers.org/go/GO:0019902" }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { IFNR'[t] == 1.0*v45 -1.0*v46, IFNRJ'[t] == 1.0*v2 +1.0*v46 -2.0*v3, IFNRJ2'[t] == 1.0*v3 +1.0*v10 +1.0*v32 +1.0*v37 +1.0*v40 -1.0*v4, IFNRJ2\[LetterSpace]star'[t] == 1.0*v4 +1.0*v6 +1.0*v41 -1.0*v5 -1.0*v7 -1.0*v9 -1.0*v29, IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2'[t] == 1.0*v9 +1.0*v39 -1.0*v10 -1.0*v43, IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1'[t] == 1.0*v34 +1.0*v43 -1.0*v35 -1.0*v39 -1.0*v40, IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]SOCS1\[LetterSpace]STAT1c'[t] == 1.0*v31 +1.0*v35 +1.0*v44 -1.0*v32 -1.0*v33, IFNRJ2\[LetterSpace]star\[LetterSpace]SHP2\[LetterSpace]STAT1c'[t] == 1.0*v33 +1.0*v36 -1.0*v37 -1.0*v44, IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1'[t] == 1.0*v29 -1.0*v30 -1.0*v34 -1.0*v41, IFNRJ2\[LetterSpace]star\[LetterSpace]SOCS1\[LetterSpace]STAT1c'[t] == 1.0*v30 +1.0*v42 -1.0*v31 -1.0*v38, IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c'[t] == 1.0*v5 +1.0*v38 -1.0*v6 -1.0*v36 -1.0*v42, IFNRJ2\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star'[t] == 1.0*v7 , JAK'[t] == -1.0*v1 -1.0*v46, PPN'[t] == 1.0*v19 +1.0*v21 -1.0*v18 -1.0*v20, PPX'[t] == 1.0*v12 +1.0*v14 -1.0*v11 -1.0*v13, R'[t] == -1.0*v1 -1.0*v45, RJ'[t] == 1.0*v1 -1.0*v2, SHP2'[t] == 1.0*v10 +1.0*v32 +1.0*v37 +1.0*v40 -1.0*v9 -1.0*v31 -1.0*v34 -1.0*v36, SOCS1'[t] == 1.0*v26 +1.0*v32 +1.0*v40 -1.0*v28 -1.0*v29 -1.0*v42 -1.0*v43 -1.0*v44, STAT1c'[t] == 1.0*v12 +1.0*v23 +1.0*v32 +1.0*v37 -1.0*v5 -1.0*v15 -1.0*v30 -1.0*v35, STAT1c\[LetterSpace]STAT1c\[LetterSpace]star'[t] == 1.0*v14 +1.0*v15 , STAT1c\[LetterSpace]star'[t] == 1.0*v6 -1.0*v7 -2.0*v8 -1.0*v11 -1.0*v15, STAT1c\[LetterSpace]star\[LetterSpace]PPX'[t] == 1.0*v11 -1.0*v12, STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star'[t] == 1.0*v8 -1.0*v13 -1.0*v16, STAT1c\[LetterSpace]star\[LetterSpace]STAT1c\[LetterSpace]star\[LetterSpace]PPX'[t] == 1.0*v13 -1.0*v14, STAT1n'[t] == 1.0*v19 -1.0*v22 -1.0*v23, STAT1n\[LetterSpace]STAT1n\[LetterSpace]star'[t] == 1.0*v21 +1.0*v22 , STAT1n\[LetterSpace]star'[t] == -2.0*v17 -1.0*v18 -1.0*v22, STAT1n\[LetterSpace]star\[LetterSpace]PPN'[t] == 1.0*v18 -1.0*v19, STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star'[t] == 1.0*v16 +1.0*v17 -1.0*v20, STAT1n\[LetterSpace]star\[LetterSpace]STAT1n\[LetterSpace]star\[LetterSpace]PPN'[t] == 1.0*v20 -1.0*v21, mRNAc'[t] == 1.0*v25 -1.0*v27, mRNAn'[t] == -1.0*v25 }; timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}]; (* Steady-state solution initialized with result of time evolution *) findRootEquations = odes /.D[_[t],t]->0; findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2]; steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100] fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations (* Plot the time evolution of the variables *) plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]