(* 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 = {
Ca\[LetterSpace]Cyt[t], Ca\[LetterSpace]ER[t], IP3\[LetterSpace]Cyt[t], O\[LetterSpace]o[t], Oc[t], S2a[t], S4[t], h[t]
};

initialValues = {
h[0] == 0.0, Ca\[LetterSpace]Cyt[0] == 0.0, Ca\[LetterSpace]ER[0] == 0.0, IP3\[LetterSpace]Cyt[0] == 0.0, O\[LetterSpace]o[0] == 0.0, Oc[0] == 0.0, S2a[0] == 0.06, S4[0] == 0.0
};

rates = {
ER\[LetterSpace]Channel, PLC, PMCA, PM\[LetterSpace]Channel, SERCA, act\[LetterSpace]apoSTIM1\[LetterSpace]deg, act\[LetterSpace]apoSTIM1\[LetterSpace]syn, apoSTIM1\[LetterSpace]oligo\[LetterSpace]deg, apoSTIM1\[LetterSpace]oligo\[LetterSpace]syn, close\[LetterSpace]CRAC\[LetterSpace]prod, closed\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]deg, deg, open\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]deg, open\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]prod, open\[LetterSpace]to\[LetterSpace]closed
};

rateEquations = {
ER\[LetterSpace]Channel -> Cytoplasm*(-Ca\[LetterSpace]Cyt[t] + Ca\[LetterSpace]ER[t])*(L + (P\[LetterSpace]IP3R*Ca\[LetterSpace]Cyt[t]^3*h[t]^3*IP3\[LetterSpace]Cyt[t]^3)/((Ka + Ca\[LetterSpace]Cyt[t])^3*(Ki + IP3\[LetterSpace]Cyt[t])^3)), PLC -> (Cytoplasm*V\[LetterSpace]PLC*Ca\[LetterSpace]Cyt[t]^2)/(K\[LetterSpace]PLC^2 + Ca\[LetterSpace]Cyt[t]^2), PMCA -> (Cytoplasm*V\[LetterSpace]PMCA*Ca\[LetterSpace]Cyt[t]^q)/(K\[LetterSpace]PMCA^q + Ca\[LetterSpace]Cyt[t]^q), PM\[LetterSpace]Channel -> Cytoplasm*(Ca\[LetterSpace]ec - Ca\[LetterSpace]Cyt[t])*(V\[LetterSpace]PMleak + k\[LetterSpace]soc*O\[LetterSpace]o[t]), SERCA -> (Cytoplasm*V\[LetterSpace]SERCA*Ca\[LetterSpace]Cyt[t]^p)/(K\[LetterSpace]SERCA^p + Ca\[LetterSpace]Cyt[t]^p), act\[LetterSpace]apoSTIM1\[LetterSpace]deg -> ER*k\[LetterSpace]i*S2a[t], act\[LetterSpace]apoSTIM1\[LetterSpace]syn -> ER*k\[LetterSpace]a*S2, apoSTIM1\[LetterSpace]oligo\[LetterSpace]deg -> ER*kd\[LetterSpace]oligo*S4[t], apoSTIM1\[LetterSpace]oligo\[LetterSpace]syn -> (ER*S2^2*Vs4)/(K2^2 + S2^2), close\[LetterSpace]CRAC\[LetterSpace]prod -> (Orai1^n\[LetterSpace]hill*PM*Vcp)/(Kc^n\[LetterSpace]hill + Orai1^n\[LetterSpace]hill), closed\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]deg -> kdc*PM*Oc[t], deg -> (Cytoplasm*kdeg*Ca\[LetterSpace]Cyt[t]^2*IP3\[LetterSpace]Cyt[t])/(K\[LetterSpace]deg^2 + Ca\[LetterSpace]Cyt[t]^2), open\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]deg -> kdo*PM*O\[LetterSpace]o[t], open\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]prod -> (kop*PM*Oc[t]*S2a[t]^l\[LetterSpace]hill)/(Ko^l\[LetterSpace]hill + S2a[t]^l\[LetterSpace]hill), open\[LetterSpace]to\[LetterSpace]closed -> kod*PM*O\[LetterSpace]o[t]
};

parameters = {
A -> 0.5, Ca\[LetterSpace]ec -> 1500.0, K1 -> 5.0, K2 -> 0.14, K\[LetterSpace]PLC -> 0.12, K\[LetterSpace]PMCA -> 0.45, K\[LetterSpace]SERCA -> 0.15, K\[LetterSpace]deg -> 0.1, Ka -> 0.4, Kc -> 2*^-05, Kd -> 0.4, Ki -> 1.0, Ko -> 0.2, L -> 0.00093, Orai1\[LetterSpace]t -> 0.001, P\[LetterSpace]IP3R -> 66.6, St -> 0.6, V\[LetterSpace]PLC -> 0.5, V\[LetterSpace]PMCA -> 1.0, V\[LetterSpace]PMleak -> 5*^-07, V\[LetterSpace]SERCA -> 1.0, Vcp -> 0.00018, Vs4 -> 0.25, k\[LetterSpace]a -> 4.0, k\[LetterSpace]i -> 6.0, k\[LetterSpace]soc -> 2.3, kd\[LetterSpace]oligo -> 0.8, kdc -> 0.5, kdeg -> 0.5, kdo -> 0.6, kod -> 1.0, kop -> 0.5, l\[LetterSpace]hill -> 1.0, n\[LetterSpace]hill -> 3.0, p -> 2.0, q -> 2.0, r\[LetterSpace]hill -> 4.0, Cytoplasm -> 1.0, ER -> 0.185, PM -> 1.0
};

assignments = {
Orai1 -> Orai1\[LetterSpace]t - r\[LetterSpace]hill*Oc[t] - r\[LetterSpace]hill*O\[LetterSpace]o[t], S2 -> (K1^2*(St - S2a[t]))/(K1^2 + Ca\[LetterSpace]ER[t]^2)
};

events = {

};

speciesAnnotations = {
Ca\[LetterSpace]Cyt[t]->"http://identifiers.org/chebi/CHEBI:29108", Ca\[LetterSpace]Cyt[t]->"http://identifiers.org/kegg.compound/C00076", Ca\[LetterSpace]ER[t]->"http://identifiers.org/chebi/CHEBI:29108", Ca\[LetterSpace]ER[t]->"http://identifiers.org/kegg.compound/C00076", IP3\[LetterSpace]Cyt[t]->"http://identifiers.org/chebi/CHEBI:16595", IP3\[LetterSpace]Cyt[t]->"http://identifiers.org/kegg.compound/C01245"
};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
Ca\[LetterSpace]Cyt'[t] == 1.0*ER\[LetterSpace]Channel +1.0*PM\[LetterSpace]Channel -1.0*SERCA -1.0*PMCA, Ca\[LetterSpace]ER'[t] == 1.0*SERCA -1.0*ER\[LetterSpace]Channel, IP3\[LetterSpace]Cyt'[t] == 1.0*PLC -1.0*deg, O\[LetterSpace]o'[t] == 1.0*open\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]prod -1.0*open\[LetterSpace]to\[LetterSpace]closed -1.0*open\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]deg, Oc'[t] == 1.0*close\[LetterSpace]CRAC\[LetterSpace]prod +1.0*open\[LetterSpace]to\[LetterSpace]closed -1.0*closed\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]deg -1.0*open\[LetterSpace]CRAC\[LetterSpace]channel\[LetterSpace]prod, S2a'[t] == 1.0*act\[LetterSpace]apoSTIM1\[LetterSpace]syn -1.0*act\[LetterSpace]apoSTIM1\[LetterSpace]deg, S4'[t] == 1.0*apoSTIM1\[LetterSpace]oligo\[LetterSpace]syn -1.0*apoSTIM1\[LetterSpace]oligo\[LetterSpace]deg, h'[t] == A*(Kd - (Kd + Ca\[LetterSpace]Cyt[t])*h[t])
};
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]}]