modelAllosteric_Hemoglobin_MWC

Monod,Wyman,Changeux (1965)

Extends from Modelica.Icons.Example (Icon for runnable examples).

Information

To understand the model is necessary to study the principles of MWC allosteric transitions first published by

[1] Monod,Wyman,Changeux (1965). "On the nature of allosteric transitions: a plausible model." Journal of molecular biology 12(1): 88-118.


In short it is about binding oxygen to hemoglobin.

Oxgen are driven by its partial pressure using clock source - from very little pressure to pressure of 10kPa.

(Partial pressure of oxygen in air is the air pressure multiplied by the fraction of the oxygen in air.)

Hemoglobin was observed (by Perutz) in two structuraly different forms R and T.

These forms are represented by blocks T0..T4 and R0..R4, where the suffexed index means the number of oxygen bounded to the form.


In equilibrated model can be four chemical reactions removed and the results will be the same, but dynamics will change a lot. ;)

If you remove the quaternaryForm1,quaternaryForm2,quaternaryForm3,quaternaryForm4 then the model in equilibrium will be exactly the same as in MWC article.


Parameters was fitted to data of Severinghaus article from 1979. (For example at pO2=26mmHg is oxygen saturation sO2 = 48.27 %).

Parameters

TypeNameDefaultDescription
Types.TemperatureT310.15Temperature
Types.FractionL7.0529*10^6=[T0]/[R0] .. dissociation constant of relaxed <-> tensed change of deoxyhemoglobin tetramer
Types.Fractionc0.00431555=KR/KT .. ration between oxygen affinities of relaxed vs. tensed subunit
Types.ConcentrationKR0.000671946oxygen dissociation on relaxed(R) hemoglobin subunit
Types.ConcentrationKTKR/coxygen dissociation on tensed(T) hemoglobin subunit
Types.AmountOfSubstancetotalAmountOfHemoglobin1

Components

TypeNameDefaultDescription
Components.SubstanceT0
Components.SubstanceT1
Components.SubstanceR1
Components.SubstanceT2
Components.SubstanceR2
Components.SubstanceT3
Components.SubstanceR3
Components.SubstanceT4
Components.SubstanceR4
Components.SubstanceR0
Components.ChemicalReactionquaternaryForm
Components.ChemicalReactionoxyR1
Components.ChemicalReactionoxyT1
Components.ChemicalReactionoxyR2
Components.ChemicalReactionoxyR3
Components.ChemicalReactionoxyR4
Components.ChemicalReactionoxyT2
Components.ChemicalReactionoxyT3
Components.ChemicalReactionoxyT4
Components.ChemicalReactionquaternaryForm1
Components.ChemicalReactionquaternaryForm2
Components.ChemicalReactionquaternaryForm3
Components.ChemicalReactionquaternaryForm4
SteadyStates.Components.MolarConservationLawhemoglobinConservationLaw
Modelica.Blocks.Math.Sumoxygen_bound
Modelica.Blocks.Math.DivisionsO2_hemoglobin oxygen saturation
Modelica.Blocks.Math.SumtHb
Components.Substanceoxygen_unbound
Modelica.Blocks.Sources.ContinuousClockclock
Sources.UnlimitedGasStorageO2_in_air
Components.GasSolubilitygasSolubility

Revisions

2013

Marek Matejak, Charles University, Prague, Czech Republic