Modulation of the Axial Water Hydrogen-Bonding Properties by

The present results augur well for detecting and characterizing changes in distal water H-bonding upon mutagenesis of residues in the distal network o...
0 downloads 0 Views 206KB Size
Published on Web 04/21/2006

Modulation of the Axial Water Hydrogen-Bonding Properties by Chemical Modification of the Substrate in Resting State, Substrate-Bound Heme Oxygenase from Neisseria meningitidis; Coupling to the Distal H-Bond Network via Ordered Water Molecules Li-Hua Ma,† Yangzhong Liu,† Xuhong Zhang,‡ Tadashi Yoshida,‡ Kevin C. Langry,† Kevin M. Smith,§ and Gerd N. La Mar*,† Contribution from the Department of Chemistry, UniVersity of California, DaVis, California 95616, Department of Biochemistry, Yamagata UniVersity School of Medicine, Yamagata 990-9585, Japan, and Department of Chemistry, Louisiana State UniVersity, Baton Rouge, Louisiana 70803 Received November 18, 2005; E-mail: [email protected]

Abstract: The hydrogen bonding of ligated water in ferric, high-spin, resting-state substrate complexes of heme oxygenase from Neisseria meningitidis has been systematically perturbed by variable electronwithdrawing substituents on the hemin periphery. The pattern of 1H NMR-detected dipolar shifts due to the paramagnetic anisotropy is strongly conserved among the four complexes, with the magnitude of dipolar shifts or anisotropy increasing in the order of substituent formyl < vinyl < methyl. The magnetic anisotropy is axial and oriented by the axial Fe-His23 bond, and while individual anisotropies have uncertainties of ∼5%, the relative values of ∆χ (and the zero-field splitting constant, D ∝ ∆χax) are defined to 1%. The unique changes in the axial field strength implied by the variable zero-field splitting are in accord with expectations for the axial water serving as a stronger H-bond donor in the order of hemin substituents formyl > vinyl > methyl. These results establish the axial anisotropy (and D) as a sensitive probe of the H-bonding properties of a ligated water in resting-state, substrate complexes of heme oxygenase. Correction of observed labile proton chemical shifts for paramagnetic influences indicates that Gln49 and His53, some ∼10 Å from the iron, sense the change in the ligated water H-bonding to the three nonligated ordered water molecules that link the two side chains to the iron ligand. The present results augur well for detecting and characterizing changes in distal water H-bonding upon mutagenesis of residues in the distal network of ordered water molecules and strong H-bonds.

Introduction

Heme oxygenase,1 HO, is a nonmetal enzyme2-6 that uses heme as substrate and cofactor to cleave it at one meso position into biliverdin, iron, and CO. The enzyme is present in mammals, plants, and some pathogenic bacteria. The mammalian role includes iron homeostasis and the generation of biliverdin as autooxidant and CO as neural messenger,5,7 while plants use HO to produce open tetrapyrroles as precursors to light harvesting antennae.8 In bacteria, the primary role appears † ‡ §

University of California. Yamagata University School of Medicine. Louisiana State University.

(1) Tenhunen, R.; Marver, H. S.; Schmid, R. J. Biol. Chem. 1969, 244, 63886394. (2) Yoshida, T.; Migita, C. T. J. Inorg. Biochem. 2000, 82, 33-41. (3) Wilks, A. Antioxidants Redox Signal 2002, 4, 603-614. (4) Frankenberg-Dinkel, N. Antioxidants Redox Signal 2004, 6, 825-834. (5) Ortiz de Montellano, P. R.; Auclair, K. In The Porphryin Handbook; Kadish, K. M., Smith, K. M., Guilard, R., Eds.; Elsevier Science: San Diego, CA, 2003; pp 175-202. (6) Rivera, M.; Zeng, Y. J. Inorg. Biochem. 2005, 99, 337-354. (7) Maines, M. D. Annual ReView of Pharmacological Toxicology 1997, 37, 517-554. (8) Beale, S. I. Ciba Foundation Symposium 1994, 180, 156-168. 10.1021/ja0578505 CCC: $33.50 © 2006 American Chemical Society

to be for “mining” the host’s iron necessary for infection.3,4 The mechanism of the HO reaction appears to be the same in the diverse HOs, where it is a ferric hydroperoxy heme,9,10 rather than the more common ferryl heme, that is the reactive species in converting heme into its first intermediate, meso-hydroxy heme. The various HOs exhibit not only significant sequence but also structural, homology, as evidenced in crystal structures for HOs from mammals,11,12 photosynthetic bacteria,13 and pathogenic bacteria.14-18 The selectivity of the ring cleavage on the R-meso position in mammalian and plant HOs, and for (9) Ortiz de Montellano, P. R. Acc. Chem. Res. 1998, 31, 543-549. (10) Davydov, R. M.; Yoshida, T.; Ikeda-Saito, M.; Hoffman, B. M. J. Am. Chem. Soc. 1999, 121, 10656-10657. (11) Lad, L.; Schuller, D. J.; Shimizu, H.; Friedman, J.; Li, H.; Ortiz de Montellano, P. R.; Poulos, T. L. J. Biol. Chem. 2003, 278, 7834-7843. (12) Sugishima, M.; Sakamoto, H.; Noguchi, M.; Fukuyama, K. Biochemistry 2003, 42, 9898-9905. (13) Sugishima, M.; Migita, C. T.; Zhang, X.; Yoshida, T.; Fukuyama, K. Eur. J. Biochem. 2004, 271, 4517-4525. (14) Schuller, D. J.; Zhu, W.; Stojiljkovic, I.; Wilks, A.; Poulos, T. L. Biochemistry 2001, 40, 11552-11558. (15) Friedman, J. M.; Lad, L.; Deshmukh, R.; Li, H. Y.; Wilks, A.; Poulos, T. L. J. Biol. Chem. 2003, 278, 34654-34659. (16) Hirotsu, S.; Chu, G. C.; Unno, M.; Lee, D.-S.; Yoshida, T.; Park, S.-Y.; Shiro, Y.; Ikeda-Saito, M. J. Biol. Chem. 2004, 279, 11937-11947. J. AM. CHEM. SOC. 2006, 128, 6391-6399

9

6391

Ma et al.

ARTICLES

all but one pathogenic bacterium,19 is controlled both by blocking three meso positions to attack by the close placement of the distal helix and by steric “tilting”/directing the exogenous ligand (and presumably, Fe3+-OOH) toward the fourth meso position by the distal helix backbone. The stabilization of the unusual reactive ferric hydroperoxy species has been attributed to three largely conserved ordered water molecules,6,15 one of which provides the primary H-bond donor/acceptor to the nascent Fe2+O2 and Fe3+-OOH species. Solution 1H NMR has shown20-26 that HOs possess a distal H-bond network, with several remarkably strong H-bonds, that extends from the distal pocket to the HO surface and that this network includes the three proposed “catalytic” water molecules, as well as numerous additional water21-23,25 molecules ordered within this network. Sequence comparison among the structurally and/or spectroscopically investigated HOs shows that human HO isozyme #1 (hHO), rat HO (rHO), and the bacterial HO from Corynebacterium diphtheriae (CdHO) possess homologous H-bond networks which, due to significant substitution, differ from those of the bacteria from Neisseria meningitidis (NmHO) or Pseudomonas aeruginosa (PaHO), which, in turn, are more similar to each other.4,6 While each of the HOs exhibit three similarly, but not identically, placed water molecules in the distal pocket, it is not yet known how strongly these H-bond networks and distal water molecules are coupled to the iron center and how this network may differentially respond to the generation of the reactive species and/or intermediates in the HO reaction. Sizable movement of heteroatoms of both water and distal residues, due to imposed perturbation introduced by mutagenesis of either the HO or substrate, can be readily detected in X-ray crystallography, but changes in H-bond strength of the ordered water/H-bond network may be unresolved and would have to rely on 1H NMR spectroscopy,21,22,25-28 where low-field bias of the labile proton chemical shift correlates with increased H-bond donor strengths. A particularly suitable HO complex for investigating distal H-bond interaction between a ligated and the distal nonligated water molecules is the high-spin, ferric, resting state complex, HO-PH-H2O (PH ) protochemin). This complex, like similar metglobin complexes,29 participates in an acid S alkaline tran(17) Friedman, J.; Lad, L.; Li, H.; Wilks, A.; Poulos, T. L. Biochemistry 2004, 43, 5239-5245. (18) Unno, M.; Matsui, T.; Chu, G. C.; Coutoure, M.; Yoshida, T.; Rousseau, D. L.; Olson, J. S.; Ikeda-Saito, M. J. Biol. Chem. 2004, 279, 2105521061. (19) Caignan, G. A.; Deshmukh, R.; Wilks, A.; Zeng, Y.; Huang, H.-w.; MoenneLoccoz, P.; Bunce, R. A.; Eastman, M. A.; Rivera, M. J. Am. Chem. Soc. 2002, 124, 14879-14892. (20) La Mar, G. N.; Asokan, A.; Espiritu, B.; Yeh, D. C.; Auclair, K.; Ortiz de Montellano, P. R. J. Biol. Chem. 2001, 276, 15676-15687. (21) Li, Y.; Syvitski, R. T.; Auclair, K.; Wilks, A.; Ortiz de Montellano, P. R.; La Mar, G. N. J. Biol. Chem. 2002, 277, 33018-33031. (22) Li, Y.; Syvitski, R. T.; Auclair, K.; Ortiz de Montellano, P. R.; La Mar, G. N. J. Am. Chem. Soc. 2003, 125, 13392-13403. (23) Syvitski, R. T.; Li, Y.; Auclair, K.; Ortiz de Montellano, P. R.; La Mar, G. N. J. Am. Chem. Soc. 2002, 124, 14296-14297. (24) Li, Y.; Syvitski, R. T.; Chu, G. C.; Ikeda-Saito, M.; La Mar, G. N. J. Biol. Chem. 2003, 279, 6651-6663. (25) Liu, Y.; Zhang, X.; Yoshida, T.; La Mar, G. N. Biochemistry 2004, 43, 10112-10126. (26) Liu, Y.; Zhang, X.; Yoshida, T.; La Mar, G. N. J. Am. Chem. Soc. 2005, 127, 6409-6422. (27) Wagner, G.; Pardi, A.; Wu¨thrich, K. J. Am. Chem. Soc. 1983, 105, 59485949. (28) Harris, T. K.; Mildvan, A. S. Proteins: Struct., Funct., Genet. 1999, 35, 275-282. (29) Antonini, E.; Brunori, M. Hemoglobin and Myoglobin and Their Reactions with Ligands; Elsevier, North-Holland Publishing: Amsterdam, 1971; pp 40-54. 6392 J. AM. CHEM. SOC.

9

VOL. 128, NO. 19, 2006

Figure 1. Depiction of a ligated water a (A) in equilibrium with an ironligated hydroxide (B). In the absence of H-bond donation, ligated water a is a weak-field ligand. Upon serving as an H-bond donor to a proteinbased (or ordered water molecule b) acceptor, the axial field strength increases until in the limit of complete transfer, as in (B), resulting in the strong-field hydroxide ligand.

sition30-35 where the axial water molecule is deprotonated and the high-spin complex is converted to a largely low-spin, hydroxide complex (see Figure 1). Differential stabilization of the H-bonding of ligated water in different HOs is indicated by the variable pK in the range 8.0 to 9.0 in WT HOs30-34 and their sensitivity to mutation for a given HO.36 However, the altered pK does not allow one to determine whether a reduced pK is due to destabilization of the H2O or stabilization of the OHby H-bonding.37 Hence it is desirable to explore the response of the distal H-bonding of the ligated water upon perturbation in the interaction with a given HO-PH-H2O complex. Both optical38,39 and NMR40 spectroscopy have demonstrated that the pK for H2O S OH- transition in metglobins, as depicted in Figure 1, can be strongly modulated by peripheral substituents on the hemin, with electron-withdrawing (donating) substituents shifting the pK to lower (higher) pH values. Hence the ligated water H-bond donor, or ligated hydroxide H-bond acceptor, properties are readily perturbed by hemin substitution. In the absence of distal interaction, a ligated water (water a in Figure 1a) is a weak-field ligand. Serving as a H-bond donor to a nonligated water b (or distal base) will increase the axial field strength, AFS, of the ligated water, which, in the limit of complete deprotonation (Figure 1B), yields the strong-field hydroxide ligand. The axial field strength, AFS, for a high-spin, ferric hemin is directly reflected41 in the zero-field splitting parameter D, and hence the AFS, with the value of D increasing in magnitude, decreases. The D values for hemins are too large (5-20 cm-1) to detect by any conventional EPR41 but are readily observed in the 1H NMR of high-spin ferric hemoproteins20,42-46 through the dipolar (30) Sun, J.; Wilks, A.; Ortiz de Montellano, P. R.; Loehr, T. M. Biochemistry 1993, 32, 14151-14157. (31) Takahashi, S.; Wang, J. L.; Rousseau, D. L.; Ishikawa, K.; Yoshida, T.; Takeuchi, N.; Ikeda-Saito, M. Biochemistry 1994, 33, 5531-5538. (32) Chu, G. C.; Tomita, T.; So¨nnichsen, F. D.; Yoshida, T.; Ikeda-Saito, M. J. Biol. Chem. 1999, 274, 24490-24496. (33) Chu, G. C.; Katakura, K.; Zhang, X.; Yoshida, T.; Ikeda-Saito, M. J. Biol. Chem. 1999, 274, 21319-21325. (34) Ratliff, M.; Zhu, M.; Deshmukh, R.; Wilks, A.; Stojiljkovic, I. J. Bacteriol. 2001, 183, 6394-6403. (35) Caignan, G. A.; Deshmukh, R.; Zeng, Y.; Wilks, A.; Bunce, R. A.; Rivera, M. J. Am. Chem. Soc. 2003, 125, 11842-11852. (36) Zeng, Y.; Deshmukh, R.; Caignan, G. A.; Bunce, R. A.; Rivera, M.; Wilks, A. Biochemistry 2004, 43, 5222-5238. (37) Caughey, W. S. In Hemes and Hemoproteins; Chance, B., Estabrook, R. W., Yonetani, T., Eds.; Academic Press: New York, 1966; pp 276-277. (38) Cowgill, R. W.; Clark, W. M. J. Biol. Chem. 1952, 198, 33-61. (39) Brunori, M.; Amiconi, G.; Antonini, E.; Wyman, J.; Zito, R.; Rossi-Fanolli, A. Biochim. Biophys. Acta 1968, 154, 315-322. (40) McGrath, T. M.; La Mar, G. N. Biochim. Biophys. Acta 1978, 534, 99111. (41) Brackett, G. C.; Richards, D. L.; Caughey, W. S. J. Chem. Phys. 1971, 54, 4383-4401. (42) Kao, Y.-H.; Lecomte, J. T. J. J. Am. Chem. Soc. 1993, 115, 9754-9762.

Ligated Water H-Bonding in N. meningitidis Heme Oxygenase

ARTICLES

Figure 2. Structure of the hemin substrates with R2 ) R4 ) CH3 (2,4dimethyldeuterohemin, DMDH); R2 ) CH3, R4 ) vinyl (2-methyl-4vinyldeuterohemin, 2M4VDH); R2 ) R4 ) vinyl (protohemin); and R2 ) formyl, R4 ) vinyl (2-formyl-4-vinyldeuterohemin; 2F4VDH) with labeling by the Fischer notation, as viewed from the proximal side of NmHO. Also shown is the position of the axial His23 imidazole plane. The reference coordinate system, x′, y′, z′ has x′, y′ in the mean heme plane and passing through NB, NA, respectively, with z′ pointing to the proximal side. The magnetic axes, x, y, z, where the axially symmetric χ is diagonal, are defined by the tilt of the major magnetic axis, z, from the heme normal, z′, by the angle β (not shown), and the direction of the tilt of the major magnetic axis z is defined by angle R, between projection of z on the x′, y′ plane and the x′ axis. For the present complex, ∆χrh ) 0 and κ becomes irrelevant.

To assess the prospects for characterizing the influence of point mutation on the distal H-bond network, we initiate here the study of the influence of “mutation” of the substrate periphery on the NMR spectral properties of high-spin ferric NmHO-PH-H2O26 (where PH is protohemin, i.e., Figure 2 with R2 ) R4 ) vinyl). We utilize three modified hemins, in addition to native protohemin, PH (R2 ) R4 ) vinyl); 2,4-dimethyldeuterohemin, DMDH (R2 ) R4 ) methyl in Figure 2); 2-methyl4-vinyldeuterohemin, 2M4VDH (R2 ) methyl, R4 ) vinyl in Figure 2); and 2-formyl-4-vinyldeuterohemin, 2F4VDH (R2 ) formyl, R4 ) vinyl in Figure 2; also called Spirographis hemin). Our initial investigation focuses on the questions: Can changes in ∆χax(D) with modified hemin be detected? How well can relative ∆χax(D) values be quantitated? How far, and how selectively, do the altered H-bond properties of the axial water, reflected in ∆χax(D) values, propagate into the distal H-bond network? Experimental Section

with B ) Bohr magneton, k ) Boltzmann constant, NA ) Avogadro’s number, g ) spectroscopic splitting factor ) ∼2.0, and θ′, R are the coordinates of a proton in a convenient reference coordinate system (x′, y′, z′ in Figure 2). Γ(R, β, γ) is the Euler rotation that provides the orientation of the unique axis in the molecular framework. The overwhelming dominance of the zero-field splitting contribution to ∆χ, manifested42,43,45-48 in strictly T-2 dependence of δdip(obsd) and the axial anisotropy, have been established26 for the NmHO-PH-H2O complex, as well as several other high-spin ferric hemoproteins.42,43,45 The necessary NMR protocols for definitively assigning the strongly relaxed dipolar-shifted residues26,48 and the conversion of these dipolar shifts to a quantitative ∆χax and D have been presented.

Enzyme Preparation: The apo-NmHO samples used in this study are the same as those described in detail previously.25 Stoichiometric amounts of the four hemins of interest, dissolved in 0.1 M KOH in 1 H2O, were added to apo-NmHO in phosphate buffer (100 mM, pH 6.8-7.2). The three modified hemins used (in addition to protohemin (R2 ) R4 ) vinyl in Figure 2), 2,4-dimethyldeuteriohemin, DMDH (R2 ) R4 ) methyl in Figure 2), 2-methyl-4-vinyldeuteriohemin, 2M4VDH (R2 ) methyl, R4 ) vinyl in Figure 2), and Spirographis or 2-formyl-4-vinyldeuteriohemin, 2F4VDH (R2 ) formyl, R4 ) vinyl in Figure 2), are the same as those reported previously.21,49 The substrate complexes were purified by column chromatography on Sephadex G25 and yielded samples ∼3 mM in NmHO-PH-H2O at pH 6.8-7.2. Samples in 1H2O were converted to 2H2O by column chromatography.50 NMR Spectroscopy: 1H NMR data were collected on Bruker AVANCE 500 and 600 spectrometers operating at 500 and 600 MHz, respectively. Reference spectra at 500 MHz were collected in 1H2O and 2H2O over the temperature range 10-35 °C at repetition rates of 1 s-1 (14 and 22 ppm bandwidths) and 5 s-1 (200 ppm bandwidth) using a standard one-pulse sequence, as well as with 3:9:19 “soft” pulse detection51 with and without saturation of the water solvent signal. Chemical shifts are referenced to 2,2-dimethyl-2-silapentane-5-sulfonate (DSS) through the water resonance calibrated at each temperature. Nonselective T1s were estimated in both 1H2O and 2H2O at 25° from the null point of a standard inversion-recovery pulse sequence. 600 MHz NOESY 52 spectra (mixing time, 40 ms) and 500 MHz CleanTOCSY (to suppress ROESY response53) spectra (spin lock 15 and 25 ms) using MLEV-1754 were recorded over a bandwidth of 15 kHz (NOESY) and 12 kHz (TOCSY) with recycle rates of 1-5 s-1, using 512 t1 blocks of 256 or 512 scans, each consisting of 2048 t2 points. 2D data sets were processed using Bruker XWIN software on a Silicon Graphics Indigo workstation and consisted of 30°- or 45°-sine-squaredbell-apodization in both dimensions and zero-filling to 2048 × 2048 data points prior to Fourier transformation. Magnetic Axes Determination: The locations of the magnetic axes were determined by finding the Euler rotation angles, Γ(R,β,γ), that rotate the crystal-structure-based, iron-centered reference coordinate

(43) Clark, K.; Dugad, L. B.; Bartsch, R. G.; Cusanovich, M. A.; La Mar, G. N. J. Am. Chem. Soc. 1996, 118, 4654-4664. (44) De´me´ne´, H.; Tsan, P.; Gans, P.; Marion, D. J. Phys. Chem. B 2000, 104, 2559-2569. (45) Asokan, A.; de Ropp, J. S.; Newmyer, S. L.; Ortiz de Montellano, P. R. J. Am. Chem. Soc. 2001, 123, 4243-4254. (46) Bertini, I.; Luchinat, C. Coord. Chem. ReV. 1996, 150, 1-296. (47) Kurland, R. J.; McGarvey, B. R. J. Magn. Reson. 1970, 2, 286-301. (48) La Mar, G. N.; Satterlee, J. D.; de Ropp, J. S. In The Porphyrins Handbook; Kadish, K. M., Smith, K. M., Guilard, R., Eds.; Academic Press: San Diego, 2000; Vol. 5, pp 185-298.

(49) Kolczak, U.; Hauksson, J. B.; Davis, N. L.; Pande, U.; de Ropp, J. S.; Langry, K. C.; Smith, K. M.; La Mar, G. N. J. Am. Chem. Soc. 1999, 121, 835-843. (50) Johnston, P. D.; Figueroa, N.; Redfield, A. G. Proc. Natl. Acad. Sci. U.S.A. 1979, 76, 3130-3134. (51) Piotto, M.; Sandek, V.; Sklenar, V. J. Biomol. NMR 1992, 2, 661-666. (52) Jeener, J.; Meier, B. H.; Bachmann, P.; Ernst, R. R. J. Chem. Phys. 1979, 71, 4546-4553. (53) Griesinger, C.; Otting, G.; Wu¨thrich, K.; Ernst, R. R. J. Am. Chem. Soc. 1988, 110, 7870-7872. (54) Bax, A.; Davis, D. G. J. Magn. Reson. 1985, 65, 355-360.

shifts for nonligated residues. With the isotopic g ) 2 value for ferric hemins, and the strong dominance of axial symmetry,20,42-46 the dipolar shifts for nonligated residue protons are given by

δdip ) -(12µoπNA)-1∆χax(3 cos2 θ′ - 1)R-3Γ(R,β,γ) (1) where ∆χax is given by46-48

∆χax ) -

21g2B2NA D/T2, 5k2

(2)

J. AM. CHEM. SOC.

9

VOL. 128, NO. 19, 2006 6393

Ma et al.

ARTICLES system, x′, y′, z′, into the magnetic coordinate system, x, y, z, where the paramagnetic susceptibility tensor, χ, is diagonal, where R, β, γ are the three Euler angles.48,55,56 The angle β dictates the tilt of the major magnetic axis, z, from the heme normal z′; R reflects the direction of this tilt and is defined as the angle between the projection of the z axis on the heme plane and the x′ axis (Figure 2); and κ ≈ R + γ is the angle between the projection of the x, y axes onto the heme plane and locates the rhombic axes (Figure 2). In the present case, the tensor is axially symmetric,26,42-45 so that ∆χrh ) 0, and γ (κ) becomes irrelevant. The magnetic axes were determined by a least-squares search for the minimum in the error function, F/n:26,48,55,56 n

F/n )

∑|δ

dip(obsd)

- δdip(calcd)|2

(3)

i)1

where the calculated dipolar shift in the reference coordinate system, x′, y′, z′, or R, θ′, is given by eq 1 where ∆χax is the axial anisotropy of the diagonal, axially symmetric paramagnetic susceptibility tensor. The observed dipolar shift, δdip(obsd), is given by

δdip(obsd) ) δDSS(obsd) - δDSS(dia)

(4)

where δDSS(obsd) and δDSS(dia) are the chemical shifts, in ppm, referenced to DSS, for the paramagnetic NmHO-PH-H2O complex and an isostructural diamagnetic complex, respectively. In the absence of an experimental δDSS(dia) for nonlabile protons, it may be reliably estimated57,58 from the available molecular structure14 and available computer programs57,58 as described previously for both NmHO-PH-CN25 and NmHO-PH-H2O.26 The calculated δDSS(dia) are estimated to be reliable to (0.5 ppm. Since the errors are random, the uncertainties in δDSS(dia) appear primarily in the deviations from the straight line expected in a plot of δdip(obsd) vs δdip(calcd) and have only a small influence on the uncertainties in the δdip(calcd) resulting from an optimized set of magnetic axes and anisotropies. The uncertainties in δdip(calcd) are estimated at (15%.

Results and Discussion

Substrate Resonances: The resolved portions of the 500 MHz 1H NMR spectra of the NmHO-hemin-H2O complexes of DMDH, 2M4VDH, PH, and 2F4VDH are illustrated in Figure 3A, B, C, and D, respectively. In the low-field portion, the positions of the heme methyls are indicated by labels M (individually assigned26 only for PH, shown as Mi, i ) position on hemin in Figure 2). Heme methyl and HRs of propionates, vinyl and formyl substituents, in the 90-30 ppm window exhibit T1’s ≈ 4-6 ms, which are typical of high-spin ferrihemoproteins48 and similar to those reported previously for NmHO-PH-H2O.26 The upfield portion of the spectra reveals the 4-vinyl Hβs for the 2M4VDH and 2F4VDH complexes. The chemical shifts for heme substituents are tabulated in the Supporting Information. The variable spread of the largely contact-shifted pyrrole substituents among the four substrate complexes reflect primarily the spin density distribution due to the asymmetric positions of the variable electron-withdrawing vinyl and formyl groups.48,59 Similar spreads of the signals are observed for the chromophore outside the protein matrix,60 such (55) Williams, G.; Clayden, N. J.; Moore, G. R.; Williams, R. J. P. J. Mol. Biol. 1985, 183, 447-460. (56) Emerson, S. D.; La Mar, G. N. Biochemistry 1990, 29, 1556-1566. (57) Neal, S.; Nip, A. M.; Zhang, H.; Wishart, D. S. J. Biomol. NMR 2003, 26, 215-240. (58) Cross, K. J.; Wright, P. E. J. Magn. Reson. 1985, 64, 220-231. (59) Walker, F. A. In The Porphyrin Handbook; Kadish, K. M., Smith, K. M., Guilard, R., Eds.; Academic Press: Boston, 2000; Vol. 5, pp 1-183. (60) Budd, D. L.; La Mar, G. N.; Langry, K. C.; Smith, K. M.; Nayyir-Mazhir, R. J. Am. Chem. Soc. 1979, 101, 6091-6096. 6394 J. AM. CHEM. SOC.

9

VOL. 128, NO. 19, 2006

that the heme peaks provide little new information on the interaction of substrate with HO. The DMDH substrate outside the protein matrix exhibits a very small spread of the heme methyl shifts49,59 since all eight pyrrole substituents are alkyl groups, rendering the heme effectively axially symmetric. The small spread of these methyls in NmHO-DMDH-H2O attests to the relatively small rhombic perturbation that the NmHO matrix exerts on the heme contact shift pattern. Axial His23 Resonances: The resonance at ∼28 ppm in NmHO-PH-H2O has been assigned26 to the axial His CβH in NmHO-PH-H2O on the basis of the characteristically shorter T1 (∼2 ms) than that for pyrrole CRHs (4-6 ms).48 Comparison of the four NMR traces in Figure 3 shows that the hyperfine shifts for this His23 CβH peak exhibits only a very small variation among the four complexes. However, the contact shift for this proton obtained via

δcon ) δDSS(obsd) - δDSS(dia) - δdip(calcd),

(5)

leads to δcon values of 35.6, 36.2, 35.2, and 32.9 ppm for the DMDH, 2M4VDH, PH, and 2F4VDH complexes, respectively. While even the δcon shifts are quite similar, it does appear that the Fe-His23 covalency is somewhat reduced upon introducing the strongly electron-withdrawing formyl substituent in the NmHO-2F4VDH-H2O complex. Such a response is consistent with the results obtained for ∆χax values below. Residue Assignments: The upfield portions of the spectra in Figure 3 exhibit one well-resolved methyl peak previously assigned26 to Ala121 CβH3 in the NmHO-PH-H2O complex. Comparison of the four complexes in Figure 3 shows that the upfield bias, and hence the magnitude of the upfield dipolar shift, decreases monotonically in the order DMDH > 2M4VDH > PH > 2F4VDH. Comparison of NOESY (see Figure 4) and TOCSY (not shown; see Supporting Information) maps among the four complexes reveals that all protons with significant dipolar shifts exhibit the same pattern of contacts among each other, dictating that active site molecular structure is conserved. Hence all residues assigned in detail for NmHO-PH-H2O are very straightforwardly transferred to each of the new target complexes using the identical 2D data described in detail for NmHO-PH-H2O.26 The remainder of the chemical shift data are provided in the Supporting Information. The chemical shifts for residues in the four complexes with at least one proton with δdip > |1.0| ppm are listed in Table 1. It is noted in Table 1 that all resonances with a significant dipolar shift exhibit monotonic decreases of δdip(obsd) in the same order as the Ala121 CβH3 signal in Figure 3, dictating that the electronic structures (magnetic properties) are systematically perturbed. This systematic effect on dipolar shifts is wellillustrated in the high-field resolved portion of the two complexes NmHO-2F4VDH-H2O (Figure 4A) and NmHODMDH-H2O (Figure 4B), the substrates with the extreme difference in the degree of electron-withdrawing power of the substituents. The trend is also obvious in the chemical shifts in Table 1. The NOESY connectivities among the upfield shifted residues in NmHO-DMDH-H2O is shown in Figure 4C-F. The lowfield portion of the NOESY spectra for NmHO-DMDH-H2O (not shown, see Supporting Information) reveals the conserved pattern of NOESY cross-peaks among the labile proton involved

Ligated Water H-Bonding in N. meningitidis Heme Oxygenase

ARTICLES

Figure 3. Resolved portions of the 500 MHz 1H NMR spectra of high-spin ferric resting state complexes in 2H2O, 100 mM in phosphate at 25 °C: (A) NmHO-DMDH-2H2O at pH 7.2; (B) NmHO-2M4VDH-2H2O at pH 7.0; (C) NmHO-PH-2H2O at pH 7.0; and (D) NmHO-2F4VDH-2H2O at pH 6.8. Heme methyl peaks are labeled M, and single protons, H, with a subscript in the Fischer notation if the peak has been assigned.26 Amino acid residue signals are identified by the residue number and position. Peak a was assigned to Ala121 CRH on the basis of the NOE to Ala121 CβH3 in NmHO-PH-H2O.26 Peak b (T1 ∼3 ms) in NmHO-PH-H2O was not detectable because it resonates under the more weakly relaxed vinyl 2Hβ peak. The presence of peak b is made obvious in the DMDH complex in part A; its T1 ∼3 ms and shift are consistent with prediction26 for Ser117 CRH. Nevertheless, both peaks a and b move to lower field in the same order as the Ala121 CβH3 (M121).

in the robust H-bonds in the distal pocket relative to that reported26 for NmHO-PH-H2O. Similar data for NmHO2F4VDH-H2O are provided in the Supporting Information. The NOESY cross-peak pattern for the Gln49, His53 side chains to Val81 and Leu142 (see Supporting Information), as previously found for NmHO-PH-CN25 and NmHO-PH-H2O,26 dictates that the Gln49 and His53 side chains in solution exhibit 180° rotation about the β-γ bond relative to those described in either crystal structure.14,15 The chemical shift data for Gln49 and His53 are included in Table 1. The data in Table 1 demonstrate clearly that even subtle perturbations such as changing the electron-withdrawing characteristic of a substrate substituent leads to detectable and systematic changes in the dipolar shift pattern for the structurally conserved protein matrix. The chemical shifts for residues with a proton with δdip > |1.0| that is assigned in at least the DMDH, PH, and 2F4VDH complexes are listed in Table 1 where we also list the δdip(calcd) for the previously reported magnetic axes of NmHO-PH-H2O. While the majority of the same residues were assigned in all four complexes, the systematically much smaller δdip for the 2FV4DH than the PH or DMDH complexes resulted in somewhat poorer resolution and, hence, provided slightly fewer definitive assignments. Similar data on other assigned residues are provided in the Supporting Information. Magnetic Axes Determination: Inspection of the chemical shifts in Table 1 shows that observed shift changes in direction (up- or downfield) and magnitude with substrate substituent (basicity) are approximately proportional to the sign and magnitude of the calculated dipolar shift for the NmHOPH-H2O complex (given for PH in the last column of Table 1), with the inferred dipolar shift magnitude increasing in the

order 2F4VDH < PH < 2M4VDH < DMDH. Hence these chemical shift changes must reflect some systematic variation in the anisotropy and/or orientation of χ with substrate. This trend is more quantitatively assessed in the magnetic anisotropies obtained by minimizing the error function in eq 3. The χ tensor for NmHO-PH-H2O was shown26 to be axial, with ∆χax ) (-2.52 ( 0.07) × 10-8 m3/mol, tilted from the heme normal by β ) 7 ( 1°, in a direction defined by R ) 210 ( 10° in Figure 2. The tilt magnitude and direction were shown to be the same as that for the Fe-N(His23) vector in the NmHO-PH-H2O crystal,14 confirming the reasonable expectation that the major magnetic axis is defined by the axial bond. Three parameter searches for the four complexes, using a variety of the available δdip(obsd), yield very good correlation between δdip(obsd) and δdip(calcd) (see Figure 5) and result in the following (∆χax, β, R) for DMDH, 2M4VDH, PH, and 2F4VDH complexes (∆χax ) (-2.71 ( 0.08) × 10-8 m3/mol, β ) 9 ( 1°, R ) 195 ( 10°), (∆χax ) (-2.59 ( 0.08) × 10-8 m3/mol, β ) 7 ( 1°, R ) 200 ( 10°), ((∆χax ) -2.52 ( 0.08) × 10-8 m3/mol, β ) 7 ( 1°, R ) 210 ( 10°), and (∆χax ) (-1.98 ( 0.09) × 10-8 m3/mol, β ) 10 ( 2°, R ) 200 ( 15°). Inspection of these resulting parameters reveals that the orientation of the χ tensor is, in fact, essentially conserved, with a mean orientation defined by a tilt β ) 8.5 ( 1.5° in a direction R ) 200 ( 15°. When we restrict R and β to the values determined by the axial His23 N-Fe vector, as determined for NmHO-PH-H2O,26 the one-parameter search results in the relative ∆χax value -2.64/-2.58/-2.52/-2.02, in 10-8 m3/mol, each with (0.08 × 10-8 m3/mol uncertainties for DMDH/2M4VDH/PH/2F4VDH, respectively. The individual values of ∆χ vary by as much as 4% upon selecting different input data, indicating that any single J. AM. CHEM. SOC.

9

VOL. 128, NO. 19, 2006 6395

Ma et al.

ARTICLES

Table 1. Chemical Shift Data for Strongly Dipolar Shifted Residues in NmHO-Hemin-H2O Complexes as a Function of Hemin Substituentsa

Thr20 Ala21 Val22

His23 Asp24 Gln49 Phe52 His53 Tyr112

Asn118 Ala121 Phe125

Leu142 Figure 4. Resolved upfield portions of the 600 MHz 1H NMR spectrum, in 2H2O, 100 mM phosphate at 25 °C of (A) NmHO-2F4VDH-2H2O at pH 6.8 and (B) NmHO-DMDH-2H2O at pH 7.2 illustrating the strong bias to lower field (smaller upfield dipolar shift) as electron-withdrawing substituents are introduced into the heme periphery. Portions of the 600 MHz NOESY spectrum (mixing time 40 ms, repetition rate 2 s-1) for NmHO-DMDH-2H2O illustrate (C) intraresidue contacts for Val22, Asp24, and Leu142; (D) intra-residue contacts for Thr20 and Ala21 and inter-residue contacts between Ala21 and Asp24, Thr20 and Ala21; (E) inter-residue contact between Leu142 and His141, and (F) inter-residue contacts for Leu142.

∆χax will not be better defined than to (0.08 × 10-8 m3/mol. However, while individual values can vary by ∼5%, the relatiVe ∆χax(D) values among the four complexes for any one set of input data (normalized to PH) are remarkably reproducible as 1.08 ( 0.01/1.04 ( 0.01/1:00/0.82 ( 0.2 for the DMDH/ 2M4VDH/PH/2F4VDH complexes. Thus relative values are much more better defined than the absolute values of ∆χax. With D ) 9.1 cm-1 reported for the NmHO-PH-H2O complex,26 the D value for the complexes with DMDH, 2M4VDH, and 2F4VDH are 9.8, 9.5, and 7.5 cm-1, respectively. On the basis of these results, we can conclude that, since the axial field strength inversely correlates with D(∆χax magnitude41), the AFS varies in the order 2F4VDH > PH > 2M4VDH > DMDH. The strongest AFS arising from the substrate with the most electron-withdrawing substituents is completely consistent with the expectation of such a perturbation in the H-bond 6396 J. AM. CHEM. SOC.

9

VOL. 128, NO. 19, 2006

CβH CγH3 CRH CβH3 NH CRH CβH Cγ1H3 Cγ2H3 CβH Cβ1H Cβ2H Nη1H Nη2H CδHs CHs CζH Cδ1H C2H N1H NH CR H Cβ1H Cβ2H CδHs CHs OηH Cβ1H Cβ2H CβH3 CRH Cβ1H Cβ2H CδHs CHs CζH NH CRH Cβ1H Cβ2H CγH Cδ1H3 Cδ2H3

DMDHb

2M4VDHc

PHd

2F4VDHe

2.52 -0.83 1.78 -0.04 6.20 2.04 0.50 -0.75 -0.23 27.6 -0.53 -0.80 8.95 5.82 7.80 8.25 8.68 6.92 6.86 10.72 9.18 5.31 4.40 3.98 8.35 8.90 14.6 0.33 0.15 -3.74 3.16 1.85 1.90 6.28 6.55 6.72 5.87 2.04 -0.15 0.55 -1.03 0.11 -1.31

2.61 -0.75 1.85 0.04 6.21 2.11 0.57 -0.65 -0.17 28.6 -0.50 -0.75 8.99 6.06

2.69g -0.60g 2.03g 0.15g 6.36 2.19 0.65 -0.63 -0.12 27.8 -0.37 -0.60 9.18 6.21 7.76 8.15 8.55 7.01 7.06 10.92 9.18 5.30 4.38 4.00 8.30 8.82 14.8 -0.41 -0.26 -3.14 3.25 1.92 2.03 6.37 6.63 6.78 6.03 2.32 0.10 0.70 -0.77 0.33 -0.85

2.60 -0.43 2.32 0.32 6.68 2.55 0.87 -0.38 0.07 26.9 0.03 -0.22 9.27 6.38 7.70 8.07 8.46

6.99 6.99 10.82 9.18 5.29 4.39 3.99 8.34 14.7 -0.36 -0.18 -3.41 3.18 1.86 1.92 6.28 6.56 6.71 5.93 2.19 -0.06 0.61 -0.91 0.18 -1.13

11.10 9.11 5.20 4.28 3.88 8.20 14.0 -2.35 3.35 2.09 2.20 6.46 6.67 6.11 2.60 0.25 0.83 -0.63 0.36 -0.70

δdip(calcd), PHf

-1.22 -2.04 -1.57 -0.93 -1.38 -1.25 -0.77 -0.91 -0.62 -8.4 -2.20 -2.51 -0.58 -0.83 0.46 0.70 1.02 -0.11 -0.60 -0.56 0.50 0.71 0.66 0.58 0.87 1.06 1.56 -1.98 -1.86 -4.50 -0.53 -0.66 -1.00 -0.54 -0.39 -0.35 -0.69 -1.51 -1.16 -0.76 -0.77 -0.53 -1.45

a

Observed chemical shift, δDSS(obsd), in ppm, referenced to DSS, in solution, 100 mM in phosphate, pH ∼7.0 at 25 °C. b 2,4-Dimethyldeuterohemin (R2 ) R4 ) methyl) in Figure 2. c 2-Methyl-4-vinyldeuterohemin (R2 ) methyl, R4 ) vinyl) in Figure 2. d Protohemin (R2 ) R4 ) vinyl) in Figure 2. e 2-Formyl-4-vinyldeuterohemin (spirographis hemin) (R2 ) formyl, R4 ) vinyl) in Figure 2. f Dipolar shift, in ppm, calculated for the NmHO-PH-H2O complex with eq 1 and the reported magnetic axes.26 g Corrections to values reported in ref 26. 1H O 2

donor strength of the axial water,37 as observed in the decrease in the pK value for the H2O > OH- conversion in metglobins.38-40 The present data, however, detect only the changes in the stabilization of the aquo complex. The trend in ∆χ correlates well with the Hammet61 σp substituent parameters that predict the basicity of the substrate, as shown in Figure 6. The fact that the ∆χ(D) values for the 2M4VDH complex are midway between those of DMDH and PH complexes demonstrates that the two vinyl groups contribute equally to the basicity of PH relative to DMDH. The strongest AFS for the 2F4VDH complex is consistent with the observation that the axial His23 CβH contact shift, and hence Fe-His covalency, is slightly reduced relative to the other three substrates, as expected from the trans(61) Hansch, C.; Leo, A.; Taft, R. W. Chem. ReV. 1991, 91, 165-195

Ligated Water H-Bonding in N. meningitidis Heme Oxygenase

ARTICLES

Figure 6. Plot of the determined axial magnetic anisotropy, ∆χax, versus the sum of the Hammett σp constant61 for substituents at positions 2 and 4 that would be related to the basicity of the porphyrin for the four substrates in high-spin ferric, resting state NmHO complexes.

Figure 5. Plot of δdip(obsd) versus δdip(calcd) for the optimized magnetic axes of (A) NmHO-DMDH-H2O (∆χax ) (-2.64 ( 0.08) × 10-8 m3/mol), (B) NmHO-2M4VDH-H2O, with ∆χax ) (-2.58 ( 0.07) × 10-8 m3/mol, (C) NmHO-PH-H2O (∆χax ) (-2.52 ( 0.08) × 10-8 m3/ mol), and (D) NmHO-2F4VDH-H2O (∆χax ) (-2.02 ( 0.09) × 10-8 m3/mol) resulting from a one-parameter search (∆χax) for the minimum in eq 3 with the χ oriented through the axial His23 N (R ) 210°, β ) 7°) as found for NmHO-PH-H2O.26

effect.62 It is expected that, as additional assignments become available on 15N- and/or 15N/13C-labeled NmHO complexes, the uncertainties in both absolute and relative values of ∆χax can be reduced. However, the uncertainties in the individual ∆χ(D) values will always be greater than those in the relative values for a given HO with a variety of substituents. H-Bond Strengths: It is noted that, because of variable H-bonding, δDSS(dia) for labile protons cannot be reliably estimated by the shift X program.25,57 Hence, comparison of H-bond strength among the four substrates requires the determination of δDSS(dia*), the chemical shifts for a labile proton in the absence of paramagnetism, which we define by

δDSS(dia*) ) δDSS(obsd) - δdip(calcd)

(6)

where δdip(calcd) for the different substrates is obtained for each set of magnetic axes. A list of the δDSS(dia*) for labile protons for peptide NHs and important side chains that depend sig-

nificantly on substituents is provided in Table 2. Also included for comparison are the range of δDSS(obsd) and ∆δDSS(obsd) and the range of the calculated ∆δDSS(dia*), as well as the δdip(calcd) for the PH complex.26 The X-Y distance in the X-H‚‚‚Y H-bond has been shown27 to correlate with the extent of low-field bias on the labile proton chemical shift, δDSS(dia*), and this distance inversely relates to H-bond strength.28 Hence increased low-field bias of labile proton chemical shift can be taken as evidence for increased H-bond strength. The reliability of the determined δDSS(dia*) is strongly supported by the fact that for the numerous peptide NHs which exhibit significantly different chemical shifts (∆δDSS(obsd*) in Table 2) with 2,4-substitution in Table 2, correction for δdip(calcd) shows that δDSS(dia*), and hence H-bonding, is essentially independent of substrates for all but one notable exception, Glu115. The Glu115 δDSS(dia*) argues for increased H-bond strength in the order DMDH < 2M4VDH < PH < 2F4VDH. The Glu115 side chain serves as an acceptor to the side chain of one of two residues, Asn118 or Arg77, with the former favored for the H-bond donor H2O ligand, and the latter favored for the H-bond acceptor NO ligands,14,15 and a change of the position of the equilibrium could selectively perturb the Glu115 NH H-bond strength. A similar anomalous difference in the Glu115 NH H-bond strength was detected for the NmHO-PH-CN25 and NmHO-PH-H2O26 complexes. Among the side chain NHs, moderate systematic differences in δDSS(obsd) for the His141 Nδ1H and Trp153 N1H are removed by correcting for δdip(calcd), to demonstrate completely conserved H-bond strength. For the Gln49 and His53 side chain, however, significant systematic differences in δDSS(obsd) are only partially decreased in δDSS(dia*), suggesting slightly increasing H-bond strength for both side chains upon introduction of electron-withdrawing substituents in the substrate. The position of three catalytically implicated, ordered water mol(62) Buchler, J. W.; Kokisch, W.; Smith, P. D. Struct. Bonding 1982, 34, 79132. J. AM. CHEM. SOC.

9

VOL. 128, NO. 19, 2006 6397

Ma et al.

ARTICLES Table 2. Calculated Diamagnetic Chemical Shifts, (δDSS(dia*)) for Labile Protons in High-Spin Ferric NmHO-Substrate-H2O Complexes as a Function of Substrate Substituentsa δDSS(dia*)

∆δDSS(obsd)f

∆δDSS- δdip(calcd)h (PH) (dia*)g

Peptide NHs 7.74 7.78 8.48 8.56 8.68 8.71 8.37 8.38 9.04 9.03 9.50 9.71 10.29 10.30 7.58 7.58 12.05 12.03 10.79 10.79 6.71 6.65 6.67 6.67 10.60 10.56 8.42 8.46 8.03 8.08 8.51 8.54

-0.52 0.28 0.11 0.30 0.16 -0.27 -0.09 -0.10 -0.09 -0.15 -0.24 -0.24 0.05 0.10 0.04 0.27

-0.12 0.06 -0.06 0.03 -0.06 -0.40 -0.05 -0.05 -0.02 0.04 -0.06 -0.04 0.12 -0.01 -0.01 0.03

-1.38 0.84 0.50 0.95 0.86 0.51 -0.14 -0.22 -0.28 -0.42 -0.69 -0.73 -0.27 0.37 0.37 0.54

Side Chain NHs 9.60 9.76 9.75 6.96 7.05 7.01 11.40 11.48 12.55 9.79 9.89 9.85 10.69 10.73 10.79 10.57 10.54 10.45

-0.32 -0.56 -0.38 0.01 -0.18 0.26

-0.17 -0.28 -0.23 0.01 -0.11 0.10

-0.58 -0.83 -0.56 -0.18 -0.36 0.89

DMDHb 2M4VDHc

Val22 Met31 Tyr112 Cys113 Ala114 Glu115 Gly138 Ala139 Arg140 His141 Leu142 Ala143 Asp147 Ala180 Phe181 Tyr184

7.66 8.62 8.65 8.35 8.97 9.31 10.25 7.53 12.01 10.75 6.60 6.63 10.68 8.45 8.01 8.57

Gln49 Nη1H 9.58 Gln49 Nη2H 6.73 His53 N1H 11.32 Arg140 NH 9.86 His141 Nδ1H 10.68 Trp153 N1H 10.55

7.63 8.58 8.67 8.35 9.02 9.41 10.27 7.56 12.02 70.76 6.65 6.63 10.66 8.44 8.02 -

PHd

2F4VDHe

a The diamagnetic chemical shifts, in ppm, referenced to DSS, as obtained via eqs 1 and 3 and with the magnetic axes for each complex described in Figure 6. b 2,4-Dimethyldeuterohemin (R2 ) R4 ) methyl) in Figure 2. c 2-Methyl-4-vinyldeuterohemin (R ) methyl, R ) vinyl) in Figure 2. 2 4 d Protohemin (R ) R ) vinyl) in Figure 2. e 2-Formyl-4-vinyldeutero2 4 hemin (spirographis hemin) (R2 ) formyl, R4 ) vinyl) in Figure 2. f The spread, with substrate, of the observed shift, ∆δDSS(obsd) ) δDSS(obsd: NmHO-DMDH-H2O) - δDSS(obsd: NmHO-2F4VDH-H2O), in ppm. g The spread, with substrate, of the calculated diamagnetic shift (eq 3), given by ∆δDSS(dia*) ) δDSS(dia*: NmHO-DMDH-H2O) - δDSS(dia*: NmHO-2F4VDH-H2O), in ppm. h Calculated δdip, in ppm, via eq 1 and the optimized magnetic axes, for the NmHO-PH-H2O complex.26

ecules (no. 37, 77, and 243 in the crystal14) (yellow spheres), the ligated water (green sphere), and residues that interact directly or indirectly with these water molecules are shown in Figure 7. The directions (H-donor f acceptor) of the H-bonds are shown as arrows in those cases where crystallography14 or solution NMR26 data provide information on the location of the proton. For other interactions where this information is not available, or ambiguous, likely H-bond interactions are shown by dashed lines. Since the ligated water no. 1 and the side chains of both Gln49 Nηs and His53 ring interact primarily as H-bond donors to ordered water molecules (Figure 7), the simplest expectation would be that increased H-bond donation of ligated H2O in the 2F4VDH (relative to the other heme) complex will result in reduced H-bond donation by the two side chains, contrary to observation. Such an expected behavior is clearly observed63 in the NmHO-PH-H2O conversion to NmHO-PH-OH. Two points may be relevant to this apparent paradox. One is that the Gln49 and His53 H-bonds and the three ordered water molecules depicted in Figure 7 are only a portion of a more extensive network of ordered water molecules and strong H-bonds, and the small changes in H-bond strength are the result of optimizing the network when increasing H-bond donor strength of the ligated water upon replacing DMDH by (63) Ma, L.-H.; Liu, Y.; Zhang, X.; Yoshida, T.; La Mar, G. N. J. Am. Chem. Soc. 2006, 128, in press. 6398 J. AM. CHEM. SOC.

9

VOL. 128, NO. 19, 2006

Figure 7. Portion of the crystal structure14 active site of NmHO-PHH2O that depicts the water molecule (no. 1 at light green sphere) ligated to the heme iron (red sphere), as well as the three catalytically implicated nonligated water molecules (no. 37, 77, and 243) (yellow spheres), and residues which directly (Gln49 (dark blue), His53 (brown), Ser117 (light blue), Ala121 (aqua), Tyr184 (magenta)) or indirectly (His141 (green), Asp84 (beige)) interact with one or more of these water molecules. Solid arrows indicate the known direction of the donor f acceptor in the implied H-bonds. The other H-bond interactions where the identities of the acceptor/ donor cannot be uniquely determined from crystallographic14 or NMR26 data are shown in dashed lines. It is noteworthy that the ligated water no. 1 and the two key distal H-bond network members Gln49 and His53 side chains all necessarily serve as donors to one or another of the three catalytically implicated water molecules.

2F4VDH. The second point is that, although the determined magnetic axes for the four complexes equally well account for the pattern of dipolar shift, the Gln49 and His53 side chains (more so the former) are placed relatively close to the node in the dipolar field, such that their δdip(calcd) values are more sensitive to small changes in structure (or uncertantities in δdip(calcd)) than those of the other residues. We have shown22 for hHO-PH-CN that peptide NH in strong H-bonds exhibit much larger 15N than 1H chemical shift changes, in ppm, upon altering H-bond strength. Since δdip(calcd) for the 15N is comparable to δdip(calcd) for 1H in a peptide NH, 15N will more unambiguously characterize the H-bond change for Gln49 and His53 than presently possible. Efficient routes to 15N labeling of NmHO are under investigation. The additional assignments accessed with 15N-labeled NmHO will also provide additional data that will more quantitatively describe the magnetic axes in the complex of interest. Conclusions

The present results demonstrate that changes in the H-bond properties in the ligated water are readily detected, and reasonably quantitated, by 1H NMR in high-spin, resting state HO-substrate complexes. The expected increase in ligated water H-bond donor strength, and hence increased effective axial field strength of the heme, with increased electron withdrawing power of substrate substituent, is observed in the decrease of the axial magnetic anisotropy (and zero-field splitting parameter D). Moreover, apparent influences of the altered H-bonding of the axial water appear to be transmitted to at least two residues

Ligated Water H-Bonding in N. meningitidis Heme Oxygenase

some 10 Å from the iron. His53 is the residue homologous to Tyr58 in mammalian HO whose mutagenesis abolishes PH cleavage.64,65 There is good reason to be optimistic that mutagenesis of the member residues of the distal H-bond network should lead to similarly perturbed ∆χ(D) values that would shed light on the link between the axial water molecule and the distal network of ordered water molecules and H-bonds. Logical extensions of these studies include the comparison of ∆χ(D) values among different HO complexes to assess differences in AFS and the relative sensitivity of ∆χax(D) on the same four substrates among complexes of a given HO. Such studies with CdHO and hHO are in progress. (64) Koenigs Lightning, L.; Huang, H.-W.; Moe¨nne-Loccoz, P.; Loehr, T. M.; Schuller, D. J.; Poulos, T. L.; Ortiz de Montellano, P. R. J. Biol. Chem. 2001, 276, 10612-10619. (65) Fujii, H.; Zhang, X.; Tomita, T.; Ikeda-Saito, M.; Yoshida, T. J. Am. Chem. Soc. 2001, 123, 6475-6484.

ARTICLES

Acknowledgment. This work was supported by the National Institutes of Health, GM62830 (GNL) and EB002064 (K.M.S.), and a Grant-in-Aid for Scientific Research (16570102) from the Ministry of Education and Sports, Science and Teaching, Japan (T.Y.). Supporting Information Available: Two tables (chemical shifts for substrate and for all assigned residues) and four figures (low-field NOESY for NmHO-DMDH-H2O, upfield TOCSY and low-field NOESY spectra for NmHO-2F4VDH-H2O, and low-field TOCSY spectra for NmHO-DMDH-H2O (11 total pages). This material is available free of charge via the Internet at http://pubs.acs.org. JA0578505

J. AM. CHEM. SOC.

9

VOL. 128, NO. 19, 2006 6399