X-ray absorption spectroscopy of chloroperoxidase compound I: Insight into the reactive intermediate of P450 chemistry

Challenges Relating to the Quantification of Ferryl(IV) Ion and Hydroxyl Radical Generation Rates Using Methyl Phenyl Sulfoxide (PMSO), Phthalhydrazide, and Benzoic Acid as Probe Compounds in the Homogeneous Fenton Reaction

Ferryl ion ([FeIVO]2+) has often been suggested to play a role in iron-based advanced oxidation processes (AOPs) with its presence commonly determined using the unique oxidation pathway from methyl phenyl sulfoxide (PMSO) to methyl phenyl sulfone (PMSO2). However, we show here that the oxidation products of PMSO, formed on reaction with hydroxyl radical, enhance PMSO2 formation as a result of their complexation with Fe(III) leading to the changes in the reactivity of Fe(III) species in the homogeneous Fenton reaction. As such, PMSO should be used with caution to investigate the role of [FeIVO]2+ in iron-based AOPs with these insights suggesting the need to reassess the findings of many previous studies in which this reagent was used. The other common target compounds, phthalhydrazide and hydroxybenzoic acids, were also found to modify the rate and extent of iron cycling as a result of complexation and/or redox reactions, either by the probe compound itself and/or oxidation products formed. Overall, this study highlights that these confounding effects of the aromatic probe compounds on the reactivity of iron species should be recognized if reliable mechanistic insights into iron-based AOPs are to be obtained.

How Is Substrate Halogenation Triggered by the Vanadium Haloperoxidase from Curvularia inaequalis?

Vanadium haloperoxidases (VHPOs) are unique enzymes in biology that catalyze a challenging halogen transfer reaction and convert a strong aromatic C-H bond into C-X (X = Cl, Br, I) with the use of a vanadium cofactor and H₂O₂. The VHPO catalytic cycle starts with the conversion of hydrogen peroxide and halide (X = Cl, Br, I) into hypohalide on the vanadate cofactor, and the hypohalide subsequently reacts with a substrate. However, it is unclear whether the hypohalide is released from the enzyme or otherwise trapped within the enzyme structure for the halogenation of organic substrates. A substrate-binding pocket has never been identified for the VHPO enzyme, which questions the role of the protein in the overall reaction mechanism. Probing its role in the halogenation of small molecules will enable further engineering of the enzyme and expand its substrate scope and selectivity further for use in biotechnological applications as an environmentally benign alternative to current organic chemistry synthesis. Using a combined experimental and computational approach, we elucidate the role of the vanadium haloperoxidase protein in substrate halogenation. Activity studies show that binding of the substrate to the enzyme is essential for the reaction of the hypohalide with substrate. Stopped-flow measurements demonstrate that the rate-determining step is not dependent on substrate binding but partially on hypohalide formation. Using a combination of molecular mechanics (MM) and molecular dynamics (MD) simulations, the substrate binding area in the protein is identified and even though the selected substrates (methylphenylindole and 2-phenylindole) have limited hydrogen-bonding abilities, they are found to bind relatively strongly and remain stable in a binding tunnel. A subsequent analysis of the MD snapshots characterizes two small tunnels leading from the vanadate active site to the surface that could fit small molecules such as hypohalide, halide, and hydrogen peroxide. Density functional theory studies using electric field effects show that a polarized environment in a specific direction can substantially lower barriers for halogen transfer. A further analysis of the protein structure indeed shows a large dipole orientation in the substrate-binding pocket that could enable halogen transfer through an applied local electric field. These findings highlight the importance of the enzyme in catalyzing substrate halogenation by providing an optimal environment to lower the energy barrier for this challenging aromatic halide insertion reaction.

Heme-Dependent Supramolecular Nanocatalysts: A Review

Enzymes fold into three-dimensional structures to distribute amino acid residues for catalysis, which inspired the supramolecular approach to construct enzyme-mimicking catalysts. A key concern in the development of supramolecular strategies is the ability to confine and orient functional groups to form enzyme-like active sites in artificial materials. This review introduces the design principles and construction of supramolecular nanomaterials exhibiting catalytic functions of heme-dependent enzymes, a large class of metalloproteins, which rely on a heme cofactor and spatially configured residues to catalyze diverse reactions via a complex multistep mechanism. We focus on the structure-activity relationship of the supramolecular catalysts and their applications in materials synthesis/degradation, biosensing, and therapeutics. The heme-free catalysts that catalyze reactions achieved by hemeproteins are also briefly discussed. Towards the end of the review, we discuss the outlook on the challenges related to catalyst design and future prospective, including the development of structure-resolving techniques and design concepts, with the aim of creating enzyme-mimicking materials that possess catalytic power rivaling that of natural enzymes..

Role of "S" Substitution on C-H Activation Reactivity of Iron(IV)-Oxo Cyclam Complexes: a Computational Investigation

A comprehensive density functional theory (DFT) investigation has been presented in this article to address the role of equatorial sulfur ligation in C-H activation. A non-heme iron-oxo compound with four nitrogen atoms constituting the equatorially connected macrocyclic framework (represented as N₄) [Fe(IV)═O(THC)(CH₃CN)]²⁺(THC = 1,4,8,11-tetrahydro1,4,8,11-tetraazacyclotetradecane) has been considered as the base compound. Other complexes have been anticipated by the sequential replacement of this nitrogen by sulfur, that is, N₄, N₃S₁, N₂S₂, N₁S₃, and S₄. Counterions, as always, have been considered to avoid the self-interaction error in DFT. Generally, the anti-conformers (with respect to equatorial N-H and Fe═O) turned out to be the most stable. It was found that with the enrichment of the equatorial sulfur atom, reactivity increases successively, that is, we get the trend N₄ < N₃S₁ < N₂S₂ < N₁S₃ < S_4. Our investigations have also verified the available experimental results where it has been reported that N₂S₂ is more reactive than N₄ in their mixed conformation. In search of insights into this typical pattern of reactivity, the interplay of several factors has been recognized, such as the distortion energy which decreases for the transition states with the addition of sulfur; the spin density on the oxygen atom which increases implying that the radical character of abstractor increases on sulfur ligation; the energy of the electron acceptor orbital (the lowest unoccupied molecular orbital (σ_z²*)) which decreases continuously with the sulfur substitution; and the triplet-quintet oxidant energy gap which decreases consistently with S enrichment in the equatorial position. The computational predictions reported here, if further validated by experiments, will definitely encourage the synthesis of sulfur-ligated bio-inspired complexes instead of the ones constituting nitrogen exclusively.

Computational Study of the Peroxygenase Mechanism Catalyzed by Hemoglobin Dehaloperoxidase Involved in the Degradation of Chlorophenols

The biochemical evidence showed that hemoglobin dehaloperoxidase (DHP B) from Amphitrite Ornata is a multifunctional hemoprotein that catalyzes both dehalogenation and hydroxylation of halophenols via the peroxidase and peroxygenase mechanism, respectively, which sets the basis for the degradation of halophenols. In the peroxygenase mechanism, the reaction was previously suggested to be triggered either by the hydrogen atom abstraction by the Fe═O center or by the proton abstraction by His55. To illuminate the peroxygenase mechanism of DHP B at the atomistic level, on the basis of the high-resolution crystal structure, computational models were constructed, and a series of quantum mechanical/molecular mechanical calculations have been performed. According to the calculation results, the pathway (Path a) initiated by the H-abstraction by the Fe═O center is feasible. In another pathway (Path b), His55 can abstract the proton from the hydroxyl group of the substrate (4-Cl-o-cresol) to initiate the reaction; however, its feasibility depends on the prior electron transfer from the substrate to the porphyrin group. The rate-limiting step of Path a is the OH-rebound, which corresponds to an energy barrier of 14.7 kcal/mol at the quartet state. His55 acts as an acid-base catalyst and directly involves in the catalysis. Our mutant study indicates that His55 can be replaced by other titratable residues. These findings may provide useful information for further understanding of the catalytic reaction of DHP B and for the design of enzymes in the degradation of pollutants, in particular, halophenols.

Tyrosyl Radical-Mediated Sequential Oxidative Decarboxylation of Coproporphyrinogen III through PCET: Theoretical Insights into the Mechanism of Coproheme Decarboxylase ChdC

The peroxide-dependent coproheme decarboxylase ChdC from Geobacillus stearothermophilus catalyzes two key steps in the synthesis of heme b, i.e., two sequential oxidative decarboxylations of coproporphyrinogen III (coproheme III) at propionate groups P2 and P4. In the binding site of coproheme III, P2 and P4 are anchored by different residues (Tyr144, Arg217, and Ser222 for P2 and Tyr113, Lys148, and Trp156 for P4); however, strong experimental evidence supports that the generated Tyr144 radical acts as an unique intermediary for hydrogen atom transfer (HAT) from both reactive propionates. So far, the reaction details are still unclear. Herein, we carried out quantum mechanics/molecular mechanics calculations to explore the decarboxylation mechanism of coproheme III. In our calculations, the coproheme Cpd I, Fe(IV) = O coupled to a porphyrin radical cation (por^•+) with four propionate groups, was used as a reactant model. Our calculations reveal that Tyr144 is directly involved in the decarboxylation of propionate group P2. First, the proton-coupled electron transfer (PCET) occurs from Tyr144 to P2, generating a Tyr144 radical, which then abstracts a hydrogen atom from the Cβ of P2. The β-H extraction was calculated to be the rate-limiting step of decarboxylation. It is the porphyrin radical cation (por^•+) that makes the PCET from Tyr144 to P2 to be quite easy to initiate the decarboxylation. Finally, the electron transfers from the Cβ^• through the porphyrin to the iron center, leading to the decarboxylation of P2. Importantly, the decarboxylation of P4 mediated by Lys148 was calculated to be very difficult, which suggests that after the P2 decarboxylation, the generated harderoheme III intermediate should rebind or rotate in the active site so that the propionate P4 occupies the binding site of P2, and Tyr144 again mediates the decarboxylation of P4. Thus, our calculations support the fact that Tyr144 is responsible for the decarboxylation of both P2 and P4.

Aspartate or arginine? Validated redox state X-ray structures elucidate mechanistic subtleties of FeIV = O formation in bacterial dye-decolorizing peroxidases

Structure determination of proteins and enzymes by X-ray crystallography remains the most widely used approach to complement functional and mechanistic studies. Capturing the structures of intact redox states in metalloenzymes is critical for assigning the chemistry carried out by the metal in the catalytic cycle. Unfortunately, X-rays interact with protein crystals to generate solvated photoelectrons that can reduce redox active metals and hence change the coordination geometry and the coupled protein structure. Approaches to mitigate such site-specific radiation damage continue to be developed, but nevertheless application of such approaches to metalloenzymes in combination with mechanistic studies are often overlooked. In this review, we summarize our recent structural and kinetic studies on a set of three heme peroxidases found in the bacterium Streptomyces lividans that each belong to the dye decolourizing peroxidase (DyP) superfamily. Kinetically, each of these DyPs has a distinct reactivity with hydrogen peroxide. Through a combination of low dose synchrotron X-ray crystallography and zero dose serial femtosecond X-ray crystallography using an X-ray free electron laser (XFEL), high-resolution structures with unambiguous redox state assignment of the ferric and ferryl (Fe^IV = O) heme species have been obtained. Experiments using stopped-flow kinetics, solvent-isotope exchange and site-directed mutagenesis with this set of redox state validated DyP structures have provided the first comprehensive kinetic and structural framework for how DyPs can modulate their distal heme pocket Asp/Arg dyad to use either the Asp or the Arg to facilitate proton transfer and rate enhancement of peroxide heterolysis.

QM-Cluster Model Study of the Guaiacol Hydrogen Atom Transfer and Oxygen Rebound with Cytochrome P450 Enzyme GcoA

The key step of the O-demethylation of guaiacol by GcoA of the cytochrome P450-reductase pair was studied with DFT using two 10-residue and three 15-residue QM-cluster models. For each model, two reaction pathways were examined, beginning with a different guaiacol orientation. Based on this study, His354, Phe349, Glu249, and Pro250 residues were found to be important for keeping the heme in a planar geometry throughout the reaction. Val241 and Gly245 residues were needed in the QM-cluster models to provide the hydrophobic pocket for an appropriate guaiacol pose in the reaction. The aromatic triad Phe75, Phe169, and Phe395 may be necessary to facilitate guaiacol migrating into the enzyme active site, but it does not qualitatively affect kinetics and thermodynamics of the proposed mechanism. All QM-cluster models created by RINRUS agree very well with previous experimental work. This study provides details for better understanding enzymatic O-demethylation of lignins to form catechol derivatives by GcoA.

A bioinspired oxoiron(IV) motif supported on a N2S2 macrocyclic ligand

A mononuclear oxoiron(iv) complex 1-trans bearing two equatorial sulfur ligations is synthesized and characterized as an active-site model of the elusive sulfur-ligated FeIV[double bond, length as m-dash]O intermediates in non-heme iron oxygenases. The introduction of sulfur ligands weakens the Fe[double bond, length as m-dash]O bond and enhances the oxidative reactivity of the FeIV[double bond, length as m-dash]O unit with a diminished deuterium kinetic isotope effect, thereby providing a compelling rationale for nature's use of the cis-thiolate ligated oxoiron(iv) motif in key metabolic transformations.

Crystals in Minutes: Instant On-Site Microcrystallisation of Various Flavours of the CYP102A1 (P450BM3) Haem Domain

Despite CYP102A1 (P450BM3) representing one of the most extensively researched metalloenzymes, crystallisation of its haem domain upon modification can be a challenge. Crystal structures are indispensable for the efficient structure-based design of P450BM3 as a biocatalyst. The abietane diterpenoid derivative N-abietoyl-l-tryptophan (AbiATrp) is an outstanding crystallisation accelerator for the wild-type P450BM3 haem domain, with visible crystals forming within 2 hours and diffracting to a near-atomic resolution of 1.22 Å. Using these crystals as seeds in a cross-microseeding approach, an assortment of P450BM3 haem domain crystal structures, containing previously uncrystallisable decoy molecules and diverse artificial metalloporphyrins binding various ligand molecules, as well as heavily tagged haem-domain variants, could be determined. Some of the structures reported herein could be used as models of different stages of the P450BM3 catalytic cycle.

X-ray snapshots reveal conformational influence on active site ligation during metalloprotein folding

Cytochrome c (cyt c) has long been utilized as a model system to study metalloprotein folding dynamics and the interplay between active site ligation and tertiary structure. However, recent reports regarding the weakness of the native Fe(ii)-S bond (Fe-Met80) call into question the role of the active site ligation in the protein folding process. In order to investigate the interplay between protein conformation and active site structures, we directly tracked the evolution of both during a photolysis-induced folding reaction using X-ray transient absorption spectroscopy and time-resolved X-ray solution scattering techniques. We observe an intermediate Fe-Met80 species appearing on ∼2 μs timescale, which should not be sustained without stabilization from the folded protein structure. We also observe the appearance of a new active site intermediate: a weakly interacting Fe-H2O state. As both intermediates require stabilization of weak metal-ligand interactions, we surmise the existence of a local structure within the unfolded protein that protects and limits the movement of the ligands, similar to the entatic state found in the native cyt c fold. Furthermore, we observe that in some of the unfolded ensemble, the local stabilizing structure is lost, leading to expansion of the unfolded protein structure and misligation to His26/His33 residues.

A Comparative Review on the Catalytic Mechanism of Nonheme Iron Hydroxylases and Halogenases

Enzymatic halogenation and haloperoxidation are unusual processes in biology; however, a range of halogenases and haloperoxidases exist that are able to transfer an aliphatic or aromatic C–H bond into C–Cl/C–Br. Haloperoxidases utilize hydrogen peroxide, and in a reaction with halides (Cl−/Br−), they react to form hypohalides (OCl−/OBr−) that subsequently react with substrate by halide transfer. There are three types of haloperoxidases, namely the iron-heme, nonheme vanadium, and flavin-dependent haloperoxidases that are reviewed here. In addition, there are the nonheme iron halogenases that show structural and functional similarity to the nonheme iron hydroxylases and form an iron(IV)-oxo active species from a reaction of molecular oxygen with α-ketoglutarate on an iron(II) center. They subsequently transfer a halide (Cl−/Br−) to an aliphatic C–H bond. We review the mechanism and function of nonheme iron halogenases and hydroxylases and show recent computational modelling studies of our group on the hectochlorin biosynthesis enzyme and prolyl-4-hydroxylase as examples of nonheme iron halogenases and hydroxylases. These studies have established the catalytic mechanism of these enzymes and show the importance of substrate and oxidant positioning on the stereo-, chemo- and regioselectivity of the reaction that takes place.

α-Oxidative decarboxylation of fatty acids catalysed by cytochrome P450 peroxygenases yielding shorter-alkyl-chain fatty acids

Shorter-alkyl-chain fatty acids such as tridecanoic acid or lauric acid were produced from myristic acid by CYP152 peroxygenases.

Control of stereoselectivity of benzylic hydroxylation catalysed by wild-type cytochrome P450BM3 using decoy molecules

The benzylic hydroxylation of non-native substrates was catalysed by cytochrome P450BM3, wherein “decoy molecules” controlled the stereoselectivity of the reactions.

Spectroscopic Investigations of Catalase Compound II: Characterization of an Iron(IV) Hydroxide Intermediate in a Non-thiolate-Ligated Heme Enzyme

We report on the protonation state of Helicobacter pylori catalase compound II. UV/visible, Mössbauer, and X-ray absorption spectroscopies have been used to examine the intermediate from pH 5 to 14. We have determined that HPC-II exists in an iron(IV) hydroxide state up to pH 11. Above this pH, the iron(IV) hydroxide complex transitions to a new species (pKa = 13.1) with Mössbauer parameters that are indicative of an iron(IV)-oxo intermediate. Recently, we discussed a role for an elevated compound II pKa in diminishing the compound I reduction potential. This has the effect of shifting the thermodynamic landscape toward the two-electron chemistry that is critical for catalase function. In catalase, a diminished potential would increase the selectivity for peroxide disproportionation over off-pathway one-electron chemistry, reducing the buildup of the inactive compound II state and reducing the need for energetically expensive electron donor molecules.

Beyond ferryl-mediated hydroxylation: 40 years of the rebound mechanism and C-H activation

Since our initial report in 1976, the oxygen rebound mechanism has become the consensus mechanistic feature for an expanding variety of enzymatic C-H functionalization reactions and small molecule biomimetic catalysts. For both the biotransformations and models, an initial hydrogen atom abstraction from the substrate (R-H) by high-valent iron-oxo species (Feⁿ=O) generates a substrate radical and a reduced iron hydroxide, [Fe^n-1-OH ·R]. This caged radical pair then evolves on a complicated energy landscape through a number of reaction pathways, such as oxygen rebound to form R-OH, rebound to a non-oxygen atom affording R-X, electron transfer of the incipient radical to yield a carbocation, R⁺, desaturation to form olefins, and radical cage escape. These various flavors of the rebound process, often in competition with each other, give rise to the wide range of C-H functionalization reactions performed by iron-containing oxygenases. In this review, we first recount the history of radical rebound mechanisms, their general features, and key intermediates involved. We will discuss in detail the factors that affect the behavior of the initial caged radical pair and the lifetimes of the incipient substrate radicals. Several representative examples of enzymatic C-H transformations are selected to illustrate how the behaviors of the radical pair [Fe^n-1-OH ·R] determine the eventual reaction outcome. Finally, we discuss the powerful potential of "radical rebound" processes as a general paradigm for developing novel C-H functionalization reactions with synthetic, biomimetic catalysts. We envision that new chemistry will continue to arise by bridging enzymatic "radical rebound" with synthetic organic chemistry.

Simulations of iron K pre-edge X-ray absorption spectra using the restricted active space method

The intensities and relative energies of metal K pre-edge features are sensitive to both geometric and electronic structures.

How Accurately Can Extended X-ray Absorption Spectra Be Predicted from First Principles? Implications for Modeling the Oxygen-Evolving Complex in Photosystem II

First principle calculations of extended X-ray absorption fine structure (EXAFS) data have seen widespread use in bioinorganic chemistry, perhaps most notably for modeling the Mn4Ca site in the oxygen evolving complex (OEC) of photosystem II (PSII). The logic implied by the calculations rests on the assumption that it is possible to a priori predict an accurate EXAFS spectrum provided that the underlying geometric structure is correct. The present study investigates the extent to which this is possible using state of the art EXAFS theory. The FEFF program is used to evaluate the ability of a multiple scattering-based approach to directly calculate the EXAFS spectrum of crystallographically defined model complexes. The results of these parameter free predictions are compared with the more traditional approach of fitting FEFF calculated spectra to experimental data. A series of seven crystallographically characterized Mn monomers and dimers is used as a test set. The largest deviations between the FEFF calculated EXAFS spectra and the experimental EXAFS spectra arise from the amplitudes. The amplitude errors result from a combination of errors in calculated S0(2) and Debye-Waller values as well as uncertainties in background subtraction. Additional errors may be attributed to structural parameters, particularly in cases where reliable high-resolution crystal structures are not available. Based on these investigations, the strengths and weaknesses of using first-principle EXAFS calculations as a predictive tool are discussed. We demonstrate that a range of DFT optimized structures of the OEC may all be considered consistent with experimental EXAFS data and that caution must be exercised when using EXAFS data to obtain topological arrangements of complex clusters.

Significantly shorter Fe-S bond in cytochrome P450-I is consistent with greater reactivity relative to chloroperoxidase

Cytochrome P450 (P450) and chloroperoxidase (CPO) are thiolate ligated heme proteins that catalyze the activation of carbon hydrogen bonds. The principal intermediate in these reactions is a ferryl radical species called compound I. P450 compound I (P450-I) is significantly more reactive than CPO-I, which only cleaves activated C-H bonds. To provide insight into the differing reactivities of these intermediates, we examined CPO-I and P450-I with variable temperature Mössbauer and X-ray absorption spectroscopies. These measurements indicate that the Fe-S bond is significantly shorter in P450-I than in CPO-I. This difference in Fe-S bond lengths can be understood in terms of variations in hydrogen bonding patterns within the “cys-pocket” (a portion of the proximal helix that encircles the thiolate ligand). Weaker hydrogen bonding in P450-I results in a shorter Fe-S bond, which enables greater electron donation from the axial-thiolate ligand. This observation may in part explain P450's greater propensity for C-H bond activation.

Enzymatic Halogenases and Haloperoxidases: Computational Studies on Mechanism and Function

Despite the fact that halogenated compounds are rare in biology, a number of organisms have developed processes to utilize halogens and in recent years, a string of enzymes have been identified that selectively insert halogen atoms into, for instance, a CH aliphatic bond. Thus, a number of natural products, including antibiotics, contain halogenated functional groups. This unusual process has great relevance to the chemical industry for stereoselective and regiospecific synthesis of haloalkanes. Currently, however, industry utilizes few applications of biological haloperoxidases and halogenases, but efforts are being worked on to understand their catalytic mechanism, so that their catalytic function can be upscaled. In this review, we summarize experimental and computational studies on the catalytic mechanism of a range of haloperoxidases and halogenases with structurally very different catalytic features and cofactors. This chapter gives an overview of heme-dependent haloperoxidases, nonheme vanadium-dependent haloperoxidases, and flavin adenine dinucleotide-dependent haloperoxidases. In addition, we discuss the S-adenosyl-l-methionine fluoridase and nonheme iron/α-ketoglutarate-dependent halogenases. In particular, computational efforts have been applied extensively for several of these haloperoxidases and halogenases and have given insight into the essential structural features that enable these enzymes to perform the unusual halogen atom transfer to substrates.

Characterizing the Intermediates Compound I and II in the Cytochrome P450 Catalytic Cycle with Nonlinear X-ray Spectroscopy: A Simulation Study

Cytochrome P450 enzymes are an important family of biocatalysts that oxidize chemically inert CH bonds. There are many unresolved questions regarding the catalytic reaction intermediates, in particular P450 Compound I (Cpd-I) and II (Cpd-II). By using simple molecular models, we simulate various X-ray spectroscopy signals, including X-ray absorption near-edge structure (XANES), resonant inelastic X-ray scattering (RIXS), and stimulated X-ray Raman spectroscopy (SXRS) of the low- and high-spin states of Cpd-I and II. Characteristic peak patterns are presented and connected to the corresponding electronic structures. These X-ray spectroscopy techniques are complementary to more conventional infrared and optical spectroscopy and they help to elucidate the evolving electronic structures of transient species along the reaction path.

Monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 enzymes

This review examines the monooxygenase, peroxidase and peroxygenase properties and reaction mechanisms of cytochrome P450 (CYP) enzymes in bacterial, archaeal and mammalian systems. CYP enzymes catalyze monooxygenation reactions by inserting one oxygen atom from O2 into an enormous number and variety of substrates. The catalytic versatility of CYP stems from its ability to functionalize unactivated carbon-hydrogen (C-H) bonds of substrates through monooxygenation. The oxidative prowess of CYP in catalyzing monooxygenation reactions is attributed primarily to a porphyrin π radical ferryl intermediate known as Compound I (CpdI) (Por•+FeIV=O), or its ferryl radical resonance form (FeIV-O•). CYP-mediated hydroxylations occur via a consensus H atom abstraction/oxygen rebound mechanism involving an initial abstraction by CpdI of a H atom from the substrate, generating a highly-reactive protonated Compound II (CpdII) intermediate (FeIV-OH) and a carbon-centered alkyl radical that rebounds onto the ferryl hydroxyl moiety to yield the hydroxylated substrate. CYP enzymes utilize hydroperoxides, peracids, perborate, percarbonate, periodate, chlorite, iodosobenzene and N-oxides as surrogate oxygen atom donors to oxygenate substrates via the shunt pathway in the absence of NAD(P)H/O2 and reduction-oxidation (redox) auxiliary proteins. It has been difficult to isolate the historically elusive CpdI intermediate in the native NAD(P)H/O2-supported monooxygenase pathway and to determine its precise electronic structure and kinetic and physicochemical properties because of its high reactivity, unstable nature (t½~2 ms) and short life cycle, prompting suggestions for participation in monooxygenation reactions of alternative CYP iron-oxygen intermediates such as the ferric-peroxo anion species (FeIII-OO-), ferric-hydroperoxo species (FeIII-OOH) and FeIII-(H2O2) complex.

How water molecules affect the catalytic activity of hydrolases--a XANES study of the local structures of peptide deformylase

Peptide deformylase (PDF) is a prokaryotic enzyme that catalyzes the deformylation of nascent peptides generated during protein synthesis and water molecules play a key role in these hydrolases. Using X-ray absorption near edge spectroscopy (XANES) and ab initio calculations we accurately probe the local atomic environment of the metal ion binding in the active site of PDF at different pH values and with different metal ions. This new approach is an effective way to monitor existing correlations among functions and structural changes. We show for the first time that the enzymatic activity depends on pH values and metal ions via the bond length of the nearest coordinating water (Wat1) to the metal ion. Combining experimental and theoretical data we may claim that PDF exhibits an enhanced enzymatic activity only when the distance of the Wat1 molecule with the metal ion falls in the limited range from 2.15 to 2.55 Å.

Setting an upper limit on the myoglobin iron(IV)hydroxide pK(a): insight into axial ligand tuning in heme protein catalysis

To provide insight into the iron(IV)hydroxide pK(a) of histidine ligated heme proteins, we have probed the active site of myoglobin compound II over the pH range of 3.9-9.5, using EXAFS, Mössbauer, and resonance Raman spectroscopies. We find no indication of ferryl protonation over this pH range, allowing us to set an upper limit of 2.7 on the iron(IV)hydroxide pK(a) in myoglobin. Together with the recent determination of an iron(IV)hydroxide pK(a) ∼ 12 in the thiolate-ligated heme enzyme cytochrome P450, this result provides insight into Nature's ability to tune catalytic function through its choice of axial ligand.

Cytochrome P450 compound I in the plane wave pseudopotential framework: GGA electronic and geometric structure of thiolate-ligated iron(IV)-oxo porphyrin

The cytochromes P450 constitute a ubiquitous family of metalloenzymes, catalyzing manifold reactions of biological and synthetic importance via a thiolate-ligated iron-oxo (IV) porphyrin radical species denoted compound I (Cpd I). Experimental investigations have implicated this intermediate in a broad spectrum of biophysically interesting phenomena, further augmenting the importance of a Cpd I model system. Ab initio molecular dynamics, including Car-Parrinello and path integral methods, conjoin electronic structure theory with finite temperature simulation, affording tools most valuable to approach such enzymes. These methods are typically driven by density functional theory (DFT) in a plane-wave pseudopotential framework; however, existing studies of Cpd I have been restricted to localized Gaussian basis sets. The appropriate choice of density functional and pseudopotential for such simulations is accordingly not obvious. To remedy this situation, a systematic benchmarking of thiolate-ligated Cpd I is performed using several generalized-gradient approximation (GGA) functionals in the Martins-Troullier and Vanderbilt ultrasoft pseudopotential schemes. The resultant electronic and structural parameters are compared to localized-basis DFT calculations using GGA and hybrid density functionals. The merits and demerits of each scheme are presented in the context of reproducing existing experimental and theoretical results for Cpd I.

Paramagnetic nuclear magnetic resonance relaxation and molecular mechanics studies of the chloroperoxidase-indole complex: insights into the mechanism of chloroperoxidase-catalyzed regioselective oxidation of indole

To unravel the mechanism of chloroperoxidase (CPO)-catalyzed regioselective oxidation of indole, we studied the structure of the CPO-indole complex using nuclear magnetic resonance (NMR) relaxation measurements and computational techniques. The dissociation constant (KD) of the CPO-indole complex was calculated to be approximately 21 mM. The distances (r) between protons of indole and the heme iron calculated via NMR relaxation measurements and molecular docking revealed that the pyrrole ring of indole is oriented toward the heme with its 2-H pointing directly at the heme iron. Both KD and r values are independent of pH in the range of 3.0-6.5. The stability and structure of the CPO-indole complex are also independent of the concentration of chloride or iodide ion. Molecular docking suggests the formation of a hydrogen bond between the NH group of indole and the carboxyl O of Glu 183 in the binding of indole to CPO. Simulated annealing of the CPO-indole complex using r values from NMR experiments as distance restraints reveals that the van der Waals interactions were much stronger than the Coulomb interactions in the binding of indole to CPO, indicating that the association of indole with CPO is primarily governed by hydrophobic rather than electrostatic interactions. This work provides the first experimental and theoretical evidence of the long-sought mechanism that leads to the "unexpected" regioselectivity of the CPO-catalyzed oxidation of indole. The structure of the CPO-indole complex will serve as a lighthouse in guiding the design of CPO mutants with tailor-made activities for biotechnological applications.

Sensitivity of X-ray core spectroscopy to changes in metal ligation: a systematic study of low-coordinate, high-spin ferrous complexes

In order to assess the sensitivity and complementarity of X-ray absorption and emission spectroscopies for determining changes in the metal ligation sphere, a systematic experimental and theoretical study of iron model complexes has been carried out. A series of high-spin ferrous complexes, in which the ligation sphere has been varied from a three-coordinate complex, [L(tBu)Fe(SPh)] (1) (where L(tBu) = bulky β-diketiminate ligand; SPh = phenyl thiolate) to four-coordinate complexes [L(tBu)Fe(SPh)(X)] (where X = CN(t)Bu (2); 1-methylimidazole (3); or N,N-dimethylformamide (DMF) (4)), has been investigated using a combination of Fe K-edge X-ray absorption (XAS) and Kβ X-ray emission (XES) spectroscopies. The Fe K XAS pre-edge and edge of all four complexes are consistent with a high-spin ferrous assignment, with the largest differences in the pre-edge intensities attributed to changes in covalency of the fourth coordination site. The X-ray emission spectra show pronounced changes in the valence to core region (V2C) as the identity of the coordinated ligand is varied. The experimental results have been correlated to density functional theory (DFT) calculations, to understand key molecular orbital contributions to the observed absorption and emission features. The calculations also have been extended to a series of hypothetical high-spin iron complexes to understand the sensitivity of XAS and XES techniques to different ligand protonation states ([L(tBu)Fe(II)(SPh)(NHn)](3-n) (n = 3, 2, 1, 0)), metal oxidation states [L(tBu)Fe(SPh)(N)](n-) (n = 3, 2, 1), and changes in the ligand identity [L(tBu)Fe(IV)(SPh)(X)](n-) (X = C(4-), N(3-), O(2-); n = 2, 1, 0). This study demonstrates that XAS pre-edge data have greater sensitivity to changes in oxidation state, while valence to core (V2C) XES data provide a more sensitive probe of ligand identity and protonation state. The combination of multiple X-ray spectroscopic methods with DFT results thus has the potential to provide for detailed characterization of complex inorganic systems in both chemical and biological catalysis.

New insight into the electronic structure of iron(IV)-oxo porphyrin compound I. A quantum chemical topological analysis

The electronic structure of iron-oxo porphyrin π-cation radical complex Por(·+) Fe(IV)=O (S-H) has been studied for doublet and quartet electronic states by means of two methods of the quantum chemical topology analysis: electron localization function (ELF) η(r) and electron density ρ(r). The formation of this complex leads to essential perturbation of the topological structure of the carbon-carbon bonds in porphyrin moiety. The double C=C bonds in the pyrrole anion subunits, represented by pair of bonding disynaptic basins V(i=1,2)(C,C) in isolated porphyrin, are replaced by single attractor V(C,C)(i=1-20) after complexation with the Fe cation. The iron-nitrogen bonds are covalent dative bonds, N→Fe, described by the disynaptic bonding basins V(Fe,N)(i=1-4), where electron density is almost formed by the lone pairs of the N atoms. The nature of the iron-oxygen bond predicted by the ELF topological analysis, shows a main contribution of the electrostatic interaction, Fe(δ+)···O(δ-), as long as no attractors between the C(Fe) and C(O) core basins were found, although there are common surfaces between the iron and oxygen basines and coupling between iron and oxygen lone pairs, that could be interpreted as a charge-shift bond. The Fe-S bond, characterized by the disynaptic bonding basin V(Fe,S), is partially a dative bond with the lone pair donated from sulfur atom. The change of electronic state from the doublet (M = 2) to quartet (M = 4) leads to reorganization of spin polarization, which is observed only for the porphyrin skeleton (-0.43e to 0.50e) and S-H bond (-0.55e to 0.52e).

A GGA+U approach to effective electronic correlations in thiolate-ligated iron-oxo (IV) porphyrin

High-valent oxo-metal complexes exhibit correlated electronic behavior on dense, low-lying electronic state manifolds, presenting challenging systems for electronic structure methods. Among these species, the iron-oxo (IV) porphyrin denoted Compound I occupies a privileged position, serving a broad spectrum of catalytic roles. The most reactive members of this family bear a thiolate axial ligand, exhibiting high activity toward molecular oxygen activation and substrate oxidation. The default approach to such systems has entailed the use of hybrid density functionals or multi-configurational/multireference methods to treat electronic correlation. An alternative approach is presented based on the GGA+U approximation to density functional theory, in which a generalized gradient approximation (GGA) functional is supplemented with a localization correction to treat on-site correlation as inspired by the Hubbard model. The electronic structure of thiolate-ligated iron-oxo (IV) porphyrin and corresponding Coulomb repulsion U are determined both empirically and self-consistently, yielding spin-distributions, state level splittings, and electronic densities of states consistent with prior hybrid functional calculations. Comparison of this detailed electronic structure with model Hamiltonian calculations suggests that the localized 3d iron moments induce correlation in the surrounding electron gas, strengthening local moment formation. This behavior is analogous to strongly correlated electronic systems such as Mott insulators, in which the GGA+U scheme serves as an effective single-particle representation for the full, correlated many-body problem.