Dioxygen activation at non-heme diiron centers: oxidation of a proximal residue in the I100W variant of toluene/o-xylene monooxygenase hydroxylase

Computational Examination on the Active Site Structure of a (Peroxo)diiron(III) Intermediate in the Amine Oxygenase AurF

In this work, we report the first computational investigation on the structure and properties of the (peroxo)diiron(III) intermediate of the AurF enzyme. Our calculations predict that, in the oxidized state of the AurF enzyme, the peroxo ligand is depicted in a μ-1,1-coordination mode with a protonated bridging ligand and is not in a μ-η(2):η(2) or μ-1,2 mode. Computed spectral data for the μ-1,1-coordination mode correlate well with experimental observations and unravel the potential of the energetics-spectroscopic approach adapted here.

Crystallographic analysis of active site contributions to regiospecificity in the diiron enzyme toluene 4-monooxygenase

Crystal structures of toluene 4-monooxygenase hydroxylase in complex with reaction products and effector protein reveal active site interactions leading to regiospecificity. Complexes with phenolic products yield an asymmetric μ-phenoxo-bridged diiron center and a shift of diiron ligand E231 into a hydrogen bonding position with conserved T201. In contrast, complexes with inhibitors p-NH(2)-benzoate and p-Br-benzoate showed a μ-1,1 coordination of carboxylate oxygen between the iron atoms and only a partial shift in the position of E231. Among active site residues, F176 trapped the aromatic ring of products against a surface of the active site cavity formed by G103, E104 and A107, while F196 positioned the aromatic ring against this surface via a π-stacking interaction. The proximity of G103 and F176 to the para substituent of the substrate aromatic ring and the structure of G103L T4moHD suggest how changes in regiospecificity arise from mutations at G103. Although effector protein binding produced significant shifts in the positions of residues along the outer portion of the active site (T201, N202, and Q228) and in some iron ligands (E231 and E197), surprisingly minor shifts (<1 Å) were produced in F176, F196, and other interior residues of the active site. Likewise, products bound to the diiron center in either the presence or absence of effector protein did not significantly shift the position of the interior residues, suggesting that positioning of the cognate substrates will not be strongly influenced by effector protein binding. Thus, changes in product distributions in the absence of the effector protein are proposed to arise from differences in rates of chemical steps of the reaction relative to motion of substrates within the active site channel of the uncomplexed, less efficient enzyme, while structural changes in diiron ligand geometry associated with cycling between diferrous and diferric states are discussed for their potential contribution to product release.

Tracking a defined route for O₂ migration in a dioxygen-activating diiron enzyme

For numerous enzymes reactive toward small gaseous compounds, growing evidence indicates that these substrates diffuse into active site pockets through defined pathways in the protein matrix. Toluene/o-xylene monooxygenase hydroxylase is a dioxygen-activating enzyme. Structural analysis suggests two possible pathways for dioxygen access through the α-subunit to the diiron center: a channel or a series of hydrophobic cavities. To distinguish which is utilized as the O(2) migration pathway, the dimensions of the cavities and the channel were independently varied by site-directed mutagenesis and confirmed by X-ray crystallography. The rate constants for dioxygen access to the diiron center were derived from the formation rates of a peroxodiiron(III) intermediate, generated upon treatment of the diiron(II) enzyme with O(2). This reaction depends on the concentration of dioxygen to the first order. Altering the dimensions of the cavities, but not the channel, changed the rate of dioxygen reactivity with the enzyme. These results strongly suggest that voids comprising the cavities in toluene/o-xylene monooxygenase hydroxylase are not artifacts of protein packing/folding, but rather programmed routes for dioxygen migration through the protein matrix. Because the cavities are not fully connected into the diiron active center in the enzyme resting state, conformational changes will be required to facilitate dioxygen access to the diiron center. We propose that such temporary opening and closing of the cavities may occur in all bacterial multicomponent monooxygenases to control O(2) consumption for efficient catalysis. Our findings suggest that other gas-utilizing enzymes may employ similar structural features to effect substrate passage through a protein matrix.

Structure and dynamics of a dizinc metalloprotein: effect of charge transfer and polarization

Structures and dynamics of a recently designed dizinc metalloprotein (DFsc) (J. Mol. Biol. 2003, 334, 1101) are studied by molecular dynamics simulation using a dynamically adapted polarized force field derived from fragment quantum calculation for protein in solvent. To properly describe the effect of charge transfer and polarization in the present approach, quantum chemistry calculation of the zinc-binding group is periodically performed (on-the-fly) to update the atomic charges of the zinc-binding group during the MD simulation. Comparison of the present result with those obtained from simulations under standard AMBER force field reveals that charge transfer and polarization are critical to maintaining the correct asymmetric metal coordination in the DFsc. Detailed analysis of the result also shows that dynamic fluctuation of the zinc-binding group facilitates solvent interaction with the zinc ions. In particular, the dynamic fluctuation of the zinc-zinc distance is shown to be an important feature of the catalytic function of the di-ion zinc-binding group. Our study demonstrates that the dynamically adapted polarization approach is computationally practical and can be used to study other metalloprotein systems.

Mechanistic studies of reactions of peroxodiiron(III) intermediates in T201 variants of toluene/o-xylene monooxygenase hydroxylase

Site-directed mutagenesis studies of a strictly conserved T201 residue in the active site of toluene/o-xylene monooxygenase hydroxylase (ToMOH) revealed that a single mutation can facilitate kinetic isolation of two distinctive peroxodiiron(III) species, designated T201(peroxo) and ToMOH(peroxo), during dioxygen activation. Previously, we characterized both oxygenated intermediates by UV-vis and Mössbauer spectroscopy, proposed structures from DFT and QM/MM computational studies, and elucidated chemical steps involved in dioxygen activation through the kinetic studies of T201(peroxo) formation. In this study, we investigate the kinetics of T201(peroxo) decay to explore the reaction mechanism of the oxygenated intermediates following O(2) activation. The decay rates of T201(peroxo) were monitored in the absence and presence of external (phenol) or internal (tryptophan residue in an I100W variant) substrates under pre-steady-state conditions. Three possible reaction models for the formation and decay of T201(peroxo) were evaluated, and the results demonstrate that this species is on the pathway of arene oxidation and appears to be in equilibrium with ToMOH(peroxo).

Insights into the different dioxygen activation pathways of methane and toluene monooxygenase hydroxylases

The methane and toluene monooxygenase hydroxylases (MMOH and TMOH, respectively) have almost identical active sites, yet the physical and chemical properties of their oxygenated intermediates, designated P*, H(peroxo), Q, and Q* in MMOH and ToMOH(peroxo) in a subclass of TMOH, ToMOH, are substantially different. We review and compare the structural differences in the vicinity of the active sites of these enzymes and discuss which changes could give rise to the different behavior of H(peroxo) and Q. In particular, analysis of multiple crystal structures reveals that T213 in MMOH and the analogous T201 in TMOH, located in the immediate vicinity of the active site, have different rotatory configurations. We study the rotational energy profiles of these threonine residues with the use of molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) computational methods and put forward a hypothesis according to which T213 and T201 play an important role in the formation of different types of peroxodiiron(III) species in MMOH and ToMOH. The hypothesis is indirectly supported by the QM/MM calculations of the peroxodiiron(III) models of ToMOH and the theoretically computed Mössbauer spectra. It also helps explain the formation of two distinct peroxodiiron(III) species in the T201S mutant of ToMOH. Additionally, a role for the ToMOD regulatory protein, which is essential for intermediate formation and protein functioning in the ToMO system, is advanced. We find that the low quadrupole splitting parameter in the Mössbauer spectrum observed for a ToMOH(peroxo) intermediate can be explained by protonation of the peroxo moiety, possibly stabilized by the T201 residue. Finally, similarities between the oxygen activation mechanisms of the monooxygenases and cytochrome P450 are discussed.

Multiple roles of component proteins in bacterial multicomponent monooxygenases: phenol hydroxylase and toluene/o-xylene monooxygenase from Pseudomonas sp. OX1

Phenol hydroxylase (PH) and toluene/o-xylene monooxygenase (ToMO) from Pseudomonas sp. OX1 require three or four protein components to activate dioxygen for the oxidation of aromatic substrates at a carboxylate-bridged diiron center. In this study, we investigated the influence of the hydroxylases, regulatory proteins, and electron-transfer components of these systems on substrate (phenol; NADH) consumption and product (catechol; H(2)O(2)) generation. Single-turnover experiments revealed that only complete systems containing all three or four protein components are capable of oxidizing phenol, a major substrate for both enzymes. Under ideal conditions, the hydroxylated product yield was ∼50% of the diiron centers for both systems, suggesting that these enzymes operate by half-sites reactivity mechanisms. Single-turnover studies indicated that the PH and ToMO electron-transfer components exert regulatory effects on substrate oxidation processes taking place at the hydroxylase actives sites, most likely through allostery. Steady state NADH consumption assays showed that the regulatory proteins facilitate the electron-transfer step in the hydrocarbon oxidation cycle in the absence of phenol. Under these conditions, electron consumption is coupled to H(2)O(2) formation in a hydroxylase-dependent manner. Mechanistic implications of these results are discussed.

Mössbauer, electron paramagnetic resonance, and density functional theory studies of synthetic S = 1/2 Fe(III)-O-Fe(IV)═O complexes. Superexchange-mediated spin transition at the Fe(IV)═O site

Previously we have characterized two high-valent complexes [LFe(IV)(μ-O)(2)Fe(III)L], 1, and [LFe(IV)(O)(μ-O)(OH) Fe(IV)L], 4. Addition of hydroxide or fluoride to 1 produces two new complexes, 1-OH and 1-F. Electron paramagnetic resonance (EPR) and Mössbauer studies show that both complexes have an S = 1/2 ground state which results from antiferromagnetic coupling of the spins of a high-spin (S(a) = 5/2) Fe(III) and a high-spin (S(b) = 2) Fe(IV) site. 1-OH can also be obtained by a 1-electron reduction of 4, which has been shown to have an Fe(IV)═O site. Radiolytic reduction of 4 at 77 K yields a Mössbauer spectrum identical to that observed for 1-OH, showing that the latter contains an Fe(IV)═O. Interestingly, the Fe(IV)═O moiety has S(b) = 1 in 4 and S(b) = 2 in 1-OH and 1-F. From the temperature dependence of the S = 1/2 signal we have determined the exchange coupling constant J (ℋ = JS(a)·S(b) convention) to be 90 ± 20 cm(-1) for both 1-OH and 1-F. Broken-symmetry density functional theory (DFT) calculations yield J = 135 cm(-1) for 1-OH and J = 104 cm(-1) for 1-F, in good agreement with the experiments. DFT analysis shows that the S(b) = 1 → S(b) = 2 transition of the Fe(IV)═O site upon reduction of the Fe(IV)-OH site to high-spin Fe(III) is driven primarily by the strong antiferromagnetic exchange in the (S(a) = 5/2, S(b) = 2) couple.

Four-electron oxidation of p-hydroxylaminobenzoate to p-nitrobenzoate by a peroxodiferric complex in AurF from Streptomyces thioluteus

The nonheme di-iron oxygenase, AurF, converts p-aminobenzoate (Ar-NH(2), where Ar = 4-carboxyphenyl) to p-nitrobenzoate (Ar-NO(2)) in the biosynthesis of the antibiotic, aureothin, by Streptomyces thioluteus. It has been reported that this net six-electron oxidation proceeds in three consecutive, two-electron steps, through p-hydroxylaminobenzoate (Ar-NHOH) and p-nitrosobenzoate (Ar-NO) intermediates, with each step requiring one equivalent of O(2) and two exogenous reducing equivalents. We recently demonstrated that a peroxodiiron(III/III) complex (peroxo- -AurF) formed by addition of O(2) to the diiron(II/II) enzyme ( -AurF) effects the initial oxidation of Ar-NH(2), generating a mu-(oxo)diiron(III/III) form of the enzyme (mu-oxo- -AurF) and (presumably) Ar-NHOH. Here we show that peroxo- -AurF also oxidizes Ar-NHOH. Unexpectedly, this reaction proceeds through to the Ar-NO(2) final product, a four-electron oxidation, and produces -AurF, with which O(2) can combine to regenerate peroxo- -AurF. Thus, conversion of Ar-NHOH to Ar-NO(2) requires only a single equivalent of O(2) and (starting from -AurF or peroxo- -AurF) is fully catalytic in the absence of exogenous reducing equivalents, by contrast to the published stoichiometry. This novel type of four-electron N-oxidation is likely also to occur in the reaction sequences of nitro-installing di-iron amine oxygenases in the biosyntheses of other natural products.

Characterization of a peroxodiiron(III) intermediate in the T201S variant of toluene/o-xylene monooxygenase hydroxylase from Pseudomonas sp. OX1

We report the observation of a novel intermediate in the reaction of a reduced toluene/o-xylene monooxygenase hydroxylase (ToMOH(red)) T201S variant, in the presence of a regulatory protein (ToMOD), with dioxygen. This species is the first oxygenated intermediate with an optical band in any toluene monooxygenase. The UV-vis and Mossbauer spectroscopic properties of the intermediate allow us to assign it as a peroxodiiron(III) species, T201S(peroxo), similar to H(peroxo) in methane monooxygenase. Although T201S generates T201S(peroxo) in addition to optically transparent ToMOH(peroxo), previously observed in wild-type ToMOH, this conservative variant is catalytically active in steady-state catalysis and single-turnover experiments and displays the same regiospecificity for toluene and slightly different regiospecificity for o-xylene oxidation.

Molecular determinants of the regioselectivity of toluene/o-xylene monooxygenase from Pseudomonas sp. strain OX1

Bacterial multicomponent monooxygenases (BMMs) are a heterogeneous family of di-iron monooxygenases which share the very interesting ability to hydroxylate aliphatic and/or aromatic hydrocarbons. Each BMM possesses defined substrate specificity and regioselectivity which match the metabolic requirements of the strain from which it has been isolated. Pseudomonas sp. strain OX1, a strain able to metabolize o-, m-, and p-cresols, produces the BMM toluene/o-xylene monooxygenase (ToMO), which converts toluene to a mixture of o-, m-, and p-cresol isomers. In order to investigate the molecular determinants of ToMO regioselectivity, we prepared and characterized 15 single-mutant and 3 double-mutant forms of the ToMO active site pocket. Using the Monte Carlo approach, we prepared models of ToMO-substrate and ToMO-reaction intermediate complexes which allowed us to provide a molecular explanation for the regioselectivities of wild-type and mutant ToMO enzymes. Furthermore, using binding energy values calculated by energy analyses of the complexes and a simple mathematical model of the hydroxylation reaction, we were able to predict quantitatively the regioselectivities of the majority of the variant proteins with good accuracy. The results show not only that the fine-tuning of ToMO regioselectivity can be achieved through a careful alteration of the shape of the active site but also that the effects of the mutations on regioselectivity can be quantitatively predicted a priori.

Structural consequences of effector protein complex formation in a diiron hydroxylase

Carboxylate-bridged diiron hydroxylases are multicomponent enzyme complexes responsible for the catabolism of a wide range of hydrocarbons and as such have drawn attention for their mechanism of action and potential uses in bioremediation and enzymatic synthesis. These enzyme complexes use a small molecular weight effector protein to modulate the function of the hydroxylase. However, the origin of these functional changes is poorly understood. Here, we report the structures of the biologically relevant effector protein-hydroxylase complex of toluene 4-monooxygenase in 2 redox states. The structures reveal a number of coordinated changes that occur up to 25 A from the active site and poise the diiron center for catalysis. The results provide a structural basis for the changes observed in a number of the measurable properties associated with effector protein binding. This description provides insight into the functional role of effector protein binding in all carboxylate-bridged diiron hydroxylases.

Characterization of the arene-oxidizing intermediate in ToMOH as a diiron(III) species

We report the generation and characterization of a diiron(III) intermediate formed during reaction with dioxygen of the reduced hydroxylase component of toluene/o-xylene monooxygenase from Pseudomonas sp. OX1. The decay rate of this species is accelerated upon mixing with phenol, a substrate for this system. Under steady-state conditions, hydrogen peroxide was generated in the absence of substrate. The oxidized hydroxylase also decomposed hydrogen peroxide to liberate dioxygen in the absence of reducing equivalents. This activity suggests that dioxygen activation may be reversible. The linear free energy relationship determined from hydroxylation of para-substituted phenols under steady-state turnover has a negative slope. A value of rho < 0 is consistent with electrophilic attack by the oxidizing intermediate on the aromatic substrates. The results from these steady and pre-steady-state experiments provide compelling evidence that the diiron(III) intermediate is the active oxidant in the toluene/o-xylene monooxygenase system and is a peroxodiiron(III) transient, despite differences between its optical and Mössbauer spectroscopic parameters and those of other peroxodiiron(III) centers.