Two-State Reactivity, Electromerism, Tautomerism, and “Surprise” Isomers in the Formation of Compound II of the Enzyme Horseradish Peroxidase from the Principal Species, Compound I

Impact of an Ionic Liquid Solution on Horseradish Peroxidase Activity

Horseradish peroxidase (HRP) is an enzyme that oxidizes pollutants from wastewater. A previous report indicated that peroxidases can have an enhancement in initial enzymatic activity in an aqueous solution of 0.26 M 1-ethyl-3-methylimidazolium ethyl sulfate ([EMIm][EtSO4]) at neutral pH. However, the atomistic details remain elusive. In the enzymatic landscape of HRP, compound II (Cpd II) plays a key role and involves a histidine (H42) residue. Cpd II exists as oxoferryl (2a) or hydroxoferryl (2b(FeIV)) forms, where 2a is the predominantly observed form in experimental studies. Intriguingly, the ferric 2b(FeIII) form seen in synthetic complexes has not been observed in HRP. Here, we have investigated the structure and dynamics of HRP in pure water and aqueous [EMIm][EtSO4] (0.26 M), as well as the reaction mechanism of 2a to 2b conversion using polarizable molecular dynamics (MD) simulations and quantum mechanics/molecular mechanics (QM/MM) calculations. When HRP is solvated in aq [EMIm][EtSO4], the catalytic water displaces, and H42 directly orients over the ferryl moiety, allowing a direct proton transfer (PT) with a significant energy barrier reduction. Conversely, in neat water, the reaction of 2a to 2b follows the previously reported mechanism. We further investigated the deprotonated form of H42. Analysis of the electric fields at the active site indicates that the aq [EMIm][EtSO4] medium facilitates the reaction by providing a more favorable environment compared with the system solvated in neat water. Overall, the atomic level supports the previous experimental observations and underscores the importance of favorable electric fields in the active site to promote catalysis.

Ferric heme b in aqueous micellar and vesicular systems: state-of-the-art and challenges

Ferric heme b (= ferric protoporphyrin IX = hemin) is an important prosthetic group of different types of enzymes, including the intensively investigated and widely applied horseradish peroxidase (HRP). In HRP, hemin is present in monomeric form in a hydrophobic pocket containing among other amino acid side chains the two imidazoyl groups of His170 and His42. Both amino acids are important for the peroxidase activity of HRP as an axial ligand of hemin (proximal His170) and as an acid/base catalyst (distal His42). A key feature of the peroxidase mechanism of HRP is the initial formation of compound I under heterolytic cleavage of added hydrogen peroxide as a terminal oxidant. Investigations of free hemin dispersed in aqueous solution showed that different types of hemin dimers can form, depending on the experimental conditions, possibly resulting in hemin crystallization. Although it has been recognized already in the 1970s that hemin aggregation can be prevented in aqueous solution by using micelle-forming amphiphiles, it remains a challenge to prepare hemin-containing micellar and vesicular systems with peroxidase-like activities. Such systems are of interest as cheap HRP-mimicking catalysts for analytical and synthetic applications. Some of the key concepts on which research in this fascinating and interdisciplinary field is based are summarized, along with major accomplishments and possible directions for further improvement. A systematic analysis of the physico-chemical properties of hemin in aqueous micellar solutions and vesicular dispersions must be combined with a reliable evaluation of its catalytic activity. Future studies should show how well the molecular complexity around hemin in HRP can be mimicked by using micelles or vesicles. Because of the importance of heme b in virtually all biological systems and the fact that porphyrins and hemes can be obtained under potentially prebiotic conditions, ideas exist about the possible role of heme-containing micellar and vesicular systems in prebiotic times.

Comparative oxidative ability of mononuclear and dinuclear high-valent iron-oxo species towards the activation of methane: does the axial/bridge atom modulate the reactivity?

Over the years, mononuclear Fe^IVO species have been extensively studied, but the presence of dinuclear Fe^IVO species in soluble methane monooxygenase (sMMO) has inspired the development of biomimic models that could activate inert substrates such as methane. There are some successful attempts; particularly the [(Por)(m-CBA) Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- species has been reported to activate methane and yield decent catalytic turnover numbers and therefore regarded as the closest to the sMMO enzyme functional model, as no mononuclear Fe^IVO analogues could achieve this feat. In this work, we have studied a series of mono and dinuclear models using DFT and ab initio DLPNO-CCSD(T) calculations to probe the importance of nuclearity in enhancing the reactivity. We have probed the catalytic activities of four complexes: [(HO)Fe^IV(O)(Por)]^- (1), [(HO)Fe^IV(O)(Por˙⁺)] (2), μ-oxo dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-O)Fe^IV(O) (Por˙⁺)]^- (3) and N-bridged dinuclear iron species [(Por)(m-CBA)Fe^IV(μ-N)Fe^IV(O)(Por˙⁺)]^- (4) towards the activation of methane. Additionally, calculations were performed on the mononuclear models [(X)Fe^IV(O)(Por˙⁺)]ⁿ {X = N 4a (n = -2), NH 4b (n = -1) and NH₂4c (n = 0)} to understand the role of nuclearity in the reactivity. DFT calculations performed on species 1-4 suggest an interesting variation among them, with species 1-3 possessing an intermediate spin (S = 1) as a ground state and species 4 possessing a high-spin (S = 2) as a ground state. Furthermore, the two Fe^IV centres in species 3 and 4 are antiferromagnetically coupled, yielding a singlet state with a distinct difference in their electronic structure. On the other hand, species 2 exhibits a ferromagnetic coupling between the Fe^IV and the Por˙⁺ moiety. Our calculations suggest that the higher barriers for the C-H bond activation of methane and the rebound step for species 1 and 3 are very high in energy, rendering them unreactive towards methane, while species 2 and 4 have lower barriers, suggesting their reactivity towards methane. Studies on the system reveal that model 4a has multiple FeN bonds facilitating greater reactivity, whereas the other two models have longer Fe-N bonds and less radical character with steeper barriers. Strong electronic cooperativity is found to be facilitated by the bridging nitride atom, and this cooperativity is suppressed by substituents such as oxygen, rendering them inactive. Thus, our study unravels that apart from enhancing the nuclearity, bridging atoms that facilitate strong cooperation between the metals are required to activate very inert substrates such as methane, and our results are broadly in agreement with earlier experimental findings.

Hemin-catalyzed oxidative oligomerization of p-aminodiphenylamine (PADPA) in the presence of aqueous sodium dodecylbenzenesulfonate (SDBS) micelles

In a previous report on the enzymatic synthesis of the conductive emeraldine salt form of polyaniline (PANI-ES) in aqueous solution using PADPA (p-aminodiphenylamine) as monomer, horseradish peroxidase isoenzyme C (HRPC) was applied as a catalyst at pH = 4.3 with H₂O₂ as a terminal oxidant. In that work, anionic vesicles were added to the reaction mixture for (i) guiding the reaction to obtain poly(PADPA) products that resemble PANI-ES, and for (ii) preventing product precipitation (known as the “template effect”). In the work now presented, instead of native HRPC, only its prosthetic group ferric heme b (= hemin) was utilized as a catalyst, and micelles formed from SDBS (sodium dodecylbenzenesulfonate) served as templates. For the elaborated optimal reaction conditions, complementary UV/vis/NIR, EPR, and Raman spectroscopy measurements clearly showed that the reaction mixture obtained after completion of the reaction contained PANI-ES-like products as dominating species, very similar to the products formed with HRPC as catalyst. HEPES (4-(2-hydroxyethyl)-1-piperazineethanesulfonate) was found to have a positive effect on the reaction rate as compared to dihydrogenphosphate. This work is the first on the template-assisted formation of PANI-ES type products under mild, environmentally friendly conditions using hemin as a cost-effective catalyst.

Polyaniline emeraldine salt-type products were synthesized under mild, environmentally friendly conditions using hemin as a cost-effective catalyst, p-aminodiphenylamine (PADPA) as a monomer, and micelles formed from SDBS as templates.

Measurement of time dependent product branching ratios indicates two-state reactivity in metal mediated chemical reactions

For several decades, the influence of Two State Reactivity (TSR) has been implicated in a host of reactions, but has lacked a stand-alone, definitive experimental kinetic signature identifying its occurrence. Here, we demonstrate that the measurement of a temporally dependent product branching ratio is indicative of spin inversion and is a kinetic signature of TSR. This is caused by products exiting different hypersurfaces with different rates and relative exothermicities. The composite measurement of product intensities with the same mass but with different multiplicities yield biexponential temporal dependences with the sampled product ratio changing in time. These measurements are made using the single photon initiated dissociative rearrangement reaction (SPIDRR) technique which identifies TSR but further determines the kinetic parameters for reaction along the original ground electronic surface in competition with spin inversion and its consequent TSR.

Histidine versus Cysteine-Bearing Heme-Dependent Halogen Peroxidases: Parallels and Differences for Cl- Oxidation

The homodimeric myeloperoxidase (MPO) features a histidine as a proximal ligand and a sulfonium linkage covalently attaching the heme porphyrin ring to the protein. MPO is able to catalyze Cl^- oxidation with about the same efficiency as chloroperoxidase at pH 7.0. In this study, we seek to explore the parallels and differences between the histidine and cysteine heme-dependent halogen peroxidases. Transition states, reaction barriers, and relevant thermodynamic properties are calculated on protein models. Together with electronic structure calculations, it gives an overview of the reaction mechanisms and of the factors that determine the selectivity between one- and two-electron paths. Conclusions point to the innate oxidizing nature of MPO with the ester and sulfonium linkages hiking up the reactivity to enable chloride oxidation. The installation of a deprotonated imidazolate as a proximal ligand does not shift the equilibrium from one- to two-electron events without influencing the chemistry of the oxidation reaction.

Metal/Metal Redox Isomerism Governed by Configuration

A pair of diastereomeric dinuclear complexes, [Tp′(CO)BrW{μ‐η²‐C,C′‐κ²‐S,P‐C₂(PPh₂)S}Ru(η⁵‐C₅H₅)(PPh₃)], in which W and Ru are bridged by a phosphinyl(thiolato)alkyne in a side‐on carbon P,S‐chelate coordination mode, were synthesized, separated and fully characterized. Even though the isomers are similar in their spectroscopic properties and redox potentials, the like‐isomer is oxidized at W while the unlike‐isomer is oxidized at Ru, which is proven by IR, NIR and EPR‐spectroscopy supported by spectro‐electrochemistry and computational methods. The second oxidation of the complexes was shown to take place at the metal left unaffected in the first redox step. Finally, the tipping point could be realized in the unlike isomer of the electronically tuned thiophenolate congener [Tp′(CO)(PhS)W{μ‐η²‐C,C′‐κ²‐S,P‐C₂(PPh₂)S}Ru(η⁵‐C₅H₅)‐(PPh₃)], in which valence trapped W^III/Ru^II and W^II/Ru^III cationic species are at equilibrium.

Computational Studies on the Mechanism and Origin of the Different Regioselectivities of Manganese Porphyrin-Catalyzed C-H Bond Hydroxylation and Amidation of Equilenin Acetate

The manganese porphyrin-catalyzed C-H bond hydroxylation and amidation of equilenin acetate developed by Breslow and his co-worker have been investigated with density functional theory (DFT) calculations. The hydroxylation of C(sp²)-H bond of equilenin acetate leading to the 6-hydroxylated product is more favorable than the hydroxylation of C(sp³)-H bond of equilenin acetate, leading to the 11β-hydroxylation product. The computational results suggest that the C(sp²)-H bond hydroxylation of equilenin acetate undergoes an oxygen-atom-transfer mechanism, which is more favorable than the C(sp³)-H bond hydroxylation undergoing the hydrogen-atom-abstraction/oxygen-rebound (HAA/OR) mechanism by 1.6 kcal/mol. That is why, the 6-hydroxylated product is the major product and the 11β-hydroxylated product is the minor product. In contrast, the 11β-amidated product is the only observed product in manganese porphyrin-catalyzed amidation reaction. The benzylic amidation undergoes a hydrogen-atom-abstraction/nitrogen-rebound (HAA/NR) mechanism, in which hydrogen atom abstraction is followed by nitrogen rebound, leading to the 11β-amidated product. The benzylic C(sp³)-H bond amidation at the C-11 position is more favorable than aromatic amidation at the C-6 position by 4.9 kcal/mol. Therefore, the DFT computational results are consistent with the experiments that manganese porphyrin-catalyzed C-H bond hydroxylation and amidation of equilenin acetate have different regioselectivities.

How Does Replacement of the Axial Histidine Ligand in Cytochrome c Peroxidase by Nδ-Methyl Histidine Affect Its Properties and Functions? A Computational Study

Heme peroxidases have important functions in nature related to the detoxification of H₂O₂. They generally undergo a catalytic cycle where, in the first stage, the iron(III)-heme-H₂O₂ complex is converted into an iron(IV)-oxo-heme cation radical species called Compound I. Cytochrome c peroxidase Compound I has a unique electronic configuration among heme enzymes where a metal-based biradical is coupled to a protein radical on a nearby Trp residue. Recent work using the engineered N_δ-methyl histidine-ligated cytochrome c peroxidase highlighted changes in spectroscopic and catalytic properties upon axial ligand substitution. To understand the axial ligand effect on structure and reactivity of peroxidases and their axially N_δ-methyl histidine engineered forms, we did a computational study. We created active site cluster models of various sizes as mimics of horseradish peroxidase and cytochrome c peroxidase Compound I. Subsequently, we performed density functional theory studies on the structure and reactivity of these complexes with a model substrate (styrene). Thus, the work shows that the N_δ-methyl histidine group has little effect on the electronic configuration and structure of Compound I and little changes in bond lengths and the same orbital occupation is obtained. However, the N_δ-methyl histidine modification impacts electron transfer processes due to a change in the reduction potential and thereby influences reactivity patterns for oxygen atom transfer. As such, the substitution of the axial histidine by N_δ-methyl histidine in peroxidases slows down oxygen atom transfer to substrates and makes Compound I a weaker oxidant. These studies are in line with experimental work on N_δ-methyl histidine-ligated cytochrome c peroxidases and highlight how the hydrogen bonding network in the second coordination sphere has a major impact on the function and properties of the enzyme.

Consecutive Marcus Electron and Proton Transfer in Heme Peroxidase Compound II-Catalysed Oxidation Revealed by Arrhenius Plots

Electron and proton transfer reactions in enzymes are enigmatic and have attracted a great deal of theoretical, experimental, and practical attention. The oxidoreductases provide model systems for testing theoretical predictions, applying experimental techniques to gain insight into catalytic mechanisms, and creating industrially important bio(electro)conversion processes. Most previous and ongoing research on enzymatic electron transfer has exploited a theoretically and practically sound but limited approach that uses a series of structurally similar ("homologous") substrates, measures reaction rate constants and Gibbs free energies of reactions, and analyses trends predicted by electron transfer theory. This approach, proposed half a century ago, is based on a hitherto unproved hypothesis that pre-exponential factors of rate constants are similar for homologous substrates. Here, we propose a novel approach to investigating electron and proton transfer catalysed by oxidoreductases. We demonstrate the validity of this new approach for elucidating the kinetics of oxidation of "non-homologous" substrates catalysed by compound II of Coprinopsis cinerea and Armoracia rusticana peroxidases. This study - using the Marcus theory - demonstrates that reactions are not only limited by electron transfer, but a proton is transferred after the electron transfer event and thus both events control the reaction rate of peroxidase-catalysed oxidation of substrates.

An unprecedented intramolecular to intermolecular mechanistic switch in 1,1-diaminoazines leading to differential product formation during the I2-induced tandem oxidative transformation

Intermolecular cyclization of 1,1-diaminoazine.

Heme isomers substantially affect heme's electronic structure and function

Inspection of heme protein structures in the protein data bank reveals four isomers of heme characterized by different relative orientations of the vinyl side chains; remarkably, all these have been reported in multiple protein structures. Density functional theory computations explain this as due to similar energy of the isomers but with a sizable (25 kJ mol^-1) barrier to interconversion arising from restricted rotation around the conjugated bonds. The four isomers, EE, EZ, ZE, and ZZ, were then investigated as 4-coordinate hemes, as 5-coordinate deoxyhemes, in 6-coordinate O₂-adducts of globins and as compound I intermediates typical of heme peroxidases. Substantial differences were observed in electronic properties relevant to heme function: notably, the spin state energy gap of O₂-heme adducts, important for fast reversible binding of O₂, depends on the isomer state, and O₂-binding enthalpies change by up to 16 kJ mol^-1; redox potentials change by up to 0.2 V depending on the isomer, and the doublet-quartet energy splitting of compound I, central to "two-state" reactivity, is affected by up to ∼15 kJ mol^-1. These effects are consistently seen with three distinct density functionals, i.e. the effects are not method-dependent. Thus, the nature of the isomer state is an important but overlooked feature of heme chemistry and function, and previous and future studies of hemes may be reconsidered in this new context.

Deciphering the chemoselectivity of nickel-dependent quercetin 2,4-dioxygenase

QM/MM calculations were performed to elucidate the reaction mechanism and chemoselectivity of 2,4-QueD. The protonation state of the first-shell ligand Glu74 plays an important role in dictating the selectivity.

Kinetic Isotope Effect Probes the Reactive Spin State, As Well As the Geometric Feature and Constitution of the Transition State during H-Abstraction by Heme Compound II Complexes

What do experimentally measured kinetic isotope effects (KIEs) tell us about H-abstraction reactions with multispin-state reactivity options? Using DFT calculations with tunneling corrections for experimentally studied H-abstraction reactions of porphyrin-Compound II species (Chem.-Eur. J. 2014, 20, 14437; Angew. Chem., Int. Ed. 2008, 47, 7321) with cyclohexane, dihydroanthracene (DHA), and xanthene (Xan), we show here that KIE is a selective probe that identifies the experimentally reactive spin state. At the same time, comparison of calculated and experimental KIE values permits us to determine the structural orientation of the transition states, as well as the presence/absence of an axial ligand, and the effect of porphyrin substituents. The studied compound II (Cpd II) species involve porphine, and porphyrin ligands with different meso-substituents, TPFPP (tetrakis(pentafluorophenyl)porphyrin dianion) and TMP (tetramesitylporphyrin dianion), with and without imidazole axial ligands. The DFT calculations reveal three potential pathways: quintet and triplet σ-pathways (⁵Hσ and ³Hσ) that possess linear transition state (TS) structures, and a triplet π -pathway (³Hπ) having a bent TS structure. Without an axial ligand, the ⁵Hσ pathways for these Cpd II complexes cross below the triplet states. The axial ligand raises the barriers for the quintet and triplet σ-pathways and quenches any chances for two-state reactivity, thus proceeding via the ³Hπ pathway. All of these pathways exhibit characteristic KIE values: very large for ³Hπ (48-200), small for ⁵Hσ (3-9), and intermediate for ³Hσ (23-51). The calculated KIEs for (TPFPP)Fe^IV═O without an axial ligand reveal that ³Hσ is the only pathway having a KIE that matches the experimental values, for the reactions with DHA and Xan (Angew. Chem., Int. Ed. 2008, 47, 7321). Indeed, theory shows that tunneling significantly lowers the ³Hσ barrier rendering it the sole reactive state for the reaction. A prediction is made for the reactivity and KIE of (TMP)Fe^IV═O complex, and a comparison is made with the analogous nonheme complexes.

Challenging Density Functional Theory Calculations with Hemes and Porphyrins

In this paper we review recent advances in computational chemistry and specifically focus on the chemical description of heme proteins and synthetic porphyrins that act as both mimics of natural processes and technological uses. These are challenging biochemical systems involved in electron transfer as well as biocatalysis processes. In recent years computational tools have improved considerably and now can reproduce experimental spectroscopic and reactivity studies within a reasonable error margin (several kcal·mol(-1)). This paper gives recent examples from our groups, where we investigated heme and synthetic metal-porphyrin systems. The four case studies highlight how computational modelling can correctly reproduce experimental product distributions, predicted reactivity trends and guide interpretation of electronic structures of complex systems. The case studies focus on the calculations of a variety of spectroscopic features of porphyrins and show how computational modelling gives important insight that explains the experimental spectra and can lead to the design of porphyrins with tuned properties.

Ferryl protonation in oxoiron(IV) porphyrins and its role in oxygen transfer

Ferryl porphyrins, P-Fe(IV)═O, are central reactive intermediates in the catalytic cycles of numerous heme proteins and a variety of model systems. There has been considerable interest in elucidating factors, such as terminal oxo basicity, that may control ferryl reactivity. Here, the sulfonated, water-soluble ferryl porphyrin complexes tetramesitylporphyrin, oxoFe(IV)TMPS (FeTMPS-II), its 2,6-dichlorophenyl analogue, oxoFe(IV)TDClPS (FeTDClPS-II), and two other analogues are shown to be protonated under turnover conditions to produce the corresponding bis-aqua-iron(III) porphyrin cation radicals. The results reveal a novel internal electromeric equilibrium, P-Fe(IV)═O ⇆ P(+)-Fe(III)(OH2)2. Reversible pKa values in the range of 4-6.3 have been measured for this process by pH-jump, UV-vis spectroscopy. Ferryl protonation has important ramifications for C-H bond cleavage reactions mediated by oxoiron(IV) porphyrin cation radicals in protic media. Both solvent O-H and substrate C-H deuterium kinetic isotope effects are observed for these reactions, indicating that hydrocarbon oxidation by these oxoiron(IV) porphyrin cation radicals occurs via a solvent proton-coupled hydrogen atom transfer from the substrate that has not been previously described. The effective FeO-H bond dissociation energies for FeTMPS-II and FeTDClPS-II were estimated from similar kinetic reactivities of the corresponding oxoFe(IV)TMPS(+) and oxoFe(IV)TDClPS(+) species to be ∼92-94 kcal/mol. Similar values were calculated from the two-proton P(+)-Fe(III)(OH2)2 pKa(obs) and the porphyrin oxidation potentials, despite a 230 mV range for the iron porphyrins examined. Thus, the iron porphyrin with the lower ring oxidation potential has a compensating higher basicity of the ferryl oxygen. The solvent-derived proton adds significantly to the driving force for C-H bond scission.

Combined experimental and theoretical study on the reactivity of compounds I and II in horseradish peroxidase biomimetics

For the exploration of the intrinsic reactivity of two key active species in the catalytic cycle of horseradish peroxidase (HRP), Compound I (HRP-I) and Compound II (HRP-II), we generated in situ [Fe(IV) O(TMP(+.) )(2-MeIm)](+) and [Fe(IV) O(TMP)(2-MeIm)](0) (TMP=5,10,15,20-tetramesitylporphyrin; 2-MeIm=2-methylimidazole) as biomimetics for HRP-I and HRP-II, respectively. Their catalytic activities in epoxidation, hydrogen abstraction, and heteroatom oxidation reactions were studied in acetonitrile at -15 °C by utilizing rapid-scan UV/Vis spectroscopy. Comparison of the second-order rate constants measured for the direct reactions of the HRP-I and HRP-II mimics with the selected substrates clearly confirmed the outstanding oxidizing capability of the HRP-I mimic, which is significantly higher than that of HRP-II. The experimental study was supported by computational modeling (DFT calculations) of the oxidation mechanism of the selected substrates with the involvement of quartet and doublet HRP-I mimics ((2,4) Cpd I) and the closed-shell triplet spin HRP-II model ((3) Cpd II) as oxidizing species. The significantly lower activation barriers calculated for the oxidation systems involving (2,4) Cpd I than those found for (3) Cpd II are in line with the much higher oxidizing efficiency of the HRP-I mimic proven in the experimental part of the study. In addition, the DFT calculations show that all three reaction types catalyzed by HRP-I occur on the doublet spin surface in an effectively concerted manner, whereas these reactions may proceed in a stepwise mechanism with the HRP-II mimic as oxidant. However, the high desaturation or oxygen rebound barriers during CH bond activation processes by the HRP-II mimic predict a sufficient lifetime for the substrate radical formed through hydrogen abstraction. Thus, the theoretical calculations suggest that the dissociation of the substrate radical may be a more favorable pathway than desaturation or oxygen rebound processes. Importantly, depending on the electronic nature of the oxidizing species, that is, (2,4) Cpd I or (3) Cpd II, an interesting region-selective conversion phenomenon between sulfoxidation and H-atom abstraction was revealed in the course of the oxidation reaction of dimethylsulfide. The combined experimental and theoretical study on the elucidation of the intrinsic reactivity patterns of the HRP-I and HRP-II mimics provides a valuable tool for evaluating the particular role of the HRP active species in biological systems.

DFT studies on the mechanism of veratryl alcohol oxidation catalyzed by Cu–phen complexes

Density functional theory (DFT) calculations have been performed to investigate the catalytic mechanism for the oxidation of veratryl alcohol to veratraldehyde by Cu–phen catalyst. In this work, the most favorable mechanism has been predicted and we hope the obtained results could provide useful insights for the reaction process.

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.

Multireference and multiconfiguration ab initio methods in heme-related systems: what have we learned so far?

This work reviews the recent applications of ab initio multireference/multiconfiguration (MR/MC) electronic structure methods to heme-related systems, involving tetra-, penta-, and hexa-coordinate species, as well as the high-valent iron-oxo species. The current accuracy of these methods in the various systems is discussed, with special attention to potential sources of systematic errors. Thus, the review summarizes and tries to rationalize the key elements of MR/MC calculations, namely, the choice of the employed active space, especially the so-called double-shell effect that has already been recognized to be important in transition-metal-containing systems, and the impact of these elements on the spin-state energetics of heme species, as well as on the bonding mechanism of small molecules to the heme. It is shown that expansion of the MC wave function into one based on localized orbitals provides a compact and insightful view on some otherwise complex electronic structures. The effects of protein environment on the MR/MC results are summarized for the few available quantum mechanical/molecular mechanical (QM/MM) studies. Comparisons with corresponding DFT results are also made wherever available. Potential future directions are proposed.

Enzymatic degradation of A2E, a retinal pigment epithelial lipofuscin bisretinoid

Some forms of blinding macular disease are associated with excessive accumulation of bisretinoid lipofuscin in retinal pigment epithelial (RPE) cells of the eye. This material is refractory to lysosomal enzyme degradation. In addition to gene and drug-based therapies, treatments that reverse the accumulation of bisretinoid would be beneficial. Thus, we have examined the feasibility of degrading the bisretinoids by delivery of exogenous enzyme. As proof of principle we report that horseradish peroxidase (HRP) can cleave the RPE bisretinoid A2E. In both cell-free and cell-based assays, A2E levels were decreased in the presence of HRP. HRP-associated cleavage products were detected by ultraperformance liquid chromatography (UPLC) coupled to electrospray ionization mass spectrometry, and the structures of the aldehyde-bearing cleavage products were elucidated by 18O-labeling and 1H NMR spectroscopy and by recording UV−vis absorbance spectra. These findings indicate that RPE bisretinoids such as A2E can be degraded by appropriate enzyme activities.

Mechanistic insight into the alcohol oxidation mediated by an efficient green [CuBr(2)(2,2'-bipy)]-TEMPO catalyst by density functional method

Density functional theory (DFT) calculations have been performed to investigate the alcohol oxidation to acetaldehyde catalyzed by [CuBr(2)(2,2'-bipy)]-TEMPO (TEMPO stands for 2,2,6,6-tetramethylpiperidinyloxy; bipy stands for bipyridine). The total charge for the studied catalytic system is +1. The catalytic cycle consists of two parts, namely, alcohol oxidation and TEMPO regeneration. In alcohol oxidation, the reaction follows the Sheldon's mechanism for the proposed two mechanisms, i.e., Semmelhack's mechanism and Sheldon's mechanism. The water participation plays minor role in the H atom abstraction step. In TEMPO regeneration, the proposed three paths are competitive in energy. By comparing with experimental observation, it is found that the path, in which alcohol provides the proton to TEMPO(-) to produce TEMPOH followed by the oxidation of TEMPOH directly to TEMPO by O(2), is favored. In TEMPO regeneration, CH(3)CN acts as the ligand to stabilize the Cu(I) species during the catalytic cycle.

Theoretical insights into heme-catalyzed oxidation of cyclohexane to adipic acid

Adipic acid is a key compound in the chemical industry, where it is mainly used in the production of polymers. The conventional process of its generation requires vast amounts of energy and, moreover, produces environmentally deleterious substances. Thus, there is interest in alternative ways to gain adequate amounts of adipic acid. Experimental reports on a one-pot iron-catalyzed conversion of cyclohexane to adipic acid motivated a theoretical investigation based on density functional theory calculations. The process investigated is interesting because it requires less energy than contemporary methods and does not produce environmentally harmful side products. The aim of the present contribution is to gain insight into the mechanism of the iron-catalyzed cyclohexane conversion to provide a basis for the further development of this process. The rate-limiting step of the process is discussed, but considering the accuracy of the calculations, it is difficult to ensure whether the rate-limiting step is in the substrate oxidation or in the generation of the catalytically active species. It is shown that the slowest step in the substrate oxidation is the conversion of cyclohexanol to cyclohexane-1,2-diol. Hydrogen-atom transfer from one of the OH groups of cyclohexane-1,2-diol makes the intradiol cleavage occur spontaneously.

How do azoles inhibit cytochrome P450 enzymes? A density functional study

To examine how azole inhibitors interact with the heme active site of the cytochrome P450 enzymes, we have performed a series of density functional theory studies on azole binding. These are the first density functional studies on azole interactions with a heme center and give fundamental insight into how azoles inhibit the catalytic function of P450 enzymes. Since azoles come in many varieties, we tested three typical azole motifs representing a broad range of azole and azole-type inhibitors: methylimidazolate, methyltriazolate, and pyridine. These structural motifs represent typical azoles, such as econazole, fluconazole, and metyrapone. The calculations show that azole binding is a stepwise mechanism whereby first the water molecule from the resting state of P450 is released from the sixth binding site of the heme to create a pentacoordinated active site followed by coordination of the azole nitrogen to the heme iron. This process leads to the breaking of a hydrogen bond between the resting state water molecule and the approaching inhibitor molecule. Although, formally, the water molecule is released in the first step of the reaction mechanism and a pentacoordinated heme is created, this does not lead to an observed spin state crossing. Thus, we show that release of a water molecule from the resting state of P450 enzymes to create a pentacoordinated heme will lead to a doublet to quartet spin state crossing at an Fe-OH(2) distance of approximately 3.0 A, while the azole substitution process takes place at shorter distances. Azoles bind heme with significantly stronger binding energies than a water molecule, so that these inhibitors block the catalytic cycle of the enzyme and prevent oxygen binding and the catalysis of substrate oxidation. Perturbations within the active site (e.g., a polarized environment) have little effect on the relative energies of azole binding. Studies with an extra hydrogen-bonded ethanol molecule in the model, mimicking the active site of the CYP121 P450, show that the resting state and azole binding structures are close in energy, which may lead to chemical equilibrium between the two structures, as indeed observed with recent protein structural studies that have demonstrated two distinct azole binding mechanisms to P450 heme.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

The influence of X-rays on the structural studies of peroxide-derived myoglobin intermediates

In recent years, the awareness of potential radiation damage of metal centers in protein crystals during crystallographic data collection has received increasing attention. The radiation damage can lead to radiation-induced changes and reduction of the metal sites. One of the research fields where these concerns have been comprehensively addressed is the study of the reaction intermediates of the heme peroxidase and oxygenase reaction cycles. For both the resting states and the high-valent intermediates, the X-rays used in the structure determination have given undesired side effects through radiation-induced changes to the trapped intermediates. However, X-rays have been used to generate and trap the peroxy/hydroperoxy state in crystals. In this review, the structural work and the influence of X-rays on these intermediates in myoglobin are summarized and viewed in light of analogous studies on similar intermediates in peroxidases and oxygenases.

A rational basis for the axial ligand effect in C-H oxidation by [MnO(porphyrin)(X)]+ (X = H2O, OH-, O2-) from a DFT study

Oxyl radical character in the MnO group of the title system is shown from a density functional theory study to be essential for efficient C-H cleavage, which is a key step in C-H oxidation. Since oxyl species have elongated Mn-O bonds relative to the more usual oxo species of type MnO, the normal expectation would be that high trans-influence ligands X should facilitate oxyl character by elongating the Mn-O bond and thus enhance both oxyl character and reactivity. Contrary to this expectation, but in line with the experimental data (Jin, N.; Ibrahim, M.; Spiro, T. G.; Groves, J. T. J. Am. Chem. Soc. 2007, 129, 12416), we find that reactivity increases along the series X = O(2-) < OH(-) < H2O for the following reasons. The ground-state singlet (S) is unreactive for all X, and only the higher-energy triplet (T) and quintet (Q) states have the oxyl character needed for reactivity, but the higher trans-influence X ligands are also shown to increase the S/T and S/Q gaps, thus making attainment of the needed T and Q states harder. The latter effect is dominant, and high trans-influence X ligands thus disfavor reaction. The higher reactivity in the presence of acid noted by Groves and co-workers is thus rationalized by the preference for having X = H2O over OH(-) or O(2-).

A density functional theory investigation on the mechanism of the second half-reaction of nitric oxide synthase

Density functional theory methods have been employed to systematically investigate the overall mechanism of the second half-reaction of nitric oxide synthases. The initial heme-bound hydrogen peroxide intermediate previously identified is found to first undergo a simple rotation about its O-O peroxide bond. Then, via a "ping-pong" peroxidase-like mechanism the -O(in)H- proton is transferred back onto the substrate's -NO oxygen then subsequently onto the outer oxygen of the resulting Fe(heme)-OOH species. As a result, O(out) is released as H2O with concomitant formation of a compound I-type (Fe(heme)-O) species. Formation of the final citrulline and NO products can then be achieved in one step via a tetrahedral transition structure resulting from direct attack of the Fe(heme)-O moiety at the substrate's guanidinium carbon center. The possible role of alternative mechanisms involving a protonated compound II-type species or an initial transfer of only the -NH- hydrogen of the =NHOH+ group of N(omega)-hydroxy-L-arginine is also discussed.