(* 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 = {
Active\[LetterSpace]R2[t], Active\[LetterSpace]Rec[t], Active\[LetterSpace]Rec\[LetterSpace]endo[t], Inactive\[LetterSpace]Rec[t], S2\[LetterSpace]c[t], S2\[LetterSpace]n[t], S4\[LetterSpace]c[t], S4\[LetterSpace]n[t], S7[t], Smad7mRNA[t], Smad7mRNA1[t], TGFb[t], TGFb\[LetterSpace]In[t], TGFb\[LetterSpace]R1\[LetterSpace]endo[t], TGFb\[LetterSpace]R1\[LetterSpace]surface[t], TGFb\[LetterSpace]R2\[LetterSpace]endo[t], TGFb\[LetterSpace]R2\[LetterSpace]surface[t], pS2\[LetterSpace]c[t], pS2\[LetterSpace]n[t], pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c[t], pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t], pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c[t], pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n[t]
};

initialValues = {
Active\[LetterSpace]R2[0] == 0.0, Active\[LetterSpace]Rec[0] == 0.0, Active\[LetterSpace]Rec\[LetterSpace]endo[0] == 0.0, Inactive\[LetterSpace]Rec[0] == 0.0, S2\[LetterSpace]c[0] == 983.8767796832187, S2\[LetterSpace]n[0] == 331.2552270775211, S4\[LetterSpace]c[0] == 321.7549460489289, S4\[LetterSpace]n[0] == 67.78337093619712, S7[0] == 5.602470383338093, Smad7mRNA[0] == 0.029984454646879696, Smad7mRNA1[0] == 0.1950853287373667, TGFb[0] == 100.0, TGFb\[LetterSpace]In[0] == 0.0, TGFb\[LetterSpace]R1\[LetterSpace]endo[0] == 17.131731480965776, TGFb\[LetterSpace]R1\[LetterSpace]surface[0] == 7.374847209795631, TGFb\[LetterSpace]R2\[LetterSpace]endo[0] == 21.632890273566286, TGFb\[LetterSpace]R2\[LetterSpace]surface[0] == 2.8485790349258178, pS2\[LetterSpace]c[0] == 0.0, pS2\[LetterSpace]n[0] == 0.0, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c[0] == 0.0, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[0] == 0.0, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c[0] == 0.0, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n[0] == 0.0
};

rates = {
reaction1, reaction10, reaction11, reaction12, reaction13, reaction14, reaction15, reaction16, reaction17, reaction18, reaction19, reaction2, reaction20, reaction21, reaction22, reaction23, reaction24, reaction25, reaction26, reaction27, reaction28, reaction29, reaction3, reaction30, reaction31, reaction32, reaction33, reaction34, reaction35, reaction36, reaction37, reaction38, reaction39, reaction4, reaction40, reaction41, reaction42, reaction43, reaction44, reaction45, reaction46, reaction47, reaction48, reaction5, reaction51, reaction52, reaction53, reaction54, reaction55, reaction56, reaction57, reaction58, reaction59, reaction6, reaction60, reaction61, reaction62, reaction63, reaction64, reaction65, reaction66, reaction67, reaction68, reaction7, reaction8, reaction9
};

rateEquations = {
reaction1 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction1\[LetterSpace]0[TGFb[t], TGFb\[LetterSpace]R2\[LetterSpace]surface[t], k\[LetterSpace]medium, kf\[LetterSpace]R2\[LetterSpace]activation], reaction10 -> cyt*kdeg\[LetterSpace]R1*Active\[LetterSpace]Rec\[LetterSpace]endo[t], reaction11 -> cyt*kin\[LetterSpace]deg\[LetterSpace]Ligand*Active\[LetterSpace]Rec\[LetterSpace]endo[t], reaction12 -> cyt*kdeg\[LetterSpace]R1*TGFb\[LetterSpace]R1\[LetterSpace]endo[t], reaction13 -> cyt*kdeg\[LetterSpace]R2*TGFb\[LetterSpace]R2\[LetterSpace]endo[t], reaction14 -> cyt*kf\[LetterSpace]Seq\[LetterSpace]S7\[LetterSpace]Rec*Active\[LetterSpace]Rec[t]*S7[t], reaction15 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction15\[LetterSpace]0[Inactive\[LetterSpace]Rec[t], index\[LetterSpace]seq\[LetterSpace]kb, kf\[LetterSpace]Seq\[LetterSpace]S7\[LetterSpace]Rec], reaction16 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction16\[LetterSpace]0[Inactive\[LetterSpace]Rec[t], index\[LetterSpace]induced\[LetterSpace]ligand\[LetterSpace]deg, kin\[LetterSpace]deg\[LetterSpace]Ligand], reaction17 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction17\[LetterSpace]0[Inactive\[LetterSpace]Rec[t], index\[LetterSpace]induced\[LetterSpace]R2\[LetterSpace]deg, kdeg\[LetterSpace]R2], reaction18 -> cyt*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4*TGFb\[LetterSpace]R1\[LetterSpace]surface[t], reaction19 -> cyt*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4*TGFb\[LetterSpace]R2\[LetterSpace]surface[t], reaction2 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction2\[LetterSpace]0[TGFb[t], TGFb\[LetterSpace]R2\[LetterSpace]surface[t], kf\[LetterSpace]R2\[LetterSpace]activation], reaction20 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction20\[LetterSpace]1[TGFb\[LetterSpace]R1\[LetterSpace]endo[t], index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4, k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4], reaction21 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction21\[LetterSpace]1[TGFb\[LetterSpace]R2\[LetterSpace]endo[t], index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4, k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4], reaction22 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction22\[LetterSpace]1[Active\[LetterSpace]Rec[t], index\[LetterSpace]active\[LetterSpace]Rec\[LetterSpace]internalize\[LetterSpace]C4, k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4], reaction23 -> cyt*k\[LetterSpace]disso\[LetterSpace]Active\[LetterSpace]Rec*Active\[LetterSpace]Rec\[LetterSpace]endo[t], reaction24 -> cyt*kin\[LetterSpace]deg\[LetterSpace]Ligand*TGFb\[LetterSpace]In[t], reaction25 -> cyt*kout\[LetterSpace]deg\[LetterSpace]Ligand\[LetterSpace]100pM*TGFb[t], reaction26 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction26\[LetterSpace]0[Active\[LetterSpace]Rec\[LetterSpace]endo[t], S2\[LetterSpace]c[t], k\[LetterSpace]phosphorylation], reaction27 -> cyt*kf\[LetterSpace]trimmer*pS2\[LetterSpace]c[t]^2*S4\[LetterSpace]c[t], reaction28 -> cyt*kb\[LetterSpace]trimmer*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c[t], reaction29 -> kf\[LetterSpace]trimmer*nuc*pS2\[LetterSpace]n[t]^2*S4\[LetterSpace]n[t], reaction3 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction3\[LetterSpace]0[TGFb[t], TGFb\[LetterSpace]R2\[LetterSpace]surface[t], kf\[LetterSpace]R2\[LetterSpace]activation], reaction30 -> kb\[LetterSpace]trimmer*nuc*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t], reaction31 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction31\[LetterSpace]2[index\[LetterSpace]kf\[LetterSpace]homotrimer, kf\[LetterSpace]trimmer, pS2\[LetterSpace]c[t]], reaction32 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction32\[LetterSpace]3[index\[LetterSpace]kb\[LetterSpace]homotrimer, kb\[LetterSpace]trimmer, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c[t]], reaction33 -> nuc*Function\[LetterSpace]for\[LetterSpace]reaction33\[LetterSpace]2[index\[LetterSpace]kf\[LetterSpace]homotrimer, kf\[LetterSpace]trimmer, pS2\[LetterSpace]n[t]], reaction34 -> nuc*Function\[LetterSpace]for\[LetterSpace]reaction34\[LetterSpace]2[index\[LetterSpace]kb\[LetterSpace]homotrimer, kb\[LetterSpace]trimmer, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n[t]], reaction35 -> Function\[LetterSpace]for\[LetterSpace]reaction35\[LetterSpace]2[S2\[LetterSpace]c[t], S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4, cyt], reaction36 -> Function\[LetterSpace]for\[LetterSpace]reaction36\[LetterSpace]2[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4, S2\[LetterSpace]n[t], nuc], reaction37 -> Function\[LetterSpace]for\[LetterSpace]reaction37\[LetterSpace]2[S4\[LetterSpace]c[t], S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4, cyt], reaction38 -> Function\[LetterSpace]for\[LetterSpace]reaction38\[LetterSpace]2[S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4, S4\[LetterSpace]n[t], nuc], reaction39 -> Function\[LetterSpace]for\[LetterSpace]reaction39\[LetterSpace]2[S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4, cyt, pS2\[LetterSpace]c[t]], reaction4 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction4\[LetterSpace]0[Active\[LetterSpace]R2[t], index\[LetterSpace]kb\[LetterSpace]R2, kf\[LetterSpace]R2\[LetterSpace]activation], reaction40 -> Function\[LetterSpace]for\[LetterSpace]reaction40\[LetterSpace]2[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4, nuc, pS2\[LetterSpace]n[t]], reaction41 -> Function\[LetterSpace]for\[LetterSpace]reaction41\[LetterSpace]2[Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc, cyt, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c[t]], reaction42 -> Function\[LetterSpace]for\[LetterSpace]reaction42\[LetterSpace]0[Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc, cyt, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c[t]], reaction43 -> k\[LetterSpace]Dephos\[LetterSpace]C4*nuc*pS2\[LetterSpace]n[t], reaction44 -> nuc*Function\[LetterSpace]for\[LetterSpace]reaction44\[LetterSpace]2[index\[LetterSpace]trimer\[LetterSpace]dephos, k\[LetterSpace]Dephos\[LetterSpace]C4, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t]], reaction45 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction45\[LetterSpace]2[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4, S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4, S2\[LetterSpace]total\[LetterSpace]C4, kdeg\[LetterSpace]S2], reaction46 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction46\[LetterSpace]2[S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4, S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4, S4\[LetterSpace]total\[LetterSpace]C4, kdeg\[LetterSpace]S4], reaction47 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction47\[LetterSpace]2[R1\[LetterSpace]total\[LetterSpace]C4, index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4, k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4, kdeg\[LetterSpace]R1], reaction48 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction48\[LetterSpace]0[R2\[LetterSpace]total\[LetterSpace]C4, index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4, k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4, kdeg\[LetterSpace]R2], reaction5 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction5\[LetterSpace]0[Active\[LetterSpace]R2[t], index\[LetterSpace]kb\[LetterSpace]R2, kf\[LetterSpace]R2\[LetterSpace]activation], reaction51 -> nuc*Function\[LetterSpace]for\[LetterSpace]reaction51\[LetterSpace]2[DRB, K\[LetterSpace]mran, S7KD, hill\[LetterSpace]fact1, k\[LetterSpace]induced\[LetterSpace]S7\[LetterSpace]production\[LetterSpace]C4, mRNA\[LetterSpace]prod\[LetterSpace]C4, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t]], reaction52 -> Function\[LetterSpace]for\[LetterSpace]reaction52\[LetterSpace]3[Smad7mRNA1[t], export\[LetterSpace]cytoplasm, nuc], reaction53 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction53\[LetterSpace]0[Smad7mRNA[t], k\[LetterSpace]S7\[LetterSpace]protein\[LetterSpace]C4], reaction54 -> cyt*kdeg\[LetterSpace]S2*S2\[LetterSpace]c[t], reaction55 -> kdeg\[LetterSpace]S2*nuc*S2\[LetterSpace]n[t], reaction56 -> kdeg\[LetterSpace]S2*nuc*pS2\[LetterSpace]n[t], reaction57 -> cyt*kdeg\[LetterSpace]S2*pS2\[LetterSpace]c[t], reaction58 -> cyt*kdeg\[LetterSpace]S4*S4\[LetterSpace]c[t], reaction59 -> kdeg\[LetterSpace]S4*nuc*S4\[LetterSpace]n[t], reaction6 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction6\[LetterSpace]0[Active\[LetterSpace]R2[t], index\[LetterSpace]kb\[LetterSpace]R2, k\[LetterSpace]medium, kf\[LetterSpace]R2\[LetterSpace]activation], reaction60 -> cyt*kdeg\[LetterSpace]S2*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c[t], reaction61 -> cyt*kdeg\[LetterSpace]S4*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c[t], reaction62 -> cyt*kdeg\[LetterSpace]S2*pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c[t], reaction63 -> kdeg\[LetterSpace]S2*nuc*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t], reaction64 -> kdeg\[LetterSpace]S4*nuc*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t], reaction65 -> kdeg\[LetterSpace]S2*nuc*pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n[t], reaction66 -> kmRNA1deg\[LetterSpace]S7*nuc*Smad7mRNA1[t], reaction67 -> cyt*kmRNAdeg\[LetterSpace]S7*Smad7mRNA[t], reaction68 -> cyt*kdeg\[LetterSpace]S7*S7[t], reaction7 -> cyt*kf\[LetterSpace]R1\[LetterSpace]activation*Active\[LetterSpace]R2[t]*TGFb\[LetterSpace]R1\[LetterSpace]surface[t], reaction8 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction8\[LetterSpace]0[Active\[LetterSpace]Rec[t], index\[LetterSpace]kb\[LetterSpace]R1, kf\[LetterSpace]R1\[LetterSpace]activation], reaction9 -> cyt*kdeg\[LetterSpace]R2*Active\[LetterSpace]Rec\[LetterSpace]endo[t]
};

parameters = {
DRB -> 1.0, K\[LetterSpace]mran -> 223.586, R1\[LetterSpace]total\[LetterSpace]C4 -> 24.5065786907614, R2\[LetterSpace]total\[LetterSpace]C4 -> 24.4814693084921, S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4 -> 3.1018054766891, S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4 -> 0.522250139330463, S2\[LetterSpace]total\[LetterSpace]C4 -> 1315.13200676074, S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4 -> 0.23373067315266, S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4 -> 0.0251255663556574, S4\[LetterSpace]total\[LetterSpace]C4 -> 389.538316985126, S7KD -> 1.0, Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc -> 0.144035002432848, export\[LetterSpace]cytoplasm -> 0.0304372961307751, hill\[LetterSpace]fact1 -> 3.98630224594983, index\[LetterSpace]active\[LetterSpace]Rec\[LetterSpace]internalize\[LetterSpace]C4 -> 0.922503805466932, index\[LetterSpace]induced\[LetterSpace]R2\[LetterSpace]deg -> 1.00037776866312, index\[LetterSpace]induced\[LetterSpace]ligand\[LetterSpace]deg -> 2.72106513634098, index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4 -> 0.28410342378491, index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4 -> 0.13081029215514, index\[LetterSpace]kb\[LetterSpace]R1 -> 1.61290474569122, index\[LetterSpace]kb\[LetterSpace]R2 -> 8.70498854136694, index\[LetterSpace]kb\[LetterSpace]homotrimer -> 0.0163921522930962, index\[LetterSpace]kf\[LetterSpace]homotrimer -> 4.4857919268227, index\[LetterSpace]seq\[LetterSpace]kb -> 0.000865131944828012, index\[LetterSpace]trimer\[LetterSpace]dephos -> 2.61013820571304, k\[LetterSpace]Dephos\[LetterSpace]C4 -> 0.0704492217625347, k\[LetterSpace]S7\[LetterSpace]protein\[LetterSpace]C4 -> 0.191333551016171, k\[LetterSpace]disso\[LetterSpace]Active\[LetterSpace]Rec -> 0.00103940411358367, k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4 -> 0.40073987230226, k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4 -> 1.1527121689746, k\[LetterSpace]induced\[LetterSpace]S7\[LetterSpace]production\[LetterSpace]C4 -> 9.74383528364265, k\[LetterSpace]medium -> 2000000.0, k\[LetterSpace]phosphorylation -> 0.0701560191772152, kb\[LetterSpace]trimmer -> 1.6077346670684, kdeg\[LetterSpace]R1 -> 0.0586584495713968, kdeg\[LetterSpace]R2 -> 0.0010003937382472, kdeg\[LetterSpace]S2 -> 0.000514084012175461, kdeg\[LetterSpace]S4 -> 0.00480202699360057, kdeg\[LetterSpace]S7 -> 0.00102401829734485, kf\[LetterSpace]R1\[LetterSpace]activation -> 4.98548548076033, kf\[LetterSpace]R2\[LetterSpace]activation -> 4.95559496701657, kf\[LetterSpace]Seq\[LetterSpace]S7\[LetterSpace]Rec -> 0.882455038106022, kf\[LetterSpace]trimmer -> 4.81424700307387, kin\[LetterSpace]deg\[LetterSpace]Ligand -> 0.720182234635314, kmRNA1deg\[LetterSpace]S7 -> 0.00523181439900629, kmRNAdeg\[LetterSpace]S7 -> 0.0990158065483902, kout\[LetterSpace]deg\[LetterSpace]Ligand\[LetterSpace]100pM -> 0.00036189622, mRNA\[LetterSpace]prod\[LetterSpace]C4 -> 0.00695852015347187, cyt -> 2.0, nuc -> 1.0
};

assignments = {
Function\[LetterSpace]for\[LetterSpace]reaction39\[LetterSpace]2[S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4_,cyt_,pS2\[LetterSpace]c_] -> cyt*pS2\[LetterSpace]c*S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4, Function\[LetterSpace]for\[LetterSpace]reaction35\[LetterSpace]2[S2\[LetterSpace]c_,S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4_,cyt_] -> cyt*S2\[LetterSpace]c*S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4, Function\[LetterSpace]for\[LetterSpace]reaction33\[LetterSpace]2[index\[LetterSpace]kf\[LetterSpace]homotrimer_,kf\[LetterSpace]trimmer_,pS2\[LetterSpace]n_] -> index\[LetterSpace]kf\[LetterSpace]homotrimer*kf\[LetterSpace]trimmer*pS2\[LetterSpace]n^3, Function\[LetterSpace]for\[LetterSpace]reaction36\[LetterSpace]2[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4_,S2\[LetterSpace]n_,nuc_] -> nuc*S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4*S2\[LetterSpace]n, Function\[LetterSpace]for\[LetterSpace]reaction38\[LetterSpace]2[S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4_,S4\[LetterSpace]n_,nuc_] -> nuc*S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4*S4\[LetterSpace]n, Function\[LetterSpace]for\[LetterSpace]reaction1\[LetterSpace]0[TGFb_,TGFb\[LetterSpace]R2\[LetterSpace]surface_,k\[LetterSpace]medium_,kf\[LetterSpace]R2\[LetterSpace]activation_] -> (kf\[LetterSpace]R2\[LetterSpace]activation*TGFb*TGFb\[LetterSpace]R2\[LetterSpace]surface)/k\[LetterSpace]medium, Function\[LetterSpace]for\[LetterSpace]reaction2\[LetterSpace]0[TGFb_,TGFb\[LetterSpace]R2\[LetterSpace]surface_,kf\[LetterSpace]R2\[LetterSpace]activation_] -> kf\[LetterSpace]R2\[LetterSpace]activation*TGFb*TGFb\[LetterSpace]R2\[LetterSpace]surface, Function\[LetterSpace]for\[LetterSpace]reaction6\[LetterSpace]0[Active\[LetterSpace]R2_,index\[LetterSpace]kb\[LetterSpace]R2_,k\[LetterSpace]medium_,kf\[LetterSpace]R2\[LetterSpace]activation_] -> (Active\[LetterSpace]R2*index\[LetterSpace]kb\[LetterSpace]R2*kf\[LetterSpace]R2\[LetterSpace]activation)/k\[LetterSpace]medium, Function\[LetterSpace]for\[LetterSpace]reaction3\[LetterSpace]0[TGFb_,TGFb\[LetterSpace]R2\[LetterSpace]surface_,kf\[LetterSpace]R2\[LetterSpace]activation_] -> kf\[LetterSpace]R2\[LetterSpace]activation*TGFb*TGFb\[LetterSpace]R2\[LetterSpace]surface, Function\[LetterSpace]for\[LetterSpace]reaction15\[LetterSpace]0[Inactive\[LetterSpace]Rec_,index\[LetterSpace]seq\[LetterSpace]kb_,kf\[LetterSpace]Seq\[LetterSpace]S7\[LetterSpace]Rec_] -> Inactive\[LetterSpace]Rec*index\[LetterSpace]seq\[LetterSpace]kb*kf\[LetterSpace]Seq\[LetterSpace]S7\[LetterSpace]Rec, Function\[LetterSpace]for\[LetterSpace]reaction4\[LetterSpace]0[Active\[LetterSpace]R2_,index\[LetterSpace]kb\[LetterSpace]R2_,kf\[LetterSpace]R2\[LetterSpace]activation_] -> Active\[LetterSpace]R2*index\[LetterSpace]kb\[LetterSpace]R2*kf\[LetterSpace]R2\[LetterSpace]activation, Function\[LetterSpace]for\[LetterSpace]reaction5\[LetterSpace]0[Active\[LetterSpace]R2_,index\[LetterSpace]kb\[LetterSpace]R2_,kf\[LetterSpace]R2\[LetterSpace]activation_] -> Active\[LetterSpace]R2*index\[LetterSpace]kb\[LetterSpace]R2*kf\[LetterSpace]R2\[LetterSpace]activation, Function\[LetterSpace]for\[LetterSpace]reaction21\[LetterSpace]1[TGFb\[LetterSpace]R2\[LetterSpace]endo_,index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4_,k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4_] -> index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4*TGFb\[LetterSpace]R2\[LetterSpace]endo, Function\[LetterSpace]for\[LetterSpace]reaction31\[LetterSpace]2[index\[LetterSpace]kf\[LetterSpace]homotrimer_,kf\[LetterSpace]trimmer_,pS2\[LetterSpace]c_] -> index\[LetterSpace]kf\[LetterSpace]homotrimer*kf\[LetterSpace]trimmer*pS2\[LetterSpace]c^3, Function\[LetterSpace]for\[LetterSpace]reaction26\[LetterSpace]0[Active\[LetterSpace]Rec\[LetterSpace]endo_,S2\[LetterSpace]c_,k\[LetterSpace]phosphorylation_] -> Active\[LetterSpace]Rec\[LetterSpace]endo*k\[LetterSpace]phosphorylation*S2\[LetterSpace]c, Function\[LetterSpace]for\[LetterSpace]reaction52\[LetterSpace]3[Smad7mRNA1_,export\[LetterSpace]cytoplasm_,nuc_] -> export\[LetterSpace]cytoplasm*nuc*Smad7mRNA1, Function\[LetterSpace]for\[LetterSpace]reaction53\[LetterSpace]0[Smad7mRNA_,k\[LetterSpace]S7\[LetterSpace]protein\[LetterSpace]C4_] -> k\[LetterSpace]S7\[LetterSpace]protein\[LetterSpace]C4*Smad7mRNA, Function\[LetterSpace]for\[LetterSpace]reaction32\[LetterSpace]3[index\[LetterSpace]kb\[LetterSpace]homotrimer_,kb\[LetterSpace]trimmer_,pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c_] -> index\[LetterSpace]kb\[LetterSpace]homotrimer*kb\[LetterSpace]trimmer*pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c, Function\[LetterSpace]for\[LetterSpace]reaction40\[LetterSpace]2[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4_,nuc_,pS2\[LetterSpace]n_] -> nuc*pS2\[LetterSpace]n*S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4, Function\[LetterSpace]for\[LetterSpace]reaction37\[LetterSpace]2[S4\[LetterSpace]c_,S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4_,cyt_] -> cyt*S4\[LetterSpace]c*S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4, Function\[LetterSpace]for\[LetterSpace]reaction45\[LetterSpace]2[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4_,S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4_,S2\[LetterSpace]total\[LetterSpace]C4_,kdeg\[LetterSpace]S2_] -> ((kdeg\[LetterSpace]S2^2 + kdeg\[LetterSpace]S2*S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4 + kdeg\[LetterSpace]S2*S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4)*S2\[LetterSpace]total\[LetterSpace]C4)/(kdeg\[LetterSpace]S2 + S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4 + 2*S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4), Function\[LetterSpace]for\[LetterSpace]reaction41\[LetterSpace]2[Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc_,cyt_,pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c_] -> cyt*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c*Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc, Function\[LetterSpace]for\[LetterSpace]reaction44\[LetterSpace]2[index\[LetterSpace]trimer\[LetterSpace]dephos_,k\[LetterSpace]Dephos\[LetterSpace]C4_,pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n_] -> index\[LetterSpace]trimer\[LetterSpace]dephos*k\[LetterSpace]Dephos\[LetterSpace]C4*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n, Function\[LetterSpace]for\[LetterSpace]reaction48\[LetterSpace]0[R2\[LetterSpace]total\[LetterSpace]C4_,index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4_,k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4_,kdeg\[LetterSpace]R2_] -> (kdeg\[LetterSpace]R2^2*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4*R2\[LetterSpace]total\[LetterSpace]C4)/(kdeg\[LetterSpace]R2^2 + kdeg\[LetterSpace]R2*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4 + index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4*kdeg\[LetterSpace]R2*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C4), Function\[LetterSpace]for\[LetterSpace]reaction51\[LetterSpace]2[DRB_,K\[LetterSpace]mran_,S7KD_,hill\[LetterSpace]fact1_,k\[LetterSpace]induced\[LetterSpace]S7\[LetterSpace]production\[LetterSpace]C4_,mRNA\[LetterSpace]prod\[LetterSpace]C4_,pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n_] -> DRB*(mRNA\[LetterSpace]prod\[LetterSpace]C4 + (k\[LetterSpace]induced\[LetterSpace]S7\[LetterSpace]production\[LetterSpace]C4*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n^hill\[LetterSpace]fact1)/(K\[LetterSpace]mran^hill\[LetterSpace]fact1 + pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n^hill\[LetterSpace]fact1))*S7KD, Function\[LetterSpace]for\[LetterSpace]reaction47\[LetterSpace]2[R1\[LetterSpace]total\[LetterSpace]C4_,index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4_,k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4_,kdeg\[LetterSpace]R1_] -> (kdeg\[LetterSpace]R1^2*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4*R1\[LetterSpace]total\[LetterSpace]C4)/(kdeg\[LetterSpace]R1^2 + kdeg\[LetterSpace]R1*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4 + index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4*kdeg\[LetterSpace]R1*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4), Function\[LetterSpace]for\[LetterSpace]reaction42\[LetterSpace]0[Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc_,cyt_,pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c_] -> cyt*pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c*Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc, Function\[LetterSpace]for\[LetterSpace]reaction34\[LetterSpace]2[index\[LetterSpace]kb\[LetterSpace]homotrimer_,kb\[LetterSpace]trimmer_,pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n_] -> index\[LetterSpace]kb\[LetterSpace]homotrimer*kb\[LetterSpace]trimmer*pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n, Function\[LetterSpace]for\[LetterSpace]reaction46\[LetterSpace]2[S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4_,S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4_,S4\[LetterSpace]total\[LetterSpace]C4_,kdeg\[LetterSpace]S4_] -> ((kdeg\[LetterSpace]S4^2 + kdeg\[LetterSpace]S4*S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4 + kdeg\[LetterSpace]S4*S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4)*S4\[LetterSpace]total\[LetterSpace]C4)/(kdeg\[LetterSpace]S4 + S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C4 + 2*S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C4), Function\[LetterSpace]for\[LetterSpace]reaction22\[LetterSpace]1[Active\[LetterSpace]Rec_,index\[LetterSpace]active\[LetterSpace]Rec\[LetterSpace]internalize\[LetterSpace]C4_,k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4_] -> Active\[LetterSpace]Rec*index\[LetterSpace]active\[LetterSpace]Rec\[LetterSpace]internalize\[LetterSpace]C4*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4, Function\[LetterSpace]for\[LetterSpace]reaction20\[LetterSpace]1[TGFb\[LetterSpace]R1\[LetterSpace]endo_,index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4_,k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4_] -> index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C4*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C4*TGFb\[LetterSpace]R1\[LetterSpace]endo, Function\[LetterSpace]for\[LetterSpace]reaction8\[LetterSpace]0[Active\[LetterSpace]Rec_,index\[LetterSpace]kb\[LetterSpace]R1_,kf\[LetterSpace]R1\[LetterSpace]activation_] -> Active\[LetterSpace]Rec*index\[LetterSpace]kb\[LetterSpace]R1*kf\[LetterSpace]R1\[LetterSpace]activation, Function\[LetterSpace]for\[LetterSpace]reaction16\[LetterSpace]0[Inactive\[LetterSpace]Rec_,index\[LetterSpace]induced\[LetterSpace]ligand\[LetterSpace]deg_,kin\[LetterSpace]deg\[LetterSpace]Ligand_] -> Inactive\[LetterSpace]Rec*index\[LetterSpace]induced\[LetterSpace]ligand\[LetterSpace]deg*kin\[LetterSpace]deg\[LetterSpace]Ligand, Function\[LetterSpace]for\[LetterSpace]reaction17\[LetterSpace]0[Inactive\[LetterSpace]Rec_,index\[LetterSpace]induced\[LetterSpace]R2\[LetterSpace]deg_,kdeg\[LetterSpace]R2_] -> Inactive\[LetterSpace]Rec*index\[LetterSpace]induced\[LetterSpace]R2\[LetterSpace]deg*kdeg\[LetterSpace]R2, NUC\[LetterSpace]CYT\[LetterSpace]SMAD2 -> (pS2\[LetterSpace]n[t] + 3*pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n[t] + 2*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t] + S2\[LetterSpace]n[t])/(pS2\[LetterSpace]c[t] + 3*pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c[t] + 2*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c[t] + S2\[LetterSpace]c[t])
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "", "metabolite" -> "", "extent" -> ""}
};

(* Time evolution *)
odes = {
Active\[LetterSpace]R2'[t] == 1.0*reaction3 +1.0*reaction8 -1.0*reaction4 -1.0*reaction7, Active\[LetterSpace]Rec'[t] == 1.0*reaction7 +1.0*reaction15 -1.0*reaction8 -1.0*reaction14 -1.0*reaction22, Active\[LetterSpace]Rec\[LetterSpace]endo'[t] == 1.0*reaction22 -1.0*reaction9 -1.0*reaction10 -1.0*reaction11 -1.0*reaction23, Inactive\[LetterSpace]Rec'[t] == 1.0*reaction14 -1.0*reaction15 -1.0*reaction16 -1.0*reaction17, S2\[LetterSpace]c'[t] == 1.0*reaction36 +1.0*reaction45 -1.0*reaction26 -1.0*reaction35 -1.0*reaction54, S2\[LetterSpace]n'[t] == 1.0*reaction35 +1.0*reaction43 +1.0*reaction44 -1.0*reaction36 -1.0*reaction55, S4\[LetterSpace]c'[t] == 1.0*reaction28 +1.0*reaction38 +1.0*reaction46 +1.0*reaction60 -1.0*reaction27 -1.0*reaction37 -1.0*reaction58, S4\[LetterSpace]n'[t] == 1.0*reaction30 +1.0*reaction37 +1.0*reaction44 +1.0*reaction63 -1.0*reaction29 -1.0*reaction38 -1.0*reaction59, S7'[t] == 1.0*reaction15 +1.0*reaction16 +1.0*reaction17 +1.0*reaction53 -1.0*reaction14 -1.0*reaction68, Smad7mRNA'[t] == 1.0*reaction52 -1.0*reaction67, Smad7mRNA1'[t] == 1.0*reaction51 -1.0*reaction52 -1.0*reaction66, TGFb'[t] == 1.0*reaction6 -1.0*reaction1 -1.0*reaction25, TGFb\[LetterSpace]In'[t] == 1.0*reaction9 +1.0*reaction10 +1.0*reaction17 +1.0*reaction23 -1.0*reaction24, TGFb\[LetterSpace]R1\[LetterSpace]endo'[t] == 1.0*reaction9 +1.0*reaction11 +1.0*reaction16 +1.0*reaction17 +1.0*reaction18 +1.0*reaction23 -1.0*reaction12 -1.0*reaction20, TGFb\[LetterSpace]R1\[LetterSpace]surface'[t] == 1.0*reaction8 +1.0*reaction20 +1.0*reaction47 -1.0*reaction7 -1.0*reaction18, TGFb\[LetterSpace]R2\[LetterSpace]endo'[t] == 1.0*reaction10 +1.0*reaction11 +1.0*reaction16 +1.0*reaction19 +1.0*reaction23 -1.0*reaction13 -1.0*reaction21, TGFb\[LetterSpace]R2\[LetterSpace]surface'[t] == 1.0*reaction5 +1.0*reaction21 +1.0*reaction48 -1.0*reaction2 -1.0*reaction19, pS2\[LetterSpace]c'[t] == 1.0*reaction26 +2.0*reaction28 +3.0*reaction32 +1.0*reaction40 +1.0*reaction60 +2.0*reaction61 +2.0*reaction62 -2.0*reaction27 -3.0*reaction31 -1.0*reaction39 -1.0*reaction57, pS2\[LetterSpace]n'[t] == 2.0*reaction30 +3.0*reaction34 +1.0*reaction39 +1.0*reaction44 +1.0*reaction63 +2.0*reaction64 +2.0*reaction65 -2.0*reaction29 -3.0*reaction33 -1.0*reaction40 -1.0*reaction43 -1.0*reaction56, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c'[t] == 1.0*reaction27 -1.0*reaction28 -1.0*reaction41 -1.0*reaction60 -1.0*reaction61, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n'[t] == 1.0*reaction29 +1.0*reaction41 -1.0*reaction30 -1.0*reaction44 -1.0*reaction63 -1.0*reaction64, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c'[t] == 1.0*reaction31 -1.0*reaction32 -1.0*reaction42 -1.0*reaction62, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n'[t] == 1.0*reaction33 +1.0*reaction42 -1.0*reaction34 -1.0*reaction65
};
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]}]