On the Catalytic Mechanism of (S)-2-Hydroxypropylphosphonic Acid Epoxidase (HppE): A Hybrid DFT Study

How Do Metalloproteins Tame the Fenton Reaction and Utilize •OH Radicals in Constructive Manners?

This Account describes the manner whereby nature controls the Fenton-type reaction of O-O homolysis of hydrogen peroxide and harnesses it to carry out various useful oxidative transformations in metalloenzymes. H₂O₂ acts as the cosubstrate for the heme-dependent peroxidases, P450BM3, P450_SPα, P450_BSβ, and the P450 decarboxylase OleT, as well as the nonheme enzymes HppE and the copper-dependent lytic polysaccharide monooxygenases (LPMOs). Whereas heme peroxidases use the Poulos-Kraut heterolytic mechanism for H₂O₂ activation, some heme enzymes prefer the alternative Fenton-type mechanism, which produces •OH radical intermediates. The fate of the •OH radical is controlled by the protein environment, using tight H-bonding networks around H₂O₂. The so-generated •OH radical is constrained by the surrounding H-bonding interactions, the orientation of which is targeted to perform H-abstraction from the Fe(III)-OH group and thereby leading to the formation of the active species, called Compound I (Cpd I), Por^+•Fe(IV)═O, which performs oxidation of the substrate. Alternatively, for the nonheme HppE enzyme, the O-O homolysis catalyzed by the resting state Fe(II) generates an Fe(III)-OH species that effectively constrains the •OH radical species by a tight H-bonding network. The so-formed H-bonded •OH radical acts directly as the oxidant, since it is oriented to perform H-abstraction from the C-H bond of the substrate (S)-2-HPP. The Fenton-type H₂O₂ activation is strongly suggested by computations to occur also in copper-dependent LPMOs and pMMO. In LPMOs, the Cu(I)-catalyzed O-O homolysis of the H₂O₂ cosubstrate generates an •OH radical that abstracts a hydrogen atom from Cu(II)-OH and forms thereby the active species of the enzyme, Cu(II)-O•. Such Fenton-type O-O activation can be shared by both the O₂-dependent activations of LPMOs and pMMOs, in which the O₂ cosubstrate may be reduced to H₂O₂ by external reductants. Our studies show that, generally, the H₂O₂ activation is highly dependent on the protein environment, as well as on the presence/absence of substrates. Since H₂O₂ is a highly flexible and hydrophilic molecule, the absence of suitable substrates may lead to unproductive binding or even to the release of H₂O₂ from the active site, as has been suggested in P450cam and LPMOs, whereas the presence of the substrate seems to play a role in steering a Fenton-type H₂O₂ activation. In the absence of a substrate, the hydrophilic active site of P450BM3 disfavors the binding and activation of H₂O₂ and protects thereby the enzyme from the damage by the Fenton reaction. Due to the distinct coordination and reaction environment, the Fenton-type H₂O₂ activation mechanism by enzymes differs from the reaction in synthetic systems. In nonenzymatic reactions, the H-bonding networks are quite dynamic and flexible and the reactivity of H₂O₂ is not strategically constrained as in the enzymatic environment. As such, our Account describes the controlled Fenton-type mechanism in metalloenzymes, and the role of the protein environment in constraining the •OH radical against oxidative damage, while directing it to perform useful oxidative transformations.

Computational Study of the C5-Hydroxylation Mechanism Catalyzed by the Diiron Monooxygenase PtmU3 as Part of the Platensimycin Biosynthesis

PtmU3 is a newly identified nonheme diiron monooxygenase, which installs a C-5 β-hydroxyl group into the C-19 CoA-ester intermediate involved in the biosynthesis of unique diterpene-derived scaffolds of platensimycin and platencin. PtmU3 possesses a noncanonical diiron active site architecture of a saturated six-coordinate iron center and lacks the μ-oxo bridge. Although the hydroxylation process is a simple reaction for nonheme mononuclear iron-dependent enzymes, how PtmU3 employs the diiron center to catalyze the H-abstraction and OH-rebound is still unknown. In particular, the electronic characteristic of diiron is also unclear. To understand the catalytic mechanism of PtmU3, we constructed two reactant models in which both the Fe1^II-Fe2^III-superoxo and Fe1^II-Fe2^IV═O are considered to trigger the H-abstraction and performed a series of quantum mechanics/molecular mechanics calculations. Our calculation results reveal that PtmU3 is a special monooxygenase, that is, both atoms of the dioxygen molecule can be incorporated into two molecules of the substrate by the successive reactions. In the first-round reaction, PtmU3 uses the Fe1^II-Fe2^III-superoxo to install a hydroxyl group into the substrate, generating the high-reactive Fe1^II-Fe2^IV═O complex. In the second-round reaction, the Fe1^II-Fe2^IV═O species is responsible for the hydroxylation of another molecule of the substrate. In the diiron center, Fe2 adopts the high spin state (S = 5/2) during the catalysis, whereas for Fe1, in addition to its structural role, it may also play an assistant role for Fe1 catalysis. In the two successive OH-installing steps, the H-abstraction is always the rate-liming step. E241 and D308 not only act as bridging ligands to connect two Fe ions but also take part in the electron reorganization. Owing to the high reactivity of Fe1^II-Fe2^IV═O compared to Fe1^II-Fe2^III-superoxo, besides the C5-hydroxylation, the C3- or C18-hydroxylation was also calculated to be feasible.

Mechanistic Insights into a Stibene Cleavage Oxygenase NOV1 from Quantum Mechanical/Molecular Mechanical Calculations

NOV1, a stilbene cleavage oxygenase, catalyzes the cleavage of the central double bond of stilbenes to two phenolic aldehydes, using a 4-His Fe(II) center and dioxygen. Herein, we use in-protein quantum mechanical/molecular mechanical (QM/MM) calculations to elucidate the reaction mechanism of the central double bond cleavage of phytoalexin resveratrol by NOV1. Our results showed that the oxygen molecule prefers to bind to the iron center in a side-on fashion, as suggested from the experiment. The quintet Fe-O2 complex with the side-on superoxo antiferromagnetic coupled to the resveratrol radical is identified as the reactive oxygen species. The QM/MM results support the dioxygenase mechanism involving a dioxetane intermediate with a rate-limiting barrier of 10.0 kcal mol-1. The alternative pathway through an epoxide intermediate is ruled out due to a larger rate-limiting barrier (26.8 kcal mol-1). These findings provide important insight into the catalytic mechanism of carotenoid cleavage oxygenases and also the dioxygen activation of non-heme enzymes.

Deciphering the chemoselectivity of nickel-dependent quercetin 2,4-dioxygenase

QM/MM calculations were performed to elucidate the reaction mechanism and chemoselectivity of 2,4-QueD. The protonation state of the first-shell ligand Glu74 plays an important role in dictating the selectivity.

Methylcobalamin-Dependent Radical SAM C-Methyltransferase Fom3 Recognizes Cytidylyl-2-hydroxyethylphosphonate and Catalyzes the Nonstereoselective C-Methylation in Fosfomycin Biosynthesis

A methylcobalamin (MeCbl)-dependent radical S-adenosyl-l-methionine (SAM) methyltransferase Fom3 was found to catalyze the C-methylation of cytidylyl-2-hydroxyethylphosphonate (HEP-CMP) to give cytidylyl-2-hydroxypropylphosphonate (HPP-CMP), although it was originally proposed to catalyze the C-methylation of 2-hydroxyethylphosphonate to give 2-hydroxypropylphosphonate in the biosynthesis of a unique C-P bond containing antibiotic fosfomycin in Streptomyces. Unexpectedly, the Fom3 reaction product from HEP-CMP was almost a 1:1 diastereomeric mixture of HPP-CMP, indicating that the C-methylation is not stereoselective. Presumably, only the CMP moiety of HEP-CMP is critical for substrate recognition; on the other hand, the enzyme does not fix the 2-hydroxy group of the substrate and either of the prochiral hydrogen atoms at the C2 position can be abstracted by the 5'-deoxyadenosyl radical generated from SAM to form the substrate radical intermediates, which react with MeCbl to afford the corresponding products. This strict substrate recognition mechanism with no stereoselectivity of a MeCbl-dependent radical SAM methyltransferase is remarkable in natural product biosynthetic chemistry, because such a hidden clue for selective substrate recognition is likely to be found in the other biosynthetic pathways.

Theoretical Study of the Mechanism of the Nonheme Iron Enzyme EgtB

EgtB is a nonheme iron enzyme catalyzing the C-S bond formation between γ-glutamyl cysteine (γGC) and N-α-trimethyl histidine (TMH) in the ergothioneine biosynthesis. Density functional calculations were performed to elucidate and delineate the reaction mechanism of this enzyme. Two different mechanisms were considered, depending on whether the sulfoxidation or the S-C bond formation takes place first. The calculations suggest that the S-O bond formation occurs first between the thiolate and the ferric superoxide, followed by homolytic O-O bond cleavage, very similar to the case of cysteine dioxygenase. Subsequently, proton transfer from a second-shell residue Tyr377 to the newly generated iron-oxo moiety takes place, which is followed by proton transfer from the TMH imidazole to Tyr377, facilitated by two crystallographically observed water molecules. Next, the S-C bond is formed between γGC and TMH, followed by proton transfer from the imidazole CH moiety to Tyr377, which was calculated to be the rate-limiting step for the whole reaction, with a barrier of 17.9 kcal/mol in the quintet state. The calculated barrier for the rate-limiting step agrees quite well with experimental kinetic data. Finally, this proton is transferred back to the imidazole nitrogen to form the product. The alternative thiyl radical attack mechanism has a very high barrier, being 25.8 kcal/mol, ruling out this possibility.

How do Enzymes Utilize Reactive OH Radicals? Lessons from Nonheme HppE and Fenton Systems

The iron(IV)-oxo (ferryl) intermediate has been amply established as the principal oxidant in nonheme enzymes and the key player in C-H bond activations and functionalizations. In contrast to this status, our present QM/MM calculations of the mechanism of fosfomycin biosynthesis (a broad range antibiotic) by the nonheme HppE enzyme rule out the iron(IV)-oxo as the reactive species in the hydrogen abstraction (H-abstraction) step of the pro-R hydrogen from the (S)-2-hydroxypropylphosphonic substrate. Moreover, the study reveals that the ferryl species is bypassed in HppE, while the actual oxidant is an HO(•) radical hydrogen-bonded to a ferric-hydroxo complex, resulting via the homolytic dissociation of the hydrogen peroxide complex, Fe(II)-H2O2. The computed energy barrier of this pathway is 12.0 kcal/mol, in fair agreement with the experimental datum of 9.8 kcal/mol. An alternative mechanism involves the iron-complexed hydroxyl radical (Fe(III)-(HO(•))) that is obtained by protonation of the iron(IV)-oxo group via the O-H group of the substrate. The barrier for this pathway, 23.0 kcal/mol, is higher than the one in the first mechanism. In both mechanisms, the HO(•) radical is highly selective; its H-abstraction leading to the final cis-fosfomycin product. It appears that HppE is prone to usage of HO(•) radicals for C-H bond activation, because the ferryl oxidant requires a specific H-abstraction trajectory (∠FeOH ∼ 180°) that cannot be met for intramolecular H-abstraction. Thus, this work broadens the landscape of nonheme iron enzymes and makes a connection to Fenton chemistry, with implications on new potential biocatalysts that may harness hydroxyl radicals for C-H bond functionalizations.

Recent application of calculations of metal complexes based on density functional theory

Recent application of density functional theory (DFT) for metal complexes is reviewed to show the achievements of DFT and the challenges for it, as well as the methods for selecting proper functionals.

Enantiopure phosphonic acids as chiral inducers: homochiral crystallization of cobalt coordination polymers showing field-induced slow magnetization relaxation

Homochiral crystallization of (M)- or (P)-Co(SO₄)(1,3-bbix)(H₂O)₃ (1M or 1P) from an achiral precursor is achieved in the presence of a catalytic amount of (S)- or (R)-3-phenyl-2-((phosphonomethyl)amino) propanoic acid.

Ferryl protonation in oxoiron(IV) porphyrins and its role in oxygen transfer

Ferryl porphyrins, P-Fe(IV)═O, are central reactive intermediates in the catalytic cycles of numerous heme proteins and a variety of model systems. There has been considerable interest in elucidating factors, such as terminal oxo basicity, that may control ferryl reactivity. Here, the sulfonated, water-soluble ferryl porphyrin complexes tetramesitylporphyrin, oxoFe(IV)TMPS (FeTMPS-II), its 2,6-dichlorophenyl analogue, oxoFe(IV)TDClPS (FeTDClPS-II), and two other analogues are shown to be protonated under turnover conditions to produce the corresponding bis-aqua-iron(III) porphyrin cation radicals. The results reveal a novel internal electromeric equilibrium, P-Fe(IV)═O ⇆ P(+)-Fe(III)(OH2)2. Reversible pKa values in the range of 4-6.3 have been measured for this process by pH-jump, UV-vis spectroscopy. Ferryl protonation has important ramifications for C-H bond cleavage reactions mediated by oxoiron(IV) porphyrin cation radicals in protic media. Both solvent O-H and substrate C-H deuterium kinetic isotope effects are observed for these reactions, indicating that hydrocarbon oxidation by these oxoiron(IV) porphyrin cation radicals occurs via a solvent proton-coupled hydrogen atom transfer from the substrate that has not been previously described. The effective FeO-H bond dissociation energies for FeTMPS-II and FeTDClPS-II were estimated from similar kinetic reactivities of the corresponding oxoFe(IV)TMPS(+) and oxoFe(IV)TDClPS(+) species to be ∼92-94 kcal/mol. Similar values were calculated from the two-proton P(+)-Fe(III)(OH2)2 pKa(obs) and the porphyrin oxidation potentials, despite a 230 mV range for the iron porphyrins examined. Thus, the iron porphyrin with the lower ring oxidation potential has a compensating higher basicity of the ferryl oxygen. The solvent-derived proton adds significantly to the driving force for C-H bond scission.

A DFT and ONIOM study of C–H hydroxylation catalyzed by nitrobenzene 1,2-dioxygenase

The DFT and ONIOM calculations show that C–H hydroxylation by nitrobenzene 1,2-dioxygenase proceeds through a HO–Fe^VO intermediate.

Not an Oxidase, But a Peroxidase

A key step in the biosynthesis of the antibiotic fosfomycin requires hydrogen peroxide, rather than molecular oxygen as previously assumed. [Also see Report by Wang et al. ]

Evidence that the fosfomycin-producing epoxidase, HppE, is a non-heme-iron peroxidase

The iron-dependent epoxidase HppE converts (S)-2-hydroxypropyl-1-phosphonate (S-HPP) to the antibiotic fosfomycin [(1R,2S)-epoxypropylphosphonate] in an unusual 1,3-dehydrogenation of a secondary alcohol to an epoxide. HppE has been classified as an oxidase, with proposed mechanisms differing primarily in the identity of the O2-derived iron complex that abstracts hydrogen (H•) from C1 of S-HPP to initiate epoxide ring closure. We show here that the preferred cosubstrate is actually H2O2 and that HppE therefore almost certainly uses an iron(IV)-oxo complex as the H• abstractor. Reaction with H2O2 is accelerated by bound substrate and produces fosfomycin catalytically with a stoichiometry of unity. The ability of catalase to suppress the HppE activity previously attributed to its direct utilization of O2 implies that reduction of O2 and utilization of the resultant H2O2 were actually operant.

Theoretical study of the catalytic mechanism of E1 subunit of pyruvate dehydrogenase multienzyme complex from Bacillus stearothermophilus

Pyruvate dehydrogenase multienzyme complex (PDHc) is a member of a family of 2-oxo acid dehydrogenase (OADH) multienzyme complexes involved in several central points of oxidative metabolism, and the E1 subunit is the most important component in the entire PDHc catalytic system, which catalyzes the reversible transfer of an acetyl group from a pyruvate to the lipoyl group of E2 subunit lipoly domain. In this article, the catalytic mechanism of the E1 subunit has been systematically studied using density functional theory (DFT). Four possible pathways with different general acid/base catalysts in decarboxylation and reductive acylation processes were explored. Our calculation results indicate that the 4'-amino pyrimidine of ThDP and residue His128 are the most likely proton donors in the decarboxylation and reductive acylation processes, respectively. During the reaction, each C-C and C-S bond formation or cleavage process, except for the liberation of CO2, is always accompanied by a proton transfer between the substrates and proton donors. The liberation of CO2 is calculated to be the rate-limiting step for the overall reaction, with an energy barrier of 13.57 kcal/mol. The decarboxylation process is endothermic by 5.32 kcal/mol, whereas the reductive acylation process is exothermic with a value of 5.74 kcal/mol. The assignment of protonation states of the surrounding residues can greatly influence the reaction. Residues His128 and His271 play roles in positioning the first substrate pyruvate and second substrate lipoyl group, respectively.