A database of thermodynamic properties of the reactions of glycolysis, the tricarboxylic acid cycle, and the pentose phosphate pathway

Table 1.

Reactions with standard reaction enthalpies forumla estimated for I = 0

EC No.	Reaction name	Reaction abbreviation	Reference reaction	(kJ/mol)
EC 2.7.1.1	Glucokinase	GLK	GLC⁰+ ATP⁴⁻= G6P²⁻+ ADP³⁻+ H⁺	−23.8^a
EC 5.3.1.9	Phosphoglucose isomerase	PGI	G6P²⁻= F6P²⁻	11.53^b
EC 2.7.1.11	Phosphofructokinase	PFK	F6P²⁻+ ATP⁴⁻= F16P⁴⁻+ ADP³⁻+ H⁺	−9.5^c
EC 4.1.2.13	Fructose-1,6-biphosphatate aldolase	FBA	F16P⁴⁻= DHAP²⁻+ GAP²⁻	48.97^b
EC 5.3.1.1	Triosphosphate isomerase	TPI	GAP²⁻= DHAP²⁻	2.73^d
EC 4.1.2.13	Fructose-1,6-biphosphatate aldolase 2	FBA2	F16P⁴⁻= 2DHAP²⁻	51.70^b
EC 1.2.1.12	Glyceraldyde-3-P dehydrogenase	GAP	GAP²⁻+ HPO₄²⁻+ NAD_ox⁻= BPG⁴⁻+ NAD_red²⁻+ H⁺	#^e
EC 2.7.2.3	Phosphoglycerate kinase	PGK	GAP²⁻+ HPO₄²⁻+ NAD_ox⁻+ ADP³⁻= PG3³⁻+ NAD_red²⁻+ ATP⁴⁻+ H⁺	#^e
EC 5.4.2.1	Phosphoglycerate mutase	PGYM	PG2³⁻= PG3³⁻	28.05^b
EC 4.2.1.11	Enolase	ENO	PG2³⁻= PEP³⁻+ H₂O⁰	15.1^b
EC 2.7.1.40	Pyruvate kinase	PYK	PYR⁻+ ATP⁴⁻= PEP³⁻+ ADP³⁻+ H⁺	5.415^a
EC 4.1.3.7	Citrate synthase	CITS	OAA²⁻+ ACoA⁰+ H₂O⁰= CIT³⁻+ COAS⁻+ 2H⁺	#^e
EC 4.2.1.3	Aconitrate hydratese	ACON	ISCIT³⁻= CIT³⁻	−20.0^b
EC 1.1.1.42	Isocitrate dehydrogenase	IDH	ISCIT³⁻+ NADP_ox³⁻+ H₂O⁰= AKG²⁻+ NADP_red⁴⁻ + CO₃²⁻+ 2H⁺	−22.17^b
EC 6.2.1.4	Auccinate-CoA ligase	SCS	GTP^4-+ SUC²⁻+ COAS⁻+ H⁺= GDP³⁻+ HPO₄²⁻ + SUCCoA⁻	−30.9^b
EC 4.2.1.2	Fumarate hydratase	FUM	FUM²⁻+ H₂O⁰= MAL²⁻	−13.18^b
EC 1.1.1.37	Malate dehydrogenase	MDH	MAL²⁻+ NAD_ox⁻= OAA²⁻+ NAD_red²⁻+ H⁺	51.29^b
EC 2.7.4.6	Nucleoside-diphosphate kinase	NDK	ATP⁴⁻+ GDP³⁻= ADP³⁻+ GTP⁴⁻	#^e
EC 1.6.1.1	NADP transhydrogenase	NPTH	NAD_ox⁻+ NADP_red⁴⁻= NAD_red²⁻+ NADP_ox³⁻	−4.1^c
EC 1.1.1.40	Malic enzyme	MLE	MAL²⁻+ NADP_ox³⁻+ H₂O⁰= PYR⁻+ NADP_red⁴⁻+ CO₃^2-+ 2H⁺	#^e
EC 1.1.1.37	Malate dehydrogenase 2	MDH2	MAL²⁻+ ACoA⁰+ NAD_ox⁻+ H₂O⁰= CIT³⁻+ COAS⁻+ NAD_red²⁻+ 3H⁺	#^e
EC 2.7.1.23	NAD⁺ kinase	NADK	ATP⁴⁻+ NAD_ox⁻= ADP³⁻+ NADP_ox³⁻+ H⁺	#^e
EC 3.6.1.32	ATPase	ATPS	ATP⁴⁻+ H₂O⁰= ADP³⁻+ HPO₄²⁻+ H⁺	−20.5^a
EC 3.1.3.1	Alkaline phosphatase/G6P hydrolysis	G6PH	G6P²⁻+ H2O⁰= GLC⁰+ Pi²⁻	0.91^a
EC 6.4.1.1	Pyruvate carboxylase	PCL	PYR⁻+ ATP⁴⁻+ CO₃²⁻= OAA²⁻+ ADP³⁻+ Pi²⁻	#^e
EC 1.1.1.49	Glucose 6-phosphate dehydrogenase	G6PD	G6P²⁻+ NADP_ox³⁻= PGLT²⁻+ NADP_red⁴⁻+ H⁺	#^e
EC 3.1.1.31	6-Phosphogluconolactonase	PGL	G6P²⁻+ NADP_ox³⁻+ H₂O⁰= PGN³⁻+ NADP_red⁴⁻+ 2H⁺	#^e
EC 1.1.1.44	6-Phosphogluconate dehydrogenase	PGD	PGN³⁻+ NADP_ox³⁻+ H₂O⁰= RU5P²⁻+ NADP_red⁴⁻ + CO₃²⁻+ 2H⁺	37.47^b
EC 5.3.1.6	Ribose-5-phosphate isomerase	R5PI	R5P²⁻= RU5P²⁻	12.86^b
EC 5.1.3.1	Ribulose-phosphate 3-epimerase	RUPE	RU5P²⁻= X5P²⁻	#^e
EC 2.2.1.1	Transketolase	TKL	S7P²⁻+ GAP²⁻= R5P²⁻+ X5P²⁻	#^e
EC 2.2.1.2	Transaldolase	TAL	S7P²⁻+ GAP²⁻= E4P²⁻+ F6P²⁻	#^e
EC 2.2.1.1	Transketolase 2	TKL2	F6P²⁻+ GAP²⁻= E4P²⁻+ X5P²⁻	#^e
EC 1.2.4.1+EC 2.3.1.12+EC 1.8.1.4	Pyruvate dehydrogenase complex	PDH	PYR⁻+ COAS⁻+ NAD_ox⁻+ H₂O = CO₃²⁻+ ACoA⁰ + NAD_red²⁻+ H⁺	#^e
EC 1.1.1.41	Isocitrate dehydrogenase	IDH2	ISCIT³⁻+ NAD_ox⁻+ H₂O⁰= AKG^2-+ NAD_red²⁻+ CO₃²⁻+ 2H⁺	−26.27^d
EC 1.2.1.52	α-Ketoglutarate dehydrogenase	AKGDH	AKG²⁻+ NAD_ox⁻ + COAS⁻+ H₂O⁰= SUCCoA⁻+ NAD_red²⁻+ CO₃^2-+ H⁺	#^e
EC 1.3.5.1	Succinate dehydrogenase	SDH	SUC²⁻+ CoQ⁰= FUM²⁻+ CoQH₂⁰	#^e

EC No.	Reaction name	Reaction abbreviation	Reference reaction	(kJ/mol)
EC 2.7.1.1	Glucokinase	GLK	GLC⁰+ ATP⁴⁻= G6P²⁻+ ADP³⁻+ H⁺	−23.8^a
EC 5.3.1.9	Phosphoglucose isomerase	PGI	G6P²⁻= F6P²⁻	11.53^b
EC 2.7.1.11	Phosphofructokinase	PFK	F6P²⁻+ ATP⁴⁻= F16P⁴⁻+ ADP³⁻+ H⁺	−9.5^c
EC 4.1.2.13	Fructose-1,6-biphosphatate aldolase	FBA	F16P⁴⁻= DHAP²⁻+ GAP²⁻	48.97^b
EC 5.3.1.1	Triosphosphate isomerase	TPI	GAP²⁻= DHAP²⁻	2.73^d
EC 4.1.2.13	Fructose-1,6-biphosphatate aldolase 2	FBA2	F16P⁴⁻= 2DHAP²⁻	51.70^b
EC 1.2.1.12	Glyceraldyde-3-P dehydrogenase	GAP	GAP²⁻+ HPO₄²⁻+ NAD_ox⁻= BPG⁴⁻+ NAD_red²⁻+ H⁺	#^e
EC 2.7.2.3	Phosphoglycerate kinase	PGK	GAP²⁻+ HPO₄²⁻+ NAD_ox⁻+ ADP³⁻= PG3³⁻+ NAD_red²⁻+ ATP⁴⁻+ H⁺	#^e
EC 5.4.2.1	Phosphoglycerate mutase	PGYM	PG2³⁻= PG3³⁻	28.05^b
EC 4.2.1.11	Enolase	ENO	PG2³⁻= PEP³⁻+ H₂O⁰	15.1^b
EC 2.7.1.40	Pyruvate kinase	PYK	PYR⁻+ ATP⁴⁻= PEP³⁻+ ADP³⁻+ H⁺	5.415^a
EC 4.1.3.7	Citrate synthase	CITS	OAA²⁻+ ACoA⁰+ H₂O⁰= CIT³⁻+ COAS⁻+ 2H⁺	#^e
EC 4.2.1.3	Aconitrate hydratese	ACON	ISCIT³⁻= CIT³⁻	−20.0^b
EC 1.1.1.42	Isocitrate dehydrogenase	IDH	ISCIT³⁻+ NADP_ox³⁻+ H₂O⁰= AKG²⁻+ NADP_red⁴⁻ + CO₃²⁻+ 2H⁺	−22.17^b
EC 6.2.1.4	Auccinate-CoA ligase	SCS	GTP^4-+ SUC²⁻+ COAS⁻+ H⁺= GDP³⁻+ HPO₄²⁻ + SUCCoA⁻	−30.9^b
EC 4.2.1.2	Fumarate hydratase	FUM	FUM²⁻+ H₂O⁰= MAL²⁻	−13.18^b
EC 1.1.1.37	Malate dehydrogenase	MDH	MAL²⁻+ NAD_ox⁻= OAA²⁻+ NAD_red²⁻+ H⁺	51.29^b
EC 2.7.4.6	Nucleoside-diphosphate kinase	NDK	ATP⁴⁻+ GDP³⁻= ADP³⁻+ GTP⁴⁻	#^e
EC 1.6.1.1	NADP transhydrogenase	NPTH	NAD_ox⁻+ NADP_red⁴⁻= NAD_red²⁻+ NADP_ox³⁻	−4.1^c
EC 1.1.1.40	Malic enzyme	MLE	MAL²⁻+ NADP_ox³⁻+ H₂O⁰= PYR⁻+ NADP_red⁴⁻+ CO₃^2-+ 2H⁺	#^e
EC 1.1.1.37	Malate dehydrogenase 2	MDH2	MAL²⁻+ ACoA⁰+ NAD_ox⁻+ H₂O⁰= CIT³⁻+ COAS⁻+ NAD_red²⁻+ 3H⁺	#^e
EC 2.7.1.23	NAD⁺ kinase	NADK	ATP⁴⁻+ NAD_ox⁻= ADP³⁻+ NADP_ox³⁻+ H⁺	#^e
EC 3.6.1.32	ATPase	ATPS	ATP⁴⁻+ H₂O⁰= ADP³⁻+ HPO₄²⁻+ H⁺	−20.5^a
EC 3.1.3.1	Alkaline phosphatase/G6P hydrolysis	G6PH	G6P²⁻+ H2O⁰= GLC⁰+ Pi²⁻	0.91^a
EC 6.4.1.1	Pyruvate carboxylase	PCL	PYR⁻+ ATP⁴⁻+ CO₃²⁻= OAA²⁻+ ADP³⁻+ Pi²⁻	#^e
EC 1.1.1.49	Glucose 6-phosphate dehydrogenase	G6PD	G6P²⁻+ NADP_ox³⁻= PGLT²⁻+ NADP_red⁴⁻+ H⁺	#^e
EC 3.1.1.31	6-Phosphogluconolactonase	PGL	G6P²⁻+ NADP_ox³⁻+ H₂O⁰= PGN³⁻+ NADP_red⁴⁻+ 2H⁺	#^e
EC 1.1.1.44	6-Phosphogluconate dehydrogenase	PGD	PGN³⁻+ NADP_ox³⁻+ H₂O⁰= RU5P²⁻+ NADP_red⁴⁻ + CO₃²⁻+ 2H⁺	37.47^b
EC 5.3.1.6	Ribose-5-phosphate isomerase	R5PI	R5P²⁻= RU5P²⁻	12.86^b
EC 5.1.3.1	Ribulose-phosphate 3-epimerase	RUPE	RU5P²⁻= X5P²⁻	#^e
EC 2.2.1.1	Transketolase	TKL	S7P²⁻+ GAP²⁻= R5P²⁻+ X5P²⁻	#^e
EC 2.2.1.2	Transaldolase	TAL	S7P²⁻+ GAP²⁻= E4P²⁻+ F6P²⁻	#^e
EC 2.2.1.1	Transketolase 2	TKL2	F6P²⁻+ GAP²⁻= E4P²⁻+ X5P²⁻	#^e
EC 1.2.4.1+EC 2.3.1.12+EC 1.8.1.4	Pyruvate dehydrogenase complex	PDH	PYR⁻+ COAS⁻+ NAD_ox⁻+ H₂O = CO₃²⁻+ ACoA⁰ + NAD_red²⁻+ H⁺	#^e
EC 1.1.1.41	Isocitrate dehydrogenase	IDH2	ISCIT³⁻+ NAD_ox⁻+ H₂O⁰= AKG^2-+ NAD_red²⁻+ CO₃²⁻+ 2H⁺	−26.27^d
EC 1.2.1.52	α-Ketoglutarate dehydrogenase	AKGDH	AKG²⁻+ NAD_ox⁻ + COAS⁻+ H₂O⁰= SUCCoA⁻+ NAD_red²⁻+ CO₃^2-+ H⁺	#^e
EC 1.3.5.1	Succinate dehydrogenase	SDH	SUC²⁻+ CoQ⁰= FUM²⁻+ CoQH₂⁰	#^e

^aGoldberg et al. (21, 33).

^bCalculated value based on experimental data at different temperatures.

^cValues obtained from Goldberg et al. (20, 21) where associated ionic strength is not reported.

^dValue calculated from sum of dependent reactions.

^eValue not available.

Reaction enthalpies (⁠⁠) for two reactions in the pentose phosphate pathway (PGD and R5PI) can be estimated using van't Hoff equation because data on apparent equilibrium constants at different temperatures are available. Neither equilibrium data at different temperatures nor prior values of are available for the other six reactions of the pentose phosphate pathway; hence the symbol ‘#’ is used to denote the absence of data. For these cases, the value is set to zero in further calculations.

Database of reactant and dissociation constants

There are seven reactants introduced to the reactant database by adding the pentose phosphate pathway into the thermodynamic database: erythrose 4-phosphate (E4P), 6-phosphoglucono-δ-lactone (PGLT), 6-phospho-d-gluconate (PGN), ribose 5-phosphate (R5P), ribulose 5-phosphate (RU5P), sedoheptulose 7-phosphate (S7P) and xylulose 5-phosphate (X5P). For the five sugar phosphates (E4P, R5P, RU5P, S7P and X5P), cation dissociation constants for only R5P can be found in NIST database. Since E4P, RU5P, S7P and X5P are structurally similar substances to R5P (all have similar near neighbors to the phosphate group), a pragmatic approach to estimate the necessary dissociation constants is to use the values for R5P (24). Specifically, E4P, RU5P, S7P and X5P are assumed to have dissociation properties equal to those for R5P in our calculations. For PGLT, we use the pK_H₁ value of 5.99 from Alberty (10). For PGN, Casazza et al. (15) report that the hydrogen ion dissociation constant for the carboxylic acid is >1.0E-4, while the constant for the phosphate ester is the same as that of glucose 6-P [pK_H₁ = 5.99 (10)]. Here we arbitrarily assign the value of 4.995 to the pK_H₁ for PGN, which is the average of 4 and 5.99. The uncertainties of these assignments are considered in the uncertainty and sensitivity analysis below.

Basic thermodynamic and ion binding data for biochemical reactants and associated reference species are listed in Table 2. Each entry in the table contains the following information: (i) detailed name of reactant; (ii) reference species abbreviation; (iii) reactant abbreviation; (iv) number of protons in reference species; and (v) the dissociation constants (provided as pK) and the corresponding dissociation enthalpies Dissociation constants and enthalpies are tabulated at 298.15 K and 0.1 M ionic strength. The symbol ‘#’ is used to indicate absence of data. For these cases, the pK's are assumed to be infinite (no binding) with corresponding dissociation constants equal to zero. In the calculations, all the pK and values are adjusted to a common reference state of T = 298.15 K and I = 0.

Table 2.

Reactant database^a

Reactant	Reference species abbreviation^b	Reactant abbreviation	N_H^b	pK_H1		pK_H2		pK_Mg1		pK_HMg		pK_Mg2		pK_K1		pK_Na1		pK_HNa		pK_Ca1		pK_HCa
Acetyl-coenzyme A	ACoA⁰	ACoA	3	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Adenosine diphosphate	ADP³⁻	ADP	12	6.496	−2	3.87	16	3.3	−15	1.59	−7.5	1.27	−11.76	1	#	1.12	#	#	#	2.86	−9.6	1.48	−6.2
Adenosine triphosphate	ATP⁴⁻	ATP	12	6.71	−2	3.99	15	4.28	−18	2.32	−9.6	1.7	−17.52	1.17	−1	1.31	0.8	#	#	3.95	−13	2.16	−7.9
1,3-Bisphosphoglycerate	BPG⁴⁻	BPG	4	7.1^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Citrate	CIT³⁻	CIT	5	5.67	−1.9	4.35	3.1	3.517	−8	1.8	#	#	#	0.6	−3.54	0.75	−1	#	#	3.54	−1.2	2.07	#
Coenzyme A-SH	COAS⁻	COAS	0	8.17^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Carbon dioxide (total)	CO₃²⁻	CO2_tot	0	9.9^c	16.1^c	6.15^c	8.27^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Dihydroxyacetone phosphate	DHAP²⁻	DHAP	5	5.9	#	#	#	1.57	#	#	#	#	#	#	#	#	#	#	#	1.38	#	#	#
d-fructose 6-phosphate	F6P²⁻	F6P	11	5.89	−0.559^d	1.1	#	1.74^e	−9.72^d	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-fructose 1,6-phosphate	F16P⁴⁻	F16P	10	6.64	#	5.92	#	2.7	#	2.12	#	#	#	#	#	#	#	#	#	#	#	#	#
Fumarate	FUM²⁻	FUM	2	4.09	−1.56	2.86	1.08	#	#	#	#	#	#	#	#	#	#	#	#	0.6	−6.44	#	#
d-glucose 6-phosphate	G6P²⁻	G6P	11	5.89^e	−0.559^d	#	#	1.74^b	−9.72^d	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-glyceraldehyde 3-phosphate	GAP²⁻	GAP	5	5.27^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Guanosine diphosphate	GDP³⁻	GDP	12	6.505	−2.14	2.8	#	3.4	−7.1	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-glucose	GLC⁰	GLC	12	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Guanosine triphosphate	GTP⁴⁻	GTP	12	6.63	−3	2.93	7.1	4.31	−17	2.31	#	#	#	#	#	#	#	#	#	3.7	#	#	#
Water	H₂O⁰	H2O	2	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Isocitrate	ISCIT³⁻	ISCIT	5	5.765	#	4.29	#	2.625	#	1.43	#	#	#	#	#	#	#	#	#	2.54	#	#	#
α-Ketoglutarate	AKG²⁻	AKG	4	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Malate	MAL²⁻	MAL	4	4.715	−0.58	3.265	3.4	1.71	−6.16	0.9^f	#	#	#	0.18	−2.86	0.28	0.4	#	#	2.005	−1.06	1.06	8
NAD_ox	NAD_ox⁻	NAD_ox	26	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NAD_red	NAD_red²⁻	NAD_red	27	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NADP_ox	NADP_ox³⁻	NADP_ox	25	6.255^k	#	3.874^k	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NADP_red	NADP_red⁴⁻	NADP_red	26	6.255^l	#	3.874^l	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Oxaloacetate	OAA²⁻	OAA	2	3.9	5.24	2.26	16.62	1.02	#	#	#	#	#	#	#	#	#	#	#	1.6	#	#	#
Orthophosphate	HPO₄²⁻	Pi	1	6.78	4.6	1.945	−8.7	1.823	−9.518	0.669	#	#	#	0.5	#	0.61	#	0.0856	#	1.745	−9.518	0.921	−10.759
2-Phospho-d-glycerate	PG2³⁻	PG2	4	7	#	3.55	#	2.45	#	#	#	#	#	1.18	#	#	#	#	#	#	#	#	#
3-Phospho-d-glycerate	PG3³⁻	PG3	4	6.89^c	#	3.64^g	#	2.21^h	#	#	#	#	#	0.87 ^h	#	#	#	#	#	#	#	#	#
Phosphoenolpyruvate	PEP³⁻	PEP	2	6.245	#	3.45	#	2.26	#	#	#	#	#	1.08	#	#	#	#	#	#	#	#	#
Pyruvate	PYR⁻	PYR	3	2.26	12.8	#	#	1.1	#	#	#	#	#	#	#	#	#	#	#	0.8	#	#	#
Succinate	SUC²⁻	SUC	4	5.275	0.41	4.02	3	1.355	#	0.62	#	#	#	0.43	−2.76	0.4212	−2.759	#	#	1.405	−8.939	0.65	−8
Succinyl-CoA	SUCCoA⁻	SUCCoA	4	3.99^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Erythrose 4-phosphate	E4P²⁻	E4P	9	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
6-Phosphoglucono-δ-lactone	PGLT²⁻	PGLT	11	5.99^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
6-Phosphogluconate	PGN³⁻	PGN	13	4.995^j	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Ribose 5-phosphate	R5P²⁻	R5P	11	6.255	−9.75	2	#	1.58	−9.71^d	#	#	#	#	#	#	#	#	#	#	1.48	#	#	#
Ribulose 5-phosphate	RU5P²⁻	RU5P	11	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
Sedoheptulose 7-phosphate	S7P²⁻	S7P	15	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
Xylulose 5-phosphate	X5P²⁻	X5P	11	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#

Reactant	Reference species abbreviation^b	Reactant abbreviation	N_H^b	pK_H1		pK_H2		pK_Mg1		pK_HMg		pK_Mg2		pK_K1		pK_Na1		pK_HNa		pK_Ca1		pK_HCa
Acetyl-coenzyme A	ACoA⁰	ACoA	3	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Adenosine diphosphate	ADP³⁻	ADP	12	6.496	−2	3.87	16	3.3	−15	1.59	−7.5	1.27	−11.76	1	#	1.12	#	#	#	2.86	−9.6	1.48	−6.2
Adenosine triphosphate	ATP⁴⁻	ATP	12	6.71	−2	3.99	15	4.28	−18	2.32	−9.6	1.7	−17.52	1.17	−1	1.31	0.8	#	#	3.95	−13	2.16	−7.9
1,3-Bisphosphoglycerate	BPG⁴⁻	BPG	4	7.1^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Citrate	CIT³⁻	CIT	5	5.67	−1.9	4.35	3.1	3.517	−8	1.8	#	#	#	0.6	−3.54	0.75	−1	#	#	3.54	−1.2	2.07	#
Coenzyme A-SH	COAS⁻	COAS	0	8.17^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Carbon dioxide (total)	CO₃²⁻	CO2_tot	0	9.9^c	16.1^c	6.15^c	8.27^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Dihydroxyacetone phosphate	DHAP²⁻	DHAP	5	5.9	#	#	#	1.57	#	#	#	#	#	#	#	#	#	#	#	1.38	#	#	#
d-fructose 6-phosphate	F6P²⁻	F6P	11	5.89	−0.559^d	1.1	#	1.74^e	−9.72^d	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-fructose 1,6-phosphate	F16P⁴⁻	F16P	10	6.64	#	5.92	#	2.7	#	2.12	#	#	#	#	#	#	#	#	#	#	#	#	#
Fumarate	FUM²⁻	FUM	2	4.09	−1.56	2.86	1.08	#	#	#	#	#	#	#	#	#	#	#	#	0.6	−6.44	#	#
d-glucose 6-phosphate	G6P²⁻	G6P	11	5.89^e	−0.559^d	#	#	1.74^b	−9.72^d	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-glyceraldehyde 3-phosphate	GAP²⁻	GAP	5	5.27^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Guanosine diphosphate	GDP³⁻	GDP	12	6.505	−2.14	2.8	#	3.4	−7.1	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-glucose	GLC⁰	GLC	12	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Guanosine triphosphate	GTP⁴⁻	GTP	12	6.63	−3	2.93	7.1	4.31	−17	2.31	#	#	#	#	#	#	#	#	#	3.7	#	#	#
Water	H₂O⁰	H2O	2	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Isocitrate	ISCIT³⁻	ISCIT	5	5.765	#	4.29	#	2.625	#	1.43	#	#	#	#	#	#	#	#	#	2.54	#	#	#
α-Ketoglutarate	AKG²⁻	AKG	4	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Malate	MAL²⁻	MAL	4	4.715	−0.58	3.265	3.4	1.71	−6.16	0.9^f	#	#	#	0.18	−2.86	0.28	0.4	#	#	2.005	−1.06	1.06	8
NAD_ox	NAD_ox⁻	NAD_ox	26	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NAD_red	NAD_red²⁻	NAD_red	27	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NADP_ox	NADP_ox³⁻	NADP_ox	25	6.255^k	#	3.874^k	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NADP_red	NADP_red⁴⁻	NADP_red	26	6.255^l	#	3.874^l	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Oxaloacetate	OAA²⁻	OAA	2	3.9	5.24	2.26	16.62	1.02	#	#	#	#	#	#	#	#	#	#	#	1.6	#	#	#
Orthophosphate	HPO₄²⁻	Pi	1	6.78	4.6	1.945	−8.7	1.823	−9.518	0.669	#	#	#	0.5	#	0.61	#	0.0856	#	1.745	−9.518	0.921	−10.759
2-Phospho-d-glycerate	PG2³⁻	PG2	4	7	#	3.55	#	2.45	#	#	#	#	#	1.18	#	#	#	#	#	#	#	#	#
3-Phospho-d-glycerate	PG3³⁻	PG3	4	6.89^c	#	3.64^g	#	2.21^h	#	#	#	#	#	0.87 ^h	#	#	#	#	#	#	#	#	#
Phosphoenolpyruvate	PEP³⁻	PEP	2	6.245	#	3.45	#	2.26	#	#	#	#	#	1.08	#	#	#	#	#	#	#	#	#
Pyruvate	PYR⁻	PYR	3	2.26	12.8	#	#	1.1	#	#	#	#	#	#	#	#	#	#	#	0.8	#	#	#
Succinate	SUC²⁻	SUC	4	5.275	0.41	4.02	3	1.355	#	0.62	#	#	#	0.43	−2.76	0.4212	−2.759	#	#	1.405	−8.939	0.65	−8
Succinyl-CoA	SUCCoA⁻	SUCCoA	4	3.99^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Erythrose 4-phosphate	E4P²⁻	E4P	9	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
6-Phosphoglucono-δ-lactone	PGLT²⁻	PGLT	11	5.99^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
6-Phosphogluconate	PGN³⁻	PGN	13	4.995^j	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Ribose 5-phosphate	R5P²⁻	R5P	11	6.255	−9.75	2	#	1.58	−9.71^d	#	#	#	#	#	#	#	#	#	#	1.48	#	#	#
Ribulose 5-phosphate	RU5P²⁻	RU5P	11	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
Sedoheptulose 7-phosphate	S7P²⁻	S7P	15	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
Xylulose 5-phosphate	X5P²⁻	X5P	11	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#

^aDissociation pK and forumla are reported for T = 298.15 K and I = 0.1 M. Unless indicated, values are the average number obtained from NIST database (27). Dissociation enthalpies are reported in units of kJ/mol. Subscripts on ‘pK’ and ‘ forumla ’ entries are defined as follows: ‘H’: hydrogen, ‘Mg’: magnesium, ‘K’: potassium, ‘Na’: sodium, ‘Ca’: calcium, ‘1’: first ion dissociation, ‘2’: second ion dissociation, ‘HMg’: hydrogen ion binds to the ligand before magnesium ion binds to the ligand. ‘#’ denotes value is not available.

^bAlberty (8).

^cAlberty (10, 34).

^dTewari et al. (35).

^eG6P and F6P are assumed to have equivalent H⁺ and Mg²⁺-dissociation properties.

^fFrom NIST database (36) at T = 293.15 K.

^gLarsson-Raźnikiewicz (37).

^hMerrill et al. (38).

ⁱE4P, RU5P, S7P and X5P are assumed to have equivalent dissociation properties as R5P (24).

^jCasazza et al. (15).

^kBriggs et al. (39).

^lNADP_red is assumed to have equivalent dissociation properties as NADP_ox.

Table 2.

Reactant database^a

Reactant	Reference species abbreviation^b	Reactant abbreviation	N_H^b	pK_H1		pK_H2		pK_Mg1		pK_HMg		pK_Mg2		pK_K1		pK_Na1		pK_HNa		pK_Ca1		pK_HCa
Acetyl-coenzyme A	ACoA⁰	ACoA	3	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Adenosine diphosphate	ADP³⁻	ADP	12	6.496	−2	3.87	16	3.3	−15	1.59	−7.5	1.27	−11.76	1	#	1.12	#	#	#	2.86	−9.6	1.48	−6.2
Adenosine triphosphate	ATP⁴⁻	ATP	12	6.71	−2	3.99	15	4.28	−18	2.32	−9.6	1.7	−17.52	1.17	−1	1.31	0.8	#	#	3.95	−13	2.16	−7.9
1,3-Bisphosphoglycerate	BPG⁴⁻	BPG	4	7.1^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Citrate	CIT³⁻	CIT	5	5.67	−1.9	4.35	3.1	3.517	−8	1.8	#	#	#	0.6	−3.54	0.75	−1	#	#	3.54	−1.2	2.07	#
Coenzyme A-SH	COAS⁻	COAS	0	8.17^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Carbon dioxide (total)	CO₃²⁻	CO2_tot	0	9.9^c	16.1^c	6.15^c	8.27^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Dihydroxyacetone phosphate	DHAP²⁻	DHAP	5	5.9	#	#	#	1.57	#	#	#	#	#	#	#	#	#	#	#	1.38	#	#	#
d-fructose 6-phosphate	F6P²⁻	F6P	11	5.89	−0.559^d	1.1	#	1.74^e	−9.72^d	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-fructose 1,6-phosphate	F16P⁴⁻	F16P	10	6.64	#	5.92	#	2.7	#	2.12	#	#	#	#	#	#	#	#	#	#	#	#	#
Fumarate	FUM²⁻	FUM	2	4.09	−1.56	2.86	1.08	#	#	#	#	#	#	#	#	#	#	#	#	0.6	−6.44	#	#
d-glucose 6-phosphate	G6P²⁻	G6P	11	5.89^e	−0.559^d	#	#	1.74^b	−9.72^d	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-glyceraldehyde 3-phosphate	GAP²⁻	GAP	5	5.27^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Guanosine diphosphate	GDP³⁻	GDP	12	6.505	−2.14	2.8	#	3.4	−7.1	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-glucose	GLC⁰	GLC	12	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Guanosine triphosphate	GTP⁴⁻	GTP	12	6.63	−3	2.93	7.1	4.31	−17	2.31	#	#	#	#	#	#	#	#	#	3.7	#	#	#
Water	H₂O⁰	H2O	2	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Isocitrate	ISCIT³⁻	ISCIT	5	5.765	#	4.29	#	2.625	#	1.43	#	#	#	#	#	#	#	#	#	2.54	#	#	#
α-Ketoglutarate	AKG²⁻	AKG	4	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Malate	MAL²⁻	MAL	4	4.715	−0.58	3.265	3.4	1.71	−6.16	0.9^f	#	#	#	0.18	−2.86	0.28	0.4	#	#	2.005	−1.06	1.06	8
NAD_ox	NAD_ox⁻	NAD_ox	26	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NAD_red	NAD_red²⁻	NAD_red	27	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NADP_ox	NADP_ox³⁻	NADP_ox	25	6.255^k	#	3.874^k	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NADP_red	NADP_red⁴⁻	NADP_red	26	6.255^l	#	3.874^l	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Oxaloacetate	OAA²⁻	OAA	2	3.9	5.24	2.26	16.62	1.02	#	#	#	#	#	#	#	#	#	#	#	1.6	#	#	#
Orthophosphate	HPO₄²⁻	Pi	1	6.78	4.6	1.945	−8.7	1.823	−9.518	0.669	#	#	#	0.5	#	0.61	#	0.0856	#	1.745	−9.518	0.921	−10.759
2-Phospho-d-glycerate	PG2³⁻	PG2	4	7	#	3.55	#	2.45	#	#	#	#	#	1.18	#	#	#	#	#	#	#	#	#
3-Phospho-d-glycerate	PG3³⁻	PG3	4	6.89^c	#	3.64^g	#	2.21^h	#	#	#	#	#	0.87 ^h	#	#	#	#	#	#	#	#	#
Phosphoenolpyruvate	PEP³⁻	PEP	2	6.245	#	3.45	#	2.26	#	#	#	#	#	1.08	#	#	#	#	#	#	#	#	#
Pyruvate	PYR⁻	PYR	3	2.26	12.8	#	#	1.1	#	#	#	#	#	#	#	#	#	#	#	0.8	#	#	#
Succinate	SUC²⁻	SUC	4	5.275	0.41	4.02	3	1.355	#	0.62	#	#	#	0.43	−2.76	0.4212	−2.759	#	#	1.405	−8.939	0.65	−8
Succinyl-CoA	SUCCoA⁻	SUCCoA	4	3.99^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Erythrose 4-phosphate	E4P²⁻	E4P	9	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
6-Phosphoglucono-δ-lactone	PGLT²⁻	PGLT	11	5.99^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
6-Phosphogluconate	PGN³⁻	PGN	13	4.995^j	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Ribose 5-phosphate	R5P²⁻	R5P	11	6.255	−9.75	2	#	1.58	−9.71^d	#	#	#	#	#	#	#	#	#	#	1.48	#	#	#
Ribulose 5-phosphate	RU5P²⁻	RU5P	11	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
Sedoheptulose 7-phosphate	S7P²⁻	S7P	15	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
Xylulose 5-phosphate	X5P²⁻	X5P	11	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#

Reactant	Reference species abbreviation^b	Reactant abbreviation	N_H^b	pK_H1		pK_H2		pK_Mg1		pK_HMg		pK_Mg2		pK_K1		pK_Na1		pK_HNa		pK_Ca1		pK_HCa
Acetyl-coenzyme A	ACoA⁰	ACoA	3	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Adenosine diphosphate	ADP³⁻	ADP	12	6.496	−2	3.87	16	3.3	−15	1.59	−7.5	1.27	−11.76	1	#	1.12	#	#	#	2.86	−9.6	1.48	−6.2
Adenosine triphosphate	ATP⁴⁻	ATP	12	6.71	−2	3.99	15	4.28	−18	2.32	−9.6	1.7	−17.52	1.17	−1	1.31	0.8	#	#	3.95	−13	2.16	−7.9
1,3-Bisphosphoglycerate	BPG⁴⁻	BPG	4	7.1^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Citrate	CIT³⁻	CIT	5	5.67	−1.9	4.35	3.1	3.517	−8	1.8	#	#	#	0.6	−3.54	0.75	−1	#	#	3.54	−1.2	2.07	#
Coenzyme A-SH	COAS⁻	COAS	0	8.17^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Carbon dioxide (total)	CO₃²⁻	CO2_tot	0	9.9^c	16.1^c	6.15^c	8.27^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Dihydroxyacetone phosphate	DHAP²⁻	DHAP	5	5.9	#	#	#	1.57	#	#	#	#	#	#	#	#	#	#	#	1.38	#	#	#
d-fructose 6-phosphate	F6P²⁻	F6P	11	5.89	−0.559^d	1.1	#	1.74^e	−9.72^d	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-fructose 1,6-phosphate	F16P⁴⁻	F16P	10	6.64	#	5.92	#	2.7	#	2.12	#	#	#	#	#	#	#	#	#	#	#	#	#
Fumarate	FUM²⁻	FUM	2	4.09	−1.56	2.86	1.08	#	#	#	#	#	#	#	#	#	#	#	#	0.6	−6.44	#	#
d-glucose 6-phosphate	G6P²⁻	G6P	11	5.89^e	−0.559^d	#	#	1.74^b	−9.72^d	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-glyceraldehyde 3-phosphate	GAP²⁻	GAP	5	5.27^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Guanosine diphosphate	GDP³⁻	GDP	12	6.505	−2.14	2.8	#	3.4	−7.1	#	#	#	#	#	#	#	#	#	#	#	#	#	#
d-glucose	GLC⁰	GLC	12	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Guanosine triphosphate	GTP⁴⁻	GTP	12	6.63	−3	2.93	7.1	4.31	−17	2.31	#	#	#	#	#	#	#	#	#	3.7	#	#	#
Water	H₂O⁰	H2O	2	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Isocitrate	ISCIT³⁻	ISCIT	5	5.765	#	4.29	#	2.625	#	1.43	#	#	#	#	#	#	#	#	#	2.54	#	#	#
α-Ketoglutarate	AKG²⁻	AKG	4	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Malate	MAL²⁻	MAL	4	4.715	−0.58	3.265	3.4	1.71	−6.16	0.9^f	#	#	#	0.18	−2.86	0.28	0.4	#	#	2.005	−1.06	1.06	8
NAD_ox	NAD_ox⁻	NAD_ox	26	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NAD_red	NAD_red²⁻	NAD_red	27	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NADP_ox	NADP_ox³⁻	NADP_ox	25	6.255^k	#	3.874^k	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
NADP_red	NADP_red⁴⁻	NADP_red	26	6.255^l	#	3.874^l	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Oxaloacetate	OAA²⁻	OAA	2	3.9	5.24	2.26	16.62	1.02	#	#	#	#	#	#	#	#	#	#	#	1.6	#	#	#
Orthophosphate	HPO₄²⁻	Pi	1	6.78	4.6	1.945	−8.7	1.823	−9.518	0.669	#	#	#	0.5	#	0.61	#	0.0856	#	1.745	−9.518	0.921	−10.759
2-Phospho-d-glycerate	PG2³⁻	PG2	4	7	#	3.55	#	2.45	#	#	#	#	#	1.18	#	#	#	#	#	#	#	#	#
3-Phospho-d-glycerate	PG3³⁻	PG3	4	6.89^c	#	3.64^g	#	2.21^h	#	#	#	#	#	0.87 ^h	#	#	#	#	#	#	#	#	#
Phosphoenolpyruvate	PEP³⁻	PEP	2	6.245	#	3.45	#	2.26	#	#	#	#	#	1.08	#	#	#	#	#	#	#	#	#
Pyruvate	PYR⁻	PYR	3	2.26	12.8	#	#	1.1	#	#	#	#	#	#	#	#	#	#	#	0.8	#	#	#
Succinate	SUC²⁻	SUC	4	5.275	0.41	4.02	3	1.355	#	0.62	#	#	#	0.43	−2.76	0.4212	−2.759	#	#	1.405	−8.939	0.65	−8
Succinyl-CoA	SUCCoA⁻	SUCCoA	4	3.99^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Erythrose 4-phosphate	E4P²⁻	E4P	9	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
6-Phosphoglucono-δ-lactone	PGLT²⁻	PGLT	11	5.99^c	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
6-Phosphogluconate	PGN³⁻	PGN	13	4.995^j	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#	#
Ribose 5-phosphate	R5P²⁻	R5P	11	6.255	−9.75	2	#	1.58	−9.71^d	#	#	#	#	#	#	#	#	#	#	1.48	#	#	#
Ribulose 5-phosphate	RU5P²⁻	RU5P	11	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
Sedoheptulose 7-phosphate	S7P²⁻	S7P	15	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#
Xylulose 5-phosphate	X5P²⁻	X5P	11	6.255ⁱ	−9.75ⁱ	2ⁱ	#	1.58ⁱ	−9.71ⁱ	#	#	#	#	#	#	#	#	#	#	1.48ⁱ	#	#	#

^aDissociation pK and forumla are reported for T = 298.15 K and I = 0.1 M. Unless indicated, values are the average number obtained from NIST database (27). Dissociation enthalpies are reported in units of kJ/mol. Subscripts on ‘pK’ and ‘ forumla ’ entries are defined as follows: ‘H’: hydrogen, ‘Mg’: magnesium, ‘K’: potassium, ‘Na’: sodium, ‘Ca’: calcium, ‘1’: first ion dissociation, ‘2’: second ion dissociation, ‘HMg’: hydrogen ion binds to the ligand before magnesium ion binds to the ligand. ‘#’ denotes value is not available.

^bAlberty (8).

^cAlberty (10, 34).

^dTewari et al. (35).

^eG6P and F6P are assumed to have equivalent H⁺ and Mg²⁺-dissociation properties.

^fFrom NIST database (36) at T = 293.15 K.

^gLarsson-Raźnikiewicz (37).

^hMerrill et al. (38).

ⁱE4P, RU5P, S7P and X5P are assumed to have equivalent dissociation properties as R5P (24).

^jCasazza et al. (15).

^kBriggs et al. (39).

^lNADP_red is assumed to have equivalent dissociation properties as NADP_ox.

Estimation of standard-state thermodynamic quantities

Given the compiled raw experimental data on the reactions and reactants, the thermodynamic model is used to estimate reference and values for the reference reactions and species. Note that in our thermodynamic model the values are adjustable parameters estimated to obtain the best fit to the biochemical equilibrium data. As in previous studies (8), values of for oxidized species of certain redox pairs are arbitrarily set to zero. Thus, these values are not true free energies of the reactions of formation for these chemical species; instead, they are parameters for which the thermodynamic model makes optimal predictions for these interdependent biochemical reactions.

Because the calculated are interrelated [e.g. in the reductive pentose phosphate pathway, five carbon sugars (X5P and R5P) are converted into three carbon (GAP) and six carbon (F6P) sugars which can then be utilized by the pathway of glycolysis], the database of may be rigorously extended only by recalculating the entire database using all of the raw data. Therefore, model fitting is based on total 686 data entries for the network of all 33 reactions (the first 33 reactions listed in Table 1). Standard free energies of formation for all the reference species are unable to be estimated independently because there are 29 stoichiometrically independent reactions and in total 40 reactants in our system. There are four reactions PDH, IDH2, AKGDH and SDH (last four reactions in Tables 1 and 4) in the TCA cycle for which direct measurements are not available in the literature. For these reactions, values are calculated using Alberty's database (8) and set as constraints to perform the optimization of the thermodynamic model. IDH2 is not an independent reaction in the overall network of 33 reactions. Therefore, values of for 32 references species may be estimated from data on the 29 reactions with the three constraints. The values of for eight species are set to fixed values which are obtained from Alberty (8), as shown in Table 3 (values of for CoQ⁰ and CoQH₂⁰ do not come into these calculations). A constrained nonlinear optimization procedure with the fmincon solver (Mathworks, Inc.) is used to analyze the whole data set. By weighting in inverse proportion to the number of data points available for a given reaction and minimizing the difference between model predictions and experimental data, a simultaneous solution of standard reaction Gibbs energies is obtained for the entire data set.

Table 3.

Values of forumla (T = 298.15 K, I = 0) used in this study and that were taken from Alberty (8) (Table 3.2)

Species	(kJ/mol)
ACoA⁰	−188.52
ADP³⁻	−1906.13
CO₃²⁻	−527.81
GTP⁴⁻	−2768.1
H₂O	−237.19
NAD_ox⁻	0^a
HPO₄²⁻	−1096.1
H⁺	0
CoQ⁰	0^a
CoQH₂⁰	−89.92

^aProperty value is based on the arbitrary assignment of zero.

Table 3.

Values of forumla (T = 298.15 K, I = 0) used in this study and that were taken from Alberty (8) (Table 3.2)

Species	(kJ/mol)
ACoA⁰	−188.52
ADP³⁻	−1906.13
CO₃²⁻	−527.81
GTP⁴⁻	−2768.1
H₂O	−237.19
NAD_ox⁻	0^a
HPO₄²⁻	−1096.1
H⁺	0
CoQ⁰	0^a
CoQH₂⁰	−89.92

^aProperty value is based on the arbitrary assignment of zero.

Table 4.

Optimal reaction free energies for reference chemical reactions (⁠ forumla ⁠) and for biochemical reactions (⁠ forumla ⁠) under physiological conditions^a

EC No.	Reaction	^b	^b			Non-standard element contributions^c
EC No.	Reaction	^b	^b			T	I	pH	P
EC 2.7.1.1	GLK	16.19	16.03	−19.22	−35.41	1.61	1.57	−41.56	2.97
EC 5.3.1.9	PGI	3.12	3.13	2.78	−0.34	−0.34	0	0	0
EC 2.7.1.11	PFK	26.79	26.79	−15.62	−42.41	1.46	−4.70	−41.56	2.40
EC 4.1.2.13	FBA	18.79	18.80	24.64	5.85	−1.21	6.26	0	0.81
EC 5.3.1.1	TPI	−7.01	−7.01	−7.57	−0.56	−0.39	0	0	−0.17
EC 4.1.2.13	FBA2	11.79	11.79	17.07	5.28	−1.61	6.26	0	0.64
EC 1.2.1.12	GAP	51.37	51.37	2.60	−48.77	2.07	−9.39	−41.56	0.12
EC 2.7.2.3	PGK	34.37	34.37	−19.00	−53.37	1.38	−9.39	−41.56	−3.79
EC 5.4.2.1	PGYM	−5.89	−5.90	−6.35	−0.46	−1.37	0	0	0.91
EC 4.2.1.11	ENO	−4.54	−4.53	−4.47	0.07	−0.79	0	0	0.86
EC 2.7.1.40	PYK	66.91	66.90	27.18	−39.73	2.48	−1.57	−41.56	0.92
EC 4.1.3.7	CITS	60.16	60.32	−36.60	−96.76	2.42	−6.26	−83.13	−9.80
EC 4.2.1.3	ACON	−5.75	−5.76	−7.58	−1.83	0.57	0	0	−2.41
EC 1.1.1.42	IDH	97.05	97.06	−3.29	−100.34	4.80	−6.26	−83.13	−15.76
EC 6.2.1.4	SCS	−56.55	−56.56	0.07	56.62	−1.03	6.26	41.56	9.82
EC 4.2.1.2	FUM	−3.38	−3.38	−3.52	−0.14	0.39	0	0	−0.53
EC 1.1.1.37	MDH	71.41	71.09	28.04	−43.37	0.81	−3.13	−41.56	0.52
EC 2.7.4.6	NDK	0.02	0.01	−0.56	−0.58	0	0	0	−0.58
EC 1.6.1.1	NPTH	−3.53	−3.58	−0.41	3.12	0.02	3.13	0	−0.03
EC 1.1.1.40	MLE	104.53	103.45	2.00	−102.53	4.21	−7.83	−83.13	−15.78
EC 1.1.1.37	MDH2	131.57	131.41	−6.50	−138.07	5.30	−9.39	−124.69	−9.28
EC 2.7.1.23	NADK	29.11	26.90	−9.95	−39.06	1.17	−1.57	−41.56	2.89
EC 3.6.1.32	ATPS	4.99	4.67	−32.42	−37.41	1.03	1.57	−41.56	1.56
EC 3.1.3.1	G6PH	−11.20	−11.36	−13.10	-1.90	−0.49	0	0	−1.41
EC 6.4.1.1	PCL	−24.60	−24.12	−4.57	20.03	−0.99	3.13	0	17.89
EC 1.1.1.49	G6PD	38.69	#^d	−7.51	−46.20	1.56	−6.26	−41.56	0.07
EC 3.1.1.31	PGL	69.16	#^d	−21.89	−91.05	2.78	−10.96	−83.13	0.25
EC 1.1.1.44	PGD	104.64	#^d	1.23	−103.41	2.70	−6.26	−83.13	−16.73
EC 5.3.1.6	R5PI	0.99	1.2^e	0.52	−0.47	−0.47	0	0	0
EC 5.1.3.1	RUPE	−1.28	−1.21^e	−1.33	−0.05	−0.05	0	0	0
EC 2.2.1.1	TKL	1.58	#^d	1.23	−0.35	0.06	0	0	−0.41
EC 2.2.1.2	TAL	2.72	#^d	2.63	−0.09	0.11	0	0	−0.20
EC 2.2.1.1	TKL2	8.98	#^d	8.71	−0.27	0.36	0	0	−0.63
EC 1.2.4.1+EC 2.3.1.12+EC 1.8.1.4	PDH	15.78^f	17.50	−39.26	−55.04	0.64	−4.70	−41.56	−9.43
EC 1.1.1.41	IDH2	93.52	93.48	−3.70	−97.22	4.82	−3.13	−83.13	−15.79
EC 1.2.1.52	AKGDH	15.85^f	15.28	−37.66	−53.51	0.64	−3.13	−41.56	−9.45
EC 1.3.5.1	SDH	−1.35^f	−3.10	−0.59	0.76	−0.05	0	0	0.81

EC No.	Reaction	^b	^b			Non-standard element contributions^c
EC No.	Reaction	^b	^b			T	I	pH	P
EC 2.7.1.1	GLK	16.19	16.03	−19.22	−35.41	1.61	1.57	−41.56	2.97
EC 5.3.1.9	PGI	3.12	3.13	2.78	−0.34	−0.34	0	0	0
EC 2.7.1.11	PFK	26.79	26.79	−15.62	−42.41	1.46	−4.70	−41.56	2.40
EC 4.1.2.13	FBA	18.79	18.80	24.64	5.85	−1.21	6.26	0	0.81
EC 5.3.1.1	TPI	−7.01	−7.01	−7.57	−0.56	−0.39	0	0	−0.17
EC 4.1.2.13	FBA2	11.79	11.79	17.07	5.28	−1.61	6.26	0	0.64
EC 1.2.1.12	GAP	51.37	51.37	2.60	−48.77	2.07	−9.39	−41.56	0.12
EC 2.7.2.3	PGK	34.37	34.37	−19.00	−53.37	1.38	−9.39	−41.56	−3.79
EC 5.4.2.1	PGYM	−5.89	−5.90	−6.35	−0.46	−1.37	0	0	0.91
EC 4.2.1.11	ENO	−4.54	−4.53	−4.47	0.07	−0.79	0	0	0.86
EC 2.7.1.40	PYK	66.91	66.90	27.18	−39.73	2.48	−1.57	−41.56	0.92
EC 4.1.3.7	CITS	60.16	60.32	−36.60	−96.76	2.42	−6.26	−83.13	−9.80
EC 4.2.1.3	ACON	−5.75	−5.76	−7.58	−1.83	0.57	0	0	−2.41
EC 1.1.1.42	IDH	97.05	97.06	−3.29	−100.34	4.80	−6.26	−83.13	−15.76
EC 6.2.1.4	SCS	−56.55	−56.56	0.07	56.62	−1.03	6.26	41.56	9.82
EC 4.2.1.2	FUM	−3.38	−3.38	−3.52	−0.14	0.39	0	0	−0.53
EC 1.1.1.37	MDH	71.41	71.09	28.04	−43.37	0.81	−3.13	−41.56	0.52
EC 2.7.4.6	NDK	0.02	0.01	−0.56	−0.58	0	0	0	−0.58
EC 1.6.1.1	NPTH	−3.53	−3.58	−0.41	3.12	0.02	3.13	0	−0.03
EC 1.1.1.40	MLE	104.53	103.45	2.00	−102.53	4.21	−7.83	−83.13	−15.78
EC 1.1.1.37	MDH2	131.57	131.41	−6.50	−138.07	5.30	−9.39	−124.69	−9.28
EC 2.7.1.23	NADK	29.11	26.90	−9.95	−39.06	1.17	−1.57	−41.56	2.89
EC 3.6.1.32	ATPS	4.99	4.67	−32.42	−37.41	1.03	1.57	−41.56	1.56
EC 3.1.3.1	G6PH	−11.20	−11.36	−13.10	-1.90	−0.49	0	0	−1.41
EC 6.4.1.1	PCL	−24.60	−24.12	−4.57	20.03	−0.99	3.13	0	17.89
EC 1.1.1.49	G6PD	38.69	#^d	−7.51	−46.20	1.56	−6.26	−41.56	0.07
EC 3.1.1.31	PGL	69.16	#^d	−21.89	−91.05	2.78	−10.96	−83.13	0.25
EC 1.1.1.44	PGD	104.64	#^d	1.23	−103.41	2.70	−6.26	−83.13	−16.73
EC 5.3.1.6	R5PI	0.99	1.2^e	0.52	−0.47	−0.47	0	0	0
EC 5.1.3.1	RUPE	−1.28	−1.21^e	−1.33	−0.05	−0.05	0	0	0
EC 2.2.1.1	TKL	1.58	#^d	1.23	−0.35	0.06	0	0	−0.41
EC 2.2.1.2	TAL	2.72	#^d	2.63	−0.09	0.11	0	0	−0.20
EC 2.2.1.1	TKL2	8.98	#^d	8.71	−0.27	0.36	0	0	−0.63
EC 1.2.4.1+EC 2.3.1.12+EC 1.8.1.4	PDH	15.78^f	17.50	−39.26	−55.04	0.64	−4.70	−41.56	−9.43
EC 1.1.1.41	IDH2	93.52	93.48	−3.70	−97.22	4.82	−3.13	−83.13	−15.79
EC 1.2.1.52	AKGDH	15.85^f	15.28	−37.66	−53.51	0.64	−3.13	−41.56	−9.45
EC 1.3.5.1	SDH	−1.35^f	−3.10	−0.59	0.76	−0.05	0	0	0.81

^aThe physiological conditions (26) are as follows: T = 310.15 K, I = 0.18 M, pH = 7, [Mg²⁺] = 0.8 mM, [K⁺] = 140 mM, [Na⁺] = 10 mM, [Ca²⁺] = 0.0001 mM. The unit of the free energy is kJ/mol.

^b forumla is the optimal reaction free energies for reference chemical reactions in this work, forumla is the optimal reaction free energies for reference chemical reaction in Li et al. (7).

^cThis column lists the actual values of each non-standard physiological condition contribution (in kJ/mol) to the difference between the physiological free energies (⁠ forumla ⁠) and the standard free energies (⁠ forumla ⁠). T denotes the temperature contribution, I denotes the ionic strength contribution, pH denotes the pH contribution and P denotes the binding polynomial contribution.

^dValue is not available.

^eCalculated from Goldberg's database (25): the forumla values for R5P, RU5P and X5P in Goldberg's database are −1582.57, −1581.37 and −1582.58 kJ/mol, respectively.

^fValues calculated from Alberty's database (8) are used as model constraints for which there are no equilibrium data in the raw-data database. Since the forumla values for CoQ and CoQH2 are not predicted here, these values are set to 0 and −89.92 kJ/mol, respectively [from Alberty (8)] to computed forumla for the SDH reaction.

Table 4.

Optimal reaction free energies for reference chemical reactions (⁠ forumla ⁠) and for biochemical reactions (⁠ forumla ⁠) under physiological conditions^a

EC No.	Reaction	^b	^b			Non-standard element contributions^c
EC No.	Reaction	^b	^b			T	I	pH	P
EC 2.7.1.1	GLK	16.19	16.03	−19.22	−35.41	1.61	1.57	−41.56	2.97
EC 5.3.1.9	PGI	3.12	3.13	2.78	−0.34	−0.34	0	0	0
EC 2.7.1.11	PFK	26.79	26.79	−15.62	−42.41	1.46	−4.70	−41.56	2.40
EC 4.1.2.13	FBA	18.79	18.80	24.64	5.85	−1.21	6.26	0	0.81
EC 5.3.1.1	TPI	−7.01	−7.01	−7.57	−0.56	−0.39	0	0	−0.17
EC 4.1.2.13	FBA2	11.79	11.79	17.07	5.28	−1.61	6.26	0	0.64
EC 1.2.1.12	GAP	51.37	51.37	2.60	−48.77	2.07	−9.39	−41.56	0.12
EC 2.7.2.3	PGK	34.37	34.37	−19.00	−53.37	1.38	−9.39	−41.56	−3.79
EC 5.4.2.1	PGYM	−5.89	−5.90	−6.35	−0.46	−1.37	0	0	0.91
EC 4.2.1.11	ENO	−4.54	−4.53	−4.47	0.07	−0.79	0	0	0.86
EC 2.7.1.40	PYK	66.91	66.90	27.18	−39.73	2.48	−1.57	−41.56	0.92
EC 4.1.3.7	CITS	60.16	60.32	−36.60	−96.76	2.42	−6.26	−83.13	−9.80
EC 4.2.1.3	ACON	−5.75	−5.76	−7.58	−1.83	0.57	0	0	−2.41
EC 1.1.1.42	IDH	97.05	97.06	−3.29	−100.34	4.80	−6.26	−83.13	−15.76
EC 6.2.1.4	SCS	−56.55	−56.56	0.07	56.62	−1.03	6.26	41.56	9.82
EC 4.2.1.2	FUM	−3.38	−3.38	−3.52	−0.14	0.39	0	0	−0.53
EC 1.1.1.37	MDH	71.41	71.09	28.04	−43.37	0.81	−3.13	−41.56	0.52
EC 2.7.4.6	NDK	0.02	0.01	−0.56	−0.58	0	0	0	−0.58
EC 1.6.1.1	NPTH	−3.53	−3.58	−0.41	3.12	0.02	3.13	0	−0.03
EC 1.1.1.40	MLE	104.53	103.45	2.00	−102.53	4.21	−7.83	−83.13	−15.78
EC 1.1.1.37	MDH2	131.57	131.41	−6.50	−138.07	5.30	−9.39	−124.69	−9.28
EC 2.7.1.23	NADK	29.11	26.90	−9.95	−39.06	1.17	−1.57	−41.56	2.89
EC 3.6.1.32	ATPS	4.99	4.67	−32.42	−37.41	1.03	1.57	−41.56	1.56
EC 3.1.3.1	G6PH	−11.20	−11.36	−13.10	-1.90	−0.49	0	0	−1.41
EC 6.4.1.1	PCL	−24.60	−24.12	−4.57	20.03	−0.99	3.13	0	17.89
EC 1.1.1.49	G6PD	38.69	#^d	−7.51	−46.20	1.56	−6.26	−41.56	0.07
EC 3.1.1.31	PGL	69.16	#^d	−21.89	−91.05	2.78	−10.96	−83.13	0.25
EC 1.1.1.44	PGD	104.64	#^d	1.23	−103.41	2.70	−6.26	−83.13	−16.73
EC 5.3.1.6	R5PI	0.99	1.2^e	0.52	−0.47	−0.47	0	0	0
EC 5.1.3.1	RUPE	−1.28	−1.21^e	−1.33	−0.05	−0.05	0	0	0
EC 2.2.1.1	TKL	1.58	#^d	1.23	−0.35	0.06	0	0	−0.41
EC 2.2.1.2	TAL	2.72	#^d	2.63	−0.09	0.11	0	0	−0.20
EC 2.2.1.1	TKL2	8.98	#^d	8.71	−0.27	0.36	0	0	−0.63
EC 1.2.4.1+EC 2.3.1.12+EC 1.8.1.4	PDH	15.78^f	17.50	−39.26	−55.04	0.64	−4.70	−41.56	−9.43
EC 1.1.1.41	IDH2	93.52	93.48	−3.70	−97.22	4.82	−3.13	−83.13	−15.79
EC 1.2.1.52	AKGDH	15.85^f	15.28	−37.66	−53.51	0.64	−3.13	−41.56	−9.45
EC 1.3.5.1	SDH	−1.35^f	−3.10	−0.59	0.76	−0.05	0	0	0.81

EC No.	Reaction	^b	^b			Non-standard element contributions^c
EC No.	Reaction	^b	^b			T	I	pH	P
EC 2.7.1.1	GLK	16.19	16.03	−19.22	−35.41	1.61	1.57	−41.56	2.97
EC 5.3.1.9	PGI	3.12	3.13	2.78	−0.34	−0.34	0	0	0
EC 2.7.1.11	PFK	26.79	26.79	−15.62	−42.41	1.46	−4.70	−41.56	2.40
EC 4.1.2.13	FBA	18.79	18.80	24.64	5.85	−1.21	6.26	0	0.81
EC 5.3.1.1	TPI	−7.01	−7.01	−7.57	−0.56	−0.39	0	0	−0.17
EC 4.1.2.13	FBA2	11.79	11.79	17.07	5.28	−1.61	6.26	0	0.64
EC 1.2.1.12	GAP	51.37	51.37	2.60	−48.77	2.07	−9.39	−41.56	0.12
EC 2.7.2.3	PGK	34.37	34.37	−19.00	−53.37	1.38	−9.39	−41.56	−3.79
EC 5.4.2.1	PGYM	−5.89	−5.90	−6.35	−0.46	−1.37	0	0	0.91
EC 4.2.1.11	ENO	−4.54	−4.53	−4.47	0.07	−0.79	0	0	0.86
EC 2.7.1.40	PYK	66.91	66.90	27.18	−39.73	2.48	−1.57	−41.56	0.92
EC 4.1.3.7	CITS	60.16	60.32	−36.60	−96.76	2.42	−6.26	−83.13	−9.80
EC 4.2.1.3	ACON	−5.75	−5.76	−7.58	−1.83	0.57	0	0	−2.41
EC 1.1.1.42	IDH	97.05	97.06	−3.29	−100.34	4.80	−6.26	−83.13	−15.76
EC 6.2.1.4	SCS	−56.55	−56.56	0.07	56.62	−1.03	6.26	41.56	9.82
EC 4.2.1.2	FUM	−3.38	−3.38	−3.52	−0.14	0.39	0	0	−0.53
EC 1.1.1.37	MDH	71.41	71.09	28.04	−43.37	0.81	−3.13	−41.56	0.52
EC 2.7.4.6	NDK	0.02	0.01	−0.56	−0.58	0	0	0	−0.58
EC 1.6.1.1	NPTH	−3.53	−3.58	−0.41	3.12	0.02	3.13	0	−0.03
EC 1.1.1.40	MLE	104.53	103.45	2.00	−102.53	4.21	−7.83	−83.13	−15.78
EC 1.1.1.37	MDH2	131.57	131.41	−6.50	−138.07	5.30	−9.39	−124.69	−9.28
EC 2.7.1.23	NADK	29.11	26.90	−9.95	−39.06	1.17	−1.57	−41.56	2.89
EC 3.6.1.32	ATPS	4.99	4.67	−32.42	−37.41	1.03	1.57	−41.56	1.56
EC 3.1.3.1	G6PH	−11.20	−11.36	−13.10	-1.90	−0.49	0	0	−1.41
EC 6.4.1.1	PCL	−24.60	−24.12	−4.57	20.03	−0.99	3.13	0	17.89
EC 1.1.1.49	G6PD	38.69	#^d	−7.51	−46.20	1.56	−6.26	−41.56	0.07
EC 3.1.1.31	PGL	69.16	#^d	−21.89	−91.05	2.78	−10.96	−83.13	0.25
EC 1.1.1.44	PGD	104.64	#^d	1.23	−103.41	2.70	−6.26	−83.13	−16.73
EC 5.3.1.6	R5PI	0.99	1.2^e	0.52	−0.47	−0.47	0	0	0
EC 5.1.3.1	RUPE	−1.28	−1.21^e	−1.33	−0.05	−0.05	0	0	0
EC 2.2.1.1	TKL	1.58	#^d	1.23	−0.35	0.06	0	0	−0.41
EC 2.2.1.2	TAL	2.72	#^d	2.63	−0.09	0.11	0	0	−0.20
EC 2.2.1.1	TKL2	8.98	#^d	8.71	−0.27	0.36	0	0	−0.63
EC 1.2.4.1+EC 2.3.1.12+EC 1.8.1.4	PDH	15.78^f	17.50	−39.26	−55.04	0.64	−4.70	−41.56	−9.43
EC 1.1.1.41	IDH2	93.52	93.48	−3.70	−97.22	4.82	−3.13	−83.13	−15.79
EC 1.2.1.52	AKGDH	15.85^f	15.28	−37.66	−53.51	0.64	−3.13	−41.56	−9.45
EC 1.3.5.1	SDH	−1.35^f	−3.10	−0.59	0.76	−0.05	0	0	0.81

^aThe physiological conditions (26) are as follows: T = 310.15 K, I = 0.18 M, pH = 7, [Mg²⁺] = 0.8 mM, [K⁺] = 140 mM, [Na⁺] = 10 mM, [Ca²⁺] = 0.0001 mM. The unit of the free energy is kJ/mol.

^b forumla is the optimal reaction free energies for reference chemical reactions in this work, forumla is the optimal reaction free energies for reference chemical reaction in Li et al. (7).

^cThis column lists the actual values of each non-standard physiological condition contribution (in kJ/mol) to the difference between the physiological free energies (⁠ forumla ⁠) and the standard free energies (⁠ forumla ⁠). T denotes the temperature contribution, I denotes the ionic strength contribution, pH denotes the pH contribution and P denotes the binding polynomial contribution.

^dValue is not available.

^eCalculated from Goldberg's database (25): the forumla values for R5P, RU5P and X5P in Goldberg's database are −1582.57, −1581.37 and −1582.58 kJ/mol, respectively.

^fValues calculated from Alberty's database (8) are used as model constraints for which there are no equilibrium data in the raw-data database. Since the forumla values for CoQ and CoQH2 are not predicted here, these values are set to 0 and −89.92 kJ/mol, respectively [from Alberty (8)] to computed forumla for the SDH reaction.

Open in new tab Download slide

Results

Estimated Gibbs free energies of reaction and formation

Figure 1 illustrates model predictions versus experimental data for all data used in the analysis. The predicted values are obtained based on accounting for the biochemical state associated with a given experimental measurement in the raw-data database. Data points in the pentose phosphate pathway are shown as filled squares. Other data points are shown as open squares. The data span almost 14 orders of magnitude, and reveal good agreement between model predictions and measured data.

Figure 1.

Model-predicted versus experimental . Model predicted apparent equilibrium constants under defined experimental conditions (T, I, [Mg2+], [Ca2+], [Na+], [K+] and pH) are plotted versus experimental measurements for all data used in the analysis. Data points in the pentose phosphate pathway are shown as filled squares.

Model-predicted forumla versus experimental forumla ⁠. Model predicted apparent equilibrium constants under defined experimental conditions (T, I, [Mg²⁺], [Ca²⁺], [Na⁺], [K⁺] and pH) are plotted versus experimental measurements for all data used in the analysis. Data points in the pentose phosphate pathway are shown as filled squares.

The optimal estimates of listed in Table 4 (⁠⁠), compared to the values from Li et al. (7) (⁠⁠). of first 25 reactions are almost the same as the results obtained previously (7). For the remaining eight reactions of the pentose phosphate pathway, since the of the reference species for E4P, PGLT, PGN and S7P are not available in the Alberty (8) or Goldberg's database (25), the symbol ‘#’ is used to denote the lack of a value. For reactions R5PI and RUPE, the absolute difference between our current model predictions and the values reported by Goldberg et al. (25) is within a reasonable margin, 0.21 and 0.07 kJ/mol, respectively. Predicted versus experimental measures are plotted in Figure 2 for these two reactions (R5PI and RUPE). The current database and the Goldberg database yield similar results. Recall that the dissociation properties of R5P, RU5P and X5P are assumed to be the same; the binding polynomials for the reactants on the left- and right-hand sides of these reactions are identical. Thus, the computed apparent equilibrium constants for these reactions do not depend on the experimental ionic composition. The validity of this assumption will be considered in the uncertainty and sensitivity analysis below. In addition, the predicted thermodynamic properties are not affected by the ionic strength for these reactions.

Figure 2.

Model-predicted versus experimental for the ribose-5-phosphate isomerase and ribuloase-phosphate 3-epimerase reactions. Open circles (GT-based ) are computed based on the Goldberg's database (25); filled squares (optimized ) are computed based on the optimized values of from Table 5.

Open in new tab Download slide

Model-predicted forumla versus experimental forumla for the ribose-5-phosphate isomerase and ribuloase-phosphate 3-epimerase reactions. Open circles (GT-based forumla ⁠) are computed based on the Goldberg's database (25); filled squares (optimized forumla ⁠) are computed based on the optimized values of forumla from Table 5.

Optimal predicted associated with the predictions are listed in Table 5, compared to the optimal values of our previous version of the database (7). After the seven reactants of the pentose phosphate pathway and three constraints (⁠ of reaction PDH, AKGDH and SDH) are introduced, most of the estimated values of are shifted substantially. However, from Table 4 we can see that these shifts do not change the predicted compared to previous predictions. This is because the shifts in do not change the optimized results (apparent or ⁠) of our model. Recall that here the estimated values represent parameters in a thermodynamic model for this set of interdependent biochemical reactions. Since the number of independent reactions for which data exist (29 reactions) is smaller than the number of values to represent the system (40 reactants), values for several reference species are set to either existing values or values reported elsewhere (Table 3). As a result, these values are not physical constants. Rather, they are parameters in a thermodynamic network model that together form a self-consistent picture of the thermodynamics of reactions of the set of reactants studied here.

Table 5.

Optimal predicted forumla for reference species (T = 298.15 K, I = 0 M)

No.	Species	^a	^a		Sensitivity
1	GLC⁰	−719.37	−916.39	197.02	0.33
2	ATP⁴⁻	−2770.03	−2769.71	0.32	56099.52
3	F6P²⁻	−1563.96	−1760.81	196.85	1.52
4	F16P⁴⁻	−2401.08	−2597.60	196.52	3.55
5	DHAP²⁻	−1194.65	−1292.91	98.26	3.51
6	GAP²⁻	−1187.64	−1285.90	98.26	3.49
7	BPG⁴⁻	−2276.28	−2354.55	78.27	19.44
8	NAD_red²⁻	43.91	23.91	20.00	0.14
9	PG3³⁻	−1429.38	−1507.96	78.58	7.83
10	PG2³⁻	−1423.49	−1502.06	78.57	7.83
11	PEP³⁻	−1190.84	−1269.40	78.56	5.54
12	PYR⁻	−393.85	−472.72	78.87	0.64
13	OAA²⁻	−714.06	−792.13	78.07	2.06
14	CIT³⁻	−1022.44	−1157.52	135.08	1.38
15	ISCIT³⁻	−1016.69	−1151.76	135.07	1.38
16	NADP_ox³⁻	−834.79	−836.68	1.89	1837.57
17	AKG²⁻	−676.45	−791.57	115.12	0.83
18	SUC²⁻	−589.56	−685.56	96.00	0.91
19	GDP³⁻	−1904.22	−1904.53	0.31	9560.62
20	FUM²⁻	−500.99	−598.74	97.75	0.67
21	MAL²⁻	−741.56	−839.30	97.74	1.44
22	G6P²⁻	−1567.08	−1763.94^b	196.86	1.52
23	COAS⁰	−57.17	0^b	57.17	0.03
24	NADP_red⁴⁻	−787.35	−809.19^b	21.84	41.94
25	SUCCoA⁻	−471.06	−509.59^b	38.53	3.47
26	E4P²⁻	−1306.62	#	#	2.39
27	PGLT²⁻	−1575.83	#	#	1.85
28	PGN³⁻	−1782.55	#	#	2.34
29	R5P²⁻	−1435.72	#	#	1.86
30	RU5P²⁻	−1434.72	#	#	1.86
31	S7P²⁻	−1685.66	#	#	1.31
32	X5P²⁻	−1436.00	#	#	1.85

No.	Species	^a	^a		Sensitivity
1	GLC⁰	−719.37	−916.39	197.02	0.33
2	ATP⁴⁻	−2770.03	−2769.71	0.32	56099.52
3	F6P²⁻	−1563.96	−1760.81	196.85	1.52
4	F16P⁴⁻	−2401.08	−2597.60	196.52	3.55
5	DHAP²⁻	−1194.65	−1292.91	98.26	3.51
6	GAP²⁻	−1187.64	−1285.90	98.26	3.49
7	BPG⁴⁻	−2276.28	−2354.55	78.27	19.44
8	NAD_red²⁻	43.91	23.91	20.00	0.14
9	PG3³⁻	−1429.38	−1507.96	78.58	7.83
10	PG2³⁻	−1423.49	−1502.06	78.57	7.83
11	PEP³⁻	−1190.84	−1269.40	78.56	5.54
12	PYR⁻	−393.85	−472.72	78.87	0.64
13	OAA²⁻	−714.06	−792.13	78.07	2.06
14	CIT³⁻	−1022.44	−1157.52	135.08	1.38
15	ISCIT³⁻	−1016.69	−1151.76	135.07	1.38
16	NADP_ox³⁻	−834.79	−836.68	1.89	1837.57
17	AKG²⁻	−676.45	−791.57	115.12	0.83
18	SUC²⁻	−589.56	−685.56	96.00	0.91
19	GDP³⁻	−1904.22	−1904.53	0.31	9560.62
20	FUM²⁻	−500.99	−598.74	97.75	0.67
21	MAL²⁻	−741.56	−839.30	97.74	1.44
22	G6P²⁻	−1567.08	−1763.94^b	196.86	1.52
23	COAS⁰	−57.17	0^b	57.17	0.03
24	NADP_red⁴⁻	−787.35	−809.19^b	21.84	41.94
25	SUCCoA⁻	−471.06	−509.59^b	38.53	3.47
26	E4P²⁻	−1306.62	#	#	2.39
27	PGLT²⁻	−1575.83	#	#	1.85
28	PGN³⁻	−1782.55	#	#	2.34
29	R5P²⁻	−1435.72	#	#	1.86
30	RU5P²⁻	−1434.72	#	#	1.86
31	S7P²⁻	−1685.66	#	#	1.31
32	X5P²⁻	−1436.00	#	#	1.85

^a forumla is the optimal free energies of formation in this work, forumla is optimal free energies of formation in Li et al. (7). The unit of the free energy is kJ/mol.

^b forumla of species is set as fixed value which is obtained from Alberty's database as in Li et al. (7).

Table 5.

Optimal predicted forumla for reference species (T = 298.15 K, I = 0 M)

No.	Species	^a	^a		Sensitivity
1	GLC⁰	−719.37	−916.39	197.02	0.33
2	ATP⁴⁻	−2770.03	−2769.71	0.32	56099.52
3	F6P²⁻	−1563.96	−1760.81	196.85	1.52
4	F16P⁴⁻	−2401.08	−2597.60	196.52	3.55
5	DHAP²⁻	−1194.65	−1292.91	98.26	3.51
6	GAP²⁻	−1187.64	−1285.90	98.26	3.49
7	BPG⁴⁻	−2276.28	−2354.55	78.27	19.44
8	NAD_red²⁻	43.91	23.91	20.00	0.14
9	PG3³⁻	−1429.38	−1507.96	78.58	7.83
10	PG2³⁻	−1423.49	−1502.06	78.57	7.83
11	PEP³⁻	−1190.84	−1269.40	78.56	5.54
12	PYR⁻	−393.85	−472.72	78.87	0.64
13	OAA²⁻	−714.06	−792.13	78.07	2.06
14	CIT³⁻	−1022.44	−1157.52	135.08	1.38
15	ISCIT³⁻	−1016.69	−1151.76	135.07	1.38
16	NADP_ox³⁻	−834.79	−836.68	1.89	1837.57
17	AKG²⁻	−676.45	−791.57	115.12	0.83
18	SUC²⁻	−589.56	−685.56	96.00	0.91
19	GDP³⁻	−1904.22	−1904.53	0.31	9560.62
20	FUM²⁻	−500.99	−598.74	97.75	0.67
21	MAL²⁻	−741.56	−839.30	97.74	1.44
22	G6P²⁻	−1567.08	−1763.94^b	196.86	1.52
23	COAS⁰	−57.17	0^b	57.17	0.03
24	NADP_red⁴⁻	−787.35	−809.19^b	21.84	41.94
25	SUCCoA⁻	−471.06	−509.59^b	38.53	3.47
26	E4P²⁻	−1306.62	#	#	2.39
27	PGLT²⁻	−1575.83	#	#	1.85
28	PGN³⁻	−1782.55	#	#	2.34
29	R5P²⁻	−1435.72	#	#	1.86
30	RU5P²⁻	−1434.72	#	#	1.86
31	S7P²⁻	−1685.66	#	#	1.31
32	X5P²⁻	−1436.00	#	#	1.85

No.	Species	^a	^a		Sensitivity
1	GLC⁰	−719.37	−916.39	197.02	0.33
2	ATP⁴⁻	−2770.03	−2769.71	0.32	56099.52
3	F6P²⁻	−1563.96	−1760.81	196.85	1.52
4	F16P⁴⁻	−2401.08	−2597.60	196.52	3.55
5	DHAP²⁻	−1194.65	−1292.91	98.26	3.51
6	GAP²⁻	−1187.64	−1285.90	98.26	3.49
7	BPG⁴⁻	−2276.28	−2354.55	78.27	19.44
8	NAD_red²⁻	43.91	23.91	20.00	0.14
9	PG3³⁻	−1429.38	−1507.96	78.58	7.83
10	PG2³⁻	−1423.49	−1502.06	78.57	7.83
11	PEP³⁻	−1190.84	−1269.40	78.56	5.54
12	PYR⁻	−393.85	−472.72	78.87	0.64
13	OAA²⁻	−714.06	−792.13	78.07	2.06
14	CIT³⁻	−1022.44	−1157.52	135.08	1.38
15	ISCIT³⁻	−1016.69	−1151.76	135.07	1.38
16	NADP_ox³⁻	−834.79	−836.68	1.89	1837.57
17	AKG²⁻	−676.45	−791.57	115.12	0.83
18	SUC²⁻	−589.56	−685.56	96.00	0.91
19	GDP³⁻	−1904.22	−1904.53	0.31	9560.62
20	FUM²⁻	−500.99	−598.74	97.75	0.67
21	MAL²⁻	−741.56	−839.30	97.74	1.44
22	G6P²⁻	−1567.08	−1763.94^b	196.86	1.52
23	COAS⁰	−57.17	0^b	57.17	0.03
24	NADP_red⁴⁻	−787.35	−809.19^b	21.84	41.94
25	SUCCoA⁻	−471.06	−509.59^b	38.53	3.47
26	E4P²⁻	−1306.62	#	#	2.39
27	PGLT²⁻	−1575.83	#	#	1.85
28	PGN³⁻	−1782.55	#	#	2.34
29	R5P²⁻	−1435.72	#	#	1.86
30	RU5P²⁻	−1434.72	#	#	1.86
31	S7P²⁻	−1685.66	#	#	1.31
32	X5P²⁻	−1436.00	#	#	1.85

^a forumla is the optimal free energies of formation in this work, forumla is optimal free energies of formation in Li et al. (7). The unit of the free energy is kJ/mol.

^b forumla of species is set as fixed value which is obtained from Alberty's database as in Li et al. (7).

The robustness of these calculations is checked by repeating the optimization with the

for a single species constrained to a value

different from the estimated optimal value. The degree to which the experimental data can be matched with one

value 10% different from the optimal values reported in Table 5 provides a measure of the sensitivity of the estimate. We define a sensitivity measure S_i for the estimate of

as

(1)

(2)

where

is the optimal value of the error function (for values listed in Table 5) and

is the error with

set to a 90% or 110% of its optimal value, M is the number of reactions, and N_j is the number of experimental measures for each reaction. Sensitivity values are listed in Table 5 for each species, revealing that estimates of

for GLC⁰, NAD_red²⁻, PYR⁻, AKG²⁻, SUC²⁻, FUM²⁻ and COAS⁰ are not highly sensitive to the data.

Predicted apparent Gibbs free energies under physiology conditions

The fifth column in Table 4 reports the predicted apparent

(⁠

⁠) at physiological conditions representative of a muscle cell (26) (T = 310.15 K, I = 0.18 M, pH = 7, [Mg²⁺] = 0.8 mM, [K⁺] = 140 mM, [Na⁺] = 10 mM, [Ca²⁺] = 0.0001 mM). Differences between the physiological free energies (⁠

⁠) and the standard free energies (⁠

⁠) are also shown in Table 4. These differences come from the effects of non-standard physiological conditions—temperature (T), ionic strength (I), pH and cation bindings. Therefore, the physiological free energy

can be expressed by the summation of the standard free energy

and the non-standard contributions

⁠:

(3)

where T₀ = 298.15 K, I₀ = 0 M,

is the apparent equilibrium constant for the associated biochemical reaction,

denotes the temperature contribution,

denotes the ionic strength contribution,

denotes the pH contribution and

denotes the binding polynomial contribution. An example calculation [for reaction GLK (EC 2.7.1.1)] is described in the Appendix A.

The contributions of non-standard elements to the differences between the standard free energy for the reference chemical reaction (⁠⁠) and the physiological energy values (⁠⁠) are listed in the last four columns of Table 4. The results can be divided to three cases: (i) if the stoichiometric coefficient of H⁺ is non-zero in the reference chemical reaction (v_H ≠ 0), the differences are substantial and pH contributes significantly, accounting for at least 73% of the difference, e.g. reactions GLK, PFK and GAP; (ii) if v_H = 0 (⁠⁠) and (⁠⁠), the differences tend to be relatively small compared with case (i) and temperature and/or binding polynomials contribute most to the differences, e.g. reactions PGI, TPI and ACON; (iii) if v_H = 0 (⁠⁠) and ionic strength and/or binding polynomials contribute most to the differences, e.g. reactions FBA, NPTH and PCL.

For the reactions of the pentose phosphate pathway, only the first three reactions G6PD, PGL and PGD show substantial difference between free energy for the reference chemical reaction (⁠⁠) and for biochemical reaction (⁠⁠) under physiological conditions. These three reactions have v_H ≠ 0. They are the oxidative portion of the pentose phosphate pathway (15, 23), which produce NADP_red and are essentially ‘irreversible’ in vivo (23).

Dissociation constants uncertainty and sensitivity analysis

The pK values listed in Table 2 are taken as the average value when there are several values (≥2) available in NIST database (27). For these pK values, the average value may not represent the best choice to be used in the model, i.e. some value among those available values may be more accurate than others. For some pK values, there exists only one estimate or no direct estimates. In order to predict the impact of uncertainty of these values on the model output, an uncertainty and sensitivity analysis is performed.

The following equation is used as a measure of uncertainty in a pK value when several independent measures are available:

(4)

where pK_max and pK_min refer to the maximum and minimum value of pK, respectively. Table 6 shows the computed uncertainties for these pKs.

Table 6.

Dissociation constants uncertainty analysis for reactants with several pK values available in NIST database (27)

Reactant abbreviation
ADP	0.0539		0.0545	0.01748
ATP	0.0849		0.1051	0.1215
CIT	0.0423		0.0597	0.0480
GDP	0.0630
GTP	0.0799		0.0835
ISCIT	0.0052		0.0571
MAL	0.0148	0.0153		0.0549
Pi	0.0251	0.0463	0.1975	0.1433
PEP	0.0336
SUC	0.0133	0.0149	0.1476
R5P	0.0016

Reactant abbreviation
ADP	0.0539		0.0545	0.01748
ATP	0.0849		0.1051	0.1215
CIT	0.0423		0.0597	0.0480
GDP	0.0630
GTP	0.0799		0.0835
ISCIT	0.0052		0.0571
MAL	0.0148	0.0153		0.0549
Pi	0.0251	0.0463	0.1975	0.1433
PEP	0.0336
SUC	0.0133	0.0149	0.1476
R5P	0.0016

Table 6.

Dissociation constants uncertainty analysis for reactants with several pK values available in NIST database (27)

Reactant abbreviation
ADP	0.0539		0.0545	0.01748
ATP	0.0849		0.1051	0.1215
CIT	0.0423		0.0597	0.0480
GDP	0.0630
GTP	0.0799		0.0835
ISCIT	0.0052		0.0571
MAL	0.0148	0.0153		0.0549
Pi	0.0251	0.0463	0.1975	0.1433
PEP	0.0336
SUC	0.0133	0.0149	0.1476
R5P	0.0016

Reactant abbreviation
ADP	0.0539		0.0545	0.01748
ATP	0.0849		0.1051	0.1215
CIT	0.0423		0.0597	0.0480
GDP	0.0630
GTP	0.0799		0.0835
ISCIT	0.0052		0.0571
MAL	0.0148	0.0153		0.0549
Pi	0.0251	0.0463	0.1975	0.1433
PEP	0.0336
SUC	0.0133	0.0149	0.1476
R5P	0.0016

When only one pK value estimate is available, the uncertainty is defined as the average number

of all calculated

⁠:

(5)

According to Table 6,

is equal to 0.0609.

The sensitivities of the computed thermodynamic database due to a 10% change of pK values are calculated (28):

(6)

where E is shown in equation (1), and x_i is the value of the i^thpK, Table 7 lists calculated sensitivities >0.1. Sensitivity values for all potassium and sodium ion dissociation constants are <0.1. Calcium ion dissociation constants are not included in these calculations because [Ca²⁺] = 0 for all reactions in our raw-data database. Thus calculated sensitivities of calcium ion dissociations are equal to 0.

Table 7.

Sensitivity >0.1 in dissociation constants sensitivity analysis

Reactant abbreviation
ADP	0.4936		0.3696	0.1709
ATP	0.2151		0.2665	0.2210
CIT			0.6029
COAS	0.5344
DHAP	0.6591		0.3772
F16P	1.1839	0.3181	1.3545
G6P	0.6372		0.4524
GDP			0.1873
GTP			0.1310
ISCIT	0.1230		0.2245
NADP_ox	1.7916
NADP_red	0.7611
Pi	0.1248		0.2015
PG2	1.2172
PG3	0.3698
PEP	0.2293
R5P	0.1818
RU5P	0.1007
S7P	0.1433
X5P	0.1597

Reactant abbreviation
ADP	0.4936		0.3696	0.1709
ATP	0.2151		0.2665	0.2210
CIT			0.6029
COAS	0.5344
DHAP	0.6591		0.3772
F16P	1.1839	0.3181	1.3545
G6P	0.6372		0.4524
GDP			0.1873
GTP			0.1310
ISCIT	0.1230		0.2245
NADP_ox	1.7916
NADP_red	0.7611
Pi	0.1248		0.2015
PG2	1.2172
PG3	0.3698
PEP	0.2293
R5P	0.1818
RU5P	0.1007
S7P	0.1433
X5P	0.1597

Table 7.

Sensitivity >0.1 in dissociation constants sensitivity analysis

Reactant abbreviation
ADP	0.4936		0.3696	0.1709
ATP	0.2151		0.2665	0.2210
CIT			0.6029
COAS	0.5344
DHAP	0.6591		0.3772
F16P	1.1839	0.3181	1.3545
G6P	0.6372		0.4524
GDP			0.1873
GTP			0.1310
ISCIT	0.1230		0.2245
NADP_ox	1.7916
NADP_red	0.7611
Pi	0.1248		0.2015
PG2	1.2172
PG3	0.3698
PEP	0.2293
R5P	0.1818
RU5P	0.1007
S7P	0.1433
X5P	0.1597

Reactant abbreviation
ADP	0.4936		0.3696	0.1709
ATP	0.2151		0.2665	0.2210
CIT			0.6029
COAS	0.5344
DHAP	0.6591		0.3772
F16P	1.1839	0.3181	1.3545
G6P	0.6372		0.4524
GDP			0.1873
GTP			0.1310
ISCIT	0.1230		0.2245
NADP_ox	1.7916
NADP_red	0.7611
Pi	0.1248		0.2015
PG2	1.2172
PG3	0.3698
PEP	0.2293
R5P	0.1818
RU5P	0.1007
S7P	0.1433
X5P	0.1597

Open in new tab Download slide

The product (US) combining uncertainty U and sensitivity S can be used to check the overlapping effect of uncertainty and sensitivity. For example, recall that we arbitrarily assign the value of 4.995 to the pK_H1 for PGN. Since the value is not available in NIST database, its uncertainty U is set to the average number 0.0609. If we consider the theoretical range of 4–5.99 discussed above, then the calculated uncertainty U is ∼0.4. For this case, because the computed pK_H1 sensitivity is 0.0225, the US product is <0.01, which is small enough that the value of pK_H1 for PGN has no substantial effect on the model output. Similarly, recall that E4P, RU5P, S7P and X5P are assumed to have dissociation properties equivalent to R5P in our calculations. Although the computed pK sensitivities of E4P, RU5P, S7P and X5P are 3.81E-6, 0.1007, 0.1433 and 0.1597, respectively, all the US products of the pK_H1 for these reactants are <0.01. Therefore, the assumption of equivalent dissociation properties for E4P, R5P, RU5P, S7P and X5P does not substantially effect our calculations.

Figure 3A illustrates that US products span eight orders of magnitude. Figure 3B illustrates the detailed distribution of the US products >0.01. All US products are <0.11. There are 23 cases for which US >0.01. These 23 US values belong to 15 reactants and four pKs as listed in Table 8. They are a subset of the pKs with sensitivity >0.1. They demonstrate the most important pKs which can make obvious impact on the model output. The largest four US values are as indicated in Figure 3B.

Figure 3.

(A) Distribution of the product of uncertainty and sensitivity (US) for all pK values; (B) detailed distribution of the product >0.01.

Table 8.

The product (⁠ forumla ⁠) >0.01 in dissociation constants uncertainty and sensitivity analysis

Reactant abbreviation
ADP
ATP
CIT
COAS
DHAP	⁠,
F16P
G6P
GDP
GTP
ISCIT
Pi
NADP_ox
NADP_red
PG2
PG3

Table 8.

The product (⁠ forumla ⁠) >0.01 in dissociation constants uncertainty and sensitivity analysis

Reactant abbreviation
ADP
ATP
CIT
COAS
DHAP	⁠,
F16P
G6P
GDP
GTP
ISCIT
Pi
NADP_ox
NADP_red
PG2
PG3

Database dissemination

ThermoML is an extensible markup language (XML)-based approach, which is an IUPAC standard for storage and exchange of thermodynamic property data (29–32). Our optimized results are stored in the standard ThermoML format with two small extensions to the current ThermoML schema (32): (i) adding ‘pseudo-Gibbs free energy of formation, kJ/mol’ in the list of ePropName in BioProperties of PureOrMixtureData; and (ii) adding ‘biochemical network calculation’ in ePredictionType of Prediction. We use the term ‘pseudo-Gibbs free energy of formation’ because in our thermodynamic model the estimated values represent adjustable parameters for the given set of interdependent biochemical reactions. In our ThermoML data files, the abbreviation name of reactant is also added in sCommonName in Compound. In the ThermoML reaction database, if the value of is not available, it is specified in sPredictionMethodDescription, and both nPropValue and nPropDigits are set to 0 in PropertyValue of NumValues. The ThermoML schema and our ThermoML data files are provided in Supplementary Data.

Discussion

We have updated our biochemical thermodynamic database by adding the reactions of the pentose phosphate pathway. To build the new database, each original publication has been studied on a case-by-case basis, case-specific assumptions and approximations determined and documented, and case-specific calculations performed and documented. Raw data and documented estimations of solution properties are made electronically available so that the updated database remains transparent and extensible.

The developed database is optimally self-consistent and consistent with the data available for reactions considered and the constraints. Theoretical predictions of apparent equilibrium constants optimally match experimental data on equilibrium constants. These apparent equilibrium constants are predicted based on the estimated species-level Gibbs free energies of formation and accounting for the effects of temperature, ionic interactions and hydrogen and metal cation binding. The new database provides thermodynamic and network-based estimates of thermodynamic properties for six reactions of the pentose phosphate pathway for which estimates are not available in the pre-existing literature.

These calculations demonstrate how network thermodynamic calculations are effectively extended. Adding raw experimental data on reactions and reactants of the pentose phosphate pathway into corresponding raw-data database, reaction database and reactant database, respectively, a previous thermodynamic network model was extended to include these elements.

Although the new optimal estimates of for reactants of glycolysis and tricarboxylic acid cycle are equal to previous estimates (7), most of the optimal predicted values which are associated with the predictions are shifted compared to the previous version. This result demonstrates that reliable and self-consistent extensions require the recalculation of entire set of species-level parameters. Furthermore, sensitivity analysis reveals that some formation energies can vary substantially without changing the optimized objective value significantly. While this set of optimal values represents a self-consistent set, combining these estimated formation energies with independently estimated parameters from other studies would require re-optimizing a combined database.

The uncertainty and sensitivity analysis of dissociation constants reveals 23 pKs most important to the model output. Additional experimental measurements of these parameter values are desirable.

Supplementary Data

Supplementary data are available at Database Online.

Acknowledgements

The authors are grateful to Robert Goldberg for advice and critical comments.

Funding

National Institutes of Health Heart Lung and Blood Institute (grant number HL072011). Funding for open access charge: National Institutes of Health Heart Lung and Blood Institute.

Conflict of interest. None declared.

References

1

Vojinović

V

,

von Stockar

U

.

Influence of uncertainties in pH, pMg, activity coefficients, metabolite concentrations, and other factors on the analysis of the thermodynamic feasibility of metabolic pathways

,

Biotechnol. Bioeng.

,

2009

, vol.

103

(pg.

780

-

795

)

2

Mavrovouniotis

ML

.

Hunter

L

,

Searls

DB

,

Shavlik

JW

.

Identification of localized and distributed bottlenecks in metabolic pathways

,

Proceedings of the First International Conference on Intelligent Systems for Molecular Biology

,

1993

Bethesda, MD

AAAI Press

(pg.

275

-

283

)

3

Yang

F

,

Beard

DA

.

Thermodynamically based profiling of drug metabolism and drug-drug metabolic interactions: a case study of acetaminophen and ethanol toxic interaction

,

Biophys. Chem.

,

2006

, vol.

120

(pg.

121

-

134

)

4

Kümmel

A

,

Panke

S

,

Heinemann

M

.

Putative regulatory sites unraveled by network-embedded thermodynamic analysis of metabolome data

,

Mol. Syst. Biol.

,

2006

, vol.

2

2006.0034

5

Henry

CS

,

Broadbelt

LJ

,

Hatzimanikatis

V

.

Thermodynamics-based metabolic flux analysis

,

Biophys. J.

,

2007

, vol.

92

(pg.

1792

-

1805

)

6

Vinnakota

KC

,

Wu

F

,

Kushmerick

MJ

, et al.

Multiple ion binding equilibria, reaction kinetics, and thermodynamics in dynamic models of biochemical pathways

,

Methods Enzymol.

,

2009

, vol.

454

(pg.

29

-

68

)

7

Li

X

,

Dash

RK

,

Pradhan

RK

, et al.

A database of thermodynamic quantities for the reactions of glycolysis and the tricarboxylic acid cycle

,

J. Phys. Chem.

,

2010

, vol.

114

(pg.

16068

-

16082

)

8

Alberty

RA

. ,

Thermodynamics of Biochemical Reactions

,

2003

Hoboken, NJ

John Wiley

9

Beard

DA

,

Qian

H

.

Conventions and calculations for biochemical systems

,

Chemical Biophysics: Quantitative Analysis of Cellular Systems

,

2008

Cambridge, UK

Cambridge University Press

(pg.

24

-

40

)

10

Alberty

RA

.

Thermodynamic properties of weak acids involved in enzyme-catalyzed reactions

,

J. Phys. Chem. B

,

2006

, vol.

110

(pg.

5012

-

5016

)

11

Clarke

ECW

,

Glew

DN

.

Evaluation of Debye–Hückel limiting slopes for water between 0 and 150°C

,

J. Chem. Soc., Faraday Trans. 1

,

1980

, vol.

76

(pg.

1911

-

1916

)

12

Alberty

RA

.

Effect of temperature on standard transformed Gibbs energies of formation of reactants at specified pH and ionic strength and apparent equilibrium constants of biochemical reactions

,

J. Phys. Chem. B

,

2001

, vol.

105

(pg.

7865

-

7870

)

13

Alberty

RA

.

Calculation of thermodynamic properties of species of biochemical reactants using the inverse legendre transform

,

J. Phys. Chem. B

,

2005

, vol.

109

(pg.

9132

-

9139

)

14

Glaser

L

,

Brown

DH

.

Purification and properties of D-glucose-6-phosphate dehydrogenase

,

J. Biol. Chem.

,

1955

, vol.

216

(pg.

67

-

79

)

15

Casazza

JP

,

Veech

RL

.

The interdependence of glycolytic and pentose cycle intermediates in ad libitum fed rats

,

J. Biol. Chem.

,

1986

, vol.

261

(pg.

690

-

698

)

16

Villet

RH

,

Dalziel

K

.

The nature of the carbon dioxide substrate and equilibrium constant of the 6-phosphogluconate dehydrogenase reaction

,

Biochem. J.

,

1969

, vol.

115

(pg.

633

-

638

)

17

Axelrod

B

,

Jang

R

.

Purification and Properties of Phosphoriboisomerase from Alfalfa

,

J. Biol. Chem.

,

1954

, vol.

209

(pg.

847

-

855

)

18

Hurwitz

J

,

Horecker

BL

.

The purification of phosphoketopentoepimerase from lactobacillus pentosus and the preparation of xylulose 5-phosphate

,

J. Biol. Chem.

,

1956

, vol.

223

(pg.

993

-

1008

)

19

Datta

AG

,

Racker

E

.

Mechanism of action of transketolase

,

J. Biol. Chem.

,

1961

, vol.

236

(pg.

617

-

623

)

20

Goldberg

RN

,

Tewari

YB

,

Bell

D

, et al.

Thermodynamics of enzyme-catalyzed reactions: Part 1. Oxidoreductases

,

J. Phys. Chem. Ref. Data

,

1993

, vol.

22

(pg.

515

-

82

)

21

Goldberg

RN

,

Tewari

YB

.

Thermodynamics of enzyme-catalyzed reactions: Part 2. Transferases

,

J. Phys. Chem. Ref. Data

,

1994

, vol.

23

(pg.

547

-

617

)

22

Goldberg

RN

,

Tewari

YB

.

Thermodynamics of enzyme-catalyzed reactions: Part 5. Isomerases and ligases

,

J. Phys. Chem. Ref. Data

,

1995

, vol.

24

(pg.

1765

-

1801

)

23

Nelson

DL

,

Cox

MM

.

Glycolysis, gluconeogenesis, and the pentose phosphate pathway

,

Lehninger Principle of Biochemistry

,

2005

New York

W. H. Freeman and Company Press

24

Perrin

DD

,

Dempsey

B

,

Serjeant

EP

. ,

pKa Prediction for Organic Acids and Bases

,

1981

London and New York

Chapman and Hall

25

Goldberg

RN

,

Tewari

YB

.

Thermodynamic and transport properties of carbohydrates and their monophosphates: the pentoses and hexoses

,

J. Phys. Chem. Ref. Data

,

1989

, vol.

18

(pg.

809

-

880

)

26

Godt

RE

,

Maughan

DW

.

On the composition of the cytosol of relaxed skeletal muscle of the frog

,

Am. J. Physiol. Cell Physiol.

,

1988

, vol.

254

(pg.

C591

-

C604

)

27

Smith

RM

,

Martell

AE

,

Motekaitis

RJ

.

NIST standard reference database 46 version 8.0: NIST critically selected stability constants of metal complexes

,

NIST Standard Reference Data

,

2004

Gaithersburg, MD

28

Wu

F

,

Yang

F

,

Vinnakota

KC

, et al.

Computer modeling of mitochondrial tricarboxylic acid cycle, oxidative phosphorylation, metabolite transport, and electrophysiology

,

J. Biol. Chem.

,

2007

, vol.

282

(pg.

24525

-

24537

)

29

Frenkel

M

,

Chirico

RD

,

Diky

VV

, et al.

ThermoML - an XML-based approach for storage and exchange of experimental and critically evaluated thermophysical and thermochemical property data. 1. Experimental data

,

J. Chem. Eng. Data

,

2002

, vol.

48

(pg.

2

-

13

)

30

Chirico

RD

,

Frenkel

M

,

Diky

VV

, et al.

ThermoML - an XML-based approach for storage and exchange of experimental and critically evaluated thermophysical and thermochemical property data. 2. Uncertainties

,

J. Chem. Eng. Data

,

2003

, vol.

48

(pg.

1344

-

1359

)

31

Frenkel

M

,

Chirico

RD

,

Diky

VV

, et al.

ThermoML - an XML-based approach for storage and exchange of experimental and critically evaluated thermophysical and thermochemical property data. 3. Critically evaluated data, predicted data, and equation representation

,

J. Chem. Eng. Data

,

2004

, vol.

49

(pg.

381

-

393

)

32

Chirico

RD

,

Frenkel

M

,

Diky

V

, et al.

ThermoML - an XML-based approach for storage and exchange of experimental and critically evaluated thermophysical and thermochemical property data. 4. Biomaterials

,

J. Chem. Eng. Data

,

2009

, vol.

55

(pg.

1564

-

1572

)

33

Goldberg

RN

,

Tewari

YB

.

Thermodynamics of enzyme-catalyzed reactions. Part 3. Hydrolases

,

J. Phys. Chem. Ref. Data

,

1994

, vol.

23

(pg.

1035

-

1103

)

34

Alberty

RA

.

Standard transformed formation properties of carbon dioxide in aqueous solutions at specified pH

,

J. Phys. Chem.

,

1995

, vol.

99

(pg.

11028

-

11034

)

35

Tewari

YB

,

Steckler

DK

,

Goldberg

RN

, et al.

Thermodynamics of hydrolysis of sugar phosphates

,

J. Biol. Chem.

,

1988

, vol.

263

(pg.

3670

-

3675

)

36

Martell

AE

,

Smith

RM

,

Motekaitis

RJ

.

NIST standard reference database 46 version 8.0: NIST critically selected stability constants of metal complexes

,

NIST Standard Reference Data

,

2004

Gaithersburg, MD

37

Larsson-Raźnikiewicz

M

.

Complex formation between 3-phospho-D-glyceric acid and divalent metal ions

,

European Journal of Biochemistry

,

1972

, vol.

30

(pg.

579

-

583

)

38

Merrill

DK

,

McAlexander

JC

,

Guynn

RW

.

Equilibrium constants under physiological conditions for the reactions of -3-phosphoglycerate dehydrogenase and -phosphoserine aminotransferase

,

Arch. Biochem. Biophys.

,

1981

, vol.

212

(pg.

717

-

729

)

39

Briggs

TN

,

Stuehr

JE

.

Simultaneous potentiometric determination of precise equivalence points and pK values of two- and three-pK systems

,

Anal. Chem.

,

1975

, vol.

47

(pg.

1916

-

1920

)