#mayya1



#reactions
PstatDimer__import:
 {1.0}PstatDimer_sol = {1.0}PstatDimer_nuc
  PstatDimer_impMax * PstatDimer_sol * (1 / (Kpsd_imp + PstatDimer_sol)) * nm
PstatDimerisation:
 {2.0}Pstat_sol = {1.0}PstatDimer_sol
  (Kf_PstatDimerisation * pow(Pstat_sol, 2) + -(Kr_PstatDimerisation * PstatDimer_sol)) * sol
PstatDimerisationNuc:
 {2.0}Pstat_nuc = {1.0}PstatDimer_nuc
  (Kf_PstatDimerisation * pow(Pstat_nuc, 2) + -(Kr_PstatDimerisation * PstatDimer_nuc)) * nuc
statDephosphorylation:
 {1.0}Pstat_nuc = {1.0}stat_nuc
  Kcat_dephos * statPhosphatase_nuc * Pstat_nuc * (1 / (Km_dephos + Pstat_nuc)) * nuc
statPhosphorylation:
 {1.0}stat_sol = {1.0}Pstat_sol + {1.0}species_test
  Kcat_phos * statKinase_sol * stat_sol * (1 / (Km_phos + stat_sol)) * sol
stat_export:
 {1.0}stat_sol = {1.0}stat_nuc
  nuc * stat_expMax * stat_nuc * (1 / (Ks_exp + stat_nuc)) * nm
stat_import:
 {1.0}stat_sol = {1.0}stat_nuc
  nuc * stat_impMax * stat_sol * (1 / (Ks_imp + stat_sol)) * nm


#parameters
PstatDimerisation_Kf_PstatDimerisation = 0.6
PstatDimerisation_Kr_PstatDimerisation = 0.03
statDephosphorylation_Kcat_dephos = 1.0
statDephosphorylation_Km_dephos = 2.0
statPhosphorylation_Kcat_phos = 1.0
statPhosphorylation_Km_phos = 4.0
PstatDimerisationNuc_Kf_PstatDimerisation = 0.6
PstatDimerisationNuc_Kr_PstatDimerisation = 0.03
PstatDimer__import_PstatDimer_impMax = 0.045
PstatDimer__import_Kpsd_imp = 0.3
stat_export_stat_expMax = -0.06
stat_export_Ks_exp = 0.6
stat_import_stat_impMax = 0.003
stat_import_Ks_imp = 3.0
nm = 1.0
nuc = 1.0
sol = 14.625


#initial values
PstatDimer_nuc = 0.0
PstatDimer_sol = 0.0
Pstat_nuc = 0.0
Pstat_sol = 0.0
species_test = 0.0
statPhosphatase_nuc = 0.05
stat_nuc = 0.0
stat_sol = 1.0


#assignment rules
!F statKinase_sol = piecewise(0.01 * sin(0.001571 * (-500.0 + t)), and(gt(t, 500.0), lt(t, 2502.54614894971)), 0.0)