Review - (2023) Volume 12, Issue 2

Equation of State for O2-Binding by Hemoglobin in Human Red Blood Cells
Francis Knowles* and Douglas Magde
 
Department of Chemistry and Biochemistry, University of California, San Diego, USA
 
*Correspondence: Francis Knowles, Department of Chemistry and Biochemistry, University of California, San Diego, USA, Email:

Received: 03-Apr-2023, Manuscript No. BABCR-22-19473; Editor assigned: 06-Apr-2023, Pre QC No. BABCR-22-19473 (PQ); Reviewed: 24-Apr-2023, QC No. BABCR-22-19473; Revised: 01-May-2023, Manuscript No. BABCR-22-19473 (R); Published: 09-May-2023, DOI: 10.35248/2161-1009.23.12.481

Abstract

O2-Equilibrium binding data of hemoglobin in whole blood under standard conditions was fitted to an equation of state comprised of three unknown quantities: Kα, the equilibrium constant for binding O2 by equivalent low affinity α-chains; KΔ, a dimensionless equilibrium constant describing the change between low-and high-affinity structures of hemoglobin, Tstate and Rstate; Kβ, the equilibrium constant for binding O2 by equivalent high affinity β-chains. Values of the unknown quantities at pH 7.4 and 37°C are: Kα=15,090 L/mol; KΔ=0.0260; Kβ=393,900 L/mol. A graph of predicted versus observed values of fractional saturation, F, is linear: FPRE=0.9998 FOBS=0.0005, R2=0.9997. The Perutz/Adair equation of state is defined as such insofar as all aspects of the stereo chemical are imposed on the earlier sequential binding model of Adair. The Perutz/Adair equation of state is general, describing: (i) the CO equilibrium binding curve of whole blood under standard conditions, Kα=4.27 × 106 L/mol, KΔ=0.05741, and Kβ=99.1 × 106 L/mol; (ii) the O2-equilibrium binding curve of purified hemoglobin in 0.100 M NaCl, 0.050 M BisTris, pH 7, 20°C, Kα=5.34 × 104 L/mol, KΔ=0.03252, and Kβ=1.81 × 106 L/mol.

Keywords

Hemoglobin, O2-Binding, Red blood cells, Stereo chemical model

Key Points

• In the presence of stoichiometric concentrations of 2,3-Bisphosphoglycerate (BPG), O2-free human hemoglobin experiences enhanced proximal strain in α-subunits and steric hindrance at the surface of the β-chains

• In the presence of stoichiometric concentrations of BPG, O2-free human hemoglobin exhibits chain heterogeneity: α-chains reacting with O2 before β-chains react with O2

• α-Chains exhibit low affinity equivalent O2-binding. β-Chains exhibit high-affinity equivalent O2-binding.

• A structural change, Tstate to Rstate, occurs after equivalent O2-binding by α-chains and before equivalent binding by β-chains

• An equation of state for the O2-equiibrium binding curve is comprised of only three unknown quantities: an equilibrium constant for O2-binding by equivalent α-chains, Kα; an equilibrium constant for O2-binding by equivalent β-chains, Kβ; an equilibrium constant for the Tstate to Rstate change, KΔ.

Introduction

Values for Kα, KΔ and Kβ address only one of three separate aspects of allosteric architecture in human hemoglobin in vivo, the other two being: (i) exothermic binding of Hb4 with BPG, converting Rstate to Tstate; and (ii) changes in pH of the interior of red blood cells (manifestation of the Bohr Effect). The Perutz/Adair equation defines the value of the equilibrium constant for the Tstate → Rstate structure change. The value of KΔ accounts for the ability of red blood cells to release O2 efficiently and provides a basis for high rates of resting metabolism accounting for warm blood. Upon release of O2 from β-chains the Rstate structure collapses to the Tstate structure. The Perutz/Adair equation does not describe O2-equilibrium binding data under all conditions. O2-binding data obtained with purified human hemoglobin in 0.050 M potassium phosphate over a range of temperatures from 4°C to 30°C: does not, for example, demonstrate equivalent binding by β-chains (Knowles, unpublished results).

The stereo chemical model describing steps involved in conversion of Hb4/BPG to (HbO2)4/BPG occurs in three stages [1,2]. In the first stage of the overall sequence of five discrete steps, both Tstate α-chains undergo identical and equivalent O2-binding reactions, remaining in the unchanged Tstate. In the first stage, β-chains are precluded from reaction with O2 by steric hindrance at the distal surface of equivalent β-chain heme moieties. The second stage in the sequence of five distinct reactions is a change in the structure of ((αO2)22/BPG, from Tstate to Rstate. The structure change (i) relaxes steric hindrance at the surface of the β-chains, opening a pathway for the third stage, and (ii) relaxes proximal strain in Tstate αO2-chains, the conversion to Rstate αO2-chains being accompanied by a negative free energy change. BPG binds to the Rstate structure. The third stage, consisting of a sequence of equivalent O2-binding reactions by Rstate β-chains, is free of changes in structure at constant pH. Equation describes the three stages of reactions leading from Hb4/BPG to HbO2)4/BPG under standard conditions of temperature, pH, and the presence of CO2.

Sequence of reactions for conversion of Hb4/BPG, species I, to (Hb4O2)4/BPG, species VI, in three stages. K1, K2, K4, and K5 are equilibrium constants for O2-binding reactions. KΔ is the equilibrium constant for the Tstate to Rstate change. The change in structure occurs between species III and IV. Conformation state for each species is indicated by a left superscript T or R. DPG is represented as in T-state conformations and o in R-state conformations.

Equation

F is fractional saturation of the binding sites of human hemoglobin with the sixth axial ligand: O2 or CO. Analytical expressions relating the concentration of each of the species, II thru VI, to the equilibrium constants, K1, K2, KΔ, K4, and K5, and concentration of O2 are as follows:

Equation

O2-binding sites of α-chains are assigned the property of being identical and equivalent. O2-binding sites of β -chains are also assigned the property of being identical and equivalent. This permits redefinition of the value of equilibrium constants for O2- binding sites: K1=2Kα; K2=Kα/2; K4=2Kβ; K5=Kβ/2. Four unknown quantities K1, K2, K4, and K5, are replaced by four expressions containing only two unknown quantities: Kα and Kβ. Substitution of these statistical equivalents into the analytical expressions for species II thru VI, leads to the Perutz/Adair equation, Equation 2.

Equation

Equation 2 expands the Adair equation to accommodate structural constraints inherent in the Perutz stereo chemical model. These steps serve to reduce the five unknown quantities of Equation 1, to three unknown quantities. These assumptions require that Equation (2) with three unknown quantities accurately predicts equilibrium binding curves for O2. Results presented below show that Equation (2) describes equilibrium binding curves for (i) both O2 and CO, in whole blood under standard conditions and for (ii) purified human hemoglobin in 0.050 M BisTris, 100 M NaCl, pH 7 with HCl, 20°C. Initial tests of Equation (2) were carried out with purified human hemoglobin in the presence in 0.050 M KPi. β-chains do not exhibit equivalent O2-binding in the presence of 0.050 M KPi, pH 7.0, 20.0°C. Four unknown quantities are required to obtain fits of O2-binding data with correlation coefficients of 0.999. Non- equivalent binding to β-chains in 0.050 M KPi, together with marked differences in the value of KΔ, provides insight into the role of E-molecules in the manifestation of allosteric properties of enzymes demonstrating cooperative properties, in general [3-6].

Literature Review

The O2-Equilibrium Binding Curve of Purified Human Hemoglobin in 0.050 M BisTris, 0.100 M NaCl, pH 7.00, 20.0°C. O2-Equilibrium binding data obtained in 0.100 M NaCl, 0.050 M BisTris, pH 7 with HCl, 20.0°C was fitted to Equation 2. The best- fitting values of Kα, Kβ, and KΔ, together with statistical information on the curve-fitting procedure, are summarized in Table 1. Binding data for O2 (Knowles, Architecture of Allosteric Structure. Equation 1 is well fitted by Equation 2, the value of R2 being 0.9997.

Kα 5.344 × 104 L/mol
KΔ 0.03252
Kβ 1.809 × 106 L/mol
R2 0.9997
RMSE 0.005455

Table 1: Curve fitting of the O2 equilibrium binding curve of purified human Hb4 in 0.050 M BisTris, 0.100 M NaCl, pH 7.00 with HCl, 20.0°C to the Perutz-Adair equation, Equation 2.

The first stage is characterized by a low value for Kα , 5.344 104 L/ mol, relative to the value of Kβ, 1.809 106 L/mol. The O2-affinity expressed by the β-chains is greater than that of the α-chains by a factor of 33.85. α-Chains, nevertheless, are titrated with CO before β-chains (Knowles). These results support the assignment of distal side steric hindrance of β-chain heme moieties of Hb4/ BPG observed by Perutz. Purified human Hb4/(Cl-1)n in 0.050 M BisTris, 0.100 M NaCl appears to behave as an ideal cooperative dimer, in which the O2-affinity of the set of equivalent β -chains is elevated from a value significantly less than that of the first set of equivalent α-chains to a value almost 34-fold higher than that expressed by the first set of equivalent α-chains, upon binding of O2. In the case of human Hb4/BPG the cooperative mechanism: (i) relaxes proximal strain in Tstate αO2-chains and (ii) provides unhindered access to the distal surface of the β-chain heme moiety to O2. O2-binding by the pair of equivalent Tstate α-chains both creates and enables the switching mechanism controlling access of pre-existing sterically blocked Rstate β-chains to O2 molecules.

The ability of Equation 2 to predict observed data was tested by comparing predicted and observed values of F of hemoglobin with O2. Predicted values of F, FPRE, are plotted against the observed values of F, FOBS, in Figure 1. This plot is expected to be linear, represented by the equation y=x, if the equation of state accurately describes the values of FOBS. The best fitting straight line has the equation FPRE=0.9998 FOBS=-0.0005, quite close to y=x. Deviation of the best fitting straight line from y=x is not visible to the unaided eye.

Analytical-Biochemistry-predicted

Figure 1: Observed and predicted values of F of the O2-Equilibrium binding curve of purified Hb4 in 0.100 M NaCl , 0.05 M BisTris, pH 7.00 with HCl, 20.0oC. Predicted values of F, FOBS, were calculated with Equation (2), using the best fitting values of Kα, Kβ, and KΔ, given in Table 1, and the corresponding O2-values. y=0.9998x-0.00065, R2=0.9997

The Perutz/Adair Equation 2 is found to accurately describe observed values of FOBS with three unknown quantities. The general form of the Adair equation should be replaced by Equation 2. Gilbert Smithson Adair was correct in recognizing the conditions justifying the use of equivalent binding constants. The results obtained by the Perutz/Adair equation, constrained to equivalent binding by identical chains are subject to a level of interpretation not possible with the results returned by curve fitting procedures using the general form of the Adair equation.

The CO Equilibrium Binding Curve of Whole Blood in Isotonic Saline, pH 7.4, PCO2=40 torr, 37°C (Data of Roughton 1970). The CO equilibrium binding curve of whole human blood was fitted to the Perutz-Adair equation. Roughton’s data (1970) presented CO concentrations as torr, in the gas phase in equilibrium with the solution containing red blood cells. Units in torr were converted to the concentration of CO in the liquid phase in units of μmol/L at 37 degrees C: 1.000 torr=1.0356 μmol CO/L. The best-fitting values of Kα, Kβ, and KΔ, together with statistical information on the curve fitting procedure, are summarized in Table 2. Roughton’s CO-equilibrium binding data is well fitted by Equation 2.

Kα 4.274 × 106 L/mol
KΔ 0.05741
Kβ 99.06 × 106 L/mol
R2 0.9997
RMSE 0.00654

Table 2: Curve fitting of the CO-equilibrium binding curve of whole blood, pH 7.4, 37°C, under standard conditions, to the Perutz-Adair equation, Equation 2.

The first stage is, again, characterized by a low value for Kα, 4.274 × 106 L/mol, relative to the value of Kβ, 99.06 × 106 L/mol. The O2- affinity expressed by the β-chains is greater than that of the α-chains by a factor of 23.17. The α-chains, nevertheless, are titrated with CO before the β-chains. These results support the assignment of distal side steric hindrance of β-chain heme moieties in Tstate structures observed by Perutz. Whole blood, under standard conditions, behaves as an ideal cooperative dimer, in which the CO-affinity of equivalent β-chains is elevated from a value significantly less than that of the first set of equivalent α-chains to a value 23-fold higher than that expressed by equivalent Tstate α-chains. In the case of hemoglobin, within RBC’s, the cooperative mechanism relies on control of steric hindrance at the distal surface of the β-chain heme moiety and this switch responds to CO-binding by the pair of equivalent T-state αO2-chains. The results obtained with whole blood, under standard conditions, are similar to those obtained with purified hemoglobin in 0.050 M BisTris, 0.100 M NaCl, pH 7.00, 20.0°C, insofar as the pattern of returned values are concerned. The actual values are, of course, quite different, insofar as CO is more tightly bound than O2. The value of KΔ, however, is quite independent of the structure of the sixth axial ligand, being 0.05741 for binding of CO by whole blood and 0.03252 for binding of O2 by purified hemoglobin in the presence of 0.100 M NaCl.

The ability of Equation 2 to describe observed data was tested by comparing predicted and observed values of fractional saturation of whole blood with CO. Predicted values of F, FPRE, are plotted against the observed values of F, FOBS, in Figure 2. Again, this plot is expected to be linear, (y=x), if the equation of state accurately predicts the values of FOBS. The best fitting straight line has the equation.

Analytical-Biochemistry-fractional

Figure 2: Observed and predicted values of fractional saturation of the CO- equilibrium binding curve of whole blood under standard conditions. Predicted values of F, FPRE, were calculated using Equation (2) and the best fitting values of Kα, Kβ, and KΔ, given in Table 2, and the corresponding CO-values. y=0.9998x-0.00065, R2=0.9997

FPRE=0.9915 FOBS+0.0052, quite close to y=x. The two lowest points deviate from the best fitting line.

The Perutz/Adair equation, Equation 2, with three unknown quantities, accurately and precisely describes Roughton’s CO-equilibrium binding curve for whole blood, under standard conditions.

The O2-Equilibrium Binding Curve of Whole Blood in Isotonic Saline, pH 7.4, p(CO2)=40 torr, 37δC [7].

The O2-equilibrium binding curve of whole blood, obtained under standard conditions, combining data collected by (1972) is presented in Figure 3. O2 concentrations expressed in units of torr in the gas phase were converted to units of concentration in the liquid phase in equilibrium with the gas phase. The concentration of O2 in the liquid phase in units of μmol O2/L is equal to p(O2) in units of torr multiplied by 1.00023. The binding curve was fitted to the Perutz-Adair equation. Best-fitting fitting values of Kα, Kβ, and KΔ, together with statistical information on the curve fitting procedure are summarized in Table 3.

Kα 15.09 × 103 L/mol
KΔ 0.02602
Kβ 393.9 × 103 L/mol
R2 0.9998
RMSE 0.003609

Table 3: Curve fitting of the O2-equilibrium binding curve of whole blood, under standard condition, to Equation 2.

Analytical-Biochemistry-Equilibrium

Figure 3: O2-Equilibrium binding curve of whole human blood under standard conditions.

Again, the first stage is characterized by the low value for Kα of 15.09 × 103 L/mol, relative to the value of Kβ, 393.9 × 103 L/mol. The O2-affinity expressed by the β-chains is greater than that of the α-chains by a factor of 26.10. Whole blood, under standard conditions, appears to behave as an ideal cooperative dimer, in which the O2-affinity of the second set of equivalent β-chains is elevated from a value significantly less than that of the first set of equivalent α-chains to a value 26-fold higher than that expressed by the first set of equivalent α-chains, upon binding of O2 [8-16]. Results obtained with whole blood for O2, under standard conditions, are similar to those obtained with (i) purified hemoglobin in 0.050 M BisTris, 0.100 M NaCl, pH 7.00, 20.0°C and (ii) whole blood for CO, under standard conditions, insofar as an identical pattern of returned values are concerned. The actual values are, of course, quite different, insofar as O2 is bound much less tightly in whole blood, under standard conditions, than by either (i) purified hemoglobin in the presence of 0.100 M NaCl or (ii) CO in whole blood. The value of KΔ is independent of the structure of the: (i) sixth axial ligand or (ii) composition of the supporting electrolyte. Direct comparison of the value of KΔ with CO and O2 in whole blood, shows only a 2.17-fold difference.

Values of FPRE are plotted against FOBS, in Figure 4. The best fitting straight line has the equation FPRE=0.9955 FOBS=0.004, quite close to y=x.

Analytical-Biochemistry-Predicted

Figure 4: Observed and predicted values of F of the O2-equilibrium binding curve of whole blood under standard conditions. Predicted values of F, FPRE, were calculated with Equation (2), using the best fitting values of Kα, Kβ, and KΔ, given in Table 3, and the corresponding O2- values. y=0.9998x-0.00065, R2=0.9997

The Perutz-Adair equation, Equation 2, with three unknown quantities, accurately and precisely describes the Roughton/ Severinghaus O2-equilibrium binding curve for whole blood, under standard conditions.

Discussion and Conclusion

The Perutz/Adair equation fits equilibrium binding data of Hb4/E for O2 for a wide variety, but not all, E-molecules. The most important E-molecule for human hemoglobin is, of course, naturally occurring BPG. The Perutz/Adair equation is the simplest possible elaboration of the Adair equation, taking into consideration significant discoveries following Adair’s publications. (i) Purple hexagonal crystals of horse Hb4/BPG will shatter and dissolve if O2 is allowed to diffuse into the crystals followed by formation of scarlet monoclinic needle crystals of (HbO2)/BPG. This observation by Felix Haurowitz is a silent argument that the transition from Tstate to Rstate is structural rather than an insignificant conformational change. (ii) The subunits of which a hemoglobin molecule is comprised are not all identical. (iii) Recognition of the role of BPG and other E-molecules on the properties of human hemoglobin was unappreciated. (iv) Formulation of the stereochemical model introduced observations leading to formulation of this manuscript. (v) Demonstration of chain heterogeneity in O2-equilibrium binding curves led to formulation. (vi) Imposition of equivalent binding on both the α-and β-chains reduces the number of unknown quantities required to describe O2-equilibrium binding reactions from 4, in the original form of the Adair equation, to 3 in the Perutz-Adair equation, Equation 2. In this instance a mathematical advance demanded the acceptance of the chemical concept that a macromolecular structure with equivalent binding sites was rigid.

Introduction of the equilibrium constant for the inter-conversion of Tstate and Rstate structures introduced a third unknown quantity, a dimensionless equilibrium constant for the Tstate to Rstate change, KΔ. Equation 2 permits, for the first time, meaningful comparisons of the values of Kα and Kβ obtained from different species, as well as a meaningful interpretation of the value of KΔ.

In each of the equilibrium binding experiments described above, the binary complex of human Hb4 and an E-molecule, behaves in a fashion similar to that of a cooperative dimer comprised of dissimilar subunits. The first subunit, or set of equivalent subunits, regulates the reactivity of the second subunit, or set of equivalent subunits. The Rstate to Tstate change induced by formation of the binary complex of human Hb4 and BPG regulates O2-affinity of α-chains by enhancing pre-existing proximal strain in the bond between an imidazole sidechain and the iron atom of the heme moiety, thereby reducing O2-affinity of α-chains. The Tstate also regulates O2-affinity of β-chains by blocking access of O2 to the distal surface of the β-chain heme moiety. Binding of O2 to β-chains requires (i) saturation of low affinity Tstate α-chains and (ii) over riding an endothermic Tstate to Rstate structure change. Reversal of the exothermic conformation change induced by the E-molecule, ΔG°=-33.7 kJ/mol, requires an energy source. Energy sources are provided by conversion of Tstate αO2-chains to Rstate αO2-chains and exothermic binding of BPG to Rstate structures. The change in energy from the change in state for the pair of Tstate αO2-chains can be estimated if it is assumed that Rstate O2-affinity is similar to that expressed by β-chains, Kα increasing from 15.09 L/mmol to 393.9 L/mmol, for an increase in ΔG° of -16.8 kJ for the pair of α-subunits. This energy would be available to uncouple BPG from Tstate conformations sufficiently for the endothermic Tstate to Rstate change to take place. These numbers alone, closely account for the values of KΔ.

The data for whole blood with O2 and CO can be compared and used to confirm the validity of Haldane’s Second Law. Haldane’s Second Law requires that the equation of state for F of hemoglobin with either O2 or CO be identical in form, the only difference being that the values of the equilibrium constants differ. Determined that equilibrium constants for CO exceed those for O2 by a factor (M-factor) of 235.15 results given in Tables 2 and 3 return an M-factor of 283 for α-chains, in the T-state, and 252 for β-chains, in the R-state. Haldane’s Second Law would appear to be substantiated by these results. The values of M are in good agreement but not as significant as is the fact that the Perutz/Adair equation fits both sets of data: CO-binding and O2-binding. The Perutz-Adair equation can be modified and used to generate O2-equilibrium binding curves in the presence of CO (unpublished data, Knowles). The Perutz-Adair equation, then, holds promise of further clarification of the functional properties of hemoglobin molecules. The Perutz/Adair equation of state, Equation 2, is incomplete. Standard conditions do not, in fact, replicate conditions experienced by red blood cells as they pass through the pulmonary circulation (increasing-pH values) and systemic circulation (decreasing pH-values). A full description of the equation of state must take into account the pH-dependence of F in red blood cells.

The 2-state model does not describe the equilibrium binding curves obtained with either O2 or CO. The 2-state model does not resemble the model presented. This work cannot make a judgment about the ability of 2-state model that may well describe the properties of other enzymes that cannot, in turn, be accommodated by the Perutz/Adair equation. Other methods of investigation may provide additional insights. Analysis of progress curves for the reaction of T((Hb4)/BPG) with O2 or reaction of R((HbO2)4/BPG) with dithionite may either support or refute the model which underlies the Perutz/Adair equation.

The following elements of allosteric structure are implicit in the Perutz/Adair equation. Two cooperative dimeric subunits, (T1α, R2β) and (T2α, R1β), exist in Hb4, in which β-chains regulate α-chains. The order of binding of O2 to a cooperative dimer is first β, then α. An Rstate to Tstate structure change arises from binding of BPG to Hb4, enhancing pre-existing proximal strain in α-chains, thereby reducing the affinity of α-chains for O2 and introducing steric hindrance at the distal surface of β-chain heme moieties. BPG acts as a non-competitive inhibitor. The order in which O2 binds to the subunits of Hb4 is inverted by formation of the Hb4/ BPG binary complex. Tstate α-chains are identical and equivalent in the presence of BPG. Binding of O2 to β-chains requires a structure change from Tstate to Rstate. The Rstate structure formed in response to oxygenation reactions of α-chains binds BPG. The significant free energy change required for the formation of an oxygenated Rstate from the oxygenated Tstate is provided by: (i) the very low value of KΔ; (ii) conversion of Tstate αO2-chains to Rstate αO2-chains; (iii) binding of BPG by the Rstate structure. No structural changes occur upon binding of O2 to Rstate structures. These elements of allosteric structure comprise the architecture of allosteric structure at constant pH.

References

Citation: Knowles F, Magde D (2023) Equation of State for O2-Binding by Hemoglobin in Human Red Blood Cells. Biochem Anal Biochem. 12:480

Copyright: © 2023 Knowles F, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.