vRBCvAK
adpRBC + mgadpRBC = ampRBC + mgatpRBC
vRBCvALD
f16p2RBC = dhapRBC + gapRBC
vRBCvATPASE
mgatpRBC = phosRBC + mgadpRBC
vRBCvBPGSP1
bpgspRBC + b13pgRBC = bpgspb13pgRBC
vRBCvBPGSP2
bpgspb13pgRBC = p3gRBC + bpgsppRBC
vRBCvBPGSP3
p3gRBC + bpgsppRBC = bpgsppp3gRBC
vRBCvBPGSP4
p2gRBC + bpgsppRBC = bpgsppp2gRBC
vRBCvBPGSP5
bpgsppp3gRBC = bpgspb23pgRBC
vRBCvBPGSP6
bpgsppp2gRBC = bpgspb23pgRBC
vRBCvBPGSP7
bpgspb23pgRBC = b23pgRBC + bpgspRBC
vRBCvBPGSP8
phosRBC + bpgsppRBC = bpgsppphosRBC
vRBCvBPGSP9
bpgsppphosRBC = {2.0}phosRBC + bpgspRBC
vRBCvENO
p2gRBC = pepRBC
vRBCvG6PDH
nadpRBC + g6pRBC = p6glRBC + nadphRBC
vRBCvGAPDH
nadRBC + phosRBC + gapRBC = b13pgRBC + nadhRBC
vRBCvGLCTRANSPORT
glcEXT = glcRBC
vRBCvGSSGR
gssgRBC + nadphRBC = nadpRBC + {2.0}gshRBC
vRBCvHBADP
adpRBC + hbRBC = hbadpRBC
vRBCvHBATP
atpRBC + hbRBC = hbatpRBC
vRBCvHBB13PG
b13pgRBC + hbRBC = hbb13pgRBC
vRBCvHBB23PG
hbRBC + b23pgRBC = hbb23pgRBC
vRBCvHBMGATP
hbRBC + mgatpRBC = hbmgatpRBC
vRBCvHK
glcRBC + mgatpRBC = mgadpRBC + g6pRBC
vRBCvLACTRANSPORT
lacRBC = lacEXT
vRBCvLDH
pyrRBC + nadhRBC = nadRBC + lacRBC
vRBCvLDHP
pyrRBC + nadphRBC = nadpRBC + lacRBC
vRBCvMGADP
adpRBC + mgRBC = mgadpRBC
vRBCvMGATP
atpRBC + mgRBC = mgatpRBC
vRBCvMGB13PG
b13pgRBC + mgRBC = mgb13pgRBC
vRBCvMGB23PG
b23pgRBC + mgRBC = mgb23pgRBC
vRBCvMGF16P2
f16p2RBC + mgRBC = mgf16p2RBC
vRBCvMGG16P2
g16p2RBC + mgRBC = mgg16p2RBC
vRBCvMGPHOS
phosRBC + mgRBC = mgphosRBC
vRBCvOX
{2.0}gshRBC = gssgRBC
vRBCvOXNADH
nadhRBC = nadRBC
vRBCvP6GDH
nadpRBC + p6gRBC = ru5pRBC + co2RBC + nadphRBC
vRBCvPFK
f6pRBC + mgatpRBC = mgadpRBC + f16p2RBC
vRBCvPGI
g6pRBC = f6pRBC
vRBCvPGK
mgadpRBC + b13pgRBC = p3gRBC + mgatpRBC
vRBCvPGLHYDROLYSIS
p6glRBC = p6gRBC
vRBCvPGM
p3gRBC = p2gRBC
vRBCvPHOSTRANSPORT
phosRBC = phosEXT
vRBCvPK
mgadpRBC + pepRBC = mgatpRBC + pyrRBC
vRBCvPYRTRANSPORT
pyrRBC = pyrEXT
vRBCvR5PI
ru5pRBC = rib5pRBC
vRBCvRu5PE
ru5pRBC = xu5pRBC
vRBCvTA
gapRBC + sed7pRBC = f6pRBC + ery4pRBC
vRBCvTIM
gapRBC = dhapRBC
vRBCvTK1
tkRBC + xu5pRBC = tkxu5pRBC
vRBCvTK2
tkxu5pRBC = gapRBC + tkgRBC
vRBCvTK3
rib5pRBC + tkgRBC = tkgrib5pRBC
vRBCvTK4
tkgrib5pRBC = sed7pRBC + tkRBC
vRBCvTK5
ery4pRBC + tkgRBC = tkgery4pRBC
vRBCvTK6
tkgery4pRBC = f6pRBC + tkRBC
pyrEXT = ConcPyrEXT * vBld
KcatrappRBCvHK = KcatrRBCvHK / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.02) + pow(10.0, -9.55) / pow(10.0, -phRBC))
KmapppyrRBCvLDH = KmpyrRBCvLDH * ((1.0 + pow(10.0, -6.8) / pow(10.0, -phRBC)) / (1.0 + pow(10.0, -6.8) / pow(10.0, -7.2)))
KmappnadhRBCvGAPDH = KmnadhRBCvGAPDH * (pow(10.0, -7.2) / pow(10.0, -phRBC))
KeqRBCvLACTRANSPORT = (1.0 + pow(10.0, phRBC - 3.73)) / (1.0 + pow(10.0, phRBC - 3.73) / RtvRBC)
KcatrappRBCvGAPDH = KcatrRBCvGAPDH / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.5) + pow(10.0, -10.0) / pow(10.0, -phRBC))
co2RBC = ConcCo2RBC * Vrbc
KiapppyrRBCvLDH = KipyrRBCvLDH * ((1.0 + pow(10.0, -6.8) / pow(10.0, -phRBC)) / (1.0 + pow(10.0, -6.8) / pow(10.0, -7.2)))
KoRBCvPHOSTRANSPORT = KiRBCvPHOSTRANSPORT / KeqRBCvPHOSTRANSPORT
KaappRBCvHBMGATP = KaRBCvHBMGATP * HbpHRBC
K2appRBCvAK = K2RBCvAK * (1.0 + hRBC * KhampRBCvAK + kRBC * KkampRBCvAK)
KiappnadhRBCvGAPDH = KinadhRBCvGAPDH * (pow(10.0, -7.2) / pow(10.0, -phRBC))
KaappRBCvHBBPG = KaRBCvHBBPG * HbpHRBC
K6appRBCvBPGSP4 = K6RBCvBPGSP4 * ((1.0 + pow(pow(10.0, -7.2) / pow(10.0, -7.17), 4.0)) / (1.0 + pow(pow(10.0, -phRBC) / pow(10.0, -7.17), 4.0)))
KaappRBCvMGG16P2 = KaRBCvMGG16P2 * CRBCvMGG16P2 * (KRBCvMGG16P2 + KhfRBCvMGG16P2 * KmghfRBCvMGG16P2 * hRBC) / (1.0 + KhfRBCvMGG16P2 * hRBC + Kh2fRBCvMGG16P2 * KhfRBCvMGG16P2 * hRBC * hRBC + kRBC * KkfRBCvMGG16P2 + kRBC * KhfRBCvMGG16P2 * KkhfRBCvMGG16P2 * hRBC)
KaappRBCvMGATP = KaRBCvMGATP * CRBCvMGATP * ((KRBCvMGATP + KhatpRBCvMGATP * KmghatpRBCvMGATP * hRBC) / (1.0 + KhatpRBCvMGATP / pHConversionFactor / pow(10.0, phRBC) + kRBC * KkatpRBCvMGATP))
K13appRBCvBPGSP7 = K13RBCvBPGSP7 * ((1.0 + pow(pow(10.0, -7.2) / pow(10.0, -7.17), 4.0)) / (1.0 + pow(pow(10.0, -phRBC) / pow(10.0, -7.17), 4.0)))
KiappgapRBCvGAPDH = KigapRBCvGAPDH / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.5) + pow(10.0, -10.0) / pow(10.0, -phRBC))
lacEXT = ConcLacEXT * vBld
phosEXT = ConcPhosEXT * vBld
KoRBCvLACTRANSPORT = KiRBCvLACTRANSPORT / KeqRBCvLACTRANSPORT
KaappRBCvHBADP = KaRBCvHBADP * HbpHRBC
KaappRBCvHBATP = KaRBCvHBATP * HbpHRBC
hRBC = pHConversionFactor / pow(10.0, phRBC)
KoRBCvPYRTRANSPORT = KiRBCvPYRTRANSPORT / RtvRBC
K4appRBCvBPGSP3 = K4RBCvBPGSP3 * ((1.0 + pow(pow(10.0, -7.2) / pow(10.0, -7.17), 4.0)) / (1.0 + pow(pow(10.0, -phRBC) / pow(10.0, -7.17), 4.0)))
KiapplacRBCvLDH = KilacRBCvLDH * ((1.0 + pow(10.0, -phRBC) / pow(10.0, -6.8)) / (1.0 + pow(10.0, -7.2) / pow(10.0, -6.8)))
KaappRBCvHBB13PG = KaRBCvHBB13PG * HbpHRBC
KeqRBCvPHOSTRANSPORT = (1.0 + pow(10.0, phRBC - 6.75)) / (1.0 / RtvRBC + pow(10.0, phRBC - 6.75) / (RtvRBC * RtvRBC))
testsum = adpRBC + atpRBC + ampRBC + mgRBC + mgadpRBC + mgatpRBC + mgb13pgRBC + mgb23pgRBC + mgf16p2RBC + mgg16p2RBC + mgphosRBC
Kiappb13pgRBCvGAPDH = Kib13pgRBCvGAPDH / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.5) + pow(10.0, -10.0) / pow(10.0, -phRBC))
LRBCvPK = pHConversionFactor * pow(10.0, -6.8) / hRBC * pow(1.0 + atpRBC / Vrbc / KtatpRBCvPK, 4.0) / (pow(1.0 + f16p2RBC / Vrbc / Krf16p2RBCvPK + g16p2RBC / Vrbc / Krg16p2RBCvPK, 4.0) * pow(1.0 + pepRBC / Vrbc / KrpepRBCvPK + pyrRBC / Vrbc / KrpyrRBCvPK, 4.0))
KaappRBCvMGPHOS = KaRBCvMGPHOS * (1.0 + pHConversionFactor * pow(10.0, -7.2) * KhphosRBCvMGPHOS + kRBC * KkphosRBCvMGPHOS) / (1.0 + hRBC * KhphosRBCvMGPHOS + kRBC * KkphosRBCvMGPHOS)
KaappRBCvMGADP = KaRBCvMGADP * CRBCvMGADP * ((KRBCvMGADP + KhadpRBCvMGADP * KmghadpRBCvMGADP * hRBC) / (1.0 + KhadpRBCvMGADP * hRBC + kRBC * KkadpRBCvMGADP))
KmapplacRBCvLDH = KmlacRBCvLDH * ((1.0 + pow(10.0, -phRBC) / pow(10.0, -6.8)) / (1.0 + pow(10.0, -7.2) / pow(10.0, -6.8)))
KcatfappRBCvHK = KcatfRBCvHK / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.02) + pow(10.0, -9.55) / pow(10.0, -phRBC))
KcatfappRBCvGAPDH = KcatfRBCvGAPDH / (1.0 + pow(10.0, -phRBC) / pow(10.0, -7.5) + pow(10.0, -10.0) / pow(10.0, -phRBC))
atpRBCsum = atpRBC + adpRBC + ampRBC
K3appRBCvBPGSP2 = K3RBCvBPGSP2 * ((1.0 + pow(pow(10.0, -7.2) / pow(10.0, -7.17), 4.0)) / (1.0 + pow(pow(10.0, -phRBC) / pow(10.0, -7.17), 4.0)))
K1appRBCvBPGSP1 = K1RBCvBPGSP1 * ((1.0 + pow(10.0, -6.8) / pow(10.0, -7.2)) / (1.0 + pow(10.0, -6.8) / pow(10.0, -phRBC)))
KaappRBCvMGF16P2 = KaRBCvMGF16P2 * CRBCvMGF16P2 * (KRBCvMGF16P2 + KhfRBCvMGF16P2 * KmghfRBCvMGF16P2 * hRBC) / (1.0 + KhfRBCvMGF16P2 * hRBC + Kh2fRBCvMGF16P2 * KhfRBCvMGF16P2 * hRBC * hRBC + kRBC * KkfRBCvMGF16P2 + kRBC * KhfRBCvMGF16P2 * KkhfRBCvMGF16P2 * hRBC)
K1appRBCvAK = K1RBCvAK * (1.0 + pHConversionFactor * pow(10.0, -phRBC) * KhadpRBCvAK + kRBC * KkadpRBCvAK)
KaappRBCvMGB13PG = KaRBCvMGB13PG * CRBCvMGB13PG * (KRBCvMGB13PG + KhbpgRBCvMGB13PG * KmghbpgRBCvMGB13PG * hRBC) / (1.0 + KhbpgRBCvMGB13PG * hRBC + Kh2bpgRBCvMGB13PG * KhbpgRBCvMGB13PG * hRBC * hRBC + kRBC * KkbpgRBCvMGB13PG + kRBC * KhbpgRBCvMGB13PG * KkhbpgRBCvMGB13PG * hRBC)
HbpHRBC = (1.0 + 2.0 * KahbRBC / (pHConversionFactor * pow(10.0, -7.2)) + pow(KahbRBC / (pHConversionFactor * pow(10.0, -7.2)), 2.0)) / (1.0 + 2.0 * KahbRBC / hRBC + KahbRBC / hRBC * KahbRBC / hRBC)
glcEXT = ConcGlcEXT*vBld
KaappRBCvMGB23PG = KaRBCvMGB23PG * CRBCvMGB23PG * (KRBCvMGB23PG + KhbpgRBCvMGB23PG * KmghbpgRBCvMGB23PG * hRBC) / (1.0 + KhbpgRBCvMGB23PG * hRBC + Kh2bpgRBCvMGB23PG * KhbpgRBCvMGB23PG * hRBC * hRBC + kRBC * KkbpgRBCvMGB23PG + kRBC * KhbpgRBCvMGB23PG * KkhbpgRBCvMGB23PG * hRBC)
Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.