Mechanistic Insights into the Decoupled Desaturation and Epoxidation Catalyzed by Dioxygenase AsqJ Involved in the Biosynthesis of Quinolone Alkaloids

Rational design improves both thermostability and activity of a new D-tagatose 3-epimerase from Kroppenstedtia eburnean to produce D-allulose

D-allulose is a naturally occurring rare sugar and beneficial to human health. However, the efficient biosynthesis of D-allulose remains a challenge. Here, we mined a new D-tagatose 3-epimerase from Kroppenstedtia eburnean (KeDt3e) with high catalytic efficiency. Initially, crucial factors contributing to the low conversion of KeDt3e were identified through crystal structure analysis, density functional theory calculations (DFT), and molecular dynamics (MD) simulations. Subsequently, based on the mechanism, combining restructuring the flexible region, proline substitution based onconsensus sequence analysis, introducing disulfide bonds, and grafting properties, and reshaping the active center, the optimal mutant M5 of KeDt3e was obtained with enhanced thermostability and activity. The optimal mutant M5 exhibited an enzyme activity of 130.8 U/mg, representing a 1.2-fold increase; Tm value increased from 52.7 °C to 71.2 °C; and half-life at 55 °C extended to 273.7 min, representing a 58.2-fold improvement, and the detailed mechanism of performance improvement was analyzed. Finally, by screening the ribosome-binding site (RBS) of the optimal mutant M5 recombinant bacterium (G01), the engineered strain G08 with higher expression levels was obtained. The engineered strain G08 catalyzed 500 g/L D-fructose to produce 172.4 g/L D-allulose, with a conversion of 34.4% in 0.5 h and productivity of 344.8 g/L/h on a 1 L scale. This study presents a promising approach for industrial-scale production of D-allulose.

Unraveling the Mechanism of the Oxidative C-C Bond Coupling Reaction Catalyzed by Deoxypodophyllotoxin Synthase

Deoxypodophyllotoxin synthase (DPS), a nonheme Fe(II)/2-oxoglutarate (2OG)-dependent oxygenase, is a key enzyme that is involved in the construction of the fused-ring system in (-)-podophyllotoxin biosynthesis by catalyzing the C-C coupling reaction. However, the mechanistic details of DPS-catalyzed ring formation remain unclear. Herein, our quantum mechanics/molecular mechanics (QM/MM) calculations reveal a novel mechanism that involves the recycling of CO2 (a product of decarboxylation of 2OG) to prevent the formation of hydroxylated byproducts. Our results show that CO2 can react with the FeIII-OH species to generate an unusual FeIII-bicarbonate species. In this way, hydroxylation is avoided by consuming the OH group. Then, the C-C coupling followed by desaturation yields the final product, deoxypodophyllotoxin. This work highlights the crucial role of the CO2 molecule, generated in the crevice between the iron active site and the substrate, in controlling the reaction selectivity.

Theoretical investigation of the reaction mechanism of THP oxidative rearrangement catalysed by BBOX

γ-Butyrobetaine hydroxylase (BBOX) is a non-heme FeII/2OG dependent enzyme that is able to perform two different kinds of catalytic reactions on 3-(2,2,2-trimethylhydrazinium) propionate (THP) to produce distinct catalytic products. Although the structure of BBOX complexed with THP has been resolved, the details of its catalytic mechanism are still elusive. In this study, by employing molecular dynamics (MD) simulations and density functional theory (DFT) calculations, the mechanism of the THP oxidative rearrangement reactions catalysed by BBOX was investigated. Our calculations revealed how the enzyme undergoes a conformational conversion to initiate the catalytic reactions. In the first catalytic step, BBOX performs hydrogen abstraction from the substrate THP as a common non-heme iron enzyme. Due to the structure of the substrate stabilizing the radical species and polarizing the adjacent N-N bond, in the next step, THP takes the pathway for N-N bond homolysis but not regular hydroxyl rebounding. The cleaved ammonium radical could either react with the hydroxyl group on the iron centre of the enzyme or recombine with the other cleaved fragment of the substrate to generate the rearranged product. This study revealed the catalytic mechanism of BBOX, detailing how the enzyme and the substrate regulated the hydroxyl rebound process to generate various products.

Structure and Ligand Based Virtual Screening and MPO Topological Analysis of Triazolo Thiadiazepine-fused Coumarin Derivatives as Anti-Parkinson Drug Candidates

Parkinson's disease (PD) is a debilitating condition that can cause locomotor problems in affected patients, such as tremors and body rigidity. PD therapy often includes the use of monoamine oxidase B (MAOB) inhibitors, particularly phenylhalogen compounds and coumarin-based semi-synthetic compounds. The objective of this study was to analyze the structural, pharmacokinetic, and pharmacodynamic profile of a series of Triazolo Thiadiazepine-fused Coumarin Derivatives (TDCDs) against MAOB, in comparison with the inhibitor safinamide. To achieve this goal, we utilized structure-based virtual screening techniques, including target prediction and absorption, distribution, metabolism, and excretion (ADME) prediction based on multi-parameter optimization (MPO) topological analysis, as well as ligand-based virtual screening techniques, such as docking and molecular dynamics. The findings indicate that the TDCDs exhibit structural similarity to other bioactive compounds containing coumarin and MAOB-binding azoles, which are present in the ChEMBL database. The topological analyses suggest that TDCD3 has the best ADME profile, particularly due to the alignment between low lipophilicity and high polarity. The coumarin and triazole portions make a strong contribution to this profile, resulting in a permeability with Papp estimated at 2.15 × 10-5 cm/s, indicating high cell viability. The substance is predicted to be metabolically stable. It is important to note that this is an objective evaluation based on the available data. Molecular docking simulations showed that the ligand has an affinity energy of - 8.075 kcal/mol with MAOB and interacts with biological substrate residues such as Pro102 and Phe103. The results suggest that the compound has a safe profile in relation to the MAOB model, making it a promising active ingredient for the treatment of PD.

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.

In Silico Study in MPO and Molecular Docking of the Synthetic Drynaran Analogues Against the Chronic Tinnitus: Modulation of the M1 Muscarinic Acetylcholine Receptor

Tinnitus is a syndrome that affects the human auditory system and is characterized by a perception of sounds in the absence of acoustic stimuli, or in total silence. Research indicates that muscarinic acetylcholine receptors (mAChRs), especially the M1 type, have a fundamental role in the alterations of auditory perceptions of tinnitus. Here, a series of computer-aided tools were used, from molecular surface analysis software to services available on the web for estimating pharmacokinetics and pharmacodynamics. The results infer that the low lipophilicity ligands, that is, the 1a-d alkyl furans, present the best pharmacokinetic profile, as compounds with an optimal alignment between permeability and clearance. However, only ligands 1a and 1b have properties that are safe for the central nervous system, the site of cholinergic modulation. These ligands showed similarity with compounds deposited in the European Molecular Biology Laboratory chemical (ChEMBL) database acting on the mAChRs M1 type, the target selected for the molecular docking test. The simulations suggest that the 1 g ligand can form the ligand-receptor complex with the best affinity energy order and that, together with the 1b ligand, they are competitive agonists in relation to the antagonist Tiotropium, in addition to acting in synergism with the drug Bromazepam in the treatment of chronic tinnitus.

QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis

Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.

Dioxygen Activation and Nδ,Nε-Dihydroxylation Mechanism Involved in the Formation of N-Nitrosourea Pharmacophore in Streptozotocin Catalyzed by Nonheme Diiron Enzyme SznF

SznF is a nonheme diiron-dependent enzyme that catalyzes the critical N-nitrosation involved in the formation of the N-nitrosourea moiety in the pancreatic cancer drug streptozotocin. The N-nitrosation contains two successive N-hydroxylation and N-nitrosation steps, which are carried out by two separate active sites, namely, the central domain and cupin domain. Recently, the crystal structure of SznF was obtained, and the central domain was proved to contain a diiron cofactor to catalyze the N-hydroxylation. In this work, to gain insights into the O₂ activation and the successive N-hydroxylation mechanism, on the basis of the high-resolution crystal structure, the enzyme-substrate complex models were constructed, and a series of combined QM/MM calculations were performed. Based on our calculations, the activation of O₂ starts from the diiron(II,III)-superoxo (S) to generate the diiron(IV)-oxo species (Q) via a diiron(III,III)-peroxo (P)-like transition state or unstable intermediate (P'), and species P' can be described as a hybridization of diiron(IV)-oxo species and diiron(III,III)-peroxo (P) owing to the long distances of Fe1-Fe2 (4.22 Å) and O1-O2 (1.89 Å), which is different from those of other nonheme diiron enzymes. In the following hydroxylation of N^δ and N^ε, the N^δ-hydroxylation was confirmed to occur first, agreeing with the experimental observations. Because the diiron(IV)-oxo species (Q) is responsible for hydroxylation, the reaction follows the H-abstraction/OH rebound mechanism, and the first abstraction occurs on the N^δ-H rather than N^ε-H, which may be attributed to the different orientation of Fe(IV)-oxo relative to N-H as well as the bond dissociation enthalpies of two N-H bonds. The hydroxylation of N-methyl-L-arginine does not employ the diiron(III,III)-hydroperoxo (P″) to trigger the electrophilic attack of the guanidine to directly form the N-O bond, as previously suggested. In addition, our calculations also revealed that the direct attack of the Fe(IV)═O unit to the N^δ of the substrate corresponds to a higher barrier than that in the H-abstraction/OH rebound mechanism. These results may provide useful information for understanding the formation of the di-hydroxylation intermediate involved in the biosynthesis of N-nitrosation.

Peroxy Intermediate Drives Carbon Bond Activation in the Dioxygenase AsqJ

Dioxygenases catalyze stereoselective oxygen atom transfer in metabolic pathways of biological, industrial, and pharmaceutical importance, but their precise chemical principles remain controversial. The α-ketoglutarate (αKG)-dependent dioxygenase AsqJ synthesizes biomedically active quinolone alkaloids via desaturation and subsequent epoxidation of a carbon-carbon bond in the cyclopeptin substrate. Here, we combine high-resolution X-ray crystallography with enzyme engineering, quantum-classical (QM/MM) simulations, and biochemical assays to describe a peroxidic intermediate that bridges the substrate and active site metal ion in AsqJ. Homolytic cleavage of this moiety during substrate epoxidation generates an activated high-valent ferryl (Fe^IV = O) species that mediates the next catalytic cycle, possibly without the consumption of the metabolically valuable αKG cosubstrate. Our combined findings provide an important understanding of chemical bond activation principles in complex enzymatic reaction networks and molecular mechanisms of dioxygenases.

How Human TET2 Enzyme Catalyzes the Oxidation of Unnatural Cytosine Modifications in Double-Stranded DNA

Methylation of cytosine bases is strongly linked to gene expression, imprinting, aging, and carcinogenesis. The Ten-eleven translocation (TET) family of enzymes, which are Fe(II)/2-oxoglutarate (2OG)-dependent enzymes, employ Fe(IV)=O species to dealkylate the lesioned bases to an unmodified cytosine. Recently, it has been shown that the TET2 enzyme can catalyze promiscuously DNA substrates containing unnatural alkylated cytosine. Such unnatural substrates of TET can be used as direct probes for measuring the TET activity or capturing TET from cellular samples. Herein, we studied the catalytic mechanisms during the oxidation of the unnatural C5-position modifications (5-ethylcytosine (5eC), 5-vinylcytosine (5vC) and 5-ethynylcytosine (5eyC)) and the demethylation of N4-methylated lesions (4-methylcytosine (4mC) and 4,4-dimethylcytosine(4dmC)) of the cytosine base by the TET2 enzyme using molecular dynamics (MD) and combined quantum mechanics and molecular mechanics (QM/MM) computational approaches. The results reveal that the chemical nature of the alkylation of the double-stranded (ds) DNA substrates induces distinct changes in the interactions in the binding site, the second coordination sphere, and long-range correlated motions of the ES complexes. The rate-determining hydrogen atom transfer (HAT) is faster in N4-methyl substituent substrates than in the C5-alkylations. Importantly, the calculations show the preference of hydroxylation over desaturation in both 5eC and 5vC substrates. The studies elucidate the post-hydroxylation rearrangements of the hydroxylated intermediates of 5eyC and 5vC to ketene and 5-formylmethylcytosine (5fmC), respectively, and hydrolysis of hemiaminal intermediate of 4mC to formaldehyde and unmodified cytosine proceed exclusively in aqueous solution outside of the enzyme environment. Overall, the studies show that the chemical nature of the unnatural alkylated cytosine substrates exercises distinct effects on the binding interactions, reaction mechanism, and dynamics of TET2.

Critical Roles of Exchange and Superexchange Interactions in Dictating Electron Transfer and Reactivity in Metalloenzymes

Electron transfer (ET) is a fundamental process in transition-metal-dependent metalloenzymes. In these enzymes, the spin-spin interactions within the same metal center and/or between different metal sites can play a pivotal role in the catalytic cycle and reactivity. This Perspective highlights that the exchange and/or superexchange interactions can intrinsically modulate the inner-sphere and long-range electron transfer, thus controlling the mechanism and activity of metalloenzymes. For mixed-valence diiron oxygenases, the spin-regulated inner-sphere ET can be dictated by exchange interactions, leading to efficient O-O bond activations. Likewise, the spin-regulated inner-sphere ET can be enhanced by both exchange and superexchange interactions in [Fe4S4]-dependent SAM enzymes, which enable the efficient cleavage of the S─C(γ) or S─C5' bond of SAM. In addition to inner-sphere ET, superexchange interactions may modulate the long-range ET between metalloenzymes. We anticipate that the exchange and superexchange enhanced reactivity could be applicable in other important metalloenzymes, such as Photosystem II and nitrogenases.

Computational Study of Aromatic Hydroxylation Catalyzed by the Iron-Dependent Hydroxylase PqqB Involved in the Biosynthesis of Redox Cofactor Pyrroloquinoline Quinone

PqqB from Methylobacterium extorquens is a unique nonheme iron-dependent hydroxylase involved in the biosynthesis of redox cofactor pyrroloquinoline quinone (PQQ). A series of recent experiments have demonstrated that PqqB catalyzes the stepwise insertions of two oxygen atoms into the tyrosine ring of the diamino acid substrate, generating the quinone moiety of PQQ; however, the reaction details have not been elucidated yet. In this paper, on the basis of the crystal structures, the enzyme-substrate complex models were constructed, and the catalytic mechanism of PqqB was explored by performing a series of combined QM/MM calculations. Our results confirmed that the first hydroxylation is performed by the highly reactive Fe^IV-oxo species and follows the typical H-abstraction/hydroxyl rebound mechanism, which is similar to the common aliphatic hydroxylation catalyzed by the α-KG enzymes. Nevertheless, the second hydroxylation is achieved by the Fe-O₂ species, and the reactant complex can be described as an intermediate radical-Fe^II-superoxide, that is, the dioxygen is activated by accepting an electron from the bidentate coordination intermediate. Since both the dioxygen and intermediate are activated by electron transfer, the distal oxygen of superoxide can directly attack the carbonyl carbon of substrate to form an alkylperoxo intermediate, then the O-O heterolysis occurs to afford the epoxide intermediate, which finally evolves into the product by rearrangement. It is the bidentate coordination of catechol moiety to iron that leads to the one-electron oxidation of the substrate by the dioxygen, which significantly activates the substrate and promotes the superoxide radical attack. During the catalysis, Asp90 and His240 in the second sphere play an important role by acting as acid-base catalysts to mediate the proton transfer and manipulate the suitable orientation of superoxide. These findings may provide useful information for understanding the unique reaction mechanism of PqqB that employs both the Fe^IV-oxo and Fe^II-superoxide to carry out the aromatic hydroxylation.

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.

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.

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.

Mechanistic Insights into the Oxidative Ring Expansion from Penicillin N to Deacetoxycephalosporin C Catalyzed by a Nonheme Iron(II) and α-KG-Dependent Oxygenase

Deacetoxycephalosporin C synthase (DAOCS) is a nonheme iron(II) and 2-oxoglutarate (α-KG)-dependent oxygenase that catalyzes the oxidative ring expansion of penicillin N (penN) to deacetoxycephalosporin C (DAOC). Earlier reported crystal structures of DAOCS indicated that the substrate penicillin binds at the same site of succinate, leading to the proposal of the unusual "ping-pong" mechanism. However, more recent data provided evidence of the formation of ternary DAOCS·α-KG·penN complex, and thus DAOCS should follow the usual consensus mechanism of α-KG-dependent nonheme iron(II) oxygenases. Nevertheless, how DAOCS catalyzes the ring expansion is unknown. In this paper, on the basis of the crystal structure, we constructed two reactant models and performed a series of combined quantum mechanics/molecular mechanics (QM/MM) calculations to illuminate the catalysis of DAOCS. The binding mode of substrate was found to be crucial in determining which hydrogen atom in two methyl groups is first abstracted and whether the second H-abstraction to be abstracted in the final desaturation step locates in a suitable orientation. The highly reactive Fe^IV-oxo species prefers to abstract a hydrogen atom from one of two methyl groups in penN to trigger the ring arrangement. After the H-abstraction, the generated methylene radical intermediate can easily initiate the ring arrangement. First, the C-S bond cleaves to generate a thiyl radical, which is in concert with the formation of the terminal C═C double bond; the newly generated thiyl radical then rapidly shifts to the more stable tertiary C atom to complete ring expansion. In the final step, the Fe^III-OH species abstracts the second hydrogen to give the desaturated DAOC product. During the catalysis, no active site residue is directly involved in the chemistry, which implies that the other pocket residues except the coordinate ones with iron play a role only in anchoring the substrate.

Epoxidation Catalyzed by the Nonheme Iron(II)- and 2-Oxoglutarate-Dependent Oxygenase, AsqJ: Mechanistic Elucidation of Oxygen Atom Transfer by a Ferryl Intermediate

Mechanisms of enzymatic epoxidation via oxygen atom transfer (OAT) to an olefin moiety is mainly derived from the studies on thiolate-heme containing epoxidases, such as cytochrome P450 epoxidases. The molecular basis of epoxidation catalyzed by nonheme-iron enzymes is much less explored. Herein, we present a detailed study on epoxidation catalyzed by the nonheme iron(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenase, AsqJ. The native substrate and analogues with different para substituents ranging from electron-donating groups (e.g., methoxy) to electron-withdrawing groups (e.g., trifluoromethyl) were used to probe the mechanism. The results derived from transient-state enzyme kinetics, Mössbauer spectroscopy, reaction product analysis, X-ray crystallography, density functional theory calculations, and molecular dynamic simulations collectively revealed the following mechanistic insights: (1) The rapid O2 addition to the AsqJ Fe(II) center occurs with the iron-bound 2OG adopting an online-binding mode in which the C1 carboxylate group of 2OG is trans to the proximal histidine (His134) of the 2-His-1-carboxylate facial triad, instead of assuming the offline-binding mode with the C1 carboxylate group trans to the distal histidine (His211); (2) The decay rate constant of the ferryl intermediate is not strongly affected by the nature of the para substituents of the substrate during the OAT step, a reactivity behavior that is drastically different from nonheme Fe(IV)-oxo synthetic model complexes; (3) The OAT step most likely proceeds through a stepwise process with the initial formation of a C(benzylic)-O bond to generate an Fe-alkoxide species, which is observed in the AsqJ crystal structure. The subsequent C3-O bond formation completes the epoxide installation.

Mechanism of Uncoupled Carbocyclization and Epimerization Catalyzed by Two Non-Heme Iron/α-Ketoglutarate Dependent Enzymes

The non-heme iron/α-ketoglutarate dependent enzymes SnoK and SnoN from Streptomyces nogalater are involved in the biosynthesis of anthracycline nogalamycin. Although they have similar active sites, SnoK is responsible for carbocyclization whereas SnoN solely catalyzes the hydroxyl epimerization. Herein, we performed docking, molecular simulations, and a series of combined quantum mechanics and molecular mechanics (QM/MM) calculations to illuminate the mechanisms of two enzymes. The catalytic reactions of two enzymes occur on the quintet state surface. For SnoK, the whole reaction includes two separated hydrogen-abstraction steps and one radical addition, and the latter step is calculated to be rate limiting with an energy barrier of 21.7 kcal/mol. Residue D106 is confirmed to participate in the construction of the hydrogen bond network, which plays a crucial role in positioning the bulky substrate in a specific orientation. Moreover, it is found that SnoN is only responsible for the hydrogen abstraction of the intermediate, and no residue was suggested to be suitable for donating a hydrogen atom to the substrate radical, which further confirms the suggestion based on experiments that either a cellular reductant or another enzyme protein could donate a hydrogen atom to the substrate. Our docking results coincide with the previous structural study that the different roles of two enzymes are achieved by minor changes in the alignment of the substrates in front of the reactive ferryl-oxo species. This work highlights the reaction mechanisms catalyzed by SnoK and SnoN, which is helpful for engineering the enzymes for the biosynthesis of anthracycline nogalamycin.

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.

Mechanism of fatty acid decarboxylation catalyzed by a non-heme iron oxidase (UndA): a QM/MM study

UndA is a non-heme iron enzyme that was recognized to catalyze the decarboxylation of medium chain (C10-C14) fatty acids to produce trace amounts of 1-alkenes. Owing to the electron imbalance during the oxidative decarboxylation of the substrate and the reduction of O₂, only single turnover reactions were obtained in UndA in vitro assays. Unlike the general non-heme iron enzymes, the catalytic efficiency of UndA is quite low. According to the previous proposal, both Fe^III-OO˙^- and Fe^IV[double bond, length as m-dash]O complexes may abstract the β-H of fatty acids to trigger the oxidative decarboxylation reaction. Herein, on the basis of the crystal structures of UndA in complex with the substrate analogues, we constructed a series of computational models and performed quantum mechanics/molecular mechanics (QM/MM) calculations to explore the UndA-catalyzed decarboxylation using lauric acid as the substrate. Our calculation results reveal that only the Fe^III-OO˙^- complex can initiate the decarboxylation, and the substrate (lauric acid) should monodentately coordinate to the Fe center to facilitate the β-H abstraction. In addition, the monodentate coordination corresponds to higher relative energy than the bidentate mode, which may explain the low efficiency of UndA. It is also revealed that as long as the β-H is extracted by the Fe^III-OO˙^-, the decarboxylation of the substrate radical is quite easy, and an electron transfer from the substrate to the iron center is the prerequisite. For the Fe^IV[double bond, length as m-dash]O complex, since the β-H is far from the O_Fe atom and the angle of ∠Fe-O-H is 53.1°, the H-abstraction is calculated to be difficult.

Insights into the Mechanism and Enantioselectivity in the Biosynthesis of Ergot Alkaloid Cycloclavine Catalyzed by Aj_EasH from Aspergillus japonicus

Cycloclavine is a complex ergot alkaloid containing an unusual cyclopropyl moiety, which has a wide range of biological activities and pharmaceutical applications. The biosynthesis of cycloclavine requires a series of enzymes, one of which is a nonheme Fe^II/α-ketoglutarate-dependent (aKG) oxidase (Aj_EasH). According to the previous proposal, the cyclopropyl ring formation catalyzed by Aj_EasH follows an unprecedented oxidative mechanism; however, the reaction details are unknown. In this article, on the basis of the recently obtained crystal structure of Aj_EasH (EasH from Aspergillus japonicas), the reactant models were built, and the reaction details were investigated by performing QM-only and combined QM and MM calculations. Our calculation results reveal that the biosynthesis of cyclopropyl moiety involves a radical intermediate rather than a carbocationic or carbanionic intermediate as in the biosynthesis of terpenoid family. The iron(IV)-oxo first abstracts a hydrogen atom from the substrate to trigger the reaction, and then the generated radical intermediate undergoes ring rearrangement to form the fused 5-3 ring system of cycloclavine. On the basis of our calculations, the absolute configuration of the cycloclavine catalyzed by Aj_EasH from Aspergillus japonicus should be (5R,8R,10R), which is different from the product isolated from Ipomoea hildebrandtii (5R,8S,10S). Residues at the active site play an important role in substrate binding, ring rearrangement, and enantioselectivity.

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.

Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment

ChemShell is a scriptable computational chemistry environment with an emphasis on multiscale simulation of complex systems using combined quantum mechanical and molecular mechanical (QM/MM) methods. Motivated by a scientific need to efficiently and accurately model chemical reactions on surfaces and within microporous solids on massively parallel computing systems, we present a major redevelopment of the ChemShell code, which provides a modern platform for advanced QM/MM embedding models. The new version of ChemShell has been re-engineered from the ground up with a new QM/MM driver module, an improved parallelization framework, new interfaces to high performance QM and MM programs, and a user interface written in the Python programming language. The redeveloped package is capable of performing QM/MM calculations on systems of significantly increased size, which we illustrate with benchmarks on zirconium dioxide nanoparticles of over 160000 atoms.

Catalytic mechanism of the PrhA (V150L/A232S) double mutant involved in the fungal meroterpenoid biosynthetic pathway: a QM/MM study

QM/MM calculations reveal the mechanism of a nonheme Fe(ii)/α-ketoglutarate-dependent oxygenase involved in the fungal meroterpenoid biosynthetic pathway.

The catalytic mechanism of S-acyltransferases: acylation is triggered on by a loose transition state and deacylation is turned off by a tight transition state

The reaction of S-acyltransferase is characterized by a loose transition state.

Conversion mechanism of enoyl thioesters into acyl thioesters catalyzed by 2-enoyl-thioester reductases from Candida Tropicalis

Enoyl thioester reductase from Candida tropicalis (Etr1p) catalyzes the NADPH-dependent conversion of enoyl thioesters into acyl thioesters, which are essential in fatty acid and second metabolite biosynthesis.

Recent examples of α-ketoglutarate-dependent mononuclear non-haem iron enzymes in natural product biosyntheses

Covering: up to 2018 α-Ketoglutarate (αKG, also known as 2-oxoglutarate)-dependent mononuclear non-haem iron (αKG-NHFe) enzymes catalyze a wide range of biochemical reactions, including hydroxylation, ring fragmentation, C-C bond cleavage, epimerization, desaturation, endoperoxidation and heterocycle formation. These enzymes utilize iron(ii) as the metallo-cofactor and αKG as the co-substrate. Herein, we summarize several novel αKG-NHFe enzymes involved in natural product biosyntheses discovered in recent years, including halogenation reactions, amino acid modifications and tailoring reactions in the biosynthesis of terpenes, lipids, fatty acids and phosphonates. We also conducted a survey of the currently available structures of αKG-NHFe enzymes, in which αKG binds to the metallo-centre bidentately through either a proximal- or distal-type binding mode. Future structure-function and structure-reactivity relationship investigations will provide crucial information regarding how activities in this large class of enzymes have been fine-tuned in nature.

On how the binding cavity of AsqJ dioxygenase controls the desaturation reaction regioselectivity: a QM/MM study

The Fe(II)/2-oxoglutarate-dependent dioxygenase AsqJ from Aspergillus nidulans catalyses two pivotal steps in the synthesis of quinolone antibiotic 4'-methoxyviridicatin, i.e., desaturation and epoxidation of a benzodiazepinedione. The previous experimental results signal that, during the desaturation reaction, hydrogen atom transfer (HAT) from the benzylic carbon atom (C10) is a more likely step to initiate the reaction than the alternative HAT from the ring moiety (C3 atom). To unravel the origins of this regioselectivity and to explain why the observed reaction is desaturation and not the "default" hydroxylation, we performed a QM/MM study on the reaction catalysed by AsqJ. Herein, we report results that complement the experimental findings and suggest that HAT at the C10 position is the preferred reaction due to favourable interactions between the substrate and the binding cavity that compensate for the relatively high intrinsic barrier associated with the process. For the resultant radical intermediate, the desaturation/hydroxylation selectivity is governed by electronic properties of the reactants, i.e., the energy gap between the orbital that hosts the unpaired electron and the sigma orbital for the C-H bond as well as the gap between the orbitals mixing in transition state structures for each elementary step. Regiospecificity of the AsqJ dehydrogenation reaction is dictated by substrate-protein interactions. 82 × 44 mm (300 × 300 dpi).

Amazing Diversity in Biochemical Roles of Fe(II)/2-Oxoglutarate Oxygenases

Since their discovery in the 1960s, the family of Fe(II)/2-oxoglutarate-dependent oxygenases has undergone a tremendous expansion to include enzymes catalyzing a vast diversity of biologically important reactions. Recent examples highlight roles in controlling chromatin modification, transcription, mRNA demethylation, and mRNA splicing. Others generate modifications in tRNA, translation factors, ribosomes, and other proteins. Thus, oxygenases affect all components of molecular biology's central dogma, in which information flows from DNA to RNA to proteins. These enzymes also function in biosynthesis and catabolism of cellular metabolites, including antibiotics and signaling molecules. Due to their critical importance, ongoing efforts have targeted family members for the development of specific therapeutics. This review provides a general overview of recently characterized oxygenase reactions and their key biological roles.

Catalytic mechanism and molecular engineering of quinolone biosynthesis in dioxygenase AsqJ

The recently discovered Fe^II/α-ketoglutarate-dependent dioxygenase AsqJ from Aspergillus nidulans stereoselectively catalyzes a multistep synthesis of quinolone alkaloids, natural products with significant biomedical applications. To probe molecular mechanisms of this elusive catalytic process, we combine here multi-scale quantum and classical molecular simulations with X-ray crystallography, and in vitro biochemical activity studies. We discover that methylation of the substrate is essential for the activity of AsqJ, establishing molecular strain that fine-tunes π-stacking interactions within the active site. To rationally engineer AsqJ for modified substrates, we amplify dispersive interactions within the active site. We demonstrate that the engineered enzyme has a drastically enhanced catalytic activity for non-methylated surrogates, confirming our computational data and resolved high-resolution X-ray structures at 1.55 Å resolution. Our combined findings provide crucial mechanistic understanding of the function of AsqJ and showcase how combination of computational and experimental data enables to rationally engineer enzymes.

The catalytic activity of dioxygenase AsqJ is strictly relying on the methylation of quinolone substrates. Here, the authors apply molecular simulations, X-ray crystallography and in vitro biochemical studies to the engineering of dioxygenase AsqJ with improved catalytic activity for modified non-methylated surrogates.

Insights into the Desaturation of Cyclopeptin and its C3 Epimer Catalyzed by a non-Heme Iron Enzyme: Structural Characterization and Mechanism Elucidation

AsqJ, an iron(II)- and 2-oxoglutarate-dependent enzyme found in viridicatin-type alkaloid biosynthetic pathways, catalyzes sequential desaturation and epoxidation to produce cyclopenins. Crystal structures of AsqJ bound to cyclopeptin and its C3 epimer are reported. Meanwhile, a detailed mechanistic study was carried out to decipher the desaturation mechanism. These findings suggest that a pathway involving hydrogen atom abstraction at the C10 position of the substrate by a short-lived Fe^IV -oxo species and the subsequent formation of a carbocation or a hydroxylated intermediate is preferred during AsqJ-catalyzed desaturation.

Unique chemistry of non-heme iron enzymes in fungal biosynthetic pathways

Reactions by non-heme iron enzymes in structurally intriguing fungal natural products pathways are summarized and discussed.