Conformational Dynamics Underlies Different Functions of Human KDM7 Histone Demethylases

Epigenetic control over the cell-intrinsic immune response antagonizes self-renewal in acute myeloid leukemia

ABSTRACT

Epigenetic modulation of the cell-intrinsic immune response holds promise as a therapeutic approach for leukemia. However, current strategies designed for transcriptional activation of endogenous transposons and subsequent interferon type-I (IFN-I) response, show limited clinical efficacy. Histone lysine methylation is an epigenetic signature in IFN-I response associated with suppression of IFN-I and IFN-stimulated genes, suggesting histone demethylation as key mechanism of reactivation. In this study, we unveil the histone demethylase PHF8 as a direct initiator and regulator of cell-intrinsic immune response in acute myeloid leukemia (AML). Site-specific phosphorylation of PHF8 orchestrates epigenetic changes that upregulate cytosolic RNA sensors, particularly the TRIM25-RIG-I-IFIT5 axis, thereby triggering the cellular IFN-I response-differentiation-apoptosis network. This signaling cascade largely counteracts differentiation block and growth of human AML cells across various disease subtypes in vitro and in vivo. Through proteome analysis of over 200 primary AML bone marrow samples, we identify a distinct PHF8/IFN-I signature in half of the patient population, without significant associations with known clinically or genetically defined AML subgroups. This profile was absent in healthy CD34+ hematopoietic progenitor cells, suggesting therapeutic applicability in a large fraction of patients with AML. Pharmacological support of PHF8 phosphorylation significantly impairs the growth in samples from patients with primary AML. These findings provide novel opportunities for harnessing the cell-intrinsic immune response in the development of immunotherapeutic strategies against AML.

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.

Lysine-specific demethylase 7A (KDM7A): A potential target for disease therapy

Histone demethylation is a kind of epigenetic modification mediated by a variety of enzymes and participates in regulating multiple physiological and pathological events. Lysine-specific demethylase 7A is a kind of α-ketoglutarate- and Fe(II)-dependent demethylase belonging to the PHF2/8 subfamily of the JmjC demethylases. KDM7A is mainly localized in the nucleus and contributes to transcriptional activation via removing mono- and di-methyl groups from the lysine residues 9 and 27 of Histone H3. Mounting studies support that KDM7A is not only necessary for normal embryonic, neural, and skeletal development, but also associated with cancer, inflammation, osteoporosis, and other diseases. Herein, the structure of KDM7A is described by comparing the similarities and differences of its amino acid sequences of KDM7A and other Histone demethylases; the functions of KDM7A in homeostasis and dyshomeostasis are summarized via documenting its content and related signaling; the currently known KDM7A-specific inhibitors and their structural relationship are listed based on their structure optimization and pharmacological activities; and the challenges and opportunities in exploring functions and developing targeted agents of KDM7A are also prospected via presenting encountered problems and potential solutions, which will provide an insight in functional exploration and drug discovery for KDM7A-related diseases.

Cyclic peptides target the aromatic cage of a PHD-finger reader domain to modulate epigenetic protein function

Plant homeodomain fingers (PHD-fingers) are a family of reader domains that can recruit epigenetic proteins to specific histone modification sites. Many PHD-fingers recognise methylated lysines on histone tails and play crucial roles in transcriptional regulation, with their dysregulation linked to various human diseases. Despite their biological importance, chemical inhibitors for targeting PHD-fingers are very limited. Here we report a potent and selective de novo cyclic peptide inhibitor (OC9) targeting the Nε-trimethyllysine-binding PHD-fingers of the KDM7 histone demethylases, developed using mRNA display. OC9 disrupts PHD-finger interaction with histone H3K4me3 by engaging the Nε-methyllysine-binding aromatic cage through a valine, revealing a new non-lysine recognition motif for the PHD-fingers that does not require cation-π interaction. PHD-finger inhibition by OC9 impacted JmjC-domain mediated demethylase activity at H3K9me2, leading to inhibition of KDM7B (PHF8) but stimulation of KDM7A (KIAA1718), representing a new approach for selective allosteric modulation of demethylase activity. Chemoproteomic analysis showed selective engagement of OC9 with KDM7s in T cell lymphoblastic lymphoma SUP T1 cells. Our results highlight the utility of mRNA-display derived cyclic peptides for targeting challenging epigenetic reader proteins to probe their biology, and the broader potential of this approach for targeting protein-protein interactions.

Targeting epigenetic regulators to overcome drug resistance in cancers

Drug resistance is mainly responsible for cancer recurrence and poor prognosis. Epigenetic regulation is a heritable change in gene expressions independent of nucleotide sequence changes. As the common epigenetic regulation mechanisms, DNA methylation, histone modification, and non-coding RNA regulation have been well studied. Increasing evidence has shown that aberrant epigenetic regulations contribute to tumor resistance. Therefore, targeting epigenetic regulators represents an effective strategy to reverse drug resistance. In this review, we mainly summarize the roles of epigenetic regulation in tumor resistance. In addition, as the essential factors for epigenetic modifications, histone demethylases mediate the histone or genomic DNA modifications. Herein, we comprehensively describe the functions of the histone demethylase family including the lysine-specific demethylase family, the Jumonji C-domain-containing demethylase family, and the histone arginine demethylase family, and fully discuss their regulatory mechanisms related to cancer drug resistance. In addition, therapeutic strategies, including small-molecule inhibitors and small interfering RNA targeting histone demethylases to overcome drug resistance, are also described.

Discovery of hit compounds for methyl-lysine reader proteins from a target class DNA-encoded library

Methyl-lysine (Kme) reader domains are prevalent in chromatin regulatory proteins which bind post-translational modification sites to recruit repressive and activating factors; therefore, these proteins play crucial roles in cellular signaling and epigenetic regulation. Proteins that contain Kme domains are implicated in various diseases, including cancer, making them attractive therapeutic targets for drug and chemical probe discovery. Herein, we report on expanding the utility of a previously reported, Kme-focused DNA-encoded library (DEL), UNCDEL003, as a screening tool for hit discovery through the specific targeting of Kme reader proteins. As an efficient method for library generation, focused DELs are designed based on structural and functional features of a specific class of proteins with the intent of novel hit discovery. To broadly assess the applicability of our library, UNCDEL003 was screened against five diverse Kme reader protein domains (53BP1 TTD, KDM7B JmjC-PHD, CDYL2 CD, CBX2 CD, and LEDGF PWWP) with varying structures and functions. From these screening efforts, we identified hit compounds which contain unique chemical scaffolds distinct from previously reported ligands. The selected hit compounds were synthesized off-DNA and confirmed using primary and secondary assays and assessed for binding selectivity. Hit compounds from these efforts can serve as starting points for additional development and optimization into chemical probes to aid in further understanding the functionality of these therapeutically relevant proteins.

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.

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.

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.

Gallic Acid Derivatives Propyl Gallate and Epigallocatechin Gallate Reduce rRNA Transcription via Induction of KDM2A Activation

We previously reported that lysine-demethylase 2A (KDM2A), a Jumonji-C histone demethylase, is activated by gallic acid to reduce H3K36me2 levels in the rRNA gene promoter and consequently inhibit rRNA transcription and cell proliferation in the breast cancer cell line MCF-7. Gallic acid activates AMP-activated protein kinase (AMPK) and increases reactive oxygen species (ROS) production to activate KDM2A. Esters of gallic acid, propyl gallate (PG) and epigallocatechin gallate (EGCG), and other chemicals, reduce cancer cell proliferation. However, whether these compounds activate KDM2A has yet to be tested. In this study, we found that PG and EGCG decreased rRNA transcription and cell proliferation through KDM2A in MCF-7 cells. The activation of both AMPK and ROS production by PG or EGCG was required to activate KDM2A. Of note, while the elevation of ROS production by PG or EGCG was limited in time, it was sufficient to activate KDM2A. Importantly, the inhibition of rRNA transcription and cell proliferation by gallic acid, PG, or EGCG was specifically observed in MCF-7 cells, whereas it was not observed in non-tumorigenic MCF10A cells. Altogether, these results suggest that the derivatization of gallic acid may be used to obtain new compounds with anti-cancer activity.

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.

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.

Domain cross-talk in regulation of histone modifications: Molecular mechanisms and targeting opportunities

Functional cross-talk between the catalytic and reader domains in chromatin-modifying enzymes and protein complexes enable coordinated regulation of chromatin modification status, and consequently impacts chromatin-associated processes. ZZ domains are a recently identified class of chromatin readers that recognize the N-terminal region of histone H3 to direct and regulate acetylation activity of several histone acetylation complexes. Cross-talk between chromatin readers sensitive to methylation, and catalytic domains of methyltransferases and demethylases impacts substrate specificity, catalytic activity, and propagation of chromatin marks. Recently described allosteric ligands that target domain communication highlight the potential of domain cross-talk in the development of the next-generation of chromatin-directed therapeutics.

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.

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.

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.