Model number

Model for Oxygen binding to myoglobin, single site first order binding.


 Binding of a ligand is a function of the affinity of 
the site for the ligand. The disociation constant is ratio of the off-rate, k_1, to 
the on-rate, k1, is Kd = k_1/k1. 
The binding of oxygen (O2) to myoglobin (Mb) is given by the reaction:

          O2 + Mb <--- MbO2 , where one O2 molecule binds to the site on Mb.

To translate from physiological measurement of the partial pressure of O2 in units of 
mmHg to molar concentration, the solubility of O2 in water for the particular temperature 
has to be used (alphO2). The PO2 at 50% saturation, i.e. 50% of Mb bound to O2, is 
called p50Mb, and is the PO2 at which the molar concentration of O2 equals Kd.
    In this model, the PO2 is raised gradually in order to run it through a 
wide range of concentrations from well below Kd to well above Kd. If the rate
of binding is slow, the binding will not be at equilibrium, and saturation 
will be below the equilibrium level when PO2 is rising fast.
Change the "rateO2" at which the pO2 is raised: at low on-rates (k1) the SMb 
(the myoglobin saturation) will lag behind the curve for instantaneous equilibration, SMbeq.
 SMbeq = 1/(1+ Kd / pO2),  where Kd is the dissociation constant.


The equations for this model may be viewed by running the JSim model applet and clicking on the Source tab at the bottom left of JSim's Run Time graphical user interface. The equations are written in JSim's Mathematical Modeling Language (MML). See the Introduction to MML and the MML Reference Manual. Additional documentation for MML can be found by using the search option at the Physiome home page.

Download JSim model project file

Help running a JSim model.

Schenkman KA, Marble DR, Burns DH, and Feigl EO: Myoglobin oxygen 
 dissociation by multiwavelength spectroscopy. J Appl Physiol 82: 86-92, 1997.

Schenkman KA, Beard DA, Ciesielski WA, and Feigl EO: Comparison of buffer 
 and red blood cell perfusion of guinea pig heart oxygenation. 
 Am J Physiol Heart Circ Physiol 285: H1819-H1825, 2003.

Schenkman KA: Cardiac performance as a function of intracellular oxygen 
 tension in buffer-perfused hearts. Am J Physiol Heart Circ Physiol 
 281: H2463-H2472, 2001.
Key terms
oxygen binding
dissociation curves
saturation curves
equilibrium dissociation
binding rates
Hill coefficient
apparent cooperativity
gas transport

Please cite in any publication for which this software is used and send one reprint to the address given below:
The National Simulation Resource, Director J. B. Bassingthwaighte, Department of Bioengineering, University of Washington, Seattle WA 98195-5061.

Model development and archiving support at provided by the following grants: NIH U01HL122199 Analyzing the Cardiac Power Grid, 09/15/2015 - 05/31/2020, NIH/NIBIB BE08407 Software Integration, JSim and SBW 6/1/09-5/31/13; NIH/NHLBI T15 HL88516-01 Modeling for Heart, Lung and Blood: From Cell to Organ, 4/1/07-3/31/11; NSF BES-0506477 Adaptive Multi-Scale Model Simulation, 8/15/05-7/31/08; NIH/NHLBI R01 HL073598 Core 3: 3D Imaging and Computer Modeling of the Respiratory Tract, 9/1/04-8/31/09; as well as prior support from NIH/NCRR P41 RR01243 Simulation Resource in Circulatory Mass Transport and Exchange, 12/1/1980-11/30/01 and NIH/NIBIB R01 EB001973 JSim: A Simulation Analysis Platform, 3/1/02-2/28/07.