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

An Active Site Tyr Residue Guides the Regioselectivity of Lysine Hydroxylation by Nonheme Iron Lysine-4-hydroxylase Enzymes through Proton-Coupled Electron Transfer

Lysine dioxygenase (KDO) is an important enzyme in human physiology involved in bioprocesses that trigger collagen cross-linking and blood pressure control. There are several KDOs in nature; however, little is known about the factors that govern the regio- and stereoselectivity of these enzymes. To understand how KDOs can selectively hydroxylate their substrate, we did a comprehensive computational study into the mechanisms and features of 4-lysine dioxygenase. In particular, we selected a snapshot from the MD simulation on KDO5 and created large QM cluster models (A, B, and C) containing 297, 312, and 407 atoms, respectively. The largest model predicts regioselectivity that matches experimental observation with rate-determining hydrogen atom abstraction from the C4-H position, followed by fast OH rebound to form 4-hydroxylysine products. The calculations show that in model C, the dipole moment is positioned along the C4-H bond of the substrate and, therefore, the electrostatic and electric field perturbations of the protein assist the enzyme in creating C4-H hydroxylation selectivity. Furthermore, an active site Tyr233 residue is identified that reacts through proton-coupled electron transfer akin to the axial Trp residue in cytochrome c peroxidase. Thus, upon formation of the iron(IV)-oxo species in the catalytic cycle, the Tyr233 phenol loses a proton to the nearby Asp179 residue, while at the same time, an electron is transferred to the iron to create an iron(III)-oxo active species. This charged tyrosyl residue directs the dipole moment along the C4-H bond of the substrate and guides the selectivity to the C4-hydroxylation of the substrate.

QM/MM Study Into the Mechanism of Oxidative C=C Double Bond Cleavage by Lignostilbene-α,β-Dioxygenase

The enzymatic biosynthesis of fragrance molecules from lignin fragments is an important reaction in biotechnology for the sustainable production of fine chemicals. In this work we investigated the biosynthesis of vanillin from lignostilbene by a nonheme iron dioxygenase using QM/MM and tested several suggested proposals via either an epoxide or dioxetane intermediate. Binding of dioxygen to the active site of the protein results in the formation of an iron(II)-superoxo species with lignostilbene cation radical. The dioxygenase mechanism starts with electrophilic attack of the terminal oxygen atom of the superoxo group on the central C=C bond of lignostilbene, and the second-coordination sphere effects in the substrate binding pocket guide the reaction towards dioxetane formation. The computed mechanism is rationalized with thermochemical cycles and valence bond schemes that explain the electron transfer processes during the reaction mechanism. Particularly, the polarity of the protein and the local electric field and dipole moments enable a facile electron transfer and an exergonic dioxetane formation pathway.

Bacteria at Work - Experimental and Theoretical Studies Reveal the Catalytic Mechanism of Ectoine Synthase

Ectoine synthase (EctC) catalyses the ultimate step of ectoine biosynthesis, a kosmotropic compound produced as compatible solute by many bacteria and some archaea or eukaryotes. EctC is an Fe2+-dependent homodimeric cytoplasmic protein. Using Mössbauer spectroscopy, molecular dynamics simulations and QM/MM calculations, we determined the most likely coordination number and geometry of the Fe2+ ion and proposed a mechanism of the EctC-catalysed reaction. Most notably, we show that apart from the three amino acids binding to the iron ion (Glu57, Tyr84 and His92), one water molecule and one hydroxide ion are required as additional ligands for the reaction to occur. They fill the first coordination sphere of the Fe2+-cofactor and act as critical proton donors and acceptors during the cyclization reaction.

Computational Study Into the Oxidative Ring-Closure Mechanism During the Biosynthesis of Deoxypodophyllotoxin

The nonheme iron dioxygenase deoxypodophyllotoxin synthase performs an oxidative ring-closure reaction as part of natural product synthesis in plants. How the enzyme enables the oxidative ring-closure reaction of (-)-yatein and avoids substrate hydroxylation remains unknown. To gain insight into the reaction mechanism and understand the details of the pathways leading to products and by-products we performed a comprehensive computational study. The work shows that substrate is bound tightly into the substrate binding pocket with the C7'-H bond closest to the iron(IV)-oxo species. The reaction proceeds through a radical mechanism starting with hydrogen atom abstraction from the C7'-H position followed by ring-closure and a final hydrogen transfer to form iron(II)-water and deoxypodophyllotoxin. Alternative mechanisms including substrate hydroxylation and an electron transfer pathway were explored but found to be higher in energy. The mechanism is guided by electrostatic perturbations of charged residues in the second-coordination sphere that prevent alternative pathways.

Catalytic divergencies in the mechanism of L-arginine hydroxylating nonheme iron enzymes

Many enzymes in nature utilize a free arginine (L-Arg) amino acid to initiate the biosynthesis of natural products. Examples include nitric oxide synthases, which generate NO from L-Arg for blood pressure control, and various arginine hydroxylases involved in antibiotic biosynthesis. Among the groups of arginine hydroxylases, several enzymes utilize a nonheme iron(II) active site and let L-Arg react with dioxygen and α-ketoglutarate to perform either C3-hydroxylation, C4-hydroxylation, C5-hydroxylation, or C4-C5-desaturation. How these seemingly similar enzymes can react with high specificity and selectivity to form different products remains unknown. Over the past few years, our groups have investigated the mechanisms of L-Arg-activating nonheme iron dioxygenases, including the viomycin biosynthesis enzyme VioC, the naphthyridinomycin biosynthesis enzyme NapI, and the streptothricin biosynthesis enzyme OrfP, using computational approaches and applied molecular dynamics, quantum mechanics on cluster models, and quantum mechanics/molecular mechanics (QM/MM) approaches. These studies not only highlight the differences in substrate and oxidant binding and positioning but also emphasize on electronic and electrostatic differences in the substrate-binding pockets of the enzymes. In particular, due to charge differences in the active site structures, there are changes in the local electric field and electric dipole moment orientations that either strengthen or weaken specific substrate C-H bonds. The local field effects, therefore, influence and guide reaction selectivity and specificity and give the enzymes their unique reactivity patterns. Computational work using either QM/MM or density functional theory (DFT) on cluster models can provide valuable insights into catalytic reaction mechanisms and produce accurate and reliable data that can be used to engineer proteins and synthetic catalysts to perform novel reaction pathways.

Identification of novel plant cysteine oxidase inhibitors from a yeast chemical genetic screen

Hypoxic responses in plants involve Plant Cysteine Oxidases (PCOs). They catalyze the N-terminal cysteine oxidation of Ethylene Response Factors VII (ERF-VII) in an oxygen-dependent manner, leading to their degradation via the cysteine N-degron pathway (Cys-NDP) in normoxia. In hypoxia, PCO activity drops, leading to the stabilization of ERF-VIIs and subsequent hypoxic gene upregulation. Thus far, no chemicals have been described to specifically inhibit PCO enzymes. In this work, we devised an in vivo pipeline to discover Cys-NDP effector molecules. Budding yeast expressing AtPCO4 and plant-based ERF-VII reporters was deployed to screen a library of natural-like chemical scaffolds and was further combined with an Arabidopsis Cys-NDP reporter line. This strategy allowed us to identify three PCO inhibitors, two of which were shown to affect PCO activity in vitro. Application of these molecules to Arabidopsis seedlings led to an increase in ERF-VII stability, induction of anaerobic gene expression, and improvement of tolerance to anoxia. By combining a high-throughput heterologous platform and the plant model Arabidopsis, our synthetic pipeline provides a versatile system to study how the Cys-NDP is modulated. Its first application here led to the discovery of at least two hypoxia-mimicking molecules with the potential to impact plant tolerance to low oxygen stress.

Melatonin Activation by Human Cytochrome P450 Enzymes: A Comparison between Different Isozymes

Cytochrome P450 enzymes in the human body play a pivotal role in both the biosynthesis and the degradation of the hormone melatonin. Melatonin plays a key role in circadian rhythms in the body, but its concentration is also linked to mood fluctuations as well as emotional well-being. In the present study, we present a computational analysis of the binding and activation of melatonin by various P450 isozymes that are known to yield different products and product distributions. In particular, the P450 isozymes 1A1, 1A2, and 1B1 generally react with melatonin to provide dominant aromatic hydroxylation at the C6-position, whereas the P450 2C19 isozyme mostly provides O-demethylation products. To gain insight into the origin of these product distributions of the P450 isozymes, we performed a comprehensive computational study of P450 2C19 isozymes and compared our work with previous studies on alternative isozymes. The work covers molecular mechanics, molecular dynamics and quantum mechanics approaches. Our work highlights major differences in the size and shape of the substrate binding pocket amongst the different P450 isozymes. Consequently, substrate binding and positioning in the active site varies substantially within the P450 isozymes. Thus, in P450 2C19, the substrate is oriented with its methoxy group pointing towards the heme, and therefore reacts favorably through hydrogen atom abstraction, leading to the production of O-demethylation products. On the other hand, the substrate-binding pockets in P450 1A1, 1A2, and 1B1 are tighter, direct the methoxy group away from the heme, and consequently activate an alternative site and lead to aromatic hydroxylation instead.

Computational Study of the Fe(II) and α-Ketoglutarate-Dependent Aryloxyalkanoate Dioxygenase (AAD-1) in the Degradation of the Herbicide 2,4-Dichlorophenoxyacetic Acid

The AAD-1 enzyme belongs to the Fe(II) and α-ketoglutarate (Fe/αKG)-dependent nonheme aryloxyalkanoate dioxygenase family (AADs), which catalyzes the breakdown of 2,4-dichlorophenoxyacetic acid (2,4-D, an active ingredient of thousands of commercial herbicides) by using the highly active Fe(IV)═O complex. Multiple species of bacteria degrade 2,4-D via a pathway initiated by AADs; however, the detail of how they promote the cleavage of the ether C-O bond to generate 2,4-dichlorophenol (2,4-DCP) and glyoxylate is still unclear, which is the prerequisite for the further degradation of these halogenated aromatics. In this work, based on the crystal structure of AAD-1, the computational models were constructed, and a series of QM/MM and QM-only calculations were performed to explore the cleavage of the ether bond in 2,4-D with the catalysis of AAD-1. Our calculations reveal that AAD-1 may be only responsible for the hydroxylation of the substrate to generate the intermediate hemiacetal, which corresponds to an overall energy barrier of 14.2 kcal/mol on the quintet state surface, and the decomposition of the hemiacetal in the active site center of AAD-1 was calculated to be rather slow, corresponding to an energy barrier of 24.5 kcal/mol. In contrast, the decomposition of the free hemiacetal molecule in a solvent was calculated to be quite easy. Whether the decomposition of the hemiacetal occurs inside or outside the activation site is still worthy of experimental verification.

Seeing the cis-Dihydroxylating Intermediate: A Mononuclear Nonheme Iron-Peroxo Complex in cis-Dihydroxylation Reactions Modeling Rieske Dioxygenases

The nature of reactive intermediates and the mechanism of the cis-dihydroxylation of arenes and olefins by Rieske dioxygenases and synthetic nonheme iron catalysts have been the topic of intense research over the past several decades. In this study, we report that a spectroscopically well characterized mononuclear nonheme iron(III)-peroxo complex reacts with olefins and naphthalene derivatives, yielding iron(III) cycloadducts that are isolated and characterized structurally and spectroscopically. Kinetics and product analysis reveal that the nonheme iron(III)-peroxo complex is a nucleophile that reacts with olefins and naphthalenes to yield cis-diol products. The present study reports the first example of the cis-dihydroxylation of substrates by a nonheme iron(III)-peroxo complex that yields cis-diol products.

Tailoring enzyme strategies and functional groups in biosynthetic pathways

Covering: 2000 to 2022Secondary metabolites are assembled by drawing off and committing some of the flux of primary metabolic building blocks to sets of enzymes that tailor the maturing scaffold to increase architectural and framework complexity, often balancing hydrophilic and hydrophobic surfaces. In this review we examine the small number of chemical strategies that tailoring enzymes employ in maturation of scaffolds. These strategies depend both on the organic functional groups present at each metabolic stage and one of two tailoring enzyme strategies. Nonoxidative tailoring enzymes typically transfer electrophilic fragments, acyl, alkyl and glycosyl groups, from a small set of thermodynamically activated but kinetically stable core metabolites. Oxidative tailoring enzymes can be oxygen-independent or oxygen-dependent. The oxygen independent oxidoreductases are often reversible nicotinamide-utilizing redox catalysts, flipping the nucleophilicity and electrophilicity of functional groups in pathway intermediates. O₂-dependent oxygenases, both mono- and dioxygenases, act by orthogonal, one electron strategies, generating carbon radical species. At sp³ substrate carbons, product alcohols may then behave as nucleophiles for subsequent waves of enzymatic tailoring. At sp² carbons in olefins, electrophilic epoxides have opposite reactivity and often function as "disappearing groups", opened by intramolecular nucleophiles during metabolite maturation. "Thwarted" oxygenases generate radical intermediates that rearrange internally and are not captured oxygenatively.

Second coordination sphere effects on the mechanistic pathways for dioxygen activation by a ferritin: involvement of a Tyr radical and the identification of a cation binding site

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp65, which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center.

The Apparently Unreactive Substrate Facilitates the Electron Transfer for Dioxygen Activation in Rieske Dioxygenases

Rieske dioxygenases belong to the non‐heme iron family of oxygenases and catalyze important cis‐dihydroxylation as well as O‐/N‐dealkylation and oxidative cyclization reactions for a wide range of substrates. The lack of substrate coordination at the non‐heme ferrous iron center, however, makes it particularly challenging to delineate the role of the substrate for productive O2 activation. Here, we studied the role of the substrate in the key elementary reaction leading to O2 activation from a theoretical perspective by systematically considering (i) the 6‐coordinate to 5‐coordinate conversion of the non‐heme Fe^II upon abstraction of a water ligand, (ii) binding of O2 , and (iii) transfer of an electron from the Rieske cluster. We systematically evaluated the spin‐state‐dependent reaction energies and structural effects at the active site for all combinations of the three elementary processes in the presence and absence of substrate using naphthalene dioxygenase as a prototypical Rieske dioxygenase. We find that reaction energies for the generation of a coordination vacancy at the non‐heme FeII center through thermoneutral H₂O reorientation and exothermic O2 binding prior to Rieske cluster oxidation are largely insensitive to the presence of naphthalene and do not lead to formation of any of the known reactive Fe‐oxygen species. By contrast, the role of the substrate becomes evident after Rieske cluster oxidation and exclusively for the 6‐coordinate non‐heme FeII sites in that the additional electron is found at the substrate instead of at the iron and oxygen atoms. Our results imply an allosteric control of the substrate on Rieske dioxygenase reactivity to happen prior to changes at the non‐heme FeII in agreement with a strategy that avoids unproductive O2 activation.

Pyrroloquinoline Quinone Chemistry, Biology, and Biosynthesis

The widely distributed, essential redox factor pyrroloquinoline quinone (PQQ, methoxatin) (1) was discovered in the mid-1960s. The breadth and depth of its biological effects are steadily being revealed, and understanding its biosynthesis at the genomic level is a continuing process. In this review, aspects of the chemistry, biology, biosynthesis, and commercial production of 1 at the gene level, and some applications, are presented from discovery through to mid-2021.

Local Charge Distributions, Electric Dipole Moments, and Local Electric Fields Influence Reactivity Patterns and Guide Regioselectivities in α-Ketoglutarate-Dependent Non-heme Iron Dioxygenases

Non-heme iron dioxygenases catalyze vital processes for human health related to the biosynthesis of essential products and the biodegradation of toxic metabolites. Often the natural product biosyntheses by these non-heme iron dioxygenases is highly regio- and chemoselective, which are commonly assigned to tight substrate-binding and positioning. However, recent high-level computational modeling has shown that substrate-binding and positioning is only part of the story and long-range electrostatic interactions can play a major additional role.In this Account, we review and summarize computational viewpoints on the high regio- and chemoselectivity of α-ketoglutarate-dependent non-heme iron dioxygenases and how external perturbations affect the catalysis. In particular, studies from our groups have shown that often a regioselectivity in enzymes can be accomplished by stabilization of the rate-determining transition state for the reaction through external charges, electric dipole moments, or local electric field effects. Furthermore, bond dissociation energies in molecules are shown to be influenced by an electric field effect, and through targeting a specific bond in an electric field, this can lead to an unusually specific reaction. For instance, in the carbon-induced starvation protein, we studied two substrate-bound conformations and showed that regardless of what C-H bond of the substrate is closest to the iron(IV)-oxo oxidant, the lowest hydrogen atom abstraction barrier is always for the pro-S C²-H abstraction due to an induced dipole moment of the protein that weakens this bond. In another example of the hygromycin biosynthesis enzyme, an oxidative ring-closure reaction in the substrate forms an ortho-δ-ester ring. Calculations on this enzyme show that the selectivity is guided by a protonated lysine residue in the active site that, through its positive charge, triggers a low energy hydrogen atom abstraction barrier. A final set of examples in this Account discuss the viomycin biosynthesis enzyme and the 2-(trimethylammonio)ethylphosphonate dioxygenase (TmpA) enzyme. Both of these enzymes are shown to possess a significant local dipole moment and local electric field effect due to charged residues surrounding the substrate and oxidant binding pockets. The protein dipole moment and local electric field strength changes the C-H bond strengths of the substrate as compared to the gas-phase triggers the regioselectivity of substrate activation. In particular, we show that in the gas phase and in a protein environment C-H bond strengths are different due to local electric dipole moments and electric field strengths. These examples show that enzymes have an intricately designed structure that enables a chemical reaction under ambient conditions through the positioning of positively and negatively charged residues that influence and enhance reaction mechanisms. These computational insights create huge possibilities in bioengineering to apply local electric field and dipole moments in proteins to achieve an unusual selectivity and specificity and trigger a fit-for-purpose biocatalyst for unique biotransformations.

Electrostatic perturbations in the substrate-binding pocket of taurine/α-ketoglutarate dioxygenase determine its selectivity

Taurine/α‐ketoglutarate dioxygenase is an important enzyme that takes part in the cysteine catabolism process in the human body and selectively hydroxylates taurine at the C¹‐position. Recent computational studies showed that in the gas‐phase the C²−H bond of taurine is substantially weaker than the C¹−H bond, yet no evidence exists of 2‐hydroxytaurine products. To this end, a detailed computational study on the selectivity patterns in TauD was performed. The calculations show that the second‐coordination sphere and the protonation states of residues play a major role in guiding the enzyme to the right selectivity. Specifically, a single proton on an active site histidine residue can change the regioselectivity of the reaction through its electrostatic perturbations in the active site and effectively changes the C¹−H and C²−H bond strengths of taurine. This is further emphasized by many polar and hydrogen bonding interactions of the protein cage in TauD with the substrate and the oxidant that weaken the pro‐R C¹−H bond and triggers a chemoselective reaction process. The large cluster models reproduce the experimental free energy of activation excellently.

Practical Enzymatic Production of Carbocycles

The Pd-catalyzed carbon-carbon bond formation pioneered by Heck in 1969 has dominated medicinal chemistry development for the ensuing fifty years. As the demand for more complex three-dimensional active pharmaceuticals continues to increase, preparative enzyme-mediated assembly, by virtue of its exquisite selectivity and sustainable nature, is poised to provide a practical and affordable alternative for accessing such compounds. In this minireview, we summarize recent state-of-the-art developments in practical enzyme-mediated assembly of carbocycles. When appropriate, background information on the enzymatic transformation is provided and challenges and/or limitations are also highlighted.

Structure and Functional Differences of Cysteine and 3-Mercaptopropionate Dioxygenases: A Computational Study

Thiol dioxygenases are important enzymes for human health; they are involved in the detoxification and catabolism of toxic thiol-containing natural products such as cysteine. As such, these enzymes have relevance to the development of Alzheimer's and Parkinson's diseases in the brain. Recent crystal structure coordinates of cysteine and 3-mercaptopropionate dioxygenase (CDO and MDO) showed major differences in the second-coordination spheres of the two enzymes. To understand the difference in activity between these two analogous enzymes, we created large, active-site cluster models. We show that CDO and MDO have different iron(III)-superoxo-bound structures due to differences in ligand coordination. Furthermore, our studies show that the differences in the second-coordination sphere and particularly the position of a positively charged Arg residue results in changes in substrate positioning, mobility and enzymatic turnover. Furthermore, the substrate scope of MDO is explored with cysteinate and 2-mercaptosuccinic acid and their reactivity is predicted.

Molecular insights into the endoperoxide formation by Fe(II)/α-KG-dependent oxygenase NvfI

Endoperoxide-containing natural products are a group of compounds with structurally unique cyclized peroxide moieties. Although numerous endoperoxide-containing compounds have been isolated, the biosynthesis of the endoperoxides remains unclear. NvfI from Aspergillus novofumigatus IBT 16806 is an endoperoxidase that catalyzes the formation of fumigatonoid A in the biosynthesis of novofumigatonin. Here, we describe our structural and functional analyses of NvfI. The structural elucidation and mutagenesis studies indicate that NvfI does not utilize a tyrosyl radical in the reaction, in contrast to other characterized endoperoxidases. Further, the crystallographic analysis reveals significant conformational changes of two loops upon substrate binding, which suggests a dynamic movement of active site during the catalytic cycle. As a result, NvfI installs three oxygen atoms onto a substrate in a single enzyme turnover. Based on these results, we propose a mechanism for the NvfI-catalyzed, unique endoperoxide formation reaction to produce fumigatonoid A.

Comparative Biology of Oxygen Sensing in Plants and Animals

Aerobic respiration is essential to almost all eukaryotes and sensing oxygen is a key determinant of survival. Analogous but mechanistically different oxygen-sensing pathways were adopted in plants and metazoan animals, and include ubiquitin-mediated degradation of transcription factors and direct sensing via non-heme iron(Fe²⁺)-dependent-dioxygenases. Key roles for oxygen sensing have been identified in both groups, with downstream signalling focussed on regulating gene transcription and chromatin modification to control development and stress responses. Components of sensing systems are promising targets for human therapeutic intervention and developing stress-resilient crops. Here, we review current knowledge about the origins, commonalities and differences between oxygen sensing in plants and animals.

Product Distributions of Cytochrome P450 OleTJE with Phenyl-Substituted Fatty Acids: A Computational Study

There are two types of cytochrome P450 enzymes in nature, namely, the monooxygenases and the peroxygenases. Both enzyme classes participate in substrate biodegradation or biosynthesis reactions in nature, but the P450 monooxygenases use dioxygen, while the peroxygenases take H₂O₂ in their catalytic cycle instead. By contrast to the P450 monooxygenases, the P450 peroxygenases do not require an external redox partner to deliver electrons during the catalytic cycle, and also no external proton source is needed. Therefore, they are fully self-sufficient, which affords them opportunities in biotechnological applications. One specific P450 peroxygenase, namely, P450 OleT_JE, reacts with long-chain linear fatty acids through oxidative decarboxylation to form hydrocarbons and, as such, has been implicated as a suitable source for the biosynthesis of biofuels. Unfortunately, the reactions were shown to produce a considerable amount of side products originating from C^α and C^β hydroxylation and desaturation. These product distributions were found to be strongly dependent on whether the substrate had substituents on the C^α and/or C^β atoms. To understand the bifurcation pathways of substrate activation by P450 OleT_JE leading to decarboxylation, C^α hydroxylation, C^β hydroxylation and C^α–C^β desaturation, we performed a computational study using 3-phenylpropionate and 2-phenylbutyrate as substrates. We set up large cluster models containing the heme, the substrate and the key features of the substrate binding pocket and calculated (using density functional theory) the pathways leading to the four possible products. This work predicts that the two substrates will react with different reaction rates due to accessibility differences of the substrates to the active oxidant, and, as a consequence, these two substrates will also generate different products. This work explains how the substrate binding pocket of P450 OleT_JE guides a reaction to a chemoselectivity.

Density Functional Theory Study into the Reaction Mechanism of Isonitrile Biosynthesis by the Nonheme Iron Enzyme ScoE

The nonheme iron enzyme ScoE catalyzes the biosynthesis of an isonitrile substituent in a peptide chain. To understand details of the reaction mechanism we created a large active site cluster model of 212 atoms that contains substrate, the active oxidant and the first- and second-coordination sphere of the protein and solvent. Several possible reaction mechanisms were tested and it is shown that isonitrile can only be formed through two consecutive catalytic cycles that both use one molecule of dioxygen and α-ketoglutarate. In both cycles the active species is an iron(IV)-oxo species that in the first reaction cycle reacts through two consecutive hydrogen atom abstraction steps: first from the N–H group and thereafter from the C–H group to desaturate the NH-CH2 bond. The alternative ordering of hydrogen atom abstraction steps was also tested but found to be higher in energy. Moreover, the electronic configurations along that pathway implicate an initial hydride transfer followed by proton transfer. We highlight an active site Lys residue that is shown to donate charge in the transition states and influences the relative barrier heights and bifurcation pathways. A second catalytic cycle of the reaction of iron(IV)-oxo with desaturated substrate starts with hydrogen atom abstraction followed by decarboxylation to give isonitrile directly. The catalytic cycle is completed with a proton transfer to iron(II)-hydroxo to generate the iron(II)-water resting state. The work is compared with experimental observation and previous computational studies on this system and put in a larger perspective of nonheme iron chemistry.

Electrostatic Perturbations from the Protein Affect C-H Bond Strengths of the Substrate and Enable Negative Catalysis in the TmpA Biosynthesis Enzyme

The nonheme iron dioxygenase 2-(trimethylammonio)-ethylphosphonate dioxygenase (TmpA) is an enzyme involved in the regio- and chemoselective hydroxylation at the C¹ -position of the substrate as part of the biosynthesis of glycine betaine in bacteria and carnitine in humans. To understand how the enzyme avoids breaking the weak C² -H bond in favor of C¹ -hydroxylation, we set up a cluster model of 242 atoms representing the first and second coordination sphere of the metal center and substrate binding pocket, and investigated possible reaction mechanisms of substrate activation by an iron(IV)-oxo species by density functional theory methods. In agreement with experimental product distributions, the calculations predict a favorable C¹ -hydroxylation pathway. The calculations show that the selectivity is guided through electrostatic perturbations inside the protein from charged residues, external electric fields and electric dipole moments. In particular, charged residues influence and perturb the homolytic bond strength of the C¹ -H and C² -H bonds of the substrate, and strongly strengthens the C² -H bond in the substrate-bound orientation.

Structure-guided insights into heterocyclic ring-cleavage catalysis of the non-heme Fe (II) dioxygenase NicX

Biodegradation of aromatic and heterocyclic compounds requires an oxidative ring cleavage enzymatic step. Extensive biochemical research has yielded mechanistic insights about catabolism of aromatic substrates; yet much less is known about the reaction mechanisms underlying the cleavage of heterocyclic compounds such as pyridine-ring-containing ones like 2,5-hydroxy-pyridine (DHP). 2,5-Dihydroxypyridine dioxygenase (NicX) from Pseudomonas putida KT2440 uses a mononuclear nonheme Fe(II) to catalyze the oxidative pyridine ring cleavage reaction by transforming DHP into N-formylmaleamic acid (NFM). Herein, we report a crystal structure for the resting form of NicX, as well as a complex structure wherein DHP and NFM are trapped in different subunits. The resting state structure displays an octahedral coordination for Fe(II) with two histidine residues (His²⁶⁵ and His³¹⁸), a serine residue (Ser³⁰²), a carboxylate ligand (Asp³²⁰), and two water molecules. DHP does not bind as a ligand to Fe(II), yet its interactions with Leu¹⁰⁴ and His¹⁰⁵ function to guide and stabilize the substrate to the appropriate position to initiate the reaction. Additionally, combined structural and computational analyses lend support to an apical dioxygen catalytic mechanism. Our study thus deepens understanding of non-heme Fe(II) dioxygenases.

What Determines the Selectivity of Arginine Dihydroxylation by the Nonheme Iron Enzyme OrfP?

The nonheme iron enzyme OrfP reacts with l-Arg selectively to form the 3R,4R-dihydroxyarginine product, which in mammals can inhibit the nitric oxide synthase enzymes involved in blood pressure control. To understand the mechanisms of dioxygen activation of l-Arg by OrfP and how it enables two sequential oxidation cycles on the same substrate, we performed a density functional theory study on a large active site cluster model. We show that substrate binding and positioning in the active site guides a highly selective reaction through C³ -H hydrogen atom abstraction. This happens despite the fact that the C³ -H and C⁴ -H bond strengths of l-Arg are very similar. Electronic differences in the two hydrogen atom abstraction pathways drive the reaction with an initial C³ -H activation to a low-energy ⁵ σ-pathway, while substrate positioning destabilizes the C⁴ -H abstraction and sends it over the higher-lying ⁵ π-pathway. We show that substrate and monohydroxylated products are strongly bound in the substrate binding pocket and hence product release is difficult and consequently its lifetime will be long enough to trigger a second oxygenation cycle.

Aerobic catabolism of sterols by microorganisms: key enzymes that open the 3-ketosteroid nucleus

Aerobic degradation of the sterol tetracyclic nucleus by microorganisms comprises the catabolism of A/B-rings, followed by that of C/D-rings. B-ring rupture at the C9,10-position is a key step involving 3-ketosteroid Δ1-dehydrogenase (KstD) and 3-ketosteroid 9α-hydroxylase (KstH). Their activities lead to the aromatization of C4,5-en-containing A-ring causing the rupture of B-ring. C4,5α-hydrogenated 3-ketosteroid could be produced by the growing microorganism containing a 5α-reductase. In this case, the microorganism synthesizes, in addition to KstD and KstH, a 3-ketosteroid Δ4-(5α)-dehydrogenase (Kst4D) in order to produce the A-ring aromatization, and consequently B-ring rupture. KstD and Kst4D are FAD-dependent oxidoreductases. KstH is composed of a reductase and a monooxygenase. This last component is the catalytic unit; it contains a Rieske-[2Fe-2S] center with a non-haem mononuclear iron in the active site. Published data regarding these enzymes are reviewed.

A novel concept for the biodegradation mechanism of dianionic catechol with homoprotocatechuate 2,3-dioxygenase: A non-proton-assisted process

The theory of "proton-assisted process" can well explain the catalytic mechanism of homoprotocatechuate 2,3-dioxygenase (2,3-HPCD) with a monoanionic substrate (homoprotocatechuate, HPCA). Here a "non-proton-assisted process" is presented to interpret catalytic mechanism of 2,3-HPCD with a dianionic substrate (4-nitrocatechol, 4NC). The ONIOM calculation is performed to investigate the reaction pathway of a wild-type 2,3-HPCD with 4NC (H200H-4NC system). The catalytic reaction is comprised of four steps: (1) A dioxygen attacks the aromatic ring to produce an alkylperoxo species. (2) O-O bond cleavage and the formation of an epoxide species occur. (3) A seven-membered O-heterocyclic compound is generated by the extinction of the epoxy structure. (4) The seven-membered ring undergoes ring opening to form the final product (C2-C3 cleavage product). The effective free energy barrier of the catalytic reaction of the H200H-4NC system is 26.2 kcal mol^-1, which is much higher than that of the H200H-HPCA system. Furthermore, two calculated electronic configurations (Fe(III)-O₂^•- and Fe(III)-SQ^•) have a high similarity to previously detected ones, which demonstrates that the Asn200 variant (H200N-4NC variant system) employs a C4 (para-carbon) pathway to produce a C4-C5 cleavage product. Our findings provide an in-depth understanding of the catalytic mechanisms of dianionic catechol and its derivatives.

C-H Amination via Nitrene Transfer Catalyzed by Mononuclear Non-Heme Iron-Dependent Enzymes

Expanding the reaction scope of natural metalloenzymes can provide new opportunities for biocatalysis. Mononuclear non-heme iron-dependent enzymes represent a large class of biological catalysts involved in the biosynthesis of natural products and catabolism of xenobiotics, among other processes. Here, we report that several members of this enzyme family, including Rieske dioxygenases as well as α-ketoglutarate-dependent dioxygenases and halogenases, are able to catalyze the intramolecular C-H amination of a sulfonyl azide substrate, thereby exhibiting a promiscuous nitrene transfer reactivity. One of these enzymes, naphthalene dioxygenase (NDO), was further engineered resulting in several active site variants that function as C-H aminases. Furthermore, this enzyme could be applied to execute this non-native transformation on a gram scale in a bioreactor, thus demonstrating its potential for synthetic applications. These studies highlight the functional versatility of non-heme iron-dependent enzymes and pave the way to their further investigation and development as promising biocatalysts for non-native metal-catalyzed transformations.

Selective Hydrogen Atom Abstraction from Dihydroflavonol by a Nonheme Iron Center Is the Key Step in the Enzymatic Flavonol Synthesis and Avoids Byproducts

The plant non-heme iron dioxygenase flavonol synthase performs a regioselective desaturation reaction as part of the biosynthesis of the signaling molecule flavonol that triggers the growing of leaves and flowers. These compounds also have health benefits for humans. Desaturation of aliphatic compounds generally proceeds through two consecutive hydrogen atom abstraction steps from two adjacent carbon atoms and in nature often is performed by a high-valent iron(IV)-oxo species. We show that the order of the hydrogen atom abstraction steps, however, is opposite of those expected from the C-H bond strengths in the substrate and determines the product distributions. As such, flavonol synthase follows a negative catalysis mechanism. Using density functional theory methods on large active-site model complexes, we investigated pathways for desaturation and hydroxylation by an iron(IV)-oxo active-site model. Contrary to thermochemical predictions, we find that the oxidant abstracts the hydrogen atom from the strong C²-H bond rather than the weaker C³-H bond of the substrate first. We analyze the origin of this unexpected selective hydrogen atom abstraction pathway and find that the alternative C³-H hydrogen atom abstraction would be followed by a low-energy and competitive substrate hydroxylation mechanism hence, should give considerable amount of byproducts. Our computational modeling studies show that substrate positioning in flavonol synthase is essential, as it guides the reactivity to a chemo- and regioselective substrate desaturation from the C²-H group, leading to desaturation products efficiently.

Expanding the roles for 2-oxoglutarate-dependent oxygenases in plant metabolism

Covering: up to 2018 2-Oxoglutarate-dependent oxygenases (2ODOs) comprise a large enzyme superfamily in plant genomes, second in size only to the cytochromes P450 monooxygenase (CYP) superfamily. 2ODOs participate in both primary and specialized plant pathways, and their occurrence across all life kingdoms points to an ancient origin. Phylogenetic evidence supports substantial expansion and diversification of 2ODOs following the split from the common ancestor of land plants. More conserved roles for these enzymes include oxidation within hormone metabolism, such as the recently described capacity of Dioxygenase for Auxin Oxidation (DAO) for governing auxin homeostasis. Conserved structural features among 2ODOs has provided a basis for continued investigation into their mechanisms, and recent structural work is expected to illuminate intriguing reactions such as that of 1-aminocyclopropane-1-carboxylic acid oxidase (ACCO). Phylogenetic radiation among this superfamily combined with neo- and subfunctionalization has enabled recruitment to highly specialized pathways, including those yielding medicines, flavours, dyes, poisons, and compounds important for plant-environment interactions. Catalytic versatility of 2ODOs in plants and across broader taxa continues to inspire biochemists tasked with the discovery of new enzymes. This highlight article summarizes recent reports up to 2018 of 2ODOs within plant metabolism. Furthermore, the respective contributions of 2ODOs and other oxidases to natural product biosynthesis are discussed as a framework for continued discovery.

Flavonol biosynthesis by nonheme iron dioxygenases: A computational study into the structure and mechanism

Plants produce flavonol compounds for vital functions regarding plant growth, fruit and flower colouring as well as fruit ripening processes. Several of these biosynthesis steps are stereo- and regioselective and are being carried out by nonheme iron enzymes. Using density functional theory calculations on a large active site model complex of flavanone-3β-hydroxylase (FHT), we established the mechanism for conversion of naringenin to its dihydroflavonol, which is a key step in the mechanism of flavonol biosynthesis. The reaction starts with dioxygen binding to the iron(II) centre and a reaction with α-ketoglutarate co-substrate gives succinate, an iron(IV)-oxo species and CO₂ with large exothermicity and small reaction barriers. The rate-determining reaction step in the mechanism; however, is hydrogen atom abstraction of an aliphatic CH bond by the iron(IV)-oxo species. We identify a large kinetic isotope effect for the replacement of the transferring hydrogen atom by deuterium. In a final step the OH and substrate radicals combine to form the alcohol product with a barrier of several kcal mol^-1. We show that the latter is the result of geometric constraints in the active site pocket. Furthermore, the calculations show that a weak tertiary CH bond is shielded from the iron(IV)-oxo species in the substrate binding position and therefore the enzyme is able to activate a stronger CH bond. As such, the flavanone-3β-hydroxylase enzyme reacts regioselectively with one specific CH bond of naringenin by avoiding activation of weaker bonds through tight substrate and oxidant positioning.

Effect of proximal ligand substitutions on the carbene and nitrene transferase activity of myoglobin

Engineered myoglobins were recently shown to be effective catalysts for abiological carbene and nitrene transfer reactions. Here, we investigated the impact of substituting the conserved heme-coordinating histidine residue with both proteinogenic (Cys, Ser, Tyr, Asp) and non-proteinogenic Lewis basic amino acids (3-(3'-pyridyl)-alanine, p-aminophenylalanine, and β-(3-thienyl)-alanine), on the reactivity of this metalloprotein toward these abiotic transformations. These studies showed that mutation of the proximal histidine residue with both natural and non-natural amino acids result in stable myoglobin variants that can function as both carbene and nitrene transferases. In addition, substitution of the proximal histidine with an aspartate residue led to a myoglobin-based catalyst capable of promoting stereoselective olefin cyclopropanation under nonreducing conditions. Overall, these studies demonstrate that proximal ligand substitution provides a promising strategy to tune the reactivity of myoglobin-based carbene and nitrene transfer catalysts and provide a first, proof-of-principle demonstration of the viability of pyridine-, thiophene-, and aniline-based unnatural amino acids for metalloprotein engineering.

Iron(ii) tetrafluoroborate complexes of new tetradentate C-scorpionates as catalysts for the oxidative cleavage of trans-stilbene with H2O2

Iron(ii) complexes of two new tetradentate C-scorpionate ligands are characterized. Both catalyze stilbene cleavage using either H₂O₂ or a O₂/photocatalyst oxidant.

Consensus Coexpression Network Analysis Identifies Key Regulators of Flower and Fruit Development in Wild Strawberry

The diploid strawberry, Fragaria vesca, is a developing model system for the economically important Rosaceae family. Strawberry fleshy fruit develops from the floral receptacle and its ripening is nonclimacteric. The external seed configuration of strawberry fruit facilitates the study of seed-to-fruit cross tissue communication, particularly phytohormone biosynthesis and transport. To investigate strawberry fruit development, we previously generated spatial and temporal transcriptome data profiling F. vesca flower and fruit development pre- and postfertilization. In this study, we combined 46 of our existing RNA-seq libraries to generate coexpression networks using the Weighted Gene Co-Expression Network Analysis package in R. We then applied a post-hoc consensus clustering approach and used bootstrapping to demonstrate consensus clustering's ability to produce robust and reproducible clusters. Further, we experimentally tested hypotheses based on the networks, including increased iron transport from the receptacle to the seed postfertilization and characterized a F. vesca floral mutant and its candidate gene. To increase their utility, the networks are presented in a web interface (www.fv.rosaceaefruits.org) for easy exploration and identification of coexpressed genes. Together, the work reported here illustrates ways to generate robust networks optimized for the mining of large transcriptome data sets, thereby providing a useful resource for hypothesis generation and experimental design in strawberry and related Rosaceae fruit crops.

Engineering plant architecture via CRISPR/Cas9-mediated alteration of strigolactone biosynthesis

BACKGROUND

Precision plant genome engineering holds much promise for targeted improvement of crop traits via unprecedented single-base level control over the genetic material. Strigolactones (SLs) are a key determinant of plant architecture, known for their role in inhibiting shoot branching (tillering).

RESULTS

We used CRISPR/Cas9 in rice (Oryza sativa) for targeted disruption of CAROTENOID CLEAVAGE DIOXYGENASE 7 (CCD7), which controls a key step in SL biosynthesis. The ccd7 mutants exhibited a striking increase in tillering, combined with a reduced height, which could be rescued by application of the synthetic SL analog GR24. Striga germination assays and liquid chromatography-mass spectrometry analysis showed that root exudates of ccd7 mutants were also SL deficient.

CONCLUSIONS

Taken together, our results show the potential and feasibility of the use of the CRISPR/Cas9 system for targeted engineering of plant architecture and for elucidating the molecular underpinnings of architecture-related traits.

Enhanced Electron-Transfer Reactivity of a Long-Lived Photoexcited State of a Cobalt-Oxygen Complex

Photodynamics and electron-transfer reactivity of an excited state derived from an earth-abundant mononuclear cobalt-oxygen complex ground state, [(TAML)Co^IV(O)]^2- (1; H₄TAML = 3,4,8,9-tetrahydro-3,3,6,6,9,9-hexamethyl-1 H-1,4,8,11-benzotetraazo-cyclotridecane-2,5,7,10-(6 H, 11 H)tetrone), prepared by electron-transfer oxidation of Li[(TAML)Co^III]·3(H₂O) (2) in a 1:1 acetonitrile/acetone solvent mixture at 5 °C, were investigated using a combination of femtosecond and nanosecond laser absorption spectroscopy. Visible light photoexcitation of 1 (λ_exc = 393 nm) resulted in generation of the excited state S₂* (lifetime: 1.4(4) ps), detected 2 ps after laser irradiation by femtosecond laser spectroscopy. The initially formed excited state S₂* converted to a lower-lying excited state, S₁* (λ_max = 580 nm), with rate constant k_c = 7(2) × 10¹¹ s^-1 (S₂* → S₁*). S₁* exhibited a 0.6(1) ns lifetime and converted to the initial ground state 1 with rate constant k_d = 1.7(3) × 10⁹ s^-1 (S₁* → 1). The same excited state dynamics was observed when 1 was generated by electron-transfer oxidation of 2 using different one-electron oxidants such as Cu(OTf)₂ (OTf^- = triflate anion), [Fe(bpy)₃]³⁺ (bpy = 2,2'-bipyridine), and tris(4-bromophenyl)ammoniumyl radical cation (TBPA^•+). The electron-transfer reactivity of S₁* was probed by nanosecond laser photoexcitation of 1 in the presence of a series of electron donors with different one-electron oxidation potentials ( E_ox vs SCE): benzene (2.35 V), toluene (2.20 V), m-xylene (2.02 V), and anisole (1.67 V). The excited state S₁* engaged in electron-transfer reactions with m-xylene and anisole to generate π-dimer radical cations of m-xylene and anisole, respectively, observed by nanosecond laser transient absorption spectroscopy, whereas no reactivity was observed toward benzene and toluene. Such differential electron-transfer reactivity depending on the E_ox values of electron donors allowed the estimation of the one-electron reduction potential of S₁* ( E_red*) as 2.1(1) V vs SCE, which is much higher than that of the ground state ( E_red = 0.86 V vs SCE).

Structure and function of a lignostilbene-α,β-dioxygenase orthologue from Pseudomonas brassicacearum

BACKGROUND

Stilbene cleaving oxygenases (SCOs), also known as lignostilbene-α,β-dioxygenases (LSDs) mediate the oxidative cleavage of the olefinic double bonds of lignin-derived intermediate phenolic stilbenes, yielding small modified benzaldehyde compounds. SCOs represent one branch of the larger carotenoid cleavage oxygenases family. Here, we describe the structural and functional characterization of an SCO-like enzyme from the soil-born, bio-control agent Pseudomonas brassicacearum.

METHODS

In vitro and in vivo assays relying on visual inspection, spectrophotometric quantification, as well as liquid-chormatographic and mass spectrometric characterization were applied for functional evaluation of the enzyme. X-ray crystallographic analyses and in silico modeling were applied for structural investigations.

RESULTS

In vitro assays demonstrated preferential cleavage of resveratrol, while in vivo analyses detected putative cleavage of the straight chain carotenoid, lycopene. A high-resolution structure containing the seven-bladed β-propeller fold and conserved 4-His-Fe unit at the catalytic site, was obtained. Comparative structural alignments, as well as in silico modelling and docking, highlight potential molecular factors contributing to both the primary in vitro activity against resveratrol, as well as the putative subsidiary activities against carotenoids in vivo, for future validation.

CONCLUSIONS

The findings reported here provide validation of the SCO structure, and highlight enigmatic points with respect to the potential effect of the enzyme's molecular environment on substrate specificities for future investigation.

Directed Evolution of Protein Catalysts

Directed evolution is a powerful technique for generating tailor-made enzymes for a wide range of biocatalytic applications. Following the principles of natural evolution, iterative cycles of mutagenesis and screening or selection are applied to modify protein properties, enhance catalytic activities, or develop completely new protein catalysts for non-natural chemical transformations. This review briefly surveys the experimental methods used to generate genetic diversity and screen or select for improved enzyme variants. Emphasis is placed on a key challenge, namely how to generate novel catalytic activities that expand the scope of natural reactions. Two particularly effective strategies, exploiting catalytic promiscuity and rational design, are illustrated by representative examples of successfully evolved enzymes. Opportunities for extending these approaches to more complex biocatalytic systems are also considered.

A global perspective on carotenoids: Metabolism, biotechnology, and benefits for nutrition and health

Carotenoids are lipophilic isoprenoid compounds synthesized by all photosynthetic organisms and some non-photosynthetic prokaryotes and fungi. With some notable exceptions, animals (including humans) do not produce carotenoids de novo but take them in their diets. In photosynthetic systems carotenoids are essential for photoprotection against excess light and contribute to light harvesting, but perhaps they are best known for their properties as natural pigments in the yellow to red range. Carotenoids can be associated to fatty acids, sugars, proteins, or other compounds that can change their physical and chemical properties and influence their biological roles. Furthermore, oxidative cleavage of carotenoids produces smaller molecules such as apocarotenoids, some of which are important pigments and volatile (aroma) compounds. Enzymatic breakage of carotenoids can also produce biologically active molecules in both plants (hormones, retrograde signals) and animals (retinoids). Both carotenoids and their enzymatic cleavage products are associated with other processes positively impacting human health. Carotenoids are widely used in the industry as food ingredients, feed additives, and supplements. This review, contributed by scientists of complementary disciplines related to carotenoid research, covers recent advances and provides a perspective on future directions on the subjects of carotenoid metabolism, biotechnology, and nutritional and health benefits.

Group VII Ethylene Response Factors in Arabidopsis: Regulation and Physiological Roles

The role of ERF-VII TFs in higher plants is to coordinate their signature response to oxygen deficiency, but additional layers of modulation of ERF-VII activity enrich their regulatory range.