mulquiney2

vRBCvAK

vRBCvAK

adpRBC + mgadpRBC = ampRBC + mgatpRBC

vRBCvALD

vRBCvALD

f16p2RBC = dhapRBC + gapRBC

vRBCvATPASE

vRBCvATPASE

mgatpRBC = phosRBC + mgadpRBC

vRBCvBPGSP1

vRBCvBPGSP1

bpgspRBC + b13pgRBC = bpgspb13pgRBC

vRBCvBPGSP2

vRBCvBPGSP2

bpgspb13pgRBC = p3gRBC + bpgsppRBC

vRBCvBPGSP3

vRBCvBPGSP3

p3gRBC + bpgsppRBC = bpgsppp3gRBC

vRBCvBPGSP4

vRBCvBPGSP4

p2gRBC + bpgsppRBC = bpgsppp2gRBC

vRBCvBPGSP5

vRBCvBPGSP5

bpgsppp3gRBC = bpgspb23pgRBC

vRBCvBPGSP6

vRBCvBPGSP6

bpgsppp2gRBC = bpgspb23pgRBC

vRBCvBPGSP7

vRBCvBPGSP7

bpgspb23pgRBC = b23pgRBC + bpgspRBC

vRBCvBPGSP8

vRBCvBPGSP8

phosRBC + bpgsppRBC = bpgsppphosRBC

vRBCvBPGSP9

vRBCvBPGSP9

bpgsppphosRBC = {2.0}phosRBC + bpgspRBC

vRBCvENO

vRBCvENO

p2gRBC = pepRBC

vRBCvG6PDH

vRBCvG6PDH

nadpRBC + g6pRBC = p6glRBC + nadphRBC

vRBCvGAPDH

vRBCvGAPDH

nadRBC + phosRBC + gapRBC = b13pgRBC + nadhRBC

vRBCvGLCTRANSPORT

vRBCvGLCTRANSPORT

glcEXT = glcRBC

vRBCvGSSGR

vRBCvGSSGR

gssgRBC + nadphRBC = nadpRBC + {2.0}gshRBC

vRBCvHBADP

vRBCvHBADP

adpRBC + hbRBC = hbadpRBC

vRBCvHBATP

vRBCvHBATP

atpRBC + hbRBC = hbatpRBC

vRBCvHBB13PG

vRBCvHBB13PG

b13pgRBC + hbRBC = hbb13pgRBC

vRBCvHBB23PG

vRBCvHBB23PG

hbRBC + b23pgRBC = hbb23pgRBC

vRBCvHBMGATP

vRBCvHBMGATP

hbRBC + mgatpRBC = hbmgatpRBC

vRBCvHK

vRBCvHK

glcRBC + mgatpRBC = mgadpRBC + g6pRBC

vRBCvLACTRANSPORT

vRBCvLACTRANSPORT

lacRBC = lacEXT

vRBCvLDH

vRBCvLDH

pyrRBC + nadhRBC = nadRBC + lacRBC

vRBCvLDHP

vRBCvLDHP

pyrRBC + nadphRBC = nadpRBC + lacRBC

vRBCvMGADP

vRBCvMGADP

adpRBC + mgRBC = mgadpRBC

vRBCvMGATP

vRBCvMGATP

atpRBC + mgRBC = mgatpRBC

vRBCvMGB13PG

vRBCvMGB13PG

b13pgRBC + mgRBC = mgb13pgRBC

vRBCvMGB23PG

vRBCvMGB23PG

b23pgRBC + mgRBC = mgb23pgRBC

vRBCvMGF16P2

vRBCvMGF16P2

f16p2RBC + mgRBC = mgf16p2RBC

vRBCvMGG16P2

vRBCvMGG16P2

g16p2RBC + mgRBC = mgg16p2RBC

vRBCvMGPHOS

vRBCvMGPHOS

phosRBC + mgRBC = mgphosRBC

vRBCvOX

vRBCvOX

{2.0}gshRBC = gssgRBC

vRBCvOXNADH

vRBCvOXNADH

nadhRBC = nadRBC

vRBCvP6GDH

vRBCvP6GDH

nadpRBC + p6gRBC = ru5pRBC + co2RBC + nadphRBC

vRBCvPFK

vRBCvPFK

f6pRBC + mgatpRBC = mgadpRBC + f16p2RBC

vRBCvPGI

vRBCvPGI

g6pRBC = f6pRBC

vRBCvPGK

vRBCvPGK

mgadpRBC + b13pgRBC = p3gRBC + mgatpRBC

vRBCvPGLHYDROLYSIS

vRBCvPGLHYDROLYSIS

p6glRBC = p6gRBC

vRBCvPGM

vRBCvPGM

p3gRBC = p2gRBC

vRBCvPHOSTRANSPORT

vRBCvPHOSTRANSPORT

phosRBC = phosEXT

vRBCvPK

vRBCvPK

mgadpRBC + pepRBC = mgatpRBC + pyrRBC

vRBCvPYRTRANSPORT

vRBCvPYRTRANSPORT

pyrRBC = pyrEXT

vRBCvR5PI

vRBCvR5PI

ru5pRBC = rib5pRBC

vRBCvRu5PE

vRBCvRu5PE

ru5pRBC = xu5pRBC

vRBCvTA

vRBCvTA

gapRBC + sed7pRBC = f6pRBC + ery4pRBC

vRBCvTIM

vRBCvTIM

gapRBC = dhapRBC

vRBCvTK1

vRBCvTK1

tkRBC + xu5pRBC = tkxu5pRBC

vRBCvTK2

vRBCvTK2

tkxu5pRBC = gapRBC + tkgRBC

vRBCvTK3

vRBCvTK3

rib5pRBC + tkgRBC = tkgrib5pRBC

vRBCvTK4

vRBCvTK4

tkgrib5pRBC = sed7pRBC + tkRBC

vRBCvTK5

vRBCvTK5

ery4pRBC + tkgRBC = tkgery4pRBC

vRBCvTK6

vRBCvTK6

tkgery4pRBC = f6pRBC + tkRBC

Global parameters

Assignment rules

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)

Function definitions

Note that constraints are not enforced in simulations. It remains the responsibility of the user to verify that simulation results satisfy these constraints.


Species:

Reactions:


Middle-click: pin/unpin nodes
Shift-click: pool/unpool species
Right-click: context menu

Apply alternate model layout to overlapping elements in current model:

log scales

y-axis min/max

x-axis min/max

Abstract
Over the last 25 years, several mathematical models of erythrocyte metabolism have been developed. Although these models have identified the key features in the regulation and control of erythrocyte metabolism, many important aspects remain unexplained. In particular, none of these models have satisfactorily accounted for 2,3-bisphosphoglycerate (2,3-BPG) metabolism. 2,3-BPG is an important modulator of haemoglobin oxygen affinity, and hence an understanding of the regulation of 2,3-BPG concentration is important for understanding blood oxygen transport. A detailed, comprehensive, and hence realistic mathematical model of erythrocyte metabolism is presented that can explain the regulation and control of 2,3-BPG concentration and turnover. The model is restricted to the core metabolic pathways, namely glycolysis, the 2,3-BPG shunt and the pentose phosphate pathway (PPP), and includes membrane transport of metabolites, the binding of metabolites to haemoglobin and Mg(2+), as well as pH effects on key enzymic reactions and binding processes. The model is necessarily complex, since it is intended to describe the regulation and control of 2,3-BPG metabolism under a wide variety of physiological and experimental conditions. In addition, since H(+) and blood oxygen tension are important external effectors of 2,3-BPG concentration, it was important that the model take into account the large array of kinetic and binding phenomena that result from changes in these effectors. Through an iterative loop of experimental and simulation analysis many values of enzyme-kinetic parameters of the model were refined to yield close conformity between model simulations and 'real' experimental data. This iterative process enabled a single set of parameters to be found which described well the metabolic behaviour of the erythrocyte under a wide variety of conditions.
A copy of the original mulquiney model but in units of fmol/min and a glucose transporter added from Philip W. Potts, Jennifer R., Kuchel. Anomeric preference of fluoroglucose exchange across human red-cell membranes. 19F-nmr studies. Biochemical Journal, 281(3):753–759, 1992.