(* 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 = {
Bid[t], CD95[t], CD95L[t], DISC[t], DISCp55[t], FADD[t], PrER[t], PrER\[LetterSpace]mGFP[t], PrNES[t], PrNES\[LetterSpace]mCherry[t], mCherry[t], mGFP[t], p18[t], p18inactive[t], p30[t], p43[t], p55free[t], tBid[t]
};

initialValues = {
Bid[0] == 224.0, CD95[0] == 12.0, CD95L[0] == 16.6, DISC[0] == 0.0, DISCp55[0] == 0.0, FADD[0] == 90.0, PrER[0] == 0.0, PrER\[LetterSpace]mGFP[0] == 3316.0, PrNES[0] == 0.0, PrNES\[LetterSpace]mCherry[0] == 1909.0, mCherry[0] == 0.0, mGFP[0] == 0.0, p18[0] == 0.0, p18inactive[0] == 0.0, p30[0] == 0.0, p43[0] == 0.0, p55free[0] == 127.0, tBid[0] == 0.0
};

rates = {
reaction\[LetterSpace]1, reaction\[LetterSpace]10, reaction\[LetterSpace]11, reaction\[LetterSpace]12, reaction\[LetterSpace]13, reaction\[LetterSpace]14, reaction\[LetterSpace]15, reaction\[LetterSpace]16, reaction\[LetterSpace]17, reaction\[LetterSpace]18, reaction\[LetterSpace]19, reaction\[LetterSpace]2, reaction\[LetterSpace]3, reaction\[LetterSpace]4, reaction\[LetterSpace]5, reaction\[LetterSpace]6, reaction\[LetterSpace]7, reaction\[LetterSpace]8, reaction\[LetterSpace]9
};

rateEquations = {
reaction\[LetterSpace]1 -> CD95act*cell*kon\[LetterSpace]FADD*FADD[t], reaction\[LetterSpace]10 -> cell*kD374*DISCp55[t], reaction\[LetterSpace]11 -> cell*kD374trans\[LetterSpace]p55*DISCp55[t]*(DISCp55[t] + p30[t]), reaction\[LetterSpace]12 -> cell*kD374trans\[LetterSpace]p43*DISCp55[t]*p43[t], reaction\[LetterSpace]13 -> cell*kD374*p30[t], reaction\[LetterSpace]14 -> cell*kD374trans\[LetterSpace]p55*p30[t]*(DISCp55[t] + p30[t]), reaction\[LetterSpace]15 -> cell*kD374trans\[LetterSpace]p43*p30[t]*p43[t], reaction\[LetterSpace]16 -> cell*kdiss\[LetterSpace]p18*p18[t], reaction\[LetterSpace]17 -> cell*kBid*Bid[t]*(p18[t] + p43[t]), reaction\[LetterSpace]18 -> cell*kD374probe*(p18[t] + p43[t])*PrNES\[LetterSpace]mCherry[t], reaction\[LetterSpace]19 -> cell*kD374probe*p18[t]*PrER\[LetterSpace]mGFP[t], reaction\[LetterSpace]2 -> cell*koff\[LetterSpace]FADD*DISC[t], reaction\[LetterSpace]3 -> cell*kDISC*DISC[t]*p55free[t], reaction\[LetterSpace]4 -> cell*kD216*DISCp55[t], reaction\[LetterSpace]5 -> cell*kD216trans\[LetterSpace]p55*DISCp55[t]*(DISCp55[t] + p30[t]), reaction\[LetterSpace]6 -> cell*kD216trans\[LetterSpace]p43*DISCp55[t]*p43[t], reaction\[LetterSpace]7 -> cell*kD216*p43[t], reaction\[LetterSpace]8 -> cell*kD216trans\[LetterSpace]p55*(DISCp55[t] + p30[t])*p43[t], reaction\[LetterSpace]9 -> cell*kD216trans\[LetterSpace]p43*p43[t]^2
};

parameters = {
KDL -> 30.0060394758199, KDR -> 57.2050013008496, kBid -> 0.00052134055139547, kD216 -> 0.00639775937416746, kD216trans\[LetterSpace]p43 -> 5.29906975294056*^-05, kD216trans\[LetterSpace]p55 -> 0.000223246421372882, kD374 -> 0.000644612643975149, kD374probe -> 0.00153710001025539, kD374trans\[LetterSpace]p43 -> 0.00413530054938906, kD374trans\[LetterSpace]p55 -> 0.000543518631342483, kDISC -> 0.000364965874405544, kdiss\[LetterSpace]p18 -> 0.064713651554491, koff\[LetterSpace]FADD -> 0.00130854998177646, kon\[LetterSpace]FADD -> 0.00108871858684363, cell -> 1.0
};

assignments = {
CD95act -> (KDL^2*CD95[t]^3*CD95L[t])/((KDL + CD95L[t])*(KDL^2*KDR + KDL^2*CD95[t]^2 + 2*KDL*KDR*CD95L[t] + KDR*CD95L[t]^2))
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
Bid'[t] ==  -1.0*reaction\[LetterSpace]17, CD95'[t] == 0.0 , CD95L'[t] == 0.0 , DISC'[t] == 1.0*reaction\[LetterSpace]1 +1.0*reaction\[LetterSpace]7 +1.0*reaction\[LetterSpace]8 +1.0*reaction\[LetterSpace]9 +1.0*reaction\[LetterSpace]13 +1.0*reaction\[LetterSpace]14 +1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]2 -1.0*reaction\[LetterSpace]3, DISCp55'[t] == 1.0*reaction\[LetterSpace]3 -1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]10 -1.0*reaction\[LetterSpace]11 -1.0*reaction\[LetterSpace]12, FADD'[t] == 1.0*reaction\[LetterSpace]2 -1.0*reaction\[LetterSpace]1, PrER'[t] == 1.0*reaction\[LetterSpace]19 , PrER\[LetterSpace]mGFP'[t] ==  -1.0*reaction\[LetterSpace]19, PrNES'[t] == 1.0*reaction\[LetterSpace]18 , PrNES\[LetterSpace]mCherry'[t] ==  -1.0*reaction\[LetterSpace]18, mCherry'[t] == 1.0*reaction\[LetterSpace]18 , mGFP'[t] == 1.0*reaction\[LetterSpace]19 , p18'[t] == 1.0*reaction\[LetterSpace]7 +1.0*reaction\[LetterSpace]8 +1.0*reaction\[LetterSpace]9 +1.0*reaction\[LetterSpace]13 +1.0*reaction\[LetterSpace]14 +1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]16, p18inactive'[t] == 1.0*reaction\[LetterSpace]16 , p30'[t] == 1.0*reaction\[LetterSpace]4 +1.0*reaction\[LetterSpace]5 +1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]13 -1.0*reaction\[LetterSpace]14 -1.0*reaction\[LetterSpace]15, p43'[t] == 1.0*reaction\[LetterSpace]10 +1.0*reaction\[LetterSpace]11 +1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]7 -1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]9, p55free'[t] ==  -1.0*reaction\[LetterSpace]3, tBid'[t] == 1.0*reaction\[LetterSpace]17 
};
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]}]