Insight into wild-type and T1372E TET2-mediated 5hmC oxidation using ab initio QM/MM calculations

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.

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.

Theoretical Insights into N-Glycoside Bond Cleavage of 5-Carboxycytosine by Thymine DNA Glycosylase: A QM/MM Study

Thymine DNA glycosylase (TDG)-mediated excision of 5-formylcytosine and 5-carboxylcytosine (5-caC) is a critical step in active DNA demethylation. Herein, we employed a combined quantum mechanics/molecular mechanics approach to investigate the reaction mechanism of TDG-catalyzed N-glycosidic bond cleavage of 5-caC. The calculated results show that TDG-catalyzed 5-caC excision follows a concerted (SN2) mechanism in which glycosidic bond dissociation is coupled with nucleophile attack. Protonation of the 5-caC anion contributes to the cleavage of the N-glycoside bond, in which the N3-protonated zwitterion and imino tautomers are more favorable than carboxyl-protonated amino tautomers. This is consistent with the experimental data. Furthermore, our results reveal that the configuration rearrangement process of the protonated 5-caC would lower the stability of the N-glycoside bond and substantially reduce the barrier height for the subsequent C1'-N1 bond cleavage. This should be attributed to the smaller electrostatic repulsion between the leaving base and the negative phosphate group as a result of the structural rearrangement.

Allosteric Modulation of the YAP/TAZ-TEAD Interaction by Palmitoylation and Small-Molecule Inhibitors

The Hippo signaling pathway is a highly conserved signaling network that plays a central role in regulating cellular growth, proliferation, and organ size. This pathway consists of a kinase cascade that integrates various upstream signals to control the activation or inactivation of YAP/TAZ proteins. Phosphorylated YAP/TAZ is sequestered in the cytoplasm; however, when the Hippo pathway is deactivated, it translocates into the nucleus, where it associates with TEAD transcription factors. This partnership is instrumental in regulating the transcription of progrowth and antiapoptotic genes. Thus, in many cancers, aberrantly hyperactivated YAP/TAZ promotes oncogenesis by contributing to cancer cell proliferation, metastasis, and therapy resistance. Because YAP and TAZ exert their oncogenic effects by binding with TEAD, it is critical to understand this key interaction to develop cancer therapeutics. Previous research has indicated that TEAD undergoes autopalmitoylation at a conserved cysteine, and small molecules that inhibit TEAD palmitoylation disrupt effective YAP/TAZ binding. However, how exactly palmitoylation contributes to YAP/TAZ-TEAD interactions and how the TEAD palmitoylation inhibitors disrupt this interaction remains unknown. Utilizing molecular dynamics simulations, our investigation not only provides detailed atomistic insight into the YAP/TAZ-TEAD dynamics but also unveils that the inhibitor studied influences the binding of YAP and TAZ to TEAD in distinct manners. This discovery has significant implications for the design and deployment of future molecular interventions targeting this interaction.

Iterative oxidation by TET1 is required for reprogramming of imprinting control regions and patterning of mouse sperm hypomethylated regions

Ten-eleven translocation (TET) enzymes iteratively oxidize 5-methylcytosine (5mC) to generate 5-hydroxymethylcytosine (5hmC), 5-formylcytosine, and 5-carboxylcytosine to facilitate active genome demethylation. Whether these bases are required to promote replication-coupled dilution or activate base excision repair during mammalian germline reprogramming remains unresolved due to the inability to decouple TET activities. Here, we generated two mouse lines expressing catalytically inactive TET1 (Tet1-HxD) and TET1 that stalls oxidation at 5hmC (Tet1-V). Tet1 knockout and catalytic mutant primordial germ cells (PGCs) fail to erase methylation at select imprinting control regions and promoters of meiosis-associated genes, validating the requirement for the iterative oxidation of 5mC for complete germline reprogramming. TET1V and TET1HxD rescue most hypermethylation of Tet1-/- sperm, suggesting the role of TET1 beyond its oxidative capability. We additionally identify a broader class of hypermethylated regions in Tet1 mutant mouse sperm that depend on TET oxidation for reprogramming. Our study demonstrates the link between TET1-mediated germline reprogramming and sperm methylome patterning.

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.

Leveraging QM/MM and Molecular Dynamics Simulations to Decipher the Reaction Mechanism of the Cas9 HNH Domain to Investigate Off-Target Effects

The clustered regularly interspaced short palindromic repeats (CRISPR) technology is an RNA-guided targeted genome-editing tool using Cas family proteins. Two magnesium-dependent nuclease domains of the Cas9 enzyme, termed HNH and RuvC, are responsible for cleaving the target DNA (t-DNA) and nontarget DNA strands, respectively. The HNH domain is believed to determine the DNA cleavage activity of both endonuclease domains and is sensitive to complementary RNA-DNA base pairing. However, the underlying molecular mechanisms of CRISPR-Cas9, by which it rebukes or accepts mismatches, are poorly understood. Thus, investigation of the structure and dynamics of the catalytic state of Cas9 with either matched or mismatched t-DNA can provide insights into improving its specificity by reducing off-target cleavages. Here, we focus on a recently discovered catalytic-active form of the Streptococcus pyogenes Cas9 (SpCas9) and employ classical molecular dynamics and coupled quantum mechanics/molecular mechanics simulations to study two possible mechanisms of t-DNA cleavage reaction catalyzed by the HNH domain. Moreover, by designing a mismatched t-DNA structure called MM5 (C to G at the fifth position from the protospacer adjacent motif region), the impact of single-guide RNA (sgRNA) and t-DNA complementarity on the catalysis process was investigated. Based on these simulations, our calculated binding affinities, minimum energy paths, and analysis of catalytically important residues provide atomic-level details of the differences between matched and mismatched cleavage reactions. In addition, several residues exhibit significant differences in their catalytic roles for the two studied systems, including K253, K263, R820, K896, and K913.

A Chemical Model of a TET Enzyme for Selective Oxidation of Hydroxymethyl Cytosine to Formyl Cytosine

Methylation/demethylation of cytosines in DNA is central to epigenetics, which plays crucial roles in the regulation of about half of all human genes. Although the methylation mechanism, which downregulates gene expression, has been sufficiently decoded; the demethylation pathway, which upregulates gene expression, still holds questions to be answered. Demethylation of 5-methylcytosine by ten-eleven translocation (TET) enzymes yields understudied but epigenetically relevant intermediates, 5-hydroxymethyl (5-hmC), 5-formyl (5-fC), and 5-carboxyl (5-caC) cytosines. Here we report an iron complex, Fe^IIITAML (TAML = tetraamido macrocyclic ligand), which can facilitate selective oxidation of 5-hmC to its oxidative derivatives by forming a high-valent Fe-oxo intermediate in the presence of H₂O₂ under physiologically relevant conditions. Detailed HPLC analyses supported by a wide reaction condition optimization for the 5-hmC → 5-fC oxidation provides us with a chemical model of the TET enzyme. This study shines light on future efforts for a better understanding of the roles of 5-hmC and the TET enzyme mechanism and potentially novel therapeutic methods.

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

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

Computational Characterization of the Inhibition Mechanism of Xanthine Oxidoreductase by Topiroxostat

Xanthine oxidase (XO) is a member of the molybdopterin-containing enzyme family. It interconverts xanthine to uric acid as the last step of purine catabolism in the human body. The high uric acid concentration in the blood directly leads to human diseases like gout and hyperuricemia. Therefore, drugs that inhibit the biosynthesis of uric acid by human XO have been clinically used for many years to decrease the concentration of uric acid in the blood. In this study, the inhibition mechanism of XO and a new promising drug, topiroxostat (code: FYX-051), is investigated by employing molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. This drug has been reported to act as both a noncovalent and covalent inhibitor and undergoes a stepwise inhibition by all its hydroxylated metabolites, which include 2-hydroxy-FYX-051, dihydroxy-FYX-051, and trihydroxy-FYX-051. However, the detailed mechanism of inhibition of each metabolite remains elusive and can be useful for designing more effective drugs with similar inhibition functions. Hence, herein we present the computational investigation of the structural and dynamical effects of FYX-051 and the calculated reaction mechanism for all of the oxidation steps catalyzed by the molybdopterin center in the active site. Calculated results for the proposed reaction mechanisms for each metabolite's inhibition reaction in the enzyme's active site, binding affinities, and the noncovalent interactions with the surrounding amino acid residues are consistent with previously reported experimental findings. Analysis of the noncovalent interactions via energy decomposition analysis (EDA) and noncovalent interaction (NCI) techniques suggests that residues L648, K771, E802, R839, L873, R880, R912, F914, F1009, L1014, and A1079 can be used as key interacting residues for further hybrid-type inhibitor development.

ALKBH5-mediated m6 A demethylation of TIRAP mRNA promotes radiation-induced liver fibrosis and decreases radiosensitivity of hepatocellular carcinoma

BACKGROUND

Radiation-induced hepatic stellate cell (HSC) activation promotes radiation-induced liver fibrosis (RILF), a complication for hepatocellular carcinoma (HCC) radiotherapy. The demethylase alpha-ketoglutarate-dependent dioxygenase alkB homolog 5 (ALKBH5) decreases N6-methyladenylate methylation (m6 A) modification of RNA, while its role in regulating RILF pathogenesis and HCC radiosensitivity remains unknown.

METHODS

Methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA-sequencing (RNA-seq) were used to screen target genes regulated by ALKBH5. HSC with altered ALKBH5 expression was used to assess irradiation-induced HSC activation and the effect of HSC on recruitment and polarisation of monocytes. Key cytokines in medium from irradiated HSC-educated monocytes were identified by cytokine array detection. The effects of blocking ALKBH5 and key cytokines on RILF and HCC radiosensitivity were also evaluated.

RESULTS

Radiation-induced ALKBH5 expression in HSC mediated m6 A demethylation of toll-interleukin 1 receptor domain containing adaptor protein (TIRAP) mRNA and activated its downstream NF-κB and JNK/Smad2 pathways to promote HSC activation. Additionally, ALKBH5 regulated CCL5 secretion by irradiated HSC to promote monocyte recruitment and M2 macrophage polarisation. Notably, polarised monocytes secreted CCL20 to up-regulate ALKBH5 expression in HSC, and reduce HCC radiosensitivity by activating ALKBH5/TIRAP axis in HCC cells. ALKBH5 knockdown-combined CCR6 (CCL20 receptor) inhibitor significantly alleviated RILF and improved HCC radiosensitivity in mice. HCC patients with high ALKBH5 and TIRAP expression were prone to radiation-induced liver injury and poor tumour response to radiotherapy.

CONCLUSIONS

Collectively, irradiation up-regulates ALKBH5 in HSC to mediate monocyte recruitment and M2 polarisation and form positive feedback to promote RILF and reduce HCC radiosensitivity. The dual roles of ALKBH5 as a microenvironmental regulator and radiosensitisation target provide new ideas for RILF prevention and radiosensitisation of HCC.

Histone tail electrostatics modulate E2-E3 enzyme dynamics: a gateway to regulate ubiquitination machinery

BRCA1 (Breast Cancer-Associated Protein 1) is a human tumor suppressor that functions as an ubiquitin (Ub) ligase enzyme (E3) and plays a key role in genomic stability and DNA repair. Heterodimerization of BRCA1 with BARD1 (BRCA1-associated RING domain protein 1) is known to increase its Ub ligase activity and is important for its stability, and cooperative activation of UbcH5c (Ub conjugating enzyme (E2)). Recent studies demonstrate the importance of ubiquitination of the nucleosomal H2A C-terminal tail by BRCA1/BARD1-UbcH5c in which its mutations inhibit ubiquitination, predispose cells to chromosomal instability and greatly increase the likelihood of breast and ovarian cancer development. Due to the lack of molecular-level insight on the flexible and dis-ordered H2A C-tail, its ubiquitination mechanism by BRCA1/BARD1-UbcH5c and its function and relationship to cancer susceptibility remain elusive. Here, we use molecular dynamics simulations to provide molecular-level insights into the dynamics of the less-studied H2A C-tail and BRCA1/BARD1-UbcH5c on the nucleosome surface and their effect on ubiquitination. Our results precisely identify the key interactions and residues that trigger conformational transitions of BRCA1/BARD1-UbcH5c, and characterize the important role of histone electrostatics in their dynamics. We provide a mechanistic basis for the H2A C-tail lysine approach to UbcH5c and show the role of H2A C-tail and UbcH5c dynamics in lysine ubiquitination. Furthermore, our data demonstrate the potential for ubiquitination based on the lysine position of the C-tail. Altogether, the findings of this study provide unrevealed insights into the mechanism of H2A C-tail ubiquitination and help us understand the communication between Ub ligase/Ub conjugating enzymes (E3/E2) and nucleosome to regulate ubiquitination machinery, paving the way for the development of effective treatments for cancer and chronic pain.

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.

Punicalagin Attenuates Neuronal Apoptosis by Upregulating 5-Hydroxymethylcytosine in the Diabetic Mouse Brain

Punicalagin exerts neuroprotective activity by improving AMP-activated kinase (AMPK) and mitochondrial Krebs cycle. AMPK and Krebs cycle metabolites regulate 5-hydroxymethylcytosine (5hmC) via acting on ten-eleven translocation (TET) enzymes. Therefore, we hypothesized that punicalagin inhibits diabetes-related neuronal apoptosis by upregulating 5hmC in the diabetic mouse brain. C57BL/6J mice aged 8 weeks were randomly separated into five groups (n = 10), normal control (NC), diabetes mellitus (DM), resveratrol (RES), low-dose punicalagin (LPU), and high-dose punicalagin (HPU). Compared with other groups, the neuronal apoptosis rate was significantly higher and the 5hmC level of the cerebral cortex was significantly lower in the DM group. The levels of TET2 and P-AMPKα/AMPKα were significantly lower in the DM group than in both LPU and HPU groups. The ratio of (succinic acid + fumaric acid)/α-ketoglutarate was significantly higher in the DM group than in other groups. The present results suggest that punicalagin upregulates 5hmC via activating AMPK and maintaining Krebs cycle homeostasis, thus inhibiting neuronal apoptosis in the diabetic mouse brain.

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.

The pH-Dependence of the Hydration of 5-Formylcytosine: an Experimental and Theoretical Study

5-Formylcytosine is an important nucleobase in epigenetic regulation, whose hydrate form has been implicated in the formation of 5-carboxycytosine as well as oligonucleotide binding events. The hydrate content of 5-formylcytosine and its uracil derivative has now been quantified using a combination of NMR and mass spectroscopic measurements as well as theoretical studies. Small amounts of hydrate can be identified for the protonated form of 5-formylcytosine and for neutral 5-formyluracil. For neutral 5-formylcytosine, however, direct detection of the hydrate was not possible due to its very low abundance. This is in full agreement with theoretical estimates.

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

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

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

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

Catalytic Space Engineering as a Strategy to Activate C-H Oxidation on 5-Methylcytosine in Mammalian Genome

Conditional remodeling of enzyme catalysis is a formidable challenge in protein engineering. Herein, we have undertaken a unique active site engineering tactic to command catalytic outcomes. With ten-eleven translocation (TET) enzyme as a paradigm, we show that variants with an expanded active site significantly enhance multistep C-H oxidation in 5-methylcytosine (5mC), whereas a crowded cavity leads to a single-step catalytic apparatus. We further identify an evolutionarily conserved residue in the TET family with a remarkable catalysis-directing ability. The activating variant demonstrated its prowess to oxidize 5mC in chromosomal DNA for potentiating expression of genes including tumor suppressors.

Combining Evolutionary Conservation and Quantum Topological Analyses To Determine Quantum Mechanics Subsystems for Biomolecular Quantum Mechanics/Molecular Mechanics Simulations

Selection of residues and other molecular fragments for inclusion in the quantum mechanics (QM) region for QM/molecular mechanics (MM) simulations is an important step for these calculations. Here, we present an approach that combines protein sequence/structure evolution and electron localization function (ELF) analyses. The combination of these two analyses allows the determination of whether a residue needs to be included in the QM subsystem or can be represented by the MM environment. We have applied this approach on two systems previously investigated by QM/MM simulations, 4-oxalocrotonate tautomerase (4OT) and ten-eleven translocation-2 (TET2), that provide examples where fragments may or may not need to be included in the QM subsystem. Subsequently, we present the use of this approach to determine the appropriate QM subsystem to calculate the minimum energy path (MEP) for the reaction catalyzed by human DNA polymerase λ (Polλ) with a third cation in the active site. Our results suggest that the combination of protein evolutionary and ELF analyses provides insights into residue/molecular fragment selection for QM/MM simulations.

TET-TDG Active DNA Demethylation at CpG and Non-CpG Sites

In mammalian genomes, cytosine methylation occurs predominantly at CG (or CpG) dinucleotide contexts. As part of dynamic epigenetic regulation, 5-methylcytosine (mC) can be erased by active DNA demethylation, whereby ten-eleven translocation (TET) enzymes catalyze the stepwise oxidation of mC to 5-hydroxymethylcytosine (hmC), 5-formylcytosine (fC), and 5-carboxycytosine (caC), thymine DNA glycosylase (TDG) excises fC or caC, and base excision repair yields unmodified cytosine. In certain cell types, mC is also enriched at some non-CG (or CH) dinucleotides, however hmC is not. To provide biochemical context for the distribution of modified cytosines observed in biological systems, we systematically analyzed the activity of human TET2 and TDG for substrates in CG and CH contexts. We find that while TET2 oxidizes mC more efficiently in CG versus CH sites, this context preference can be diminished for hmC oxidation. Remarkably, TDG excision of fC and caC is only modestly dependent on CG context, contrasting its strong context dependence for thymine excision. We show that collaborative TET-TDG oxidation-excision activity is only marginally reduced for CA versus CG contexts. Our findings demonstrate that the TET-TDG-mediated demethylation pathway is not limited to CG sites and suggest a rationale for the depletion of hmCH in genomes rich in mCH.

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.

Comparison of DNA and RNA substrate effects on TET2 structure

Ten-eleven translocation (TET) enzymes can perform the stepwise oxidation of 5-methylcytosine (5mC) to 5-carboxylcytosine on both single-stranded (ss) and double-stranded (ds) DNA and RNA. It has been established that TET2 has a preference for ds DNA substrates, but it can catalyze the oxidation reaction on both ssDNA and RNA. The reasons for this substrate preference have been investigated for only a substrate 5mC ribonucleotide in a DNA strand, but not other nucleic acid configurations (Biochemistry58 (2019) 411). We performed molecular dynamics simulations on TET2 with various ss and ds substrates in order to better understand the structural and dynamical reasons for TET2's preference to act on ds DNA. Our simulations show that substrates that have a ribonucleotide experience several disruptions in their overall backbone shape, hydrogen bonding character, and non-bonded interactions. These differences appear to lead to the instability of ribonucleotide in the active site, and provide further rational for TET2's experimental behavior.

LICHEM 1.1: Recent Improvements and New Capabilities

The QM/MM method has become a useful tool to investigate various properties of complex systems. We previously introduced the layered interacting chemical models (LICHEM) package to enable QM/MM simulations with advanced potentials by combining various (unmodified) QM and MM codes ( J. Comp. Chem., 2016, 37, 1019). LICHEM provides several capabilities such as the ability to use polarizable force fields, such as AMOEBA, for the MM environment. Here, we describe an updated version of LICHEM (v1.1), which includes several new functionalities including a new method to account for long-range electrostatic effects in QM/MM (QM/MM-LREC), a new implementation for QM/MM with the Gaussian electrostatic model (GEM), and new capabilities for path optimizations using the quadratic string model (QSM) coupled with restrained MM environment optimization.

Polarizable ab initio QM/MM Study of the Reaction Mechanism of N-tert-Butyloxycarbonylation of Aniline in [EMIm][BF₄]

N-t e r t-butoxycarbonylation of amines in solution (water, organic solvents, or ionic liquids) is a common reaction for the preparation of drug molecules. To understand the reaction mechanism and the role of the solvent, quantum mechanical/molecular mechanical simulations using a polarizable multipolar force field with long⁻range electrostatic corrections were used to optimize the minimum energy paths (MEPs) associated with various possible reaction mechanisms employing the nudged elastic band (NEB) and the quadratic string method (QSM). The calculated reaction energies and energy barriers were compared with the corresponding gas-phase and dichloromethane results. Complementary Electron Localization Function (ELF)/NCI analyses provide insights on the critical structures along the MEP. The calculated results suggest the most likely path involves a sequential mechanism with the rate⁻limiting step corresponding to the nucleophilic attack of the aniline, followed by proton transfer and the release of CO 2 without the direct involvement of imidazolium cations as catalysts.