4-Hydroxyphenylpyruvate Dioxygenase: A Hybrid Density Functional Study of the Catalytic Reaction Mechanism

Design, Synthesis, and Bioactivity of Novel 2-(Arylformyl)cyclohexane-1,3-dione Derivatives as HPPD Inhibitors

4-Hydroxyphenylpyruvate dioxygenase inhibitors (Echinochloa crus-galli 1.13.11.27, HPPD) have gained significant popularity as one of the best-selling herbicides worldwide. To identify highly effective HPPD inhibitors, a rational design approach utilizing bioisosterism was employed to create a series of 2-(arylformyl)cyclohexane-1,3-dione derivatives. A total of 29 novel compounds were synthesized and characterized through various techniques, including IR, 1H NMR, 13C NMR, and HRMS. Evaluation of their inhibitory activity against Arabidopsis thaliana HPPD (AtHPPD) revealed that certain derivatives exhibited superior potency compared to mesotrione (IC50 = 0.204 μM). Initial herbicidal activity tests demonstrated that compounds 27 and 28 were comparable to mesotrione in terms of weed control and crop safety, with compound 28 exhibiting enhanced safety in canola crops. Molecular docking analyses indicated that the quinoline rings of compounds 27 and 28 formed more stable π-π interactions with the amino acid residues Phe-360 and Phe-403 in the active cavity of AtHPPD, surpassing the benzene ring of mesotrione. Molecular dynamics simulations and molecular structure comparisons confirmed the robust binding capabilities of compounds 27 and 28 to AtHPPD. This study provides a valuable reference for the development of novel triketone herbicide structures, serving as a blueprint for future advancements in this field.

Discovery of novel HPPD inhibitors: Virtual screening, molecular design, structure modification and biological evaluation

4-Hydroxyphenylpyruvate dioxygenase (EC 1.13.11.27, HPPD, a Fe(II)/α-ketoglutarate dependent oxygenases), is a popular herbicide target. In this work, two pharmacophore models based on common molecular characteristics (HipHop) and receptor-ligand complex (CBP) were generated for virtual screening for HPPD inhibitors. About 1,000,000 molecules containing diketone structure from PubChem were filtered by Lipinski's rules to build a 3D database. Then the database was screened through combining HipHop model, CBP model, ADMET (absorption, distribution, metabolism, excretion and toxicity) prediction and molecular docking. Subsequently, based on the specific binding mode and affinity of HPPD inhibitors, 4 molecules with high -CDOCKER energy, good aqueous solubility and human safety predicative properties values were screened. From the screening results and combined with previous work, three novel HPPD inhibitors were designed and synthesized through fragment splicing and bioisosterism strategies. Compound IV-a exhibited similar inhibition of Arabidopsis thaliana HPPD (AtHPPD) and herbicidal activity as mesotrione. Crop selectivity showed that compound IV-a had better crop safety than mesotrione. Comparing the molecular properties, ADMET and molecular docking studies indicated that compounds IV-a exhibited better properties than mesotrione, which could be further modified as novel HPPD inhibitor herbicides.

Insights into 4-hydroxyphenylpyruvate dioxygenase-inhibitor interactions from comparative structural biology

4-Hydroxyphenylpyruvate dioxygenase (HPPD) plays a key role in tyrosine metabolism and has been identified as a promising target for herbicide and drug discovery. The structures of HPPD complexed with different types of inhibitors have been determined previously. We summarize the structures of HPPD complexed with structurally diverse molecules, including inhibitors, natural products, substrates, and catalytic intermediates; from these structures, the detailed inhibitory mechanisms of different inhibitors were analyzed and compared, and the key structural factors determining the slow-binding behavior of inhibitors were identified. Further, we propose four subpockets that accommodate different inhibitor substructures. We believe that these analyses will facilitate in-depth understanding of the enzymatic reaction mechanism and enable the design of new inhibitors with higher potency and selectivity.

New insights and advances on pyomelanin production: from microbial synthesis to applications

Pyomelanin is a brown-black phenolic polymer and results from the oxidation of homogentisic acid (HGA) in the L-tyrosine pathway. As part of the research for natural and active ingredients issued from realistic bioprocesses, this work re-evaluates the HGA pigment and makes an updated inventory of its syntheses, microbial pathways, and properties, with tracks and recent advances for its large-scale production. The mechanism of the HGA polymerization is also well documented. In alkaptonuria, pyomelanin formation leads to connective tissue damages and arthritis, most probably due to the ROS issued from HGA oxidation. While UV radiation on human melanin may generate degradation products, pyomelanin is not photodegradable, is hyperthermostable, and has other properties better than the L-Dopa melanin. This review aims to raise awareness about the potential of this pigment for various applications, not only for skin coloring and protection but also for other cells, materials, and as a promising (semi)conductor for bioelectronics and energy.

Rational Redesign of Enzyme via the Combination of Quantum Mechanics/Molecular Mechanics, Molecular Dynamics, and Structural Biology Study

Increasing demands for efficient and versatile chemical reactions have prompted innovations in enzyme engineering. A major challenge in engineering α-ketoglutarate-dependent oxygenases is to develop a rational strategy which can be widely used for directly evolving the desired mutant to generate new products. Herein, we report a strategy for rational redesign of a model enzyme, 4-hydroxyphenylpyruvate dioxygenase (HPPD), based on quantum mechanics/molecular mechanics (QM/MM) calculation and molecular dynamic simulations. This strategy enriched our understanding of the HPPD catalytic reaction pathway and led to the discovery of a series of HPPD mutants producing hydroxyphenylacetate (HPA) as the alternative product other than the native product homogentisate. The predicted HPPD-Fe(IV)═O-HPA intermediate was further confirmed by the crystal structure of Arabidopsis thaliana HPPD/S267W complexed with HPA. These findings not only provide a good understanding of the structure-function relationship of HPPD but also demonstrate a generally applicable platform for the development of biocatalysts.

Functional role of residues involved in substrate binding of human 4-hydroxyphenylpyruvate dioxygenase

4-Hydroxylphenylpyruvate dioxygenase (HPPD) catalyzes the conversion of 4-hydroxylphenylpyruvate (HPP) to homogentisate, the important step for tyrosine catabolism. Comparison of the structure of human HPPD with the substrate-bound structure of A. thaliana HPPD revealed notably different orientations of the C-terminal helix. This helix performed as a closed conformation in human enzyme. Simulation revealed a different substrate-binding mode in which the carboxyl group of HPP interacted by a H-bond network formed by Gln334, Glu349 (the metal-binding ligand), and Asn363 (in the C-terminal helix). The 4-hydroxyl group of HPP interacted with Gln251 and Gln265. The relative activity and substrate-binding affinity were preserved for the Q334A mutant, implying the alternative role of Asn363 for HPP binding and catalysis. The reduction in kcat/Km of the Asn363 mutants confirmed the critical role in catalysis. Compared to the N363A mutant, the dramatic reduction in the Kd and thermal stability of the N363D mutant implies the side-chain effect in the hinge region rotation of the C-terminal helix. The activity and binding affinity were not recovered by double mutation; however, the 4-hydroxyphenylacetate intermediate formation by the uncoupled reaction of Q334N/N363Q and Q334A/N363D mutants indicated the importance of the H-bond network in the electrophilic reaction. These results highlight the functional role of the H-bond network in a closed conformation of the C-terminal helix to stabilize the bound substrate. The extremely low activity and reduction in Q251E's Kd suggest that interaction coupled with the H-bond network is crucial to locate the substrate for nucleophilic reaction.

Discovery and structure optimization of a novel corn herbicide, tolpyralate

Tolpyralate, a new selective postemergence herbicide developed for the weed control in corn, possesses a unique chemical structure with a 1-alkoxyethyl methyl carbonate group on the N-ethyl pyrazole moiety. This compound shows high herbicidal activity against many weed species, including glyphosate-resistant Amaranthus tuberculatus. Tolpyralate targets 4-hydroxyphenylpyruvate dioxygenase (4-HPPD), which is involved in the tyrosine degradation pathway. Inhibition of the enzyme destroys the chlorophyll, thereby killing the susceptible weeds. Details of tolpyralate discovery, structure optimization, and biological activities are described.

Computational studies of DNA base repair mechanisms by nonheme iron dioxygenases: selective epoxidation and hydroxylation pathways

DNA base repair mechanisms of alkylated DNA bases is an important reaction in chemical biology and particularly in the human body. It is typically catalyzed by an α-ketoglutarate-dependent nonheme iron dioxygenase named the AlkB repair enzyme. In this work we report a detailed computational study into the structure and reactivity of AlkB repair enzymes with alkylated DNA bases. In particular, we investigate the aliphatic hydroxylation and C[double bond, length as m-dash]C epoxidation mechanisms of alkylated DNA bases by a high-valent iron(iv)-oxo intermediate. Our computational studies use quantum mechanics/molecular mechanics methods on full enzymatic structures as well as cluster models on active site systems. The work shows that the iron(iv)-oxo species is rapidly formed after dioxygen binding to an iron(ii) center and passes a bicyclic ring structure as intermediate. Subsequent cluster models explore the mechanism of substrate hydroxylation and epoxidation of alkylated DNA bases. The work shows low energy barriers for substrate activation and consequently energetically feasible pathways are predicted. Overall, the work shows that a high-valent iron(iv)-oxo species can efficiently dealkylate alkylated DNA bases and return them into their original form.

Discovery of Novel Pyrazole-Quinazoline-2,4-dione Hybrids as 4-Hydroxyphenylpyruvate Dioxygenase Inhibitors

4-Hydroxyphenylpyruvate dioxygenase (HPPD, EC 1.13.11.27) has been identified as one of the most significant targets in herbicide discovery for resistant weed control. In a continuing effort to discover potent novel HPPD inhibitors, we adopted a ring-expansion strategy to design a series of novel pyrazole-quinazoline-2,4-dione hybrids based on the previously discovered pyrazole-isoindoline-1,3-dione scaffold. One compound, 3-(2-chlorophenyl)-6-(5-hydroxy-1,3-dimethyl-1H-pyrazole-4-carbonyl)-1,5-dimethylquinazoline-2,4(1H,3H)-dione (9bj), displayed excellent potency against AtHPPD, with an IC₅₀ value of 84 nM, which is approximately 16-fold more potent than pyrasulfotole (IC₅₀ = 1359 nM) and 2.7-fold more potent than mesotrione (IC₅₀ = 226 nM). Furthermore, the co-crystal structure of the AtHPPD-9bj complex (PDB ID 6LGT) was determined at a resolution of 1.75 Å. Similar to the existing HPPD inhibitors, compound 9bj formed a bidentate chelating interaction with the metal ion and a π-π stacking interaction with Phe381 and Phe424. In contrast, o-chlorophenyl at the N3 position of quinazoline-2,4-dione with a double conformation was surrounded by hydrophobic residues (Met335, Leu368, Leu427, Phe424, Phe392, and Phe381). Remarkably, the greenhouse assay indicated that most compounds displayed excellent herbicidal activity (complete inhibition) against at least one of the tested weeds at the application rate of 150 g of active ingredient (ai)/ha. Most promisingly, compounds 9aj and 9bi not only exhibited prominent weed control effects with a broad spectrum but also showed very good crop safety to cotton, peanuts, and corn at the dose of 150 g of ai/ha.

Role of the N-terminus in human 4-hydroxyphenylpyruvate dioxygenase activity

4-Hydroxyphenylpyruvate dioxygenase (HPPD) is a key enzyme in tyrosine catabolism, catalysing the oxidation of 4-hydroxyphenylpyruvate to homogentisate. Genetic deficiency of this enzyme causes type III tyrosinaemia. The enzyme comprises two barrel-shaped domains formed by the N- and C-termini, with the active site located in the C-terminus. This study investigated the role of the N-terminus, located at the domain interface, in HPPD activity. We observed that the kcat/Km decreased ∼8-fold compared with wild type upon removal of the 12 N-terminal residues (ΔR13). Interestingly, the wild-type level of activity was retained in a mutant missing the 17 N-terminal residues, with a kcat/Km 11-fold higher than that of the ΔR13 mutant; however, the structural stability of this mutant was lower than that of wild type. A 2-fold decrease in catalytic efficiency was observed for the K10A and E12A mutants, indicating synergism between these residues in the enzyme catalytic function. A molecular dynamics simulation showed large RMS fluctuations in ΔR13 suggesting that conformational flexibility at the domain interface leads to lower activity in this mutant. These results demonstrate that the N-terminus maintains the stability of the domain interface to allow for catalysis at the active site of HPPD.

Synthesis and herbicidal activities of aryloxyacetic acid derivatives as HPPD inhibitors

A series of aryloxyacetic acid derivatives were designed and synthesized as 4-hydoxyphenylpyruvate dioxygenase (HPPD) inhibitors. Preliminary bioassay results reveal that these derivatives are promising Arabidopsis thaliana HPPD (AtHPPD) inhibitors, in particular compounds I12 (K_i = 0.011 µM) and I23 (K_i = 0.012 µM), which exhibit similar activities to that of mesotrione, a commercial HPPD herbicide (K_i = 0.013 µM). Furthermore, the newly synthesized compounds show significant greenhouse herbicidal activities against tested weeds at dosages of 150 g ai/ha. In particular, II4 exhibited high herbicidal activity for pre-emergence treatment that was slightly better than that of mesotrione. In addition, compound II4 was safe for weed control in maize fields at a rate of 150 g ai/ha, and was identified as the most potent candidate for a novel HPPD inhibitor herbicide. The compounds described herein may provide useful guidance for the design of new HPPD inhibiting herbicides and their modification.

Reaction Mechanism of Histone Demethylation in αKG-dependent Non-Heme Iron Enzymes

Histone demethylases (KDMs) catalyze histone lysine demethylation, an important epigenetic process that controls gene expression in eukaryotes, and represent important cancer drug targets for cancer treatment. Demethylation of histone is comprised of sequential reaction steps including oxygen activation, decarboxylation, and demethylation. The initial oxygen binding and activation steps have been studied. However, the information on the complete catalytic reaction cycle is limited, which has impeded the structure-based design of inhibitors targeting KDMs. Here we report the mechanism of the complete reaction steps catalyzed by a representative nonheme iron αKG-dependent KDM, PHF8 using QM/MM approaches. The atomic-level understanding on the complete reaction mechanism of PHF8 would shed light on the structure-based design of selective inhibitors targeting KDMs to intervene in cancer epigenetics.

Mechanistic insights into a non-heme 2-oxoglutarate-dependent ethylene-forming enzyme: selectivity of ethylene-formation versusl-Arg hydroxylation

The ethylene-forming enzyme (EFE) is a unique member of the Fe(ii)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenases. It converts 2OG into ethylene plus three CO2 molecules (ethylene-forming reaction) and also catalyzes the C5 hydroxylation of l-arginine coupled to the oxidative decarboxylation of 2OG (l-Arg hydroxylation reaction). To uncover the mechanisms of the dual transformations by EFE, quantum mechanical/molecular mechanical (QM/MM) calculations were carried out. Based on the results, a branched mechanism was proposed. An FeII-peroxysuccinate complex with a dissociated CO2 generated through the nucleophilic attack of the superoxo moiety of the Fe-O2 species on the keto carbon of 2OG is the key common intermediate in both reactions. A competition between the subsequent CO2 insertion (a key step in the ethylene-forming pathway) and the O-O bond cleavage (leading to the formation of succinate) governs the product selectivity. The calculated reaction barriers suggested that the CO2 insertion is favored over the O-O bond cleavage. This is consistent with the product preference observed in experiments. By comparison with the results of AsqJ (an Fe/2OG oxygenase that leads to substrate oxidation exclusively), the protein environment was found to be crucial for the selectivity. Further calculations demonstrated that the local electric field of the protein environment in EFE promotes ethylene formation by acting as a charge template, exemplifying the importance of the electrostatic interaction in enzyme catalysis. These findings offer mechanistic insights into the EFE catalysis and provide important clues for better understanding the unique ethylene-forming capability of EFE compared with other Fe/2OG oxygenases.

Advances in Sustainable Catalysis: A Computational Perspective

The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.

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.

A Comparative Review on the Catalytic Mechanism of Nonheme Iron Hydroxylases and Halogenases

Enzymatic halogenation and haloperoxidation are unusual processes in biology; however, a range of halogenases and haloperoxidases exist that are able to transfer an aliphatic or aromatic C–H bond into C–Cl/C–Br. Haloperoxidases utilize hydrogen peroxide, and in a reaction with halides (Cl−/Br−), they react to form hypohalides (OCl−/OBr−) that subsequently react with substrate by halide transfer. There are three types of haloperoxidases, namely the iron-heme, nonheme vanadium, and flavin-dependent haloperoxidases that are reviewed here. In addition, there are the nonheme iron halogenases that show structural and functional similarity to the nonheme iron hydroxylases and form an iron(IV)-oxo active species from a reaction of molecular oxygen with α-ketoglutarate on an iron(II) center. They subsequently transfer a halide (Cl−/Br−) to an aliphatic C–H bond. We review the mechanism and function of nonheme iron halogenases and hydroxylases and show recent computational modelling studies of our group on the hectochlorin biosynthesis enzyme and prolyl-4-hydroxylase as examples of nonheme iron halogenases and hydroxylases. These studies have established the catalytic mechanism of these enzymes and show the importance of substrate and oxidant positioning on the stereo-, chemo- and regioselectivity of the reaction that takes place.

Novel HPPD inhibitors: triketone 2H-benzo[b][1,4]oxazin-3(4H)-one analogs

BACKGROUND

Herbicides that inhibit 4-hydroxyphenylpyruvate dioxygenase (HPPD, EC 1.13.11.27) are very important for grass weed control. In order to discover novel HPPD herbicides, a series of triketone 2H-benzo[b] oxazin-3(4H)-one analogs was designed and synthesized.

RESULTS

In comparison with the commercial triketone HPPD herbicide mesotrione (IC₅₀ = 0.252 μM), some of these new triketone analogs displayed excellent HPPD inhibitory potency in vitro, for example B39 (IC₅₀ = 0.172 μM) and B41 (IC₅₀ = 0.156 μM). In addition, some of these compounds exhibited pre- and post-emergence herbicidal activity similar to mesotrione when applied at 375 g/ha.

CONCLUSION

Many of the title compounds described in this paper could be important lead structures for the further development of novel HPPD herbicides. © 2017 Society of Chemical Industry.

Combination of Virtual Screening Protocol by in Silico toward the Discovery of Novel 4-Hydroxyphenylpyruvate Dioxygenase Inhibitors

4-Hydroxyphenylpyruvate dioxygenase (EC 1.13.11.27, HPPD) is a potent new bleaching herbicide target. Therefore, in silico structure-based virtual screening was performed in order to speed up the identification of promising HPPD inhibitors. In this study, an integrated virtual screening protocol by combining 3D-pharmacophore model, molecular docking and molecular dynamics (MD) simulation was established to find novel HPPD inhibitors from four commercial databases. 3D-pharmacophore Hypo1 model was applied to efficiently narrow potential hits. The hit compounds were subsequently submitted to molecular docking studies, showing four compounds as potent inhibitor with the mechanism of the Fe(II) coordination and interaction with Phe360, Phe403, and Phe398. MD result demonstrated that nonpolar term of compound 3881 made great contributions to binding affinities. It showed an IC₅₀ being 2.49 μM against AtHPPD in vitro. The results provided useful information for developing novel HPPD inhibitors, leading to further understanding of the interaction mechanism of HPPD inhibitors.

4-Hydroxyphenylpyruvate Dioxygenase Inhibitors: From Chemical Biology to Agrochemicals

The development of new herbicides is receiving considerable attention to control weed biotypes resistant to current herbicides. Consequently, new enzymes are always desired as targets for herbicide discovery. 4-Hydroxyphenylpyruvate dioxygenase (HPPD, EC 1.13.11.27) is an enzyme engaged in photosynthetic activity and catalyzes the transformation of 4-hydroxyphenylpyruvic acid (HPPA) into homogentisic acid (HGA). HPPD inhibitors constitute a promising area of discovery and development of innovative herbicides with some advantages, including excellent crop selectivity, low application rates, and broad-spectrum weed control. HPPD inhibitors have been investigated for agrochemical interests, and some of them have already been commercialized as herbicides. In this review, we mainly focus on the chemical biology of HPPD, discovery of new potential inhibitors, and strategies for engineering transgenic crops resistant to current HPPD-inhibiting herbicides. The conclusion raises some relevant gaps for future research directions.

3D Pharmacophore-Based Virtual Screening and Docking Approaches toward the Discovery of Novel HPPD Inhibitors

p-Hydroxyphenylpyruvate dioxygenase (HPPD) is not only the useful molecular target in treating life-threatening tyrosinemia type I, but also an important target for chemical herbicides. A combined in silico structure-based pharmacophore and molecular docking-based virtual screening were performed to identify novel potential HPPD inhibitors. The complex-based pharmacophore model (CBP) with 0.721 of ROC used for screening compounds showed remarkable ability to retrieve known active ligands from among decoy molecules. The ChemDiv database was screened using CBP-Hypo2 as a 3D query, and the best-fit hits subjected to molecular docking with two methods of LibDock and CDOCKER in Accelrys Discovery Studio 2.5 (DS 2.5) to discern interactions with key residues at the active site of HPPD. Four compounds with top rankings in the HipHop model and well-known binding model were finally chosen as lead compounds with potential inhibitory effects on the active site of target. The results provided powerful insight into the development of novel HPPD inhibitors herbicides using computational techniques.

Atom Tunneling in the Hydroxylation Process of Taurine/α-Ketoglutarate Dioxygenase Identified by Quantum Mechanics/Molecular Mechanics Simulations

Taurine/α-ketoglutarate dioxygenase is one of the most studied α-ketoglutarate-dependent dioxygenases (αKGDs), involved in several biotechnological applications. We investigated the key step in the catalytic cycle of the αKGDs, the hydrogen transfer process, by a quantum mechanics/molecular mechanics approach (B3LYP/CHARMM22). Analysis of the charge and spin densities during the reaction demonstrates that a concerted mechanism takes place, where the H atom transfer happens simultaneously with the electron transfer from taurine to the Fe═O cofactor. We found the quantum tunneling of the hydrogen atom to increase the rate constant by a factor of 40 at 5 °C. As a consequence, a quite high kinetic isotope effect close to 60 is obtained, which is consistent with the experimental value.

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.

4-Hydroxyphenylpyruvate Dioxygenase and Its Inhibition in Plants and Animals: Small Molecules as Herbicides and Agents for the Treatment of Human Inherited Diseases

This review mainly focuses on the physiological function of 4-hydroxyphenylpyruvate dioxygenase (HPPD), as well as on the development and application of HPPD inhibitors of several structural classes. Among them, one illustrative example is represented by compounds belonging to the class of triketone compounds. They were discovered by serendipitous observations on weed growth and were developed as bleaching herbicides. Informed reasoning on nitisinone (NTBC, 14), a triketone that failed to reach the final steps of the herbicidal design and development process, allowed it to become a curative agent for type I tyrosinemia (T1T) and to enter clinical trials for alkaptonuria. These results boosted the research of new compounds able to interfere with HPPD activity to be used for the treatment of the tyrosine metabolism-related diseases.

Go it alone: four-electron oxidations by mononuclear non-heme iron enzymes

This review discusses the current mechanistic understanding of a group of mononuclear non-heme iron-dependent enzymes that catalyze four-electron oxidation of their organic substrates without the use of any cofactors or cosubstrates. One set of enzymes acts on α-ketoacid-containing substrates, coupling decarboxylation to oxygen activation. This group includes 4-hydroxyphenylpyruvate dioxygenase, 4-hydroxymandelate synthase, and CloR involved in clorobiocin biosynthesis. A second set of enzymes acts on substrates containing a thiol group that coordinates to the iron. This group is comprised of isopenicillin N synthase, thiol dioxygenases, and enzymes involved in the biosynthesis of ergothioneine and ovothiol. The final group of enzymes includes HEPD and MPnS that both carry out the oxidative cleavage of the carbon-carbon bond of 2-hydroxyethylphosphonate but generate different products. Commonalities amongst many of these enzymes are discussed and include the initial substrate oxidation by a ferric-superoxo-intermediate and a second oxidation by a ferryl species.

Catalytic strategies of the non-heme iron dependent oxygenases and their roles in plant biology

•

Current evidence for iron-oxo reactive intermediates is reviewed.

•

In crystallo intermediates detected in a native extradiol dioxygenase reaction.

•

Carotenoid cleavage dioxygenases catalyse strigolactone biosynthesis.

•

Identification of plant cysteine oxidases involved in the plant hypoxic response.

•

Applications of enzyme manipulation to plant biology and agriculture are discussed.

Non-heme iron-dependent oxygenases catalyse the incorporation of O₂ into a wide range of biological molecules and use diverse strategies to activate their substrates. Recent kinetic studies, including in crystallo, have provided experimental support for some of the intermediates used by different subclasses of this enzyme family. Plant non-heme iron-dependent oxygenases have diverse and important biological roles, including in growth signalling, stress responses and secondary metabolism. Recently identified roles include in strigolactone biosynthesis, O-demethylation in morphine biosynthesis and regulating the stability of hypoxia-responsive transcription factors. We discuss current structural and mechanistic understanding of plant non-heme iron oxygenases, and how their chemical/genetic manipulation could have agricultural benefit, for example, for improved yield, stress tolerance or herbicide development.

Mechanism of O2 Activation by α-Ketoglutarate Dependent Oxygenases Revisited. A Quantum Chemical Study

Four mechanisms previously proposed for dioxygen activation catalyzed by α-keto acid dependent oxygenases (α-KAO) were studied with dispersion-corrected DFT methods employing B3LYP and TPSSh functionals in combination with triple-ζ basis set (cc-pVTZ). The aim of this study was to revisit mechanisms suggested in the past decade and resolve remaining issues related to dioxygen activation. Mechanism A, which runs on the quintet potential energy surface (PES) and includes formation of an Fe(III)-superoxide radical anion complex, subsequent oxidative decarboxylation, and O-O bond cleavage, was found to be most likely. However, mechanism B taking place on the septet PES involves a rate limiting barrier comparable to the one found for mechanism A, and thus it cannot be excluded, though two other mechanisms (C and D) were ruled out. Mechanism C is a minor variation of mechanism A, whereas mechanism D proceeds through formation of a triplet Fe(IV)-alkyl peroxo bridged intermediate. The study covered also full optimization of relevant minimum energy crossing points (MECPs). The relative energy of critical intermediates was also studied with the CCSD(T) method in order to benchmark TPSSh and B3LYP functionals with respect to their credibility in predicting relative energies of septet and triplet spin states of the ternary enzyme-Fe-α-keto glutarate (α-KG)-O2 complex.

Quantum Mechanics/Molecular Mechanics Modeling of Enzymatic Processes: Caveats and Breakthroughs

Nature has developed large groups of enzymatic catalysts with the aim to transfer substrates into useful products, which enables biosystems to perform all their natural functions. As such, all biochemical processes in our body (we drink, we eat, we breath, we sleep, etc.) are governed by enzymes. One of the problems associated with research on biocatalysts is that they react so fast that details of their reaction mechanisms cannot be obtained with experimental work. In recent years, major advances in computational hardware and software have been made and now large (bio)chemical systems can be studied using accurate computational techniques. One such technique is the quantum mechanics/molecular mechanics (QM/MM) technique, which has gained major momentum in recent years. Unfortunately, it is not a black-box method that is easily applied, but requires careful set-up procedures. In this work we give an overview on the technical difficulties and caveats of QM/MM and discuss work-protocols developed in our groups for running successful QM/MM calculations.

4-Hydroxyphenylpyruvate Dioxygenase and Hydroxymandelate Synthase: 2-Oxo Acid-Dependent Oxygenases of Importance to Agriculture and Medicine

Despite a separate evolutionary lineage, 4-hydroxyphenylpyruvate dioxygenase (HPPD) and hydroxymandelate synthase (HMS) are appropriately grouped with the 2-oxo acid-dependent oxygenase (2OADO) family of enzymes. HPPD and HMS accomplish highly similar overall chemistry to that observed in the majority of 2OADOs but require only two substrates rather than three. 2OADOs typically use the 2-oxo acid of 2-oxoglutarate (2OG) as a source of electrons to reduce and activate dioxygen in order to oxidize a third specific substrate. HPPD and HMS use instead the pyruvate substituent of 4-hydroxyphenylpyruvate to activate dioxygen and then proceed to also hydroxylate this substrate, each yielding a distinctly different aromatic product. HPPD catalyses the second and committed step of tyrosine catabolism, a pathway common to nearly all aerobes. Plants require the HPPD reaction to biosynthesize plastoquinones and therefore HPPD inhibitors can have potent herbicidal activity. The ubiquity of the HPPD reaction, however, has meant that HPPD-specific molecules developed as herbicides have other uses in different forms of life. In humans herbicidal HPPD inhibitors can be used therapeutically to alleviate specific inborn defects and also to retard the progress of certain bacterial and fungal infections. This review is intended as a concise overview of the contextual and catalytic chemistries of HPPD and HMS.

Dioxygen activation and two consecutive oxidative decarboxylations of phenylpyruvate by nonheme iron(II) complexes: functional models of hydroxymandelate synthase (HMS) and CloR

Two mononuclear iron(ii)-phenylpyruvate complexes of monoanionic facial N3 ligands are reported to react with dioxygen to undergo two consecutive oxidative decarboxylation steps via an iron-mandelate complex mimicking the function of HMS and CloR.

Synthesis and herbicidal evaluation of triketone-containing quinazoline-2,4-diones

Exploring novel 4-hydroxyphenylpyruvate dioxygenase (EC 1.13.11.27, HPPD) inhibitors is one of the most promising research directions in herbicide discovery. To discover new triketone herbicides with broad-spectrum weed control as well as excellent crop selectivity, a series of (total 52) novel triketone-containing quinazoline-2,4-dione derivatives were synthesized and further bioevaluated. The greenhouse testing indicated that many of the newly synthesized compounds showed better or excellent herbicidal activity against broadleaf and monocotyledonous weeds at the dosages of 37.5-150 g of active ingredient (ai)/ha. The structure and activity relationship in this study indicated that the triketone-containing quinazoline-2,4-dione motif has possessed great impact on herbicide activity and may be used for further optimization. Among the new compounds, III-b and VI-a-VI-d displayed a broader spectrum of weed control than mesotrione. In addition, the compound III-b also demonstrated comparatively superior crop selectivity to mesotrione, thus possessing great potential for weed control in the field.

Insights into enzymatic halogenation from computational studies

The halogenases are a group of enzymes that have only come to the fore over the last 10 years thanks to the discovery and characterization of several novel representatives. They have revealed the fascinating variety of distinct chemical mechanisms that nature utilizes to activate halogens and introduce them into organic substrates. Computational studies using a range of approaches have already elucidated many details of the mechanisms of these enzymes, often in synergistic combination with experiment. This Review summarizes the main insights gained from these studies. It also seeks to identify open questions that are amenable to computational investigations. The studies discussed herein serve to illustrate some of the limitations of the current computational approaches and the challenges encountered in computational mechanistic enzymology.

4-Hydroxyphenylpyruvate dioxygenase and hydroxymandelate synthase: exemplars of the α-keto acid dependent oxygenases

4-Hydroxyphenylpyruvate dioxygenase (HPPD) and hydroxymandelate synthase (HMS) are outliers within the α-keto acid dependent oxygenase (αKAO) family. HPPD and HMS catalyze the chemistry of the majority of enzymes within the αKAO family but are clearly mechanistically convergent, having a grossly different structural topology. Some of the unique characteristics of HPPD and HMS have elucidated select parts of the catalytic cycle that are obscured in other family members. Moreover, the inhibitory chemistry of HPPD is a phenomenon with ever-expanding relevance across all kingdoms of life. This review is a synopsis of the literature pertaining to HPPD and HMS. It is not intended as an exhaustive compilation of all observations made for these enzymes but rather a condensed narrative that connects those studies that have advanced the understanding of the chemistry of both enzymes.

Chiral hydroxylation at the mononuclear nonheme Fe(II) center of 4-(S) hydroxymandelate synthase--a structure-activity relationship analysis

(S)-Hydroxymandelate synthase (Hms) is a nonheme Fe(II) dependent dioxygenase that catalyzes the oxidation of 4-hydroxyphenylpyruvate to (S)-4-hydroxymandelate by molecular oxygen. In this work, the substrate promiscuity of Hms is characterized in order to assess its potential for the biosynthesis of chiral α-hydroxy acids. Enzyme kinetic analyses, the characterization of product spectra, quantitative structure activity relationship (QSAR) analyses and in silico docking studies are used to characterize the impact of substrate properties on particular steps of catalysis. Hms is found to accept a range of α-oxo acids, whereby the presence of an aromatic substituent is crucial for efficient substrate turnover. A hydrophobic substrate binding pocket is identified as the likely determinant of substrate specificity. Upon introduction of a steric barrier, which is suspected to obstruct the accommodation of the aromatic ring in the hydrophobic pocket during the final hydroxylation step, the racemization of product is obtained. A steady state kinetic analysis reveals that the turnover number of Hms strongly correlates with substrate hydrophobicity. The analysis of product spectra demonstrates high regioselectivity of oxygenation and a strong coupling efficiency of C-C bond cleavage and subsequent hydroxylation for the tested substrates. Based on these findings the structural basis of enantioselectivity and enzymatic activity is discussed.

Applications of Potential Energy Surfaces in the Study of Enzymatic Reactions

From a generated PES, one can determine the relative energies of species involved, the sequence in which they occur, and the activation barrier(s) associated with individual steps or the overall mechanism. Furthermore, they can provide more insights than a simple indication of a path of sequential mechanistic structures and their energetic relationships. The investigation into the activation of O2 by alpha-ketoglutarate-dependent dioxygenase (AlkB) clearly shows the opportunity for spin inversion, where one can see that the lowest energy product may be formed via several possible routes. In the investigation of uroporphyrinogen decarboxylase III (UROD), the use of QM/MM methods allowed for the inclusion of the anisotropic protein environment providing greater insight into the rate-limiting barrier. Lastly, the mechanism of 6-phospho-α-glucosidase (GlvA) was discussed using different active site models. In particular, a continuum model PES was compared to the gas-phase PES.

Mechanism of benzylic hydroxylation by 4-hydroxymandelate synthase. A computational study

Hydroxymandelate synthase (HMS) and 4-hydroxyphenylpyruvate dioxygenase (HPPD) are highly related enzymes using the same substrates but catalyzing hydroxylation reactions yielding different products. The first steps of the HMS and HPPD catalytic reactions are believed to proceed in the same way and lead to an Fe(IV)═O-hydroxyphenylacetate (HPA) intermediate. Further down the catalytic cycles, HMS uses Fe(IV)═O to perform hydroxylation of the benzylic carbon, whereas in HPPD, the reactive oxoferryl intermediate attacks the aromatic ring of HPA. This study focuses on this part of the HMS catalytic cycle that starts from the oxoferryl intermediate and aims to identify interactions within the active site that are responsible for enzyme specificity. To this end, a HMS-Fe(IV)═O-HPA complex was modeled with molecular dynamics simulations. On the basis of the molecular dynamics-equilibrated structure, an active site model suitable for quantum chemical investigations was constructed and used for density functional theory (B3LYP) calculations of the mechanism of the native reaction of HMS, i.e., benzylic hydroxylation, and the alternative electrophilic attack on the ring, which is a step of the HPPD catalytic cycle. The most important result of this study is the finding that the conformation of the Ser201 side chain in the second coordination shell has a key role in directing the reaction of Fe(IV)═O into either the HMS or the HPPD channel.