Titration behavior and tautomeric states of individual histidine residues of myoglobins. Application of natural abundance carbon 13 nuclear magnetic resonance spectroscopy

The Biophysical Probes 2-fluorohistidine and 4-fluorohistidine: Spectroscopic Signatures and Molecular Properties

Fluorinated amino acids serve as valuable biological probes, by reporting on local protein structure and dynamics through ¹⁹F NMR chemical shifts. 2-fluorohistidine and 4-fluorohistidine, studied here with DFT methods, have even more capabilities for biophysical studies, as their altered pK_a values, relative to histidine, allow for studies of the role of proton transfer and tautomeric state in enzymatic mechanisms. Considering the two tautomeric forms of histidine, it was found that 2-fluorohistidine primarily forms the common (for histidine) τ-tautomer at neutral pH, while 4-fluorohistidine exclusively forms the less common π-tautomer. This suggests the two isomers of fluorohistidine can also serve as probes of tautomeric form within biomolecules, both by monitoring NMR chemical shifts and by potential perturbation of the tautomeric equilibrium within biomolecules. Fluorine also enables assignment of tautomeric states in crystal structures. The differences in experimental pK_a values between the isomers was found to arise from solvation effects, providing insight into the polarization and molecular properties of each isomer. Results also encompass ¹³C and ¹⁹F NMR chemical shifts, from both tautomers of 2-fluorohistidine and 4-fluorohistidine in a number of different environments. This work can serve as a guide for interpretation of spectroscopic results in biophysical studies employing 2-fluorohistidine and 4-fluorohistidine.

A Solid State 13C NMR, Crystallographic, and Quantum Chemical Investigation of Chemical Shifts and Hydrogen Bonding in Histidine Dipeptides

We report the first solid-state NMR, crystallographic, and quantum chemical investigation of the origins of the 13C NMR chemical shifts of the imidazole group in histidine-containing dipeptides. The chemical shift ranges for Cgamma and Cdelta2 seen in eight crystalline dipeptides were very large (12.7-13.8 ppm); the shifts were highly correlated (R2= 0.90) and were dominated by ring tautomer effects and intermolecular interactions. A similar correlation was found in proteins, but only for buried residues. The imidazole 13C NMR chemical shifts were predicted with an overall rms error of 1.6-1.9 ppm over a 26 ppm range, by using quantum chemical methods. Incorporation of hydrogen bond partner molecules was found to be essential in order to reproduce the chemical shifts seen experimentally. Using AIM (atoms in molecules) theory we found that essentially all interactions were of a closed shell nature and the hydrogen bond critical point properties were highly correlated with the N...H...O (average R2= 0.93) and Nepsilon2...H...N (average R2= 0.98) hydrogen bond lengths. For Cepsilon1, the 13C chemical shifts were also highly correlated with each of these properties (at the Nepsilon2 site), indicating the dominance of intermolecular interactions for Cepsilon1. These results open up the way to analyzing 13C NMR chemical shifts, tautomer states (from Cdelta2, Cepsilon1 shifts), and hydrogen bond properties (from Cepsilon1 shifts) of histidine residue in proteins and should be applicable to imidazole-containing drug molecules bound to proteins, as well.

Structural investigations of the active-site mutant Asn156Ala of outer membrane phospholipase A: function of the Asn-His interaction in the catalytic triad

Outer membrane phospholipase A (OMPLA) from Escherichia coli is an integral-membrane enzyme with a unique His-Ser-Asn catalytic triad. In serine proteases and serine esterases usually an Asp occurs in the catalytic triad; its role has been the subject of much debate. Here the role of the uncharged asparagine in the active site of OMPLA is investigated by structural characterization of the Asn156Ala mutant. Asparagine 156 is not involved in maintaining the overall active-site configuration and does not contribute significantly to the thermal stability of OMPLA. The active-site histidine retains an active conformation in the mutant notwithstanding the loss of the hydrogen bond to the asparagine side chain. Instead, stabilization of the correct tautomeric form of the histidine can account for the observed decrease in activity of the Asn156Ala mutant.

Crystal structures of myoglobin-ligand complexes at near-atomic resolution

We have used x-ray crystallography to determine the structures of sperm whale myoglobin (Mb) in four different ligation states (unligated, ferric aquomet, oxygenated, and carbonmonoxygenated) to a resolution of better than 1.2 A. Data collection and analysis were performed in as much the same way as possible to reduce model bias in differences between structures. The structural differences among the ligation states are much smaller than previously estimated, with differences of <0.25 A root-mean-square deviation among all atoms. One structural parameter previously thought to vary among the ligation states, the proximal histidine (His-93) azimuthal angle, is nearly identical in all the ferrous complexes, although the tilt of the proximal histidine is different in the unligated form. There are significant differences, however, in the heme geometry, in the position of the heme in the pocket, and in the distal histidine (His-64) conformations. In the CO complex the majority conformation of ligand is at an angle of 18 +/- 3 degrees with respect to the heme plane, with a geometry similar to that seen in encumbered model compounds; this angle is significantly smaller than reported previously by crystallographic studies on monoclinic Mb crystals, but still significantly larger than observed by photoselection. The distal histidine in unligated Mb and in the dioxygenated complex is best described as having two conformations. Two similar conformations are observed in MbCO, in addition to another conformation that has been seen previously in low-pH structures where His-64 is doubly protonated. We suggest that these conformations of the distal histidine correspond to the different conformational substates of MbCO and MbO(2) seen in vibrational spectra. Full-matrix refinement provides uncertainty estimates of important structural parameters. Anisotropic refinement yields information about correlated disorder of atoms; we find that the proximal (F) helix and heme move approximately as rigid bodies, but that the distal (E) helix does not.

Theory of electrostatic interactions in macromolecules

In the past year, substantial progress has been made in the modeling of electrostatic interactions in biomolecules. This review highlights advances in the following areas: first, the efficient computation of long-range electrostatic interactions in detailed molecular simulations; second, the application of the Poisson-Boltzmann electrostatic model in conformational analysis; third, the application of the Poisson-Boltzmann model in quantum chemistry calculations; fourth, the development of atomic parameters; and finally, the modeling of ionization equilibria in proteins.

Are the electrospray mass spectra of proteins related to their aqueous solution chemistry?

Institute for Environmental Chemistry, National Research Council of Canada, Ottawa, Ontario, Canada The shape of the profile described by the relative abundances of multiply charged ions of proteins in the electrospray mass spectrum can be described by means of one or more Gaussian functions. An aqueous solution equilibrium model of the distribution of multiply charged ions of equine cytochrome c and myoglobin has been shown to match qualitatively the shape of the distribution of these ions in an electrospray mass spectrum. Monotonic functions such as the quadrupole mass spectrometric transmission efficiency may alter the centroid of the profile, but the shape of the ion abundance pattern appears to be controlled by the aqueous solution chemistry of the proteins. NRCC No. 32938.

Distal and proximal ligand interactions in heme proteins: correlations between C-O and Fe-C vibrational frequencies, oxygen-17 and carbon-13 nuclear magnetic resonance chemical shifts, and oxygen-17 nuclear quadrupole coupling constants in C17O- and 13CO-labeled species

We have obtained the oxygen-17 nuclear magnetic resonance (NMR) spectra of a variety of C17O-labeled heme proteins, including sperm whale (Physeter catodon) myoglobin, two synthetic sperm whale myoglobin mutants (His E7----Val E7; His E7----Phe E7), adult human hemoglobin, rabbit (Oryctolagus cuniculus) hemoglobin, horseradish (Cochlearia armoracia) peroxidase (E.C. 1.11.1.7) isoenzymes A and C, and Caldariomyces fumago chloroperoxidase (E.C. 1.11.1.10), in some cases as a function of pH, and have determined their isotropic 17O NMR chemical shifts, delta i, and spin-lattice relaxation times, T1. We have also obtained similar results on a picket fence prophyrin, [5,10,15,20-tetrakis(alpha, alpha, alpha, alpha, alpha-pivalamidophenyl)porphyrinato]iron(II) (1-MeIm)CO, both in solution and in the solid state. Our results show an excellent correlation between the infrared C-O vibrational frequencies, v(C-O), and delta i, between v(C-O) and the 17O nuclear quadrupole coupling constant (e2qQ/h, derived from T1), and as expected between e2qQ/h and delta i. Taken together with the work of others on the 13C NMR of 13CO-labeled proteins, where we find an excellent correlation between delta i(13C) and v(Fe-C), our results suggest that IR and NMR measurements reflect the same interaction, which is thought to be primarily the degree of pi-back-bonding from Fe d to CO pi* orbitals, as outlined previously [Li, X.-Y., & Spiro, T.G. (1988) J. Am. Chem. Soc. 110, 6024]. The modulation of this interaction by the local charge field of the distal heme residue (histidine, glutamine, arginine, and possibly lysine) in a variety of species and mutants, as reflected in the NMR and IR measurements, is discussed, as is the effect of cysteine as the proximal heme ligand.

Resonance Raman evidence that distal histidine protonation removes the steric hindrance to upright binding of carbon monoxide by myoglobin

The resonance Raman band assigned to Fe--CO stretching in the sperm whale myoglobin CO adduct shifts from 507 cm-1 at neutral pH to 488 cm-1 at low pH, in concert with a shift of the C-O stretching infrared band from 1947 to 1967 cm-1 (Fuchsman & Appleby, 1979), while the 575-cm-1 Fe-C-O bending RR band loses intensity. The pKa that characterizes these changes is approximately 4.4. The vibrational frequencies at low pH are well modeled by the protein-free CO, imidazole adduct of protoheme in a nonpolar solvent while those at high pH are modeled by the adduct of a heme with a covalent strap (Yu et al., 1983) which inhibits upright CO binding. It is inferred that the Fe-C-O unit changes from a tilted to an upright geometry when the distal histidine is protonated, because its side chain swings out of the heme pocket due to electrostatic repulsion with a nearby arginine residue. A different protonation step (pKa = 5.7), which has been shown to modulate the CO rebinding kinetics (Doster et al., 1982) as well as the optical spectrum (Fuchsman & Appleby, 1979), is suggested to involve a global structure change associated with protonation of histidine residues distant from the heme.

Conformational differences between various myoglobin ligated states as monitored by 1H NMR spectroscopy

In paramagnetic metmyoglobin, cyanomyoglobin (CNMb), and deoxymyoglobin, His-36 has a high pK (approximately 8), and the NMR titration behavior of the H-2 resonance is perturbed, due to the presence at low pH of a hydrogen bond with Glu-38, which is broken at high pH. The His-36 H-4 resonance shows no shift with pK approximately 8 because of two opposing chemical shift effects but monitors the titration of nearby Glu-36 (pK = 5.6). In diamagnetic derivatives [(carbon monoxy)myoglobin (COMb) and oxymyoglobin (oxyMb)], the titration behavior of His-36 H-2 and H-4 resonances is normalized (pK approximately 6.8). The very slight alkaline Bohr effect in sperm whale myoglobin (Mb) is interpreted in terms of the pK change of His-36 from deoxyMb to oxyMb and compensating pK changes in the opposite direction of other unspecified groups. In sperm whale COMb at 40 degrees C, the distal histidine (His-64) and His-97 have pK values of 5.0 and 5.9. The meso proton resonances remote from these groups do not show a titration shift, but the nearby gamma-meso proton (pK = 5.3) responds to titration of both histidines, and the upfield Val-68 methyl at -2.3 ppm (pK = 4.7) witnesses the titration of nearby His-64. At 20 degrees C, the latter resonance is reduced in size, and a second resonance occurs at -2.8 ppm, which is insensitive to pH and, hence, more remote from His-64. Both resonances arise from two conformations of Val-68 in slow equilibrium. In oxyMb at 20 degrees C, only the latter resonance is observed, presumably because of the steric restrictions imposed by the hydrogen bond between ligand and His-64 in oxyMb, which is not present in COMb. In oxyMb the pK of His-97 (5.6) is similar to that of the meso proton resonances (5.5) and to the pK of other pH-dependent processes, including the very small acid Bohr effect. It is likely that these processes are controlled by the titration of His-97.(ABSTRACT TRUNCATED AT 250 WORDS)

Assignment of 1H NMR resonances of histidine and other aromatic residues in met-, cyano-, oxy-, and (carbon monoxy)myoglobins

The resolved 1H NMR resonances of the aromatic region in the 270-MHz NMR spectrum of sperm whale, horse, and pig metmyoglobin (metMb) have been assigned, including the observable H-2 and H-4 histidine resonances, the tryptophan H-2 resonances, and upfield-shifted resonances from one tyrosine residue. The use of different Mb species, carboxymethylation, and matching of pK values allows the assignment of the H-4 resonances, which agree in only three cases out of seven with scalar-correlated two-dimensional NMR spectroscopy assignments by others. The conversion to hydroxymyoglobin at high pH involves rearrangements throughout the molecule and is observed by many assigned residues. In sperm whale ferric cyanomyoglobin, nine H-2 and eight H-4 histidine resonances have been assigned, including the His-97 H-2 resonance and tyrosine resonances from residues 103 and 146. The hyperfine-shifted resonances from heme and near-heme protons observe a shift with a pK = 5.3 +/- 0.3 (probably due to deprotonation of His-97, pK = 5.6) and another shift at pK = 10.8 +/- 0.3. The spectrum of high-spin ferrous sperm whale deoxymyoglobin is very similar to that of metMb, which allows the assignment of seven surface histidine H-2 and H-4 resonances and also resonances from the two tryptophan residues and one tyrosine. In diamagnetic sperm whale (carbon monoxy)myoglobin (COMb), 10 His H-2 and 11 His H-4 resonances are observed, and 8 H-2 and 9 H-4 resonances are assigned, including His-64 H-4, the distal histidine. This important resonance is not observed in sperm whale oxymyoglobin, which in general shows very similar titration curves to COMb. Histidine-36 shows unusual titration behavior in the paramagnetic derivatives but normal behavior in the diamagnetic derivatives, which is discussed in the accompanying paper [Bradbury, J. H., & Carver, J. A. (1984) Biochemistry (following paper in this issue)].

Resonance Raman and absorption spectroscopic detection of distal histidine--fluoride interactions in human methemoglobin fluoride and sperm whale metmyoglobin fluoride: measurements of distal histidine ionization constants

The pH dependence of the resonance Raman and absorption spectra of human methemoglobin fluoride (HbIIIF) and sperm whale metmyoglobin fluoride (MbIIIF) has been examined. Both the Raman and absorption spectra of HbIIIF and MbIIIF indicate the existence at alkaline pH of an equilibrium between the hydroxide and fluoride complexes. The absorption maxima of HbIIIF and MbIIIF solutions shift to longer wavelengths as the pH is decreased from neutrality. The Raman data show a corresponding shift of the 461- and 468-cm-1 Fe-F vibrational stretching peaks at pH 7.0 [Asher, S. A., & Schuster, T. M. (1979) Biochemistry 18, 5377] to 399 and 407 cm-1 at acid pH in MbIIIF and HbIIIF, respectively. These shifts are interpreted to result from protonation of the distal histidine and the formation of a hydrogen bond to the fluoride ligand. Measurements of the pH dependence of the absorption and resonance Raman spectra give distal histidine ionization constants (apparent) corresponding to pK = 5.1 (+/- 0.1) for HbIIIF and pK = 5.5 (+/- 0.1) for MbIIIF. An examination of the distal histidine pK values and the frequency of the hydrogen-bonded Fe--F stretching vibration at pH 5.0 of HbIIIF with and without inositol hexaphosphate indicates little difference in the distal histidine--heme distance between the so-called R and T quaternary forms of HbIIIF. These results indicate that the changes in the electronic spectrum of HbIIIF that occur upon switching from the R to the T form do not result from alterations in (1) the iron--fluoride bond distance, (2) the iron out-of-heme plane distance, or (3) the distal histidine--fluoride distance.

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.

CO and O2 complexes of soybean leghemoglobins: pH effects upon infrared and visible spectra. Comparisons with CO and O2 complexes of myoglobin and hemoglobin

The effects of pH upon infrared spectra [CO stretching frequency (vco) region] and visible spectra of the CO complexes of soybean leghemoglobins a, c1, and c2, sperm whale myoglobin, and human hemoglobin A are reported. The vco for leghemoglobin--CO complexes was 1947.5 cm-1 at neutral pH. At acid pH myoglobin-- and hemoglobin--CO complexes developed vco bands at 1966--1968 cm-1, whereas leghemoglobin--CO complexes developed vco bands at approximately 1957 cm-1. All pKapp co values determined by pH-dependent variation of vco fell in the range 4.0--4.6. The pKapp co values determined from visible spectra were consistent with vco-determined values except for that of myoglobin--CO (visible pKapp co = 5.8). The pKapp co values in the 4.0--4.6 range appear to be pK values of the distal histidines, while the visible pKapp co of myoglobin--CO appears to be the pK of a group other than the distal and proximal histidines. The data are consistent with a model in which protonation of the distal histidine permits protein-free heme FeCO geometry in leghemoglobin--CO complexes but not in myoglobin-- or hemoglobin--CO complexes. Thus the heme pockets of leghemoglobins appear to be more flexible than the heme pockets of myoglobin and hemoglobin. The effects of pH upon visible spectra of the O2 complexes of soybean leghemoglobins a, c1, and c2, sperm whale myoglobin, and human hemoglobin A also are reported. pKapp o2 values of approximately 5.5 (leghemoglobins) and 4.4 (hemoglobin) are probably the pK values of the distal histidines. Comparisons of pKapp o2 values with pKapp co values indicate a more flexible heme pocket in leghemoglobins than in hemoglobin. The O2 complex of leghemoglobin c2 differed significantly from the O2 complexes of leghemoglobins a and c1 in visible spectra and titration behavior. These differences might be associated with the small structural differences in the region between the E and F helixes of leghemoglobins.