Quantum mechanics/molecular mechanics study on the oxygen binding and substrate hydroxylation step in AlkB repair enzymes

Effects of Clinical Mutations in the Second Coordination Sphere and Remote Regions on the Catalytic Mechanism of Non-Heme Fe(II)/2-Oxoglutarate-Dependent Aspartyl Hydroxylase AspH

Aspartyl/asparaginyl hydroxylase (AspH) catalyzes the post-translational hydroxylations of vital human proteins, playing an essential role in maintaining their biological functions. Single-point mutations in the Second Coordination Sphere (SCS) and long-range (LR) residues of AspH have been linked to pathological conditions such as the ophthalmologic condition Traboulsi syndrome and chronic kidney disease (CKD). Although the clinical impacts of these mutations are established, there is a critical knowledge gap regarding their specific atomistic effects on the catalytic mechanism of AspH. In this study, we report integrated computational investigations on the potential mechanistic implications of four mutant forms of human AspH with clinical importance: R735W, R735Q, R688Q, and G434V. All the mutant forms exhibited altered binding interactions with the co-substrate 2-oxoglutarate (2OG) and the main substrate in the ferric-superoxo and ferryl complexes, which are critical for catalysis, compared to the wild-type (WT). Importantly, the mutations strongly influence the energetics of the frontier molecular orbitals (FMOs) and, thereby, the activation energies for the hydrogen atom transfer (HAT) step compared to the WT AspH. Insights from our study can contribute to enzyme engineering and the development of selective modulators for WT and mutants of AspH, ultimately aiding in treating cancers, Traboulsi syndrome and, CKD.

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.

The Unique Role of the Second Coordination Sphere to Unlock and Control Catalysis in Nonheme Fe(II)/2-Oxoglutarate Histone Demethylase KDM2A

Nonheme Fe(II) and 2-oxoglutarate (2OG)-dependent histone lysine demethylases 2A (KDM2A) catalyze the demethylation of the mono- or dimethylated lysine 36 residue in the histone H3 peptide (H3K36me1/me2), which plays a crucial role in epigenetic regulation and can be involved in many cancers. Although the overall catalytic mechanism of KDMs has been studied, how KDM2 catalysis takes place in contrast to other KDMs remains unknown. Understanding such differences is vital for enzyme redesign and can help in enzyme-selective drug design. Herein, we employed molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) to explore the complete catalytic mechanism of KDM2A, including dioxygen diffusion and binding, dioxygen activation, and substrate oxidation. Our study demonstrates that the catalysis of KDM2A is controlled by the conformational change of the second coordination sphere (SCS), specifically by a change in the orientation of Y222, which unlocks the 2OG rearrangement from off-line to in-line mode. The study demonstrates that the variant Y222A makes the 2OG rearrangement more favorable. Furthermore, the study reveals that it is the size of H3K36me3 that prevents the 2OG rearrangement, thus rendering the enzyme inactivity with trimethylated lysine. Calculations show that the SCS and long-range interacting residues that stabilize the HAT transition state in KDM2A differ from those in KDM4A, KDM7B, and KDM6A, thus providing the basics for the enzyme-selective redesign and modulation of KDM2A without influencing other KDMs.

Unusual catalytic strategy by non-heme Fe(ii)/2-oxoglutarate-dependent aspartyl hydroxylase AspH

Biocatalytic C-H oxidation reactions are of important synthetic utility, provide a sustainable route for selective synthesis of important organic molecules, and are an integral part of fundamental cell processes. The multidomain non-heme Fe(ii)/2-oxoglutarate (2OG) dependent oxygenase AspH catalyzes stereoselective (3R)-hydroxylation of aspartyl- and asparaginyl-residues. Unusually, compared to other 2OG hydroxylases, crystallography has shown that AspH lacks the carboxylate residue of the characteristic two-His-one-Asp/Glu Fe-binding triad. Instead, AspH has a water molecule that coordinates Fe(ii) in the coordination position usually occupied by the Asp/Glu carboxylate. Molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies reveal that the iron coordinating water is stabilized by hydrogen bonding with a second coordination sphere (SCS) carboxylate residue Asp721, an arrangement that helps maintain the six coordinated Fe(ii) distorted octahedral coordination geometry and enable catalysis. AspH catalysis follows a dioxygen activation-hydrogen atom transfer (HAT)-rebound hydroxylation mechanism, unusually exhibiting higher activation energy for rebound hydroxylation than for HAT, indicating that the rebound step may be rate-limiting. The HAT step, along with substrate positioning modulated by the non-covalent interactions with SCS residues (Arg688, Arg686, Lys666, Asp721, and Gln664), are essential in determining stereoselectivity, which likely proceeds with retention of configuration. The tetratricopeptide repeat (TPR) domain of AspH influences substrate binding and manifests dynamic motions during catalysis, an observation of interest with respect to other 2OG oxygenases with TPR domains. The results provide unique insights into how non-heme Fe(ii) oxygenases can effectively catalyze stereoselective hydroxylation using only two enzyme-derived Fe-ligating residues, potentially guiding enzyme engineering for stereoselective biocatalysis, thus advancing the development of non-heme Fe(ii) based biomimetic C-H oxidation catalysts, and supporting the proposal that the 2OG oxygenase superfamily may be larger than once perceived.

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.

Catalysis by KDM6 Histone Demethylases - A Synergy between the Non-Heme Iron(II) Center, Second Coordination Sphere, and Long-Range Interactions

KDM6A (UTX) and KDM6B (JMJD3) are human non-heme Fe(II) and 2-oxoglutarate (2OG) dependent JmjC oxygenases that catalyze the demethylation of trimethylated lysine 27 in the N-terminal tail of histone H3, a post-translational modification that regulates transcription. A Combined Quantum Mechanics/ Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) study on the catalytic mechanism of KDM6A/B reveals that the transition state for the rate-limiting hydrogen atom transfer (HAT) reaction in KDM6A catalysis is stabilized by polar (Asn217) and aromatic (Trp369)/non-polar (Pro274) residues in contrast to KDM4, KDM6B and KDM7 demethylases where charged residues (Glu, Arg, Asp) are involved. KDM6A employs both σ- and π-electron transfer pathways for HAT, whereas KDM6B employs the σ-electron pathway. Differences in hydrogen bonding of the Fe-chelating Glu252(KDM6B) contribute to the lower energy barriers in KDM6B vs. KDM6A. The study reveals a dependence of the activation barrier of the rebound hydroxylation on the Fe-O-C angle in the transition state of KDM6A. Anti-correlation of the Zn-binding domain with the active site residues is a key factor distinguishing KDM6A/B from KDM7/4s. The results reveal the importance of communication between the Fe center, second coordination sphere, and long-range interactions in catalysis by KDMs and, by implication, other 2OG oxygenases.

Kinetic Studies on the 2-Oxoglutarate/Fe(II)-Dependent Nucleic Acid Modifying Enzymes from the AlkB and TET Families

Nucleic acid methylations are important genetic and epigenetic biomarkers. The formation and removal of these markers is related to either methylation or demethylation. In this review, we focus on the demethylation or oxidative modification that is mediated by the 2-oxoglutarate (2-OG)/Fe(II)-dependent AlkB/TET family enzymes. In the catalytic process, most enzymes oxidize 2-OG to succinate, in the meantime oxidizing methyl to hydroxymethyl, leaving formaldehyde and generating demethylated base. The AlkB enzyme from Escherichia coli has nine human homologs (ALKBH1-8 and FTO) and the TET family includes three members, TET1 to 3. Among them, some enzymes have been carefully studied, but for certain enzymes, few studies have been carried out. This review focuses on the kinetic properties of those 2-OG/Fe(II)-dependent enzymes and their alkyl substrates. We also provide some discussions on the future directions of this field.

Quantum based effects of therapeutic nuclear magnetic resonance persistently reduce glycolysis

Electromagnetic fields are known to induce the clock protein cryptochrome to modulate intracellular reactive oxygen species (ROS) via the quantum based radical pair mechanism (RPM) in mammalian cells. Recently, therapeutic Nuclear Magnetic Resonance (tNMR) was shown to alter protein levels of the circadian clock associated Hypoxia Inducible Factor-1α (HIF-1α) in a nonlinear dose response relationship. Using synchronized NIH3T3 cells, we show that tNMR under normoxia and hypoxia persistently modifies cellular metabolism. After normoxic tNMR treatment, glycolysis is reduced, as are lactate production, extracellular acidification rate, the ratio of ADP/ATP and cytosolic ROS, whereas mitochondrial and extracellular ROS, as well as cellular proliferation are increased. Remarkably, these effects are even more pronounced after hypoxic tNMR treatment, driving cellular metabolism to a reduced glycolysis while mitochondrial respiration is kept constant even during reoxygenation. Hence, we propose tNMR as a potential therapeutic tool in ischemia driven diseases like inflammation, infarct, stroke and cancer.

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.

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.

Properties of the Reactants and Their Interactions within and with the Enzyme Binding Cavity Determine Reaction Selectivities. The Case of Fe(II)/2-Oxoglutarate Dependent Enzymes

Fe(II)/2-oxoglutarate dependent dioxygenases (ODDs) share a double stranded beta helix (DSBH) fold and utilise a common reactive intermediate, ferryl species, to catalyse oxidative transformations of substrates. Despite the structural similarities, ODDs accept a variety of substrates and facilitate a wide range of reactions, that is hydroxylations, desaturations, (oxa)cyclisations and ring rearrangements. In this review we present and discuss the factors contributing to the observed (regio)selectivities of ODDs. They span from inherent properties of the reactants, that is, substrate molecule and iron cofactor, to the interactions between the substrate and the enzyme's binding cavity; the latter can counterbalance the effect of the former. Based on results of both experimental and computational studies dedicated to ODDs, we also line out the properties of the reactants which promote reaction outcomes other than the "default" hydroxylation. It turns out that the reaction selectivity depends on a delicate balance of interactions between the components of the investigated system.

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 C1 -position. Recent computational studies showed that in the gas-phase the C2 -H bond of taurine is substantially weaker than the C1 -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 C1 -H and C2 -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 C1 -H bond and triggers a chemoselective reaction process. The large cluster models reproduce the experimental free energy of activation excellently.

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.

Biodegradation of herbicides by a plant nonheme iron dioxygenase: mechanism and selectivity of substrate analogues

Aryloxyalkanoate dioxygenases are unique herbicide biodegrading nonheme iron enzymes found in plants and hence, from environmental and agricultural point of view they are important and valuable. However, they often are substrate specific and little is known on the details of the mechanism and the substrate scope. To this end, we created enzyme models and calculate the mechanism for 2,4-dichlorophenoxyacetic acid biodegradation and 2-methyl substituted analogs by density functional theory. The work shows that the substrate binding is tight and positions the aliphatic group close to the metal center to enable a chemoselective reaction mechanism to form the C 2 -hydroxy products, whereas the aromatic hydroxylation barriers are well higher in energy. Subsequently, we investigated the metabolism of R - and S -methyl substituted inhibitors and show that these do not react as efficiently as 2,4-dichlorophenoxyacetic acid substrate due to stereochemical clashes in the active site and particularly for the R -isomer give high rebound barriers.

Computational investigations of selected enzymes from two iron and α-ketoglutarate-dependent families

DNA alkylation is used as the key epigenetic mark in eukaryotes, however, most alkylation in DNA can result in deleterious effects. Therefore, this process needs to be tightly regulated. The enzymes of the AlkB and Ten-Eleven Translocation (TET) families are members of the Fe and alpha-ketoglutarate-dependent superfamily of enzymes that are tasked with dealkylating DNA and RNA in cells. Members of these families span all species and are an integral part of transcriptional regulation. While both families catalyze oxidative dealkylation of various bases, each has specific preference for alkylated base type as well as distinct catalytic mechanisms. This perspective aims to provide an overview of computational work carried out to investigate several members of these enzyme families including AlkB, ALKB Homolog 2, ALKB Homolog 3 and Ten-Eleven Translocate 2. Insights into structural details, mutagenesis studies, reaction path analysis, electronic structure features in the active site, and substrate preferences are presented and discussed.

What Is the Catalytic Mechanism of Enzymatic Histone N-Methyl Arginine Demethylation and Can It Be Influenced by an External Electric Field?

Arginine methylation is an important mechanism of epigenetic regulation. Some Fe(II) and 2-oxoglutarate dependent Jumonji-C (JmjC) Nϵ-methyl lysine histone demethylases also have N-methyl arginine demethylase activity. We report combined molecular dynamic (MD) and Quantum Mechanical/Molecular Mechanical (QM/MM) studies on the mechanism of N-methyl arginine demethylation by human KDM4E and compare the results with those reported for N-methyl lysine demethylation by KDM4A. At the KDM4E active site, Glu191, Asn291, and Ser197 form a conserved scaffold that restricts substrate dynamics; substrate binding is also mediated by an out of active site hydrogen-bond between the substrate Ser1 and Tyr178. The calculations imply that in either C-H or N-H potential bond cleaving pathways for hydrogen atom transfer (HAT) during N-methyl arginine demethylation, electron transfer occurs via a σ-channel; the transition state for the N-H pathway is ∼10 kcal/mol higher than for the C-H pathway due to the higher bond dissociation energy of the N-H bond. The results of applying external electric fields (EEFs) reveal EEFs with positive field strengths parallel to the Fe=O bond have a significant barrier-lowering effect on the C-H pathway, by contrast, such EEFs inhibit the N-H activation rate. The overall results imply that KDM4 catalyzed N-methyl arginine demethylation and N-methyl lysine demethylation occur via similar C-H abstraction and rebound mechanisms leading to methyl group hydroxylation, though there are differences in the interactions leading to productive binding of intermediates.

Thebaine is Selectively Demethylated by Thebaine 6-O-Demethylase and Codeine-3-O-demethylase at Distinct Binding Sites: A Computational Study

Two homologous 2-oxoglutarate-dependent (ODD) nonheme enzymes thebaine 6-O-demethylase (T6ODM) and codeine-3-O-demethylase (CODM), are involved in the morphine biosynthesis pathway from thebaine, catalyzing the O-demethylation reaction with precise regioselectivity at C6 and C3 positions of thebaine respectively. We investigated the origin of the regioselectivity of these enzymes by combined molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations and found that Thebaine binds at the two distinct sites of T6ODM and CODM, which determines the regioselectivity of the enzymes. A remarkable oxo rotation is observed in the decarboxylation process. Starting from the closed pentacoordinate configuration, the C-terminal lid adopts an open conformation in the octahedral Fe(IV) = O complex to facilitate the subsequent demethylation. Phe241 and Phe311 stabilize the substrate in the binding pocket, while Arg219 acts as a gatekeeper residue to stabilize the substrate. Our results unravel the regioselectivity in 2-OG dependent nonheme enzymes and may shed light for exploring the substrate scope of these enzymes and developing novel biotechnology for morphine biosynthesis.

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.

How external perturbations affect the chemoselectivity of substrate activation by cytochrome P450 OleTJE

Cytochrome P450 enzymes are versatile biocatalysts found in most forms of life. Generally, the cytochrome P450s react with dioxygen and hence are haem-based mono-oxygenases; however, in specific isozymes, H2O2 rather than O2 is used and these P450s act as peroxygenases. The P450 OleTJE is a peroxygenase that binds long to medium chain fatty acids and converts them to a range of products originating from Cα-hydroxylation, Cβ-hydroxylation, Cα-Cβ desaturation and decarboxylation of the substrate. There is still controversy regarding the details of the reaction mechanism of P450 OleTJE; how the products are formed and whether the product distributions can be influenced by external perturbations. To gain further insights into the structure and reactivity of P450 OleTJE, we set up a range of large active site model complexes as well as full enzymatic structures and did a combination of density functional theory studies and quantum mechanics/molecular mechanics calculations. In particular, the work focused on the mechanisms leading to these products under various reaction conditions. Thus, for a small cluster model, we find a highly selective Cα-hydroxylation pathway that is preferred over Cβ-H hydrogen atom abstraction by at least 10 kcal mol-1. Introduction of polar residues to the model, such as an active site protonated histidine residue or through external electric field effects, lowers the Cβ-H hydrogen atom abstraction barriers are lowered, while a full QM/MM model brings the Cα-H and Cβ-H hydrogen atom abstraction barriers within 1 kcal mol-1. Our studies; therefore, implicate that environmental effects in the second-coordination sphere can direct and guide selectivities in enzymatic reaction mechanisms.

The molecular mechanism of P450-catalyzed amination of the pyrrolidine derivative of lidocaine: insights from multiscale simulations

Nitrogen heterocycles are key and prevalent motifs in drugs. Evolved variants of cytochrome P450_BM3 (CYP102A1) from Bacillus megaterium employ high-valent oxo-iron(iv) species to catalyze the synthesis of imidazolidine-4-ones via an intramolecular C–H amination. Herein, we use multi-scale simulations, including classical molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical (QM/MM) calculations and QM calculations, to reveal the molecular mechanism of the intramolecular C–H amination of the pyrrolidine derivative of lidocaine bearing cyclic amino moieties catalyzed by the variant RP/FV/EV of P450_BM3, which bears five mutations compared to wild type. Our calculations show that overall catalysis includes both the enzymatic transformation in P450 and non-enzymatic transformation in water solution. The enzymatic transformation involves the exclusive hydroxylation of the C–H bond of the pyrrolidine derivative of lidocaine, leading to the hydroxylated intermediate, during which the substrate radical would be bypassed. The following dehydration and C–N coupling reactions are found to be much favored in aqueous situation compared to that in the non-polar protein environment. The present findings expand our understanding of the P450-catalyzed C(sp₃)–H amination reaction.

Nitrogen heterocycles are key and prevalent motifs in drugs.

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.

Regioselectivity of hyoscyamine 6β-hydroxylase-catalysed hydroxylation as revealed by high-resolution structural information and QM/MM calculations

Hyoscyamine 6β-hydroxylase (H6H) is a bifunctional non-heme 2-oxoglutarate/Fe2+-dependent dioxygenase that catalyzes the two final steps in the biosynthesis of scopolamine. Based on high resolution crystal structures of H6H from Datura metel, detailed information on substrate binding was obtained that provided insights into the onset of the enzymatic process. In particular, the role of two prominent residues was revealed - Glu-116 that interacts with the tertiary amine located on the hyoscyamine tropane moiety and Tyr-326 that forms CH-π hydrogen bonds with the hyoscyamine phenyl ring. The structures were used as the basis for QM/MM calculations that provided an explanation for the regioselectivity of the hydroxylation reaction on the hyoscyamine tropane moiety (C6 vs. C7) and quantified contributions of active site residues to respective barrier heights.

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

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

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.

Catalysis by the JmjC histone demethylase KDM4A integrates substrate dynamics, correlated motions and molecular orbital control

The N^ε-methyl lysine status of histones is important in the regulation of eukaryotic transcription. The Fe(ii) and 2-oxoglutarate (2OG) -dependent JmjC domain enzymes are the largest family of histone N^ε-methyl lysine demethylases (KDMs). The human KDM4 subfamily of JmjC KDMs is linked with multiple cancers and some of its members are medicinal chemistry targets. We describe the use of combined molecular dynamics (MD) and Quantum Mechanical/Molecular Mechanical (QM/MM) methods to study the mechanism of KDM4A, which catalyzes demethylation of both tri- and di-methylated forms of histone H3 at K9 and K36. The results show that the oxygen activation at the active site of KDM4A is optimized towards the generation of the reactive Fe(iv)-oxo intermediate. Factors including the substrate binding mode, correlated motions of the protein and histone substrates, and molecular orbital control synergistically contribute to the reactivity of the Fe(iv)-oxo intermediate. In silico substitutions were performed to investigate the roles of residues (Lys241, Tyr177, and Asn290) in substrate orientation. The Lys241Ala substitution abolishes activity due to altered substrate orientation consistent with reported experimental studies. Calculations with a macrocyclic peptide substrate analogue reveal that induced conformational changes/correlated motions in KDM4A are sequence-specific in a manner that influences substrate binding affinity. Second sphere residues, such as Ser288 and Thr289, may contribute to KDM4A catalysis by correlated motions with active site residues. Residues that stabilize key intermediates, and which are predicted to be involved in correlated motions with other residues in the second sphere and beyond, are shown to be different in KDM4A compared to those in another JmjC KDM (PHF8), which acts on H3K9 di- and mono-methylated forms, suggesting that allosteric type inhibition is of interest from the perspective of developing selective JmjC KDM inhibitors.

The second sphere residues and regions of the protein in histone demethylase enzymes that makes correlated motion with the active site contribute to efficient catalysis.

Role of Structural Dynamics in Selectivity and Mechanism of Non-heme Fe(II) and 2-Oxoglutarate-Dependent Oxygenases Involved in DNA Repair

AlkB and its human homologue AlkBH2 are Fe(II)- and 2-oxoglutarate (2OG)-dependent oxygenases that repair alkylated DNA bases occurring as a consequence of reactions with mutagenic agents. We used molecular dynamics (MD) and combined quantum mechanics/molecular mechanics (QM/MM) methods to investigate how structural dynamics influences the selectivity and mechanisms of the AlkB- and AlkBH2-catalyzed demethylation of 3-methylcytosine (m3C) in single (ssDNA) and double (dsDNA) stranded DNA. Dynamics studies reveal the importance of the flexibility in both the protein and DNA components in determining the preferences of AlkB for ssDNA and of AlkBH2 for dsDNA. Correlated motions, including of a hydrophobic β-hairpin, are involved in substrate binding in AlkBH2-dsDNA. The calculations reveal that 2OG rearrangement prior to binding of dioxygen to the active site Fe is preferred over a ferryl rearrangement to form a catalytically productive Fe(IV)=O intermediate. Hydrogen atom transfer proceeds via a σ-channel in AlkBH2-dsDNA and AlkB-dsDNA; in AlkB-ssDNA, there is a competition between σ- and π-channels, implying that the nature of the complexed DNA has potential to alter molecular orbital interactions during the substrate oxidation. Our results reveal the importance of the overall protein-DNA complex in determining selectivity and how the nature of the substrate impacts the mechanism.

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

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

Catalysis by the Non-Heme Iron(II) Histone Demethylase PHF8 Involves Iron Center Rearrangement and Conformational Modulation of Substrate Orientation

PHF8 (KDM7B) is a human non-heme 2-oxoglutarate (2OG) JmjC domain oxygenase that catalyzes the demethylation of the di/mono-Nε-methylated K9 residue of histone H3. Altered PHF8 activity is linked to genetic diseases and cancer; thus, it is an interesting target for epigenetic modulation. We describe the use of combined quantum mechanics/molecular mechanics (QM/MM) and molecular dynamics (MD) simulations to explore the mechanism of PHF8, including dioxygen activation, 2OG binding modes, and substrate demethylation steps. A PHF8 crystal structure manifests the 2OG C-1 carboxylate bound to iron in a nonproductive orientation, i.e., trans to His247. A ferryl-oxo intermediate formed by activating dioxygen bound to the vacant site in this complex would be nonproductive, i.e., "off-line" with respect to reaction with Nε-methylated K9. We show rearrangement of the "off-line" ferryl-oxo intermediate to a productive "in-line" geometry via a solvent exchange reaction (called "ferryl-flip") is energetically unfavorable. The calculations imply that movement of the 2OG C-1 carboxylate prior to dioxygen binding at a five-coordination stage in catalysis proceeds with a low barrier, suggesting that two possible 2OG C-1 carboxylate geometries can coexist at room temperature. We explored alternative mechanisms for hydrogen atom transfer and show that second sphere interactions orient the Nε-methylated lysine in a conformation where hydrogen abstraction from a methyl C-H bond is energetically more favorable than hydrogen abstraction from the N-H bond of the protonated Nε-methyl group. Using multiple HAT reaction path calculations, we demonstrate the crucial role of conformational flexibility in effective hydrogen transfer. Subsequent hydroxylation occurs through a rebound mechanism, which is energetically preferred compared to desaturation, due to second sphere interactions. The overall mechanistic insights reveal the crucial role of iron-center rearrangement, second sphere interactions, and conformational flexibility in PHF8 catalysis and provide knowledge useful for the design of mechanism-based PHF8 inhibitors.

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.

On the reaction mechanism of an endoperoxide ring formation by fumitremorgin B endoperoxidase. The right arrangement makes a difference

Fumitremorgin B endoperoxidase (FtmF) belongs to 2-oxoglutarate dependent dioxygenases and catalyzes an unusual oxidative reaction of endoperoxide formation at the final stage of biosynthesis of verruculogen - a mycotoxin produced by Aspergillus and Penicillinum strains. The published crystal structure of FtmF (PDB: ), which is of overall good quality, contains a model of the substrate bound in the active site, which, however, has very low occupancy and its conformation does not comply with the small molecule crystal structure. Moreover, a previous computational study that employed a model based on this crystal structure revealed a substantial reaction barrier, which might indicate that the model of FtmF/substrate complex can have serious errors. The purpose of this work was to model with computational methods the structure of the enzyme-substrate complex and to investigate the mechanisms of the enzymatic reaction. Docking, molecular dynamics simulation and DFT results, all indicate the substrate most likely binds in the active site in a configuration very different from that originally suggested. Moreover, for this newly proposed structure of the enzyme-substrate complex, the reaction energy profile is characterised exclusively by low barriers and it successfully explains the observed regiospecificity of the enzymatic process. Finally, a plausible binding site for ascorbate was found and it is suggested that ascorbate is involved in the final step of the FtmF reaction.

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.

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.

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.

Bridging the Homogeneous-Heterogeneous Divide: Modeling Spin for Reactivity in Single Atom Catalysis

Single atom catalysts (SACs) are emergent catalytic materials that have the promise of merging the scalability of heterogeneous catalysts with the high activity and atom economy of homogeneous catalysts. Computational, first-principles modeling can provide essential insight into SAC mechanism and active site configuration, where the sub-nm-scale environment can challenge even the highest-resolution experimental spectroscopic techniques. Nevertheless, the very properties that make SACs attractive in catalysis, such as localized d electrons of the isolated transition metal center, make them challenging to study with conventional computational modeling using density functional theory (DFT). For example, Fe/N-doped graphitic SACs have exhibited spin-state dependent reactivity that remains poorly understood. However, spin-state ordering in DFT is very sensitive to the nature of the functional approximation chosen. In this work, we develop accurate benchmarks from correlated wavefunction theory (WFT) for relevant octahedral complexes. We use those benchmarks to evaluate optimal DFT functional choice for predicting spin state ordering in small octahedral complexes as well as models of pyridinic and pyrrolic nitrogen environments expected in larger SACs. Using these guidelines, we determine Fe/N-doped graphene SAC model properties and reactivity as well as their sensitivities to DFT functional choice. Finally, we conclude with broad recommendations for computational modeling of open-shell transition metal single-atom catalysts.

Advances in Sustainable Catalysis: A Computational Perspective

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

Conformational flexibility influences structure-function relationships in nucleic acid N-methyl demethylases

N-Methylation of DNA/RNA bases can be regulatory or damaging and is linked to diseases including cancer and genetic disorders. Bacterial AlkB and human FTO are DNA/RNA demethylases belonging to the Fe(ii) and 2-oxoglutarate oxygenase superfamily. Modelling studies reveal conformational dynamics influence structure-function relationships of AlkB and FTO, e.g. why 1-methyladenine is a better substrate for AlkB than 6-methyladenine. Simulations show that the flexibility of the double stranded DNA substrate in AlkB influences correlated motions, including between the core jelly-roll fold and an active site loop involved in substrate binding. The FTO N- and C-terminal domains move in respect to one another in a manner likely important for substrate binding. Substitutions, including clinically observed ones, influencing catalysis contribute to the network of correlated motions in AlkB and FTO. Overall, the calculations highlight the importance of the overall protein environment and its flexibility to the geometry of the reactant complexes.

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.