Electrostatic effects in hemoglobin: Bohr effect and ionic strength dependence of individual groups.

The electrostatic treatment applied in the preceding paper in this issue [Matthew, J. B., Hanania, G.I.H., & Gurd, F.R.N. (1979) Biochemistry (preceding paper in this issue)] to the titration behavior of individual groups in human deoxyhemoglobin and oxyhemoglobin was applied to the computation of the alkaline Bohr effect at various values of ionic strength. The enhanced proton binding of deoxyhemoglobin in the pH range of 6--9 was accounted for at ionic strength 0.01 M by the effects of the unique charge distributions of ionizable groups in the two quaternary states. At ionic strength 0.10 M the effects of 2--4 bound anions had to be considered in addition in the deoxyhemoglobin charge configuration. At the higher ionic strength 10 groups per tetramer contributed to the Bohr effect, whereas 28 groups were contributory at the lower ionic strength. The ionic strength dependence of individual groups in the two tetrameric structures as well as in the alpha-chain monomer was explained in terms of the electrostatic treatment. This examination showed that the differences in electrostatic behavior of deoxy- and oxyhemoglobin follow from particular dissymmetries in their configurations with respect to charge and static solvent accessibility.     

Discrete charge calculations of potentiometric titrations for globular proteins: sperm whale myoglobin, hemoglobin alpha chain, cytochrome c.

The modified Tanford-Kirkwood theory of Shire et al. for intramolecular electrostatic interactions has been applied to hydrogen ion equilibria of sperm whale ferrimyoglobin, human hemoglobin α-chain and horse cytochrome c. The model employs two sets of parameters derived from the crystalline protein structures, first, the atomic coordinates of charged amino acid residues and, second, static accessibility factors to reflect their solvent exposure. In addition, a consistent set of intrinsic pK values (pK_int) for the individual groups is employed. The theoretical pK values at half-titration for individual groups in each protein correspond to the available observed pK values, and the theoretical titration curves compare closely with experimental potentiometric curves.     

Electrostatic effects in myoglobin. Application of the modified Tanford-Kirkwood theory to myoglobins from horse, California grey whale, harbor seal, and California sea lion.

The modified Tanford-Kirkwood electrostatic theory (Shire et al., 1974a) was applied to ferrimyoglobins from the following animal species: sperm whale (Physeter catodon), horse, California grey whale (Eschrichtius gibbosus), harbor seal (Phoca vitulina), and California sea lion (Zalophus californianus). Computations were made of the overall hydrogen ion titration curves of the proteins, and of pH and ionic strength variations of ionization equilibria for individual groups in the protein, with particular reference to the hemic acid ionization of the iron bound water molecule. Coordinates and static solvent accessibility were estimated in terms of the sperm whale myoglobin structure. Where possible, theoretical results and experimental data are compared. Some comparative features of charge and ionization properties among the various myoglobins are presented.     

Electrochemical reduction of methemoglobin either directly or with flavin mononucleotide as a mediator

Electrochemical reduction of methemoglobin on a platinum electrode is studied by means of thin layer spectroelectrochemistry. For methemoglobin alone in solution, direct reduction is very slow even for potentials close to those of the reduction of the solvent. The reduction of a methemoglobin-oxyhemoglobin mixture with an imposed potential causes the electrochemical reduction of oxygen, the conversion of oxyhemoglobin into deoxyhemoglobin, and a simultaneous transformation of part of the molecules into methemoglobin. When fixed oxygen has disappeared, reduction of methemoglobin takes place. The reduction of methemoglobin and deoxyhemoglobin is catalyzed by the presence of flavin mononucleotide (FMN). For the oxyhemoglobin-methemoglobin mixture, flavin makes a fast deoxygenation of oxyhemoglobin without a change in the oxidation state of the iron. It also allows the rapid reduction of methemoglobin. In each case, the resulting deoxyhemoglobin solutions do not show any electrolysis-induced modification of the equilibrium curves for oxygen binding.     

Calculation of the energy difference between the quaternary structures of deoxy- and oxyhemoglobin.

The energy difference between the quaternary structures of deoxy- and oxyhemoglobin is evaluated on the basis of the atomic coordinates determined by X-ray diffraction analysis. Calculation of the van der Waals interaction between subunits shows that in a hemoglobin molecule as a whole, the interaction is more attractive in the oxy form than in the deoxy form by about 8 kcal/mol, and that in each pair of two subunits except the pair alpha1alpha2, the interaction energy varies by about 15 kcal/mol. The electrostatic interactions originating in the partial charges on all constituent atoms of hemoglobin and in the polar residues on the surface of hemoglobin make only a small contribution to the energy difference between the quaternary structures of deoxy- and oxyhemoglobin. Thus, the contribution of the clusters of the polar residues in the internal cavity between like subunits and also of the freedom of rotation of the C-terminal of each subunit in oxyhemoglobin may be important energetically in the transition from deoxy to oxy quaternary structure. In this point, the present calculation supports Perutz' model, but suggests necessity of further investigations on the transitional characteristics of the quaternary structure in the intermediate steps of oxygenation. The discussion on the transitional characteristics is given in the last section.     

Electrostatic interactions in sperm whale myoglobin. Site specificity, roles in structural elements, and external electrostatic potential distributions

The electrostatic free energy contribution to the stability of sperm whale ferrimyoglobin was evaluated according to the static accessibility modified Tanford-Kirkwood model. The electrostatic free energy contribution of each distinct structural element was divided into one term arising from interactions between it and other elements (interelemental) and another from interactions within the particular element itself (intraelemental). At pH 7 the majority of the terms were found to be stabilizing. The interelemental terms are the dominant ones for most structural elements. The small interelemental terms of the C and D helices are compensated by large intraelemental interactions which stabilize these short helices. Perturbations in pH can be accommodated by the structural elements through a redistribution of stabilizing and destabilizing interactions. The electrostatic potentials calculated at the surface of the protein indicate that the internal compensation of local potentials achieved during folding results in a generally neutral protein-solvent interface save for two distinct areas of nonzero potential. The accessibility of each charged atom to solvent was analyzed in terms of the surface area lost to charged, polar and nonpolar atoms separately. The net solvent accessibility lost parallels closely that lost to nonpolar atoms alone, indicating a specific role for nonpolar atoms in defining dielectric shielding of charged atoms, aside from their participation in the well-known hydrophobic interactions.     

The mutation beta 99 Asp-Tyr stabilizes Y--a new, composite quaternary state of human hemoglobin

Carbonmonoxy hemoglobin Ypsilanti (beta 99 Asp-Tyr) exhibits a quaternary form distinctly different from any structures previously observed for human hemoglobins. The relative orientation of alpha beta dimers in the new quaternary form lies well outside the range of values observed for normal unliganded and liganded tetramers (Baldwin, J., Chothia, C., J. Mol. Biol. 129:175-220, 1979). Despite this large quaternary structural difference between carbonmonoxy hemoglobin Ypsilanti and the two canonical structures, the new quaternary structure's hydrogen bonding interactions in the "switch" region, and packing interactions in the "flexible joint" region, show noncovalent interactions characteristic of the alpha 1 beta 2 contacts of both unliganded and liganded normal hemoglobins. In contrast to both canonical structures, the beta 97 histidine residue in carbonmonoxy hemoglobin Ypsilanti is disengaged from quaternary packing interactions that are generally believed to enforce two-state behavior in ligand binding. These features of the new quaternary structure, denoted Y, may therefore be representative of quaternary states that occur transiently along pathways between the normal unliganded, T, and liganded, R, hemoglobin structures.     

Observation of allosteric transition in hemoglobin

Two conclusions have been drawn from NMR studies of mixed state hemoglobins. First the α and β subunits in hemoglobin are not equivalent in their conformational properties. Second the mixed state hemoglobin (α^IIICN β^II)₂ can take two different quaternary structures without changing the degree of ligation. One of the two structures is similar to that of deoxyhemoglobin and the other to that of oxyhemoglobin.     

Electrostatic stabilization in myoglobin. Interactive free energies between individual sites.

The pattern of electrostatic interactions between pairs of charge sites in sperm whale ferrimyoglobin was examined as a function of pH in terms of proton site occupancy, static solvent accessibility, and distance of separation. By grouping all examples of the most stabilizing interactions and all examples of the most destabilizing interactions, we can easily show that at pH 7.50 the former is much stronger; that is, the negative contributions to electrostatic free energy far outweigh the positive contributions. Much of the electrostatic energy of stabilization in native myoglobin is provided by specific charge-pair partners that are very highly conserved among 53 mammalian myoglobin species and is invariant substantially from pH 8.5 to 3.5. Destablizing interactions that become most significant, but not actually dominant, near the acid unfolding pH range can be recognized in emerging clusters of uncompensated positive charges. Binding of azide ion by the heme iron effectively reduces the most prominent destabilizing set of such interactions. In general, thoe charged residues that experience the largest summed stabilizing interactions with other groups are the most conserved between species. The histidine residues, however, show their best correlation of conservation with low values of static accessibility. Although histidine residue 64 has an effective pK corresponding to the midpoint of the unfolding transition near pH 4.2 at an ionic strength of 0.10 M and so might be called a "trigger group", its interactions contribute only a modest fraction of the overall pH-dependent free energy change. An examination of the primary stabilizing interactions represented by the charge-pair partners indicates a probably major role of electrostatic interactions in the nucleation and docking stages of the condensation of the polypeptide chain into the compact native structure.     

The assignment of carbon monoxide association rate constants to the alpha and beta subunits in native and mutant human deoxyhemoglobin tetramers

The association kinetics of CO binding to site-directed mutants of human deoxyhemoglobin were measured by stopped-flow rapid mixing techniques at pH 7.0, 20 degrees C. Hemoglobin tetramers were constructed from one set of native subunits and one set of mutated partners containing His(E7) to Gly, Val(E11) to Ala, or Val(E11) to Ile substitutions. The reactivity of beta Cys93 with p-hydroxymercuribenzoate was measured to ensure that the mutant deoxyhemoglobins were capable of forming T-state quaternary conformations. Time courses for the complete binding of CO were measured by mixing the deoxygenated proteins with a 5-fold excess of ligand in the absence and presence of inositol hexaphosphate. Association rate constants for the individual alpha and beta subunits in the T-state conformation were assigned by measuring time courses for the reaction of a small, limiting amount of CO with a 20-fold excess of deoxyhemoglobin (i.e. Hb4 + CO----Hb4(CO)). The effects of the E7 and E11 mutations in T-state alpha subunits were qualitatively similar to those observed for the same subunit in the R-state (Mathews, A.J., Rohlfs, R.J., Olson, J.S., Tame, J., Renaud, J-P., and Nagai, K. (1989) J. Biol. Chem. 264, 16573-16583). The alpha His58(E7) to Gly and Val62(E11) to Ala substitutions caused 80- and 3-fold increases, respectively, in k'CO for T-state alpha subunits, and the alpha Val62(E11) to Ile mutation caused a 3-fold decrease. The beta His63(E7) to Gly and Val67(E11) to Ala substitutions produced 70- and 8-fold increases, respectively, in k'CO for T-state beta subunits whereas these same mutations caused little effect on the rate of CO binding to R-state beta subunits. The beta Val67(E11) to Ile mutation produced the same large effect, a 23-fold reduction in k'CO, in both quaternary conformations of beta subunits. These kinetic results can be interpreted qualitatively in terms of differences between the alpha and beta subunits in the deoxy and liganded crystal structures of human hemoglobin (Perutz, M.F. (1990) Annu. Rev. Physiol. 52, 1-25). Both the structural and functional data suggest that the distal portion of the beta heme pocket is tightly packed in deoxyhemoglobin whereas the CO binding site in R-state beta subunits is much more open. In contrast, the distal portion of the alpha heme pocket is restricted sterically in both quaternary states.(ABSTRACT TRUNCATED AT 400 WORDS)     

Mechanism of surface peptide proton exchange in bovine pancreatic trypsin inhibitor. Salt effects and O-protonation

The acid-catalyzed hydrogen exchange rate constants kH, and the base-catalyzed rate constants kOH, have been determined (in the preceding paper) for the 25 most rapidly exchanging NH groups of bovine pancreatic trypsin inhibitor. Most of these NH groups are at the protein-solvent interface. The correlation of kH, but not kOH, with the static accessibility and hydrogen bonding of the peptide carbonyl O atom indicates that the mechanism of acid catalysis in proteins involves O-protonation. Agreement between the ionic strength dependence observed for kH and kOH and the ionic strength dependence calculated for an O-protonation mechanism supports this conclusion. N-protonation for acid catalysis, as well as N-deprotonation for base catalysis, have traditionally been assumed in the mechanism of the chemical step in peptide amide proton exchange. A preference for the alternative O-protonation mechanism has far-reaching implications in the interpretation of protein hydrogen exchange kinetics. With an O-protonation mechanism, acid-catalyzed rates of surface NH groups are primarily a function of the average solvent accessibility of the carbonyl O atoms in the dynamic solution structure, while base-catalyzed rates of surface NH groups measure solvent accessibility of the peptide N. The relative dynamic accessibilities of peptide O atoms, as measured by relative values of kH (corrected for electrostatic effects), correlate with O static accessibilities in the crystal structure. A lower correlation of static accessibility of N atoms with kOH is observed for surface NH groups in peptide groups in which the carbonyl O is not hydrogen bonded. For some surface NH groups, the observed pH of minimum rate, pHmin, deviates widely from the pHmin of model compounds. This is explained as the combined result of electrostatic effects and of the differences in accessibility of the carbonyl O and N atoms that result in a change in the relative values of kH and kOH as compared to those of model peptides. A mechanism whereby exchange of interior sites is catalyzed by interactions of catalysis ions with protein surface atoms via charge transfer is suggested.     

Quaternary structure of carbonmonoxyhemoglobins in solution: structural changes induced by the allosteric effector inositol hexaphosphate

We have applied the residual dipolar coupling (RDC) method to investigate the solution quaternary structures of (2)H- and (15)N-labeled human normal adult recombinant hemoglobin (rHb A) and a low-oxygen-affinity mutant recombinant hemoglobin, rHb(alpha96Val-->Trp), both in the carbonmonoxy form, in the absence and presence of an allosteric effector, inositol hexaphosphate (IHP), using a stretched polyacrylamide gel as the alignment medium. Our recent RDC results [Lukin, J. A., Kontaxis, G., Simplaceanu, V., Yuan, Y., Bax, A., and Ho, C. (2003) Proc. Natl. Acad. Sci. U.S.A. 100, 517-520] indicate that the quaternary structure of HbCO A in solution is a dynamic ensemble between two previously determined crystal structures, R (crystals grown under high-salt conditions) and R2 (crystals grown under low-salt conditions). On the basis of a comparison of the geometric coordinates of the T, R, and R2 structures, it has been suggested that the oxygenation of Hb A follows the transition pathway from T to R and then to R2, with R being the intermediate structure [Srinivasan, R., and Rose, G. D. (1994) Proc. Natl. Acad. Sci. U.S.A. 91, 11113-11117]. The results presented here suggest that IHP can shift the solution quaternary structure of HbCO A slightly toward the R structure. The solution quaternary structure of rHbCO(alpha96Val-->Trp) in the absence of IHP is similar to that of HbCO A in the presence of IHP, consistent with rHbCO(alpha96Val-->Trp) having an affinity for oxygen lower than that of Hb A. Moreover, IHP has a much stronger effect in shifting the solution quaternary structure of rHbCO(alpha96Val-->Trp) toward the R structure and toward the T structure, consistent with IHP causing a more pronounced decrease in its oxygen affinity. The results presented in this work, as well as other results recently reported in the literature, clearly indicate that there are multiple quaternary structures for the ligated form of hemoglobin. These results also provide new insights regarding the roles of allosteric effectors in regulating the structure and function of hemoglobin. The classical two-state/two-structure allosteric mechanism for the cooperative oxygenation of hemoglobin cannot account for the structural and functional properties of this protein and needs to be revised.     

A proton nuclear magnetic resonance investigation of the anion Bohr effect of human normal adult hemoglobin

High-resolution proton nuclear magnetic resonance spectroscopy has been used to investigate the molecular mechanism of the Bohr effect of human normal adult hemoglobin in the presence of two allosteric effectors, i.e., chloride and inorganic phosphate ions. The individual hydrogen ion equilibria of 22-26 histidyl residues of hemoglobin have been measured in anion-free 0.1 M HEPES buffer and in the presence of 0.18 M chloride or 0.1 M inorganic phosphate ions in both deoxy and carbonmonoxy forms. The results indicate that the beta 2-histidyl residues are strong binding sites for chloride and inorganic phosphate ions in hemoglobin. The affinity of the beta 2-histidyl residues for these anions is larger in the deoxy than in the carbonmonoxy form. Nevertheless, the contribution of these histidyl residues to the anion Bohr effect is small due to their low pK value in deoxyhemoglobin in anion-free solvents. The interactions of chloride and inorganic phosphate ions with the hemoglobin molecule also result in lower pK values and/or changes in the shapes of the hydrogen ion binding curves for several other surface histidyl residues. These results suggest that long-range electrostatic interactions between individual ionizable sites in hemoglobin could play an important role in the molecular mechanism of the anion Bohr effect.     

A thermodynamic model for the cooperative functional properties of the tetraheme cytochrome c3 from Desulfovibrio gigas.

A thermodynamic model is presented to describe the redox behaviour of the tetraheme cytochrome c3 from Desulfovibrio gigas. This molecule displays different intrinsic redox potentials for the four hemes and during the redox titration process, interactions among different hemes occur, thus altering the values of redox potentials according to which of the hemes are oxidized [Santos, H., Moura, J.J.G., Moura, I., LeGall, J. & Xavier, A.V. (1984) Eur. J. Biochem. 141, 283-296]. This complex cooperative behaviour [Xavier, A.V. (1986) J. Inorg. Biochem. 28, 239-243] has been analyzed here using an I2H4-interaction network [Cornish-Bowden, A. & Koshland, D.E. Jr (1970) J. Biol. Chem. 245, 6241-6250] coupled to a proton-linked equilibrium between two tertiary structures. Such a formalism, which requires a reduced number of parameters, is able to fully account quantitatively for the pH dependence of the NMR redox-titration curves. The 'redox-Bohr' effect is discussed in terms of the available structure and thermodynamic data and a functional mechanism is proposed.     

Distance-dependent dielectric constants and their application to double-helical DNA

Detailed studies of structures of biological macromolecules, even in simplified models, involve many costly and time-consuming calculations. Any thorough methods require sampling of an extremely large conformation and momentum space. Calculations of electrostatic interactions, which depend on many physical factors, such as the details of solvent, solvent accessibility in macromolecules, and molecular polarizability, are always developed in a compromise between more rigorous, detailed models and the need for immediate application to complicated biological systems. In this paper, a middle ground is taken between the more exact theoretical models and the simplest constant values for the dielectric constant. The effects of solvent, counterions, and molecular polarizability are incorporated through a set of adjustable parameters that should be determined from experimental conditions. Several previous forms for the dielectric function are compared with the new ones. The present methods use Langevin functions to span the region of dielectric constant between bulk solvent and cavity values. Application of such dielectric models to double-helical DNA is important because base-stacking preferences were previously demonstrated [A. Sarai, J. Mazur, R. Nussinov, and R. L. Jernigan (1988) Biochemistry, vol. 27, pp. 8498-8502] to be sensitive to the electrostatic formulation. Here we find that poly(dG).poly(dC) can be A form for high screening and B form for low screening. By contrast, poly(dA).poly(dT) can only take the B form. Base stacking is more sensitive to the form of the dielectric function than are the sugar-phosphate backbone conformations. Also in B form, the backbone conformations are not so affected by the base types as in A form.     

Comparison of reduced and oxidized yeast iso-1-cytochrome c using proton paramagnetic shifts.

Dipolar paramagnetic shifts for protons of yeast iso-1-cytochrome c have been calculated by using an optimized g-tensor and the X-ray crystallographic coordinates of the reduced form of yeast iso-1-cytochrome c [Louie, G. V., & Brayer, G. D. (1990) J. Mol. Biol. 214, 527-555]. The calculated values are compared with the observed paramagnetic shift determined from over 450 nonequivalent protons that have been assigned in both oxidation states [Gao, Y., Boyd, J., Williams, R. J. P., & Pielak, G. J. (1990) Biochemistry 29, 6994-7003]. There is good agreement between the calculated and the experimental data with a few exceptions. This indicates that, overall, the solution structures must be very similar in both the reduced and oxidized states in solution as is the case in crystals. The differences between observed and calculated shift values for the molecule in solution are most readily explained by slight movement of the heme and certain changes in diamagnetic shift due to small rearrangements of a few residues and some considerable changes in a few hydrogen bonds. It is also known that small differences exist between the structures of the two oxidation states in crystals but the hydrogen-bond changes are not so easily observed there. Structural changes from nuclear magnetic resonance data are in reasonable agreement with those deduced from crystallography, but additional information is clearly available concerning changes in hydrogen bonding.     

The three-dimensional structure of Bacillus amyloliquefaciens subtilisin at 1.8 A and an analysis of the structural consequences of peroxide inactivation

The three-dimensional structure of the subtilisin from Bacillus amyloliquefaciens (BAS) has been refined to 1.8 A using the amino acid sequence deduced from the DNA coding sequence. The structure is essentially the same as the previously reported structures of subtilisin BPN' (Wright, C.S., Alden, R.A., and Kraut, J. (1969) Nature 221, 235-242) and Novo (Drenth, J., Hol, W. G. J., Jansonius, J. N., and Koekoek, R. (1972) Eur. J. Biochem. 26, 177-181) determined in different crystal forms, at 2.5 and 2.8 A resolution, respectively. The largest differences in the three crystallographic models are seen in regions where the amino acid sequence used in the fit to the electron density maps of BPN' and Novo differs from the gene sequence of BAS (Wells, J. A., Ferrari, E., Henner, D. J., Estell, D. A., and Chen, E. Y. (1983) Nucleic Acids Res. 11, 7911-7925). The refined BAS model shows new features of cation binding, hydrogen bonding, and internal solvent structure. The refined BAS model has served as a basis for the analysis of stereochemical factors involved in the peroxide inactivation of the enzyme. Methionine 222, which is adjacent to the catalytic Ser221, is quantitatively oxidized to the sulfoxide by hydrogen peroxide as had been previously shown for the related Bacillus licheniformis enzyme (Stauffer, C. E., and Etson, D. (1969) J. Biol. Chem. 244, 5333-5338). In addition to this site of modification, we observe partial to full oxidation of two of the four remaining methionines. The oxidation of the methionines does not correlate well with their solvent accessibility calculated from the x-ray structure coordinates; in addition, only one of the two possible stereoisomers of methionine sulfoxide is formed. We also detect hydrogen peroxide-induced modification of the hydroxyl groups of two tyrosines. Modeling suggests that most of the observed effect of oxidation on the enzyme's catalytic efficiency can be attributed to unfavorable interactions at the oxyanion binding site between the sulfoxide group at 222 and the carbonyl oxygen of the scissile peptide bond of the bound substrate.     

Oxygenation of hemoglobin: Correspondence of crystal and solution properties using diffusion coefficient measurements

Quasi-elastic light scattering has been used to measure the change in the translational diffusion coefficient of hemoglobin upon oxygenation and the difference in the diffusion coefficients of oxy- and methemoglobin. The diffusion coefficients of oxy- and methemoglobin were found to be the same within the experimental accuracy of 0.2%, while the diffusion coefficient of oxyhemoglobin tetramers in solution at 13 mg/ml was found to be 0.8% smaller than that of deoxyhemoglobin at the same concentration, when the reversible dissociation of oxyhemoglobin tetramers into dimers was taken into account. In the limit of zero concentration, the oxyhemoglobin diffusion coefficient was found to be 1.5% ± 1.0% smaller than that of deoxyhemoglobin. This result is in very good agreement with what we predict using atomic coordinates to model the liganded and unliganded hemoglobin molecules as ellipsoids of revolution.     

The effect of potassium chloride on the Bohr effect of human hemoglobin.

The normal and differential titration curves of liganded and unliganded hemoglobin were measured at various KCl concentrations (0.1 to 2.0 M). In this range of KCl concentrations, the curves for deoxyhemoglobin showed no salt-induced pK changes of titratable groups. In the same salt concentration range oxyhemoglobin showed a marked change in titration behavior which could only be accounted for by a salt-induced increase in pK of some titratable groups. These results show that the suppression of the alkaline Bohr effect by high concentrations of neutral univalent salt is not caused by a weakening of the salt bridges in deoxyhemoglobin but is due to an interaction of chloride ions with oxyhemoglobin. Measurements of the Bohr effect at various KCl concentrations showed that at low chloride ion concentration (5 times 10-3 M) the alkaline Bohr effect is smaller than at a concentration of 0.1 M. This observation indicates that at a chloride ion concentration of 0.1 M, part of the alkaline Bohr effect is due to an interaction of chloride ions with hemoglobin. Furthermore, at low concentrations of chloride ions the acid Bohr effect has almost vanished. This result suggests that part of the acid Bohr effect arises from an interaction of chloride ions with oxyhemoglobin. The dependence of the Bohr effect upon the chloride ion concentration can be explained by assuming specific binding of chloride ions to both oxy- and deoxyhemoglobin, with deoxyhemoglobin having the highest affinity.     

Relationship between electrostatics and redox function in human thioredoxin: characterization of pH titration shifts using two-dimensional homo- and heteronuclear NMR.

The electrostatic behavior of potentially titrating groups in reduced human thioredoxin was investigated using two-dimensional (2D) 1H and 15N nuclear magnetic resonance (NMR) spectroscopy. A total of 241 chemical shift titration curves were measured over the pH range of 2.1-10.6 from homonuclear 1H-1H Hartmann-Hahn (HOHAHA) and heteronuclear 1H-15N Overbodenhausen correlation spectra. Nonlinear least-squares fits of the data to simple relationships derived from the Henderson-Hasselbalch equation led to the determination of pKas for certain isolated ionizable groups, including the single histidine residue at position 43 (pKa = 5.5 +/- 0.1) and a number of aspartic and glutamic acid carboxylate groups. Many of the titration curves demonstrate complex behavior due to the effects of interacting titrating groups, the long range of electrostatic interactions through the protein interior, and, perhaps, pH-induced conformational changes on the chemical shifts. Unambiguous assignment of the pKas for most of the 38 potentially ionizing groups of human thioredoxin could therefore not be made. In addition, there was no clear evidence that Asp-26 titrates in a manner corresponding to that observed in the Escherichia coli protein [Dyson, H. J., Tennant, L. L., & Holmgren, A. (1991) Biochemistry 30, 4262-4268]. The pKas of the active site cysteines were measured, however, with Cys-32 having an anomalously low value of 6.3 +/- 0.1 and that of Cys-35 between 7.5 and 8.6. These pKas are in agreement with proposed mechanisms for redox catalysis of thioredoxin and previously measured pKas within the active site of E. coli thioredoxin [Kallis, G. B., & Holmgren, A. (1980) J. Biol. Chem. 255, 10261-10265]. The stabilization of a thiolate anion at physiological pH can be explained by the interaction of the S gamma of Cys-32 with the amide of Cys-35 observed in the previously determined high-resolution solution structure of reduced human thioredoxin [Forman-Kay, J. D., Clore, G. M., Wingfield, P. T., & Gronenborn, A. M. (1991) Biochemistry 30, 2685-2698].