Spectroscopic and Crystallographic Evidence for the Role of a Water-Containing H-Bond Network in Oxidase Activity of an Engineered Myoglobin.

Heme-copper oxidases (HCOs) catalyze efficient reduction of oxygen to water in biological respiration. Despite progress in studying native enzymes and their models, the roles of non-covalent interactions in promoting this activity are still not well understood. Here we report EPR spectroscopic studies of cryoreduced oxy-F33Y-CuBMb, a functional model of HCOs engineered in myoglobin (Mb). We find that cryoreduction at 77 K of the O2-bound form, trapped in the conformation of the parent oxyferrous form, displays a ferric-hydroperoxo EPR signal, in contrast to the cryoreduced oxy-wild-type (WT) Mb, which is unable to deliver a proton and shows a signal from the peroxo-ferric state. Crystallography of oxy-F33Y-CuBMb reveals an extensive H-bond network involving H2O molecules, which is absent from oxy-WTMb. This H-bonding proton-delivery network is the key structural feature that transforms the reversible oxygen-binding protein, WTMb, into F33Y-CuBMb, an oxygen-activating enzyme that reduces O2 to H2O. These results provide direct evidence of the importance of H-bond networks involving H2O in conferring enzymatic activity to a designed protein. Incorporating such extended H-bond networks in designing other metalloenzymes may allow us to confer and fine-tune their enzymatic activities.

Activation of O2 is important in biology because O2, despite its high abundance in the Earth’s atmosphere and high oxidizing power, is kinetically inert. A prominent example of O2 activation is aerobic respiration, in which the oxidizing potential of O2 is used to generate a proton gradient, driving synthesis of adenosine triphosphate, the energy source for most biological processes. This reaction is efficiently catalyzed by terminal oxidases, e.g., heme-copper oxidases (HCOs) containing a heme and CuB center[1] and cytochrome (cyt) bd-oxygen reductases with only hemes at the active site.[2] While these enzymes and their model systems have been extensively studied and mechanisms have been proposed,[1−3] an understanding of structural features responsible for O2 activation is still being sought. For example, while a number reports have elucidated H-bond networks in the channels leading to the CuB site in HCOs,[1,3a,3b] experimental support for structural features, e.g., a putative H-bond network involving H2O within the CuB site, proposed by computational modeling,[4] remains elusive.

To elucidate structural features responsible for efficient O2 activation, biochemical and biophysical studies of the oxidases and related enzymes have identified potentially important amino acid residues in the reaction, which has been confirmed by loss of activity in site-directed mutagenesis (SDM) studies.[1,3a,3b,5,6] However, these mutations often cause structural perturbations in addition to those directly being investigated, which may cause loss of function. Moreover, SDM studies reveal only necessary, but not minimally sufficient structural features for function. To overcome this limitation, synthetic models of native enzymes have been prepared, and their studies have offered deeper insights into the structural features responsible for oxidase activity.[3g,3h,7]

As a complementary approach to the studies of native enzymes and their synthetic models, use of small designed metalloproteins and -peptides has shown remarkable success.[8] For example, we have employed a “biosynthetic modeling” approach, based on a small natural protein scaffold, to gain fundamental understanding of sufficient features conferring oxidase activity, i.e., 4e– reduction of O2 to H2O.[8a,9] In contrast to native enzymes, such as HCOs, which are membrane proteins that are not easy to purify to homogeneity with high yields and contain other accessory metal sites that can complicate spectroscopic studies, biosynthetic models are simpler to synthesize and contain only the metal-binding site central to the enzymatic activity. While most HCO models use organic molecules as ligands, our biosynthetic approach uses small, stable, and well-characterized proteins, such as sperm whale myoglobin (Mb), as the scaffold. Thanks to recent advances in biology, biosynthetic models of complex metal-binding sites can be prepared more easily and with higher yields than using organic molecule ligands. Furthermore, the structurally rigid scaffold protein allows defining long-range non-covalent interactions and probing their roles more precisely. Using this approach, we previously reported the design of an HCO mimic by introducing two His and one Tyr residues in the active site of wild-type (WT) Mb (named F33Y-CuBMb) and found that, unlike WTMb, which only binds O2 reversibly, this protein achieves O2 reduction to H2O with >500 turnovers.[9d] Thus, this minimal functional mimic, much simpler than native oxidases while still functioning under physiological conditions, provides a platform for elucidating structural features responsible for O2 reduction that bridges the gap between native enzymes and small molecular mimics of native enzymes.

We hypothesize that F33Y-CuBMb operates by heterolytic O–O bond cleavage (Scheme ), analogous to the accepted mechanism for native oxidases and monooxygenases, but the detailed structural features and interactions that promote this reaction in F33Y-CuBMb are still unclear. Here we use cryogenic EPR spectroscopy and high-resolution X-ray crystallography to show that an extensive H-bond network involving H2O molecules exists in the active site of F33Y-CuBMb and plays an important role in facilitating proton delivery to the oxygen and activating it for reduction, thus, imparting oxidase activity.

Scheme 1

Proposed Mechanism of O2 Reduction by F33Y-CuBMb in Comparison with WTMb

EPR has been widely used to probe the environment of O2-activating iron enzymes.[10] Although oxy-heme is an electronically closed-shell system with no EPR signal, EPR spectroscopy has been used to offer insight into the properties of oxy-heme, e.g., its interactions with nearby residues and its intermediates upon reduction, by a method called cryoreduction.[11] Briefly, the oxy-heme protein, trapped at 77 K, is exposed to γ-radiation, releasing “free” electrons into the sample matrix to efficiently reduce EPR-silent oxy-heme into EPR-active peroxo-heme. Because cryogenic trapping suppresses thermal relaxation and equilibration of the protein and its active site, EPR spectra of the resulting sample, maintained at 77 K, provide valuable structural information on the oxy-heme precursor (Scheme ).[10a,10f−10l,11] Furthermore, the peroxo-heme products are intermediates along the pathway of O2 reduction, and subsequent decay, achieved by stepwise warming of the samples, followed by rapid recooling to 77 K after each step to trap intermediates-a procedure called annealing, provides further insight into the reaction (Scheme ).[10a,11,12] As a result, EPR spectroscopy in combination with cryoreduction and annealing has been used over the past two decades to probe the structure and interactions of O2 in many heme proteins, including myoglobins and monooxygenases.[10a,10f−10q,11,12] Because WTMb, the scaffold for F33Y-CuBMb, was among the first and most studied proteins by the above method, the abundance of results facilitates understanding of the latter by direct comparison. In this study, F33Y-CuBMb contained heme but no metal ion in the CuB site, because a previous report showed that this protein exhibits high oxidase activity without a metal ion at CuB,[9d] mimicking cyt bd-oxygen reductases, which lack this site.[2]

Scheme 2

Formation and Trapping of (Hydro)peroxo Intermediates by Cryoreduction

Figure shows EPR spectra of oxy-WTMb and oxy-F33Y-CuBMb cryoreduced at 77 K, obtained under the same conditions. The WTMb spectrum shows the signal from a single species with g = [2.22, 2.13, 1.96] (see Figure S1 for spectra at other temperatures) whose EPR parameters (Table ) identify it as a ferric-peroxo state, with the peroxo ligand H-bonded to the distal His-64.[10f−10i,11] In contrast, the cryoreduced oxy-F33Y-CuBMb exhibits two rhombic EPR signals (Figures , S1, and S2). The first, with g = [2.24, 2.13, 1.95], is similar to that of the ferric-peroxo species in WTMb. The second, more rhombic EPR signal, has g = [2.29, 2.17, 1.96], in the range characteristic of ferric-hydroperoxo intermediates.[10a,10h,10k−10m] This assignment is supported by 1H ENDOR measurements (Figure S6), where the gmax = 2.24 species has been removed by annealing, showing a maximum proton coupling of ∼9.5 MHz. Ferric-peroxo intermediates show strongly coupled exchangeable proton signals with Amax ≈ 14–18 MHz,[10a,10i,10o] while hydroperoxo- exhibit a smaller coupling, Amax ≤ 13 MHz.[10j,10k,10q] The absence of signal with Amax > 13 MHz confirms that the more rhombic species observed in F33Y-CuBMb is hydroperoxo, not peroxo.

Figure 1

EPR spectra of oxy-WTMb and oxy-F33Y CuBMb after radiolytic reduction with a ∼3 Mrad dose of γ radiation from 60Co, and after subsequent stepwise annealing for 1 min at indicated temperatures (sharp signal marked by asterisk is due to radiolytically generated H-atoms at 77 K in quartz EPR tube).

Table 1

g-Values of Cryoreduced Oxy-ferrous Myoglobins

protein		g1	g2	g3
WTMb	peroxo	2.22	2.12	1.965
	hydroperoxoa	2.32	2.18	1.943
F33y CuBMb	peroxo	2.24	2.13	1.964
	hydroperoxob	2.29	2.16	1.950
	hydroperoxoa	2.34	2.19	1.937

Hydroperoxo product formed upon annealing.

Primary hydroperoxo product formed at 77 K.

To gain further insight into O2 reduction by F33Y-CuBMb, we monitored the decay of the above cryoreduced species by annealing. As previously reported, annealing the cryoreduced oxy-WTMb at 190 K converts the peroxo signal with gmax = 2.22 to the more rhombic EPR signal of a ferric-hydroperoxo intermediate, with g = [2.32, 2.18, 1.93] (Figures and S1), through protonation of the peroxo moiety.[13] Upon annealing of cryoreduced oxy-F33Y-CuBMb at 173–175 K for at least 1 min, a signal with g = [2.34, 2.19, 1.94] begins to appear and continues to grow in at 190 K (Figures , S1, and S2). These EPR parameters are very similar to those of ferric-hydroperoxo-WTMb. Appearance of this signal during annealing is interpreted to arise from both protonation of the ferric-peroxo centers with gmax = 2.24 to form ferric-hydroperoxo, and relaxation/repositioning of the hydroperoxo ligand of the gmax = 2.29 species of the heme iron (Scheme ).

Based on the known reactivity of the (hydro)peroxo-hemes, these intermediates are expected to thermally decompose to compound II (CpdII), a ferryl species that is EPR silent. Consistent with this expectation, further annealing both samples at 210 and 220 K (Figures and S1) results in obvious decay of the (hydro)peroxo signals to an EPR-silent state. It has been shown that further cryoreduction of the EPR-silent CpdII (Fe(IV)=O2–) at 77 K generates a low-spin rhombic EPR signal from the Fe(III)–O2– or Fe(III)–OH– species, depending on the protonation state of ferryl precursor. Subsequent cryoreduction of the EPR-silent species in WTMb resulted in a signal with g = [2.43, −, 1.93] (Figure S3), characteristic of Fe(III)–O2– (ferric oxo) species observed in cryoreduced CpdII (g = [2.43, 2.12, 1.93]),[12d] confirming the EPR-silent state arising after annealing of (hydro)peroxo-WTMb as the Fe(IV)=O2– CpdII. A minor signal with g = [2.52, −, 1.90] (Figure S3) was also observed, previously assigned as protonated ferric-oxo, i.e., ferric-hydroxo (Fe(III)–OH–).[12d] In the case of F33Y-CuBMb, further cryoreduction of the EPR-silent annealed species yielded a dominant EPR signal with g = [2.52, 2.15, 1.90] (Figure S3), likewise assigned as ferric-hydroxo. This result supports the hypothesis that the hydroperoxo-F33Y-CuBMb converts to CpdII. This assignment is confirmed by comparison with the same signals observed in cryoreduced samples of CpdII, generated in both proteins by reaction with H2O2 at 21–25 °C (Figure S4).

Observation of two distinct signals in cryoreduced oxy-F33Y-CuBMb at 77 K suggests the presence of two structural states of the oxy-ferrous precursor. The peroxo signal with gmax = 2.24 is nearly identical to the signal observed in WTMb and well within the range of those of other oxygen-binding proteins such as hemoglobin,[10h,10i] suggesting that the oxy-ferrous population leading to this state has similar interactions and behaves in a similar manner. On the other hand, the observation of a hydroperoxo species with gmax = 2.29 at 77 K requires a configuration in which the oxygen can be easily protonated at 77 K (Scheme ). The observation of protonated Fe(III)–OH– in F33Y-CuBMb upon subsequent cryoreduction of annealed species provides further support for the facile protonation of intermediates. Previous cryoreduction studies of a variety of oxy-heme proteins[10a,10j−10o,12a,12c] showed that proton transfer to a cryoreduced oxy-heme at 77 K requires the presence of an extended H-bonded proton delivery network in the oxy-ferrous protein that includes protic residues and at least one ordered H2O molecule H-bonded to the terminal O-atom of the O2 ligand. Based on these observations, we hypothesized that the F33Y-CuBMb may contain a similar H-bond network for faster proton delivery to the oxy-heme than in WTMb. To test this, we obtained the crystal structure of oxy-F33Y-CuBMb.

The oxy-F33Y-CuBMb crystal was prepared by soaking ferric-Mb crystals in dithionite, followed by exposure to O2,[14] and the oxy state was confirmed by single-crystal UV/vis spectroscopy before diffraction (Figure S5). The structure, refined at 1.27 Å resolution, is compared with the structure of oxy-WTMb in Figure .[15] An omit-map showing the electron density of the oxygen ligand is shown in Figure S6. Due to reduced oxygen affinity of F33Y-CuBMb, refinement of this structure yielded an incomplete occupancy for the oxygen ligand, but it nevertheless supports the above interpretation of the EPR data. As suggested by the EPR results, an extended H-bond network involving two H2O molecules is observed, with W1 stabilized by H-bonds with the engineered His43 and Tyr33 residues, and W2 stabilized by H-bonds to W1 and the engineered His29 (Figure ). Together, this H-bond network links the designed residues and the distal O-atom of bound O–O. As a result of these changes in the active-site pocket, the O–O in F33Y-CuBMb is rotated from that in WTMb by ∼13° about the Fe–O bond, apparently to maximize its ability to accept H-bonds from His64 and W2. Such an interaction of the oxygen with the W2 is expected to enable the facile protonation, by W2, of the peroxo intermediate formed on cryoreduction at 77 K, yielding the hydroperoxo observed by the cryogenic EPR.

Figure 2

Crystal structure of oxy-F33Y-CuBMb determined at 1.27 Å resolution (PDB: 5HAV), compared with that of oxy-WTMb (1A6M).[16]

Given the spectroscopic and crystallographic evidence for an H-bond network that interacts with the O2 in oxy-F33Y-CuBMb, and its absence in oxy-WTMb, we conclude that this H-bond network is the important structural feature that transforms WTMb, a protein that can only bind O2 reversibly, into F33Y-CuBMb, a protein that can activate O2 and reduce it to H2O. We propose that this activation is achieved by polarization of the O–O bond by the H-bonded H2O, allowing facile reduction and protonation of the oxygen and intermediates. While H-bond networks involving water have been observed in native heme enzymes such as cyt P450[10l,10m,17] and a similar network has been proposed to occur in HCOs-based computational and isotope studies,[4] design of metalloenzymes incorporating this important structural feature has been difficult to perform and confirm. The combination of cryoreduction EPR spectroscopy and high-resolution crystallography presented in this work provides a clear example of the importance of H-bond networks involving water in conferring activity to designed metalloenzymes. We believe this structural feature will be critical to enhancing the future success of high-activity metalloenzyme design.