// Created by libAntimony v2.8.0
model *Conant2007_glycolysis_2C()

  // Compartments and Species:
  compartment cyto, mito;
  species GLCi in cyto, ATP in cyto, G6P in cyto, ADP in cyto, F6P in cyto;
  species F16bP in cyto, $F26bP in cyto, AMP in cyto, DHAP in cyto, GAP in cyto;
  species NAD in cyto, BPG in cyto, NADH in cyto, P3G in cyto, P2G in cyto;
  species PEP in cyto, PYR in cyto, PYRmito in mito, AcAld in cyto, $GLCo in cyto;
  species $CO2 in cyto, $EtOH in cyto, $Glycerol in cyto, $Glycogen in cyto;
  species $Trehalose in cyto, $Succinate in cyto, $AcCoA in mito, $CO2mito in mito;

  // Reactions:
  HXT: $GLCo -> GLCi; (HXT_Vmax_1*WGD_E*(GLCo - GLCi))/HXT_Kglc_1/(1 + (GLCo + GLCi)/HXT_Kglc_1 + (HXT_Ki_1*GLCo*GLCi)/HXT_Kglc_1^2);
  HXK: GLCi + ATP -> G6P + ADP; (cyto*WGD_E*HXK_Vmax_2*((GLCi*ATP)/(HXK_Kglc_2*HXK_Katp_2) - (G6P*ADP)/(HXK_Kglc_2*HXK_Katp_2*HXK_Keq_2)))/((1 + GLCi/HXK_Kglc_2 + G6P/HXK_Kg6p_2)*(1 + ATP/HXK_Katp_2 + ADP/HXK_Kadp_2));
  PGI: G6P -> F6P; (cyto*PGI_Vmax_3*WGD_E*(G6P/PGI_Kg6p_3 - F6P/(PGI_Kg6p_3*PGI_Keq_3)))/(1 + G6P/PGI_Kg6p_3 + F6P/PGI_Kf6p_3);
  PFK: ATP + F6P => ADP + F16bP; cyto*PFK_Vmax_4*WGD_E*((PFK_gR_4*(F6P/PFK_Kf6p_4)*(ATP/PFK_Katp_4)*(1 + F6P/PFK_Kf6p_4 + ATP/PFK_Katp_4 + ((PFK_gR_4*F6P)/PFK_Kf6p_4*ATP)/PFK_Katp_4))/((1 + F6P/PFK_Kf6p_4 + ATP/PFK_Katp_4 + ((PFK_gR_4*F6P)/PFK_Kf6p_4*ATP)/PFK_Katp_4)^2 + PFK_L0_4*((1 + (PFK_Ciatp_4*ATP)/PFK_Kiatp_4)/(1 + ATP/PFK_Kiatp_4))^2*((1 + (PFK_Camp_4*AMP)/PFK_Kamp_4)/(1 + AMP/PFK_Kamp_4))^2*((1 + (PFK_Cf26_4*F26bP)/PFK_Kf26_4 + (PFK_Cf16_4*F16bP)/PFK_Kf16_4)/(1 + F26bP/PFK_Kf26_4 + F16bP/PFK_Kf16_4))^2*(1 + (PFK_Catp_4*ATP)/PFK_Katp_4)^2));
  FBA: F16bP -> DHAP + GAP; (cyto*FBA_Vmax_5*WGD_E*(F16bP/FBA_Kf16bp_5 - (DHAP*GAP)/(FBA_Kf16bp_5*FBA_Keq_5)))/(1 + F16bP/FBA_Kf16bp_5 + DHAP/FBA_Kdhap_5 + GAP/FBA_Kgap_5 + (F16bP*GAP)/(FBA_Kf16bp_5*FBA_Kigap_5) + (DHAP*GAP)/(FBA_Kdhap_5*FBA_Kgap_5));
  TPI: DHAP -> GAP; cyto*(TPI_k1_6*DHAP - TPI_k2_6*GAP);
  TDH: GAP + NAD -> BPG + NADH; (cyto*TDH_C_7*((TDH_Vmaxf_7*WGD_E*GAP*NAD)/(TDH_Kgap_7*TDH_Knad_7) - (TDH_Vmaxr_7*WGD_E*BPG*NADH)/(TDH_Kbpg_7*TDH_Knadh_7)))/((1 + GAP/TDH_Kgap_7 + BPG/TDH_Kbpg_7)*(1 + NAD/TDH_Knad_7 + NADH/TDH_Knadh_7));
  PGK: ADP + BPG -> ATP + P3G; (cyto*PGK_Vmax_8*WGD_E*((PGK_Keq_8*BPG*ADP - P3G*ATP)/(PGK_Kp3g_8*PGK_Katp_8)))/((1 + BPG/PGK_Kbpg_8 + P3G/PGK_Kp3g_8)*(1 + ADP/PGK_Kadp_8 + ATP/PGK_Katp_8));
  PGM: P3G -> P2G; (cyto*PGM_Vmax_9*WGD_E*(P3G/PGM_Kp3g_9 - P2G/(PGM_Kp3g_9*PGM_Keq_9)))/(1 + P3G/PGM_Kp3g_9 + P2G/PGM_Kp2g_9);
  ENO: P2G -> PEP; (cyto*ENO_Vmax_10*WGD_E*(P2G/ENO_Kp2g_10 - PEP/(ENO_Kp2g_10*ENO_Keq_10)))/(1 + P2G/ENO_Kp2g_10 + PEP/ENO_Kpep_10);
  PYK: ADP + PEP -> ATP + PYR; (cyto*PYK_Vmax_11*WGD_E*((PEP*ADP)/(PYK_Kpep_11*PYK_Kadp_11) - (PYR*ATP)/(PYK_Kpep_11*PYK_Kadp_11*PYK_Keq_11)))/((1 + PEP/PYK_Kpep_11 + PYR/PYK_Kpyr_11)*(1 + ADP/PYK_Kadp_11 + ATP/PYK_Katp_11));
  PDC: PYR => AcAld + $CO2; (cyto*PDC_Vmax_12*WGD_E*(PYR/PDC_Kpyr_12)^PDC_nH_12)/(1 + (PYR/PDC_Kpyr_12)^PDC_nH_12);
  ADH: NAD + $EtOH -> NADH + AcAld; (cyto*ADH_Vmax_13*WGD_E*((EtOH*NAD)/(ADH_Ketoh_13*ADH_Kinad_13) - (AcAld*NADH)/(ADH_Ketoh_13*ADH_Kinad_13*ADH_Keq_13)))/(1 + NAD/ADH_Kinad_13 + (EtOH*ADH_Knad_13)/(ADH_Kinad_13*ADH_Ketoh_13) + (AcAld*ADH_Knadh_13)/(ADH_Kinadh_13*ADH_Kacald_13) + NADH/ADH_Kinadh_13 + (EtOH*NAD)/(ADH_Kinad_13*ADH_Ketoh_13) + (NAD*AcAld*ADH_Knadh_13)/(ADH_Kinad_13*ADH_Kinadh_13*ADH_Kacald_13) + (EtOH*NADH*ADH_Knad_13)/(ADH_Kinad_13*ADH_Kinadh_13*ADH_Ketoh_13) + (AcAld*NADH)/(ADH_Kacald_13*ADH_Kinadh_13) + (EtOH*NAD*AcAld)/(ADH_Kinad_13*ADH_Kiacald_13*ADH_Ketoh_13) + (EtOH*AcAld*NADH)/(ADH_Kietoh_13*ADH_Kinadh_13*ADH_Kacald_13));
  ATPase: ATP => ADP; cyto*ATPase_Katpase_14*ATP;
  AK: 2ADP -> ATP + AMP; cyto*(AK_k1_15*ADP*ADP - AK_k2_15*ATP*AMP);
  G3PDH: DHAP + NADH -> NAD + $Glycerol; (cyto*G3PDH_Vmax_16*WGD_E*(DHAP/G3PDH_Kdhap_16*(NADH/G3PDH_Knadh_16) - Glycerol/G3PDH_Kdhap_16*(NAD/G3PDH_Knadh_16)*(1/G3PDH_Keq_16)))/((1 + DHAP/G3PDH_Kdhap_16 + Glycerol/G3PDH_Kglycerol_16)*(1 + NADH/G3PDH_Knadh_16 + NAD/G3PDH_Knad_16));
  glycogen_branch: ATP + G6P => ADP + $Glycogen; cyto*glycogen_branch_KGLYCOGEN_17;
  trehalose_synth: ATP + 2G6P => ADP + $Trehalose; cyto*trehalose_synth_Ktrehalose_18;
  succinate_syn: 3NAD + 2AcAld => 3NADH + $Succinate; cyto*succinate_syn_k_19*AcAld;
  PYR_shut: PYR -> PYRmito; t_m*(PYR_shut_k1*PYR - PYR_shut_k2*PYRmito);
  PDH: PYRmito => $AcCoA + $CO2mito; (mito*WGD_E*PDH_Vmax_PDH*PYRmito*(PDH_NADX_tot - PDH_NADX_tot/(1 + NAD/NADH)))/((PDH_NADX_tot*PDH_Ki_PYR*PDH_K_NAD)/PDH_Ki_NADH/(1 + NAD/NADH) + PDH_K_PYR*(PDH_NADX_tot - PDH_NADX_tot/(1 + NAD/NADH)) + PDH_K_NAD*PYRmito + ((PDH_NADX_tot*PDH_K_NAD)/PDH_Ki_NADH*PYRmito)/(1 + NAD/NADH) + (PDH_NADX_tot - PDH_NADX_tot/(1 + NAD/NADH))*PYRmito);

  // Species initializations:
  GLCi = 0.097652231064563;
  ATP = 2.52512746499271;
  G6P = 2.67504014044787;
  ADP = 1.28198768168719;
  F6P = 0.624976405532373;
  F16bP = 6.22132076069411;
  F26bP = 0.02;
  AMP = 0.292884853320091;
  DHAP = 1.00415254899644;
  GAP = 0.0451809175780963;
  NAD = 1.50329030201531;
  BPG = 0.000736873499865602;
  NADH = 0.0867096979846952;
  P3G = 0.885688538360659;
  P2G = 0.127695817386632;
  PEP = 0.0632352144936527;
  PYR = 1.81531251192736;
  PYRmito = 0;
  AcAld = 0.178140579850657;
  GLCo = 50;
  CO2 = 1;
  EtOH = 50;
  Glycerol = 0.15;
  Glycogen = 0;
  Trehalose = 0;
  Succinate = 0;
  AcCoA = 1;
  CO2mito = 1;

  // Compartment initializations:
  cyto = 1;
  mito = 1;

  // Variable initializations:
  t_m = 1;
  t_m has dimensionless;
  WGD_E = 0.65;
  WGD_E has dimensionless;
  HXT_Vmax_1 = 97.24;
  HXT_Vmax_1 has mmpmin;
  HXT_Kglc_1 = 1.1918;
  HXT_Kglc_1 has mM;
  HXT_Ki_1 = 0.91;
  HXT_Ki_1 has dimensionless;
  HXK_Vmax_2 = 236.7;
  HXK_Vmax_2 has mMpmin;
  HXK_Kglc_2 = 0.08;
  HXK_Kglc_2 has mM;
  HXK_Katp_2 = 0.15;
  HXK_Katp_2 has mM;
  HXK_Keq_2 = 2000;
  HXK_Keq_2 has dimensionless;
  HXK_Kg6p_2 = 30;
  HXK_Kg6p_2 has mM;
  HXK_Kadp_2 = 0.23;
  HXK_Kadp_2 has mM;
  PGI_Vmax_3 = 1056;
  PGI_Vmax_3 has mMpmin;
  PGI_Kg6p_3 = 1.4;
  PGI_Kg6p_3 has mM;
  PGI_Keq_3 = 0.29;
  PGI_Keq_3 has dimensionless;
  PGI_Kf6p_3 = 0.3;
  PGI_Kf6p_3 has mM;
  PFK_Vmax_4 = 110;
  PFK_Vmax_4 has mMpmin;
  PFK_gR_4 = 5.12;
  PFK_gR_4 has dimensionless;
  PFK_Kf6p_4 = 0.1;
  PFK_Kf6p_4 has mM;
  PFK_Katp_4 = 0.71;
  PFK_Katp_4 has mM;
  PFK_L0_4 = 0.66;
  PFK_L0_4 has dimensionless;
  PFK_Ciatp_4 = 100;
  PFK_Ciatp_4 has dimensionless;
  PFK_Kiatp_4 = 0.65;
  PFK_Kiatp_4 has mM;
  PFK_Camp_4 = 0.0845;
  PFK_Camp_4 has dimensionless;
  PFK_Kamp_4 = 0.0995;
  PFK_Kamp_4 has mM;
  PFK_Cf26_4 = 0.0174;
  PFK_Cf26_4 has dimensionless;
  PFK_Kf26_4 = 0.000682;
  PFK_Kf26_4 has mM;
  PFK_Cf16_4 = 0.397;
  PFK_Cf16_4 has dimensionless;
  PFK_Kf16_4 = 0.111;
  PFK_Kf16_4 has mM;
  PFK_Catp_4 = 3;
  PFK_Catp_4 has dimensionless;
  FBA_Vmax_5 = 94.69;
  FBA_Vmax_5 has mMpmin;
  FBA_Kf16bp_5 = 0.3;
  FBA_Kf16bp_5 has mM;
  FBA_Keq_5 = 0.069;
  FBA_Keq_5 has mM;
  FBA_Kdhap_5 = 2;
  FBA_Kdhap_5 has mM;
  FBA_Kgap_5 = 2.4;
  FBA_Kgap_5 has mM;
  FBA_Kigap_5 = 10;
  FBA_Kigap_5 has mM;
  TPI_k1_6 = 450000;
  TPI_k1_6 has pmin;
  TPI_k2_6 = 10000000;
  TPI_k2_6 has pmin;
  TDH_C_7 = 1;
  TDH_C_7 has dimensionless;
  TDH_Vmaxf_7 = 1152;
  TDH_Vmaxf_7 has mMpmin;
  TDH_Kgap_7 = 0.21;
  TDH_Kgap_7 has mM;
  TDH_Knad_7 = 0.09;
  TDH_Knad_7 has mM;
  TDH_Vmaxr_7 = 6719;
  TDH_Vmaxr_7 has mMpmin;
  TDH_Kbpg_7 = 0.0098;
  TDH_Kbpg_7 has mM;
  TDH_Knadh_7 = 0.06;
  TDH_Knadh_7 has mM;
  PGK_Vmax_8 = 1288;
  PGK_Vmax_8 has mMpmin;
  PGK_Keq_8 = 3200;
  PGK_Keq_8 has dimensionless;
  PGK_Kp3g_8 = 0.53;
  PGK_Kp3g_8 has mM;
  PGK_Katp_8 = 0.3;
  PGK_Katp_8 has mM;
  PGK_Kbpg_8 = 0.003;
  PGK_Kbpg_8 has mM;
  PGK_Kadp_8 = 0.2;
  PGK_Kadp_8 has mM;
  PGM_Vmax_9 = 2585;
  PGM_Vmax_9 has mMpmin;
  PGM_Kp3g_9 = 1.2;
  PGM_Kp3g_9 has mM;
  PGM_Keq_9 = 0.19;
  PGM_Keq_9 has dimensionless;
  PGM_Kp2g_9 = 0.08;
  PGM_Kp2g_9 has mM;
  ENO_Vmax_10 = 201.6;
  ENO_Vmax_10 has mMpmin;
  ENO_Kp2g_10 = 0.04;
  ENO_Kp2g_10 has mM;
  ENO_Keq_10 = 6.7;
  ENO_Keq_10 has dimensionless;
  ENO_Kpep_10 = 0.5;
  ENO_Kpep_10 has mM;
  PYK_Vmax_11 = 1000;
  PYK_Vmax_11 has mMpmin;
  PYK_Kpep_11 = 0.14;
  PYK_Kpep_11 has mM;
  PYK_Kadp_11 = 0.53;
  PYK_Kadp_11 has mM;
  PYK_Keq_11 = 6500;
  PYK_Keq_11 has dimensionless;
  PYK_Kpyr_11 = 21;
  PYK_Kpyr_11 has mM;
  PYK_Katp_11 = 1.5;
  PYK_Katp_11 has mM;
  PDC_Vmax_12 = 857.8;
  PDC_Vmax_12 has mMpmin;
  PDC_Kpyr_12 = 4.33;
  PDC_Kpyr_12 has mM;
  PDC_nH_12 = 1.9;
  PDC_nH_12 has dimensionless;
  ADH_Vmax_13 = 209.5;
  ADH_Vmax_13 has mMpmin;
  ADH_Ketoh_13 = 17;
  ADH_Ketoh_13 has mM;
  ADH_Kinad_13 = 0.92;
  ADH_Kinad_13 has mM;
  ADH_Keq_13 = 6.9e-005;
  ADH_Keq_13 has dimensionless;
  ADH_Knad_13 = 0.17;
  ADH_Knad_13 has mM;
  ADH_Knadh_13 = 0.11;
  ADH_Knadh_13 has mM;
  ADH_Kinadh_13 = 0.031;
  ADH_Kinadh_13 has mM;
  ADH_Kacald_13 = 1.11;
  ADH_Kacald_13 has mM;
  ADH_Kiacald_13 = 1.1;
  ADH_Kiacald_13 has mM;
  ADH_Kietoh_13 = 90;
  ADH_Kietoh_13 has mM;
  ATPase_Katpase_14 = 39.5;
  ATPase_Katpase_14 has pmin;
  AK_k1_15 = 45;
  AK_k1_15 has pmMpmin;
  AK_k2_15 = 100;
  AK_k2_15 has pmMpmin;
  G3PDH_Vmax_16 = 47.11;
  G3PDH_Vmax_16 has mMpmin;
  G3PDH_Kdhap_16 = 0.4;
  G3PDH_Kdhap_16 has mM;
  G3PDH_Knadh_16 = 0.023;
  G3PDH_Knadh_16 has mM;
  G3PDH_Keq_16 = 4300;
  G3PDH_Keq_16 has dimensionless;
  G3PDH_Kglycerol_16 = 1;
  G3PDH_Kglycerol_16 has mM;
  G3PDH_Knad_16 = 0.93;
  G3PDH_Knad_16 has mM;
  glycogen_branch_KGLYCOGEN_17 = 6;
  glycogen_branch_KGLYCOGEN_17 has mMpmin;
  trehalose_synth_Ktrehalose_18 = 2.4;
  trehalose_synth_Ktrehalose_18 has mMpmin;
  succinate_syn_k_19 = 21.4;
  succinate_syn_k_19 has pmin;
  PYR_shut_k1 = 1;
  PYR_shut_k1 has lpmin;
  PYR_shut_k2 = 1;
  PYR_shut_k2 has lpmin;
  PDH_Vmax_PDH = 379.2;
  PDH_Vmax_PDH has mMpmin;
  PDH_NADX_tot = 8.01;
  PDH_NADX_tot has mM;
  PDH_K_PYR = 70;
  PDH_K_PYR has mM;
  PDH_K_NAD = 160;
  PDH_K_NAD has mM;
  PDH_Ki_NADH = 50;
  PDH_Ki_NADH has mM;
  PDH_Ki_PYR = 20;
  PDH_Ki_PYR has mM;

  // Other declarations:
  const cyto, mito, t_m, WGD_E;

  // Unit definitions:
  unit time_unit = 60 second;
  unit substance = 1e-3 mole;
  unit mM = 1e-3 mole / litre;
  unit mmpmin = 1e-3 mole / 60 second;
  unit mMpmin = 1e-3 mole / (60 second * litre);
  unit pmin = 1 / 60 second;
  unit lpmin = litre / 60 second;
  unit pmMpmin = litre / (1e-3 mole * 60 second);

  // Display Names:
  time_unit is "minute";
  substance is "mmol";
  mmpmin is "mmolepermin";
  mMpmin is "mMpermin";
  pmin is "permin";
  lpmin is "permin";
  pmMpmin is "permMpermin";
  cyto is "cytoplasm";
  mito is "mitochondrion";
  F26bP is "Fru2,6-P2";
  t_m is "toggle_transport";
  WGD_E is "WGD enzyme conc change";
  AK is "adenylate_kinase";
  G3PDH is "glycerol-3-phosphate_dehydrogenase";
  glycogen_branch is "glycogenbranch";
  trehalose_synth is "trehalose_synthesis";
  succinate_syn is "succinate_synthesis";
  PYR_shut is "pyruvate_shuttling";
  PDH is "pyruvate dehydrogenase";
end