Protein Environment Facilitates O2 Binding in Non-Heme Iron Enzyme. An Insight from ONIOM Calculations on Isopenicillin N Synthase (IPNS)

X-ray free-electron laser studies reveal correlated motion during isopenicillin N synthase catalysis

Isopenicillin N synthase (IPNS) catalyzes the unique reaction of l-δ-(α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) with dioxygen giving isopenicillin N (IPN), the precursor of all natural penicillins and cephalosporins. X-ray free-electron laser studies including time-resolved crystallography and emission spectroscopy reveal how reaction of IPNS:Fe(II):ACV with dioxygen to yield an Fe(III) superoxide causes differences in active site volume and unexpected conformational changes that propagate to structurally remote regions. Combined with solution studies, the results reveal the importance of protein dynamics in regulating intermediate conformations during conversion of ACV to IPN. The results have implications for catalysis by multiple IPNS-related oxygenases, including those involved in the human hypoxic response, and highlight the power of serial femtosecond crystallography to provide insight into long-range enzyme dynamics during reactions presently impossible for nonprotein catalysts.

Isopenicillin N Synthase: Crystallographic Studies

Isopenicillin N synthase (IPNS) is a non-heme iron oxidase (NHIO) that catalyses the cyclisation of tripeptide δ-(l-α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) to bicyclic isopenicillin N (IPN). Over the last 25 years, crystallography has shed considerable light on the mechanism of IPNS catalysis. The first crystal structure, for apo-IPNS with Mn bound in place of Fe at the active site, reported in 1995, was also the first structure for a member of the wider NHIO family. This was followed by the anaerobic enzyme-substrate complex IPNS-Fe-ACV (1997), this complex plus nitric oxide as a surrogate for co-substrate dioxygen (1997), and an enzyme product complex (1999). Since then, crystallography has been used to probe many aspects of the IPNS reaction mechanism, by crystallising the protein with a diversity of substrate analogues and triggering the oxidative reaction by using elevated oxygen pressures to force the gaseous co-substrate throughout protein crystals and maximise synchronicity of turnover in crystallo. In this way, X-ray structures have been elucidated for a range of complexes closely related to and/or directly derived from key intermediates in the catalytic cycle, thereby answering numerous mechanistic questions that had arisen from solution-phase experiments, and posing many new ones. The results of these crystallographic studies have, in turn, informed computational experiments that have brought further insight. These combined crystallographic and computational investigations augment and extend the results of earlier spectroscopic analyses and solution phase studies of IPNS turnover, to enrich our understanding of this important protein and the wider NHIO enzyme family.

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.

Mechanistic insights into a non-heme 2-oxoglutarate-dependent ethylene-forming enzyme: selectivity of ethylene-formation versusl-Arg hydroxylation

The ethylene-forming enzyme (EFE) is a unique member of the Fe(ii)- and 2-oxoglutarate-dependent (Fe/2OG) oxygenases. It converts 2OG into ethylene plus three CO2 molecules (ethylene-forming reaction) and also catalyzes the C5 hydroxylation of l-arginine coupled to the oxidative decarboxylation of 2OG (l-Arg hydroxylation reaction). To uncover the mechanisms of the dual transformations by EFE, quantum mechanical/molecular mechanical (QM/MM) calculations were carried out. Based on the results, a branched mechanism was proposed. An FeII-peroxysuccinate complex with a dissociated CO2 generated through the nucleophilic attack of the superoxo moiety of the Fe-O2 species on the keto carbon of 2OG is the key common intermediate in both reactions. A competition between the subsequent CO2 insertion (a key step in the ethylene-forming pathway) and the O-O bond cleavage (leading to the formation of succinate) governs the product selectivity. The calculated reaction barriers suggested that the CO2 insertion is favored over the O-O bond cleavage. This is consistent with the product preference observed in experiments. By comparison with the results of AsqJ (an Fe/2OG oxygenase that leads to substrate oxidation exclusively), the protein environment was found to be crucial for the selectivity. Further calculations demonstrated that the local electric field of the protein environment in EFE promotes ethylene formation by acting as a charge template, exemplifying the importance of the electrostatic interaction in enzyme catalysis. These findings offer mechanistic insights into the EFE catalysis and provide important clues for better understanding the unique ethylene-forming capability of EFE compared with other Fe/2OG oxygenases.

Roles of 2-oxoglutarate oxygenases and isopenicillin N synthase in β-lactam biosynthesis

Covering: up to 2017 2-Oxoglutarate (2OG) dependent oxygenases and the homologous oxidase isopenicillin N synthase (IPNS) play crucial roles in the biosynthesis of β-lactam ring containing natural products. IPNS catalyses formation of the bicyclic penicillin nucleus from a tripeptide. 2OG oxygenases catalyse reactions that diversify the chemistry of β-lactams formed by both IPNS and non-oxidative enzymes. Reactions catalysed by the 2OG oxygenases of β-lactam biosynthesis not only involve their typical hydroxylation reactions, but also desaturation, epimerisation, rearrangement, and ring-forming reactions. Some of the enzymes involved in β-lactam biosynthesis exhibit remarkable substrate and product selectivities. We review the roles of 2OG oxygenases and IPNS in β-lactam biosynthesis, highlighting opportunities for application of knowledge of their roles, structures, and mechanisms.

Unique chemistry of non-heme iron enzymes in fungal biosynthetic pathways

Reactions by non-heme iron enzymes in structurally intriguing fungal natural products pathways are summarized and discussed.

Multistructural microiteration technique for geometry optimization and reaction path calculation in large systems

We propose a multistructural microiteration (MSM) method for geometry optimization and reaction path calculation in large systems. MSM is a simple extension of the geometrical microiteration technique. In conventional microiteration, the structure of the non-reaction-center (surrounding) part is optimized by fixing atoms in the reaction-center part before displacements of the reaction-center atoms. In this method, the surrounding part is described as the weighted sum of multiple surrounding structures that are independently optimized. Then, geometric displacements of the reaction-center atoms are performed in the mean field generated by the weighted sum of the surrounding parts. MSM was combined with the QM/MM-ONIOM method and applied to chemical reactions in aqueous solution or enzyme. In all three cases, MSM gave lower reaction energy profiles than the QM/MM-ONIOM-microiteration method over the entire reaction paths with comparable computational costs. © 2017 Wiley Periodicals, Inc.

Artificial Force Induced Reaction Method for Systematic Determination of Complex Reaction Mechanisms

Nowadays, computational studies are very important for the elucidation of reaction mechanisms and selectivity of complex reactions. However, traditional computational methods usually require an estimated reaction path, mainly driven by limited experimental implications, intuition, and assumptions of stationary points. However, the artificial force induced reaction (AFIR) method in the global reaction route mapping (GRRM) strategy can be used for unbiased and automatic reaction path searches for complex reactions. In this account, we highlight applications of the AFIR method to a variety of reactions (organic, organometallic, enzymatic, and photochemical) of complex molecular systems. In addition, the AFIR method has been successfully used to rationalise the origin of stereo- and regioselectivity. The AFIR method can be applied from small to large molecular systems, and will be a very useful tool for the study of complex molecular problems in many areas of chemistry, biology, and material sciences.

Spectroscopic Evidence for the Two C-H-Cleaving Intermediates of Aspergillus nidulans Isopenicillin N Synthase

The enzyme isopenicillin N synthase (IPNS) installs the β-lactam and thiazolidine rings of the penicillin core into the linear tripeptide l-δ-aminoadipoyl-l-Cys-d-Val (ACV) on the pathways to a number of important antibacterial drugs. A classic set of enzymological and crystallographic studies by Baldwin and co-workers established that this overall four-electron oxidation occurs by a sequence of two oxidative cyclizations, with the β-lactam ring being installed first and the thiazolidine ring second. Each phase requires cleavage of an aliphatic C-H bond of the substrate: the pro-S-CCys,β-H bond for closure of the β-lactam ring, and the CVal,β-H bond for installation of the thiazolidine ring. IPNS uses a mononuclear non-heme-iron(II) cofactor and dioxygen as cosubstrate to cleave these C-H bonds and direct the ring closures. Despite the intense scrutiny to which the enzyme has been subjected, the identities of the oxidized iron intermediates that cleave the C-H bonds have been addressed only computationally; no experimental insight into their geometric or electronic structures has been reported. In this work, we have employed a combination of transient-state-kinetic and spectroscopic methods, together with the specifically deuterium-labeled substrates, A[d2-C]V and AC[d8-V], to identify both C-H-cleaving intermediates. The results show that they are high-spin Fe(III)-superoxo and high-spin Fe(IV)-oxo complexes, respectively, in agreement with published mechanistic proposals derived computationally from Baldwin's founding work.

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

JMJD2A catalyses the demethylation of di- and trimethylated lysine residues in histone tails and is a target for the development of new anticancer medicines. Mechanistic details of demethylation are yet to be elucidated and are important for the understanding of epigenetic processes. We have evaluated the initial step of histone demethylation by JMJD2A and demonstrate the dramatic effect of the protein environment upon oxygen binding using quantum mechanics/molecular mechanics (QM/MM) calculations. The changes in electronic structure have been studied for possible spin states and different conformations of O2 , using a combination of quantum and classical simulations. O2 binding to this histone demethylase is computed to occur preferentially as an end-on superoxo radical bound to a high-spin ferric centre, yielding an overall quintet ground state. The favourability of binding is strongly influenced by the surrounding protein: we have quantified this effect using an energy decomposition scheme into electrostatic and dispersion contributions. His182 and the methylated lysine assist while Glu184 and the oxoglutarate cofactor are deleterious for O2 binding. Charge separation in the superoxo-intermediate benefits from the electrostatic stabilization provided by the surrounding residues, stabilizing the binding process significantly. This work demonstrates the importance of the extended protein environment in oxygen binding, and the role of energy decomposition in understanding the physical origin of binding/recognition.

The Importance of the MM Environment and the Selection of the QM Method in QM/MM Calculations: Applications to Enzymatic Reactions

In this chapter, we discuss the influence of an anisotropic protein environment on the reaction mechanisms of saccharopine reductase and uroporphyrinogen decarboxylase, respectively, via the use of a quantum mechanical and molecular mechanical (QM/MM) approach. In addition, we discuss the importance of selecting a suitable DFT functional to be used in a QM/MM study of a key intermediate in the mechanism of 8R-lipoxygenase, a nonheme iron enzyme. In the case of saccharopine reductase, while the enzyme utilizes a substrate-assisted catalytic pathway, it was found that only through treating the polarizing effect of the active site, via the use of an electronic embedding formalism, was agreement with experimental kinetic data obtained. Similarly, in the case of uroporphyrinogen decarboxylase, the effect of the protein environment on the catalytic mechanism was found to be such that the calculated rate-limiting barrier is in good agreement with related experimentally determined values for the first decarboxylation of the substrate. For 8R-lipoxygenase, it was found that the geometries and energies of the multicentered open-shell intermediate complexes formed during the mechanism are quite sensitive to the choice of the density functional theory method. Thus, while density functional theory has become the method of choice in QM/MM studies, care must be taken in the selection of a particular high-level method.

A DFT and ONIOM study of C–H hydroxylation catalyzed by nitrobenzene 1,2-dioxygenase

The DFT and ONIOM calculations show that C–H hydroxylation by nitrobenzene 1,2-dioxygenase proceeds through a HO–Fe^VO intermediate.

Binding affinity of substituted ureido-benzenesulfonamide ligands to the carbonic anhydrase receptor: a theoretical study of enzyme inhibition

The binding properties of a series of benzenesulfonamide inhibitors (4-substituted-ureido-benzenesulfonamides, UBSAs) of human carbonic anhydrase II (hCA II) enzyme with active site residues have been studied using a hybrid quantum mechanical/molecular mechanical (QM/MM) model. To account for the important docking interactions between the UBSAs ligand and hCA II enzyme, a molecular docking program AutoDock Vina is used. The molecular docking results obtained by AutoDock Vina revealed that the docked conformer has root mean square deviation value less than 1.50 Å compared to X-ray crystal structures. The inhibitory activity of UBSA ligands against hCA II is found to be in good agreement with the experimental results. The thermodynamic parameters for inhibitor binding show that hydrogen bonding, hydrophilic, and hydrophobic interactions play a major role in explaining the diverse inhibitory range of these derivatives. Additionally, natural bond orbital analysis is performed to characterize the ligand-metal charge transfer stability. The insights gained from this study have great potential to design new hCA-II inhibitor, 4-[3-(1-p-Tolyl-4-trifluoromethyl-1H-pyrazol-3-yl)-ureido]-benzenesulfonamide, which belongs to the family of UBSA inhibitors and shows similar type of inhibitor potency with hCA II. This work also reveals that a QM/MM model and molecular docking method are computationally feasible and accurate for studying substrate-protein inhibition.

Elucidating the catalytic mechanism of β-secretase (BACE1): a quantum mechanics/molecular mechanics (QM/MM) approach

In this quantum mechanics/molecular mechanics (QM/MM) study, the mechanisms of the hydrolytic cleavage of the Met2-Asp3 and Leu2-Asp3 peptide bonds of the amyloid precursor protein (WT-substrate) and its Swedish mutant (SW) respectively catalyzed by β-secretase (BACE1) have been investigated by explicitly including the electrostatic and steric effects of the protein environment in the calculations. BACE1 catalyzes the rate-determining step in the generation of Alzheimer amyloid beta peptides and is widely acknowledged as a promising therapeutic target. The general acid-base mechanism followed by the enzyme proceeds through the following two steps: (1) formation of the gem-diol intermediate and (2) cleavage of the peptide bond. The formation of the gem-diol intermediate occurs with the barriers of 19.6 and 16.1 kcal/mol for the WT- and SW-substrate respectively. The QM/MM energetics predict that with the barriers of 21.9 and 17.2 kcal/mol for the WT- and SW-substrate respectively the cleavage of the peptide bond occurs in the rate-determining step. The computed barriers are in excellent agreement with the measured barrier of ∼18.0 kcal/mol for the SW-substrate and in line with the experimental observation that the cleavage of this substrate is sixty times more efficient than the WT-substrate.

Combined QM/MM (ONIOM) and QSAR approach to the study of complex formation of matrix metalloproteinase‑9 with a series of biphenylsulfonamides−LERE-QSAR analysis (V)

We previously proposed a novel QSAR (quantitative structure-activity relationship) procedure called LERE (linear expression by representative energy terms)-QSAR involving molecular calculations such as an ab initio fragment molecular orbital ones. In the present work, we applied LERE-QSAR to complex formation of matrix metalloproteinase-9 (MMP-9) with a series of substituted biphenylsulfonamides. The results shows that the overall free-energy change accompanying complex formation is due to predominantly the contribution from the electrostatic interaction with the zinc atom in the active site of MMP-9. Carbonic anhydrase (CA) belongs to the zinc-containing protease family. In contrast to the current case of MMP-9, the overall free-energy change during complex formation of CA with a series of benzenesulfonamides is due to the contributions from the solvation and dissociation free-energy changes, as previously reported. Comparison of the two sets of results indicates quantitative differences in the relative contributions of free-energy components to the overall free-energy change between the two data sets, corresponding with those in the respective classical QSAR equations. The LERE-QSAR procedure was demonstrated to quantitatively reveal differences in the binding mechanisms between the two cases involving similar but different zinc-containing proteins at the electronic and atomic levels.

[Recent progress in the theoretical studies of structure, function, and reaction of biological molecules]

Essential biomolecular functions often involve electron-related events such as chemical reactions and photoluminescence phenomena. Theoretical description of such electronic processes requires the use of quantum mechanics (QM), but the number of atoms that can be handled with QM is usually smaller than the number of atoms present in a single protein. A reasonable strategy is therefore to give priority to a few tens or hundreds of atoms in the system and deal with them quantum mechanically. Lower-priority atoms influence the event occurring in the higher-priority area; therefore, their effect should also be taken into account. Under these circumstances, a reasonable approach is to apply two or more different theoretical methods to differently prioritized subsystems. QM can be combined, for example, with less accurate yet much less demanding molecular mechanics (MM). Our own N-layered integrated molecular orbital and molecular mechanics (ONIOM) method allows for such hybrid calculations, and our group has been applying it to a wide range of biology-related problems. In this paper, we briefly explain the theoretical background and the procedure for the theoretical investigation of biological systems. Subsequently, we provide an overview of some of our recent studies of metalloenzymes and photobiology-related problems.

How molecular oxygen binds to bis[trifluoroacetylacetonato(-1)]cobalt(II)--ab initio and density functional theory studies

Cobalt(II) diketonate complexes, such as bis[trifluoroacetylacetonato(-1)]cobalt(II) [Co(tfa)(2)], catalyze the aerobic oxidation of alkenols into functionalized tetrahydrofurans. To gain insight into activation of triplet dioxygen by Co(tfa)(2) in a protic solvent, as used in oxidation catalysis, the electronic structure of aquabis[trifluoroacetylacetonato(-1)]cobalt(II)--Co(tfa)(2)(H(2)O)--and the derived dioxygen adduct were characterized using ab initio (CASSCF, NEVPT2) and density functional theory (BP86, TPSSh, B3LYP) methods. The ground state of Co(tfa)(2)(H(2)O) is a high-spin, quartet state. As dioxygen approaches the cobalt atom, the quartet state couples with a triplet dioxygen molecule and forms a sextet, a quartet, and a doublet spin state with the high-spin state being the lowest in energy. At the equilibrium Co-O(2) distance of 1.9 Å, Co(tfa)(2)(H(2)O)(O(2)) has a doublet superoxo Co(III) ground state with the unpaired electron residing on the oxygen moiety, in a nearly unchanged O(2)π* orbital.

Formation of the iron-oxo hydroxylating species in the catalytic cycle of aromatic amino acid hydroxylases

The first part of the catalytic cycle of the pterin-dependent, dioxygen-using nonheme-iron aromatic amino acid hydroxylases, leading to the Fe(IV)=O hydroxylating intermediate, has been investigated by means of density functional theory. The starting structure in the present investigation is the water-free Fe-O(2) complex cluster model that represents the catalytically competent form of the enzymes. A model for this structure was obtained in a previous study of water-ligand dissociation from the hexacoordinate model complex of the X-ray crystal structure of the catalytic domain of phenylalanine hydroxylase in complex with the cofactor (6R)-L-erythro-5,6,7,8-tetrahydrobiopterin (BH(4)) (PAH-Fe(II)-BH(4)). The O-O bond rupture and two-electron oxidation of the cofactor are found to take place via a Fe-O-O-BH(4) bridge structure that is formed in consecutive radical reactions involving a superoxide ion, O(2)(-). The overall effective free-energy barrier to formation of the Fe(IV)=O species is calculated to be 13.9 kcal  mol(-1), less than 2 kcal  mol(-1) lower than that derived from experiment. The rate-limiting step is associated with a one-electron transfer from the cofactor to dioxygen, whereas the spin inversion needed to arrive at the quintet state in which the O-O bond cleavage is finalized, essentially proceeds without activation.

Oxygen cleavage with manganese and iron in ribonucleotide reductase from Chlamydia trachomatis

The oxygen cleavage in Chlamydia trachomatis ribonucleotide reductase (RNR) has been studied using B3LYP* hybrid density functional theory. Class Ic C. trachomatis RNR lacks the radical-bearing tyrosine, crucial for activity in conventional class I (subclass a and b) RNR. Instead of the Fe(III)Fe(III)-Tyr(rad) active state, C. trachomatis RNR has a mixed Mn(IV)Fe(III) metal center in subunit II (R2). A mixed MnFe metal center has never been observed as a radical cofactor before. The active state is generated by reductive oxygen cleavage at the metal site. On the basis of calculated barriers for oxygen cleavage in C. trachomatis R2 and R2 from Escherichia coli with a diiron, a mixed manganese-iron, and a dimanganese center, conclusions can be drawn about the effect of changing metals in R2. The oxygen cleavage is found to be governed by two factors: the redox potentials of the metals and the relative stability of the different peroxides. Mn(IV) has higher stability than Fe(IV), and the barrier is therefore lower with a mixed metal center than with a diiron center. With a dimanganese center, an asymmetric peroxide is more stable than the symmetric peroxide, and the barrier therefore becomes too high. Calculated proton-coupled redox potentials are compared to identify three possible R2 active states, the Fe(III)Fe(III)-Tyr(rad) state, the Mn(IV)Fe(III) state, and the Mn(IV)Mn(IV) state. A tentative energy profile of the thermodynamics of the radical transfer from R2 to subunit I is constructed to illustrate how the stability of the active states can be understood from a thermodynamical point of view.

Correlation analyses on binding affinity of substituted benzenesulfonamides with carbonic anhydrase using ab initio MO calculations on their complex structures

Quantitative structure-activity relationship analyses on the free energy change during complex formation between substituted benzenesulfonamides (BSAs) and bovine carbonic anhydrase II (bCA II) were performed using generilized Born/surface area (GB/SA) and ab initio fragment molecular orbital (FMO) calculations for the whole complex structures. The result shows that the overall free energy change is governed by the contribution from solvation and dissociation free energy changes accompanying by complex formation. The FMO-IFIE (interfragment interaction energy) analysis quantitatively revealed that the intrinsic interaction energy of bCA II with BSAs is mostly from interactions with amino acid residues in the active site of bCA II. The "Zn block" (Zn(2+) and three histidine residues coordinated to Zn(2+)) in the active site shows the lowest interaction energy and the greatest variance of interaction energy with BSAs through their coordination interaction. The proposed procedure was demonstrated to provide a quantitative basis for understanding a ligand-protein interaction at electronic and atomic levels.

Production and isolation of ligated metal(IV)-oxo ions by tandem mass spectrometry

High valent metal(IV)-oxo species, [M(==O)(MeIm)(n)(OAc)](+) (M = Mn-Ni, MeIm = 1-methylimidazole, n = 1-2), which are relevant to biology and oxidative catalysis, were produced and isolated in gas-phase reactions of the metal(II) precursor ions [M(MeIm)(n)(OAc)](+) (M = Mn-Zn, n = 1-3) with ozone. The precursor ions [M(MeIm)(OAc)](+) and [M(MeIm)(2)(OAc)](+) were generated via collision-induced dissociation of the corresponding [M(MeIm)(3)(OAc)](+) ion. The dependence of ozone reactivity on metal and coordination number is discussed.

Reaction coordinate of isopenicillin N synthase: oxidase versus oxygenase activity

Isopenicillin N synthase (IPNS) can have both oxidase and oxygenase activity depending on the substrate. For the native substrate, ACV, oxidase activity exists; however, for the substrate analogue ACOV, which lacks an amide nitrogen, IPNS exhibits oxygenase activity. The potential energy surfaces for the O-O bond elongation and cleavage were calculated for three different reactions: homolytic cleavage via traditional Fenton chemistry, heterolytic cleavage, and nucleophilic attack. These surfaces show that the hydroperoxide-ferrous intermediate, formed by O(2)-activated H atom abstraction from the substrate, can exploit different reaction pathways and that interactions with the substrate govern the pathway. The hydrogen bonds from hydroperoxide to the amide nitrogen of ACV polarize the sigma* orbital of the peroxide toward the proximal oxygen, facilitating heterolytic cleavage. For the substrate analogue ACOV, this hydrogen bond is no longer present, leading to nucleophilic attack on the substrate intermediate C-S bond. After cleavage of the hydroperoxide, the two reaction pathways proceed with minimal barriers, resulting in the closure of the beta-lactam ring for the oxidase activity (ACV) or formation of the thiocarboxylate for oxygenase activity (ACOV).

Insights into the (superoxo)Fe(III)Fe(III) intermediate and reaction mechanism of myo-inositol oxygenase: DFT and ONIOM(DFT:MM) study

The (superoxo)Fe(III)Fe(III) reactive species and the catalytic reaction mechanism of a diiron enzyme, myo-inositol oxygenase (MIOX), were theoretically investigated by means of density functional theory (DFT) and ONIOM quantum mechanical/molecular mechanical (QM/MM) approaches. The ground state of the (superoxo)Fe(III)Fe(III) intermediate was shown to have a side-on coordination geometry and an S = 1/2 spin state, wherein the two iron sites are antiferromagnetically coupled while the superoxide site and the nearest iron are ferromagnetically coupled. A full reaction pathway leading to a D-glucuronate product from myo-inositol was proposed based on ONIOM computational results. Two major roles of the enzyme surrounding during the catalytic reaction were identified. One is to facilitate the initial H-abstraction step, and the other is to restrict the movement of the substrate via H-bonding interactions in order to avoid unwanted side reactions. In our proposed mechanism, O-O bond cleavage has the highest barrier, thus constituting the rate-limiting step of the reaction. The unique role of the bridging hydroxide ligand as a catalytic base was also identified.

C-X...H contacts in biomolecular systems: how they contribute to protein-ligand binding affinity

The hydrogen bond acceptor capability of halogens has long been underappreciated in the field of biology. In this work, we have surveyed structures of protein complexes with halogenated ligands to characterize geometrical preferences of C-X...H contacts and contributions of such interactions to protein-ligand binding affinity. Notably, F...H interactions in biomolecules exhibit a remarkably different behavior as compared to three other kinds of X...H (X = Cl, Br, I) interactions, which has been rationalized by means of ab initio calculations using simple model systems. The C-X...H contacts in biological systems are characterized as weak hydrogen bonding interactions. Furthermore, the electrophile "head on" and nucleophile "side on" interactions of halogens have been extensively investigated through the examination of interactions in protein structures and a two-layer ONIOM-based QM/MM method. In biomolecular systems, C-X...H contacts are recognized as secondary interaction contributions to C-X...O halogen bonds that play important roles in conferring specificity and affinity for halogenated ligands. The results presented here are within the context of their potential applications in drug design, including relevance to the development of accurate force fields for halogens.

Origin of the correlation of the rate constant of substrate hydroxylation by nonheme iron(IV)-oxo complexes with the bond-dissociation energy of the C-H bond of the substrate

Mononuclear nonheme iron containing systems are versatile and vital oxidants of substrate hydroxylation reactions in many biosystems, whereby the rate constant of hydroxylation correlates with the strength of the C-H bond that is broken in the process. The thermodynamic reason behind these correlations, however, has never been established. In this work results of a series of density functional theory calculations of substrate hydroxylation by a mononuclear nonheme iron(IV)-oxo oxidant with a 2 His/1 Asp structural motif analogous to alpha-ketoglutarate dependent dioxygenases are presented. The calculations show that these oxidants are very efficient and able to hydroxylate strong C-H bonds, whereby the hydrogen abstraction barriers correlate linearly with the strength of the C-H bond of the substrate that is broken. These trends have been rationalized using a valence bond (VB) curve-crossing diagram, which explains the correlation using electron transfer mechanisms in the hydrogen abstraction processes. We also rationalized the subsequent reaction step for radical rebound and show that the barrier is proportional to the electron affinity of the iron(III)-hydroxo intermediate complex. It is shown that nonheme iron(IV)-hydroxo complexes have a larger electron affinity than heme iron(IV)-hydroxo complexes and therefore also experience larger radical rebound barriers, which may have implications for product distributions and rearrangement reactions. Thus, detailed comparisons between heme and nonheme iron(IV)-oxo oxidants reveal the fundamental differences in monoxygenation capabilities of these important classes of oxidants in biosystems and synthetic analogues for the first time and enable us to make predictions of experimental processes.

A DFT study of nucleobase dealkylation by the DNA repair enzyme AlkB

Oxidative dealkylation is a unique mechanistic pathway found in the alpha-ketoglutarate-Fe(II)-dependent AlkB family of enzymes to remove the alkylation damage to DNA bases and regenerate nucleobases to their native state. The B3LYP density functional combined with a self-consistent reaction field was used to explore the triplet, quintet, and septet spin-state potential energy surfaces of the multistep catalytic mechanism of AlkB. The mechanism was found to consist of four stages. First, binding of dioxygen to iron in the active-site complex occurs concerted with electron transfer, thereby yielding a ferric-superoxido species. Second, competing initiation for the activation of oxygen to generate the high-valent iron-oxygen intermediates (ferryl-oxo Fe(IV)O and ferric-oxyl Fe(III)O(*) species) was found to occur on the quintet and septet surfaces. Then, conformational reorientation of the activated iron-oxygen ligand was found to be nearly thermoneutral with a barrier of ca. 50 kJ mol(-1). The final stage is the oxidative dealkylation of the damaged nucleobase with the rate-controlling step being the abstraction of a hydrogen atom from the damaging methyl group by the ferryl-oxo ligand. For this step, the calculated barrier of 87.4 kJ mol(-1) is in good agreement with the experimental activation energy of ca. 83 kJ mol(-1) for the enzyme-catalyzed reaction.