(* 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] == 1269.5683411637406, S2\[LetterSpace]n[0] == 482.4252172700195, S4\[LetterSpace]c[0] == 632.1908785904236, S4\[LetterSpace]n[0] == 159.99335806446143, S7[0] == 44.316321578389555, Smad7mRNA[0] == 0.10950532934545645, Smad7mRNA1[0] == 0.7124652899456714, TGFb[0] == 100.0, TGFb\[LetterSpace]In[0] == 0.0, TGFb\[LetterSpace]R1\[LetterSpace]endo[0] == 11.105668789597434, TGFb\[LetterSpace]R1\[LetterSpace]surface[0] == 6.412297617178062, TGFb\[LetterSpace]R2\[LetterSpace]endo[0] == 15.542698664608169, TGFb\[LetterSpace]R2\[LetterSpace]surface[0] == 1.9569147929681294, 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[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[Inactive\[LetterSpace]Rec[t], index\[LetterSpace]seq\[LetterSpace]kb, kf\[LetterSpace]Seq\[LetterSpace]S7\[LetterSpace]Rec], reaction16 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction16[Inactive\[LetterSpace]Rec[t], index\[LetterSpace]induced\[LetterSpace]ligand\[LetterSpace]deg, kin\[LetterSpace]deg\[LetterSpace]Ligand], reaction17 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction17[Inactive\[LetterSpace]Rec[t], index\[LetterSpace]induced\[LetterSpace]R2\[LetterSpace]deg, kdeg\[LetterSpace]R2], reaction18 -> cyt*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1*TGFb\[LetterSpace]R1\[LetterSpace]surface[t], reaction19 -> cyt*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1*TGFb\[LetterSpace]R2\[LetterSpace]surface[t], reaction2 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction2[TGFb[t], TGFb\[LetterSpace]R2\[LetterSpace]surface[t], kf\[LetterSpace]R2\[LetterSpace]activation], reaction20 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction20\[LetterSpace]5[TGFb\[LetterSpace]R1\[LetterSpace]endo[t], index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1, k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1], reaction21 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction21\[LetterSpace]5[TGFb\[LetterSpace]R2\[LetterSpace]endo[t], index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1, k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1], reaction22 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction22\[LetterSpace]5[Active\[LetterSpace]Rec[t], index\[LetterSpace]active\[LetterSpace]Rec\[LetterSpace]internalize\[LetterSpace]C1, k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1], 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[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[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\[LetterSpace]0[index\[LetterSpace]kf\[LetterSpace]homotrimer, kf\[LetterSpace]trimmer, pS2\[LetterSpace]c[t]], reaction32 -> nuc*Function\[LetterSpace]for\[LetterSpace]reaction32\[LetterSpace]1\[LetterSpace]0[index\[LetterSpace]kf\[LetterSpace]homotrimer, kf\[LetterSpace]trimmer, pS2\[LetterSpace]n[t]], reaction33 -> nuc*Function\[LetterSpace]for\[LetterSpace]reaction33\[LetterSpace]4[index\[LetterSpace]kb\[LetterSpace]homotrimer, kb\[LetterSpace]trimmer, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n[t]], reaction34 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction34\[LetterSpace]4[index\[LetterSpace]kb\[LetterSpace]homotrimer, kb\[LetterSpace]trimmer, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c[t]], reaction35 -> Function\[LetterSpace]for\[LetterSpace]reaction35\[LetterSpace]6[S2\[LetterSpace]c[t], S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1, cyt], reaction36 -> Function\[LetterSpace]for\[LetterSpace]reaction36\[LetterSpace]6[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1, S2\[LetterSpace]n[t], nuc], reaction37 -> Function\[LetterSpace]for\[LetterSpace]reaction37\[LetterSpace]6[S4\[LetterSpace]c[t], S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1, cyt], reaction38 -> Function\[LetterSpace]for\[LetterSpace]reaction38\[LetterSpace]6[S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1, S4\[LetterSpace]n[t], nuc], reaction39 -> Function\[LetterSpace]for\[LetterSpace]reaction39\[LetterSpace]6[S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1, cyt, pS2\[LetterSpace]c[t]], reaction4 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction4[Active\[LetterSpace]R2[t], index\[LetterSpace]kb\[LetterSpace]R2, kf\[LetterSpace]R2\[LetterSpace]activation], reaction40 -> Function\[LetterSpace]for\[LetterSpace]reaction40\[LetterSpace]6[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1, nuc, pS2\[LetterSpace]n[t]], reaction41 -> Function\[LetterSpace]for\[LetterSpace]reaction41\[LetterSpace]2\[LetterSpace]0[Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc, cyt, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]c[t]], reaction42 -> Function\[LetterSpace]for\[LetterSpace]reaction42[Trimer\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc, cyt, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]c[t]], reaction43 -> k\[LetterSpace]Dephos\[LetterSpace]C1*nuc*pS2\[LetterSpace]n[t], reaction44 -> nuc*Function\[LetterSpace]for\[LetterSpace]reaction44\[LetterSpace]6[index\[LetterSpace]trimer\[LetterSpace]dephos, k\[LetterSpace]Dephos\[LetterSpace]C1, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t]], reaction45 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction45\[LetterSpace]6[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1, S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1, S2\[LetterSpace]total\[LetterSpace]C1, kdeg\[LetterSpace]S2], reaction46 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction46\[LetterSpace]6[S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1, S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1, S4\[LetterSpace]total\[LetterSpace]C1, kdeg\[LetterSpace]S4], reaction47 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction47\[LetterSpace]6[R1\[LetterSpace]total\[LetterSpace]C1, index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1, k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1, kdeg\[LetterSpace]R1], reaction48 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction48\[LetterSpace]3[R2\[LetterSpace]total\[LetterSpace]C1, index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1, k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1, kdeg\[LetterSpace]R2], reaction5 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction5[Active\[LetterSpace]R2[t], index\[LetterSpace]kb\[LetterSpace]R2, kf\[LetterSpace]R2\[LetterSpace]activation], reaction51 -> nuc*Function\[LetterSpace]for\[LetterSpace]reaction51\[LetterSpace]5[K\[LetterSpace]mran, hill\[LetterSpace]fact1, k\[LetterSpace]induced\[LetterSpace]S7\[LetterSpace]production\[LetterSpace]C1, mRNA\[LetterSpace]prod\[LetterSpace]C1, pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n[t]], reaction52 -> Function\[LetterSpace]for\[LetterSpace]reaction52\[LetterSpace]3\[LetterSpace]0[Smad7mRNA1[t], export\[LetterSpace]cytoplasm, nuc], reaction53 -> cyt*Function\[LetterSpace]for\[LetterSpace]reaction53\[LetterSpace]3[Smad7mRNA[t], k\[LetterSpace]S7\[LetterSpace]protein\[LetterSpace]C1], 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[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[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 = {
K\[LetterSpace]mran -> 223.586, R1\[LetterSpace]total\[LetterSpace]C1 -> 17.5179664067755, R2\[LetterSpace]total\[LetterSpace]C1 -> 17.4996134575763, S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1 -> 2.104009010142, S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1 -> 0.399850475956441, S2\[LetterSpace]total\[LetterSpace]C1 -> 1751.99355843376, S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1 -> 0.258219489840727, S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1 -> 0.033282428730513, S4\[LetterSpace]total\[LetterSpace]C1 -> 792.184236654885, 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]C1 -> 0.5069458643727, 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]C1 -> 0.347249922630908, index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1 -> 0.124917902398145, 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]C1 -> 0.0876347503049301, k\[LetterSpace]S7\[LetterSpace]protein\[LetterSpace]C1 -> 0.414415667607617, k\[LetterSpace]disso\[LetterSpace]Active\[LetterSpace]Rec -> 0.00103940411358367, k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1 -> 0.254881953552674, k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1 -> 1.01272286568778, k\[LetterSpace]induced\[LetterSpace]S7\[LetterSpace]production\[LetterSpace]C1 -> 9.21283556540176, 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]C1 -> 0.0254130031757049, cyt -> 2.0, nuc -> 1.0
};

assignments = {
Function\[LetterSpace]for\[LetterSpace]reaction35\[LetterSpace]6[S2\[LetterSpace]c_,S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1_,cyt_] -> cyt*S2\[LetterSpace]c*S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1, Function\[LetterSpace]for\[LetterSpace]reaction38\[LetterSpace]6[S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1_,S4\[LetterSpace]n_,nuc_] -> nuc*S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1*S4\[LetterSpace]n, Function\[LetterSpace]for\[LetterSpace]reaction36\[LetterSpace]6[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1_,S2\[LetterSpace]n_,nuc_] -> nuc*S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1*S2\[LetterSpace]n, Function\[LetterSpace]for\[LetterSpace]reaction1[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]reaction17[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, Function\[LetterSpace]for\[LetterSpace]reaction6[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[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[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[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[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]reaction31\[LetterSpace]2\[LetterSpace]0[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[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]reaction44\[LetterSpace]6[index\[LetterSpace]trimer\[LetterSpace]dephos_,k\[LetterSpace]Dephos\[LetterSpace]C1_,pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n_] -> index\[LetterSpace]trimer\[LetterSpace]dephos*k\[LetterSpace]Dephos\[LetterSpace]C1*pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n, Function\[LetterSpace]for\[LetterSpace]reaction33\[LetterSpace]4[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]reaction22\[LetterSpace]5[Active\[LetterSpace]Rec_,index\[LetterSpace]active\[LetterSpace]Rec\[LetterSpace]internalize\[LetterSpace]C1_,k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1_] -> Active\[LetterSpace]Rec*index\[LetterSpace]active\[LetterSpace]Rec\[LetterSpace]internalize\[LetterSpace]C1*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1, Function\[LetterSpace]for\[LetterSpace]reaction39\[LetterSpace]6[S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1_,cyt_,pS2\[LetterSpace]c_] -> cyt*pS2\[LetterSpace]c*S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1, Function\[LetterSpace]for\[LetterSpace]reaction46\[LetterSpace]6[S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1_,S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1_,S4\[LetterSpace]total\[LetterSpace]C1_,kdeg\[LetterSpace]S4_] -> ((kdeg\[LetterSpace]S4^2 + kdeg\[LetterSpace]S4*S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1 + kdeg\[LetterSpace]S4*S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1)*S4\[LetterSpace]total\[LetterSpace]C1)/(kdeg\[LetterSpace]S4 + S4\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1 + 2*S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1), Function\[LetterSpace]for\[LetterSpace]reaction53\[LetterSpace]3[Smad7mRNA_,k\[LetterSpace]S7\[LetterSpace]protein\[LetterSpace]C1_] -> k\[LetterSpace]S7\[LetterSpace]protein\[LetterSpace]C1*Smad7mRNA, Function\[LetterSpace]for\[LetterSpace]reaction37\[LetterSpace]6[S4\[LetterSpace]c_,S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1_,cyt_] -> cyt*S4\[LetterSpace]c*S4\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1, Function\[LetterSpace]for\[LetterSpace]reaction40\[LetterSpace]6[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1_,nuc_,pS2\[LetterSpace]n_] -> nuc*pS2\[LetterSpace]n*S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1, Function\[LetterSpace]for\[LetterSpace]reaction52\[LetterSpace]3\[LetterSpace]0[Smad7mRNA1_,export\[LetterSpace]cytoplasm_,nuc_] -> export\[LetterSpace]cytoplasm*nuc*Smad7mRNA1, Function\[LetterSpace]for\[LetterSpace]reaction42[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]reaction41\[LetterSpace]2\[LetterSpace]0[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]reaction51\[LetterSpace]5[K\[LetterSpace]mran_,hill\[LetterSpace]fact1_,k\[LetterSpace]induced\[LetterSpace]S7\[LetterSpace]production\[LetterSpace]C1_,mRNA\[LetterSpace]prod\[LetterSpace]C1_,pS2\[LetterSpace]pS2\[LetterSpace]S4\[LetterSpace]n_] -> mRNA\[LetterSpace]prod\[LetterSpace]C1 + (k\[LetterSpace]induced\[LetterSpace]S7\[LetterSpace]production\[LetterSpace]C1*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), Function\[LetterSpace]for\[LetterSpace]reaction47\[LetterSpace]6[R1\[LetterSpace]total\[LetterSpace]C1_,index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1_,k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1_,kdeg\[LetterSpace]R1_] -> (kdeg\[LetterSpace]R1^2*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1*R1\[LetterSpace]total\[LetterSpace]C1)/(kdeg\[LetterSpace]R1^2 + kdeg\[LetterSpace]R1*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1 + index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1*kdeg\[LetterSpace]R1*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1), Function\[LetterSpace]for\[LetterSpace]reaction45\[LetterSpace]6[S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1_,S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1_,S2\[LetterSpace]total\[LetterSpace]C1_,kdeg\[LetterSpace]S2_] -> ((kdeg\[LetterSpace]S2^2 + kdeg\[LetterSpace]S2*S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1 + kdeg\[LetterSpace]S2*S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1)*S2\[LetterSpace]total\[LetterSpace]C1)/(kdeg\[LetterSpace]S2 + S2\[LetterSpace]export\[LetterSpace]from\[LetterSpace]nuc\[LetterSpace]C1 + 2*S2\[LetterSpace]import\[LetterSpace]to\[LetterSpace]nuc\[LetterSpace]C1), Function\[LetterSpace]for\[LetterSpace]reaction20\[LetterSpace]5[TGFb\[LetterSpace]R1\[LetterSpace]endo_,index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1_,k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1_] -> index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]1\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1*k\[LetterSpace]in\[LetterSpace]R1\[LetterSpace]C1*TGFb\[LetterSpace]R1\[LetterSpace]endo, Function\[LetterSpace]for\[LetterSpace]reaction21\[LetterSpace]5[TGFb\[LetterSpace]R2\[LetterSpace]endo_,index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1_,k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1_] -> index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1*TGFb\[LetterSpace]R2\[LetterSpace]endo, Function\[LetterSpace]for\[LetterSpace]reaction34\[LetterSpace]4[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]reaction48\[LetterSpace]3[R2\[LetterSpace]total\[LetterSpace]C1_,index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1_,k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1_,kdeg\[LetterSpace]R2_] -> (kdeg\[LetterSpace]R2^2*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1*R2\[LetterSpace]total\[LetterSpace]C1)/(kdeg\[LetterSpace]R2^2 + kdeg\[LetterSpace]R2*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1 + index\[LetterSpace]k\[LetterSpace]out\[LetterSpace]2\[LetterSpace]relative\[LetterSpace]speed\[LetterSpace]C1*kdeg\[LetterSpace]R2*k\[LetterSpace]in\[LetterSpace]R2\[LetterSpace]C1), Function\[LetterSpace]for\[LetterSpace]reaction32\[LetterSpace]1\[LetterSpace]0[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]reaction8[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[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]reaction2[TGFb_,TGFb\[LetterSpace]R2\[LetterSpace]surface_,kf\[LetterSpace]R2\[LetterSpace]activation_] -> kf\[LetterSpace]R2\[LetterSpace]activation*TGFb*TGFb\[LetterSpace]R2\[LetterSpace]surface, 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*reaction18 +1.0*reaction23 -1.0*reaction12 -1.0*reaction20, TGFb\[LetterSpace]R1\[LetterSpace]surface'[t] == 1.0*reaction8 +1.0*reaction16 +1.0*reaction17 +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*reaction19 +1.0*reaction23 -1.0*reaction13 -1.0*reaction21, TGFb\[LetterSpace]R2\[LetterSpace]surface'[t] == 1.0*reaction5 +1.0*reaction16 +1.0*reaction21 +1.0*reaction48 -1.0*reaction2 -1.0*reaction19, pS2\[LetterSpace]c'[t] == 1.0*reaction26 +2.0*reaction28 +3.0*reaction34 +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*reaction33 +1.0*reaction39 +1.0*reaction44 +1.0*reaction63 +2.0*reaction64 +2.0*reaction65 -2.0*reaction29 -3.0*reaction32 -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*reaction34 -1.0*reaction42 -1.0*reaction62, pS2\[LetterSpace]pS2\[LetterSpace]pS2\[LetterSpace]n'[t] == 1.0*reaction32 +1.0*reaction42 -1.0*reaction33 -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]}]