Dioxygen Binding in the Active Site of Histone Demethylase JMJD2A and the Role of the Protein Environment

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.

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.

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.

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.

Reaction Mechanism of Histone Demethylation in αKG-dependent Non-Heme Iron Enzymes

Histone demethylases (KDMs) catalyze histone lysine demethylation, an important epigenetic process that controls gene expression in eukaryotes, and represent important cancer drug targets for cancer treatment. Demethylation of histone is comprised of sequential reaction steps including oxygen activation, decarboxylation, and demethylation. The initial oxygen binding and activation steps have been studied. However, the information on the complete catalytic reaction cycle is limited, which has impeded the structure-based design of inhibitors targeting KDMs. Here we report the mechanism of the complete reaction steps catalyzed by a representative nonheme iron αKG-dependent KDM, PHF8 using QM/MM approaches. The atomic-level understanding on the complete reaction mechanism of PHF8 would shed light on the structure-based design of selective inhibitors targeting KDMs to intervene in cancer epigenetics.

Conformational Dynamics Underlies Different Functions of Human KDM7 Histone Demethylases

The human KDM7 subfamily histone H3 Nϵ-methyl lysine demethylases PHF8 (KDM7B) and KIAA1718 (KDM7A) have different substrate selectivities and are linked to genetic diseases and cancer. We describe experimentally based computational studies revealing that flexibility of the region linking the PHD finger and JmjC domains in PHF8 and KIAA1718 regulates interdomain interactions, the nature of correlated motions, and ultimately H3 binding and demethylation site selectivity. F279S an X-linked mental retardation mutation in PHF8 is involved in correlated motions with the iron ligands and second sphere residues. The calculations reveal key roles of a flexible protein environment in productive formation of enzyme-substrate complexes and suggest targeting the flexible KDM7 linker region is of interest from a medicinal chemistry perspective.

Hydrogen-Bond-Dependent Conformational Switching: A Computational Challenge from Experimental Thermochemistry

We have compiled an experimental data set (SWITCH10) of equilibrium constants for a series of hydrogen-bond-dependent conformational switches. These organic molecules possess common functionalities and are representative in terms of size and composition of systems routinely studied computationally. They exist as two well-defined conformations which serve as a useful tool to benchmark computational estimates of experimental Gibbs energy differences. We examine the performance of HF theory and a variety of density functionals (B3LYP, B3LYP-D3, CAM-B3LYP, ωB97X-D, M06-2X) against these experimental benchmarks. Surprisingly, despite a strong similarity between the two switch conformations, the average errors (0.4-1.7 kcal·mol^-1) obtained across the data set for all methods are larger than obtained with HF calculations. B3LYP was found to outperform implicitly and explicitly dispersion-corrected functionals, with an average error smaller by 1 kcal·mol^-1. Unsystematic errors in the optimized structures were found to contribute to the relatively poor performance obtained, while quasi-rigid rotor harmonic oscillator thermal contributions are important in improving the accuracy of computed Gibbs energy differences. These results emphasize the challenge of quantitative accuracy in computing solution-phase thermochemistry for flexible systems and caution against the often used (but unstated) assumption of favorable error cancellation in comparing conformers or stereoisomers.

Density Functional Theory Study on the Demethylation Reaction between Methylamine, Dimethylamine, Trimethylamine, and Tamoxifen Catalyzed by a Fe(IV)-Oxo Porphyrin Complex

In this work, we studied computationally the N-demethylation reaction of methylamine, dimethylamine, and trimethylamine as archetypal examples of primary, secondary, and tertiary amines catalyzed by high-field low-spin Fe-containing enzymes such as cytochromes P450. Using DFT calculations, we found that the expected C-H hydroxylation process was achieved for trimethylamine. When dimethylamine and methylamine were studied, two different reaction mechanisms (C-H hydroxylation and a double hydrogen atom transfer) were computed to be energetically accessible and both are equally preferred. Both processes led to the formation of formaldehyde and the N-demethylated substrate. Finally, as an illustrative example, the relative contribution of the three primary oxidation routes of tamoxifen was rationalized through energetic barriers obtained from density functional calculations and docking experiments involving CYP3A4 and CYP2D6 isoforms. We found that the N-demethylation process was the intrinsically favored one, whereas other oxidation reactions required most likely preorganization imposed by the residues close to the active sites.

Multiscale Modelling of Lytic Polysaccharide Monooxygenases

Lytic polysaccharide monooxygenase (LPMO) enzymes have attracted considerable attention owing to their ability to enhance polysaccharide depolymerization, making them interesting with respect to production of biofuel from cellulose. LPMOs are metalloenzymes that contain a mononuclear copper active site, capable of activating dioxygen. However, many details of this activation are unclear. Some aspects of the mechanism have previously been investigated from a computational angle. Yet, either these studies have employed only molecular mechanics (MM), which are inaccurate for metal active sites, or they have described only the active site with quantum mechanics (QM) and neglected the effect of the protein. Here, we employ hybrid QM and MM (QM/MM) methods to investigate the first steps of the LPMO mechanism, which is reduction of Cu^II to Cu^I and the formation of a Cu^II-superoxide complex. In the latter complex, the superoxide can bind either in an equatorial or an axial position. For both steps, we obtain structures that are markedly different from previous suggestions, based on small QM-cluster calculations. Our calculations show that the equatorial isomer of the superoxide complex is over 60 kJ/mol more stable than the axial isomer because it is stabilized by interactions with a second-coordination-sphere glutamine residue, suggesting a possible role for this residue. The coordination of superoxide in this manner agrees with recent experimental suggestions.