Rapid Freeze-Quench ENDOR Study of Chloroperoxidase Compound I: The Site of the Radical

Axial Ligation Impedes Proton-Coupled Electron-Transfer Reactivity of a Synthetic Compound-I Analogue

The nature of the axial ligand in high-valent iron-oxo heme enzyme intermediates and related synthetic catalysts is a critical structural element for controlling proton-coupled electron-transfer (PCET) reactivity of these species. Herein, we describe the generation and characterization of three new 6-coordinate, iron(IV)-oxo porphyrinoid-π-cation-radical complexes and report their PCET reactivity together with a previously published 5-coordinate analogue, FeIV(O)(TBP8Cz+•) (TBP8Cz = octakis(p-tert-butylphenyl)corrolazinato3-) (2) (Cho, K. A high-valent iron-oxo corrolazine activates C-H bonds via hydrogen-atom transfer. J. Am. Chem. Soc. 2012, 134, 7392-7399). The new complexes FeIV(O)(TBP8Cz+•)(L) (L = 1-methyl imidazole (1-MeIm) (4a), 4-dimethylaminopyridine (DMAP) (4b), cyanide (CN-)(4c)) can be generated from either oxidation of the ferric precursors or by addition of L to the Compound-I (Cpd-I) analogue at low temperatures. These complexes were characterized by UV-vis, electron paramagnetic resonance (EPR), and Mössbauer spectroscopies, and cryospray ionization mass spectrometry (CSI-MS). Kinetic studies using 4-OMe-TEMPOH as a test substrate indicate that coordination of a sixth axial ligand dramatically lowers the PCET reactivity of the Cpd-I analogue (rates up to 7000 times slower). Extensive density functional theory (DFT) calculations together with the experimental data show that the trend in reactivity with the axial ligands does not correlate with the thermodynamic driving force for these reactions or the calculated strengths of the O-H bonds being formed in the FeIV(O-H) products, pointing to non-Bell-Evans-Polanyi behavior. However, the PCET reactivity does follow a trend with the bracketed reduction potential of Cpd-I analogues and calculated electron affinities. The combined data suggest a concerted mechanism (a concerted proton electron transfer (CPET)) and an asynchronous movement of the electron/proton pair in the transition state.

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.

Intramolecular electron transfer from biopterin to FeII-O2 complex in nitric oxide synthases occurs at very different rates between bacterial and mammalian enzymes: Direct observation of a catalytically active intermediate

Nitric oxide synthase (NOS) is a cytochrome P450-type mono‑oxygenase that catalyzes the oxidation of L-arginine to nitric oxide. We previously observed that intramolecular electron transfer from biopterin to Fe²⁺-O₂ in Deinococcus radiodurans NOS (DrNOS) using pulse radiolysis. However, the rate of electron transfer in DrNOS (2.2 × 10³ s^-1) contrasts with a reported corresponding rate (11 s^-1) in a mammalian NOS determined using rapid freeze-quench (RFQ) EPR. We applied pulse radiolysis to Bacillus subtilis NOS (bsNOS) and to rat neural NOS oxygenase domain NOS (mNOS). Concurrently, RFQ EPR was used to trap a pterin radical during single-turnover enzyme reactions of the enzymes. By using the pulse radiolysis method, hydrated electrons (e_aq^-) reduced the heme iron of NOS enzymes. Subsequently, ferrous heme reacted with O₂ to form a Fe²⁺-O₂ intermediate. In the presence of pterin, the intermediate of bsNOS was found to convert to other intermediate in the time range of milliseconds. A similar process was determined to have occurred after pulse radiolysis of the pterin-bound mNOS, though the rate was much slower. The intermediates of all of the NOS enzymes further converted to the original ferric form in the time range of seconds. When using the RFQ method, pterin radicals were formed very rapidly in both DrNOS and bsNOS in the time range of milliseconds. In contrast, the pterin radical in mNOS was observed to form slowly, at a rate of ∼20 s^-1.

17O Electron Nuclear Double Resonance Analysis of Compound I: Inverse Correlation between Oxygen Spin Population and Electron Donation

Although the activation of inert C-H bonds by metal-oxo complexes has been widely studied, important questions remain, particularly regarding the role of oxygen spin population (i.e., unpaired electrons on the oxo ligand) in facilitating C-H bond cleavage. In order to shed light on this issue, we have utilized ¹⁷O electron nuclear double resonance spectroscopy to measure the oxygen spin populations of three compound I intermediates in heme enzymes with different reactivities toward C-H bonds: chloroperoxidase, cytochrome P450, and a selenolate (selenocysteinyl)-ligated cytochrome P450. The experimental data suggest an inverse correlation between oxygen spin population and electron donation from the axial ligand. We have explored the implications of this result using a Hückel-type molecular orbital model and constrained density functional theory calculations. These investigations have allowed us to examine the relationship between oxygen spin population, oxygen charge, electron donation from the axial ligand, and reactivity.

What Drives Radical Halogenation versus Hydroxylation in Mononuclear Nonheme Iron Complexes? A Combined Experimental and Computational Study

Nonheme iron halogenases are unique enzymes in nature that selectively activate an aliphatic C–H bond of a substrate to convert it into C–X (X = Cl/Br, but not F/I). It is proposed that they generate an Fe^III(OH)(X) intermediate in their catalytic cycle. The analogous Fe^III(OH) intermediate in nonheme iron hydroxylases transfers OH^• to give alcohol product, whereas the halogenases transfer X^• to the carbon radical substrate. There remains significant debate regarding what factors control their remarkable selectivity of the halogenases. The reactivity of the complexes Fe^III(BNPA^Ph2O)(OH)(X) (X = Cl, Br) with a secondary carbon radical (R^•) is described. It is found that X^• transfer occurs with a secondary carbon radical, as opposed to OH^• transfer with tertiary radicals. Comprehensive computational studies involving density functional theory were carried out to examine the possible origins of this selectivity. The calculations reproduce the experimental findings, which indicate that halogen transfer is not observed for the tertiary radicals because of a nonproductive equilibrium that results from the endergonic nature of these reactions, despite a potentially lower reaction barrier for the halogenation pathway. In contrast, halogen transfer is favored for secondary carbon radicals, for which the halogenated product complex is thermodynamically more stable than the reactant complex. These results are rationalized by considering the relative strengths of the C–X bonds that are formed for tertiary versus secondary carbon centers. The computational analysis also shows that the reaction barrier for halogen transfer is significantly dependent on secondary coordination sphere effects, including steric and H-bonding interactions.

Radical SAM Enzyme Spore Photoproduct Lyase: Properties of the Ω Organometallic Intermediate and Identification of Stable Protein Radicals Formed during Substrate-Free Turnover

Spore photoproduct lyase is a radical S-adenosyl-l-methionine (SAM) enzyme with the unusual property that addition of SAM to the [4Fe-4S]¹⁺ enzyme absent substrate results in rapid electron transfer to SAM with accompanying homolytic S-C5' bond cleavage. Herein, we demonstrate that this unusual reaction forms the organometallic intermediate Ω in which the unique Fe atom of the [4Fe-4S] cluster is bound to C5' of the 5'-deoxyadenosyl radical (5'-dAdo^•). During catalysis, homolytic cleavage of the Fe-C5' bond liberates 5'-dAdo^• for reaction with substrate, but here, we use Ω formation without substrate to determine the thermal stability of Ω. The reaction of Geobacillus thermodenitrificans SPL (GtSPL) with SAM forms Ω within ∼15 ms after mixing. By monitoring the decay of Ω through rapid freeze-quench trapping at progressively longer times we find an ambient temperature decay time of the Ω Fe-C5' bond of τ ≈ 5-6 s, likely shortened by enzymatic activation as is the case with the Co-C5' bond of B₁₂. We have further used hand quenching at times up to 10 min, and thus with multiple SAM turnovers, to probe the fate of the 5'-dAdo^• radical liberated by Ω. In the absence of substrate, Ω undergoes low-probability conversion to a stable protein radical. The WT enzyme with valine at residue 172 accumulates a Val^•; mutation of Val172 to isoleucine or cysteine results in accumulation of an Ile^• or Cys^• radical, respectively. The structures of the radical in WT, V172I, and V172C variants have been established by detailed EPR/DFT analyses.

Does Substrate Positioning Affect the Selectivity and Reactivity in the Hectochlorin Biosynthesis Halogenase?

In this work we present the first computational study on the hectochlorin biosynthesis enzyme HctB, which is a unique three-domain halogenase that activates non-amino acid moieties tethered to an acyl-carrier, and as such may have biotechnological relevance beyond other halogenases. We use a combination of small cluster models and full enzyme structures calculated with quantum mechanics/molecular mechanics methods. Our work reveals that the reaction is initiated with a rate-determining hydrogen atom abstraction from substrate by an iron (IV)-oxo species, which creates an iron (III)-hydroxo intermediate. In a subsequent step the reaction can bifurcate to either halogenation or hydroxylation of substrate, but substrate binding and positioning drives the reaction to optimal substrate halogenation. Furthermore, several key residues in the protein have been identified for their involvement in charge-dipole interactions and induced electric field effects. In particular, two charged second coordination sphere amino acid residues (Glu223 and Arg245) appear to influence the charge density on the Cl ligand and push the mechanism toward halogenation. Our studies, therefore, conclude that nonheme iron halogenases have a chemical structure that induces an electric field on the active site that affects the halide and iron charge distributions and enable efficient halogenation. As such, HctB is intricately designed for a substrate halogenation and operates distinctly different from other nonheme iron halogenases.

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.

How the Proximal Pocket May Influence the Enantiospecificities of Chloroperoxidase-Catalyzed Epoxidations of Olefins

Chloroperoxidase-catalyzed enantiospecific epoxidations of olefins are of significant biotechnological interest. Typical enantiomeric excesses are in the range of 66%-97% and translate into free energy differences on the order of 1 kcal/mol. These differences are generally attributed to the effect of the distal pocket. In this paper, we show that the influence of the proximal pocket on the electron transfer mechanism in the rate-limiting event may be just as significant for a quantitatively accurate account of the experimentally-measured enantiospecificities.

Radical SAM catalysis via an organometallic intermediate with an Fe-[5'-C]-deoxyadenosyl bond

Radical S-adenosylmethionine (SAM) enzymes use a [4Fe-4S] cluster to cleave SAM to initiate diverse radical reactions. These reactions are thought to involve the 5'-deoxyadenosyl radical intermediate, which has not yet been detected. We used rapid freeze-quenching to trap a catalytically competent intermediate in the reaction catalyzed by the radical SAM enzyme pyruvate formate-lyase activating enzyme. Characterization of the intermediate by electron paramagnetic resonance and (13)C, (57)Fe electron nuclear double-resonance spectroscopies reveals that it contains an organometallic center in which the 5' carbon of a SAM-derived deoxyadenosyl moiety forms a bond with the unique iron site of the [4Fe-4S] cluster. Discovery of this intermediate extends the list of enzymatic bioorganometallic centers to the radical SAM enzymes, the largest enzyme superfamily known, and reveals intriguing parallels to B12 radical enzymes.

Rapid Freeze-Quench EPR Spectroscopy: Improved Collection of Frozen Particles

Rapid freeze-quench (RFQ) in combination with electron paramagnetic resonance (EPR) spectroscopy at X-band is a proven technique to trap and characterize paramagnetic intermediates of biochemical reactions. Preparation of suitable samples is still cumbersome, despite many attempts to remedy this problem, and limits the wide applicability of RFQ EPR. We present a method, which improves the collection of freeze-quench particles from isopentane and their packing in an EPR tube. The method is based on sucking the particle suspension into an EPR tube with a filter at the bottom. This procedure results in a significant reduction of the required volume of reactants, which allows the economical use of valuable reactants such as proteins. The approach also enables the successful collection of smaller frozen particles, which are generated at higher flow rates. The method provides for a reproducible, efficient and fast collection of the freeze-quench particles and can be easily adapted to RFQ EPR at higher microwave frequencies than X-band.

Chloroperoxidase-Catalyzed Epoxidation of Cis-β-Methylstyrene: NH-S Hydrogen Bonds and Proximal Helix Dipole Change the Catalytic Mechanism and Significantly Lower the Reaction Barrier

Proximal hydrogen bonding of the axial sulfur with the backbone amides (NH-S) is a conserved feature of heme-thiolate enzymes such as chloroperoxidase (CPO) and cytochrome P450 (P450). In CPO, the effect of NH-S bonds is amplified by the dipole moment of the proximal helix. Our gas-phase DFT studies show that the proximal pocket effect significantly enhances CPO's reactivity toward the epoxidation of olefinic substrates. Comparison of models with and without proximal pocket residues shows that with them, the barrier for Cβ-O bond formation is lowered by about ∼4.6 kcal/mol, while Cα-O-Cβ ring closure becomes barrierless. The dipole moment of the proximal helix was estimated to contribute 1/3 of the decrease, while the rest is attributed to the effect of NH-S bonds. The decrease of the reaction barrier correlates with increased electron density transfer to residues of the proximal pocket. The effect is most pronounced on the doublet spin surface and involves a change in the electron-transfer mechanism. A full enzyme QMMM study on the doublet spin surface gives about the same barrier as the gas-phase DFT study. The free-energy barrier was estimated to be in agreement with the experimental results for the CPO-catalyzed epoxidation of styrene.

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.

Spectroscopic studies reveal that the heme regulatory motifs of heme oxygenase-2 are dynamically disordered and exhibit redox-dependent interaction with heme

Heme oxygenase (HO) catalyzes a key step in heme homeostasis: the O2- and NADPH-cytochrome P450 reductase-dependent conversion of heme to biliverdin, Fe, and CO through a process in which the heme participates both as a prosthetic group and as a substrate. Mammals contain two isoforms of this enzyme, HO2 and HO1, which share the same α-helical fold forming the catalytic core and heme binding site, as well as a membrane spanning helix at their C-termini. However, unlike HO1, HO2 has an additional 30-residue N-terminus as well as two cysteine-proline sequences near the C-terminus that reside in heme regulatory motifs (HRMs). While the role of the additional N-terminal residues of HO2 is not yet understood, the HRMs have been proposed to reversibly form a thiol/disulfide redox switch that modulates the affinity of HO2 for ferric heme as a function of cellular redox poise. To further define the roles of the N- and C-terminal regions unique to HO2, we used multiple spectroscopic techniques to characterize these regions of the human HO2. Nuclear magnetic resonance spectroscopic experiments with HO2 demonstrate that, when the HRMs are in the oxidized state (HO2(O)), both the extra N-terminal and the C-terminal HRM-containing regions are disordered. However, protein NMR experiments illustrate that, under reducing conditions, the C-terminal region gains some structure as the Cys residues in the HRMs undergo reduction (HO2(R)) and, in experiments employing a diamagnetic protoporphyrin, suggest a redox-dependent interaction between the core and the HRM domains. Further, electron nuclear double resonance and X-ray absorption spectroscopic studies demonstrate that, upon reduction of the HRMs to the sulfhydryl form, a cysteine residue from the HRM region ligates to a ferric heme. Taken together with EPR measurements, which show the appearance of a new low-spin heme signal in reduced HO2, it appears that a cysteine residue(s) in the HRMs directly interacts with a second bound heme.

Ab initio dynamics of the cytochrome P450 hydroxylation reaction

The iron(IV)-oxo porphyrin π-cation radical known as Compound I is the primary oxidant within the cytochromes P450, allowing these enzymes to affect the substrate hydroxylation. In the course of this reaction, a hydrogen atom is abstracted from the substrate to generate hydroxyiron(IV) porphyrin and a substrate-centered radical. The hydroxy radical then rebounds from the iron to the substrate, yielding the hydroxylated product. While Compound I has succumbed to theoretical and spectroscopic characterization, the associated hydroxyiron species is elusive as a consequence of its very short lifetime, for which there are no quantitative estimates. To ascertain the physical mechanism underlying substrate hydroxylation and probe this timescale, ab initio molecular dynamics simulations and free energy calculations are performed for a model of Compound I catalysis. Semiclassical estimates based on these calculations reveal the hydrogen atom abstraction step to be extremely fast, kinetically comparable to enzymes such as carbonic anhydrase. Using an ensemble of ab initio simulations, the resultant hydroxyiron species is found to have a similarly short lifetime, ranging between 300 fs and 3600 fs, putatively depending on the enzyme active site architecture. The addition of tunneling corrections to these rates suggests a strong contribution from nuclear quantum effects, which should accelerate every step of substrate hydroxylation by an order of magnitude. These observations have strong implications for the detection of individual hydroxylation intermediates during P450 catalysis.

ATPase site configuration of the RNA helicase DbpA probed by ENDOR spectroscopy

Electron-nuclear double resonance (ENDOR) is a method that probes the local structure of paramagnetic centers via their hyperfine interactions with nearby magnetic nuclei. Here we describe the use of this technique to structurally characterize the ATPase active site of the RNA helicase DbpA, where Mg(2+)-ATP binds. This is achieved by substituting the EPR (electron paramagnetic resonance) silent Mg(2+) ion with paramagnetic, EPR active, Mn(2+) ion. (31)P ENDOR provides the interaction of the Mn(2+) with the nucleotide (ADP, ATP and its analogs) through the phosphates. The ENDOR spectra clearly distinguish between ATP- and ADP-binding modes. In addition, by preparing (13)C-enriched DbpA, (13)C ENDOR is used to probe the interaction of the Mn(2+) with protein residues. This combination allows tracking structural changes in the Mn(2+) coordination shell, in the ATPase site, in different states of the protein, namely with and without RNA and with different ATP analogs. Here, a detailed description of sample preparation and the ENDOR measurement methodology is provided, focusing on measurements at W-band (95 GHz) where sensitivity is high and spectral interpretations are relatively simple.

Electron paramagnetic resonance and electron-nuclear double resonance studies of the reactions of cryogenerated hydroperoxoferric-hemoprotein intermediates

The fleeting ferric peroxo and hydroperoxo intermediates of dioxygen activation by hemoproteins can be readily trapped and characterized during cryoradiolytic reduction of ferrous hemoprotein–O₂ complexes at 77 K. Previous cryoannealing studies suggested that the relaxation of cryogenerated hydroperoxoferric intermediates of myoglobin (Mb), hemoglobin, and horseradish peroxidase (HRP), either trapped directly at 77 K or generated by cryoannealing of a trapped peroxo-ferric state, proceeds through dissociation of bound H₂O₂ and formation of the ferric heme without formation of the ferryl porphyrin π-cation radical intermediate, compound I (Cpd I). Herein we have reinvestigated the mechanism of decays of the cryogenerated hydroperoxyferric intermediates of α- and β-chains of human hemoglobin, HRP, and chloroperoxidase (CPO). The latter two proteins are well-known to form spectroscopically detectable quasistable Cpds I. Peroxoferric intermediates are trapped during 77 K cryoreduction of oxy Mb, α-chains, and β-chains of human hemoglobin and CPO. They convert into hydroperoxoferric intermediates during annealing at temperatures above 160 K. The hydroperoxoferric intermediate of HRP is trapped directly at 77 K. All studied hydroperoxoferric intermediates decay with measurable rates at temperatures above 170 K with appreciable solvent kinetic isotope effects. The hydroperoxoferric intermediate of β-chains converts to the S = 3/2 Cpd I, which in turn decays to an electron paramagnetic resonance (EPR)-silent product at temperature above 220 K. For all the other hemoproteins studied, cryoannealing of the hydroperoxo intermediate directly yields an EPR-silent majority product. In each case, a second follow-up 77 K γ-irradiation of the annealed samples yields low-spin EPR signals characteristic of cryoreduced ferrylheme (compound II, Cpd II). This indicates that in general the hydroperoxoferric intermediates relax to Cpd I during cryoanealing at low temperatures, but when this state is not captured by reaction with a bound substrate, it is reduced to Cpd II by redox-active products of radiolysis.

Multiple, disparate redox pathways exhibited by a tris(pyrrolido)ethane iron complex

Iron(III) complexes of the tris(pyrrolide)ethane trianion have been synthesized by reaction of one- and two-electron oxidants with [(tpe)Fe(THF)][Li(THF)4] (tpe = tris(5-mesitylpyrrolyl)ethane). X-ray crystallography, (57)Fe Mössbauer, (1)H NMR and EPR spectroscopy, SQUID magnetometry, and density functional theory calculations were employed to rigorously establish the iron 3+ oxidation state. All oxidants employed are proposed to operate via an inner-sphere electron transfer mechanism. Dialkyl peroxides and dibenzyldisulfide served to oxidize iron by one electron, and group transfer of an aryl nitrene unit to the Fe(2+) starting material resulted in formation of Fe(3+) amido species following H-atom abstraction by a presumed nitrenoid intermediate. Single electron transfer to and from diphenyldiazoalkane was also observed to yield a diphenyldiazomethanyl radical anion antiferromagnetically coupled to the S = 5/2 Fe(3+). Isolation of Fe(3+) complexes of tpe, in comparison with previous results wherein the tpe ligand was the redox active moiety, presents an unusual juxtaposition of two noncommunicating redox reservoirs, each accessible via different reaction pathways (namely, inner- and outer-sphere electron transfer).

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.

A novel microfluidic rapid freeze-quench device for trapping reactions intermediates for high field EPR analysis

Rapid freeze quench electron paramagnetic resonance (RFQ)-EPR is a method for trapping short lived intermediates in chemical reactions and subjecting them to EPR spectroscopy investigation for their characterization. Two (or more) reacting components are mixed at room temperature and after some delay the mixture is sprayed into a cold trap and transferred into the EPR tube. A major caveat in using commercial RFQ-EPR for high field EPR applications is the relatively large amount of sample needed for each time point, a major part of which is wasted as the dead volume of the instrument. The small sample volume (∼2μl) needed for high field EPR spectrometers, such as W-band (∼3.5T, 95GHz), that use cavities calls for the development of a microfluidic based RFQ-EPR apparatus. This is particularly important for biological applications because of the difficulties often encountered in producing large amounts of intrinsically paramagnetic proteins and spin labeled nucleic acid and proteins. Here we describe a dedicated microfluidic based RFQ-EPR apparatus suitable for small volume samples in the range of a few μl. The device is based on a previously published microfluidic mixer and features a new ejection mechanism and a novel cold trap that allows collection of a series of different time points in one continuous experiment. The reduction of a nitroxide radical with dithionite, employing the signal of Mn(2+) as an internal standard was used to demonstrate the performance of the microfluidic RFQ apparatus.

Reaction of ferric Caldariomyces fumago chloroperoxidase with meta-chloroperoxybenzoic acid: sequential formation of compound I, compound II and regeneration of the ferric state using one reactant

The mechanism of the reaction between ferric Caldariomyces fumago chloroperoxidase (CCPO) and meta-chloroperoxybenzoic acid (mCPBA) has been examined. It has previously been established that an Fe(IV) -oxo porphyrin radical species known as Compound I (Cpd I) is formed by two-electron oxidation of the native ferric enzyme by a variety of oxidants including organic peracids and hydroperoxides. Cpd I can return to the ferric state either by direct oxygen insertion into an organic substrate (e.g. a P450 oxygenase-like reaction), or by two consecutive one-electron additions, the first resulting in an intermediate Fe(IV) -oxo species known as Compound II (Cpd II). There has been much debate over the role of Cpd II and the details of its structure. In the present study, both CCPO Fe(IV) -oxo intermediates are formed, but unlike most CCPO reactions, Cpd I and Cpd II are formed using the same reactant, mCPBA. Thus, the peracid is used as an oxo donor to produce Cpd I and then as a reductant to reduce Cpd I to Cpd II, and finally, Cpd II to the ferric state. The observation of saturation kinetics with respect to mCPBA concentration for each step is consistent with the formation of CCPO-mCPBA complexes in each phase of the reaction. The original reaction mechanism for ferric CCPO with mCPBA was hypothesized to involve a scrambling mechanism with a unique Fe -OOO-C(O)R intermediate formed with no observed Cpd II intermediate. The data reported herein clearly demonstrate the formation of Cpd II in returning the oxidized enzyme back to its native ferric state.

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.

Chloroperoxidase-catalyzed epoxidation of cis-β-methylstyrene: distal pocket flexibility tunes catalytic reactivity

Chloroperoxidase, the most versatile heme protein, has a hybrid active site pocket that shares structural features with peroxidases and cytochrome P450s. The simulation studies presented here show that the enzyme possesses a remarkable ability to efficiently utilize its hybrid structure, assuming structurally different peroxidase-like and P450-like distal pocket faces and thereby enhancing the inherent catalytic capability of the active center. We find that, during epoxidation of cis-β-methylstyrene (CBMS), the native peroxidase-like aspect of the distal pocket is diminished as the polar Glu183 side chain is displaced away from the active center and the distal pocket takes on a more hydrophobic, P450-like, aspect. The P450-like distal pocket provides a significant enthalpic stabilization of ∼4 kcal/mol of the 14 kcal/mol reaction barrier for gas-phase epoxidation of CMBS by an oxyferryl heme-thiolate species. This stabilization comes from breathing of the distal pocket. As until recently the active site of chloroperoxidase was postulated to be inflexible, these results suggest a new conceptual understanding of the enzyme's versatility: catalytic reactivity is tuned by flexibility of the distal pocket.

Catalases versus peroxidases: DFT investigation of H₂O₂ oxidation in models systems and implications for heme protein engineering

Catalases and peroxidases are ubiquitous heme enzymes that catalyze the removal of hydrogen peroxide (H(2)O(2)). Both enzymes use one molecule of hydrogen peroxide to form a high valent iron intermediate named Compound I (Cpd I). However, whereas catalase Cpd I oxidizes a second H(2)O(2) molecule to oxygen, peroxidases use this intermediate to oxidize other substrates rather than H(2)O(2). The origin of the different reactivity of peroxidases and catalases is not known, but it is likely to be related to structural differences between the two heme active sites. Recent modeling studies suggest that the oxidation of H(2)O(2) by catalase Cpd I may take place by two hydrogen atom transfer steps. In this work, we investigate how catalases and peroxidases compare along the same hydrogen transfer steps to give hints into the question why peroxidases cannot efficiently oxidize H(2)O(2). The use of simplified models allows us to probe the direct effect of the proximal ligand (tyrosinate in catalases and histidine in peroxidases) without masking from the protein environment. We show that the nature of the fifth ligand (His in peroxidase and Tyr in catalase) has little effect on the energy barriers of the hydrogen transfer steps. On the contrary, the Cpd I-hydrogen peroxide (O(Fe)-O(peroxide)) distance affects significantly the reaction barriers. We propose that the distal side architecture of peroxidases do not allow to attain short O(Cpd I)-O(peroxide) distances, thus resulting in a lower efficiency towards H(2)O(2) oxidation.

Enantiospecificity of chloroperoxidase-catalyzed epoxidation: biased molecular dynamics study of a cis-β-methylstyrene/chloroperoxidase-compound I complex

Molecular dynamics simulations of an explicitly solvated cis-β-methylstyrene/chloroperoxidase-Compound I complex are performed to determine the cause of the high enantiospecificity of epoxidation. From the simulations, a two-dimensional free energy potential is calculated to distinguish binding potential wells from which reaction to 1S2R and 1R2S epoxide products may occur. Convergence of the free energy potential is accelerated with an adaptive biasing potential. Analysis of binding is followed by analysis of 1S2R and 1R2S reaction precursor structures in which the substrate, having left the binding wells, places its reactive double bond in steric proximity to the oxyferryl heme center. Structural analysis of binding and reaction precursor conformations is presented. We find that 1), a distortion of Glu(183) is important for CPO-catalyzed epoxidation as was postulated previously based on experimental results; 2), the free energy of binding does not provide significant differentiation between structures leading to the respective epoxide enantiomers; and 3), CPO's enantiospecificity toward cis-β-methylstyrene is likely to be caused by a specific group of residues which form a hydrophobic core surrounding the oxyferryl heme center.

Design, Implementation, Simulation, and Visualization of a Highly Efficient RIM Microfluidic Mixer for Rapid Freeze-Quench of Biological Samples

Rapid freeze-quench (RFQ) trapping of short-lived reaction intermediates for spectroscopic study plays an important role in the characterization of biological reactions. Recently there has been considerable effort to achieve submillisecond reaction deadtimes. We present here a new, robust, high-velocity microfluidic mixer that enables such rapid freeze-quenching. It is a based on the mixing method of two impinging jets commonly used in reaction injection molding (RIM) of plastics. This method achieves efficient mixing by inducing chaotic flow at relatively low Reynolds numbers (Re =140). We present the first mathematical simulation and microscopic visualization of mixing in such RFQ micromixers, the results of which show that the impinging solutions efficiently mix within the mixing chamber. These tests, along with a practical demonstration in a RFQ setup that involves copper wheels, show this new mixer can in practice provide reaction deadtimes as low as 100 microseconds.

Spectroscopic features of cytochrome P450 reaction intermediates

Cytochromes P450 constitute a broad class of heme monooxygenase enzymes with more than 11,500 isozymes which have been identified in organisms from all biological kingdoms [1]. These enzymes are responsible for catalyzing dozens chemical oxidative transformations such as hydroxylation, epoxidation, N-demethylation, etc., with very broad range of substrates [2,3]. Historically these enzymes received their name from 'pigment 450' due to the unusual position of the Soret band in UV-vis absorption spectra of the reduced CO-saturated state [4,5]. Despite detailed biochemical characterization of many isozymes, as well as later discoveries of other 'P450-like heme enzymes' such as nitric oxide synthase and chloroperoxidase, the phenomenological term 'cytochrome P450' is still commonly used as indicating an essential spectroscopic feature of the functionally active protein which is now known to be due to the presence of a thiolate ligand to the heme iron [6]. Heme proteins with an imidazole ligand such as myoglobin and hemoglobin as well as an inactive form of P450 are characterized by Soret maxima at 420nm [7]. This historical perspective highlights the importance of spectroscopic methods for biochemical studies in general, and especially for heme enzymes, where the presence of the heme iron and porphyrin macrocycle provides rich variety of specific spectroscopic markers available for monitoring chemical transformations and transitions between active intermediates of catalytic cycle.

Cytochrome P450 compound I: capture, characterization, and C-H bond activation kinetics

Cytochrome P450 enzymes are responsible for the phase I metabolism of approximately 75% of known pharmaceuticals. P450s perform this and other important biological functions through the controlled activation of C-H bonds. Here, we report the spectroscopic and kinetic characterization of the long-sought principal intermediate involved in this process, P450 compound I (P450-I), which we prepared in approximately 75% yield by reacting ferric CYP119 with m-chloroperbenzoic acid. The Mössbauer spectrum of CYP119-I is similar to that of chloroperoxidase compound I, although its electron paramagnetic resonance spectrum reflects an increase in |J|/D, the ratio of the exchange coupling to the zero-field splitting. CYP119-I hydroxylates the unactivated C-H bonds of lauric acid [D(C-H) ~ 100 kilocalories per mole], with an apparent second-order rate constant of k(app) = 1.1 × 10(7) per molar per second at 4°C. Direct measurements put a lower limit of k ≥ 210 per second on the rate constant for bound substrate oxidation, whereas analyses involving kinetic isotope effects predict a value in excess of 1400 per second.

What factors influence the rate constant of substrate epoxidation by compound I of cytochrome P450 and analogous iron(IV)-oxo oxidants?

The cytochromes P450 are a versatile range of mono-oxygenase enzymes that catalyze a variety of different chemical reactions, of which the key reactions include aliphatic hydroxylation and C=C double bond epoxidation. To establish the fundamental factors that govern substrate epoxidation by these enzymes we have done a systematic density functional theory study on substrate epoxidation by the active species of P450 enzymes, namely the iron(IV)-oxo porphyrin cation radical oxidant or Compound I. We show here, for the first time, that the rate constant of substrate epoxidation, and hence the activation energy, correlates with the ionization potential of the substrate as well as with intrinsic electronic properties of the active oxidant such as the polarizability volume. To explain these findings we present an electron-transfer model for the reaction mechanism that explains the factors that determine the barrier heights and developed a valence bond (VB) curve crossing mechanism to rationalize the observed trends. In addition, we have found a correlation for substrate epoxidation reactions catalyzed by a range of heme and nonheme iron(IV)-oxo oxidants with the strength of the O-H bond in the iron-hydroxo complex, i.e. BDE(OH), which is supported by the VB model. Finally, the fundamental factors that determine the regioselectivity change between substrate hydroxylation and epoxidation are discussed. It is shown that the regioselectivity of aliphatic hydroxylation versus double bond epoxidation is not influenced by the choice of the oxidant but is purely substrate dependent.

Effects of substrate, protein environment, and proximal ligand mutation on compound I and compound 0 of chloroperoxidase

This paper investigates the enzyme chloroperoxidase (CPO) by means of hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. The effects of anionic substrate, protein environment, and proximal ligand mutation on the high-valent iron-oxo species, compound I (Cpd I), and the ferric hydroperoxide complex, compound 0 (Cpd 0), are studied. The results indicate that the presence of an anionic substrate (acetate) and the protonation state of one critical residue (Glu104) have a considerable impact on the relative stabilities of Cpd I and Cpd 0. In the absence of the substrate or when the substrate is protonated, Cpd I is considerably more stable, and its formation barrier is smaller than in the case where the substrate is in its anionic state and when Glu104 is deprotonated. This trend, which is shown to be a simple manifestation of the Hammond principle, reproduces the experimental observation that the working pH of the enzyme is acidic. Furthermore, in the absence of substrate (or when it is protonated), the relative Cpd 0/Cpd I energies are found to be a good index of Cpd I stability in heme enzymes and to follow the experimental order: horseradish peroxidase (HRP) > CPO > P450. In silico mutation of the proximal ligand from cysteine to selenocysteine was found to have no effect at all on the properties of Cpd I (e.g., spin density on the chalcogen, Mössbauer parameters, etc.) and its relative stability to Cpd 0 or on the corresponding barrier for formation. This surprising finding shows that the polar CPO pocket applies a leveling effect that stabilizes the anionic forms of the proximal ligands (CysS(-) and CysSe(-)). This in turn means that the Se-Cpd I of the mutant CPO is observable.