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

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.

An Unusual Ferryl Intermediate and Its Implications for the Mechanism of Oxacyclization by the Loline-Producing Iron(II)- and 2-Oxoglutarate-Dependent Oxygenase, LolO

N-Acetylnorloline synthase (LolO) is one of several iron(II)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenases that catalyze sequential reactions of different types in the biosynthesis of valuable natural products. LolO hydroxylates C2 of 1-exo-acetamidopyrrolizidine before coupling the C2-bonded oxygen to C7 to form the tricyclic loline core. Each reaction requires cleavage of a C-H bond by an oxoiron(IV) (ferryl) intermediate; however, different carbons are targeted, and the carbon radicals have different fates. Prior studies indicated that the substrate-cofactor disposition (SCD) controls the site of H· abstraction and can affect the reaction outcome. These indications led us to determine whether a change in SCD from the first to the second LolO reaction might contribute to the observed reactivity switch. Whereas the single ferryl complex in the C2 hydroxylation reaction was previously shown to have typical Mössbauer parameters, one of two ferryl complexes to accumulate during the oxacyclization reaction has the highest isomer shift seen to date for such a complex and abstracts H· from C7 ∼ 20 times faster than does the first ferryl complex in its previously reported off-pathway hydroxylation of C7. The detectable hydroxylation of C7 in competition with cyclization by the second ferryl complex is not enhanced in 2H2O solvent, suggesting that the C2 hydroxyl is deprotonated prior to C7-H cleavage. These observations are consistent with the coordination of the C2 oxygen to the ferryl complex, which may reorient its oxo ligand, the substrate, or both to positions more favorable for C7-H cleavage and oxacyclization.

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.

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.

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.

Synthesis, Structure and Reactivity of a Mononuclear N,N,O-Bound Fe(II) α-Keto-Acid Complex

A bulky, tridentate phenolate ligand (ImPh2 NNOtBu ) was used to synthesise the first example of a mononuclear, facial, N,N,O-bound iron(II) benzoylformate complex, [Fe(ImPh2 NNOtBu )(BF)] (2). The X-ray crystal structure of 2 reveals that the iron centre is pentacoordinate (τ=0.5), with a vacant site located cis to the bidentate BF ligand. The Mössbauer parameters of 2 are consistent with high-spin iron(II), and are very close to those reported for α-ketoglutarate-bound non-heme iron enzyme active sites. According to NMR and UV-vis spectroscopies, the structural integrity of 2 is retained in both coordinating and non-coordinating solvents. Cyclic voltammetry studies show that the iron centre has a very low oxidation potential and is more prone to electrochemical oxidation than the redox-active phenolate ligand. Complex 2 reacts with NO to form a S=3 /2 {FeNO}7 adduct in which NO binds directly to the iron centre, according to EPR, UV-vis, IR spectroscopies and DFT analysis. Upon O2 exposure, 2 undergoes oxidative decarboxylation to form a diiron(III) benzoate complex, [Fe2 (ImPh2 NNOtBu )2 (μ2 -OBz)(μ2 -OH)2 ]+ (3). A small amount of hydroxylated ligand was also observed by ESI-MS, hinting at the formation of a high-valent iron(IV)-oxo intermediate. Initial reactivity studies show that 2 is capable of oxygen atom transfer reactivity with O2 , converting methyl(p-tolyl)sulfide to sulfoxide.

Why Nonheme Iron Halogenases Do Not Fluorinate C-H Bonds: A Computational Investigation

Selective halogenation is necessary for a range of fine chemical applications, including the development of therapeutic drugs. While synthetic processes to achieve C-H halogenation require harsh conditions, enzymes such as nonheme iron halogenases carry out some types of C-H halogenation, i.e., chlorination or bromination, with ease, while others, i.e., fluorination, have never been observed in natural or engineered nonheme iron enzymes. Using density functional theory and correlated wave function theory, we investigate the differences in structural and energetic preferences of the smaller fluoride and the larger chloride or bromide intermediates throughout the catalytic cycle. Although we find that the energetics of rate-limiting hydrogen atom transfer are not strongly impacted by fluoride substitution, the higher barriers observed during the radical rebound reaction for fluoride relative to chloride and bromide contribute to the difficulty of C-H fluorination. We also investigate the possibility of isomerization playing a role in differences in reaction selectivity, and our calculations reveal crucial differences in terms of isomer energetics of the key ferryl intermediate between fluoride and chloride/bromide intermediates. While formation of monodentate isomers believed to be involved in selective catalysis is shown for chloride and bromide intermediates, we find that formation of the fluoride monodentate intermediate is not possible in our calculations, which lack additional stabilizing interactions with the greater protein environment. Furthermore, the shorter Fe-F bonds are found to increase isomerization reaction barriers, suggesting that incorporation of residues that form a halogen bond with F and elongate Fe-F bonds could make selective C-H fluorination possible in nonheme iron halogenases. Our work highlights the differences between the fluoride and chloride/bromide intermediates and suggests potential steps toward engineering nonheme iron halogenases to enable selective C-H fluorination.

Impacts of QM region sizes and conformation numbers on modelling enzyme reactions: a case study of polyethylene terephthalate hydrolase

A quantum mechanics/molecular mechanics (QM/MM) approach is a broadly used tool in computational enzymology. Treating the QM region with a high-level DFT method is one of the important branches. Here, taking leaf-branch compost cutinase-catalyzed polyethylene terephthalate depolymerization as an example, the convergence behavior of energy barriers as well as key structural and charge features with respect to the size of the QM region (up to 1000 atoms) is systematically investigated. BP86/6-31G(d)//CHARMM and M06-2X/6-311G(d,p)//CHARMM level of theories were applied for geometry optimizations and single-point energy calculations, respectively. Six independent enzyme conformations for all the four catalytic steps (steps (i)-(iv)) were considered. Most of the twenty-four cases show that at least 500 QM atoms are needed while only two rare cases show that ∼100 QM atoms are sufficient for convergence when only a single conformation was considered. This explains why most previous studies showed that 500 or more QM atoms are required while a few others showed that ∼100 QM atoms are sufficient for DFT/MM calculations. More importantly, average energy barriers and key structural/charge features from six conformations show an accelerated convergence than that in a single conformation. For instance, to reach energy barrier convergence (within 2.0 kcal mol-1) for step (ii), only ∼100 QM atoms are required if six conformations are considered while 500 or more QM atoms are needed with a single conformation. The convergence is accelerated to be more rapid if hundreds and thousands of conformations were considered, which aligns with previous findings that only several dozens of QM atoms are required for convergence with semi-empirical QM/MM MD simulations.

From random to rational: improving enzyme design through electric fields, second coordination sphere interactions, and conformational dynamics

Enzymes are versatile and efficient biological catalysts that drive numerous cellular processes, motivating the development of enzyme design approaches to tailor catalysts for diverse applications. In this perspective, we investigate the unique properties of natural, evolved, and designed enzymes, recognizing their strengths and shortcomings. We highlight the challenges and limitations of current enzyme design protocols, with a particular focus on their limited consideration of long-range electrostatic and dynamic effects. We then delve deeper into the impact of the protein environment on enzyme catalysis and explore the roles of preorganized electric fields, second coordination sphere interactions, and protein dynamics for enzyme function. Furthermore, we present several case studies illustrating successful enzyme-design efforts incorporating enzyme strategies mentioned above to achieve improved catalytic properties. Finally, we envision the future of enzyme design research, spotlighting the challenges yet to be overcome and the synergy of intrinsic electric fields, second coordination sphere interactions, and conformational dynamics to push the state-of-the-art boundaries.

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.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

Kinetic and inhibition studies on human Jumonji-C (JmjC) domain-containing protein 5

Jumonji-C (JmjC) domain-containing protein 5 (JMJD5) is a human 2-oxoglutarate (2OG) and Fe(ii)-dependent oxygenase which catalyses the post-translational C3 hydroxylation of arginyl-residues and which is linked to the circadian rhythm and to cancer biology through as yet unidentified mechanisms. We report robust solid phase extraction coupled to mass spectrometry (SPE-MS)-based JMJD5 assays which enable kinetic and high-throughput inhibition studies. The kinetic studies reveal that some synthetic 2OG derivatives, notably including a 2OG derivative with a cyclic carbon backbone (i.e. (1R)-3-(carboxycarbonyl)cyclopentane-1-carboxylic acid), are efficient alternative cosubstrates of JMJD5 and of factor inhibiting hypoxia-inducible transcription factor HIF-α (FIH), but not of the Jumonji-C (JmjC) histone Nε-methyl lysine demethylase KDM4E, apparently reflecting the closer structural similarity of JMJD5 and FIH. The JMJD5 inhibition assays were validated by investigating the effect of reported 2OG oxygenase inhibitors on JMJD5 catalysis; the results reveal that broad-spectrum 2OG oxygenase inhibitors are also efficient JMJD5 inhibitors (e.g. N-oxalylglycine, pyridine-2,4-dicarboxylic acid, ebselen) whereas most 2OG oxygenase inhibitors that are in clinical use (e.g. roxadustat) do not inhibit JMJD5. The SPE-MS assays will help enable the development of efficient and selective JMJD5 inhibitors for investigating the biochemical functions of JMJD5 in cellular studies.

Dioxygen Binding Is Controlled by the Protein Environment in Non-heme FeII and 2-Oxoglutarate Oxygenases: A Study on Histone Demethylase PHF8 and an Ethylene-Forming Enzyme

This study investigates dioxygen binding and 2-oxoglutarate (2OG) coordination by two model non-heme FeII /2OG enzymes: a class 7 histone demethylase (PHF8) that catalyzes the hydroxylation of its H3K9me2 histone substrate leading to demethylation reactivity and the ethylene-forming enzyme (EFE), which catalyzes two competing reactions of ethylene generation and substrate l-Arg hydroxylation. Although both enzymes initially bind 2OG by using an off-line 2OG coordination mode, in PHF8, the substrate oxidation requires a transition to an in-line mode, whereas EFE is catalytically productive for ethylene production from 2OG in the off-line mode. We used classical molecular dynamics (MD), quantum mechanics/molecular mechanics (QM/MM) MD and QM/MM metadynamics (QM/MM-MetD) simulations to reveal that it is the dioxygen binding process and, ultimately, the protein environment that control the formation of the in-line FeIII -OO⋅- intermediate in PHF8 and the off-line FeIII -OO⋅- intermediate in EFE.

Spin-Forbidden Addition of Molecular Oxygen to Stable Enol Intermediates-Decarboxylation of 2-Methyl-1-tetralone-2-carboxylic Acid

The deprotonation of an organic substrate is a common preactivation step for the enzymatic cofactorless addition of O2 to this substrate, as it promotes charge-transfer between the two partners, inducing intersystem crossing between the triplet and singlet states involved in the process. Nevertheless, the spin-forbidden addition of O2 to uncharged ligands has also been observed in the laboratory, and the detailed mechanism of how the system circumvents the spin-forbiddenness of the reaction is still unknown. One of these examples is the cofactorless peroxidation of 2-methyl-3,4-dihydro-1-naphthol, which will be studied computationally using single and multi-reference electronic structure calculations. Our results show that the preferred mechanism is that in which O2 picks a proton from the substrate in the triplet state, and subsequently hops to the singlet state in which the product is stable. For this reaction, the formation of the radical pair is associated with a higher barrier than that associated with the intersystem crossing, even though the absence of the negative charge leads to relatively small values of the spin-orbit coupling.

Development of JmjC-domain-containing histone demethylase (KDM2-7) inhibitors for cancer therapy

Histone methylation is the most common histone modification and a highly dynamic regulator of gene transcription. Methylation of lysine residues can alter the structure of chromatin, helping to regulate DNA-based nuclear activities. Lysine demethylases control and maintain epigenetic factors that affect chromatin structure and cell characteristics. A variety of diseases, including malignant tumors, are connected to their dysregulation. Advances in biochemistry and pathogenesis have prompted the discovery of small molecule inhibitors and tool compounds that disrupt lysine demethylation. In this review, we focus on JmjC-domain-containing histone lysine demethylases (KDM2-7), discussing their structures and biological roles, representative inhibitors, and therapeutic potential in cancer therapy, and aiming to provide unique insights into the development of JmjC-KDM inhibitors.

JMJD8 Functions as a Novel AKT1 Lysine Demethylase

JMJD8 is a protein from the JMJD family that only has the JmjC domain. Studies on the function of JMJD8 indicate that JMJD8 is involved in signaling pathways, including AKT/NF-κB, and thus affects cell proliferation and development. Here, we reported the activity of JMJD8 as a non-histone demethylase. We investigated the demethylation of JMJD8 on trimethylated lysine of AKT1 in vivo and in vitro using trimethylated AKT1 short peptide and AKT1 protein, and we tracked the regulation of JMJD8 on AKT1 activity at the cellular level. The results showed that JMJD8, a mini lysine demethylase, altered AKT1 protein function via changing its degree of methylation.

Can Second Coordination Sphere and Long-Range Interactions Modulate Hydrogen Atom Transfer in a Non-Heme Fe(II)-Dependent Histone Demethylase?

Fe(II)-dependent oxygenases employ hydrogen atom transfer (HAT) to produce a myriad of products. Understanding how such enzymes use dynamic processes beyond the immediate vicinity of the active site to control the selectivity and efficiency of HAT is important for metalloenzyme engineering; however, obtaining such knowledge by experiments is challenging. This study develops a computational framework for identifying second coordination sphere (SCS) and especially long-range (LR) residues relevant for catalysis through dynamic cross-correlation analysis (DCCA) using the human histone demethylase PHF8 (KDM7B) as a model oxygenase. Furthermore, the study explores the mechanistic pathways of influence of the SCS and LR residues on the HAT reaction. To demonstrate the plausibility of the approach, we investigated the effect of a PHF8 F279S clinical mutation associated with X-linked intellectual disability, which has been experimentally shown to ablate PHF8-catalyzed demethylation. In agreement, the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) studies showed a change in the H3_1-14K9me2 substrate orientation and an increased HAT barrier. We systematically analyzed the pathways by which the identified SCS and LR residues may influence HAT by exploring changes in H3K9me2 substrate orientation, interdomain correlated motions, HAT transition state stabilization, reaction energetics, electron transfer mechanism, and alterations in the intrinsic electric field of PHF8. Importantly, SCS and LR variations decrease key motions of α9-α12 of the JmjC domain toward the Fe(IV)-center that are associated with tighter binding of the H3_1-14K9me2 substrate. SCS and LR residues alter the intrinsic electric field of the enzyme along the reaction coordinate and change the individual energetic contributions of residues toward TS stabilization. The overall results suggest that DCCA can indeed identify non-active-site residues relevant for catalysis. The substitutions of such dynamically correlated residues might be used as a tool to tune HAT in non-heme Fe(II)- and 2OG-dependent enzymes.

Spectroscopic studies reveal details of substrate-induced conformational changes distant from the active site in isopenicillin N synthase

Isopenicillin N synthase (IPNS) catalyzes formation of the β-lactam and thiazolidine rings of isopenicillin N from its linear tripeptide l-δ-(α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) substrate in an iron- and dioxygen (O₂)-dependent four-electron oxidation without precedent in current synthetic chemistry. Recent X-ray free-electron laser studies including time-resolved serial femtosecond crystallography show that binding of O₂ to the IPNS–Fe(II)–ACV complex induces unexpected conformational changes in α-helices on the surface of IPNS, in particular in α3 and α10. However, how substrate binding leads to conformational changes away from the active site is unknown. Here, using detailed ¹⁹F NMR and electron paramagnetic resonance experiments with labeled IPNS variants, we investigated motions in α3 and α10 induced by binding of ferrous iron, ACV, and the O₂ analog nitric oxide, using the less mobile α6 for comparison. ¹⁹F NMR studies were carried out on singly and doubly labeled α3, α6, and α10 variants at different temperatures. In addition, double electron–electron resonance electron paramagnetic resonance analysis was carried out on doubly spin-labeled variants. The combined spectroscopic and crystallographic results reveal that substantial conformational changes in regions of IPNS including α3 and α10 are induced by binding of ACV and nitric oxide. Since IPNS is a member of the structural superfamily of 2-oxoglutarate-dependent oxygenases and related enzymes, related conformational changes may be of general importance in nonheme oxygenase catalysis.

Use of Noncanonical Tyrosine Analogues to Probe Control of Radical Intermediates during Endoperoxide Installation by Verruculogen Synthase (FtmOx1)

Important bioactive natural products, including prostaglandin H2 and artemisinin, contain reactive endoperoxides. Known enzymatic pathways for endoperoxide installation require multiple hydrogen-atom transfers (HATs). For example, iron(II)- and 2-oxoglutarate-dependent verruculogen synthase (FtmOx1; EC 1.14.11.38) mediates HAT from aliphatic C21 of fumitremorgin B, capture of O2 by the C21 radical (C21•), addition of the peroxyl radical (C21-O-O•) to olefinic C27, and HAT to the resultant C26•. Recent studies proposed conflicting roles for FtmOx1 tyrosine residues, Tyr224 and Tyr68, in the HATs from C21 and to C26•. Here, analysis of variant proteins bearing a ring-halogenated tyrosine or (amino)phenylalanine in place of either residue establishes that Tyr68 is the hydrogen donor to C26•, while Tyr224 has no essential role. The radicals that accumulate rapidly in FtmOx1 variants bearing a HAT-competent tyrosine analog at position 68 exhibit hypsochromically shifted absorption and, in cases of fluorine substitution, 19F-coupled electron-paramagnetic-resonance (EPR) spectra. By contrast, functional Tyr224-substituted variants generate radicals with unaltered light-absorption and EPR signatures as they produce verruculogen. The alternative major product of the Tyr68Phe variant, which forms competitively with verruculogen also in wild-type FtmOx1 in 2H2O and in the variant with the less readily oxidized 2,3-F2Tyr at position 68, is identified by mass spectrometry and isotopic labeling as the 26-hydroxy-21,27-endoperoxide compound formed after capture of another equivalent of O2 by the longer lived C26•. The results highlight the considerable chemical challenges the enzyme must navigate in averting both oxygen rebound and a second O2 coupling to obtain verruculogen selectively over other possible products.

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.

Structural Insight into the Catalytic Mechanism of the Endoperoxide Synthase FtmOx1

The 2-oxoglutarate (2OG)-dependent non-heme enzyme FtmOx1 catalyzes the endoperoxide biosynthesis of verruculogen. Although several mechanistic studies have been carried out, the catalytic mechanism of FtmOx1 is not well determined owing to the lack of a reliable complex structure of FtmOx1 with fumitremorgin B. Herein we provide the X-ray crystal structure of the ternary complex FtmOx1⋅2OG⋅fumitremorgin B at a resolution of 1.22 Å. Our structures show that the binding of fumitremorgin B induces significant compression of the active pocket and that Y68 is in close proximity to C26 of the substrate. Further MD simulation and QM/MM calculations support a CarC-like mechanism, in which Y68 acts as the H atom donor for quenching the C26-centered substrate radical. Our results are consistent with all available experimental data and highlight the importance of accurate complex structures in the mechanistic study of enzymatic catalysis.

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.

Can the isonitrile biosynthesis enzyme ScoE assist with the biosynthesis of isonitrile groups in drug molecules? A computational study

Computational studies show that the isonitrile synthesizing enzyme ScoE can catalyse the conversion of γ-Gly substituents in substrates to isonitrile. This enables efficient isonitrile substitution into target molecules such as axisonitrile-1.

Effects of the Nature of the Metal Ion, Protein and Substrate on the Catalytic Center in Matrix Metalloproteinase-1: Insights from Multilevel MD, QM/MM and QM Studies

Matrix metalloproteinase-1 (MMP-1) is a Zn(II) dependent endopeptidase involved in the degradation of collagen, the most abundant structural protein in the extracellular matrix of connective tissues and the human body. Herein we performed a multilevel computational analysis including molecular dynamics (MD), combined quantum mechanics/molecular mechanics (QM/MM), and quantum mechanics (QM) calculations to characterize the structure and geometry of the catalytic Zn(II) within the MMP-1 protein environment in comparison to crystallographic and spectroscopic data. The substrate's removal fine-tuned impact on the conformational dynamics and geometry of the catalytic Zn(II) center was also explored. Finally, the study examined the effect of substituting catalytic Zn(II) by Co(II) on the overall structure and dynamics of the MMP-1 THP complex and specifically on the geometry of the catalytic metal center. Overall our QM/MM and QM studies were in good agreement with the MM description of the Zn(II) centers in the MD simulations.

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.

Nonadiabatic Effects in the Molecular Oxidation of Subnanometric Cu5 Clusters

The electronic structure of subnanometric clusters, far off the bulk regime, is still dominated by molecular characteristics. The spatial arrangement of the notoriously undercoordinated metal atoms is strongly coupled to the electronic properties of the system, which makes this class of materials particularly interesting for applications including luminescence, sensing, bioimaging, theranostics, energy conversion, catalysis, and photocatalysis. Opposing a common rule of thumb that assumes an increasing chemical reactivity with smaller cluster size, Cu₅ clusters have proven to be exceptionally resistant to irreversible oxidation, i.e., the dissociative chemisorption of molecular oxygen. Besides providing reasons for this behavior in the case of heavy loading with molecular oxygen, we investigate the competition between physisorption and molecular chemisorption from the perspective of nonadiabatic effects. Landau–Zener theory is applied to the Cu₅(O₂)₃ complex to estimate the probability for a switching between the electronic states correlating the neutral O₂ + Cu₅(O₂)₂ and the ionic O₂^– + (Cu₅(O₂)₂)⁺ fragments in a diabatic representation. Our work demonstrates the involvement of strong nonadiabatic effects in the associated charge transfer process, which might be a common motive in reactions involving subnanometric metal structures.

Active site characterization and activity of the human aspartyl (asparaginyl) β-hydroxylase

Human aspartyl/asparaginyl beta-hydroxylase (HAAH) is a member of the superfamily of nonheme Fe2+/α-ketoglutarate (αKG) dependent oxygenase enzymes with a noncanonical active site. HAAH hydroxylates epidermal growth factor (EGF) like domains to form the β-hydroxylated product from substrate asparagine or aspartic acid and has been suggested to have a negative impact in a variety of cancers. In addition to iron, HAAH also binds divalent calcium, although the role of the latter is not understood. Herein, the metal binding chemistry and influence on enzyme stability and activity have been evaluated by a combined biochemical and biophysical approach. Metal binding parameters for the HAAH active site were determined by use of isothermal titration calorimetry, demonstrating a high-affinity regulatory binding site for Ca2+ in the catalytic domain in addition to the catalytic Fe2+ cofactor. We have analyzed various active site derivatives, utilizing LC-MS and a new HPLC technique to determine the role of metal binding and the second coordination sphere in enzyme activity, discovering a previously unreported residue as vital for HAAH turnover. This analysis of the in vitro biochemical function of HAAH furthers the understanding of its importance to cellular biochemistry and metabolic pathways.

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.

N 6-methyladenosine binding induces a metal-centered rearrangement that activates the human RNA demethylase Alkbh5

Alkbh5 catalyzes demethylation of the N ⁶-methyladenosine (m⁶A), an epigenetic mark that controls several physiological processes including carcinogenesis and stem cell differentiation. The activity of Alkbh5 comprises two coupled reactions. The first reaction involves decarboxylation of α-ketoglutarate (αKG) and formation of a Fe⁴⁺═O species. This oxyferryl intermediate oxidizes the m⁶A to reestablish the canonical base. Despite coupling between the two reactions being required for the correct Alkbh5 functioning, the mechanisms linking dioxygen activation to m⁶A binding are not fully understood. Here, we use solution NMR to investigate the structure and dynamics of apo and holo Alkbh5. We show that binding of m⁶A to Alkbh5 induces a metal-centered rearrangement of αKG that increases the exposed area of the metal, making it available for binding O₂ Our study reveals the molecular mechanisms underlying activation of Alkbh5, therefore opening new perspectives for the design of novel strategies to control gene expression and cancer progression.

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

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

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.

Comparative molecular dynamics simulations identify a salt-sensitive loop responsible for the halotolerant activity of GH5 cellulases

Halotolerant glycoside hydrolases (GH) have broad application potentials in biorefinery industries. Elucidating the structure-activity relationship underlying the halotolerant catalysis is essential to design superior biocatalysts. Here, we performed molecular dynamics simulations to investigate the structural dynamics of two GH5 cellulases, namely the halotolerant Cel5R and non-halotolerant TfCel5A. Through characterizing the physical properties at different salt concentrations, the results revealed that the overall structures of Cel5R and TfCel5A were marginally affected by the increase in salt concentrations. However, a salt-sensitive loop was identified from both Cel5R and TfCel5A based on its significantly increased flexibility at high salt concentrations. Importantly, compared to TfCel5A the salt-sensitive loop of Cel5R engaged more sodium ions and water molecules around the active site of the enzyme. Besides, the unique residue motif of the salt-sensitive loop in Cel5R formed more intramolecular hydrogen bonds, stabilizing the active architecture of Cel5R at high salt concentrations. Collectively, the structural and dynamic differences may contribute to the various catalytic halotolerance of Cel5R and TfCel5A. These findings provide mechanistic insight into the halotolerant catalysis and will guide the ration design of GH5 cellulases with improved catalytic properties.Communicated by Ramaswamy H. Samy.

Multi- and single-reference methods for the analysis of multi-state peroxidation of enolates

In spite of being spin-forbidden, some enzymes are capable of catalyzing the incorporation of O₂(Σg-3) to organic substrates without needing any cofactor. It has been established that the process followed by these enzymes starts with the deprotonation of the substrate forming an enolate. In a second stage, the peroxidation of the enolate formation occurs, a process in which the system changes its spin multiplicity from a triplet state to a singlet state. In this article, we study the addition of O₂ to enolates using state-of-the-art multi-reference and single-reference methods. Our results confirm that intersystem crossing is promoted by stabilization of the singlet state along the reaction path. When multi-reference methods are used, large active spaces are required, and in this situation, semistochastic heat-bath configuration interaction emerges as a powerful method to study these multi-configurational systems and is in good agreement with PNO-LCCSD(T) when the system is well-represented by a single-configuration.

Vitamin C-Sources, Physiological Role, Kinetics, Deficiency, Use, Toxicity, and Determination

Vitamin C (L-ascorbic acid) has been known as an antioxidant for most people. However, its physiological role is much larger and encompasses very different processes ranging from facilitation of iron absorption through involvement in hormones and carnitine synthesis for important roles in epigenetic processes. Contrarily, high doses act as a pro-oxidant than an anti-oxidant. This may also be the reason why plasma levels are meticulously regulated on the level of absorption and excretion in the kidney. Interestingly, most cells contain vitamin C in millimolar concentrations, which is much higher than its plasma concentrations, and compared to other vitamins. The role of vitamin C is well demonstrated by miscellaneous symptoms of its absence-scurvy. The only clinically well-documented indication for vitamin C is scurvy. The effects of vitamin C administration on cancer, cardiovascular diseases, and infections are rather minor or even debatable in the general population. Vitamin C is relatively safe, but caution should be given to the administration of high doses, which can cause overt side effects in some susceptible patients (e.g., oxalate renal stones). Lastly, analytical methods for its determination with advantages and pitfalls are also discussed in this review.

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.

DpgC-Catalyzed Peroxidation of 3,5-Dihydroxyphenylacetyl-CoA (DPA-CoA): Insights into the Spin-Forbidden Transition and Charge Transfer Mechanisms*

Despite being a very strong oxidizing agent, most organic molecules are not oxidized in the presence of O₂ at room temperature because O₂ is a diradical whereas most organic molecules are closed-shell. Oxidation then requires a change in the spin state of the system, which is forbidden according to non-relativistic quantum theory. To overcome this limitation, oxygenases usually rely on metal or redox cofactors to catalyze the incorporation of, at least, one oxygen atom into an organic substrate. However, some oxygenases do not require any cofactor, and the detailed mechanism followed by these enzymes remains elusive. To fill this gap, here the mechanism for the enzymatic cofactor-independent oxidation of 3,5-dihydroxyphenylacetyl-CoA (DPA-CoA) is studied by combining multireference calculations on a model system with QM/MM calculations. Our results reveal that intersystem crossing takes place without requiring the previous protonation of molecular oxygen. The characterization of the electronic states reveals that electron transfer is concomitant with the triplet-singlet transition. The enzyme plays a passive role in promoting the intersystem crossing, although spontaneous reorganization of the water wire connecting the active site with the bulk presets the substrate for subsequent chemical transformations. The results show that the stabilization of the singlet radical-pair between dioxygen and enolate is enough to promote spin-forbidden reaction without the need for neither metal cofactors nor basic residues in the active site.

Production of ROS by Gallic Acid Activates KDM2A to Reduce rRNA Transcription

Metformin, which is suggested to have anti-cancer effects, activates KDM2A to reduce rRNA transcription and proliferation of cancer cells. Thus, the specific activation of KDM2A may be applicable to the treatment of cancers. In this study, we screened a food-additive compound library to identify compounds that control cell proliferation. We found that gallic acid activated KDM2A to reduce rRNA transcription and cell proliferation in breast cancer MCF-7 cells. Gallic acid accelerated ROS production and activated AMPK. When ROS production or AMPK activity was inhibited, gallic acid did not activate KDM2A. These results suggest that both ROS production and AMPK activation are required for activation of KDM2A by gallic acid. Gallic acid did not reduce the succinate level, which was required for KDM2A activation by metformin. Metformin did not elevate ROS production. These results suggest that the activation of KDM2A by gallic acid includes mechanisms distinct from those by metformin. Therefore, signals from multiple intracellular conditions converge in KDM2A to control rRNA transcription. Gallic acid did not induce KDM2A-dependent anti-proliferation activity in non-tumorigenic MCF10A cells. These results suggest that the mechanism of KDM2A activation by gallic acid may be applicable to the treatment of breast cancers.

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.

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.