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

  // Variable initializations:
  time_ = ;
end

model leloup_2004__MP(MP, vstot, n, BN, time_)

  // Assignment Rules:
  vsP := vstot;

  // Rate Rules:
  MP' = vsP * (power(BN, n) / (power(KAP, n) + power(BN, n))) - (vmP * (MP / (KmP + MP)) + kdmp * MP);

  // Variable initializations:
  MP = 0.5;
  vmP = 1.1;
  kdmp = 0.01;
  KAP = 0.7;
  KmP = 0.3;
  vstot = ;
  n = ;
  BN = ;
  time_ = ;
end

model leloup_2004__MC(MC, vstot, n, BN, time_)

  // Assignment Rules:
  vsC := 0.8 * vstot;

  // Rate Rules:
  MC' = vsC * (power(BN, n) / (power(KAC, n) + power(BN, n))) - (vmC * (MC / (KmC + MC)) + kdmc * MC);

  // Variable initializations:
  MC = 0.3;
  vmC = 1.0;
  kdmc = 0.01;
  KAC = 1.0;
  KmC = 0.4;
  vstot = ;
  n = ;
  BN = ;
  time_ = ;
end

model leloup_2004__MB(MB, vstot, m, BN, time_)

  // Assignment Rules:
  vsB := 0.7 * vstot;

  // Rate Rules:
  MB' = vsB * (power(KIB, m) / (power(KIB, m) + power(BN, m))) - (vmB * (MB / (KmB + MB)) + kdmb * MB);

  // Variable initializations:
  MB = 3.1;
  vmB = 0.2;
  kdmb = 0.01;
  KIB = 0.8;
  KmB = 0.4;
  vstot = ;
  m = ;
  BN = ;
  time_ = ;
end

model leloup_2004__PC(PC, ksP, Kp, Kdp, k3, k4, kdn, V1P, V2P, PCP, PCC, CC, MP, time_)

  // Rate Rules:
  PC' = (ksP * MP + V2P * (PCP / (Kdp + PCP)) + k4 * PCC) - (V1P * (PC / (Kp + PC)) + k3 * PC * CC + kdn * PC);

  // Variable initializations:
  PC = 0.1;
  ksP = ;
  Kp = ;
  Kdp = ;
  k3 = ;
  k4 = ;
  kdn = ;
  V1P = ;
  V2P = ;
  PCP = ;
  PCC = ;
  CC = ;
  MP = ;
  time_ = ;
end

model leloup_2004__CC(CC, ksC, Kp, Kdp, k3, k4, kdnc, V1C, V2C, CCP, PC, PCC, MC, time_)

  // Rate Rules:
  CC' = (ksC * MC + V2C * (CCP / (Kdp + CCP)) + k4 * PCC) - (V1C * (CC / (Kp + CC)) + k3 * PC * CC + kdnc * CC);

  // Variable initializations:
  CC = 0.1;
  ksC = ;
  Kp = ;
  Kdp = ;
  k3 = ;
  k4 = ;
  kdnc = ;
  V1C = ;
  V2C = ;
  CCP = ;
  PC = ;
  PCC = ;
  MC = ;
  time_ = ;
end

model leloup_2004__PCP(PCP, vdPC, Kp, Kdp, Kd, kdn, V1P, V2P, PC, time_)

  // Rate Rules:
  PCP' = V1P * (PC / (Kp + PC)) - (V2P * (PCP / (Kdp + PCP)) + vdPC * (PCP / (Kd + PCP)) + kdn * PCP);

  // Variable initializations:
  PCP = 0.1;
  vdPC = ;
  Kp = ;
  Kdp = ;
  Kd = ;
  kdn = ;
  V1P = ;
  V2P = ;
  PC = ;
  time_ = ;
end

model leloup_2004__CCP(CCP, vdCC, Kp, Kdp, Kd, kdn, V1C, V2C, CC, time_)

  // Rate Rules:
  CCP' = V1C * (CC / (Kp + CC)) - (V2C * (CCP / (Kdp + CCP)) + vdCC * (CCP / (Kd + CCP)) + kdn * CCP);

  // Variable initializations:
  CCP = 0.1;
  vdCC = ;
  Kp = ;
  Kdp = ;
  Kd = ;
  kdn = ;
  V1C = ;
  V2C = ;
  CC = ;
  time_ = ;
end

model leloup_2004__PCC(PCC, Kp, Kdp, k1, k2, k3, k4, kdn, V1PC, V2PC, PCCP, PCN, PC, CC, time_)

  // Rate Rules:
  PCC' = (V2PC * (PCCP / (Kdp + PCCP)) + k3 * PC * CC + k2 * PCN) - (V1PC * (PCC / (Kp + PCC)) + k4 * PCC + k1 * PCC + kdn * PCC);

  // Variable initializations:
  PCC = 0.1;
  Kp = ;
  Kdp = ;
  k1 = ;
  k2 = ;
  k3 = ;
  k4 = ;
  kdn = ;
  V1PC = ;
  V2PC = ;
  PCCP = ;
  PCN = ;
  PC = ;
  CC = ;
  time_ = ;
end

model leloup_2004__PCN(PCN, Kp, Kdp, k1, k2, k7, k8, kdn, V3PC, V4PC, PCNP, PCC, BN, IN_, time_)

  // Rate Rules:
  PCN' = (V4PC * (PCNP / (Kdp + PCNP)) + k1 * PCC + k8 * IN_) - (V3PC * (PCN / (Kp + PCN)) + k2 * PCN + k7 * BN * PCN + kdn * PCN);

  // Variable initializations:
  PCN = 0.1;
  Kp = ;
  Kdp = ;
  k1 = ;
  k2 = ;
  k7 = ;
  k8 = ;
  kdn = ;
  V3PC = ;
  V4PC = ;
  PCNP = ;
  PCC = ;
  BN = ;
  IN_ = ;
  time_ = ;
end

model leloup_2004__PCCP(PCCP, Kp, Kdp, Kd, kdn, V1PC, V2PC, vdPCC, PCC, time_)

  // Rate Rules:
  PCCP' = V1PC * (PCC / (Kp + PCC)) - (V2PC * (PCCP / (Kdp + PCCP)) + vdPCC * (PCCP / (Kd + PCCP)) + kdn * PCCP);

  // Variable initializations:
  PCCP = 0.1;
  Kp = ;
  Kdp = ;
  Kd = ;
  kdn = ;
  V1PC = ;
  V2PC = ;
  vdPCC = ;
  PCC = ;
  time_ = ;
end

model leloup_2004__PCNP(PCNP, Kp, Kdp, Kd, kdn, V3PC, V4PC, vdPCN, PCN, time_)

  // Rate Rules:
  PCNP' = V3PC * (PCN / (Kp + PCN)) - (V4PC * (PCNP / (Kdp + PCNP)) + vdPCN * (PCNP / (Kd + PCNP)) + kdn * PCNP);

  // Variable initializations:
  PCNP = 0.1;
  Kp = ;
  Kdp = ;
  Kd = ;
  kdn = ;
  V3PC = ;
  V4PC = ;
  vdPCN = ;
  PCN = ;
  time_ = ;
end

model leloup_2004__BC(BC, Kp, Kdp, ksB, k5, k6, kdn, V1B, V2B, MB, BCP, BN, time_)

  // Rate Rules:
  BC' = (V2B * (BCP / (Kdp + BCP)) + k6 * BN + ksB * MB) - (V1B * (BC / (Kp + BC)) + k5 * BC + kdn * BC);

  // Variable initializations:
  BC = 0.1;
  Kp = ;
  Kdp = ;
  ksB = ;
  k5 = ;
  k6 = ;
  kdn = ;
  V1B = ;
  V2B = ;
  MB = ;
  BCP = ;
  BN = ;
  time_ = ;
end

model leloup_2004__BCP(BCP, Kp, Kdp, Kd, kdn, V1B, V2B, vdBC, BC, time_)

  // Rate Rules:
  BCP' = V1B * (BC / (Kp + BC)) - (V2B * (BCP / (Kdp + BCP)) + vdBC * (BCP / (Kd + BCP)) + kdn * BCP);

  // Variable initializations:
  BCP = 0.1;
  Kp = ;
  Kdp = ;
  Kd = ;
  kdn = ;
  V1B = ;
  V2B = ;
  vdBC = ;
  BC = ;
  time_ = ;
end

model leloup_2004__BN(BN, Kp, Kdp, k5, k6, k7, k8, kdn, V3B, V4B, BNP, BC, PCN, IN_, time_)

  // Rate Rules:
  BN' = (V4B * (BNP / (Kdp + BNP)) + k5 * BC + k8 * IN_) - (V3B * (BN / (Kp + BN)) + k6 * BN + k7 * BN * PCN + kdn * BN);

  // Variable initializations:
  BN = 0.1;
  Kp = ;
  Kdp = ;
  k5 = ;
  k6 = ;
  k7 = ;
  k8 = ;
  kdn = ;
  V3B = ;
  V4B = ;
  BNP = ;
  BC = ;
  PCN = ;
  IN_ = ;
  time_ = ;
end

model leloup_2004__BNP(BNP, Kp, Kdp, Kd, kdn, V3B, V4B, vdBN, BN, time_)

  // Rate Rules:
  BNP' = V3B * (BN / (Kp + BN)) - (V4B * (BNP / (Kdp + BNP)) + vdBN * (BNP / (Kd + BNP)) + kdn * BNP);

  // Variable initializations:
  BNP = 0.1;
  Kp = ;
  Kdp = ;
  Kd = ;
  kdn = ;
  V3B = ;
  V4B = ;
  vdBN = ;
  BN = ;
  time_ = ;
end

model leloup_2004__IN(IN_, Kd, kdn, k7, k8, vdIN, PCN, BN, time_)

  // Rate Rules:
  IN_' = k7 * BN * PCN - (k8 * IN_ + vdIN * (IN_ / (Kd + IN_)) + kdn * IN_);

  // Variable initializations:
  IN_ = 0.1;
  Kd = ;
  kdn = ;
  k7 = ;
  k8 = ;
  vdIN = ;
  PCN = ;
  BN = ;
  time_ = ;
end

model leloup_2004__model_parameters(k1, k2, k3, k4, k5, k6, k7, k8, kdnc, kdn, ksB, ksC, ksP, m, n, vstot, Kd, Kdp, Kp, V1B, V1C, V1P, V1PC, V2B, V2C, V2P, V2PC, V3B, V3PC, V4B, V4PC, vdBC, vdBN, vdCC, vdIN, vdPC, vdPCC, vdPCN)

  // Assignment Rules:
  ksB := kstot;
  ksC := kstot;
  ksP := 0.5 * kstot;
  V1P := Vphos;
  V1PC := Vphos;
  V3PC := Vphos;

  // Variable initializations:
  k1 = 0.8;
  k2 = 0.2;
  k3 = 0.8;
  k4 = 0.2;
  k5 = 0.4;
  k6 = 0.2;
  k7 = 0.5;
  k8 = 0.1;
  kdnc = 0.01;
  kdn = 0.01;
  kstot = 1.0;
  m = 4.0;
  n = 4.0;
  vstot = 1.0;
  Kd = 0.3;
  Kdp = 0.3;
  Kp = 0.1;
  Vphos = 0.6;
  V1B = 1.0;
  V1C = 0.6;
  V2B = 0.1;
  V2C = 0.1;
  V2P = 0.3;
  V2PC = 0.1;
  V3B = 1.0;
  V4B = 0.2;
  V4PC = 0.1;
  vdBC = 1.0;
  vdBN = 0.5;
  vdCC = 0.7;
  vdIN = 0.8;
  vdPC = 0.7;
  vdPCC = 1.0;
  vdPCN = 1.0;
end

model *leloup_2004____main()

  // Sub-modules, and any changes to those submodules:
  environment: leloup_2004__environment(time_);
  MP: leloup_2004__MP(MP0, vstot, n, BN0, time_);
  MC: leloup_2004__MC(MC0, vstot, n, BN0, time_);
  MB: leloup_2004__MB(MB0, vstot, m, BN0, time_);
  PC: leloup_2004__PC(PC0, ksP, Kp, Kdp, k3, k4, kdn, V1P, V2P, PCP0, PCC0, CC0, MP0, time_);
  CC: leloup_2004__CC(CC0, ksC, Kp, Kdp, k3, k4, kdnc, V1C, V2C, CCP0, PC0, PCC0, MC0, time_);
  PCP: leloup_2004__PCP(PCP0, vdPC, Kp, Kdp, Kd, kdn, V1P, V2P, PC0, time_);
  CCP: leloup_2004__CCP(CCP0, vdCC, Kp, Kdp, Kd, kdn, V1C, V2C, CC0, time_);
  PCC: leloup_2004__PCC(PCC0, Kp, Kdp, k1, k2, k3, k4, kdn, V1PC, V2PC, PCCP0, PCN0, PC0, CC0, time_);
  PCN: leloup_2004__PCN(PCN0, Kp, Kdp, k1, k2, k7, k8, kdn, V3PC, V4PC, PCNP0, PCC0, BN0, IN_0, time_);
  PCCP: leloup_2004__PCCP(PCCP0, Kp, Kdp, Kd, kdn, V1PC, V2PC, vdPCC, PCC0, time_);
  PCNP: leloup_2004__PCNP(PCNP0, Kp, Kdp, Kd, kdn, V3PC, V4PC, vdPCN, PCN0, time_);
  BC: leloup_2004__BC(BC0, Kp, Kdp, ksB, k5, k6, kdn, V1B, V2B, MB0, BCP0, BN0, time_);
  BCP: leloup_2004__BCP(BCP0, Kp, Kdp, Kd, kdn, V1B, V2B, vdBC, BC0, time_);
  BN: leloup_2004__BN(BN0, Kp, Kdp, k5, k6, k7, k8, kdn, V3B, V4B, BNP0, BC0, PCN0, IN_0, time_);
  BNP: leloup_2004__BNP(BNP0, Kp, Kdp, Kd, kdn, V3B, V4B, vdBN, BN0, time_);
  IN_: leloup_2004__IN(IN_0, Kd, kdn, k7, k8, vdIN, PCN0, BN0, time_);
  model_parameters: leloup_2004__model_parameters(k1, k2, k3, k4, k5, k6, k7, k8, kdnc, kdn, ksB, ksC, ksP, m, n, vstot, Kd, Kdp, Kp, V1B, V1C, V1P, V1PC, V2B, V2C, V2P, V2PC, V3B, V3PC, V4B, V4PC, vdBC, vdBN, vdCC, vdIN, vdPC, vdPCC, vdPCN);
end