//Created by libAntimony v2.4
model fujita_2010__environment(time_)

  // Variable initializations:
  time_ = ;
end

model fujita_2010__EGF(EGF)

  // Variable initializations:
  EGF = 30.0;
end

model fujita_2010__inhibitor(inhibitor)

  // Variable initializations:
  inhibitor = 15.0;
end

model fujita_2010__pro_EGFR(pro_EGFR)

  // Variable initializations:
  pro_EGFR = 6.81902e-4;
end

model fujita_2010__EGFR(EGFR, EGF_EGFR, EGF, pro_EGFR, inhibitor, EGFR_i, EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, EGFR_turnover, time_)

  // Rate Rules:
  EGFR' = (EGF_binding_kb * EGF_EGFR + inhibitor_binding_kb * EGFR_i + EGFR_turnover * (pro_EGFR - EGFR)) - (EGF_binding_kf * EGF * EGFR + inhibitor_binding_kf * inhibitor * EGFR);

  // Variable initializations:
  EGFR = 68190.000000002;
  EGF_EGFR = ;
  EGF = ;
  pro_EGFR = ;
  inhibitor = ;
  EGFR_i = ;
  EGF_binding_kb = ;
  EGF_binding_kf = ;
  inhibitor_binding_kb = ;
  inhibitor_binding_kf = ;
  EGFR_turnover = ;
  time_ = ;
end

model fujita_2010__EGF_EGFR(EGF_EGFR, EGFR, EGF, inhibitor, EGF_EGFR_i, EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, k1_EGFR_phosphorylation, time_)

  // Rate Rules:
  EGF_EGFR' = (EGF_binding_kf * EGF * EGFR + inhibitor_binding_kb * EGF_EGFR_i) - (EGF_binding_kb * EGF_EGFR + inhibitor_binding_kf * inhibitor * EGF_EGFR + k1_EGFR_phosphorylation * EGF_EGFR);

  // Variable initializations:
  EGF_EGFR = 0.0;
  EGFR = ;
  EGF = ;
  inhibitor = ;
  EGF_EGFR_i = ;
  EGF_binding_kb = ;
  EGF_binding_kf = ;
  inhibitor_binding_kb = ;
  inhibitor_binding_kf = ;
  k1_EGFR_phosphorylation = ;
  time_ = ;
end

model fujita_2010__EGFR_i(EGFR_i, EGF_EGFR_i, EGF, EGFR, inhibitor, EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, EGFR_turnover, time_)

  // Rate Rules:
  EGFR_i' = (EGF_binding_kb * EGF_EGFR_i + inhibitor_binding_kf * inhibitor * EGFR) - (EGF_binding_kf * EGF * EGFR_i + inhibitor_binding_kb * EGFR_i + EGFR_turnover * EGFR_i);

  // Variable initializations:
  EGFR_i = 0.0;
  EGF_EGFR_i = ;
  EGF = ;
  EGFR = ;
  inhibitor = ;
  EGF_binding_kb = ;
  EGF_binding_kf = ;
  inhibitor_binding_kb = ;
  inhibitor_binding_kf = ;
  EGFR_turnover = ;
  time_ = ;
end

model fujita_2010__EGF_EGFR_i(EGF_EGFR_i, EGFR_i, EGF, EGF_EGFR, inhibitor, EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, time_)

  // Rate Rules:
  EGF_EGFR_i' = (EGF_binding_kf * EGF * EGFR_i + inhibitor_binding_kf * inhibitor * EGF_EGFR) - (EGF_binding_kb * EGF_EGFR_i + inhibitor_binding_kb * EGF_EGFR_i);

  // Variable initializations:
  EGF_EGFR_i = 0.0;
  EGFR_i = ;
  EGF = ;
  EGF_EGFR = ;
  inhibitor = ;
  EGF_binding_kb = ;
  EGF_binding_kf = ;
  inhibitor_binding_kb = ;
  inhibitor_binding_kf = ;
  time_ = ;
end

model fujita_2010__pEGFR(pEGFR, EGF_EGFR, pEGFR_Akt, Akt, k1_pEGFR_degradation, k1_Akt_phosphorylation, k1_EGFR_phosphorylation, k1_pEGFR_Akt, k2_pEGFR_Akt, time_)

  // Rate Rules:
  pEGFR' = (k1_EGFR_phosphorylation * EGF_EGFR + k1_Akt_phosphorylation * pEGFR_Akt + k2_pEGFR_Akt * pEGFR_Akt) - (k1_pEGFR_degradation * pEGFR + k1_pEGFR_Akt * pEGFR * Akt);

  // Variable initializations:
  pEGFR = 0.0;
  EGF_EGFR = ;
  pEGFR_Akt = ;
  Akt = ;
  k1_pEGFR_degradation = ;
  k1_Akt_phosphorylation = ;
  k1_EGFR_phosphorylation = ;
  k1_pEGFR_Akt = ;
  k2_pEGFR_Akt = ;
  time_ = ;
end

model fujita_2010__pEGFR_Akt(pEGFR_Akt, pEGFR, Akt, k1_Akt_phosphorylation, k1_pEGFR_Akt, k2_pEGFR_Akt, time_)

  // Rate Rules:
  pEGFR_Akt' = k1_pEGFR_Akt * pEGFR * Akt - (k2_pEGFR_Akt * pEGFR_Akt + k1_Akt_phosphorylation * pEGFR_Akt);

  // Variable initializations:
  pEGFR_Akt = 0.0;
  pEGFR = ;
  Akt = ;
  k1_Akt_phosphorylation = ;
  k1_pEGFR_Akt = ;
  k2_pEGFR_Akt = ;
  time_ = ;
end

model fujita_2010__Akt(Akt, pEGFR, pEGFR_Akt, pAkt, k1_pAkt_dephosphorylation, k1_pEGFR_Akt, k2_pEGFR_Akt, time_)

  // Rate Rules:
  Akt' = (k2_pEGFR_Akt * pEGFR_Akt + k1_pAkt_dephosphorylation * pAkt) - k1_pEGFR_Akt * pEGFR * Akt;

  // Variable initializations:
  Akt = 4.33090e-2;
  pEGFR = ;
  pEGFR_Akt = ;
  pAkt = ;
  k1_pAkt_dephosphorylation = ;
  k1_pEGFR_Akt = ;
  k2_pEGFR_Akt = ;
  time_ = ;
end

model fujita_2010__pAkt(pAkt, S6, pEGFR_Akt, pAkt_S6, k1_Akt_phosphorylation, k1_pAkt_S6, k2_pAkt_S6, k1_S6_phosphorylation, k1_pAkt_dephosphorylation, time_)

  // Rate Rules:
  pAkt' = (k1_Akt_phosphorylation * pEGFR_Akt + k2_pAkt_S6 * pAkt_S6 + k1_S6_phosphorylation * pAkt_S6) - (k1_pAkt_dephosphorylation * pAkt + k1_pAkt_S6 * pAkt * S6);

  // Variable initializations:
  pAkt = 0.0;
  S6 = ;
  pEGFR_Akt = ;
  pAkt_S6 = ;
  k1_Akt_phosphorylation = ;
  k1_pAkt_S6 = ;
  k2_pAkt_S6 = ;
  k1_S6_phosphorylation = ;
  k1_pAkt_dephosphorylation = ;
  time_ = ;
end

model fujita_2010__pAkt_S6(pAkt_S6, S6, pAkt, k1_pAkt_S6, k2_pAkt_S6, k1_S6_phosphorylation, time_)

  // Rate Rules:
  pAkt_S6' = k1_pAkt_S6 * pAkt * S6 - (k2_pAkt_S6 * pAkt_S6 + k1_S6_phosphorylation * pAkt_S6);

  // Variable initializations:
  pAkt_S6 = 0.0;
  S6 = ;
  pAkt = ;
  k1_pAkt_S6 = ;
  k2_pAkt_S6 = ;
  k1_S6_phosphorylation = ;
  time_ = ;
end

model fujita_2010__S6(S6, pAkt_S6, pS6, pAkt, k1_pAkt_S6, k2_pAkt_S6, k1_pS6_dephosphorylation, time_)

  // Rate Rules:
  S6' = (k2_pAkt_S6 * pAkt_S6 + k1_pS6_dephosphorylation * pS6) - k1_pAkt_S6 * pAkt * S6;

  // Variable initializations:
  S6 = 3.54317e0;
  pAkt_S6 = ;
  pS6 = ;
  pAkt = ;
  k1_pAkt_S6 = ;
  k2_pAkt_S6 = ;
  k1_pS6_dephosphorylation = ;
  time_ = ;
end

model fujita_2010__pS6(pS6, pAkt_S6, k1_S6_phosphorylation, k1_pS6_dephosphorylation, time_)

  // Rate Rules:
  pS6' = k1_S6_phosphorylation * pAkt_S6 - k1_pS6_dephosphorylation * pS6;

  // Variable initializations:
  pS6 = 0.0;
  pAkt_S6 = ;
  k1_S6_phosphorylation = ;
  k1_pS6_dephosphorylation = ;
  time_ = ;
end

model fujita_2010__pEGFR_total(pEGFR, pEGFR_Akt)

  // Assignment Rules:
  pEGFR_total := (pEGFR + pEGFR_Akt) * pEGFR_scalefactor;

  // Variable initializations:
  pEGFR_scalefactor = 1.81735e-4;
  pEGFR = ;
  pEGFR_Akt = ;
end

model fujita_2010__pAkt_total(pAkt, pAkt_S6)

  // Assignment Rules:
  pAkt_total := (pAkt + pAkt_S6) * pAkt_scalefactor;

  // Variable initializations:
  pAkt_scalefactor = 6.00588e1;
  pAkt = ;
  pAkt_S6 = ;
end

model fujita_2010__pS6_total(pS6)

  // Assignment Rules:
  pS6_total := pS6 * pS6_scalefactor;

  // Variable initializations:
  pS6_scalefactor = 4.98862e4;
  pS6 = ;
end

model fujita_2010__model_parameters(EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, EGFR_turnover, k1_S6_phosphorylation, k1_pS6_dephosphorylation, k1_pAkt_S6, k2_pAkt_S6, k1_Akt_phosphorylation, k1_pAkt_dephosphorylation, k1_pEGFR_Akt, k2_pEGFR_Akt, k1_pEGFR_degradation, k1_EGFR_phosphorylation)

  // Variable initializations:
  EGF_binding_kb = 4.07490e-2;
  EGF_binding_kf = 6.73816e-3;
  inhibitor_binding_kb = 5.25097e-5;
  inhibitor_binding_kf = 2.434663e-5;
  EGFR_turnover = 1.06386129269658e-4;
  k1_S6_phosphorylation = 1.21498e-3;
  k1_pS6_dephosphorylation = 1.13102e-3;
  k1_pAkt_S6 = 2.10189e-6;
  k2_pAkt_S6 = 5.17940e-15;
  k1_Akt_phosphorylation = 3.05684e-2;
  k1_pAkt_dephosphorylation = 3.27962e-2;
  k1_pEGFR_Akt = 1.55430e-5;
  k2_pEGFR_Akt = 5.17473e-3;
  k1_pEGFR_degradation = 9.97194e-2;
  k1_EGFR_phosphorylation = 1.92391e-2;
end

model *fujita_2010____main()

  // Sub-modules, and any changes to those submodules:
  environment: fujita_2010__environment(time_);
  EGF: fujita_2010__EGF(EGF0);
  inhibitor: fujita_2010__inhibitor(inhibitor0);
  pro_EGFR: fujita_2010__pro_EGFR(pro_EGFR0);
  EGFR: fujita_2010__EGFR(EGFR0, EGF_EGFR0, EGF0, pro_EGFR0, inhibitor0, EGFR_i0, EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, EGFR_turnover, time_);
  EGF_EGFR: fujita_2010__EGF_EGFR(EGF_EGFR0, EGFR0, EGF0, inhibitor0, EGF_EGFR_i0, EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, k1_EGFR_phosphorylation, time_);
  EGFR_i: fujita_2010__EGFR_i(EGFR_i0, EGF_EGFR_i0, EGF0, EGFR0, inhibitor0, EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, EGFR_turnover, time_);
  EGF_EGFR_i: fujita_2010__EGF_EGFR_i(EGF_EGFR_i0, EGFR_i0, EGF0, EGF_EGFR0, inhibitor0, EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, time_);
  pEGFR: fujita_2010__pEGFR(pEGFR0, EGF_EGFR0, pEGFR_Akt0, Akt0, k1_pEGFR_degradation, k1_Akt_phosphorylation, k1_EGFR_phosphorylation, k1_pEGFR_Akt, k2_pEGFR_Akt, time_);
  pEGFR_Akt: fujita_2010__pEGFR_Akt(pEGFR_Akt0, pEGFR0, Akt0, k1_Akt_phosphorylation, k1_pEGFR_Akt, k2_pEGFR_Akt, time_);
  Akt: fujita_2010__Akt(Akt0, pEGFR0, pEGFR_Akt0, pAkt0, k1_pAkt_dephosphorylation, k1_pEGFR_Akt, k2_pEGFR_Akt, time_);
  pAkt: fujita_2010__pAkt(pAkt0, S60, pEGFR_Akt0, pAkt_S60, k1_Akt_phosphorylation, k1_pAkt_S6, k2_pAkt_S6, k1_S6_phosphorylation, k1_pAkt_dephosphorylation, time_);
  pAkt_S6: fujita_2010__pAkt_S6(pAkt_S60, S60, pAkt0, k1_pAkt_S6, k2_pAkt_S6, k1_S6_phosphorylation, time_);
  S6: fujita_2010__S6(S60, pAkt_S60, pS60, pAkt0, k1_pAkt_S6, k2_pAkt_S6, k1_pS6_dephosphorylation, time_);
  pS6: fujita_2010__pS6(pS60, pAkt_S60, k1_S6_phosphorylation, k1_pS6_dephosphorylation, time_);
  pEGFR_total: fujita_2010__pEGFR_total(pEGFR0, pEGFR_Akt0);
  pAkt_total: fujita_2010__pAkt_total(pAkt0, pAkt_S60);
  pS6_total: fujita_2010__pS6_total(pS60);
  model_parameters: fujita_2010__model_parameters(EGF_binding_kb, EGF_binding_kf, inhibitor_binding_kb, inhibitor_binding_kf, EGFR_turnover, k1_S6_phosphorylation, k1_pS6_dephosphorylation, k1_pAkt_S6, k2_pAkt_S6, k1_Akt_phosphorylation, k1_pAkt_dephosphorylation, k1_pEGFR_Akt, k2_pEGFR_Akt, k1_pEGFR_degradation, k1_EGFR_phosphorylation);
end