Computations Reveal a Rich Mechanistic Variation of Demethylation of N-Methylated DNA/RNA Nucleotides by FTO

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.

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.

Flavin-N5OOH Functions as both a Powerful Nucleophile and a Base in the Superfamily of Flavoenzymes

Flavoenzymes can mediate a large variety of oxidation reactions through the activation of oxygen. However, the O2 activation chemistry of flavin enzymes is not yet fully exploited. Normally, the O2 activation occurs at the C4a site of the flavin cofactor, yielding the flavin C4a-(hydro)hydroperoxyl species in monooxygenases or oxidases. Using extensive MD simulations, QM/MM calculations and QM calculations, our studies reveal the formation of the common nucleophilic species, Flavin-N5OOH, in two distinct flavoenzymes (RutA and EncM). Our studies show that Flavin-N5OOH acts as a powerful nucleophile that promotes C-N cleavage of uracil in RutA, and a powerful base in the deprotonation of substrates in EncM. We reason that Flavin-N5OOH can be a common reactive species in the superfamily of flavoenzymes, which accomplish generally selective general base catalysis and C-X (X=N, S, Cl, O) cleavage reactions that are otherwise challenging with solvated hydroxide ion base. These results expand our understanding of the chemistry and catalysis of flavoenzymes.

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.

Understanding the initial events of the oxidative damage and protection mechanisms of the AA9 lytic polysaccharide monooxygenase family

Lytic polysaccharide monooxygenase (LPMO) is a new class of oxidoreductases that boosts polysaccharide degradation employing a copper active site. This boost may facilitate the cost-efficient production of biofuels and high-value chemicals from polysaccharides such as lignocellulose. Unfortunately, self-oxidation of the active site inactivates LPMOs. Other oxidoreductases employ hole-hopping mechanisms as protection against oxidative damage, but little is generally known about the details of these mechanisms. Herein, we employ highly accurate theoretical models based on density functional theory (DFT) molecular mechanics (MM) hybrids to understand the initial steps in LPMOs' protective measures against self-oxidation; we identify several intermediates recently proposed from experiment, and quantify which are important for protective hole-hopping pathways. Investigations on two different LPMOs show consistently that a tyrosine residue close to copper is crucial for protection: this explains recent experiments, showing that LPMOs without this tyrosine are more susceptible to self-oxidation.

Impact of the T296S mutation in P450 GcoA for aryl-O-demethylation: a QM/MM study

Lignin, a complex plant cell wall component, holds promise as a renewable aromatic carbon feedstock. p-Vanillin is a key product of lignin depolymerization and a precursor of protocatechuic acid (PCA) that has tremendous potential for biofuel production. While the GcoAB enzyme, native to Amycolatopsis sp., naturally catalyzes aryl-O-demethylation toward guaiacol, recent research introduced a single mutation, T296S, into the GcoAP450 enzyme, enabling it to catalyze aryl-O-demethylation of p-vanillin. This structural modification increases the efficiency of GcoAP450 for the natural substrate while being active for p-vanillin. This study reveals the increased flexibility of p-vanillin and its ability to adapt a favorable conformation by aligning the methoxy group in close proximity to Fe(IV) = O of Cpd I in the active site of the T296S variant. The QM/MM calculations in accordance with the experimental data validated that the rate-limiting step for the oxidation of p-vanillin is hydrogen atom abstraction and provided a detailed geometric structure of stationary and saddle points for the oxidation of p-vanillin.

Deciphering Deoxynybomycin Biosynthesis Reveals Fe(II)/α-Ketoglutarate-Dependent Dioxygenase-Catalyzed Oxazoline Ring Formation and Decomposition

The antibacterial agents deoxynybomycin (DNM) and nybomycin (NM) have a unique tetracyclic structure featuring an angularly fused 4-oxazoline ring. Here, we report the identification of key enzymes responsible for forming the 4-oxazoline ring in Embleya hyalina NBRC 13850 by comparative bioinformatics analysis of the biosynthetic gene clusters encoding structurally similar natural products DNM, deoxynyboquinone (DNQ), and diazaquinomycins (DAQs). The N-methyltransferase DnmS plays a crucial role in catalyzing the N-dimethylation of a tricyclic precursor prenybomycin to generate NM D; subsequently, the Fe(II)/α-ketoglutarate-dependent dioxygenase (Fe/αKGD) DnmT catalyzes the formation of a 4-oxazoline ring from NM D to produce DNM; finally, a second Fe/αKGD DnmU catalyzes the C-12 hydroxylation of DNM to yield NM. Strikingly, DnmT is shown to display unexpected functions to also catalyze the decomposition of the 4-oxazoline ring and the N-demethylation, thereby converting DNM back to prenybomycin, to putatively serve as a manner to control the intracellular yield of DNM. Structure modeling, site-directed mutagenesis, and quantum mechanics calculations provide mechanistic insights into the DnmT-catalyzed reactions. This work expands our understanding of the functional diversity of Fe/αKGDs in natural product biosynthesis.

Characterization and structure-based protein engineering of a regiospecific saponin acetyltransferase from Astragalus membranaceus

Acetylation contributes to the bioactivity of numerous medicinally important natural products. However, little is known about the acetylation on sugar moieties. Here we report a saponin acetyltransferase from Astragalus membranaceus. AmAT7-3 is discovered through a stepwise gene mining approach and characterized as the xylose C3'/C4'-O-acetyltransferse of astragaloside IV (1). To elucidate its catalytic mechanism, complex crystal structures of AmAT7-3/1 and AmAT7-3A310G/1 are obtained, which reveal a large active pocket decided by a specific sequence AADAG. Combining with QM/MM computation, the regiospecificity of AmAT7-3 is determined by sugar positioning modulated by surrounding amino acids including #A310 and #L290. Furthermore, a small mutant library is built using semi-rational design, where variants A310G and A310W are found to catalyze specific C3'-O and C4'-O acetylation, respectively. AmAT7-3 and its variants are also employed to acetylate other bioactive saponins. This work expands the understanding of saponin acetyltransferases, and provide efficient catalytic tools for saponin acetylation.

QM/MM and MM MD Simulations on Enzymatic Degradation of the Nerve Agent VR by Phosphotriesterase

V-type nerve agents are hardly degraded by phosphotriesterase (PTE). Interestingly, the PTE variant of BHR-73MNW can effectively improve the hydrolytic efficiency of VR, especially for its Sp-enantiomer. Here, the whole enzymatic degradation of both Sp and Rp enantiomers of VR by the wild-type PTE and its variant BHR-73MNW was investigated by quantum mechanics/molecular mechanics (QM/MM) calculations and MM molecular dynamics simulations. Present results indicate that the degradation of VR can be initiated by the nucleophilic attack of the bridging OH- and the zinc-bound water molecule. The QM/MM-predicted energy barriers for the hydrolytic process of Sp-VR are 19.8 kcal mol-1 by the variant with water as a nucleophile and 22.0 kcal mol-1 by the wild-type PTE with OH- as a nucleophile, and corresponding degraded products are bound to the dinuclear metal site in monodentate and bidentate coordination modes, respectively. The variant effectively increases the volume of the large pocket, allowing more water molecules to enter the active pocket and resulting in the improvement of the degradation efficiency of Sp-VR. The hydrolysis of Rp-VR is triggered only by the hydroxide with an energy span of 20.6 kcal mol-1 for the wild-type PTE and 20.7 kcal mol-1 for the variant BHR-73-MNW PTE. Such mechanistic insights into the stereoselective degradation of VR by PTE and the role of water may inspire further studies to improve the catalytic efficiency of PTE toward the detoxification of nerve agents.

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.

Endowing homodimeric carbamoyltransferase GdmN with iterative functions through structural characterization and mechanistic studies

Iterative enzymes, which catalyze sequential reactions, have the potential to improve the atom economy and diversity of industrial enzymatic processes. Redesigning one-step enzymes to be iterative biocatalysts could further enhance these processes. Carbamoyltransferases (CTases) catalyze carbamoylation, an important modification for the bioactivity of many secondary metabolites with pharmaceutical applications. To generate an iterative CTase, we determine the X-ray structure of GdmN, a one-step CTase involved in ansamycin biosynthesis. GdmN forms a face-to-face homodimer through unusual C-terminal domains, a previously unknown functional form for CTases. Structural determination of GdmN complexed with multiple intermediates elucidates the carbamoylation process and identifies key binding residues within a spacious substrate-binding pocket. Further structural and computational analyses enable multi-site enzyme engineering, resulting in an iterative CTase with the capacity for successive 7-O and 3-O carbamoylations. Our findings reveal a subclade of the CTase family and exemplify the potential of protein engineering for generating iterative 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.

Molecular Mechanism of the Mononuclear Copper Complex-Catalyzed Water Oxidation from Cluster-Continuum Model Calculations

Cluster-continuum model calculations were conducted to decipher the mechanism of water oxidation catalyzed by a mononuclear copper complex. Among various O-O bond formation mechanisms investigated in this study, the most favorable pathway involved the nucleophilic attack of OH^- onto the ^.+ L-Cu^II -OH^- intermediate. During such process, the initial binding of OH^- to the proximity of ^.+ L-Cu^II -OH^- would result in the spontaneous oxidation of OH^- , leading to OH⋅ radical and Cu^II -OH^- species. The further O-O coupling between OH⋅ radical and Cu^II -OH^- was associated with a barrier of 14.8 kcal mol^-1 , leading to the formation of H₂ O₂ intermediate. Notably, the formation of "Cu^III -O^.- " species, a widely proposed active species for O-O bond formation, was found to be thermodynamically unfavorable and could be bypassed during the catalytic reactions. On the basis the present calculations, a catalytic cycle of the mononuclear copper complex-catalyzed water oxidation was proposed.

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.

Crystal structure and catalytic mechanism of the MbnBC holoenzyme required for methanobactin biosynthesis

Methanobactins (Mbns) are a family of copper-binding peptides involved in copper uptake by methanotrophs, and are potential therapeutic agents for treating diseases characterized by disordered copper accumulation. Mbns are produced via modification of MbnA precursor peptides at cysteine residues catalyzed by the core biosynthetic machinery containing MbnB, an iron-dependent enzyme, and MbnC. However, mechanistic details underlying the catalysis of the MbnBC holoenzyme remain unclear. Here, we present crystal structures of MbnABC complexes from two distinct species, revealing that the leader peptide of the substrate MbnA binds MbnC for recruitment of the MbnBC holoenzyme, while the core peptide of MbnA resides in the catalytic cavity created by the MbnB-MbnC interaction which harbors a unique tri-iron cluster. Ligation of the substrate sulfhydryl group to the tri-iron center achieves a dioxygen-dependent reaction for oxazolone-thioamide installation. Structural analysis of the MbnABC complexes together with functional investigation of MbnB variants identified a conserved catalytic aspartate residue as a general base required for MbnBC-mediated MbnA modification. Together, our study reveals the similar architecture and function of MbnBC complexes from different species, demonstrating an evolutionarily conserved catalytic mechanism of the MbnBC holoenzymes.

Identification of Potent and Selective Inhibitors of Fat Mass Obesity-Associated Protein Using a Fragment-Merging Approach

Fat mass obesity-associated protein (FTO) is a DNA/RNA demethylase involved in the epigenetic regulation of various genes and is considered a therapeutic target for obesity, cancer, and neurological disorders. Here, we aimed to design novel FTO-selective inhibitors by merging fragments of previously reported FTO inhibitors. Among the synthesized analogues, compound 11b, which merges key fragments of Hz (3) and MA (4), inhibited FTO selectively over alkylation repair homologue 5 (ALKBH5), another DNA/RNA demethylase. Treatment of acute monocytic leukemia NOMO-1 cells with a prodrug of 11b decreased the viability of acute monocytic leukemia cells, increased the level of the FTO substrate N⁶-methyladenosine in mRNA, and induced upregulation of MYC and downregulation of RARA, which are FTO target genes. Thus, Hz (3)/MA (4) hybrid analogues represent an entry into a new class of FTO-selective inhibitors.

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

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

Mechanistic Insights into the P450 TleB-Catalyzed Unusual Intramolecular C-N Bond Formation Involved in the Biosynthesis of Indolactam V

Indolactam V, a known biosynthetic precursor of indolactam alkaloids, is the main pharmacophore of this family and exhibits potential protein kinase C activation. A key step in the biosynthesis of indolactam V is the formation of an indole-fused nine-membered lactam core by intramolecular C-N bond formation. In this work, we report a computational study of the unique cytochrome P450 TleB enzyme-catalyzed direct and selective C-H bond amination reaction that can generate indolactam V from the dipeptide N-methylvalyl-tryptophanol. By performing molecular dynamics simulations and quantum-mechanical/molecular-mechanical calculations, we revealed that the C-H bond amination involves one step of proton transfer from N1-H of the indole ring to the Fe^IV═O unit, one step of hydrogen abstraction of N13-H in the side chain of the substrate by the Fe^IV-OH unit, and diradical coupling, in which two conformational changes of the side chain of the substrate are necessary. In the enzyme-substrate complex of TleB, the N-H bond of the indole ring of the substrate forms a strong hydrogen bond with the Fe^IV═O unit in compound I, and the porphyrin radical cation accepts an electron from the substrate to form the closed-shell electronic configuration. Thus, compound I in the enzyme-substrate complex cannot be described as Fe^IV═O coupled to a porphyrin radical cation, which is different from those of other P450 enzymes. Besides, two stages of conformational changes of the side chains of the substrate may increase the relative energies of reaction intermediates by 10-12 kcal/mol. From the structure point of view, it is the rotatable long side chain of the substrate and the large flexible active pocket of TleB that make the intramolecular diradical coupling feasible. Our findings may provide useful information to further understand the Tleb-catalyzed intramolecular C-H bond amination and the other bio-catalyzed intramolecular diradical coupling.

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

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

Nitrogen heterocycles are key and prevalent motifs in drugs.

The role of the active site tyrosine in the mechanism of lytic polysaccharide monooxygenase

Catalytic breakdown of polysaccharides can be achieved more efficiently by means of the enzymes lytic polysaccharide monooxygenases (LPMOs). However, the LPMO mechanism has remained controversial, preventing full exploitation of their potential. One of the controversies has centered around an active site tyrosine, present in most LPMO classes. Recent investigations have for the first time obtained direct (spectroscopic) evidence for the possibility of chemical modification of this tyrosine. However, the spectroscopic features obtained in the different investigations are remarkably different, with absorption maximum at 420 and 490 nm, respectively. In this paper we use density functional theory (DFT) in a QM/MM formulation to reconcile these (apparently) conflicting results. By modeling the spectroscopy as well as the underlying reaction mechanism we can show how formation of two isomers (both involving deprotonation of tyrosine) explains the difference in the observed spectroscopic features. Both isomers have a [TyrO–Cu–OH]⁺ moiety with the OH in either the cis- or trans-position to a deprotonated tyrosine. Although the cis-[TyrO–Cu–OH]⁺ moiety is well positioned for oxidation of the substrate, preliminary calculations with the substrate reveal that the reactivity is at best moderate, making a protective role of tyrosine more likely.

With QM/MM, we investigate the mechanism of tyrosine deprotonation in lytic polysaccharide monooxygenases. Our results support deprotonation and our calculated UV-vis spectra show that two isomers must be formed to match recent experiments.

Mechanistic Insights into the Oxidative Rearrangement Catalyzed by the Unprecedented Dioxygenase ChaP Involved in Chartreusin Biosynthesis

ChaP is a non-heme iron-dependent dioxygenase belonging to the vicinal oxygen chelate (VOC) enzyme superfamily that catalyzes the final α-pyrone ring formation in the biosynthesis of chartreusin. In contrast to other common dioxygenases, for example, 2,3-catechol dioxygenase which uses the dioxygen molecule as the oxidant, ChaP requires the flavin-activated oxygen (O₂^2-) as the equivalent. Previous experiments showed that the ChaP-catalyzed ring rearrangement contains two successive C-C bond cleavages and one lactonization; however, the detailed reaction mechanism is unknown. In this work, on the basis of the recently obtained crystal structure of ChaP, the computational model was constructed and the catalytic mechanism of ChaP was explored by performing quantum mechanical/molecular mechanical (QM/MM) calculations. Our calculation results reveal that ChaP uses the proximal oxygen in iron-coordinated HOO^- to attack the carbonyl carbon of the substrate, whereas the previous proposal that Asp49 acts as a base to deprotonate the iron-coordinated HOO^- to generate O₂^2- is unlikely. In the first stage reaction, owing to the coordination of the substrate with iron, the substrate is activated by accepting an electron from iron and the resulting oxy-radical intermediate formed by O-O cleavage can easily trigger the ring rearrangement. In the final decarboxylation, the phenolic anion of the substrate cooperatively accepts the proton of iron-coordinated HOO^- to facilitate the attack of the distal oxygen, and the proton-coupled electron transfer (PCET) from the substrate to the Fe^IV═O plays a key role for the decarboxylation. These findings may provide useful information for understanding the ChaP-catalyzed oxidative rearrangement and other flavin-dependent non-heme dioxygenases.

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.

Mechanism of fatty acid decarboxylation catalyzed by a non-heme iron oxidase (UndA): a QM/MM study

UndA is a non-heme iron enzyme that was recognized to catalyze the decarboxylation of medium chain (C10-C14) fatty acids to produce trace amounts of 1-alkenes. Owing to the electron imbalance during the oxidative decarboxylation of the substrate and the reduction of O₂, only single turnover reactions were obtained in UndA in vitro assays. Unlike the general non-heme iron enzymes, the catalytic efficiency of UndA is quite low. According to the previous proposal, both Fe^III-OO˙^- and Fe^IV[double bond, length as m-dash]O complexes may abstract the β-H of fatty acids to trigger the oxidative decarboxylation reaction. Herein, on the basis of the crystal structures of UndA in complex with the substrate analogues, we constructed a series of computational models and performed quantum mechanics/molecular mechanics (QM/MM) calculations to explore the UndA-catalyzed decarboxylation using lauric acid as the substrate. Our calculation results reveal that only the Fe^III-OO˙^- complex can initiate the decarboxylation, and the substrate (lauric acid) should monodentately coordinate to the Fe center to facilitate the β-H abstraction. In addition, the monodentate coordination corresponds to higher relative energy than the bidentate mode, which may explain the low efficiency of UndA. It is also revealed that as long as the β-H is extracted by the Fe^III-OO˙^-, the decarboxylation of the substrate radical is quite easy, and an electron transfer from the substrate to the iron center is the prerequisite. For the Fe^IV[double bond, length as m-dash]O complex, since the β-H is far from the O_Fe atom and the angle of ∠Fe-O-H is 53.1°, the H-abstraction is calculated to be difficult.

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.

Identification of Flavin Mononucleotide as a Cell-Active Artificial N6 -Methyladenosine RNA Demethylase

N⁶ -Methyladenosine (m⁶ A) represents a common and highly dynamic modification in eukaryotic RNA that affects various cellular pathways. Natural dioxygenases such as FTO and ALKBH5 are enzymes that demethylate m⁶ A residues in mRNA. Herein, the first identification of a small-molecule modulator that functions as an artificial m⁶ A demethylase is reported. Flavin mononucleotide (FMN), the metabolite produced by riboflavin kinase, mediates substantial photochemical demethylation of m⁶ A residues of RNA in live cells. This study provides a new perspective to the understanding of demethylation of m⁶ A residues in mRNA and sheds light on the development of powerful small molecules as RNA demethylases and new probes for use in RNA biology.

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

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

Open-Source, Python-Based Redevelopment of the ChemShell Multiscale QM/MM Environment

ChemShell is a scriptable computational chemistry environment with an emphasis on multiscale simulation of complex systems using combined quantum mechanical and molecular mechanical (QM/MM) methods. Motivated by a scientific need to efficiently and accurately model chemical reactions on surfaces and within microporous solids on massively parallel computing systems, we present a major redevelopment of the ChemShell code, which provides a modern platform for advanced QM/MM embedding models. The new version of ChemShell has been re-engineered from the ground up with a new QM/MM driver module, an improved parallelization framework, new interfaces to high performance QM and MM programs, and a user interface written in the Python programming language. The redeveloped package is capable of performing QM/MM calculations on systems of significantly increased size, which we illustrate with benchmarks on zirconium dioxide nanoparticles of over 160000 atoms.

Mechanical insights into the oxidative cleavage of resveratrol catalyzed by dioxygenase NOV1 from Novosphingobium aromaticivorans: confirmation of dioxygenase mechanism by QM/MM calculations

QM/MM calculations confirm that the oxidative cleavage of resveratrol catalyzed by dioxygenase NOV1 follows the dioxygenase mechanism.

Chemo- and Regioselective Dihydroxylation of Benzene to Hydroquinone Enabled by Engineered Cytochrome P450 Monooxygenase

Hydroquinone (HQ) is produced commercially from benzene by multi-step Hock-type processes with equivalent amounts of acetone as side-product. We describe an efficient biocatalytic alternative using the cytochrome P450-BM3 monooxygenase. Since the wildtype enzyme does not accept benzene, a semi-rational protein engineering strategy was developed. Highly active mutants were obtained which transform benzene in a one-pot sequence first into phenol and then regioselectively into HQ without any overoxidation. A computational study shows that the chemoselective oxidation of phenol by the P450-BM3 variant A82F/A328F leads to the regioselective formation of an epoxide intermediate at the C3=C4 double bond, which departs from the binding pocket and then undergoes fragmentation in aqueous medium with exclusive formation of HQ. As a practical application, an E. coli designer cell system was constructed, which enables the cascade transformation of benzene into the natural product arbutin, which has anti-inflammatory and anti-bacterial activities.

Oxidative Rearrangement Mechanism of Pentalenolactone F Catalyzed by Cytochrome P450 CYP161C2 (PntM)

The CYP161C2 (PntM) from Streptomyces arenae is a member of the cytochrome P450 enzymes, which catalyzes the unusual oxidative rearrangement of pentalenolactone F (1) to the sesquiterpenoid antibiotic pentalenolactone (3). On the basis of the crystal structure of PntM bound with substrate, quantum mechanical/molecular mechanics (QM/MM) calculations have been performed to explore the detailed mechanism of PntM-catalyzed oxidative rearrangement. The conversion from pentalenolactone F (1) to pentalenolactone (3) involves the stereospecific removal of the H-1 si from 1, the syn-1,2-migration of the 2 si methyl group, and the antarafacial loss of H-3 re. The abstraction of H-1 si by Cpd I is calculated to be rate limiting with an energy barrier of 20.3 kcal/mol, which basically agrees with the estimated free energy barrier from experiments (18.6 kcal/mol). It is the unfavorable geometry of Fe-OH-C1 that blocks the oxygen rebound reaction, and the subsequent intramolecular syn-1,2-methyl migration is accompanied by an electron transfer from the substrate to the porphyrin ring via an Fe-OH group, generating the carbocation intermediate. Owing to the positive charge, the intermediate can easily lose a proton to form the final products. Our calculation results indicate that both the carboxyl group of porphyrin and Fe-OH can act as bases to accept the proton of the substrate. The target product pentalenolactone and the three isomeric byproducts correspond to four different modes of deprotonation.

On how the binding cavity of AsqJ dioxygenase controls the desaturation reaction regioselectivity: a QM/MM study

The Fe(II)/2-oxoglutarate-dependent dioxygenase AsqJ from Aspergillus nidulans catalyses two pivotal steps in the synthesis of quinolone antibiotic 4'-methoxyviridicatin, i.e., desaturation and epoxidation of a benzodiazepinedione. The previous experimental results signal that, during the desaturation reaction, hydrogen atom transfer (HAT) from the benzylic carbon atom (C10) is a more likely step to initiate the reaction than the alternative HAT from the ring moiety (C3 atom). To unravel the origins of this regioselectivity and to explain why the observed reaction is desaturation and not the "default" hydroxylation, we performed a QM/MM study on the reaction catalysed by AsqJ. Herein, we report results that complement the experimental findings and suggest that HAT at the C10 position is the preferred reaction due to favourable interactions between the substrate and the binding cavity that compensate for the relatively high intrinsic barrier associated with the process. For the resultant radical intermediate, the desaturation/hydroxylation selectivity is governed by electronic properties of the reactants, i.e., the energy gap between the orbital that hosts the unpaired electron and the sigma orbital for the C-H bond as well as the gap between the orbitals mixing in transition state structures for each elementary step. Regiospecificity of the AsqJ dehydrogenation reaction is dictated by substrate-protein interactions. 82 × 44 mm (300 × 300 dpi).

Mechanistic insights into dioxygen activation, oxygen atom exchange and substrate epoxidation by AsqJ dioxygenase from quantum mechanical/molecular mechanical calculations

Herein, we use in-protein quantum mechanical/molecular mechanical (QM/MM) calculations to elucidate the mechanism of dioxygen activation, oxygen atom exchange and substrate epoxidation processes by AsqJ, an Fe^II/α-ketoglutarate-dependent dioxygenase (α-KGD) using a 2-His-1-Asp facial triad. Our results demonstrated that the whole reaction proceeds through a quintet surface. The dioxygen activation by AsqJ leads to a quintet penta-coordinated Fe^IV-oxo species, which has a square pyramidal geometry with the oxo group trans to His134. This penta-coordinated Fe^IV-oxo species is not the reactive one in the substrate epoxidation reaction since its oxo group is pointing away from the target C[double bond, length as m-dash]C bond. Instead, it can undergo the oxo group isomerization followed by water binding or the water binding followed by oxygen atom exchange to form the reactive hexa-coordinated Fe^IV-oxo species with the oxo group trans to His211. The calculated parameters of Mössbauer spectra for this hexa-coordinated Fe^IV-oxo intermediate are in excellent agreement with the experimental values, suggesting that it is most likely the experimentally trapped species. The calculated energetics indicated that the rate-limiting step is the substrate C[double bond, length as m-dash]C bond activation. This work improves our understanding of the dioxygen activation by α-KGD and provides important structural information about the reactive Fe^IV-oxo species.

Insights into the unprecedented epoxidation mechanism of fumitremorgin B endoperoxidase (FtmOx1) from Aspergillus fumigatus by QM/MM calculations

QM/MM calculations indicate that the quintet of the Fe^IVO complex firstly abstracts the hydrogen from Tyr228 to initiate the reaction, then the generated Tyr228 radical extracts the hydrogen from C₂₁ to form the C₂₁ radical, which binds the second dioxygen to complete the epoxidation.

Cleavage mechanism of the aliphatic C–C bond catalyzed by 2,4′-dihydroxyacetophenone dioxygenase from Alcaligenes sp. 4HAP: a QM/MM study

Calculations suggest that the reactant complex may firstly undergo a triplet–quintet crossing to initiate the reaction and then the subsequent chemistry occurs on the multiple-states surfaces. The key C–C bond cleavage is accompanied by an insertion reaction of oxygen radical.

Visible-light-mediated oxidative demethylation of N6-methyl adenines

A novel solution for the chemical modulation of RNA epigenetics via biologically compatible photo-catalysis under visible light irradiation has been described.