ONIOM(DFT:MM) Study of 2-Hydroxyethylphosphonate Dioxygenase: What Determines the Destinies of Different Substrates?

QM/MM Study Into the Mechanism of Oxidative C=C Double Bond Cleavage by Lignostilbene-α,β-Dioxygenase

The enzymatic biosynthesis of fragrance molecules from lignin fragments is an important reaction in biotechnology for the sustainable production of fine chemicals. In this work we investigated the biosynthesis of vanillin from lignostilbene by a nonheme iron dioxygenase using QM/MM and tested several suggested proposals via either an epoxide or dioxetane intermediate. Binding of dioxygen to the active site of the protein results in the formation of an iron(II)-superoxo species with lignostilbene cation radical. The dioxygenase mechanism starts with electrophilic attack of the terminal oxygen atom of the superoxo group on the central C=C bond of lignostilbene, and the second-coordination sphere effects in the substrate binding pocket guide the reaction towards dioxetane formation. The computed mechanism is rationalized with thermochemical cycles and valence bond schemes that explain the electron transfer processes during the reaction mechanism. Particularly, the polarity of the protein and the local electric field and dipole moments enable a facile electron transfer and an exergonic dioxetane formation pathway.

Which is the real oxidant in the competitive ligand self-hydroxylation and substrate oxidation, a biomimetic iron(II)-hydroperoxo species or an oxo-iron(IV)-hydroxy one?

Nonheme iron(II)-hydroperoxo species (FeII-(η2-OOH)) 1 and the concomitant oxo-iron(IV)-hydroxyl one 2 are proposed as the key intermediates of a large class of 2-oxoglutarate dependent dioxygenases (e.g., isopenicillin N synthase). Extensive...

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.

Fenton-Derived OH Radicals Enable the MPnS Enzyme to Convert 2-Hydroxyethylphosphonate to Methylphosphonate: Insights from Ab Initio QM/MM MD Simulations

The mechanism for dioxygen activation represents one of the core issues in metalloenzymes. In most cases, the activation of the O₂ molecule requires additional electrons from an external reducant. However, nonheme hydroxyethylphosphonate dioxygenase (HEPD) and methylphosphonate synthase (MPnS) are exceptional C-H oxygenases. Both enzymes do not utilize reductants, rather they employ directly iron(III)-superoxide species to initiate H-abstraction reactions and lead thereby to catalysis of the C-C cleavage in 2-hydroxyethylphosphonate (2-HEP). Using the recently characterized MPnS structure and QM(B3LYP)/MM-based metadynamics simulations, we deciphered the chemical mechanism for MPnS. Our simulations demonstrate O₂ activation in MPnS is mediated by an adjacent Lysine residue (Lys28) in the active site, leading to an unusual H ₂ O ₂ intermediate in the reductant-independent nonheme MPnS enzyme. Furthermore, the so-generated H ₂ O ₂ intermediate is subsequently employed in a Fenton-type reaction, leading to a locked •OH radical that spontaneously attaches to the substrate carbonyl group. Meanwhile, the proton from the Fe(III)-OH is shuttled back to the deprotonated Lys28, affording the Fe(IV)-oxo species that is identified by experiment in HEPD. Thus, our calculations demonstrate an unusual proton-shuttle mechanism for O ₂ activation in metalloenzymes.

On the encapsulation of Olsalazine by β-cyclodextrin: A DFT-based computational and spectroscopic investigations

In this work, the supramolecular host-guest interaction of the prodrug Olsalazine (OLZ) and β-Cyclodextrin (β-CD) was examined experimentally and computationally. Experimentally, employing the UV-Vis spectroscopic method in aqueous media at various pH's, results obtained using the Benesi-Hilderbrand approach demonstrated that OLZ can form supramolecular inclusion complex with β-CD with stoichiometric ratio of 1:1. Furthermore, these results revealed that the formation of OLZ: β-CD complexes exhibited insignificant pH dependency in the range 5-8 with an average binding constant (K_b) of approximately 1×10³M^-1. Computationally, geometry optimization of 1:1 OLZ: β-CD complexes was performed employing the ONIOM (DFT((ωB97XB)/6-31+G(d)),SQM(PM3)) approach. Obtained results demonstrated that OLZ: β-CD complex is stabilized by the formation of intermolecular hydrogen bonds with an average length of approximately 1.8Å. Additionally, the stability of OLZ: β-CD complex was demonstrated employing ADMP molecular dynamic simulations over a timeframe of 500fs. The molecularity of the supramolecular host-guest interaction between OLZ and β-CD is presented and interpreted in the essence of TD-DFT and molecular orbitals analyses.

Mechanistic Insights into a Stibene Cleavage Oxygenase NOV1 from Quantum Mechanical/Molecular Mechanical Calculations

NOV1, a stilbene cleavage oxygenase, catalyzes the cleavage of the central double bond of stilbenes to two phenolic aldehydes, using a 4-His Fe(II) center and dioxygen. Herein, we use in-protein quantum mechanical/molecular mechanical (QM/MM) calculations to elucidate the reaction mechanism of the central double bond cleavage of phytoalexin resveratrol by NOV1. Our results showed that the oxygen molecule prefers to bind to the iron center in a side-on fashion, as suggested from the experiment. The quintet Fe-O2 complex with the side-on superoxo antiferromagnetic coupled to the resveratrol radical is identified as the reactive oxygen species. The QM/MM results support the dioxygenase mechanism involving a dioxetane intermediate with a rate-limiting barrier of 10.0 kcal mol-1. The alternative pathway through an epoxide intermediate is ruled out due to a larger rate-limiting barrier (26.8 kcal mol-1). These findings provide important insight into the catalytic mechanism of carotenoid cleavage oxygenases and also the dioxygen activation of non-heme enzymes.

Mechanistic insights into ring cleavage of hydroquinone by PnpCD from quantum mechanical/molecular mechanical calculations

QM/MM calculations for ring cleavage of hydroquinone by PnpCD show that Asn258 loses coordination to the iron when the reaction begins. The first-sphere Glu262 can act as an acid–base catalyst to lower the rate-limiting barrier.

Mechanistic Insight on the Activity and Substrate Selectivity of Nonheme Iron Dioxygenases

Nonheme iron dioxygenases catalyze vital reactions for human health particularly related to aging processes. They are involved in the biosynthesis of amino acids, but also the biodegradation of toxic compounds. Typically they react with their substrate(s) through oxygen atom transfer, although often with the assistance of a co-substrate like α-ketoglutarate that is converted to succinate and CO₂ . Many reaction processes catalyzed by the nonheme iron dioxygenases are stereoselective or regiospecific and hence understanding the mechanism and protein involvement in the selectivity is important for the design of biotechnological applications of these enzymes. To this end, I will review recent work of our group on nonheme iron dioxygenases and include background information on their general structure and catalytic cycle. Examples of stereoselective and regiospecific reaction mechanisms we elucidated are for the AlkB repair enzyme, prolyl-4-hydroxylase and the ergothioneine biosynthesis enzyme. Finally, I cover an example where we bioengineered S-p-hydroxymandelate synthase into the R-p-hydroxymandelate synthase.

Sulfoxide Synthase versus Cysteine Dioxygenase Reactivity in a Nonheme Iron Enzyme

The sulfoxide synthase EgtB represents a unique family of nonheme iron enzymes that catalyze the formation of a C-S bond between N-α-trimethyl histidine and γ-glutamyl cysteine, which is the key step in the biosynthesis of ergothioneine, an important amino acid related to aging. A controversy has arisen regarding its catalytic mechanism related to the function of the active-site Tyr₃₇₇ residue. The biosynthesis of ergothioneine in EgtB shows structural similarities to cysteine dioxygenase which transfers two oxygen atoms to the thiolate group of cysteine. The question, therefore, is how do EgtB enzymes catalyze the C-S bond-formation reaction, while also preventing a dioxygenation of its cysteinate substrate? In this work we present a quantum mechanics/molecular mechanics study into the mechanism of sulfoxide synthase enzymes as compared to cysteine dioxygenase enzymes and present pathways for both reaction channels in EgtB. We show that EgtB contains a conserved tyrosine residue that reacts via proton-coupled electron transfer with the iron(III)-superoxo species and creates an iron(III)-hydroperoxo intermediate, thereby preventing the possible thiolate dioxygenation side reaction. The nucleophilic C-S bond-formation step happens subsequently concomitant to relay of the proton of the iron(II)-hydroperoxo back to Tyr₃₇₇. This is the rate-determining step in the reaction cycle and is followed by hydrogen-atom transfer from the CE1-H group of trimethyl histidine substrate to iron(II)-superoxo. In the final step, a quick and almost barrierless sulfoxidation leads to the sulfoxide product complexes. The work highlights a unique machinery and active-site setup of the enzyme that drives the sulfoxide synthase reaction.

O-H Activation by an Unexpected Ferryl Intermediate during Catalysis by 2-Hydroxyethylphosphonate Dioxygenase

Activation of O–H bonds by inorganic metal-oxo complexes has been documented, but no cognate enzymatic process is known. Our mechanistic analysis of 2-hydroxyethylphosphonate dioxygenase (HEPD), which cleaves the C1–C2 bond of its substrate to afford hydroxymethylphosphonate on the biosynthetic pathway to the commercial herbicide phosphinothricin, uncovered an example of such an O–H-bond-cleavage event. Stopped-flow UV–visible absorption and freeze-quench Mössbauer experiments identified a transient iron(IV)-oxo (ferryl) complex. Maximal accumulation of the intermediate required both the presence of deuterium in the substrate and, importantly, the use of ²H₂O as solvent. The ferryl complex forms and decays rapidly enough to be on the catalytic pathway. To account for these unanticipated results, a new mechanism that involves activation of an O–H bond by the ferryl complex is proposed. This mechanism accommodates all available data on the HEPD reaction.

Phosphonate Biochemistry

Organophosphonic acids are unique as natural products in terms of stability and mimicry. The C-P bond that defines these compounds resists hydrolytic cleavage, while the phosphonyl group is a versatile mimic of transition-states, intermediates, and primary metabolites. This versatility may explain why a variety of organisms have extensively explored the use organophosphonic acids as bioactive secondary metabolites. Several of these compounds, such as fosfomycin and bialaphos, figure prominently in human health and agriculture. The enzyme reactions that create these molecules are an interesting mix of chemistry that has been adopted from primary metabolism as well as those with no chemical precedent. Additionally, the phosphonate moiety represents a source of inorganic phosphate to microorganisms that live in environments that lack this nutrient; thus, unusual enzyme reactions have also evolved to cleave the C-P bond. This review is a comprehensive summary of the occurrence and function of organophosphonic acids natural products along with the mechanisms of the enzymes that synthesize and catabolize these molecules.

Go it alone: four-electron oxidations by mononuclear non-heme iron enzymes

This review discusses the current mechanistic understanding of a group of mononuclear non-heme iron-dependent enzymes that catalyze four-electron oxidation of their organic substrates without the use of any cofactors or cosubstrates. One set of enzymes acts on α-ketoacid-containing substrates, coupling decarboxylation to oxygen activation. This group includes 4-hydroxyphenylpyruvate dioxygenase, 4-hydroxymandelate synthase, and CloR involved in clorobiocin biosynthesis. A second set of enzymes acts on substrates containing a thiol group that coordinates to the iron. This group is comprised of isopenicillin N synthase, thiol dioxygenases, and enzymes involved in the biosynthesis of ergothioneine and ovothiol. The final group of enzymes includes HEPD and MPnS that both carry out the oxidative cleavage of the carbon-carbon bond of 2-hydroxyethylphosphonate but generate different products. Commonalities amongst many of these enzymes are discussed and include the initial substrate oxidation by a ferric-superoxo-intermediate and a second oxidation by a ferryl species.

Quantum Mechanics/Molecular Mechanics Modeling of Enzymatic Processes: Caveats and Breakthroughs

Nature has developed large groups of enzymatic catalysts with the aim to transfer substrates into useful products, which enables biosystems to perform all their natural functions. As such, all biochemical processes in our body (we drink, we eat, we breath, we sleep, etc.) are governed by enzymes. One of the problems associated with research on biocatalysts is that they react so fast that details of their reaction mechanisms cannot be obtained with experimental work. In recent years, major advances in computational hardware and software have been made and now large (bio)chemical systems can be studied using accurate computational techniques. One such technique is the quantum mechanics/molecular mechanics (QM/MM) technique, which has gained major momentum in recent years. Unfortunately, it is not a black-box method that is easily applied, but requires careful set-up procedures. In this work we give an overview on the technical difficulties and caveats of QM/MM and discuss work-protocols developed in our groups for running successful QM/MM calculations.

Implicit and explicit solvent effects on the selectivity of the cycloaddition reaction of cyclopentadiene and methyl acrylate; a theoretical study

The product selectivity of the Diels–Alder reaction between cyclopentadiene and methyl acrylate has been studied through a solvent effect analysis by implicit and explicit solvent models. Two paths for this reaction have been proposed and thermodynamic and kinetic parameters for each path were calculated by the B3LYP functional with 6-311++G(d,p) basis set at 298.15 K. In the explicit solvent model, hydrogen bonds were investigated by the QM/MM method. According to the results, the proportion of the exo product in the gas phase is more than for the endo one, which is in contrast to the solvent phase. HOMO–LUMO analyses in the gas phase and different solvents have been done and solvent effects on the energies of the HOMO–LUMO orbitals were investigated. According to the results, the band gap of the TSs was increased in the presence of the solvents. Finally, natural population analysis has been performed to calculate the electronic charges of some atoms that are involved in the centre of the reaction. The data obtained showed that different solvents, including polar and nonpolar solvents, have different effects on the electronic charges of the atoms. Topological analysis of the structures during the reaction by the quantum theory of atoms in molecules (QTAIM) method in the gas phase and different solvents confirmed a cyclic structure for the transition state.

Molecular Dynamics Simulation of Nitrobenzene Dioxygenase Using AMBER Force Field

Molecular dynamics simulation of the oxygenase component of nitrobenzene dioxygenase (NBDO) system, a member of the naphthalene family of Rieske nonheme iron dioxygenases, has been carried out using the AMBER force field combined with a new set of parameters for the description of the mononuclear nonheme iron center and iron-sulfur Rieske cluster. Simulation results provide information on the structure and dynamics of nitrobenzene dioxygenase in an aqueous environment and shed light on specific interactions that occur in its catalytic center. The results suggest that the architecture of the active site is stabilized by key hydrogen bonds, and Asn258 positions the substrate for oxidation. Analysis of protein-water interactions reveal the presence of a network of solvent molecules at the entrance to the active site, which could be of potential catalytic importance.

Applications of density functional theory to iron-containing molecules of bioinorganic interest

The past decades have seen an explosive growth in the application of density functional theory (DFT) methods to molecular systems that are of interest in a variety of scientific fields. Owing to its balanced accuracy and efficiency, DFT plays particularly useful roles in the theoretical investigation of large molecules. Even for biological molecules such as proteins, DFT finds application in the form of, e.g., hybrid quantum mechanics and molecular mechanics (QM/MM), in which DFT may be used as a QM method to describe a higher prioritized region in the system, while a MM force field may be used to describe remaining atoms. Iron-containing molecules are particularly important targets of DFT calculations. From the viewpoint of chemistry, this is mainly because iron is abundant on earth, iron plays powerful (and often enigmatic) roles in enzyme catalysis, and iron thus has the great potential for biomimetic catalysis of chemically difficult transformations. In this paper, we present a brief overview of several recent applications of DFT to iron-containing non-heme synthetic complexes, heme-type cytochrome P450 enzymes, and non-heme iron enzymes, all of which are of particular interest in the field of bioinorganic chemistry. Emphasis will be placed on our own work.

Standard role for a conserved aspartate or more direct involvement in deglycosylation? An ONIOM and MD investigation of adenine-DNA glycosylase

8-Oxoguanine (OG) is one of the most frequently occurring forms of DNA damage and is particularly deleterious since it forms a stable Hoogsteen base pair with adenine (A). The repair of an OG:A mispair is initiated by adenine-DNA glycosylase (MutY), which hydrolyzes the sugar-nucleobase bond of the adenine residue before the lesion is processed by other proteins. MutY has been proposed to use a two-part chemical step involving protonation of the adenine nucleobase, followed by SN1 hydrolysis of the glycosidic bond. However, differences between a recent (fluorine recognition complex, denoted as the FLRC) crystal structure and the structure on which most mechanistic conclusions have been based to date (namely, the lesion recognition complex or LRC) raise questions regarding the mechanism used by MutY and the discrete role of various active-site residues. The present work uses both molecular dynamics (MD) and quantum mechanical (ONIOM) models to compare the active-site conformational dynamics in the two crystal structures, which suggests that only the understudied FLRC leads to a catalytically competent reactant. Indeed, all previous computational studies on MutY have been initiated from the LRC structure. Subsequently, for the first time, various mechanisms are examined with detailed ONIOM(M06-2X:PM6) reaction potential energy surfaces (PES) based on the FLRC structure, which significantly extends the mechanistic picture. Specifically, our work reveals that the reaction proceeds through a different route than the commonly accepted mechanism and the catalytic function of various active-site residues (Geobacillus stearothermophilus numbering). Specifically, contrary to proposals based on the LRC, E43 is determined to solely be involved in the initial adenine protonation step and not the deglycosylation reaction as the general base. Additionally, a novel catalytic role is proposed for Y126, whereby this residue plays a significant role in stabilizing the highly charged active site, primarily through interactions with E43. More importantly, D144 is found to explicitly catalyze the nucleobase dissociation step through partial nucleophilic attack. Although this is a more direct role than previously proposed for any other DNA glycosylase, comparison to previous work on other glycosylases justifies the larger contribution in the case of MutY and allows us to propose a unified role for the conserved Asp/Glu in the DNA glycosylases, as well as other enzymes that catalyze nucleotide deglycosylation reactions.