TY - JOUR AU - Yamamoto, Yasuhiko AB - Abstract Human adult haemoglobin (Hb A), a tetrameric oxygen transfer haemoprotein, has been recognized as an excellent model for investigating the structure–function relationships in allosteric proteins, and has been characterized exhaustively from both experimental and theoretical aspects. Despite the detailed structural and spectroscopic information available for the protein, functional properties have not been as fully elucidated as expected, and hence have remained unexplored. A major drawback for the functional characterization of Hb A is the lack of experimental techniques which enable quantitative characterization of functional properties of the individual subunits of the intact protein. In this study, we have developed techniques for determining the equilibrium constant of the acid–alkaline transition, usually represented as the ‘pKa’ value, in the individual subunits of the met-forms of Hb A (metHb A) and human foetal haemoglobin (metHb F). The pKa values of the individual subunits of metHb A and metHb F have been shown to constitute novel and highly sensitive probes for characterizing the effects of structural changes of not only the interfaces between the subunits within the protein, but also the contact between haem and the protein in the haem pocket. In addition, haem replacement studies of the proteins revealed that the contact between the haem peripheral vinyl side chain and the protein in the haem pocket is important for maintaining the non-equivalence in the haem environment between the subunits of Hb A and Hb F, which could be relevant to the cooperative ligand binding of the proteins. acid–alkaline transition, 19F NMR, haemoglobin, myoglobin, thermodynamics Human adult haemoglobin (Hb A), a multimeric oxygen transport haemoprotein, is the best, albeit incompletely, understood allosteric protein in terms of its structure and function (1). Hb A has served as an excellent model for investigating the structure–function relationships in allosteric proteins. Hb A is a tetrameric protein consisting of two α subunits of 141 amino acid residues and two β subunits of 146 residues, respectively (2). The α and β subunits each contain one haem [Fe-protoporphyrin IX complex (protohaem) (Fig. 1A)] as a prosthetic group, to which oxygen molecule binds. The structures of the α and β subunits are quite similar to each other and also to that of myoglobin (Mb), a monomeric oxygen storage haemoprotein, even though the sequence homology among them is only ∼18% (3). Fig. 1 View largeDownload slide Molecular structure of haem and orientation of haem relative to His F8. (A) The structures and numbering system for protohaem (R2 = CH3, R3 = R8 = CH=CH2, R7 = CH3), mesohaem (R2 = R7 = CH3, R3 = R8 = C2H5) and 7-PF (R2 = CH3, R3 = R8 = C2H5, R7 = CF3). (B) Two possible orientations of the haem relative to His F8, N and R forms. The Φ value represents the angle between the projection of the axial His imidazole onto the haem plane and the NII–Fe–NIV axis. Fig. 1 View largeDownload slide Molecular structure of haem and orientation of haem relative to His F8. (A) The structures and numbering system for protohaem (R2 = CH3, R3 = R8 = CH=CH2, R7 = CH3), mesohaem (R2 = R7 = CH3, R3 = R8 = C2H5) and 7-PF (R2 = CH3, R3 = R8 = C2H5, R7 = CF3). (B) Two possible orientations of the haem relative to His F8, N and R forms. The Φ value represents the angle between the projection of the axial His imidazole onto the haem plane and the NII–Fe–NIV axis. In general, the functional properties of such proteins are characterized spectrophotometrically (4), since haem exhibits a characteristic absorption spectrum which sharply reflects changes in the haem Fe coordination state (5). Despite the high sensitivity of absorption spectra to the local haem environment, however, the line widths of the absorption spectra are too large to differentiate the haem environments of the individual subunits of Hb A. There have been a number of studies attempting to characterize the functional properties of the individual subunits of Hb A, most of which have utilized artificial proteins, in which one of the subunits is chemically modified, e.g. metal ion substitution (6–16), to ruin its function in order to selectively evaluate the properties of the other intact subunit. Nevertheless questions must remain about the utilization of such chemically-modified proteins to isolate the functional properties of one of the subunits because of the presence of indirect effects of the modification. Hence, the development of techniques which enable quantitative characterization of the functional properties of the individual subunits of intact Hb A has been highly demanded, because they are expected to provide experimental data which directly reflect the subunit interactions responsible for the cooperative oxygen binding of Hb A, and hence important clues for elucidating the structure–function relationship in the protein. In fact, taking advantage of subunit-specific isotope labelling techniques applicable to Hb A (17–20), the preparation of the so-called valency hybrid Hb A possessing ferrohaem in one subunit and ferrihaem in the other, (15, 21) or the high sensitivity of paramagnetic NMR in distinguishing the individual subunits of the protein (22, 23), some attempts have been successful in the characterization of the functional properties of the individual subunits of the intact protein. In this study, we develop techniques which allow characterization of ‘the acid–alkaline transition’ in the subunits of the met form of Hb A (metHb A). The haem active sites in metHb A and the met forms of Mb (metMb) exhibit characteristic pH-dependent structure changes collectively known as the acid–alkaline transition (5). The subunits of metHb A as well as metMb possessing highly conserved distal His E7 [E7 represents an alphanumeric code referring to the positions of amino acid residues in the helices and turns of Hb and Mb, E7 represents the seventh residue in the E helix (see the amino acid sequences of the proteins in Supplementary Material SM1)] have H2O and OH− as coordinated external ligands under low and high pH conditions, respectively (5, 24–31). As illustrated in Scheme 1 (32), the acid–alkaline transition in the proteins comprises three distinct reactions, i.e. interconversion of the coordinated ligand between H2O and OH−, tautomerism of the His E7 imidazole and deprotonation/protonation of His E7 NδH (31). Furthermore, since the transition is associated with the deprotonation/protonation process, its equilibrium constant is usually represented as the ‘pKa’ value in analogy with the ionization equilibrium of an acid. The pKa value has been shown to remarkably vary with the protein, ranging from 7.2 to 10 (5, 24–31), indicating that the acid–alkaline transition in the protein sharply reflects characteristics of the structural features of the haem active site (24–31). In addition, we have recently demonstrated both experimentally and theoretically that the pKa value in Mb is highly sensitive to electronic nature of haem peripheral side chains (33). Furthermore, the sensitivity of the pKa value to the haem environment has also been confirmed by the observation that the value is affected by the well-documented haem orientational disorder (34) resulting from the incorporation of haem into the apoprotein in two orientations differing by 180° rotation about the 5,10-meso axis (Fig. 1B) (31). Scheme 1 View largeDownload slide The acid–alkaline transition in metMb, and the subunits of metHb A and metHb F. Scheme 1 View largeDownload slide The acid–alkaline transition in metMb, and the subunits of metHb A and metHb F. One of our techniques is based on quantitative fitting of the pH-dependent change in the absorption spectra of the protein to the sum of two Henderson–Hasselbalch equations, with the assumption of non-equivalent pKa values for the constituent α and β subunits of metHb A. The other two techniques take advantage of the preparation of the valency hybrid metHb A (35), in which ferrihaem is accommodated in one subunit, i.e. either the α or β subunit, and ferrohaem in the other, and 19F NMR in combination with fluorinated haem. For the former technique, the acid–alkaline transition in only the subunit possessing ferrihaem is reflected in the pH-profile of the absorption spectra of the prepared hybrid protein. For the latter one, as reported previously (36–41), NMR signals arising from fluorine atoms introduced into the haem as peripheral side chain(s) are extremely sensitive to the haem electronic structure. In particular, fluorinated haems possessing CF3 group(s) such as 13,17-bis(2-carboxylatoethyl)-3,8-diethyl-2,12,18-trimethyl-7-trifluoro-methyl-porphyrinatoiron (III) (36) (7-PF: Fig. 1A) exhibit relatively narrow 19F NMR signals, even though the haem is incorporated into an apoprotein with a large molecular weight, mainly because of its relatively small chemical shift anisotropy, which largely contributes to the 19F relaxation (42, 43). In addition, the Curie spin relaxation mechanism effect (44, 45) is also expected to be small for the 19F relaxation of CF3 due to its rapid rotation around the C3 axis. These three techniques have been applied to determine the pKa values of the subunits of not only metHb A, but also the met form of human foetal haemoglobin (46) (metHb F) composed of two α subunits and two γ subunits of 146 amino acid residues, respectively, in order to characterize the haem active sites in the individual subunits of the proteins. Furthermore, the pKa values of the α, β and γ subunits of metHb A and metHb F were compared with those of the corresponding subunits in their isolated states to gain an insight into the effect of the hetero-tetramer assembly on the haem environments of the subunits. Considering that the isolated α and β subunits, which exist as a dimer (47) and a tetramer (47–49), respectively, do not exhibit cooperative ligand binding at all, the non-equivalence in the haem environment between the constituent subunits within a tetrameric Hb, together with the molecular mechanisms which allow transmission of structural changes among the haem active sites of the subunits, is thought to be indispensable for the cooperativity of Hb. The high sensitivity of the pKa value to the nature of the haem active site is expected to permit detailed characterization of the haem environments of the subunits of metHb A and metHb F. Finally, the pKa values of Hb A and Hb F reconstituted with mesohaem (Fig. 1A) and 7-PF were determined and compared with each other and with those of native Hbs in order to determine the effect of the haem-protein interaction on the haem environments of the constituent subunits. We report herein the development of techniques for determination of the pKa values of the subunits of metHb A and metHb F. The subunit interactions in the proteins was clearly manifested in the pKa values of metHb A, metHb F, the isolated subunits and valency hybrid proteins. In addition, since the pKa value has been shown to be highly affected by the electronic nature of the haem peripheral side chains, we have analysed the relationship between the pKa value and the conformation of the haem vinyl side chains through the analysis of well-known Mulliken charge analysis using density functional theory (DFT) calculations. The calculations suggested that the pKa value is largely affected through the electronic interaction between the porphyrin and vinyl π-systems. Finally, comparison of the pKa values between metHb A and the protein reconstituted with mesohaem or 7-PF revealed that the contact between the haem peripheral vinyl side chains and the protein is important for maintaining the non-equivalence in the haem environment between the α and β subunits of the protein. Materials and Methods Protein Samples Hb A was prepared from blood obtained from the Medical Center of the University of Tsukuba using the reported procedure (50). Hb F was isolated and purified in the carbonmonoxy (CO) form from the blood withdrawn from the umbilical cord of a patient, who agreed to the donation, in the Medical Center according to the method previously described (51). α, β and γ subunits were isolated from Hb A and Hb F according to the method previously described (52, 53). The quality of each Hb subunit was confirmed with a Mass Spectrometer, QStar/Pulsar i (Applied Biosystems). Valency hybrid α-metHb A (α-metHb F) possessing ferrihaem in the α subunit and ferrohaem in the β (γ) one was prepared from metHb A (metHb F) using the reported procedure (54). Valency hybrid β-metHb A (γ-metHb F) possessing ferrohaem in the α subunit and ferrihaem in the β (γ) one was prepared through reassembly of the isolated α, β and γ subunits according to the method previously described (55). In the valency hybrid proteins prepared, CO is coordinated as an external ligand to haem Fe of ferrohaem. The quality of each hybrid protein was confirmed by analysis of 1H NMR spectra of their azide adducts (Fig. 2). MetHb A and metHb F were prepared from the CO forms of the proteins under a stream of O2 gas with strong illumination in the presence of a 5-fold molar excess of potassium ferricyanide (Wako Chemical Co.). The protein was separated from the residual chemicals using a Sephadex G-50 (Sigma Chemical Co.) column equilibrated with 50 mM Bis–Tris (Sigma Chemical Co.), pH 6.5. 7-PF was synthesized as previously described (31). Protohaem was purchased from Sigma Chemical Co. Mesohaem (Fig. 1A) was prepared from mesoporphyrin IX dimethyl ester purchased from Aldrich Chemical Co. Sperm whale Mb was purchased as a lyophilized powder from Biozyme and used without further purification. The apoprotein of Hb was prepared at 4°C according to the procedure of Teale (56), and reconstitution of the apoprotein with haem was carried out by the standard procedure (50). The reconstituted Hb A and Hb F were concentrated to ∼1 mM in an ultra-filtration cell (Amicon). The pH of each sample was measured with a Horiba F-22 pH meter equipped with a Horiba type 6069-10c electrode. The pH of the sample was adjusted using 0.2 M NaOH or HCl. Fig. 2 View largeDownload slide 400 MHz 1H NMR spectra of met-azido Hbs. 1H NMR spectra (400 MHz) of the met-azido forms of valency hybrid α-metHb A (A), metHb A (B), valency hybrid β-metHb A (C), valency hybrid α-metHb F (D), metHb F (E) and valency hybrid γ-metHb F (F) at pH 7.0. (A–C) and (D–F) were recorded at 25 and 45°C, respectively. The corresponding haem methyl proton signals are connected by dotted lines. In the spectra of valency hybrid Hbs, signals due to the subunit possessing a ferrohaem with a coordinated CO are not observed in the indicated chemical shift region, because of the absence of unpaired electrons. Fig. 2 View largeDownload slide 400 MHz 1H NMR spectra of met-azido Hbs. 1H NMR spectra (400 MHz) of the met-azido forms of valency hybrid α-metHb A (A), metHb A (B), valency hybrid β-metHb A (C), valency hybrid α-metHb F (D), metHb F (E) and valency hybrid γ-metHb F (F) at pH 7.0. (A–C) and (D–F) were recorded at 25 and 45°C, respectively. The corresponding haem methyl proton signals are connected by dotted lines. In the spectra of valency hybrid Hbs, signals due to the subunit possessing a ferrohaem with a coordinated CO are not observed in the indicated chemical shift region, because of the absence of unpaired electrons. UV–vis absorption and NMR spectroscopies UV–vis absorption spectra were recorded at 25°C with a Beckman DU 640 spectrophotometer and a protein concentration of 10 µM in 20 mM phosphate buffer. 1H NMR and 19F NMR spectra were recorded on a Bruker AVANCE-400 spectrometer operating at the 1H and 19F frequencies of 400 and 376 MHz, respectively. Typical 1H and 19F NMR spectra consisted of ∼20 k transients with a 100 kHz spectral width and 16 k data points. The signal-to-noise ratio of each spectrum was improved by apodization, which introduced 20–100 Hz line broadening. The chemical shifts of 1H and 19F NMR spectra are given in ppm downfield from the residual 1H2HO, as an internal reference, and from trifluoroacetic acid, as an external reference, respectively. DFT calculation Evaluation of the electron density of the haem Fe atom (ρFe) value was carried out by Mulliken charge analysis through the DFT calculations carried out using the Gaussian 03 programme package (57). Restricted spin orbital approach involving the B3LYP method, together with electron basis sets of Pople’s 6-31G(d), was employed. For simplification, the calculations were performed for the 3,8-divinylporphinatoiron(II) complex as a model for the protohaem (Supplementary Material SM2B). Geometry optimization of the model compounds was carried out in the gas phase, with a fixed vinyl group conformation possessing the dihedral angles C2–C3–Cα–Cβ (Ψ3-vinyl) and C7–C8–Cα–Cβ (Ψ8-vinyl) of 0° (Fig. 3), i.e. all vinyl group atoms were in the porphyrin plane (Supplementary Material SM2B), and then single point calculations of the model compound were carried out with changing the Ψ3-vinyl and Ψ8-vinyl values by 5° steps between 0° and 90° in order to analyse the relationship between the ρFe value and dihedral angle (Fig. 3). The oxidation, spin and coordination states of the haem Fe atoms of the model compound were assumed to be Fe(II), S = 0, and a six-coordinated structure with CO ligands, respectively (Supplementary Material SM2B). Fig. 3 View largeDownload slide Plots of calculated Mulliken charges of haem Fe atom against the Ψ8-vinyl value, with a fixed 3-vinyl side chain conformation at Ψ3-vinyl = 0°, (open circle) and the Ψ3-vinyl and Ψ8-vinyl values, which changed simultaneously by 5° steps (filled circle). The definition of the dihedral angles C2–C3–Cα–Cβ (Ψ3-vinyl) and C7–C8–Cα–Cβ (Ψ8-vinyl) is illustrated in the inset; s-cis vinyl side chain conformation with all vinyl group atoms being in the porphyrin plane is defined to be 0°. Fig. 3 View largeDownload slide Plots of calculated Mulliken charges of haem Fe atom against the Ψ8-vinyl value, with a fixed 3-vinyl side chain conformation at Ψ3-vinyl = 0°, (open circle) and the Ψ3-vinyl and Ψ8-vinyl values, which changed simultaneously by 5° steps (filled circle). The definition of the dihedral angles C2–C3–Cα–Cβ (Ψ3-vinyl) and C7–C8–Cα–Cβ (Ψ8-vinyl) is illustrated in the inset; s-cis vinyl side chain conformation with all vinyl group atoms being in the porphyrin plane is defined to be 0°. Results Determination of the pKa values of the constituent subunits of metHb A and metHb F The pH dependence of the absorption spectra of metHb A at 25°C and plots of 575-nm absorption against pH are presented in Fig. 4A and B, respectively. Quantitative fitting of the plots to the Henderson–Hasselbalch equation yielded an ‘apparent pKa value’ of 7.90 ± 0.05 (Supplementary Material SM3). Careful scrutiny revealed that the fitting with a single Henderson–Hasselbalch equation resulted in systematic errors as to the observed and theoretical values, i.e. the observed values were larger and smaller than the theoretical ones in lower and higher pH regions relative to the apparent pKa value, respectively (Supplementary Material SM3). With the assumption of non-equivalence in the pKa value between the α and β subunits of metHb A, the plots were significantly better fitted with the sum of two Henderson–Hasselbalch equations, yielding pKa values of 7.48 ± 0.05 and 8.32 ± 0.05 (Supplementary Material SM2). Fig. 4 View largeDownload slide Analysis of the acid-alkaline transition of metHb A. pH-dependent optical spectra, 450–700 nm, of metHb A (A), and pKa values of 7.48 ± 0.05 and 8.32 ± 0.05 obtained from plots of absorbance at 575 nm against pH (B). Fig. 4 View largeDownload slide Analysis of the acid-alkaline transition of metHb A. pH-dependent optical spectra, 450–700 nm, of metHb A (A), and pKa values of 7.48 ± 0.05 and 8.32 ± 0.05 obtained from plots of absorbance at 575 nm against pH (B). In order to assign the obtained pKa values to the subunits of metHb A, the pKa value of valency hybrid α-metHb A (β-metHb A) possessing ferrihaem in the α (β) subunit and ferrohaem in the β (α) one was similarly determined. In these hybrid proteins, the acid–alkaline transition in only the subunit possessing ferrihaem is reflected in the pH-profile of the absorption spectra. The pKa values of 8.27 ± 0.10 and 7.61 ± 0.10 were obtained for α-metHb A and β-metHb A, respectively (Supplementary Material SM4). These results unequivocally indicated that the pKa values of 8.32 ± 0.05 and 7.48 ± 0.05, obtained on fitting of the metHb A plots using the sum of two Henderson–Hasselbalch equations, are due to the α and β subunits of metHb A, respectively. Similarly, on fitting of the pH-dependent 575-nm absorption of metHb F with the sum of two Henderson–Hasselbalch equations, the pKa values of 7.76 ± 0.05 and 8.48 ± 0.05 were obtained (Supplementary Material SM5). Considering the pKa values of 8.27 ± 0.1 and 7.89 ± 0.10 obtained for valency hybrid α-metHb F and γ-metHb F, respectively (Supplementary Material SM6), the 8.48 ± 0.05 and 7.76 ± 0.05 yielded for metHb F could be due to the pKa values of the α and γ subunits of metHb F, respectively. The appreciable differences in the pKa values between the corresponding subunits of metHb A (or metHb F) and valency hybrid metHb A (or metHb F) indicated that the structural properties of one subunit in Hb A (or Hb F) are affected by the haem Fe oxidation and coordination states of the other. Determination of the pKa values of the isolated α, β and γ subunits We next similarly determined the pKa values of the isolated α, β and γ subunits through analysis of the pH dependence of the optical spectra. pKa values of 8.17 ± 0.05, 7.75 ± 0.05 and 8.26 ± 0.05 were obtained for the isolated α, β and γ subunits, respectively (Fig. 5 and Table I), which were different from those of the corresponding subunits of metHb A and metHb F in Table I, indicating that the pKa values of the subunits were significantly affected by the hetero-tetramer assembly. Fig. 5 View largeDownload slide Analysis of the acid-alkaline transition of met-forms of isolated subunits. pH-dependent optical spectra, 450–750 nm, of the isolated α subunit (A), β subunit (B) and γ subunit (C), and pKa values of 8.17 ± 0.05, 7.75 ± 0.05 and 8.26 ± 0.05 obtained from plots of 575 nm-absorbance for the isolated α, β and γ subunits, respectively (D). Fig. 5 View largeDownload slide Analysis of the acid-alkaline transition of met-forms of isolated subunits. pH-dependent optical spectra, 450–750 nm, of the isolated α subunit (A), β subunit (B) and γ subunit (C), and pKa values of 8.17 ± 0.05, 7.75 ± 0.05 and 8.26 ± 0.05 obtained from plots of 575 nm-absorbance for the isolated α, β and γ subunits, respectively (D). Table I. The pKa values of metHb A, metHb F, Hybrid Hbs and isolated subunit at 25°C.   metHb A  metHb F  Hybrid Hb  Isolated subunit  α  8.32 ± 0.05  8.48 ± 0.05  8.27 ± 0.10  8.17 ± 0.05  β  7.48 ± 0.05  –  7.61 ± 0.10  7.75 ± 0.05  γ  –  7.76 ± 0.05  7.89 ± 0.10  8.26 ± 0.05    metHb A  metHb F  Hybrid Hb  Isolated subunit  α  8.32 ± 0.05  8.48 ± 0.05  8.27 ± 0.10  8.17 ± 0.05  β  7.48 ± 0.05  –  7.61 ± 0.10  7.75 ± 0.05  γ  –  7.76 ± 0.05  7.89 ± 0.10  8.26 ± 0.05  View Large Absorption spectra of Hb A, Hb F and their constituent subunits We measured the absorption spectra of the CO forms of Hb A, Hb F and their constituent subunits, and the absorption maxima (λmax) are summarized in Table II. The λmax values observed for Hb A and isolated α and β subunits were similar to those previously reported (58). Assuming that the λmax values of Hb A and Hb F represent those of the constituent subunits, comparison of the λmax values among the Hbs and isolated subunits indicated that the hetero-tetramer assembly results in bathochrimic and hypochromic shifts of the absorption of α subunit, and β and γ subunits, respectively. Since the Soret absorption has been shown to exhibit bathochromic shift with increasing the polarity of solvent (59), the present results supported that the haem environments of the α subunit, and the β and γ subunits become more polar and less polar upon the hetero-tetramer assembly. Table II. The absorption maximum wavelength (λmax) of CO forms of Hb A, Hb F and their constituent subtunit at pH 7.00 and 25°C.   Soret (nm)   Visible (nm)     Hb A or Hb F  Isolated subunits  Hb A or Hb F  Isolated subunits  α  419.5  418.5  537.0  538.5a  568.0  568.5  β  419.5  420.0  540.0  538.5  569.5  568.5  γ  419.5  420.8  539.2  538.0  569.4  568.5    Soret (nm)   Visible (nm)     Hb A or Hb F  Isolated subunits  Hb A or Hb F  Isolated subunits  α  419.5  418.5  537.0  538.5a  568.0  568.5  β  419.5  420.0  540.0  538.5  569.5  568.5  γ  419.5  420.8  539.2  538.0  569.4  568.5  aTaken from the value of Hb A. The value of Hb F is 538.0 nm. View Large Determination of the pKa value of metHb A reconstituted with mesohaem We also determined the pKa value of metHb A reconstituted with mesohaem [metHb A(Meso)]. Fitting of the pH-dependent 575-nm absorption to the Henderson–Hasselbalch equation yielded the apparent pKa value of 8.50 ± 0.05, and, in contrast to the case of native metHb A, the use of two Henderson–Hasselbalch equations for the fitting did not improve the agreement between the observed and theoretical values (Supplementary Material SM7). These results indicated that the pKa values of both the α and β subunits of metHb A(Meso) are ∼8.5. Considering that the pKa values for the α and β subunits of native metHb A are 8.32 ± 0.05 and 7.48 ± 0.05, respectively, the substitution of haem peripheral vinyl side chains by ethyl groups, as a result of the haem replacement, was found to reduce the non-equivalence of the pKa values between the α and β subunits of metHb A, and also to increase the pKa values of the α and β subunits by ∼0.3 and 1.1 pH units, respectively. Similarly, 8.50 ± 0.05 was determined for the pKa value of metHb F reconstituted with mesohaem [metHb F(Meso)] (Supplementary Material SM7). Consequently, as in the case of Hb A, the replacement of protohaem by mesohaem resulted in the loss of the non-equivalence in the pKa value between the subunits of metHb F, and increases in the pKa values of the α and γ subunits by ∼1.0 and 0.7 pH units, respectively. Determination of the pKa values of the subunits of metHb A and metHb F reconstituted with 7-PF 19F NMR spectra of metHb A and metHb F reconstituted with 7-PF [metHb A(7-PF) and metHb F(7-PF), respectively] at 25°C and basic pH, together with that of metMb(7-PF) at pH 10.91, are shown in Fig. 6A. The line widths of the Hb signals, i.e. 600–800 Hz, were larger than those of the Mb ones, i.e. ∼500 Hz, possibly due to the difference in size between the proteins. Two signals due to the haem orientational disorder were observed in the spectrum of metMb(7-PF), the ratio being 1.0 : 2.2 for N form: R form at equilibrium (Fig. 6A) (31). On the other hand, in the spectra of metHb A(7-PF) and metHb F(7-PF), three signals, peaks ai–ci (i = A or F), were resolved. These signals could be assigned through analysis of their intensities, and spectral comparison between metHb A(7-PF) and metHb F(7-PF). First of all, the sum of the intensities of peaks ai and bi was almost the same as that of peak ci, indicating that the former signals arise from one subunit and hence the latter ones from the other. In addition, the shifts of peaks aA and bA, i.e. 71.2 and 63.0 ppm, respectively, are similar to those of peaks aF and bF, i.e. 71.0 and 62.9 ppm, respectively, suggesting that peaks ai and bi arise from the α subunit, which is common to both Hb A and Hb F. Hence, peaks cA and cF observed at 51.3 and 53.1 ppm, respectively, are assigned to the β and γ subunits of the proteins, respectively. Consequently, the shift differences between the two signals due to the haem orientational disorder were found to be larger for Hbs, i.e. 8.06 and 9.23 ppm for metHb A(7-PF) and metHb F(7-PF), respectively, than metMb(7-PF), i.e. 2.30 ppm. Furthermore, assuming that the R form dominates over the N one, as in the case of metMb(7-PF), peaks ai and bi are tentatively assigned to the N and R forms, respectively. Thus, peaks aA and bA could be assigned to the N and R forms of the α subunit of metHb A(7-PF), respectively [αA(N) and αA(R), respectively], and aF and bF to the N and R forms of the α subunit of metHb F(7-PF), respectively [αF(N) and αF(R), respectively]. Analysis of the signal intensities revealed that the ratios for N form: R form for the α subunits of metHb A(7-PF) and metHb F(7-PF) are ∼1 : 3, and, on the other hand, the haem orientational disorder in both the β subunit of metHb A(7-PF) and the γ subunit of metHb F(7-PF) is negligible, if any. Fig. 6 View largeDownload slide 470 MHz 19F NMR spectra of metHb A(7-PF). (A) 470 MHz 19F NMR spectra of metHb A(7-PF) at pH 9.0, metHb F(7-PF) at pH 9.0, and metMb(7-PF) at pH 10.91. Three peaks, ai–ci (i = A or F), were observed in the spectrum of each of metHb A(7-PF) and metHb F(7-PF). In the spectrum of metHb A(7-PF), peaks aA and bA are assigned to the N and R forms of the α subunit, and cA to either the N or R form of the β subunit (see text). Similarly, in the spectrum of metHb F(7-PF), peaks aF and bF are assigned to the N and R forms of the α subunit, and cF to either the N or R form of the γ subunit (see text). The assignments of two signals, due to the haem orientational disorder, in the spectrum of metMb(7-PF)14 are indicated in the spectrum. Structure of 7-PF is also illustrated. (B) The spectra of metHb A(7-PF) at 25°C and the indicated pH values. Fig. 6 View largeDownload slide 470 MHz 19F NMR spectra of metHb A(7-PF). (A) 470 MHz 19F NMR spectra of metHb A(7-PF) at pH 9.0, metHb F(7-PF) at pH 9.0, and metMb(7-PF) at pH 10.91. Three peaks, ai–ci (i = A or F), were observed in the spectrum of each of metHb A(7-PF) and metHb F(7-PF). In the spectrum of metHb A(7-PF), peaks aA and bA are assigned to the N and R forms of the α subunit, and cA to either the N or R form of the β subunit (see text). Similarly, in the spectrum of metHb F(7-PF), peaks aF and bF are assigned to the N and R forms of the α subunit, and cF to either the N or R form of the γ subunit (see text). The assignments of two signals, due to the haem orientational disorder, in the spectrum of metMb(7-PF)14 are indicated in the spectrum. Structure of 7-PF is also illustrated. (B) The spectra of metHb A(7-PF) at 25°C and the indicated pH values. With decreasing pH, the signals exhibited progressive downfield shifts of ∼40–50 ppm. The acid–alkaline transition in the protein was manifested in the large pH-induced shift changes attributable to a change in the spin state between the acidic (essentially S = 5/2) and alkaline (mainly S = 1/2) forms (Fig. 6B). Quantitative fitting of their pH-dependent shifts to the Henderson–Hasselbach equation yielded pKa values of 7.55 ± 0.05, 7.59 ± 0.05 and 7.38 ± 0.05 for the N and R forms of the α subunit, and the β subunit of metHb A(7-PF), respectively (Fig. 7). On the other hand, an optical study of metHb A(7-PF) yielded a pKa value of 7.52 ± 0.05 (see Supplementary Material SM8). Similarly, pKa values of 7.50 ± 0.05, 7.60 ± 0.05 and 7.61 ± 0.05 were obtained for the N and R forms of the α subunit, and the γ subunit of metHb F(7-PF), respectively (Fig. 7). The obtained pKa values are summarized in Table III. Fig. 7 View largeDownload slide Plots of the shifts of the 19F NMR signals of metHb A(7-PF) (top) and metHb F(7-PF) (bottom) at 25°C as a function of pH. pKa values of 7.55 ± 0.05, 7.59 ± 0.05 and 7.38 ± 0.05 were obtained for the N and R forms of the α subunit, and the β subunit of metHb A(7-PF), respectively; and those of 7.50 ± 0.05, 7.60 ± 0.05 and 7.61 ± 0.05 for the N and R forms of the α subunit, and the γ subunit of metHb F(7-PF), respectively. Fig. 7 View largeDownload slide Plots of the shifts of the 19F NMR signals of metHb A(7-PF) (top) and metHb F(7-PF) (bottom) at 25°C as a function of pH. pKa values of 7.55 ± 0.05, 7.59 ± 0.05 and 7.38 ± 0.05 were obtained for the N and R forms of the α subunit, and the β subunit of metHb A(7-PF), respectively; and those of 7.50 ± 0.05, 7.60 ± 0.05 and 7.61 ± 0.05 for the N and R forms of the α subunit, and the γ subunit of metHb F(7-PF), respectively. DFT calculations We have recently shown that the pKa value is closely related to the ρFe value which is determined by electronic nature of the haem peripheral side chains (33). Since the electronic effect of the haem vinyl side chains on the porphyrin π-system is though to be affected by the vinyl side chain conformation, we have carried out DFT calculations using a model compound (see ‘Materials and Methods’ section, and Supplementary Material SM2 for details) in order to quantitatively assess the relationship between the Ψ3-vinyl and Ψ8-vinyl values and the ρFe value. The ρFe value was estimated on the basis of the Mulliken charge of the haem Fe atom. The obtained Mulliken charge decreased with increasing the Ψ3-vinyl and Ψ8-vinyl values from 0° to 90°, and the Ψ3-vinyl and Ψ8-vinyl values were found to alter the Mulliken charge in an additive manner (Fig. 3). Although the Mulliken charges were calculated for the ferrous model compounds, a similar relationship between the ρFe value and the Ψ3-vinyl and Ψ8-vinyl values is also expected for the ferrihaems in the active sites of the proteins. Discussion pKa values of the subunits of metHb A and metHb F We first compared the pKa values of the subunits of metHb A and metHb F in order to determine the non-equivalence in the haem environment among them. As shown in Table I, there is an appreciable difference in the pKa value of the α subunit between metHb A and metHb F. This finding indicated that the structural properties of the haem active site in the α subunit are influenced by interactions between the constituent subunits of the protein (‘inter-subunit interactions’). The pKa value is affected by haem environment and the electronic nature of the haem peripheral side chains. Ferrihaem in metHb carries a net positive-charge, and hence needs to be stabilized by neutralization through electrostatic interaction with near-by polar groups in a hydrophobic haem environment. Consequently, as far as the stability of the ferrihaem in the protein is concerned, the alkaline form is more stable than the acidic one because of neutralization of the cationic character of ferrihaem by the coordinated OH− (Scheme 1). On the other hand, as shown in Fig. 3, the Mulliken charge is affected by the Ψ3-vinyl and Ψ8-vinyl values. Since a decrease in the Mulliken charge is thought to represent an increase in the ρFe value, the plots in Fig. 3 demonstrated that the ρFe value increases with increasing the Ψ3-vinyl and/or Ψ8-vinyl values. The pKa value has been shown to be closely related to the ρFe value (33). The backward reaction of the acid–alkaline transition (Scheme 1) can be considered as protonation of the Fe-bound OH−. Since a decrease in the ρFe value results in removal of an electron from the Fe-bound OH−, the proton affinity of the Fe-bound OH− is thought to decrease with decreasing the ρFe value, leading to a decrease in the pKa value. Hence the pKa value decreases with decreasing the Ψ3-vinyl and/or Ψ8-vinyl values. X-ray structure of Hb F is available only for deoxy form, and X-ray structural comparison between the α subunits of deoxy Hb A and Hb F indicated that there are sizable differences in the Ψ3-vinyl and Ψ8-vinyl values between them (Table IV) (3, 46, 49, 60), although their protein structures are highly alike (61, 46). Assuming that the larger Ψ3-vinyl and Ψ8-vinyl values for the α subunit of Hb F relative to the corresponding values of Hb A, as observed in the X-ray structures of their deoxy forms, are retained in their met-forms, the pKa value of the α subunit of metHb F is expected to be higher than that of metHb A. Thus the higher pKa value of the α subunit of metHb F than that of metHb A could be attributed to the differences in the Ψ3-vinyl and Ψ8-vinyl values between the α subunits of the two metHbs. Table III. The pKa values of individual subunits of native metHb A, metHb F and the proteins reconstituted with meso and 7-PF at 25°C.   metHb A   metHb F       7-PF       7-PF   protoa  meso  Nb  Rb  protoa  meso  Nb  Rb  α  8.32 ± 0.05  8.52 ± 0.05  7.55 ± 0.05  7.59 ± 0.05  8.48 ± 0.05  8.50 ± 0.05  7.50 ± 0.05  7.60 ± 0.05  β  7.48 ± 0.05  8.52 ± 0.05  7.38 ± 0.05c  –  –  –  –  γ  –  –  –  –  7.76 ± 0.05  8.50 ± 0.05  7.61 ± 0.05c    metHb A   metHb F       7-PF       7-PF   protoa  meso  Nb  Rb  protoa  meso  Nb  Rb  α  8.32 ± 0.05  8.52 ± 0.05  7.55 ± 0.05  7.59 ± 0.05  8.48 ± 0.05  8.50 ± 0.05  7.50 ± 0.05  7.60 ± 0.05  β  7.48 ± 0.05  8.52 ± 0.05  7.38 ± 0.05c  –  –  –  –  γ  –  –  –  –  7.76 ± 0.05  8.50 ± 0.05  7.61 ± 0.05c  aTaken from native proteins. bN and R represent the normal and reversed haem orientations, respectively (Fig. 1B). cThis subunit predominantly possesses a single heme orientation and the determination of the haem orientation could not be made from the present study (see text). View Large On the other hand, the difference in the pKa value between the subunits of metHb A could be interpreted in terms of their haem environments, because the Ψ3-vinyl and Ψ8-vinyl values of the α and β subunits of metHb A are almost identical to each other (Table IV). In the case of metHb A, the pKa value of the α subunit is higher by ∼0.8 pH units than that of the β subunit. This finding could be attributed to a more polar haem environment in the α subunit than the latter. This interpretation is supported from the study of autoxidation of oxy Hb A. Yasuda et al. (62) revealed that the autoxidation rate for the α subunit is larger than that for the β one. The autoxidation rate is expected to increase with increasing polarity of the haem pocket. Therefore, the higher polarity of the haem environment in the α subunit than that in the β one, as manifested in the larger autoxidation rate, is consistent with our interpretation, although the haem Fe oxidation and ligation states of the two systems are different from each other. Finally, considering the amino acid sequence homology between the β and γ subunits (see Supplementary Material SM1), the higher pKa value for the α subunit than the γ one in metHb F might also be attributed to a more polar haem environment in the former than the latter (see below). Table IV. Ψ3-vinyl and Ψ8-vinyl values of Hbs and isolated subunits.   |Ψ3-vinyl|a(degree)   |Ψ8-vinyl|a(degree)   Protein  α  β  γ  α  β  γ  Deoxy form               Hb Ab  33.2  32.0  –  31.0  32.8  –  36.5  32.8  37.5  33.1      Hb Fc  69.2  –  81.4  81.5  –  15.7  CO form               Hb Ab  32.7  6.0  –  38.3  23.9  –  30.5  26.8   β subunitd  –  48.2  –  –  30.8  –  61.4  31.3  71.6  31.6   γ subunite  –  –  49.3  –  –  39.6  50.8  60.2    |Ψ3-vinyl|a(degree)   |Ψ8-vinyl|a(degree)   Protein  α  β  γ  α  β  γ  Deoxy form               Hb Ab  33.2  32.0  –  31.0  32.8  –  36.5  32.8  37.5  33.1      Hb Fc  69.2  –  81.4  81.5  –  15.7  CO form               Hb Ab  32.7  6.0  –  38.3  23.9  –  30.5  26.8   β subunitd  –  48.2  –  –  30.8  –  61.4  31.3  71.6  31.6   γ subunite  –  –  49.3  –  –  39.6  50.8  60.2  aThe Ψ3-vinyl and Ψ8-vinyl values are indicated in the absolute values. In addition, since the orientation of the vinyl side chain π system with respect to the porphyrin π one is considered in the study, the values are indicated in the range of 0°–90°. Multiple entry represents disorder in the conformation of the vinyl side chains. bObtained from ref. 3. cObtained from ref. 46. dObtained from ref. 49. eObtained from ref. 60. View Large Effects of hetero-tetramer assembly on the pKa values of the subunits We next compared the pKa values of the isolated subunits and metHbs in order to elucidate the effect of the hetero-tetramer assembly on the haem environments in the individual subunits. As indicated in Table I, the pKa values of the isolated α, β and γ subunits were different from each other. Interestingly, the value of the isolated γ subunit was closer to that of the isolated α subunit than that of the isolated β subunit, although the β and γ subunits exhibit high sequence homology of 73% (51). This finding clearly indicated that the pKa value is determined solely by the haem local environment. The difference in the pKa values between the isolated β and γ subunits may be interpreted in terms of amino acid substitutions in them (Supplementary Material SM1). Among the 39 amino acid residue substitutions in the two subunits, two replacement sites, i.e. E14 and E15, are located in close proximity to the haem, and Ala E14 and Phe E15 in the β subunit are replaced by Ser and Leu in the γ subunit, respectively. Although E14 and E15 residues are located at least ∼0.8 nm away from haem iron, the replacements of the amino acid residues at these sites have been shown to alter the haem chemical environment largely enough to be manifested in the shifts of haem methyl proton NMR signals (63). Therefore, with these replacements, the polarity of the haem active site in the γ subunit is assumed to be higher than that in the β one, and hence the pKa value of the former is expected to be higher than that of the latter. The pKa value of the γ subunit was indeed larger than that of the β one (Table I). Furthermore, the λmax value of Soret absorption of the isolated γ subunit was larger than that of the isolated β one (Table II), also supporting that the polarity of the haem environment in the γ subunit is higher than that in the β one, because the Soret absorption has been shown to exhibit bathochromic shift with increasing the polarity of solvent (59). In addition, the increase and decrease in the polarity of the haem environments of the α subunit, and the β and γ subunits, respectively, upon the hetero-tetramer assembly, as suggested from the analysis of the λmax values (Table II), are consistent with the observation that, for the α subunit, the pKa value of the isolated state is smaller than that of the Hb state, while, for the β and γ subunis, the values of the isolated states were larger than those of the Hb states (Table I). The pKa values of the α subunits in metHb A and metHb F were higher by 0.15 and 0.31 pH unit, respectively, relative to that of the isolated α subunit. On the other hand, the values of the β and γ subunits in metHb A and metHb F, respectively, were lower by 0.27 and 0.50 pH units, respectively, relative to those of the corresponding isolated subunits. As a result, the difference in the pKa values between the α and β subunits of metHb A [ΔpKa(A)] is larger than that between the isolated α and β subunits [ΔpKa(α/β)], and similarly the difference in the pKa values between the α and γ subunits of metHb F [ΔpKa(F)] is larger than that between the isolated α and γ subunits [ΔpKa(α/γ)]. Thus, the non-equivalence in the pKa value between the constituent subunits, and hence the difference in the haem local environment between them, was found to be enhanced by the hetero-tetramer assembly. Such an asymmetric nature of active sites might be essential for allosteric proteins composed of multiple subunits. Effects of hetero-tetramer assembly on the pKa values of the constituent subunits of valency hybrid Hbs The effect of the hetero-tetramer assembly on the haem local environment was also reflected in the pKa values of valency hybrid proteins. For all the subunits, the valency hybrid proteins exhibited a pKa value between those of the corresponding subunits in the isolated and tetrameric Hb states, i.e. for the α subunits of both metHb A and metHb F, the pKa value of the valency hybrid state is larger and smaller than those of the isolated and Hb states, respectively, and, on the other hand, for the β and γ subunits of metHb A and metHb F, the value of the valency hybrid state is smaller and larger than those of the isolated and Hb states, respectively. These results demonstrated that the pKa value of a subunit is affected by not only the hetero-tetramer assembly, but also the haem Fe oxidation and ligation states of the partner subunit. The subunit interactions in Hb A and Hb F enable communication between the haem active sites of the constituent subunits of the proteins. Hence the molecular mechanisms responsible for the subunit interactions in the protein are likely to involve structural changes in not only the sites of inter-subunits interaction, but also the contact between haem and the protein at the active site (‘haem–protein interactions’). As described above, the effect of the hetero-tetramer assembly on the haem environment was clearly manifested in the pKa values of the subunits. Furthermore, the effect of the haem–protein interactions on the inter-subunit interaction was also manifested in the differences in the pKa values of the subunits between the valency hybrid and Hb states. These results confirmed the presence of a pathway for transmission of a structural change of one subunit to the other ones within tetrameric Hb, although its molecular mechanism is remained to be elucidated. Haem active site structures of metHb A(7-PF) and metHb F(7-PF) Because of its wide spectral range, 19F NMR is remarkably sensitive as to the local magnetic environment, exhibiting not only extremely high resolution of signals, but also a large dynamic range (41). Well-resolved signals were observed in the 19F NMR spectra of metHb A(7-PF) and metHb F(7-PF) (Fig. 6A). The shift differences between the signals due to the haem orientational disorder in the α subunits of metHb A(7-PF) and metHb F(7-PF) are similar to each other, and are larger than that in the case of metMb(7-PF). The shifts of the signals have been shown to sharply reflect the unpaired electron delocalization of the haem, which is predominantly determined by the nature of the axial ligands (64, 65). Consequently, the present results suggested that the Φ angle between the projection of the axial His imidazole onto the haem plane and the NII-Fe-NIV axis (Fig. 1B) of the α subunits of metHb A(7-PF) and metHb F(7-PF) are similar to each other, and are smaller than that in the case of metMb(7-PF), although the X-ray crystal structures of the native proteins indicated that the Φ angles of the α subunits of native Hb A and Hb F, and met Mb are similar to each other (3, 47). Haem Orientational Disorder in metHb A(7-PF) and metHb F(7-PF) Analysis of the signal intensities in 19F NMR spectra revealed that the N form: R form ratios for the α subunits of metHb A(7-PF) and metHb F(7-PF) are ∼1 : 3, and essentially no haem orientational disorder was detected in the β and γ subunits of either protein. These ratios were different from that for metMb(7-PF), i.e. N form: R form = 1.0 : 2.2 (31). Considering the symmetry of the 7-PF molecular structure, it is obvious that steric contact between the haem peripheral side chains attached at positions 2, 3, 7 and 8, and near-by amino acid residues is crucial for the thermodynamic stabilities of the N and R forms. The haem side chains at positions 2, 3, 7 and 8 are in close contact with the amino acid residues at E10:E14:F4, E15:G5:G8:G12:H15, B13:FG5:G4:G8 and B13:CD1:C4:C7:FG5:G4, respectively (3, 46). Among these residues, amino acid replacements were found at the B13 and C7 positions between the α and β (or γ) subunits, i.e. Met B13 and Tyr C7 in the α subunit are replaced by Leu and Phe in the β (or γ) subunit, respectively (Supplementary Material SM1) (3, 46). Furthermore, in the case of Mb, Leu and Lys occupy B13 and C7, respectively (see Supplementary Material SM1) (2). Consequently, the amino acid residue at C7 may play a role in determination of the thermodynamics of the haem orientational disorder of 7-PF in the subunits of Hbs and Mb. According to the X-ray crystal structures of native Hb A and Mb (3, 66), the side chains of Tyr C7 and Lys C7 in the α subunit and Mb, respectively, point somewhat away from the haem, and are hydrogen bonded to the side chain NH2 group of Asn G4 and carbonyl oxygen of Lys FG4, respectively. Hence the contact between the haem side chain at position 8 and the C7 residue in the α subunit and Mb would not be so tight. On the other hand, the side chain of Phe C7 of the β subunit is in close contact with haem, leading to tight contact with the haem side chain at position 8. Such steric hindrance might render the energy difference between the two haem orientations, i.e. the N and R forms, large enough to almost completely inhibit one of the orientations. pKa values of reconstituted metHbs We finally analysed the effect of the haem substitution on the pKa values of the subunits of metHbs in order to determine the effect of alteration of the haem–protein interaction on the haem environments of the proteins. In the case of the metHb A system, both the α and β subunits of metHb A(Meso) exhibited pKa values of ∼8.5, and values of 7.55 ± 0.05, 7.59 ± 0.05 and 7.38 ± 0.05 were obtained for the N and R forms of the α subunit, and the β subunit of metHb A(7-PF), respectively (Table I). Since mesohaem and 7-PF differ from each other only in the substituent at position 7 (Fig. 1A), the decrease in the pKa value by ∼1 pH unit upon the haem replacement from the former to the latter demonstrated the significant effect of the CF3 substitution on the acid–alkaline transition. The substitution of the electron-withdrawing CF3 group for the haem peripheral CH3 side chain decreases the the ρFe value, leading to a decrease in the pKa value (33). The haem substitution also influenced the non-equivalence of the pKa values of the constituent subunits of Hbs. In the case of metHb A, the ΔpKa(A) value was significantly decreased by the haem substitution, i.e. the ΔpKa(A) values were ∼0.8, 0 and 0.2 pH unit for metHb A, metHb A(Meso) and metHb A(7-PF), respectively (Table III). Similarly, since the values of 7.50 ± 0.05, 7.60 ± 0.05 and 7.61 ± 0.05 were obtained for the N and R forms of the α subunit, and the γ subunit of metHb F(7-PF), respectively (Table III), the ΔpKa(F) value was also decreased by the haem substitution, i.e. the ΔpKa(F) values were ∼0.7 and 0 pH units for metHb F and metHb F(7-PF), respectively. Thus, the substitution of the native haem, i.e. protohaem, by mesohaem or 7-PF was found to reduce the non-equivalence of the haem environment of the constituent subunits of the proteins. As described earlier, since the pKa values of the subunits of the proteins have been shown to be affected by the inter-subunit interactions and the haem–protein interaction, the decrease in the ΔpKa(A) [or ΔpKa(F)] value on the haem replacement could be due to alteration of the haem–protein interaction. The substitution of the haem peripheral vinyl side chains by ethyl groups, as a result of the haem replacement from protohaem to mesohaem (or 7-PF), should significantly alter the haem–protein interaction. Furthermore, the cooperative ligand binding of Hb A has been shown to be affected by the haem–protein interaction, especially the contacts between the haem peripheral side chains attached to pyrroles I and II, and the protein (67). Consequently, the present results demonstrated that the interaction between the haem vinyl side chain and protein contributes significantly to maintain not only the non-equivalence of the haem active sites of the constituent subunits of the proteins, but also the inter-subunit interaction essential for the cooperative oxygen binding of the protein. In fact, the oxygen binding study on Hb A demonstrated that the cooperativity, as measured as Hill’s n-value, decreased from 3.2 to 1.6 with the replacement of protohaem by mesohaem (68), indicating the importance of the haem vinyl groups for the subunit interaction. Funding Grant-in-aid for Scientific Research on Innovative Areas (No. 21108505, ‘π-Space’) from the Ministry of Education, Culture, Sports, Science and Technology, Japan; the Yazaki Memorial Foundation for Science and Technology, and the NOVARTIS Foundation (Japan) for the Promotion of Science. Conflict of interest None declared. Acknowledgements The 19F NMR spectra were recorded on a Bruker AVANCE-500 spectrometer at the Chemical Analysis Center, University of Tsukuba. Abbreviations Abbreviations 7-PF 13,17-bis(2-carboxylatoethyl)-3,8-diethyl-2,12,18-trimethyl-7-trifluoromethyl-porphyrinatoiron Bis-Tris Bis(2-hydroxyethyl)iminotris(hydroxymethyl)methane DFT density functional theory Hb A human adult haemoglobin Hb F human foetal haemoglobin Mb myoglobin mesohaem 13,17-bis(2-carboxylatoethyl)-3,8-diethyl-2,7,12,18-tetramethylporphyrinatoiron meso mesohaem metHb A met-form of human adult haemoglobin metHb F met-form of human foetal haemoglobin protohaem Fe-protoporphyrin IX complex TFA trifluoroacetic acid References 1 Perutz MF,  Wilkinson AJ,  Paoli M,  Dodson GG.  The stereochemical mechanism of the cooperative effects in hemoglobin revisited,  Annu. Rev. Biophys. Biomol. Struct. ,  1998, vol.  27 (pg.  1- 34) Google Scholar CrossRef Search ADS PubMed  2 Dickerson RE,  Geis I. ,  Hemoglobin: Structure, Function, Evolution, and Pathology ,  1983 Menlo Park, CA Benjamin/Cummings 3 Park SY,  Yokoyama T,  Shibayama N,  Shiro Y,  Tame J.RH.  1.25 Å resolution crystal structures of human haemoglobin in the oxy, deoxy and carbonmonoxy forms,  J. Mol. Biol. ,  2006, vol.  360 (pg.  690- 701) Google Scholar CrossRef Search ADS PubMed  4 Imai K.  Antonini E,  Rossi-Bernardi L,  Chiancone E.  Measurement of accurate oxygen equilibrium curves by an automatic oxygenation apparatus,  Methods in Enzymology ,  1981, vol.  Vol. 76  New York Academic Press(pg.  438- 449) 5 Antonini E,  Brunori M. ,  Hemoglobins and Myoglobins in their Reactions with Ligands ,  1971, vol.  Chapters 2 and 3  Amsterdam North Holland Publishing 6 Blough NV,  Zemel H,  Hoffman BM,  Lee T.CK,  Gibson QH.  Kinetics of CO binding to manganese, zinc, and cobalt hybrid hemoglobins,  J. Am. Chem. Soc. ,  1980, vol.  102 (pg.  5683- 5685) Google Scholar CrossRef Search ADS   7 Inubushi T,  Ikeda-Saito M,  Yonetani T.  Isotropcally shifted NMR resonances for the proximal histidyl imidazole NH protons in cobalt hemoglobin and iron-cobalt hybrid hemoglobins. Binding of the proximal histidine toward porphyrin metal ion in the intermediate state of cooperative ligand binding,  Biochemistry ,  1983, vol.  22 (pg.  2904- 2907) Google Scholar CrossRef Search ADS PubMed  8 Blough NV,  Hoffman BM.  Carbon monoxide binding to the ferrous chains of [Mn, Fe(II)] hybrid hemoglobins: pH dependence of the chain affinity constants associated with specific hemoglobin ligation pathways,  Biochemistry ,  1984, vol.  23 (pg.  2875- 2882) Google Scholar CrossRef Search ADS PubMed  9 Shibayama N,  Morimoto H,  Miyazaki G.  Oxygen equilibrium study and light absorption spectra of Ni(II)-Fe(II) hybrid hemoglobins,  J. Mol. Biol. ,  1986, vol.  192 (pg.  323- 329) Google Scholar CrossRef Search ADS PubMed  10 Inubushi T,  D’Ambtosio C,  Ileda-Saito M,  Yonetani T.  NMR studies of monoliganded Fe-Co hybrid hemoglobins: their quaternary structure and proximal histidine coordination,  J. Am. Chem. Soc. ,  1986, vol.  108 (pg.  3799- 3803) Google Scholar CrossRef Search ADS   11 Shibayama N,  Ikeda-Saito M,  Hori H,  Itaroku K,  Morimoto H,  Saigo S.  Oxygen equilibrium and electron paramagnetic resonance studies on copper(II)-iron(II) hybrid hemoglobins at room temperature,  FEBS Lett. ,  1995, vol.  372 (pg.  126- 130) Google Scholar CrossRef Search ADS PubMed  12 Bowen JH,  Shokhirev NV,  Raitsimring AM,  Buttlaire DH,  Walker FA.  EPR studies of the dynamics of rotation of dioxygen in model cobalt(II) hemes and cobalt-containing hybrid hemoglobins,  J. Phys. Chem. B ,  1997, vol.  101 (pg.  8683- 8691) Google Scholar CrossRef Search ADS   13 Unzai S,  Eich R,  Shibayama N,  Olson JS,  Morimoto H.  Rate constants for O2 and CO binding to the α and β subunits within the R and T states of human hemoglobin,  J. Biol. Chem. ,  1998, vol.  273 (pg.  23150- 23159) Google Scholar CrossRef Search ADS PubMed  14 Miyazaki G,  Morimoto H,  Yun K.-M,  Park S.-Y,  Nakagawa A,  Minagawa H,  Shibayama M.  Magnesium(II) and zinc(II)-protoporphyrin IX’s stabilize the lowest oxygen affinity state of human hemoglobin even more strongly than deoxyheme,  J. Mol. Biol. ,  1999, vol.  292 (pg.  1121- 1136) Google Scholar CrossRef Search ADS PubMed  15 Nagatomo S,  Nagai M,  Shibayama N,  Kitagawa T.  Differences in changes of the α1-β2 subunit contacts between ligand binding to the α and β subunits of hemoglobin A: UV resonance Raman analysis using Ni-Fe hybrid hemoglobin,  Biochemistry ,  2002, vol.  41 (pg.  10010- 10020) Google Scholar CrossRef Search ADS PubMed  16 Samuni U,  Juszczak L,  Dantsker D,  Khan I,  Friedman AJ,  Pérez-González-de-Apodaca J,  Bruno S,  Hui HL,  Colby JE,  Karasik E,  Kwiatkowski LD,  Mozzarelli A,  Noble R,  Friedman JM.  Functional and spectroscopic characterization of half-liganded iron−zinc hybrid hemoglobin: evidence for conformational plasticity within the T state,  Biochemistry ,  2003, vol.  42 (pg.  8272- 8288) Google Scholar CrossRef Search ADS PubMed  17 Simplaceanu V,  Lukin JA,  Fang T.-Y,  Zou M,  Ho NT,  Ho C.  Chain-selective isotopic labeling for NMR studies of large multimeric proteins: application to hemoglobin,  Biophys. J. ,  2000, vol.  79 (pg.  1146- 1154) Google Scholar CrossRef Search ADS PubMed  18 Song X.-J,  Simplaceanu V,  Lukin JA,  Ho NT,  Ho C.  Effector-induced structural fluctuation regulates the ligand affinity of an allosteric protein: binding of inositol hexaphosphate has distinct dynamic consequences for the T and R states of hemoglobin,  Biochemistry ,  2008, vol.  47 (pg.  4907- 4915) Google Scholar CrossRef Search ADS PubMed  19 Hu X,  Spiro TG.  Tyrosine and tryptophan structure markers in hemoglobin ultraviolet resonance Raman spectra: mode assignments via subunit-specific isotope labeling of recombinant protein,  Biochemistry ,  1997, vol.  36 (pg.  15701- 15712) Google Scholar CrossRef Search ADS PubMed  20 Balakrishnan G,  Zhao X,  Podstawska E,  Proniewicz LM,  Kincaid JR,  Spiro TG.  Subunit-selective interrogation of CO recombination in carbonmonoxy hemoglobin by isotope-edited time-resolved resonance Raman spectroscopy,  Biochemistry ,  2009, vol.  48 (pg.  3120- 3126) Google Scholar CrossRef Search ADS PubMed  21 Nagai K,  Kitagawa T.  Differences in Fe(II)-Nε(His-F8) stretching frequencies between deoxyhemoglobins in the two alternative quaternary structures,  Proc. Natl Acad. Sci. USA ,  1980, vol.  77 (pg.  2033- 2037) Google Scholar CrossRef Search ADS   22 Jue T,  La Mar GN,  Han K.-H,  Yamamoto Y.  NMR study of the exchange rates of allosterically responsive labile protons in the heme pockets of hemoglobin A,  Biophys. J. ,  1984, vol.  46 (pg.  117- 120) Google Scholar CrossRef Search ADS PubMed  23 Han K.-H,  La Mar GN.  Nuclear magnetic resonance study of the isotope exchange of the proximal histidyl ring labile protons in hemoglobin A: the exchange rates and mechanisms of individual subunits in deoxy and oxy-hemoglobin,  J. Mol. Biol. ,  1986, vol.  189 (pg.  541- 552) Google Scholar CrossRef Search ADS PubMed  24 Giacometti GM,  Da Ros A,  Antonini E,  Brunori M.  Equilibrium and kinetics of the reaction of Aplysia myoglobin with azide,  Biochemistry ,  1975, vol.  14 (pg.  1584- 1588) Google Scholar CrossRef Search ADS PubMed  25 Iizuka T,  Morishima I.  NMR studies of hemoproteins. VI. Acid-base transitions of ferric myoglobin and its imidazole complex,  Biochim. Biophys. Acta ,  1975, vol.  400 (pg.  143- 153) Google Scholar CrossRef Search ADS PubMed  26 McGrath TM,  La Mar GN.  Proton NMR study of the thermodynamics and kinetics of the acid in equilibrium base transitions in reconstituted metmyoglobins,  Biochim. Biophys. Acta ,  1978, vol.  534 (pg.  99- 111) Google Scholar CrossRef Search ADS PubMed  27 Pande U,  La Mar GN,  Lecomte J.TL,  Ascoli F,  Brunori M,  Smith KM,  Pandey RK,  Parish DW,  Thanabal V.  NMR study of the molecular and electronic structure of the heme cavity of Aplysia metmyoglobin. Resonance assignments based on isotope labeling and proton nuclear Overhauser effect measurements,  Biochemistry ,  1986, vol.  25 (pg.  5638- 5646) Google Scholar CrossRef Search ADS PubMed  28 Yamamoto Y,  Osawa A,  Inoue Y,  Chûjô R,  Suzuki T.  A 1H-NMR study of electronic structure of the active site of Galeorhinus japonicus metmyoglobin,  Eur. J. Biochem. ,  1990, vol.  192 (pg.  225- 229) Google Scholar CrossRef Search ADS PubMed  29 Yamamoto Y,  Chûjô R,  Inoue Y,  Suzuki T.  Kinetic characterization of the acid-alkaline transition in Dolabella auricularia ferric myoglobin using 1H-NMR saturation transfer experiments,  FEBS Lett. ,  1992, vol.  310 (pg.  71- 74) Google Scholar CrossRef Search ADS PubMed  30 Yamamoto Y,  Suzuki T,  Hori H.  Dynamics and thermodynamics of acid-alkaline transitions in metmyoglobins lacking the usual distal histidine residue,  Biochim. Biophys. Acta ,  1993, vol.  1203 (pg.  267- 275) Google Scholar CrossRef Search ADS PubMed  31 Nagao S,  Hirai Y,  Suzuki A,  Yamamoto Y.  19F NMR characterization of the thermodynamics and dynamics of the acid-alkaline transition in a reconstituted sperm whale metmyoglobin,  J. Am. Chem. Soc. ,  2005, vol.  127 (pg.  4146- 4147) Google Scholar CrossRef Search ADS PubMed  32 Shu F,  Ramakrishnan V,  Schoenborn BP.  Enhanced visibility of hydrogen atoms by neutron crystallography on fully deuterated myoglobin,  Proc. Natl Acad. Sci. USA ,  2000, vol.  97 (pg.  3872- 3877) Google Scholar CrossRef Search ADS   33 Shibata T,  Nagao S,  Fukaya M,  Tai H,  Nagatomo S,  Morihashi K,  Matsuo T,  Hirota S,  Suzuki A,  Imai K,  Yamamoto Y.  Effect of heme modification on oxygen affinity of myoglobin and equilibrium of the acid-alkaline transition in metmyoglobin,  J. Am. Chem. Soc. ,  2010, vol.  132 (pg.  6091- 6098) Google Scholar CrossRef Search ADS PubMed  34 La Mar GN,  Budd DL,  Viscio DB,  Smith KM,  Langry LC.  Proton nuclear magnetic resonance characterization of heme disorder in hemoproteins,  Proc. Natl Acad. Sci. USA ,  1978, vol.  75 (pg.  5755- 5759) Google Scholar CrossRef Search ADS   35 Cassoly R.  Antonini E,  Rossi-Bernardi L,  Chiancone E.  Preparation of hemoglobin hybrid carrying different ligands: valency hybrids and related compounds,  Methods in Enzymology ,  1981, vol.  Vol. 76  New York Academic Press(pg.  106- 113) 36 Toi H,  Homma M,  Suzuki A,  Ogoshi H.  Paramagnetic 19F n.m.r. spectra of iron(III) porphyrins substituted with CF3 groups and reconstituted myoglobin,  J. Chem. Soc. Chem. Commun. ,  1985(pg.  1791- 1792) 37 Yamamoto Y,  Hirai Y,  Suzuki A.  19F NMR study of protein-induced rhombic perturbations on the electronic structure of the active site of myoglobin,  J. Biol. Inorg. Chem. ,  2000, vol.  5 (pg.  455- 462) Google Scholar CrossRef Search ADS PubMed  38 Hirai Y,  Yamamoto Y,  Suzuki A.  19F NMR study of the heme orientation and electronic structure in a myoglobin reconstituted with a ring-fluorinated heme,  Bull. Chem. Soc. Jpn. ,  2000, vol.  73 (pg.  2309- 2316) Google Scholar CrossRef Search ADS   39 Yamamoto Y,  Nagao S,  Hirai Y,  Inose T,  Terui N,  Mita H,  Suzuki A.  NMR investigation of the heme electronic structure in deoxymyoglobin possessing a fluorinated heme,  J. Biol. Inorg. Chem. ,  2004, vol.  9 (pg.  152- 160) Google Scholar CrossRef Search ADS PubMed  40 Hirai Y,  Nagao S,  Mita H,  Suzuki A,  Yamamoto Y.  19F NMR study on the heme electronic structure in oxy and carbonmonoxy reconstituted myoglobins,  Bull. Chem. Soc. Jpn. ,  2004, vol.  77 (pg.  1485- 1486) Google Scholar CrossRef Search ADS   41 Yamamoto Y.  NMR study of active sites in paramagnetic hemoproteins,  Annu. Rep. NMR Spectrosc. ,  1998, vol.  36 (pg.  1- 77) 42 Hull WE,  Sykes BD.  Fluorotyrosine alkaline phosphatase. 19F nuclear magnetic resonance relaxation times and molecular motion of the individual fluorotyrosines,  Biochemistry ,  1974, vol.  13 (pg.  3431- 3437) Google Scholar CrossRef Search ADS PubMed  43 Hull WE,  Sykes BD.  Fluorotyrosine alkaline phosphatase: internal mobility of individual tyrosines and the role of chemical shift anisotropy as a 19F nuclear spin relaxation mechanism in proteins,  J. Mol. Biol. ,  1975, vol.  98 (pg.  121- 153) Google Scholar CrossRef Search ADS PubMed  44 Gueron M.  Nuclear relaxation in macromolecules by paramagnetic ions. Novel mechanism,  J. Magn. Reson. ,  1975, vol.  19 (pg.  58- 66) 45 Vega AJ,  Fiat D.  Nuclear relaxation processes of paramagnetic complexes. The slow-motion case,  Mol. Phys. ,  1976, vol.  31 (pg.  347- 355) Google Scholar CrossRef Search ADS   46 Frier JA,  Perutz MF.  Structure of human foetal deoxyhaemoglobin,  J. Mol. Biol. ,  1977, vol.  112 (pg.  97- 112) Google Scholar CrossRef Search ADS PubMed  47 Valdes R.Jr,  Ackers GK.  Thermodynamic studies on subunit assembly in human hemoglobin. Self-association of oxygenated chains (alphaSH and betaSH): determination of stoichiometries and equilibrium constants as a function of temperature,  J. Biol. Chem. ,  1977, vol.  252 (pg.  74- 81) Google Scholar PubMed  48 Borgstahl G.EO,  Rogers PH,  Arnone A.  The 1.9 Å structure of deoxy beta 4 hemoglobin. Analysis of the partitioning of quaternary-associated and ligand-induced changes in tertiary structure,  J. Mol. Biol. ,  1994, vol.  236 (pg.  831- 843) Google Scholar CrossRef Search ADS PubMed  49 Borgstahl G.EO,  Rogers PH,  Arnone A.  The 1.8 Å structure of carbonmonoxy-beta 4 hemoglobin. Analysis of a homotetramer with the R quaternary structure of liganded alpha 2 beta 2 hemoglobin,  J. Mol. Biol. ,  1994, vol.  236 (pg.  817- 830) Google Scholar CrossRef Search ADS PubMed  50 La Mar GN,  Yamamoto Y,  Jue T,  Smith KM,  Pandey RK.  1H NMR characterization of metastable and equilibrium heme orientational heterogeneity in reconstituted and native human hemoglobin,  Biochemistry ,  1985, vol.  24 (pg.  3826- 3831) Google Scholar CrossRef Search ADS PubMed  51 Prins HK.  Separation of diferent types of human hemoglobin,  J. Chromatogr. ,  1956, vol.  2 (pg.  445- 486) Google Scholar CrossRef Search ADS   52 Bucci E,  Fronticelli C.  A new method for the preparation of alpha and beta subunits of human hemoglobin,  J. Biol. Chem. ,  1965, vol.  240 (pg.  551- 552) 53 Noble RW.  The effect of p-hydroxymercuribenzoate on the reactions of the isolated gamma chains of human hemoglobin with ligands,  J. Biol. Chem. ,  1971, vol.  246 (pg.  2972- 2976) Google Scholar PubMed  54 Tomoda A,  Takeshita M,  Yoneyama Y.  Characterization of intermediate hemoglobin produced during methemoglobin reduction by ascorbic acid,  J. Biol. Chem. ,  1978, vol.  253 (pg.  7415- 7419) Google Scholar PubMed  55 Banerjee R,  Cassoly R.  Oxygen equilibria of human hemoglobin valency hybrids. Discussion on the intrinsic properties of alpha and beta chains in the native protein,  J. Mol. Chem. ,  1969, vol.  42 (pg.  351- 361) 56 Teale F.WJ.  Cleavage of the haem-protein link by acid methylethylketone,  Biochim. Biophys. Acta ,  1959, vol.  35 pg.  543  Google Scholar CrossRef Search ADS PubMed  57 Frisch MJ,  Trucks GW,  Schlegel HB,  Scuseria GE,  Robb MA,  Cheeseman JR,  Montgomery J.A.Jr,  Vreven T,  Kudin KN,  Burant JC,  Millam JM,  Iyengar SS,  Tomasi J,  Barone V,  Mennucci B,  Cossi M,  Scalmani G,  Rega N,  Petersson GA,  Nakatsuji H,  Hada M,  Ehara M,  Toyota K,  Fukuda R,  Hasegawa J,  Ishida M,  Nakajima T,  Honda Y,  Kitao O,  Nakai H,  Klene M,  Li X,  Knox JE,  Hratchian HP,  Cross JB,  Bakken V,  Adamo C,  Jaramillo J,  Gomperts R,  Stratmann RE,  Yazyev O,  Austin AJ,  Cammi R,  Pomelli C,  Ochterski JW,  Ayala PY,  Morokuma K,  Voth GA,  Salvador P,  Dannenberg JJ,  Zakrzewski VG,  Dapprich S,  Daniels AD,  Strain MC,  Farkas O,  Malick DK,  Rabuck AD,  Raghavachari K,  Foresman JB,  Ortiz JV,  Cui Q,  Baboul AG,  Clifford S,  Cioslowski J,  Stefanov BB,  Liu G,  Liashenko A,  Piskorz P,  Komaromi I,  Martin RL,  Fox DJ,  Keith T,  Al-Laham MA,  Peng CY,  Nanayakkara A,  Challacombe M,  Gill P.MW,  Johnson B,  Chen W,  Wong MW,  Gonzalez C,  Pople JA. ,  Gaussian 03 ,  2004 Wallingford, CT Revision E.01, Gaussian, Inc 58 Banerjee R,  Alpert Y,  Leterrier F,  Williams R.JP.  Visible absorption and electron spin resonance spectra of the isolated chains of human hemoglobin. Discussion of chain-mediated heme-heme interaction,  Biochemistry ,  1969, vol.  8 (pg.  2862- 2867) Google Scholar CrossRef Search ADS PubMed  59 Falk JE. ,  Porphyrin and Metalloporphyrins ,  1994 Amsterdam Elsevier 60 Kidd RD,  Baker HM,  Mathews AJ,  Brittain T,  Baker EN.  Oligomerization and ligand binding in a homotetrameric hemoglobin: two high-resolution crystal structures of hemoglobin Bart's (gamma(4)), a marker for alpha-thalassemia,  Protein Sci. ,  2001, vol.  10 (pg.  1739- 1749) Google Scholar CrossRef Search ADS PubMed  61 Liddington R,  Derewenda Z,  Dodson E,  Hubbard R,  Dodson G.  High resolution crystal structures and comparisons of T-state deoxyhaemoglobin and two liganded T-state haemoglobins: T(α-oxy)haemoglobin and T(met)haemoglobin,  J. Mol. Biol. ,  1992, vol.  228 (pg.  551- 579) Google Scholar CrossRef Search ADS PubMed  62 Yasuda J,  Ichikawa T,  Tsuruga M,  Matsuoka A,  Sugawara Y,  Shikama K.  The alpha 1 beta 1 contact of human hemoglobin plays a key role in stabilizing the bound dioxygen,  Eur. J. Biochem. ,  2002, vol.  269 (pg.  202- 211) Google Scholar CrossRef Search ADS PubMed  63 Yamamoto Y,  Nagaoka T.  A 1H NMR comparative study of human adult and fetal hemoglobins,  FEBS Lett. ,  1998, vol.  424 (pg.  169- 172) Google Scholar CrossRef Search ADS PubMed  64 Yamamoto Y,  Nanai N,  Chujo R,  Suzuki T.  Heme methyl hyperfine shift pattern as a probe for determining the orientation of the functionally relevant proximal histidyl imidazole with respect to the heme in hemoproteins,  FEBS Lett. ,  1990, vol.  264 (pg.  113- 116) Google Scholar CrossRef Search ADS PubMed  65 Shokhirev NV,  Walker FA.  The effect of axial ligand plane orientation on the contact and pseudocontact shifts of low-spin ferriheme proteins,  J. Biol. Inorg. Chem. ,  1998, vol.  3 (pg.  581- 594) Google Scholar CrossRef Search ADS   66 Vojtěchovský J,  Chu K,  Berendzen J,  Sweet RM,  Schlichting I.  Crystal structures of myoglobin-ligand complexes at near-atomic resolution,  Biophys. J. ,  1999, vol.  77 (pg.  2153- 2174) Google Scholar CrossRef Search ADS PubMed  67 Nagai M,  Nagai Y,  Aki Y,  Imai K,  Wada Y,  Nagatomo S,  Yamamoto Y.  Effect of reversed heme orientation on circular dichroism and cooperative oxygen binding of human adult hemoglobin,  Biochemistry ,  2008, vol.  47 (pg.  517- 525) Google Scholar CrossRef Search ADS PubMed  68 Sugita Y,  Yoneyama Y.  Oxygen equilibrium of hemoglobins containing unnatural hemes. Effect of modification of heme carboxyl groups and side chains at positions 2 and 4,  J. Biol. Chem. ,  1971, vol.  246 (pg.  389- 394) Google Scholar PubMed  © The Authors 2010. Published by Oxford University Press on behalf of the Japanese Biochemical Society. All rights reserved TI - Characterization of the acid–alkaline transition in the individual subunits of human adult and foetal methaemoglobins JF - The Journal of Biochemistry DO - 10.1093/jb/mvq055 DA - 2010-06-02 UR - https://www.deepdyve.com/lp/oxford-university-press/characterization-of-the-acid-alkaline-transition-in-the-individual-3Oo70TCjp7 SP - 217 EP - 229 VL - 148 IS - 2 DP - DeepDyve ER -