Oxidation of heme to beta- and delta-biliverdin by Pseudomonas aeruginosa heme oxygenase as a consequence of an unusual seating of the heme

The origin of the unusual regioselectivity of heme oxygenation, i.e. the oxidation of heme to delta-biliverdin (70%) and beta-biliverdin (30%), that is exhibited by heme oxygenase from Pseudomonas aeruginosa (pa-HO) has been studied by (1)H NMR, (13)C NMR, and resonance Raman spectroscopies. Whereas resonance Raman indicates that the heme-iron ligation in pa-HO is homologous to that observed in previously studied alpha-hydroxylating heme oxygenases, the NMR spectroscopic studies suggest that the heme in this enzyme is seated in a manner that is distinct from that observed for all other alpha-hydroxylating heme oxygenase enzymes for which a structure is known. In pa-HO, the heme is rotated in-plane approximately 110 degrees, so the delta-meso-carbon of the major orientational isomer is located within the HO-fold in the place where the alpha-hydroxylating enzymes typically place the alpha-meso-carbon. The unusual heme seating displayed by pa-HO places the heme propionates so that these groups point in the direction of the solvent-exposed heme edge and appears to originate in large part from the absence of stabilizing interactions between the polypeptide and the heme propionates, which are typically found in alpha-hydroxylating heme oxygenase enzymes. These interactions typically involve Lys-16 and Tyr-112, in Neisseriae meningitidis HO, and Lys-16 and Tyr-134, in human and rat HO-1. The corresponding residues in pa-HO are Asn-19 and Phe-117, respectively. In agreement with this hypothesis, we found that the Asn-19 Lys/Phe-117 Tyr double mutant of pa-HO exists as a mixture of molecules exhibiting two distinct heme seatings; one seating is identical to that exhibited by wild-type pa-HO, whereas the alternative seating is very similar to that typical of alpha-hydroxylating heme oxygenase enzymes and is related to the wild-type seating by approximately 110 degrees in-plane rotation of the heme. Furthermore, each of these heme seatings in the pa-HO double mutant gives rise to a subset of two heme isomeric orientations that are related to each other by 180 degrees rotation about the alpha-gamma-meso-axis. The coexistence of these molecules in solution, in the proportions suggested by the corresponding area under the peaks in the (1)H NMR spectrum, explains the unusual regioselectivity of heme oxygenation observed with the double mutant, which we found produces alpha- (55%), delta- (35%), and beta-biliverdin (10%). Alpha-biliverdin is obtained by oxidation of the heme seated similar to that of alpha-hydroxylating enzymes, whereas beta- and delta-biliverdin are formed from the oxidation of heme seated as in wild-type pa-HO.     

Millisecond dynamics in glutaredoxin during catalytic turnover is dependent on substrate binding and absent in the resting states

Conformational dynamics is important for enzyme function. Which motions of enzymes determine catalytic efficiency and whether the same motions are important for all enzymes, however, are not well understood. Here we address conformational dynamics in glutaredoxin during catalytic turnover with a combination of NMR magnetization transfer, R(2) relaxation dispersion, and ligand titration experiments. Glutaredoxins catalyze a glutathione exchange reaction, forming a stable glutathinoylated enzyme intermediate. The equilibrium between the reduced state and the glutathionylated state was biochemically tuned to exchange on the millisecond time scale. The conformational changes of the protein backbone during catalysis were followed by (15)N nuclear spin relaxation dispersion experiments. A conformational transition that is well described by a two-state process with an exchange rate corresponding to the glutathione exchange rate was observed for 23 residues. Binding of reduced glutathione resulted in competitive inhibition of the reduced enzyme having kinetics similar to that of the reaction. This observation couples the motions observed during catalysis directly to substrate binding. Backbone motions on the time scale of catalytic turnover were not observed for the enzyme in the resting states, implying that alternative conformers do not accumulate to significant concentrations. These results infer that the turnover rate in glutaredoxin is governed by formation of a productive enzyme-substrate encounter complex, and that catalysis proceeds by an induced fit mechanism rather than by conformer selection driven by intrinsic conformational dynamics.     

Catalytic mechanism of heme oxygenase through EPR and ENDOR of cryoreduced oxy-heme oxygenase and its Asp 140 mutants

Heme oxygenase (HO) catalyzes the O(2)- and NADPH-cytochrome P450 reductase-dependent conversion of heme to biliverdin, Fe, and CO through a process in which the heme participates both as a prosthetic group and as a substrate. In the present study, we have generated a detailed reaction cycle for the first monooxygenation step of HO catalysis, conversion of the heme to alpha-meso-hydroxyheme. We employed EPR (using both (16)O(2) and (17)O(2)) and (1)H, (14)N ENDOR spectroscopies to characterize the intermediates generated by 77 K radiolytic cryoreduction and subsequent annealing of wild-type oxy-HO and D140A, F mutants. One-electron cryoreduction of oxy-HO yields a hydroperoxoferri-HO with g-tensor, g = [2.37, 2.187, 1.924]. Annealing of this species to 200 K is accompanied by spectroscopic changes that include the appearance of a new (1)H ENDOR signal, reflecting rearrangements in the active site. Kinetic measurements at 214 K reveal that the annealed hydroperoxoferri-HO species, denoted R, generates the ferri-alpha-meso-hydroxyheme product in a first-order reaction. Disruption of the H-bonding network within the distal pocket of HO by the alanine and phenylalanine mutations of residue D140 prevents product formation. The hydroperoxoferri-HO (D140A) instead undergoes heterolytic cleavage of the O-O bond, ultimately yielding an EPR-silent compound II-like species that does not form product. These results, which agree with earlier suggestions, establish that hydroperoxoferri-HO is indeed the reactive species, directly forming the alpha-meso-hydroxyheme product by attack of the distal OH of the hydroperoxo moiety at the heme alpha-carbon.     

Resonance Raman evidence for protein-induced out-of-plane distortion of the heme prosthetic group of mammalian lactoperoxidase

Resonance Raman spectra have been acquired for resting state mammalian lactoperoxidase, LPO(N), and its six-coordinate, low-spin (6CLS) cyanide complex, LPO(CN), as well as for various heme l containing fragments resulting from partial or complete proteolytic digestion. These proteolytic fragments provide a useful set of reference compounds for analysis of the LPO(N) and LPO(CN) enzymes, using various ligands to generate well-defined five-coordinate and six-coordinate high-spin (5CHS and 6CHS) species. In addition, these model compounds, which contain zero, one, or two covalently attached ester linkages to polypeptide chains, are quite useful for determining the extent to which the presence of the ester linkages at the heme periphery affects the characteristic heme resonance Raman marker bands. The spectral results not only provide strong evidence for the formulation of the resting state enzyme as a 6CHS species, but also confirm the previously documented anomalous intensities of several low-frequency resonance Raman bands, which are most reasonably interpreted to arise from a protein-induced out-of-plane distortion of the heme l macrocycle mediated by the covalent ester linkages to the associated polypeptide residues of the intact protein.     

Kinetic isotope effects on the rate-limiting step of heme oxygenase catalysis indicate concerted proton transfer/heme hydroxylation

Heme oxygenase (HO) catalyzes the O2 and NADPH/cytochrome P450 reductase-dependent conversion of heme to biliverdin, free iron ion, and CO through a process in which the heme participates as both dioxygen-activating prosthetic group and substrate. We earlier confirmed that the first step of HO catalysis is a monooxygenation in which the addition of one electron and two protons to the HO oxy-ferroheme produces ferric-alpha-meso-hydroxyheme (h). Cryoreduction/EPR and ENDOR measurements further showed that hydroperoxo-ferri-HO converts directly to h in a single kinetic step without formation of a Compound I. We here report details of that rate-limiting step. One-electron 77 K cryoreduction of human oxy-HO and annealing at 200 K generates a structurally relaxed hydroperoxo-ferri-HO species, denoted R. We here report the cryoreduction/annealing experiments that directly measure solvent and secondary kinetic isotope effects (KIEs) of the rate-limiting R --> h conversion, using enzyme prepared with meso-deuterated heme and in H2O/D2O buffers to measure the solvent KIE (solv-KIE), and the secondary KIE (sec-KIE) associated with the conversion. This approach is unique in that KIEs measured by monitoring the rate-limiting step are not susceptible to masking by KIEs of other processes, and these results represent the first direct measurement of the KIEs of product formation by a kinetically competent reaction intermediate in any dioxygen-activating heme enzyme.The observation of both solv-KIE(298) = 1.8 and sec-KIE(298) = 0.8 (inverse) indicates that the rate-limiting step for formation of h by HO is a concerted process: proton transfer to the hydroperoxo-ferri-heme through the distal-pocket H-bond network, likely from a carboxyl group acting as a general acid catalyst, occurring in synchrony with bond formation between the terminal hydroperoxo-oxygen atom and the alpha-meso carbon to form a tetrahedral hydroxylated-heme intermediate. Subsequent rearrangement and loss of H2O then generates h.     

Ultrafast heme dynamics in ferrous versus ferric cytochrome c studied by time-resolved resonance Raman and transient absorption spectroscopy

Cytochrome c (Cyt c) is a heme protein involved in electron transfer and also in apoptosis. Its heme iron is bisaxially ligated to histidine and methionine side chains and both ferric and ferrous redox states are physiologically relevant, as well as a ligand exchange between internal residue and external diatomic molecule. The photodissociation of internal axial ligand was observed for several ferrous heme proteins including Cyt c, but no time-resolved studies have been reported on ferric Cyt c. To investigate how the oxidation state of the heme influences the primary photoprocesses, we performed a comprehensive comparative study on horse heart Cyt c by subpicosecond time-resolved resonance Raman and femtosecond transient absorption spectroscopy. We found that in ferric Cyt c, in contrast to ferrous Cyt c, the photodissociation of an internal ligand does not take place, and relaxation dynamics is dominated by vibrational cooling in the ground electronic state of the heme. The intermolecular vibrational energy transfer was found to proceed in a single phase with a temperature decay of approximately 7 ps in both ferric and ferrous Cyt c. For ferrous Cyt c, the instantaneous photodissociation of the methionine side chain from the heme iron is the dominant event, and its rebinding proceeds in two phases, with time constants of approximately 5 and approximately 16 ps. A mechanism of this process is discussed, and the difference in photoinduced coordination behavior between ferric and ferrous Cyt c is explained by an involvement of the excited electronic state coupled with conformational relaxation of the heme.     

The hydroxide complex of Pseudomonas aeruginosa heme oxygenase as a model of the low-spin iron(III) hydroperoxide intermediate in heme catabolism: 13C NMR spectroscopic studies suggest the active participation of the heme in macrocycle hydroxylation

13C NMR spectroscopic studies have been conducted with the hydroxide complex of Pseudomonas aeruginosa heme oxygenase (Fe(III)-OH), where OH(-) has been used as a model of the OOH(-) ligand to gain insights regarding the elusive ferric hydroperoxide (Fe(III)-OOH) intermediate in heme catabolism at ambient temperatures. Analysis of the heme core carbon resonances revealed that the coordination of hydroxide in the distal site of the enzyme results in the formation of at least three populations of Fe(III)-OH complexes with distinct electronic configurations and nonplanar ring distortions that are in slow exchange relative to the NMR time scale. The most abundant population exhibits a spin crossover between S = (1)/(2) and S = (3)/(2) spin states, and the two less abundant populations exhibit pure, S = (3)/(2) and S = (1)/(2), (d(xy)())(1) electronic configurations. We propose that the highly organized network of water molecules in the distal pocket of heme oxygenase, by virtue of donating a hydrogen bond to the coordinated hydroxide ligand, lowers its ligand field strength, thereby increasing the field strength of the porphyrin (equatorial) ligand, which results in nonplanar deformations of the macrocycle. This tendency to deform from planarity, which is imparted by the ligand field strength of the coordinated OH(-), is likely reinforced by the flexibility of the distal pocket in HO. These findings suggest that if the ligand field strength of the coordinated OOH(-) in heme oxygenase is modulated in a similar manner, the resultant large spin density at the meso carbons and nonplanar deformations of the pophyrin ring prime the macrocycle to actively participate in its own hydroxylation.     

Azide-inhibited bacterial heme oxygenases exhibit an S = 3/2 (dxz,dyz)3(dxy)1(dz2)1 spin state: mechanistic implications for heme oxidation

The azide complexes of heme oxygenase from Pseudomonas aeruginosa (pa-HO) and Neisseriae meningitidis (nm-HO) have been studied with the aid of (1)H and (13)C NMR spectroscopy. These complexes have been shown to exist as an equilibrium mixture of two populations, one exhibiting an S = (1)/(2), (d(xy))(2)(d(xz), d(yz))(3) electron configuration and planar heme and a second with a novel S = (3)/(2), (d(xz), d(yz))(3)(d(xy))(1)(d(z)(2))(1) spin state and nonplanar heme. At physiologically relevant temperatures, the equilibrium shifts in the direction of the population exhibiting the latter electron configuration and nonplanar heme, whereas at temperatures approaching the freezing point of water, the equilibrium shifts in the direction of the population with the former electronic structure and planar heme. These findings indicate that the microenvironment of the distal pocket in heme oxygenase is unique among heme-containing proteins in that it lowers the sigma-donating (field strength) ability of the distal ligand and, therefore, promotes the attainment of heme electronic structures thus far only observed in heme oxygenase. When the field strength of the distal ligand is slightly lower than that of azide, such as OH(-) (J. Am. Chem. Soc. 2003, 125, 11842), the corresponding complex exists as a mixture of populations with nonplanar hemes and electronic structures that place significant spin density at the meso positions. The ease with which these unusual heme electronic structures are attained by heme oxygenase is likely related to activation of meso carbon reactivity which, in turn, facilitates hydroxylation of a meso carbon by the obligatory ferric hydroperoxide intermediate.     

Ultraviolet resonance Raman evidence for utilization of the heme 6-propionate hydrogen-bond network in signal transmission from heme to protein in Ec DOS protein

The direct oxygen sensor protein from Escherichia coli (Ec DOS) is a heme-based signal transducer protein responsible for phosphodiesterase (PDE) activity. Binding of either O2 or CO molecule to a reduced heme enhances the PDE activity toward 3',5'-cyclic diguanylic acid. We report ultraviolet resonance Raman (UVRR) spectroscopic investigations of the reduced, O2- and CO-bound forms of heme-bound PAS domain of Ec DOS. The UVRR results show that heme discriminates different ligands, resulting in altered conformations in the protein moiety. Specifically, the environment around Trp53 that contacts the 2-vinyl group of heme, is changed to a more hydrophobic environment by O2 binding, whereas it is changed to a more hydrophilic environment by CO-binding. In addition, the PDE activity of the O2- and CO-bound forms for the Trp53Phe mutant is significantly decreased compared with that of the wild type (WT), demonstrating the importance of Trp53 for the catalytic reaction. On the other hand, the binding of O2 or CO to the heme produces drastic changes in the Tyr126 of Ibeta-strand at the surface of the sensor domain. Furthermore, we found that Asn84 forms a hydrogen bond with Tyr126 either in the O2- or CO-bound forms but not in the reduced form. Finally, the PDE activities of the ligand-bound forms for Asn84Val and Tyr126Phe mutants are significantly reduced as compared with that of WT, suggesting the importance of the hydrogen-bonding network from heme 6-propionate to Tyr126 through Asn84 in signal transmission.     

Theory favors a stepwise mechanism of porphyrin degradation by a ferric hydroperoxide model of the active species of heme oxygenase

The report uses density functional theory to address the mechanism of heme degradation by the enzyme heme oxygenase (HO) using a model ferric hydroperoxide complex. HO is known to trap heme molecules and degrade them to maintain iron homeostasis in the biosystem. The degradation is initiated by complexation of the heme, then formation of the iron-hydroperoxo species, which subsequently oxidizes the meso position of the porphyrin by hydroxylation, thereby enabling eventually the cleavage of the porphyrin ring. Kinetic isotope effect studies indicate that the mechanism is assisted by general acid catalysis, via a chain of water molecules, and that all the events occur in concert. However, previous theoretical treatments indicated that the concerted mechanism has a high barrier, much higher than an alternative mechanism that is initiated by O-O bond homolysis of iron-hydroperoxide. The present contribution studies the stepwise and concerted acid-catalyzed mechanisms using H(3)O(+)(H(2)O)(n)(), n = 0-2. The effect of the acid strength is tested using the H(4)N(+)(H(2)O)(2) cluster and a fully protonated ferric hydroperoxide. All the calculations show that a stepwise mechanism that involves proton relay and O-O homolysis, in the rate-determining step, has a much lower barrier (>10 kcal/mol) than the corresponding fully concerted mechanism. The best fit of the calculated solvent kinetic isotope effect, to the experimental data, is obtained for the H(3)O(+)(H(2)O)(2) cluster. The calculated alpha-deuterium secondary kinetic isotope effect is inverse (0.95-0.98), but much less so than the experimental value (0.7). Possible reasons for this quantitative difference are discussed. Some probes are suggested that may enable experiment to distinguish the stepwise from the concerted mechanism.     

Freezing a single distal motion in dihydrofolate reductase

Constraining a single motion between distal residues separated by approximately 28 A in hybrid quantum/classical molecular dynamics simulations is found to increase the free energy barrier for hydride transfer in dihydrofolate reductase by approximately 3 kcal/mol. Our analysis indicates that a single distal constraint alters equilibrium motions throughout the enzyme on a wide range of time scales. This alteration of the conformational sampling of the entire system is sufficient to significantly increase the free energy barrier and decrease the rate of hydride transfer. Despite the changes in conformational sampling introduced by the constraint, the system assumes a similar transition state conformation with a donor-acceptor distance of approximately 2.72 A to enable the hydride transfer reaction. The modified thermal sampling leads to a substantial increase in the average donor-acceptor distance for the reactant state, however, thereby decreasing the probability of sampling the transition state conformations with the shorter distances required for hydride transfer. These simulations indicate that fast thermal fluctuations of the enzyme, substrate, and cofactor lead to conformational sampling of configurations that facilitate hydride transfer. The fast thermal motions are in equilibrium as the reaction progresses along the collective reaction coordinate, and the overall average equilibrium conformational changes occur on the slower time scale measured experimentally. Recent single molecule experiments suggest that at least some of these thermally averaged equilibrium conformational changes occur on the millisecond time scale of the hydride transfer reaction. Thus, introducing a constraint that modifies the conformational sampling of an enzyme could significantly impact its catalytic activity.     

Design of a five-coordinate heme protein maquette: a spectroscopic model of deoxymyoglobin

The substitution of 1-methyl-l-histidine for the histidine heme ligands in a de novo designed four-alpha-helix bundle scaffold results in conversion of a six-coordinate cytochrome maquette into a self-assembled five-coordinate mono-(1-methyl-histidine)-ligated heme as an initial maquette for the dioxygen carrier protein myoglobin. UV-vis, magnetic circular dichroism, and resonance Raman spectroscopies demonstrate the presence of five-coordinate mono-(1-methyl-histidine) ligated ferrous heme spectroscopically similar to deoxymyoglobin. Thermodynamic analysis of the ferric and ferrous heme dissociation constants indicates greater destabilization of the ferric state than the ferrous state. The ferrous heme protein reacts with carbon monoxide to form a (1-methyl-histidine)-Fe(II)(heme)-CO complex; however, reaction with dioxygen leads to autoxidation and ferric heme dissociation. These results indicate that negative protein design can be used to generate a five-coordinate heme within a maquette scaffold.     

Protein-based radicals in the catalase-peroxidase of synechocystis PCC6803: a multifrequency EPR investigation of wild-type and variants on the environment of the heme active site

Catalase-peroxidases are bifunctional heme enzymes with a high structural homology to peroxidases from prokaryotic origin and a catalatic activity comparable to monofunctional catalases. These unique features of catalase-peroxidases make them good systems to study and understand the role of alternative electron pathways both in catalases and peroxidases. In particular, it is of interest to study the poorly understood role of tyrosyl and tryptophanyl radicals as alternative cofactors in the catalytic cycle of catalases and peroxidases. In this work, we have used a powerful combination of multifrequency EPR spectroscopy, isotopic labeling of tryptophan and tyrosine residues, and site-directed mutagenesis to unequivocally identify the reactive intermediates formed by the wild-type Synechocystis PCC6803 catalase-peroxidase. Selected variants of the heme distal and proximal sides of the Synechocystis enzyme were investigated. Variants on the aromatic residues of the short stretch located relatively close to the heme and spanning the distal and proximal sides were also investigated. In the wild-type enzyme, the EPR signal of the catalases and peroxidases (typical) Compound I intermediate [Fe(IV)=O por.+] was observed. Two protein-based radical intermediates were also detected and identified as a Tyr. and a Trp. . The site of Trp. is proposed to be Trp 106, a residue belonging to the conserved short stretch in catalase-peroxidases and located at a 7-8 A distance to the heme propionate groups. An extensive hydrogen-bonding network on the heme distal side, involving Trp122, His123, Arg119, seven structural waters, the heme 6-propionate group, and Trp106, is proposed to have a key role on the formation of the tryptophanyl radical. We used high-field EPR spectroscopy (95-285 GHz) to resolve the g-anisotropy of the protein-based radicals in Synechocystis catalase-peroxidase. The broad gx component of the HF EPR spectrum of the Tyr. in Synechocystis catalase-peroxidase was consistent with a distributed electropositive protein environment to the tyrosyl radical.     

High catalytic activities of artificial peroxidases based on supramolecular hydrogels that contain heme models

Composed of a supramolecular hydrogel and a heme model compound, a new type of artificial peroxidase shows high catalytic activity in organic media. The activity of this new type of artificial enzyme is significantly higher than that of the heme model compounds alone. Changes in the distal substituents above the coordinated-metal centers of the model compounds directly modulate catalytic activity. This supramolecular-hydrogel-based artificial enzyme is most active in toluene, reaching about 90% of the nascent activity of horseradish peroxidase. Moreover, this study confirms that the incorporation of the heme models into the nanofibers of gelators accounts for most of the enhancement of catalytic activity.     

Oxygen activation and electron transfer in flavocytochrome P450 BM3

In flavocytochrome P450 BM3, there is a conserved phenylalanine residue at position 393 (Phe393), close to Cys400, the thiolate ligand to the heme. Substitution of Phe393 by Ala, His, Tyr, and Trp has allowed us to modulate the reduction potential of the heme, while retaining the structural integrity of the enzyme's active site. Substrate binding triggers electron transfer in P450 BM3 by inducing a shift from a low- to high-spin ferric heme and a 140 mV increase in the heme reduction potential. Kinetic analysis of the mutants indicated that the spin-state shift alone accelerates the rate of heme reduction (the rate determining step for overall catalysis) by 200-fold and that the concomitant shift in reduction potential is only responsible for a modest 2-fold rate enhancement. The second step in the P450 catalytic cycle involves binding of dioxygen to the ferrous heme. The stabilities of the oxy-ferrous complexes in the mutant enzymes were also analyzed using stopped-flow kinetics. These were found to be surprisingly stable, decaying to superoxide and ferric heme at rates of 0.01-0.5 s(-)(1). The stability of the oxy-ferrous complexes was greater for mutants with higher reduction potentials, which had lower catalytic turnover rates but faster heme reduction rates. The catalytic rate-determining step of these enzymes can no longer be the initial heme reduction event but is likely to be either reduction of the stabilized oxy-ferrous complex, i.e., the second flavin to heme electron transfer or a subsequent protonation event. Modulating the reduction potential of P450 BM3 appears to tune the two steps in opposite directions; the potential of the wild-type enzyme appears to be optimized to maximize the overall rate of turnover. The dependence of the visible absorption spectrum of the oxy-ferrous complex on the heme reduction potential is also discussed.     

The ferrous verdoheme-heme oxygenase complex is six-coordinate and low-spin

A biosynthetic and enzymatic method was developed for the preparation of 13C-labeled verdoheme, which permits the 13C NMR spectroscopic characterization of this elusive intermediate in the heme oxidation path catalyzed by the enzyme heme oxygenase. The 13C NMR data indicate that the ferrous verdoheme complex of Neisseria meningitides heme oxygenase is hexacoordinate and diamagnetic, with a proximal histidine and likely a distal hydroxide as axial ligands. The coordination number and spin state of the ferrous verdoheme-heme oxygenase complex is in stark contrast to the pentacoordinate and paramagnetic nature of the heme-heme oxygenase complex and heme centers in general.     

A role for highly conserved carboxylate, aspartate-140, in oxygen activation and heme degradation by heme oxygenase-1.

Heme oxygenase (HO) catalyzes the oxygen-dependent degradation of heme to biliverdinIXalpha, CO, and free iron ion via three sequential monooxygenase reactions. Although the distinct active-site structure of HO from cytochrome P450 families suggests unique distal protein machinery to activate molecular oxygen, the mechanism and the key amino acid for the oxygen activation have not been clear. To investigate the functionality of highly conserved polar amino acids in the distal helix of HO-1, we have prepared alanine mutants: T135A, R136A, D140A, and S142A, and found drastic changes in the heme degradation reactions of D140A. In this paper, we report the first evidence that D140 is involved in the oxygen activation mechanism in HO-1. The heme complexes of HO mutants examined in this study fold and bind heme normally. The pK(a) values of the iron-bound water and autoxidation rates of the oxy-form are increased with R136A, D140A, and S142A mutations, but are not changed with T135A mutation. As the wild-type, T135A, R136A, and S142A degrade heme to verdohemeIXalpha with H(2)O(2) and to biliverdinIXalpha with the NADPH reductase system. On the other hand, D140A heme complex forms compound II with H(2)O(2), and no heme degradation occurs. For the NADPH reductase system, the oxy-form of D140A heme complex is accumulated in the reaction, and only 50% of heme is degraded. The stopped flow experiments suggest that D140A cannot activate iron-bound dioxygen and hydroperoxide properly. To investigate the carboxylate functionality of D140, we further replaced D140 with glutamic acid (D140E), phenylalanine (D140F), and asparagine (D140N). D140E degrades heme normally, but D140N shows reactivity similar to that of D140A. D140F loses heme degradation activity completely. All of these results indicate that the carboxylate at position 140 is essential to activate the iron-bound dioxygen and hydroperoxide. On the basis of the present findings, we propose an oxygen activation mechanism involving the hydrogen-bonding network through the bridging water and D140 side chain.     

How do azoles inhibit cytochrome P450 enzymes? A density functional study

To examine how azole inhibitors interact with the heme active site of the cytochrome P450 enzymes, we have performed a series of density functional theory studies on azole binding. These are the first density functional studies on azole interactions with a heme center and give fundamental insight into how azoles inhibit the catalytic function of P450 enzymes. Since azoles come in many varieties, we tested three typical azole motifs representing a broad range of azole and azole-type inhibitors: methylimidazolate, methyltriazolate, and pyridine. These structural motifs represent typical azoles, such as econazole, fluconazole, and metyrapone. The calculations show that azole binding is a stepwise mechanism whereby first the water molecule from the resting state of P450 is released from the sixth binding site of the heme to create a pentacoordinated active site followed by coordination of the azole nitrogen to the heme iron. This process leads to the breaking of a hydrogen bond between the resting state water molecule and the approaching inhibitor molecule. Although, formally, the water molecule is released in the first step of the reaction mechanism and a pentacoordinated heme is created, this does not lead to an observed spin state crossing. Thus, we show that release of a water molecule from the resting state of P450 enzymes to create a pentacoordinated heme will lead to a doublet to quartet spin state crossing at an Fe-OH(2) distance of approximately 3.0 A, while the azole substitution process takes place at shorter distances. Azoles bind heme with significantly stronger binding energies than a water molecule, so that these inhibitors block the catalytic cycle of the enzyme and prevent oxygen binding and the catalysis of substrate oxidation. Perturbations within the active site (e.g., a polarized environment) have little effect on the relative energies of azole binding. Studies with an extra hydrogen-bonded ethanol molecule in the model, mimicking the active site of the CYP121 P450, show that the resting state and azole binding structures are close in energy, which may lead to chemical equilibrium between the two structures, as indeed observed with recent protein structural studies that have demonstrated two distinct azole binding mechanisms to P450 heme.     

How heme metabolism occurs in heme oxygenase: computational study of oxygen-donation ability of the oxo and hydroperoxo species

We report a density functional theory study on the heme metabolism in heme oxygenase using iron-hydroperoxo and -oxo models. The activation energies for heme oxidation at the alpha-carbon by the iron-hydroperoxo and -oxo species are calculated to be 42.9 and 39.9 kcal/mol, respectively. These high activation barriers lead us to reconsider the catalytic mechanism of heme oxygenase