modelHemoglobinMultipleAllosteryCO
Extends from Modelica.Icons.Example (Icon for runnable examples).
Information
Oxygen dissociation curve of hemoglobin.
M. Mateják, T. Kulhánek, and S. Matoušek, "Adair-based hemoglobin equilibrium with oxygen, carbon dioxide and hydrogen ion activity," Scandinavian Journal of Clinical & Laboratory Investigation, pp. 1-8, 2015.
J. W. Severinghaus, "Simple, accurate equations for human blood O2 dissociation computations," Journal of Applied Physiology, vol. 46, pp. 599-602, 1979.
pO2 .. partial pressure of oxygen in gas
pCO2 .. partial pressure of carbon dioxide
sO2 .. oxygen saturation of hemoglobin
pH = log10(aH), where aH is mole fraction based activity of hydrogen ions
R. B. Reeves, "The effect of temperature on the oxygen equilibrium curve of human blood," Respiration physiology, vol. 42, pp. 317-328, 1980.
T .. temperature
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| Modelica.Units.SI.AmountOfSubstance | THb | 0.001 | Total amount of hemoglobin |
| Real | RT | Modelica.Constants.R*298.15 | |
| Modelica.Units.SI.Volume | OneLiter | 0.001 | |
| Real | L_old | 7.0529*10^6 | =[T0]/[R0] .. dissociation constant of relaxed <-> tensed change of deoxyhemoglobin tetramer |
| Real | c | 0.00431555 | =KR/KT .. ration between oxygen affinities of relaxed vs. tensed subunit |
| Modelica.Units.SI.Concentration | KR | 0.000671946*(55.508/38.7) | oxygen dissociation on relaxed(R) hemoglobin subunit |
| Modelica.Units.SI.Concentration | KT | KR/c | oxygen dissociation on tensed(T) hemoglobin subunit |
| Modelica.Units.SI.MoleFraction | KRo37 | KR*OneLiter | |
| Modelica.Units.SI.MoleFraction | KTo37 | KT*OneLiter | |
| Modelica.Units.SI.ChemicalPotential | DfG_O2 | -RT*log(0.0013) | |
| Modelica.Units.SI.ChemicalPotential | DfG_CO2 | -RT*log(0.034) - 394400 | |
| Modelica.Units.SI.ChemicalPotential | DfG_tT | 0 | |
| Modelica.Units.SI.ChemicalPotential | DfG_tR | DfG_tT + RT*log(L) | |
| Real | KC | 1e-8 | Slow down factor |
| Modelica.Units.SI.MoleFraction | initialO2 | 1.9594e-07 | Initial O2 at 37degC, pO2=100Pa |
| Modelica.Units.SI.MoleFraction | initialH | 10^(-7.2) | |
| Modelica.Units.SI.MoleFraction | initialCO2 | 2.4217e-10 | Initial CO2 at 37degC, pCO2=40mmHg |
| Modelica.Units.SI.MoleFraction | initialCO | 1e-12 | Initial CO at 37degC, pCO=0mmHg |
| Modelica.Units.SI.MoleFraction | KRh37 | 10^(-6.89) | |
| Modelica.Units.SI.MoleFraction | KTh37 | 10^(-7.52) | |
| Modelica.Units.SI.MoleFraction | KRz37 | 10^(-7.25) | |
| Modelica.Units.SI.MoleFraction | KTz37 | 10^(-7.73) | |
| Modelica.Units.SI.MoleFraction | KRc37 | (10^(-8.35))/(OneLiter) | |
| Modelica.Units.SI.MoleFraction | KTc37 | (10^(-7.54))/(OneLiter) | |
| Real | L | L_old*(((KTh37/((10^(-7.2)) + KTh37))/(KRh37/((10^(-7.2)) + KRh37)))^4)*(((KTz37*((10^(-7.2))^2 + KRz37*(10^(-7.2)) + KRz37*KRc37*(2.4217e-5)))/(KRz37*((10^(-7.2))^2 + KTz37*(10^(-7.2)) + KTz37*KTc37*(2.4217e-5))))^4) | =[T0]/[R0] .. dissociation constant of relaxed <-> tensed change of deoxyhemoglobin tetramer |
Components
Revisions
2013-2018
Marek Matejak, Charles University, Prague, Czech Republic