modelHemoglobinCarboxylation

Multiple-ligand allosteric hemoglobin model
Diagram of HemoglobinCarboxylation

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

TypeNameDefaultDescription
Modelica.Units.SI.AmountOfSubstanceTHb0.001Total amount of hemoglobin
RealRTModelica.Constants.R*298.15
Modelica.Units.SI.VolumeOneLiter0.001
RealL_old7.0529*10^6=[T0]/[R0] .. dissociation constant of relaxed <-> tensed change of deoxyhemoglobin tetramer
Realc0.00431555=KR/KT .. ration between oxygen affinities of relaxed vs. tensed subunit
Modelica.Units.SI.ConcentrationKR0.000671946*(55.508/38.7)oxygen dissociation on relaxed(R) hemoglobin subunit
Modelica.Units.SI.ConcentrationKTKR/coxygen dissociation on tensed(T) hemoglobin subunit
Modelica.Units.SI.MoleFractionKRo37KR*OneLiter
Modelica.Units.SI.MoleFractionKTo37KT*OneLiter
Modelica.Units.SI.ChemicalPotentialDfG_O2-RT*log(0.0013)
Modelica.Units.SI.ChemicalPotentialDfG_CO2-RT*log(0.034) - 394400
Modelica.Units.SI.ChemicalPotentialDfG_tT0
Modelica.Units.SI.ChemicalPotentialDfG_tRDfG_tT + RT*log(L)
RealKC1e-8Slow down factor
Modelica.Units.SI.MoleFractioninitialO20.000132233Initial O2 at 37degC, pO2=100Pa
Modelica.Units.SI.MoleFractioninitialH10^(-6.9)
Modelica.Units.SI.MoleFractioninitialCO20.00136212Initial CO2 at 37degC, pCO2=40mmHg
Modelica.Units.SI.MoleFractioninitialCO1e-11Initial CO at 37degC, pCO=0mmHg
Modelica.Units.SI.MoleFractionKRh3710^(-6.89)
Modelica.Units.SI.MoleFractionKTh3710^(-7.52)
Modelica.Units.SI.MoleFractionKRz3710^(-7.25)
Modelica.Units.SI.MoleFractionKTz3710^(-7.73)
Modelica.Units.SI.MoleFractionKRc37(10^(-8.35))/(OneLiter)
Modelica.Units.SI.MoleFractionKTc37(10^(-7.54))/(OneLiter)
RealLL_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

TypeNameDefaultDescription
Chemical.Obsolete.Components.Solutionsolution
Chemical.Obsolete.Components.ReactionquaternaryForm
Chemical.Obsolete.Components.SubstanceO2_free
Chemical.Obsolete.Sources.ExternalIdealGasSubstanceoxygen_in_air
Chemical.Obsolete.Components.GasSolubilitypartialPressure1
Obsolete.Components.SubstanceH2O
HemoglobinQuaternaryFormrelaxed
HemoglobinQuaternaryFormtensed
Chemical.Obsolete.Sources.ExternalMoleFractionH
Chemical.Obsolete.Components.SubstanceCO2_free
Chemical.Obsolete.Sources.ExternalIdealGasSubstanceCO2_gas
Chemical.Obsolete.Components.GasSolubilitypartialPressure2
RealsO2Hemoglobin oxygen saturation
RealsCO2Hemoglobin carbon dioxide saturation
RealdHHemoglobin charge change caused by binding of Bohr's protons
Modelica.Blocks.Sources.ContinuousClockcarbonDioxideSource

Revisions

2013-2018

Marek Matejak, Charles University, Prague, Czech Republic