Location: Metabolic Pathways @ eeb465a0d792 / model / CarbohydrateMetabolism.cellml
- Author:
- Weiwei Ai <wai484@aucklanduni.ac.nz>
- Date:
- 2024-08-21 08:44:37+12:00
- Desc:
- The 1st version from Prof. Hunter.
- Permanent Source URI:
- http://models.cellml.org/workspace/b91/rawfile/eeb465a0d7920d2d188ec1b9ec6f1a5295e0d2aa/model/CarbohydrateMetabolism.cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="my_model" xmlns="http://www.cellml.org/cellml/1.1#" xmlns:cellml="http://www.cellml.org/cellml/1.1#">
<!-- Carbohydrate metabolism-->
<units name="dim">
<unit units="dimensionless"/>
</units>
<units name="per_s">
<unit exponent="-1" units="second"/>
</units>
<units name="C_per_mol">
<unit units="coulomb"/>
<unit exponent="-1" units="mole"/>
</units>
<units name="J_per_mol">
<unit units="joule"/>
<unit exponent="-1" units="mole"/>
</units>
<units name="J_per_C">
<unit units="joule"/>
<unit exponent="-1" units="coulomb"/>
</units>
<component name="main">
<variable initial_value="0" name="t" units="second"/>
<!-- Glycolytic pwy fluxes -->
<!-- Hexokinase: C6H12O6(Glucose) + PO3(ATP-ADP) = C6H11O9P(Glucose-6-Pi) + H -->
<variable name="v_GLYr1" units="per_s"/>
<variable initial_value="10" name="A_GLYf1" units="per_s"/>
<variable initial_value="10" name="A_GLYr1" units="per_s"/>
<variable initial_value="1" name="B_GLYf1" units="dim"/>
<variable initial_value="1" name="B_GLYr1" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr1</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf1</ci>
<ci>q_Glc</ci>
<ci>q_ATP_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr1</ci>
<ci>q_G6P</ci>
<ci>q_ADP_i</ci>
<ci>q_H_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf1</ci>
<ci>q_Glc</ci>
<ci>q_ATP_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr1</ci>
<ci>q_G6P</ci>
<ci>q_ADP_i</ci>
<ci>q_H_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Glucose-6-Pi isomerase: C6H11O9P(Glucose-6-Pi) = C6H11O9P(Fructose-6-Pi)-->
<variable name="v_GLYr2" units="per_s"/>
<variable initial_value="10" name="A_GLYf2" units="per_s"/>
<variable initial_value="1" name="A_GLYr2" units="per_s"/>
<variable initial_value="0" name="A_GLYf2" units="per_s"/>
<variable initial_value="0" name="A_GLYr2" units="per_s"/>
<variable initial_value="1" name="B_GLYf2" units="dim"/>
<variable initial_value="1" name="B_GLYr2" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr2</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf2</ci>
<ci>q_G6P</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr2</ci>
<ci>q_F6P</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf2</ci>
<ci>q_G6P</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr2</ci>
<ci>q_F6P</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Phospho-fructokinase: C6H11O9P(Fructose-6-Pi) + PO3(ATP-ADP) = C6H10O12P2(Fructose-1,6-bisphosphate) + H-->
<variable name="v_GLYr3" units="per_s"/>
<variable initial_value="10" name="A_GLYf3" units="per_s"/>
<variable initial_value="1" name="A_GLYr3" units="per_s"/>
<!-- var A_GLYf3: per_s {init: 0};
var A_GLYr3: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf3" units="dim"/>
<variable initial_value="1" name="B_GLYr3" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr3</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf3</ci>
<ci>q_F6P</ci>
<ci>q_ATP_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr3</ci>
<ci>q_F16bP</ci>
<ci>q_ADP_i</ci>
<ci>q_H_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf3</ci>
<ci>q_F6P</ci>
<ci>q_ATP_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr3</ci>
<ci>q_F16bP</ci>
<ci>q_ADP_i</ci>
<ci>q_H_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Aldolase: C6H10O12P2(Fructose-1,6-bisphosphate) = C3H5O6P(Glyceraldehyde-3-Pi) + C3H5O6P(Dihydroxyacetone-Pi)-->
<variable name="v_GLYr4" units="per_s"/>
<variable initial_value="10" name="A_GLYf4" units="per_s"/>
<variable initial_value="1" name="A_GLYr4" units="per_s"/>
<!-- var A_GLYf4: per_s {init: 0};
var A_GLYr4: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf4" units="dim"/>
<variable initial_value="1" name="B_GLYr4" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr4</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf4</ci>
<ci>q_F16bP</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr4</ci>
<ci>q_GA3P</ci>
<ci>q_DHAP</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf4</ci>
<ci>q_F16bP</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr4</ci>
<ci>q_GA3P</ci>
<ci>q_DHAP</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Triose-Pi-isomerase: C3H5O6P(Glyceraldehyde-3-Pi) = C3H5O6P(Dihydroxyacetone-Pi)-->
<variable name="v_GLYr5" units="per_s"/>
<variable initial_value="10" name="A_GLYf5" units="per_s"/>
<variable initial_value="1" name="A_GLYr5" units="per_s"/>
<!-- var A_GLYf5: per_s {init: 0};
var A_GLYr5: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf5" units="dim"/>
<variable initial_value="1" name="B_GLYr5" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr5</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf5</ci>
<ci>q_GA3P</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr5</ci>
<ci>q_DHAP</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf5</ci>
<ci>q_GA3P</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr5</ci>
<ci>q_DHAP</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Glyceraldehyde-3-Pi-dehydrogenase: C3H5O6P(Glyceraldehyde-3-Pi) + HPO4(Phosphate) = C3H4O10P2(1-3-bisphosphoglycerate) + H + H(NADH-NAD)-->
<variable name="v_GLYr6" units="per_s"/>
<variable initial_value="10" name="A_GLYf6" units="per_s"/>
<variable initial_value="1" name="A_GLYr6" units="per_s"/>
<!-- var A_GLYf6: per_s {init: 0};
var A_GLYr6: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf6" units="dim"/>
<variable initial_value="1" name="B_GLYr6" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr6</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf6</ci>
<ci>q_GA3P</ci>
<ci>q_P_i</ci>
<ci>q_NAD_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr6</ci>
<ci>q_13bPG</ci>
<ci>q_H_i</ci>
<ci>q_NADH_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf6</ci>
<ci>q_GA3P</ci>
<ci>q_P_i</ci>
<ci>q_NAD_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr6</ci>
<ci>q_13bPG</ci>
<ci>q_H_i</ci>
<ci>q_NADH_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Phosphoglycerate kinase: C3H4O10P2(1-3-bisphosphoglycerate) = C3H4O7P(3-phosphoglycerate) + PO3(ATP-ADP)-->
<variable name="v_GLYr7" units="per_s"/>
<variable initial_value="10" name="A_GLYf7" units="per_s"/>
<variable initial_value="1" name="A_GLYr7" units="per_s"/>
<!-- var A_GLYf7: per_s {init: 0};
var A_GLYr7: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf7" units="dim"/>
<variable initial_value="1" name="B_GLYr7" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr7</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf7</ci>
<ci>q_13bPG</ci>
<ci>q_ADP_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr7</ci>
<ci>q_3PG</ci>
<ci>q_ATP_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf7</ci>
<ci>q_13bPG</ci>
<ci>q_ADP_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr7</ci>
<ci>q_3PG</ci>
<ci>q_ATP_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Phosphoglyceromutase: C3H4O7P(3-phosphoglycerate) = C3H4O7P(2-phosphoglycerate)-->
<variable name="v_GLYr8" units="per_s"/>
<variable initial_value="10" name="A_GLYf8" units="per_s"/>
<variable initial_value="1" name="A_GLYr8" units="per_s"/>
<!-- var A_GLYf8: per_s {init: 0};
var A_GLYr8: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf8" units="dim"/>
<variable initial_value="1" name="B_GLYr8" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr8</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf8</ci>
<ci>q_3PG</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr8</ci>
<ci>q_2PG</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf8</ci>
<ci>q_3PG</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr8</ci>
<ci>q_2PG</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Enolase: C3H4O7P(2-phosphoglycerate) = C3H2O6P(Phosphoenolpyruvate) + H2O-->
<variable name="v_GLYr9" units="per_s"/>
<variable initial_value="10" name="A_GLYf9" units="per_s"/>
<variable initial_value="1" name="A_GLYr9" units="per_s"/>
<!-- var A_GLYf9: per_s {init: 0};
var A_GLYr9: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf9" units="dim"/>
<variable initial_value="1" name="B_GLYr9" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr9</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf9</ci>
<ci>q_2PG</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr9</ci>
<ci>q_PEP_i</ci>
<ci>q_H2O_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf9</ci>
<ci>q_2PG</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr9</ci>
<ci>q_PEP_i</ci>
<ci>q_H2O_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Pyruvate kinase: C3H2O6P(Phosphoenolpyruvate) + H = C3H3O3(Pyruvate) + PO3(ATP-ADP) -->
<variable name="v_GLYr10" units="per_s"/>
<variable initial_value="10" name="A_GLYf10" units="per_s"/>
<variable initial_value="1" name="A_GLYr10" units="per_s"/>
<!-- var A_GLYf10: per_s {init: 0};
var A_GLYr10: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf10" units="dim"/>
<variable initial_value="1" name="B_GLYr10" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr10</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf10</ci>
<ci>q_PEP_i</ci>
<ci>q_ADP_i</ci>
<ci>q_H_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr10</ci>
<ci>q_Pyr_i</ci>
<ci>q_ATP_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf10</ci>
<ci>q_PEP_i</ci>
<ci>q_ADP_i</ci>
<ci>q_H_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr10</ci>
<ci>q_Pyr_i</ci>
<ci>q_ATP_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Lactate dehydrogenase: C3H3O3(Pyruvate) + H + H(NADH-NAD) = C3H5O3(Lactate)-->
<variable name="v_GLYr11" units="per_s"/>
<variable initial_value="10" name="A_GLYf11" units="per_s"/>
<variable initial_value="1" name="A_GLYr11" units="per_s"/>
<!-- var A_GLYf11: per_s {init: 0};
var A_GLYr11: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf11" units="dim"/>
<variable initial_value="1" name="B_GLYr11" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr11</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf11</ci>
<ci>q_Pyr_i</ci>
<ci>q_NADH_i</ci>
<ci>q_H_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr11</ci>
<ci>q_Lactate</ci>
<ci>q_NAD_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf11</ci>
<ci>q_Pyr_i</ci>
<ci>q_NADH_i</ci>
<ci>q_H_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr11</ci>
<ci>q_Lactate</ci>
<ci>q_NAD_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Phospho-glucomutase: C6H11O9P(Glucose-1-Pi) = C6H11O9P(Glucose-6-P)-->
<variable name="v_GLYr12" units="per_s"/>
<variable initial_value="10" name="A_GLYf12" units="per_s"/>
<variable initial_value="1" name="A_GLYr12" units="per_s"/>
<!-- var A_GLYf12: per_s {init: 0};
var A_GLYr12: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf12" units="dim"/>
<variable initial_value="1" name="B_GLYr12" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr12</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf12</ci>
<ci>q_G1P</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr12</ci>
<ci>q_G6P</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf12</ci>
<ci>q_G1P</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr12</ci>
<ci>q_G6P</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Glucose-1-Pi uridyl transferase: C6H11O9P(Glucose-1-Pi) + C9H15N2O15P3(UTP) = C15H24N2O16P2(UDP-glucose) + H2P2O8(Diphosphate) (Oxygen imbalance)-->
<variable name="v_GLYr13" units="per_s"/>
<variable initial_value="10" name="A_GLYf13" units="per_s"/>
<variable initial_value="1" name="A_GLYr13" units="per_s"/>
<!-- var A_GLYf13: per_s {init: 0};
var A_GLYr13: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf13" units="dim"/>
<variable initial_value="1" name="B_GLYr13" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr13</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf13</ci>
<ci>q_G1P</ci>
<ci>q_UTP</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr13</ci>
<ci>q_PPi</ci>
<ci>q_UDPglc</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf13</ci>
<ci>q_G1P</ci>
<ci>q_UTP</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr13</ci>
<ci>q_PPi</ci>
<ci>q_UDPglc</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Glycogen synthetase: C15H24N2O17P2(UDP-glucose) = C9H14N2O12P2(UDP) + C6H10O5(add to Glycogen)-->
<variable name="v_GLYr14" units="per_s"/>
<variable initial_value="10" name="A_GLYf14" units="per_s"/>
<variable initial_value="1" name="A_GLYr14" units="per_s"/>
<!-- var A_GLYf14: per_s {init: 0};
var A_GLYr14: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf14" units="dim"/>
<variable initial_value="1" name="B_GLYr14" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr14</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf14</ci>
<ci>q_UDPglc</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr14</ci>
<ci>q_Gly</ci>
<ci>q_UDP</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf14</ci>
<ci>q_UDPglc</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr14</ci>
<ci>q_Gly</ci>
<ci>q_UDP</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Glycogen phosphorylase: C6H10O5(Glycogen) + HPO4(Phosphate) = C6H11O9P(Glucose-1-Pi)-->
<variable name="v_GLYr15" units="per_s"/>
<variable initial_value="10" name="A_GLYf15" units="per_s"/>
<variable initial_value="1" name="A_GLYr15" units="per_s"/>
<!-- var A_GLYf15: per_s {init: 0};
var A_GLYr15: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf15" units="dim"/>
<variable initial_value="1" name="B_GLYr15" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr15</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf15</ci>
<ci>q_Gly</ci>
<ci>q_P_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr15</ci>
<ci>q_G1P</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf15</ci>
<ci>q_Gly</ci>
<ci>q_P_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr15</ci>
<ci>q_G1P</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Pyrophosphatase: H2P2O8(Diphosphate) = H2P2O8(2xPhosphate)-->
<variable name="v_GLYr16" units="per_s"/>
<variable initial_value="10" name="A_GLYf16" units="per_s"/>
<variable initial_value="1" name="A_GLYr16" units="per_s"/>
<!-- var A_GLYf16: per_s {init: 0};
var A_GLYr16: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf16" units="dim"/>
<variable initial_value="1" name="B_GLYr16" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr16</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf16</ci>
<ci>q_PPi</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr16</ci>
<apply>
<power/>
<ci>q_P_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf16</ci>
<ci>q_PPi</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr16</ci>
<apply>
<power/>
<ci>q_P_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Fructose-1,6-biphosphatase: C6H10O12P2(Fructose-1,6-bisphosphate) + H2O = C6H11O9P(Fructose-6-Pi) + HPO4(Phosphate)-->
<variable name="v_GLYr17" units="per_s"/>
<variable initial_value="1" name="A_GLYf17" units="per_s"/>
<variable initial_value="0.1" name="A_GLYr17" units="per_s"/>
<!-- var A_GLYf17: per_s {init: 0};
var A_GLYr17: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf17" units="dim"/>
<variable initial_value="1" name="B_GLYr17" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr17</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf17</ci>
<ci>q_F16bP</ci>
<ci>q_H2O_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr17</ci>
<ci>q_F6P</ci>
<ci>q_P_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf17</ci>
<ci>q_F16bP</ci>
<ci>q_H2O_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr17</ci>
<ci>q_F6P</ci>
<ci>q_P_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Glucose-6-P-ase: C6H11O9P(Glucose-6-Pi) + H2O = C6H12O6(Glucose) + HPO4(Phosphate)-->
<variable name="v_GLYr18" units="per_s"/>
<variable initial_value="1" name="A_GLYf18" units="per_s"/>
<variable initial_value="0.1" name="A_GLYr18" units="per_s"/>
<!-- var A_GLYf18: per_s {init: 0};
var A_GLYr18: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf18" units="dim"/>
<variable initial_value="1" name="B_GLYr18" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr18</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf18</ci>
<ci>q_G6P</ci>
<ci>q_H2O_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr18</ci>
<ci>q_Glc</ci>
<ci>q_P_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf18</ci>
<ci>q_G6P</ci>
<ci>q_H2O_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr18</ci>
<ci>q_Glc</ci>
<ci>q_P_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Glycerol-3-Pi dehydrogenase: C3H5O6P(Dihydroxyacetone-Pi) + H2 = C3H7O6P(Glycerol-3-Pi)-->
<variable name="v_GLYr19" units="per_s"/>
<variable initial_value="10" name="A_GLYf19" units="per_s"/>
<variable initial_value="1" name="A_GLYr19" units="per_s"/>
<!-- var A_GLYf19: per_s {init: 0};
var A_GLYr19: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf19" units="dim"/>
<variable initial_value="1" name="B_GLYr19" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr19</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf19</ci>
<ci>q_DHAP</ci>
<ci>q_NADH_i</ci>
<ci>q_H_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr19</ci>
<ci>q_Gly3P</ci>
<ci>q_NAD_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf19</ci>
<ci>q_DHAP</ci>
<ci>q_NADH_i</ci>
<ci>q_H_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr19</ci>
<ci>q_Gly3P</ci>
<ci>q_NAD_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Glycerol kinase: C3H7O6P(Glycerol-3-Pi) + H = C3H8O3(Glycerol) + PO3(ATP-ADP)-->
<variable name="v_GLYr20" units="per_s"/>
<variable initial_value="10" name="A_GLYf20" units="per_s"/>
<variable initial_value="1" name="A_GLYr20" units="per_s"/>
<!-- var A_GLYf20: per_s {init: 0};
var A_GLYr20: per_s {init: 0};-->
<variable initial_value="1" name="B_GLYf20" units="dim"/>
<variable initial_value="1" name="B_GLYr20" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_GLYr20</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_GLYf20</ci>
<ci>q_Gly3P</ci>
<ci>q_ADP_i</ci>
<ci>q_H_i</ci>
</apply>
<apply>
<times/>
<ci>A_GLYr20</ci>
<ci>q_Glycerol</ci>
<ci>q_ATP_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_GLYf20</ci>
<ci>q_Gly3P</ci>
<ci>q_ADP_i</ci>
<ci>q_H_i</ci>
</apply>
<apply>
<times/>
<ci>B_GLYr20</ci>
<ci>q_Glycerol</ci>
<ci>q_ATP_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Pyruvate transport (VDAC1): C3H3O3(Pyruvate) = C3H3O3(Pyruvate)-->
<variable name="v_PYRr1" units="per_s"/>
<variable initial_value="10" name="A_PYRf1" units="per_s"/>
<!-- var A_PYRf1: per_s {init: 0};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_PYRr1</ci>
<apply>
<times/>
<ci>A_PYRf1</ci>
<apply>
<minus/>
<ci>q_Pyr_i</ci>
<ci>q_Pyr_m</ci>
</apply>
</apply>
</apply>
</math>
<!-- PEP transport (SLC25A1): C3H2O6P(Phosphoenolpyruvate) = C3H2O6P(Phosphoenolpyruvate)-->
<variable name="v_PYRr2" units="per_s"/>
<variable initial_value="10" name="A_PYRf2" units="per_s"/>
<variable initial_value="1" name="A_PYRr2" units="per_s"/>
<!-- var A_PYRf2: per_s {init: 0};
var A_PYRr2: per_s {init: 0};-->
<variable initial_value="1" name="B_PYRf2" units="dim"/>
<variable initial_value="1" name="B_PYRr2" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_PYRr2</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_PYRf2</ci>
<ci>q_PEP_m</ci>
</apply>
<apply>
<times/>
<ci>A_PYRr2</ci>
<ci>q_PEP_i</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_PYRf2</ci>
<ci>q_PEP_m</ci>
</apply>
<apply>
<times/>
<ci>B_PYRr2</ci>
<ci>q_PEP_i</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- TCA cycle fluxes -->
<!-- Citrate synthase: C4H2O5(Oxaloacetate) + C2H3OS(AcetlyCoA) + H2O = C6H5O7(citrate) + SH(CoA-SH) + H-->
<variable name="v_TCAr1" units="per_s"/>
<variable initial_value="10" name="A_TCAf1" units="per_s"/>
<variable initial_value="1" name="A_TCAr1" units="per_s"/>
<!-- var A_TCAf1: per_s {init: 0};
var A_TCAr1: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf1" units="dim"/>
<variable initial_value="1" name="B_TCAr1" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr1</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf1</ci>
<ci>q_AcetylCoA</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_H2O_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr1</ci>
<ci>q_CoASH</ci>
<ci>q_Citrate</ci>
<ci>q_H_m</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf1</ci>
<ci>q_AcetylCoA</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_H2O_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr1</ci>
<ci>q_CoASH</ci>
<ci>q_Citrate</ci>
<ci>q_H_m</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Aconitase: C6H5O7(citrate) = C6H5O7(isocitrate)-->
<variable name="v_TCAr2" units="per_s"/>
<variable initial_value="10" name="A_TCAf2" units="per_s"/>
<variable initial_value="1" name="A_TCAr2" units="per_s"/>
<!-- var A_TCAf2: per_s {init: 0};
var A_TCAr2: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf2" units="dim"/>
<variable initial_value="1" name="B_TCAr2" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr2</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf2</ci>
<ci>q_Citrate</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr2</ci>
<ci>q_Isocitrate</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf2</ci>
<ci>q_Citrate</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr2</ci>
<ci>q_Isocitrate</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Isocitrate dehydrogenase: C6H5O7(isocitrate) = C5H4O5(2-Oxogluterate) + H(NADH-NAD) + CO2-->
<variable name="v_TCAr3" units="per_s"/>
<variable initial_value="10" name="A_TCAf3" units="per_s"/>
<variable initial_value="1" name="A_TCAr3" units="per_s"/>
<!-- var A_TCAf3: per_s {init: 0};
var A_TCAr3: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf3" units="dim"/>
<variable initial_value="1" name="B_TCAr3" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr3</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf3</ci>
<ci>q_Isocitrate</ci>
<ci>q_NAD_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr3</ci>
<ci>q_2OG</ci>
<ci>q_NADH_m</ci>
<ci>q_CO2</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf3</ci>
<ci>q_Isocitrate</ci>
<ci>q_NAD_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr3</ci>
<ci>q_2OG</ci>
<ci>q_NADH_m</ci>
<ci>q_CO2</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- 2-oxogluterate dehydrogenase: C5H4O5(2-Oxogluterate) + SH(CoA-SH) = C4H4O3S(Succinyl-CoA) + H(NADH-NAD) + CO2-->
<variable name="v_TCAr4" units="per_s"/>
<variable initial_value="10" name="A_TCAf4" units="per_s"/>
<variable initial_value="1" name="A_TCAr4" units="per_s"/>
<!-- var A_TCAf4: per_s {init: 0};
var A_TCAr4: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf4" units="dim"/>
<variable initial_value="1" name="B_TCAr4" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr4</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf4</ci>
<ci>q_2OG</ci>
<ci>q_CoASH</ci>
<ci>q_NAD_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr4</ci>
<ci>q_SuccinylCoA</ci>
<ci>q_NADH_m</ci>
<ci>q_CO2</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf4</ci>
<ci>q_2OG</ci>
<ci>q_CoASH</ci>
<ci>q_NAD_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr4</ci>
<ci>q_SuccinylCoA</ci>
<ci>q_NADH_m</ci>
<ci>q_CO2</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Succinyl-CoA synthase: C4H4O3S(Succinyl-CoA) + HPO4(Phosphate) = C4H4O4(Succinate) + SH(CoA-SH) + PO3(GTP-GDP)-->
<variable name="v_TCAr5" units="per_s"/>
<variable initial_value="10" name="A_TCAf5" units="per_s"/>
<variable initial_value="1" name="A_TCAr5" units="per_s"/>
<!-- var A_TCAf5: per_s {init: 0};
var A_TCAr5: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf5" units="dim"/>
<variable initial_value="1" name="B_TCAr5" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr5</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf5</ci>
<ci>q_SuccinylCoA</ci>
<ci>q_GDP</ci>
<ci>q_P_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr5</ci>
<ci>q_Succinate</ci>
<ci>q_CoASH</ci>
<ci>q_GTP</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf5</ci>
<ci>q_SuccinylCoA</ci>
<ci>q_GDP</ci>
<ci>q_P_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr5</ci>
<ci>q_Succinate</ci>
<ci>q_CoASH</ci>
<ci>q_GTP</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Succinate dehydrogenase: C4H4O4(Succinate) = C4H2O4(Fumarate) + H2(FADH2-FAD)-->
<variable name="v_TCAr6" units="per_s"/>
<variable initial_value="10" name="A_TCAf6" units="per_s"/>
<variable initial_value="1" name="A_TCAr6" units="per_s"/>
<!-- var A_TCAf6: per_s {init: 0};
var A_TCAr6: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf6" units="dim"/>
<variable initial_value="1" name="B_TCAr6" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr6</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf6</ci>
<ci>q_Succinate</ci>
<ci>q_FAD</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr6</ci>
<ci>q_Fumarate</ci>
<ci>q_FADH2</ci>
<ci>q_e</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf6</ci>
<ci>q_Succinate</ci>
<ci>q_FAD</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr6</ci>
<ci>q_Fumarate</ci>
<ci>q_FADH2</ci>
<ci>q_e</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Fumarate hydratase: C4H2O4(Fumarate) + H2O = C4H4O5(Malate)-->
<variable name="v_TCAr7" units="per_s"/>
<variable initial_value="10" name="A_TCAf7" units="per_s"/>
<variable initial_value="1" name="A_TCAr7" units="per_s"/>
<!-- var A_TCAf7: per_s {init: 0};
var A_TCAr7: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf7" units="dim"/>
<variable initial_value="1" name="B_TCAr7" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr7</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf7</ci>
<ci>q_Fumarate</ci>
<ci>q_H2O_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr7</ci>
<ci>q_Malate</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf7</ci>
<ci>q_Fumarate</ci>
<ci>q_H2O_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr7</ci>
<ci>q_Malate</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Malate dehydrogenase: C4H4O5(Malate) = C4H2O5(Oxoloacetate) + H(NADH-NAD) + H-->
<variable name="v_TCAr8" units="per_s"/>
<variable initial_value="10" name="A_TCAf8" units="per_s"/>
<variable initial_value="1" name="A_TCAr8" units="per_s"/>
<!-- var A_TCAf8: per_s {init: 0};
var A_TCAr8: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf8" units="dim"/>
<variable initial_value="1" name="B_TCAr8" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr8</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf8</ci>
<ci>q_Malate</ci>
<ci>q_NAD_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr8</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_NADH_m</ci>
<ci>q_H_m</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf8</ci>
<ci>q_Malate</ci>
<ci>q_NAD_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr8</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_NADH_m</ci>
<ci>q_H_m</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Pyruvate kinase: C3H2O6P(PEP) + H = C3H3O3(Pyruvate) + PO3(ATP-ADP)-->
<variable name="v_TCAr9" units="per_s"/>
<variable initial_value="10" name="A_TCAf9" units="per_s"/>
<variable initial_value="1" name="A_TCAr9" units="per_s"/>
<!-- var A_TCAf9: per_s {init: 0};
var A_TCAr9: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf9" units="dim"/>
<variable initial_value="1" name="B_TCAr9" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr9</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf9</ci>
<ci>q_PEP_m</ci>
<ci>q_ADP_m</ci>
<ci>q_H_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr9</ci>
<ci>q_Pyr_m</ci>
<ci>q_ATP_m</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf9</ci>
<ci>q_PEP_m</ci>
<ci>q_ADP_m</ci>
<ci>q_H_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr9</ci>
<ci>q_Pyr_m</ci>
<ci>q_ATP_m</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Pyruvate dehydrogenase: C3H3O3(Pyruvate) + SH(CoA-SH) + H = C2H3OS(Acetyl-CoA) + H(NADH-NAD) + H + CO2 -->
<variable name="v_TCAr10" units="per_s"/>
<variable initial_value="10" name="A_TCAf10" units="per_s"/>
<variable initial_value="1" name="A_TCAr10" units="per_s"/>
<!-- var A_TCAf10: per_s {init: 0};
var A_TCAr10: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf10" units="dim"/>
<variable initial_value="1" name="B_TCAr10" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr10</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf10</ci>
<ci>q_Pyr_m</ci>
<ci>q_CoASH</ci>
<ci>q_NAD_m</ci>
<ci>q_H_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr10</ci>
<ci>q_AcetylCoA</ci>
<ci>q_NADH_m</ci>
<ci>q_H_m</ci>
<ci>q_CO2</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf10</ci>
<ci>q_Pyr_m</ci>
<ci>q_CoASH</ci>
<ci>q_NAD_m</ci>
<ci>q_H_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr10</ci>
<ci>q_AcetylCoA</ci>
<ci>q_NADH_m</ci>
<ci>q_H_m</ci>
<ci>q_CO2</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Pyruvate carboxylase: C3H3O3(Pyruvate) + CO2 + PO3(ATP-ADP) = C4H2O5(Oxoloacetate) + HPO4
Pyruvate carboxylase: C3H3O3(Pyruvate) + HCO3 + PO3(ATP-ADP) = C4H2O5(Oxoloacetate) + HPO4 + H-->
<variable name="v_TCAr11" units="per_s"/>
<variable initial_value="10" name="A_TCAf11" units="per_s"/>
<variable initial_value="1" name="A_TCAr11" units="per_s"/>
<!-- var A_TCAf11: per_s {init: 0};
var A_TCAr11: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf11" units="dim"/>
<variable initial_value="1" name="B_TCAr11" units="dim"/>
<!-- v_TCAr11 = (A_TCAf11*q_Pyr_m*q_CO2*q_ATP_m-A_TCAr11*q_Oxaloacetate*q_ADP_m*q_P_m)/(1{dim}+B_TCAf11*q_Pyr_m*q_CO2*q_ATP_m+B_TCAr11*q_Oxaloacetate*q_ADP_m*q_P_m);-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr11</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf11</ci>
<ci>q_Pyr_m</ci>
<ci>q_HCO3_m</ci>
<ci>q_ATP_m</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr11</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_ADP_m</ci>
<ci>q_P_m</ci>
<ci>q_H_m</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf11</ci>
<ci>q_Pyr_m</ci>
<ci>q_HCO3_m</ci>
<ci>q_ATP_m</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr11</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_ADP_m</ci>
<ci>q_P_m</ci>
<ci>q_H_m</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Phosphoenolpyruvate carboxykinase: C4H2O5(Oxoloacetate) + PO3(GTP-GDP) = C3H2O6P(PEP) + CO2-->
<variable name="v_TCAr12" units="per_s"/>
<variable initial_value="10" name="A_TCAf12" units="per_s"/>
<variable initial_value="1" name="A_TCAr12" units="per_s"/>
<!-- var A_TCAf12: per_s {init: 0};
var A_TCAr12: per_s {init: 0};-->
<variable initial_value="1" name="B_TCAf12" units="dim"/>
<variable initial_value="1" name="B_TCAr12" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_TCAr12</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_TCAf12</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_GTP</ci>
</apply>
<apply>
<times/>
<ci>A_TCAr12</ci>
<ci>q_PEP_m</ci>
<ci>q_GDP</ci>
<ci>q_CO2</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_TCAf12</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_GTP</ci>
</apply>
<apply>
<times/>
<ci>B_TCAr12</ci>
<ci>q_PEP_m</ci>
<ci>q_GDP</ci>
<ci>q_CO2</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Electron transport chain fluxes -->
<!-- NADH-ubiquinone oxido-reductase (CI): H5 = H(NADH-NAD) + H4-->
<variable name="v_ETCr1" units="per_s"/>
<!-- var A_ETCf1: per_s {init: 10};
var A_ETCr1: per_s {init: 1};-->
<variable initial_value="0" name="A_ETCf1" units="per_s"/>
<variable initial_value="0" name="A_ETCr1" units="per_s"/>
<variable initial_value="1" name="B_ETCf1" units="dim"/>
<variable initial_value="1" name="B_ETCr1" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ETCr1</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_ETCf1</ci>
<ci>q_NADH_m</ci>
<apply>
<power/>
<ci>q_H_m</ci>
<cn cellml:units="dim">5</cn>
</apply>
</apply>
<apply>
<times/>
<ci>A_ETCr1</ci>
<ci>q_NAD_m</ci>
<ci>q_e</ci>
<apply>
<power/>
<ci>q_H_p</ci>
<cn cellml:units="dim">4</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_ETCf1</ci>
<ci>q_NADH_m</ci>
<apply>
<power/>
<ci>q_H_m</ci>
<cn cellml:units="dim">5</cn>
</apply>
</apply>
<apply>
<times/>
<ci>B_ETCr1</ci>
<ci>q_NAD_m</ci>
<ci>q_e</ci>
<apply>
<power/>
<ci>q_H_p</ci>
<cn cellml:units="dim">4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Ubiquinol cytochrome c oxidoreductase (CIII): H4 = H4 -->
<variable name="v_ETCr2" units="per_s"/>
<!-- var A_ETCf2: per_s {init: 10};
var A_ETCr2: per_s {init: 1};-->
<variable initial_value="0" name="A_ETCf2" units="per_s"/>
<variable initial_value="0" name="A_ETCr2" units="per_s"/>
<variable initial_value="1" name="B_ETCf2" units="dim"/>
<variable initial_value="1" name="B_ETCr2" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ETCr2</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_ETCf2</ci>
<apply>
<power/>
<ci>q_H_m</ci>
<cn cellml:units="dim">4</cn>
</apply>
</apply>
<apply>
<times/>
<ci>A_ETCr2</ci>
<apply>
<power/>
<ci>q_H_p</ci>
<cn cellml:units="dim">4</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_ETCf2</ci>
<apply>
<power/>
<ci>q_H_m</ci>
<cn cellml:units="dim">4</cn>
</apply>
</apply>
<apply>
<times/>
<ci>B_ETCr2</ci>
<apply>
<power/>
<ci>q_H_p</ci>
<cn cellml:units="dim">4</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Cytochrome c oxidase (CIV): O2 + H6 = H4O2(2xH2O) + H2-->
<variable name="v_ETCr3" units="per_s"/>
<!-- var A_ETCf3: per_s {init: 10};
var A_ETCr3: per_s {init: 1};-->
<variable initial_value="0" name="A_ETCf3" units="per_s"/>
<variable initial_value="0" name="A_ETCr3" units="per_s"/>
<variable initial_value="1" name="B_ETCf3" units="dim"/>
<variable initial_value="1" name="B_ETCr3" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ETCr3</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_ETCf3</ci>
<ci>q_Oxy</ci>
<apply>
<power/>
<ci>q_H_m</ci>
<cn cellml:units="dim">6</cn>
</apply>
<apply>
<power/>
<ci>q_e</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<ci>A_ETCr3</ci>
<apply>
<power/>
<ci>q_H2O_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>q_H_p</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_ETCf3</ci>
<ci>q_Oxy</ci>
<apply>
<power/>
<ci>q_H_m</ci>
<cn cellml:units="dim">6</cn>
</apply>
<apply>
<power/>
<ci>q_e</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<ci>B_ETCr3</ci>
<apply>
<power/>
<ci>q_H2O_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<apply>
<power/>
<ci>q_H_p</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Inorganic phosphate carrier (PiC) (SLC25A3): HPO4 + H = HPO4 + H -->
<variable name="v_ETCr4" units="per_s"/>
<!-- var A_ETCf4: per_s {init: 10};
var A_ETCr4: per_s {init: 1};-->
<variable initial_value="0" name="A_ETCf4" units="per_s"/>
<variable initial_value="0" name="A_ETCr4" units="per_s"/>
<variable initial_value="1" name="B_ETCf4" units="dim"/>
<variable initial_value="1" name="B_ETCr4" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ETCr4</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_ETCf4</ci>
<ci>q_P_p</ci>
<ci>q_H_p</ci>
</apply>
<apply>
<times/>
<ci>A_ETCr4</ci>
<ci>q_P_m</ci>
<ci>q_H_m</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_ETCf4</ci>
<ci>q_P_p</ci>
<ci>q_H_p</ci>
</apply>
<apply>
<times/>
<ci>B_ETCr4</ci>
<ci>q_P_m</ci>
<ci>q_H_m</ci>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Adenosine nucleotide translocase (ANT1) (SLC25A4): ADP + ATP = ATP + ADP -->
<variable name="v_ETCr5" units="per_s"/>
<!-- var V_ETCr5max: per_s {init: 10};-->
<variable initial_value="0" name="V_ETCr5max" units="per_s"/>
<variable initial_value="2" name="kappa_f" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ETCr5</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>V_ETCr5max</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>q_ATP_m</ci>
<ci>q_ADP_p</ci>
</apply>
<apply>
<times/>
<ci>q_ATP_p</ci>
<ci>q_ADP_m</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<ci>q_ATP_m</ci>
<apply>
<times/>
<ci>kappa_f</ci>
<ci>q_ADP_m</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<ci>q_ATP_p</ci>
<ci>q_ADP_p</ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<ci>q_ATP_p</ci>
<apply>
<times/>
<ci>kappa_f</ci>
<ci>q_ADP_p</ci>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<ci>q_ATP_m</ci>
<ci>q_ADP_m</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- ATP-synthase: HPO4 + H9 = PO3(ATP-ADP) + H2O + H8-->
<variable name="v_ETCr6" units="per_s"/>
<!-- var A_ETCf6: per_s {init: 10};
var A_ETCr6: per_s {init: 1};-->
<variable initial_value="0" name="A_ETCf6" units="per_s"/>
<variable initial_value="0" name="A_ETCr6" units="per_s"/>
<variable initial_value="1" name="B_ETCf6" units="dim"/>
<variable initial_value="1" name="B_ETCr6" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ETCr6</ci>
<apply>
<divide/>
<apply>
<minus/>
<apply>
<times/>
<ci>A_ETCf6</ci>
<ci>q_ADP_m</ci>
<ci>q_P_m</ci>
<apply>
<power/>
<ci>q_H_p</ci>
<cn cellml:units="dim">9</cn>
</apply>
</apply>
<apply>
<times/>
<ci>A_ETCr6</ci>
<ci>q_ATP_m</ci>
<ci>q_H2O_m</ci>
<apply>
<power/>
<ci>q_H_m</ci>
<cn cellml:units="dim">8</cn>
</apply>
</apply>
</apply>
<apply>
<plus/>
<cn cellml:units="dim">1</cn>
<apply>
<times/>
<ci>B_ETCf6</ci>
<ci>q_ADP_m</ci>
<ci>q_P_m</ci>
<apply>
<power/>
<ci>q_H_p</ci>
<cn cellml:units="dim">9</cn>
</apply>
</apply>
<apply>
<times/>
<ci>B_ETCr6</ci>
<ci>q_ATP_m</ci>
<ci>q_H2O_m</ci>
<apply>
<power/>
<ci>q_H_m</ci>
<cn cellml:units="dim">8</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
</math>
<!-- Cytosol substrates-->
<variable initial_value="10" name="q_UTP" units="dim"/>
<variable initial_value="0" name="q_UDP" units="dim"/>
<variable initial_value="10" name="q_ATP_i" units="dim"/>
<variable initial_value="0" name="q_ADP_i" units="dim"/>
<variable initial_value="0" name="q_P_i" units="dim"/>
<variable initial_value="10" name="q_H2O_i" units="dim"/>
<variable initial_value="0" name="q_H_i" units="dim"/>
<variable initial_value="20" name="q_Glc" units="dim"/>
<variable initial_value="0" name="q_G6P" units="dim"/>
<variable initial_value="0" name="q_G1P" units="dim"/>
<variable initial_value="0" name="q_Gly" units="dim"/>
<variable initial_value="0" name="q_UDPglc" units="dim"/>
<variable initial_value="0" name="q_PPi" units="dim"/>
<variable initial_value="0" name="q_F6P" units="dim"/>
<variable initial_value="0" name="q_F16bP" units="dim"/>
<variable initial_value="0" name="q_GA3P" units="dim"/>
<variable initial_value="0" name="q_13bPG" units="dim"/>
<variable initial_value="0" name="q_3PG" units="dim"/>
<variable initial_value="0" name="q_2PG" units="dim"/>
<variable initial_value="0" name="q_PEP_i" units="dim"/>
<variable initial_value="0" name="q_Lactate" units="dim"/>
<variable initial_value="0" name="q_DHAP" units="dim"/>
<variable initial_value="0" name="q_Gly3P" units="dim"/>
<variable initial_value="0" name="q_Glycerol" units="dim"/>
<variable initial_value="0" name="q_Pyr_i" units="dim"/>
<variable initial_value="0.5" name="q_NADH_i" units="dim"/>
<variable initial_value="0.5" name="q_NAD_i" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_UTP</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr13</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_UDP</ci>
</apply>
<ci>v_GLYr14</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ATP_i</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_GLYr1</ci>
</apply>
<ci>v_GLYr3</ci>
</apply>
<ci>v_GLYr7</ci>
<ci>v_GLYr10</ci>
<ci>v_GLYr20</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ADP_i</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_GLYr1</ci>
<ci>v_GLYr3</ci>
</apply>
<ci>v_GLYr7</ci>
</apply>
<ci>v_GLYr10</ci>
</apply>
<ci>v_GLYr20</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_P_i</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_GLYr6</ci>
</apply>
<ci>v_GLYr15</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<ci>v_GLYr16</ci>
</apply>
<ci>v_GLYr17</ci>
<ci>v_GLYr18</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_H2O_i</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_GLYr9</ci>
<ci>v_GLYr17</ci>
</apply>
<ci>v_GLYr18</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_H_i</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_GLYr1</ci>
<ci>v_GLYr3</ci>
<ci>v_GLYr6</ci>
</apply>
<ci>v_GLYr10</ci>
</apply>
<ci>v_GLYr11</ci>
</apply>
<ci>v_GLYr19</ci>
</apply>
<ci>v_GLYr20</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Glc</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr18</ci>
<ci>v_GLYr1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_G6P</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_GLYr1</ci>
<ci>v_GLYr2</ci>
</apply>
<ci>v_GLYr12</ci>
</apply>
<ci>v_GLYr18</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_G1P</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_GLYr12</ci>
</apply>
<ci>v_GLYr13</ci>
</apply>
<ci>v_GLYr15</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Gly</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr14</ci>
<ci>v_GLYr15</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_UDPglc</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr13</ci>
<ci>v_GLYr14</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_PPi</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr13</ci>
<ci>v_GLYr16</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_F6P</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_GLYr2</ci>
<ci>v_GLYr3</ci>
</apply>
<ci>v_GLYr17</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_F16bP</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_GLYr3</ci>
<ci>v_GLYr4</ci>
</apply>
<ci>v_GLYr17</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_GA3P</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_GLYr4</ci>
<ci>v_GLYr5</ci>
</apply>
<ci>v_GLYr6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_13bPG</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr6</ci>
<ci>v_GLYr7</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_3PG</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr7</ci>
<ci>v_GLYr8</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_2PG</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr8</ci>
<ci>v_GLYr9</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_PEP_i</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_GLYr9</ci>
<ci>v_GLYr10</ci>
</apply>
<ci>v_PYRr2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Lactate</ci>
</apply>
<ci>v_GLYr11</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_DHAP</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_GLYr4</ci>
<ci>v_GLYr5</ci>
</apply>
<ci>v_GLYr19</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Gly3P</ci>
</apply>
<apply>
<minus/>
<ci>v_GLYr19</ci>
<ci>v_GLYr20</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Glycerol</ci>
</apply>
<ci>v_GLYr20</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Pyr_i</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_GLYr10</ci>
<ci>v_GLYr11</ci>
</apply>
<ci>v_PYRr1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_NADH_i</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_GLYr6</ci>
<ci>v_GLYr11</ci>
</apply>
<ci>v_GLYr19</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_NAD_i</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_GLYr6</ci>
</apply>
<ci>v_GLYr11</ci>
<ci>v_GLYr19</ci>
</apply>
</apply>
</math>
<!-- Matrix substrates-->
<variable initial_value="0" name="q_Pyr_m" units="dim"/>
<variable initial_value="10" name="q_AcetylCoA" units="dim"/>
<variable initial_value="10" name="q_CoASH" units="dim"/>
<variable initial_value="0" name="q_Citrate" units="dim"/>
<variable initial_value="0" name="q_Isocitrate" units="dim"/>
<variable initial_value="0" name="q_2OG" units="dim"/>
<variable initial_value="0" name="q_SuccinylCoA" units="dim"/>
<variable initial_value="0" name="q_Succinate" units="dim"/>
<variable initial_value="10" name="q_FAD" units="dim"/>
<variable initial_value="0" name="q_FADH2" units="dim"/>
<variable initial_value="0" name="q_Fumarate" units="dim"/>
<variable initial_value="0" name="q_Malate" units="dim"/>
<variable initial_value="1" name="q_Oxaloacetate" units="dim"/>
<variable initial_value="0" name="q_PEP_m" units="dim"/>
<variable initial_value="0" name="q_CO2" units="dim"/>
<variable initial_value="10" name="q_HCO3_m" units="dim"/>
<variable initial_value="10" name="q_GTP" units="dim"/>
<variable initial_value="0" name="q_GDP" units="dim"/>
<variable initial_value="0" name="q_P_m" units="dim"/>
<variable initial_value="0.5" name="q_H_m" units="dim"/>
<variable initial_value="10" name="q_ATP_m" units="dim"/>
<variable initial_value="0" name="q_ADP_m" units="dim"/>
<variable initial_value="10" name="q_H2O_m" units="dim"/>
<variable initial_value="5" name="q_NADH_m" units="dim"/>
<variable initial_value="5" name="q_NAD_m" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Pyr_m</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_PYRr1</ci>
<ci>v_TCAr9</ci>
</apply>
<ci>v_TCAr10</ci>
</apply>
<ci>v_TCAr11</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_AcetylCoA</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_TCAr1</ci>
</apply>
<ci>v_TCAr10</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_CoASH</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_TCAr1</ci>
<ci>v_TCAr4</ci>
</apply>
<ci>v_TCAr5</ci>
</apply>
<ci>v_TCAr10</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Citrate</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr1</ci>
<ci>v_TCAr2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Isocitrate</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr2</ci>
<ci>v_TCAr3</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_2OG</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr3</ci>
<ci>v_TCAr4</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_SuccinylCoA</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr4</ci>
<ci>v_TCAr5</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Succinate</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr5</ci>
<ci>v_TCAr6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_FAD</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_FADH2</ci>
</apply>
<ci>v_TCAr6</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Fumarate</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr6</ci>
<ci>v_TCAr7</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Malate</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr7</ci>
<ci>v_TCAr8</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Oxaloacetate</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_TCAr1</ci>
</apply>
<ci>v_TCAr8</ci>
<ci>v_TCAr11</ci>
</apply>
<ci>v_TCAr12</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_PEP_m</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_TCAr9</ci>
</apply>
<ci>v_TCAr12</ci>
</apply>
<ci>v_PYRr2</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_CO2</ci>
</apply>
<apply>
<plus/>
<ci>v_TCAr3</ci>
<ci>v_TCAr4</ci>
<ci>v_TCAr10</ci>
<ci>v_TCAr12</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_HCO3_m</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr11</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_GTP</ci>
</apply>
<apply>
<minus/>
<ci>v_TCAr5</ci>
<ci>v_TCAr12</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_GDP</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_TCAr5</ci>
</apply>
<ci>v_TCAr12</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_P_m</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_TCAr5</ci>
</apply>
<ci>v_TCAr11</ci>
<ci>v_ETCr4</ci>
</apply>
<ci>v_ETCr6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_H_m</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_TCAr1</ci>
<ci>v_TCAr8</ci>
</apply>
<ci>v_TCAr9</ci>
</apply>
<ci>v_TCAr11</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">5</cn>
<ci>v_ETCr1</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dim">4</cn>
<ci>v_ETCr2</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dim">6</cn>
<ci>v_ETCr3</ci>
</apply>
</apply>
<ci>v_ETCr4</ci>
<apply>
<times/>
<cn cellml:units="dim">8</cn>
<ci>v_ETCr6</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ATP_m</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_TCAr9</ci>
<ci>v_TCAr11</ci>
</apply>
<ci>v_ETCr5</ci>
</apply>
<ci>v_ETCr6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ADP_m</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_TCAr9</ci>
</apply>
<ci>v_TCAr11</ci>
<ci>v_ETCr5</ci>
</apply>
<ci>v_ETCr6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_H2O_m</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_TCAr1</ci>
</apply>
<ci>v_TCAr7</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<ci>v_ETCr3</ci>
</apply>
<ci>v_ETCr6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_NADH_m</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<ci>v_TCAr3</ci>
<ci>v_TCAr4</ci>
<ci>v_TCAr8</ci>
<ci>v_TCAr10</ci>
</apply>
<ci>v_ETCr1</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_NAD_m</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<ci>v_TCAr3</ci>
</apply>
<ci>v_TCAr4</ci>
</apply>
<ci>v_TCAr8</ci>
</apply>
<ci>v_TCAr10</ci>
</apply>
<ci>v_ETCr1</ci>
</apply>
</apply>
</math>
<!-- Substrates associated with inner mitochondrial membrane -->
<variable initial_value="10" name="q_H_p" units="dim"/>
<variable initial_value="0" name="q_ATP_p" units="dim"/>
<variable initial_value="10" name="q_ADP_p" units="dim"/>
<variable initial_value="10" name="q_P_p" units="dim"/>
<variable initial_value="0" name="q_e" units="dim"/>
<variable initial_value="50" name="q_Oxy" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_H_p</ci>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dim">4</cn>
<ci>v_ETCr1</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">4</cn>
<ci>v_ETCr2</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<ci>v_ETCr3</ci>
</apply>
</apply>
<ci>v_ETCr4</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dim">9</cn>
<ci>v_ETCr6</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ATP_p</ci>
</apply>
<ci>v_ETCr5</ci>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_ADP_p</ci>
</apply>
<apply>
<minus/>
<ci>v_ETCr5</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_P_p</ci>
</apply>
<apply>
<minus/>
<ci>v_ETCr4</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_e</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<ci>v_ETCr1</ci>
<ci>v_ETCr3</ci>
</apply>
<ci>v_TCAr6</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>t</ci>
</bvar>
<ci>q_Oxy</ci>
</apply>
<apply>
<minus/>
<ci>v_ETCr3</ci>
</apply>
</apply>
</math>
<!-- Check conservation of chemical species-->
<!-- Carbon -->
<variable name="q_C1" units="dim"/>
<variable name="q_C2" units="dim"/>
<variable name="q_C3" units="dim"/>
<variable name="q_C4" units="dim"/>
<variable name="q_C5" units="dim"/>
<variable name="q_C6" units="dim"/>
<variable name="q_C9" units="dim"/>
<variable name="q_C15" units="dim"/>
<variable name="q_Carbon" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_C1</ci>
<apply>
<times/>
<cn cellml:units="dim">1</cn>
<apply>
<plus/>
<ci>q_CO2</ci>
<ci>q_HCO3_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_C2</ci>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<ci>q_AcetylCoA</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_C3</ci>
<apply>
<times/>
<cn cellml:units="dim">3</cn>
<apply>
<plus/>
<ci>q_GA3P</ci>
<ci>q_13bPG</ci>
<ci>q_3PG</ci>
<ci>q_2PG</ci>
<ci>q_DHAP</ci>
<ci>q_Gly3P</ci>
<ci>q_Glycerol</ci>
<ci>q_Lactate</ci>
<ci>q_Pyr_i</ci>
<ci>q_Pyr_m</ci>
<ci>q_PEP_i</ci>
<ci>q_PEP_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_C4</ci>
<apply>
<times/>
<cn cellml:units="dim">4</cn>
<apply>
<plus/>
<ci>q_SuccinylCoA</ci>
<ci>q_Succinate</ci>
<ci>q_Fumarate</ci>
<ci>q_Malate</ci>
<ci>q_Oxaloacetate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_C5</ci>
<apply>
<times/>
<cn cellml:units="dim">5</cn>
<ci>q_2OG</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_C6</ci>
<apply>
<times/>
<cn cellml:units="dim">6</cn>
<apply>
<plus/>
<ci>q_Glc</ci>
<ci>q_G6P</ci>
<ci>q_G1P</ci>
<ci>q_F6P</ci>
<ci>q_F16bP</ci>
<ci>q_Gly</ci>
<ci>q_Citrate</ci>
<ci>q_Isocitrate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_C9</ci>
<apply>
<times/>
<cn cellml:units="dim">9</cn>
<apply>
<plus/>
<ci>q_UTP</ci>
<ci>q_UDP</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_C15</ci>
<apply>
<times/>
<cn cellml:units="dim">15</cn>
<ci>q_UDPglc</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Carbon</ci>
<apply>
<plus/>
<ci>q_C1</ci>
<ci>q_C2</ci>
<ci>q_C3</ci>
<ci>q_C4</ci>
<ci>q_C5</ci>
<ci>q_C6</ci>
<ci>q_C9</ci>
<ci>q_C15</ci>
</apply>
</apply>
</math>
<!-- Oxygen -->
<variable name="q_O1" units="dim"/>
<variable name="q_O2" units="dim"/>
<variable name="q_O3" units="dim"/>
<variable name="q_O4" units="dim"/>
<variable name="q_O5" units="dim"/>
<variable name="q_O6" units="dim"/>
<variable name="q_O7" units="dim"/>
<variable name="q_O8" units="dim"/>
<variable name="q_O9" units="dim"/>
<variable name="q_O10" units="dim"/>
<variable name="q_O11" units="dim"/>
<variable name="q_O12" units="dim"/>
<variable name="q_O15" units="dim"/>
<variable name="q_O16" units="dim"/>
<variable name="q_Oxygen" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_O1</ci>
<apply>
<times/>
<cn cellml:units="dim">1</cn>
<apply>
<plus/>
<ci>q_H2O_i</ci>
<ci>q_H2O_m</ci>
<ci>q_AcetylCoA</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O2</ci>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<apply>
<plus/>
<ci>q_Oxy</ci>
<ci>q_CO2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O3</ci>
<apply>
<times/>
<cn cellml:units="dim">3</cn>
<apply>
<plus/>
<ci>q_HCO3_m</ci>
<ci>q_Pyr_i</ci>
<ci>q_Pyr_m</ci>
<ci>q_Lactate</ci>
<ci>q_Glycerol</ci>
<ci>q_ATP_i</ci>
<ci>q_ATP_m</ci>
<ci>q_ATP_p</ci>
<ci>q_GTP</ci>
<ci>q_SuccinylCoA</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O4</ci>
<apply>
<times/>
<cn cellml:units="dim">4</cn>
<apply>
<plus/>
<ci>q_P_i</ci>
<ci>q_P_m</ci>
<ci>q_P_p</ci>
<ci>q_Succinate</ci>
<ci>q_Fumarate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O5</ci>
<apply>
<times/>
<cn cellml:units="dim">5</cn>
<apply>
<plus/>
<ci>q_Gly</ci>
<ci>q_2OG</ci>
<ci>q_Malate</ci>
<ci>q_Oxaloacetate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O6</ci>
<apply>
<times/>
<cn cellml:units="dim">6</cn>
<apply>
<plus/>
<ci>q_Glc</ci>
<ci>q_GA3P</ci>
<ci>q_DHAP</ci>
<ci>q_PEP_i</ci>
<ci>q_PEP_m</ci>
<ci>q_Gly3P</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O7</ci>
<apply>
<times/>
<cn cellml:units="dim">7</cn>
<apply>
<plus/>
<ci>q_3PG</ci>
<ci>q_2PG</ci>
<ci>q_Citrate</ci>
<ci>q_Isocitrate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O8</ci>
<apply>
<times/>
<cn cellml:units="dim">8</cn>
<ci>q_PPi</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O9</ci>
<apply>
<times/>
<cn cellml:units="dim">9</cn>
<apply>
<plus/>
<ci>q_G6P</ci>
<ci>q_G1P</ci>
<ci>q_F6P</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O10</ci>
<apply>
<times/>
<cn cellml:units="dim">10</cn>
<ci>q_13bPG</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O11</ci>
<apply>
<times/>
<cn cellml:units="dim">11</cn>
<ci>q_UDP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O12</ci>
<apply>
<times/>
<cn cellml:units="dim">12</cn>
<ci>q_F16bP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O15</ci>
<apply>
<times/>
<cn cellml:units="dim">15</cn>
<ci>q_UTP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_O16</ci>
<apply>
<times/>
<cn cellml:units="dim">16</cn>
<ci>q_UDPglc</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Oxygen</ci>
<apply>
<plus/>
<ci>q_O1</ci>
<ci>q_O2</ci>
<ci>q_O3</ci>
<ci>q_O4</ci>
<ci>q_O5</ci>
<ci>q_O6</ci>
<ci>q_O7</ci>
<ci>q_O8</ci>
<ci>q_O9</ci>
<ci>q_O10</ci>
<ci>q_O11</ci>
<ci>q_O12</ci>
<ci>q_O15</ci>
<ci>q_O16</ci>
</apply>
</apply>
</math>
<!-- Phosphorus -->
<variable name="q_P1" units="dim"/>
<variable name="q_P2" units="dim"/>
<variable name="q_P3" units="dim"/>
<variable name="q_Phosphorus" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_P1</ci>
<apply>
<times/>
<cn cellml:units="dim">1</cn>
<apply>
<plus/>
<ci>q_P_i</ci>
<ci>q_P_m</ci>
<ci>q_P_p</ci>
<ci>q_G6P</ci>
<ci>q_G1P</ci>
<ci>q_F6P</ci>
<ci>q_GA3P</ci>
<ci>q_3PG</ci>
<ci>q_2PG</ci>
<ci>q_DHAP</ci>
<ci>q_Gly3P</ci>
<ci>q_PEP_i</ci>
<ci>q_PEP_m</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_P2</ci>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<apply>
<plus/>
<ci>q_ADP_i</ci>
<ci>q_ADP_m</ci>
<ci>q_GDP</ci>
<ci>q_UDPglc</ci>
<ci>q_UDP</ci>
<ci>q_PPi</ci>
<ci>q_F16bP</ci>
<ci>q_13bPG</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_P3</ci>
<apply>
<times/>
<cn cellml:units="dim">3</cn>
<apply>
<plus/>
<ci>q_ATP_i</ci>
<ci>q_ATP_m</ci>
<ci>q_GTP</ci>
<ci>q_UTP</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Phosphorus</ci>
<apply>
<plus/>
<ci>q_P1</ci>
<ci>q_P2</ci>
<ci>q_P3</ci>
</apply>
</apply>
</math>
<!-- Sulphur -->
<variable name="q_Sulphur" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_Sulphur</ci>
<apply>
<plus/>
<ci>q_AcetylCoA</ci>
<ci>q_CoASH</ci>
<ci>q_SuccinylCoA</ci>
</apply>
</apply>
</math>
<!-- Nitrogen -->
<variable name="q_Nitrogen" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_Nitrogen</ci>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<apply>
<plus/>
<ci>q_UTP</ci>
<ci>q_UDP</ci>
<ci>q_UDPglc</ci>
</apply>
</apply>
</apply>
</math>
<!-- Hydrogen -->
<variable name="q_H1" units="dim"/>
<variable name="q_H2" units="dim"/>
<variable name="q_H3" units="dim"/>
<variable name="q_H4" units="dim"/>
<variable name="q_H5" units="dim"/>
<variable name="q_H7" units="dim"/>
<variable name="q_H8" units="dim"/>
<variable name="q_H10" units="dim"/>
<variable name="q_H11" units="dim"/>
<variable name="q_H12" units="dim"/>
<variable name="q_H14" units="dim"/>
<variable name="q_H15" units="dim"/>
<variable name="q_H24" units="dim"/>
<variable name="q_Hydrogen" units="dim"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>q_H1</ci>
<apply>
<times/>
<cn cellml:units="dim">1</cn>
<apply>
<plus/>
<ci>q_H_i</ci>
<ci>q_H_m</ci>
<ci>q_H_p</ci>
<ci>q_HCO3_m</ci>
<ci>q_P_i</ci>
<ci>q_P_m</ci>
<ci>q_NADH_i</ci>
<ci>q_NADH_m</ci>
<ci>q_CoASH</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H2</ci>
<apply>
<times/>
<cn cellml:units="dim">2</cn>
<apply>
<plus/>
<ci>q_H2O_i</ci>
<ci>q_H2O_m</ci>
<ci>q_PEP_i</ci>
<ci>q_PEP_m</ci>
<ci>q_PPi</ci>
<ci>q_Fumarate</ci>
<ci>q_Oxaloacetate</ci>
<ci>q_FADH2</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H3</ci>
<apply>
<times/>
<cn cellml:units="dim">3</cn>
<apply>
<plus/>
<ci>q_Pyr_i</ci>
<ci>q_Pyr_m</ci>
<ci>q_AcetylCoA</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H4</ci>
<apply>
<times/>
<cn cellml:units="dim">4</cn>
<apply>
<plus/>
<ci>q_13bPG</ci>
<ci>q_3PG</ci>
<ci>q_2PG</ci>
<ci>q_2OG</ci>
<ci>q_SuccinylCoA</ci>
<ci>q_Succinate</ci>
<ci>q_Malate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H5</ci>
<apply>
<times/>
<cn cellml:units="dim">5</cn>
<apply>
<plus/>
<ci>q_GA3P</ci>
<ci>q_DHAP</ci>
<ci>q_Lactate</ci>
<ci>q_Citrate</ci>
<ci>q_Isocitrate</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H7</ci>
<apply>
<times/>
<cn cellml:units="dim">7</cn>
<ci>q_Gly3P</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H8</ci>
<apply>
<times/>
<cn cellml:units="dim">8</cn>
<ci>q_Glycerol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H10</ci>
<apply>
<times/>
<cn cellml:units="dim">10</cn>
<apply>
<plus/>
<ci>q_F16bP</ci>
<ci>q_Gly</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H11</ci>
<apply>
<times/>
<cn cellml:units="dim">11</cn>
<apply>
<plus/>
<ci>q_G6P</ci>
<ci>q_F6P</ci>
<ci>q_G1P</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H12</ci>
<apply>
<times/>
<cn cellml:units="dim">12</cn>
<ci>q_Glc</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H14</ci>
<apply>
<times/>
<cn cellml:units="dim">14</cn>
<ci>q_UDP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H15</ci>
<apply>
<times/>
<cn cellml:units="dim">15</cn>
<ci>q_UTP</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_H24</ci>
<apply>
<times/>
<cn cellml:units="dim">24</cn>
<ci>q_UDPglc</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>q_Hydrogen</ci>
<apply>
<plus/>
<ci>q_H1</ci>
<ci>q_H2</ci>
<ci>q_H3</ci>
<ci>q_H4</ci>
<ci>q_H5</ci>
<ci>q_H7</ci>
<ci>q_H8</ci>
<ci>q_H10</ci>
<ci>q_H11</ci>
<ci>q_H12</ci>
<ci>q_H14</ci>
<ci>q_H15</ci>
<ci>q_H24</ci>
</apply>
</apply>
</math>
</component>
</model>