Unique Structural and Electronic Features of Perferryl–Oxo Oxidant in Cytochrome P450

Low-Lying Electromeric States in Chloro-Ligated Iron(IV)-Oxo Porphyrin as a Model for Compound I, Studied with Second-Order Perturbation Theory Based on Density Matrix Renormalization Group

Employing second-order perturbation theory based on the density matrix renormalization group (DMRG-CASPT2), this work aims at providing a quantitative description of the spin state energetics of a chloro-ligated iron(IV)-oxo porphyrin as a model for the cytochromes P450 active species, also known as compound I (Cpd I). We explored DMRG-CASPT2 to its full extent with an extensive active space (up to 31 active orbitals) as well as a large number of renormalized states m (up to 10000). Different flavors of DMRG-CASPT2, using either the costly exact 4-particle reduced density matrix (4-RDM) or the cheaper cumulant approximated 4-RDM (cu(4)), were analyzed. All flavors essentially converge to similar relative energies between different spin states. Including a correction for the protein environment, we found a quartet Fe^IVO ground state and, more importantly, a thermally accessible doublet Fe^VO excited state that might directly contribute to the reactivity of this iron-oxo species. Our results also showed that cheaper approaches, such as CASPT2 based on a smaller active space or the cumulant approximation DMRG-cu(4)-CASPT2, are capable of accurately describing the spin state energetics of this species.

Revitalizing Spin Natural Orbital Analysis: Electronic Structures of Mixed-Valence Compounds, Singlet Biradicals, and Antiferromagnetically Coupled Systems

Chemical systems with open-shell electronic structure have been gaining attention these days. Their potential applications in first-row transition metal catalysis, molecular wires, photovoltaics and other potential applications have urged the adoption of a simple analysis tool to better understand their open-shell electronic structures, especially the role played by the unpaired electrons. Despite its lack of popularity, spin natural orbital (SNO) analysis is a tool we found to well-suit this purpose. We have therefore re-examined how the SNO could help us analyze some interesting open-shell systems, including mixed-valence compounds, singlet biradicals, and antiferromagnetically coupled systems. We found that some interesting patterns emerge from SNO analysis, especially those associated with exchange interaction. © 2019 Wiley Periodicals, Inc.

Conceptual Insights into DFT Spin-State Energetics of Octahedral Transition-Metal Complexes through a Density Difference Analysis

In this study, an intuitive concept is derived, which explains the characteristic dependence of spin-state energetics on the exact exchange admixture of DFT functionals in the case of octahedral transition metal complexes. The change in electron density distributions upon varying the admixture, c₃ , in the B3LYP functional is analyzed for archetype ionic and covalent systems as well as for the Fe²⁺ ion in an ideal octahedral field. An understanding of how the DFT description of the electronic structure of octahedral complexes changes as a function of c₃ is sought. A systematic spin-state energy analysis of 50 octahedral complexes of various metals and ligands with consistent experimental data is presented, allowing the derivation, in theory, of an optimal c₃ value for each system. The notion that the admixture dependence of DFT spin-state energetics stems from the treatment of nondynamic electrons arising from the mixing of (M-Lz2 )⁰ (dz2 )² and (M-Lx2-y2 )⁰ (dx2-y2 )² configurations into the dominant (M-Lx2-y2 )² (dx2-y2 )⁰ and (M-Lx2-y2 )² (dx2-y2 )⁰ ones in the low(er) spin states is put forward. That is, in the effort to mimic such electron-electron interactions, E_x^LDA overestimates, whereas exact exchange downplays the contribution of this type of electron correlation to the stability of low(er) spin states, leading to the widespread practical observation that the higher the exact exchange admixture, the more stable the high-spin-state configuration.

Role of the Proximal Cysteine Hydrogen Bonding Interaction in Cytochrome P450 2B4 Studied by Cryoreduction, Electron Paramagnetic Resonance, and Electron-Nuclear Double Resonance Spectroscopy

Crystallographic studies have shown that the F429H mutation of cytochrome P450 2B4 introduces an H-bond between His429 and the proximal thiolate ligand, Cys436, without altering the protein fold but sharply decreases the enzymatic activity and stabilizes the oxyferrous P450 2B4 complex. To characterize the influence of this hydrogen bond on the states of the catalytic cycle, we have used radiolytic cryoreduction combined with electron paramagnetic resonance (EPR) and (electron-nuclear double resonance (ENDOR) spectroscopy to study and compare their characteristics for wild-type (WT) P450 2B4 and the F429H mutant. (i) The addition of an H-bond to the axial Cys436 thiolate significantly changes the EPR signals of both low-spin and high-spin heme-iron(III) and the hyperfine couplings of the heme-pyrrole (14)N but has relatively little effect on the (1)H ENDOR spectra of the water ligand in the six-coordinate low-spin ferriheme state. These changes indicate that the H-bond introduced between His and the proximal cysteine decreases the extent of S → Fe electron donation and weakens the Fe(III)-S bond. (ii) The added H-bond changes the primary product of cryoreduction of the Fe(II) enzyme, which is trapped in the conformation of the parent Fe(II) state. In the wild-type enzyme, the added electron localizes on the porphyrin, generating an S = (3)/2 state with the anion radical exchange-coupled to the Fe(II). In the mutant, it localizes on the iron, generating an S = (1)/2 Fe(I) state. (iii) The additional H-bond has little effect on g values and (1)H-(14)N hyperfine couplings of the cryogenerated, ferric hydroperoxo intermediate but noticeably slows its decay during cryoannealing. (iv) In both the WT and the mutant enzyme, this decay shows a significant solvent kinetic isotope effect, indicating that the decay reflects a proton-assisted conversion to Compound I (Cpd I). (v) We confirm that Cpd I formed during the annealing of the cryogenerated hydroperoxy intermediate and that it is the active hydroxylating species in both WT P450 2B4 and the F429H mutant. (vi) Our data also indicate that the added H-bond of the mutation diminishes the reactivity of Cpd I.

F429 Regulation of Tunnels in Cytochrome P450 2B4: A Top Down Study of Multiple Molecular Dynamics Simulations

The root causes of the outcomes of the single-site mutation in enzymes remain by and large not well understood. This is the case of the F429H mutant of the cytochrome P450 (CYP) 2B4 enzyme where the substitution, on the proximal surface of the active site, of a conserved phenylalanine 429 residue with histidine seems to hamper the formation of the active species, Compound I (porphyrin cation radical-Fe^(IV) = O, Cpd I) from the ferric hydroperoxo (Fe^(III)OOH^-, Cpd 0) precursor. Here we report a study based on extensive molecular dynamic (MD) simulations of 4 CYP-2B4 point mutations compared to the WT enzyme, having the goal of better clarifying the importance of the proximal Phe429 residue on CYP 2B4 catalytic properties. To consolidate the huge amount of data coming from five simulations and extract the most distinct structural features of the five species studied we made an extensive use of cluster analysis. The results show that all studied single polymorphisms of F429, with different side chain properties: i) drastically alter the reservoir of conformations accessible by the protein, perturbing global dynamics ii) expose the thiolate group of residue Cys436 to the solvent, altering the electronic properties of Cpd0 and iii) affect the various ingress and egress channels connecting the distal sites with the bulk environment, altering the reversibility of these channels. In particular, it was observed that the wild type enzyme exhibits unique structural features as compared to all mutant species in terms of weak interactions (hydrogen bonds) that generate a completely different dynamical behavior of the complete system. Albeit not conclusive, the current computational investigation sheds some light on the subtle and critical effects that proximal single-site mutations can exert on the functional mechanisms of human microsomal CYPs which should go rather far beyond local structure characterization.

A comprehensive test set of epoxidation rate constants for iron(iv)-oxo porphyrin cation radical complexes

Trends in oxygen atom transfer to Compound I of the P450 models with an extensive test set have been studied and show a preferred regioselectivity of epoxidation over hydroxylation in the gas-phase for the first time.

Cytochrome P450 enzymes are heme based monoxygenases that catalyse a range of oxygen atom transfer reactions with various substrates, including aliphatic and aromatic hydroxylation as well as epoxidation reactions. The active species is short-lived and difficult to trap and characterize experimentally, moreover, it reacts in a regioselective manner with substrates leading to aliphatic hydroxylation and epoxidation products, but the origin of this regioselectivity is poorly understood. We have synthesized a model complex and studied it with low-pressure Fourier transform-ion cyclotron resonance (FT-ICR) mass spectrometry (MS). A novel approach was devised using the reaction of [Fe^III(TPFPP)]⁺ (TPFPP = meso-tetrakis(pentafluorophenyl)porphinato dianion) with iodosylbenzene as a terminal oxidant which leads to the production of ions corresponding to [Fe^IV(O)(TPFPP⁺˙)]⁺. This species was isolated in the gas-phase and studied in its reactivity with a variety of olefins. Product patterns and rate constants under Ideal Gas conditions were determined by FT-ICR MS. All substrates react with [Fe^IV(O)(TPFPP⁺˙)]⁺ by a more or less efficient oxygen atom transfer process. In addition, substrates with low ionization energies react by a charge-transfer channel, which enabled us to determine the electron affinity of [Fe^IV(O)(TPFPP⁺˙)]⁺ for the first time. Interestingly, no hydrogen atom abstraction pathways are observed for the reaction of [Fe^IV(O)(TPFPP⁺˙)]⁺ with prototypical olefins such as propene, cyclohexene and cyclohexadiene and also no kinetic isotope effect in the reaction rate is found, which suggests that the competition between epoxidation and hydroxylation – in the gas-phase – is in favour of substrate epoxidation. This notion further implies that P450 enzymes will need to adapt their substrate binding pocket, in order to enable favourable aliphatic hydroxylation over double bond epoxidation pathways. The MS studies yield a large test-set of experimental reaction rates of iron(iv)–oxo porphyrin cation radical complexes, so far unprecedented in the gas-phase, providing a benchmark for calibration studies using computational techniques. Preliminary computational results presented here confirm the observed trends excellently and rationalize the reactivities within the framework of thermochemical considerations and valence bond schemes.

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

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

Revisiting the role of exact exchange in DFT spin-state energetics of transition metal complexes

Sensitivity of DFT spin-state energetics to exact exchange is rooted in the description of metal–ligand bonding.

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).

Similarities of artificial photosystems by ruthenium oxo complexes and native water splitting systems

The nature of chemical bonds of ruthenium(Ru)-quinine(Q) complexes, mononuclear [Ru(trpy)(3,5-t-Bu(2)Q)(OH(2))](ClO(4))(2) (trpy = 2,2':6',2''-terpyridine, 3,5-di-tert-butyl-1,2-benzoquinone) (1), and binuclear [Ru(2)(btpyan)(3,6-di-Bu(2)Q)(2)(OH(2))](2+) (btpyan = 1,8-bis(2,2':6',2''-terpyrid-4'-yl)anthracene, 3,6-t-Bu(2)Q = 3,6-di-tert-butyl-1,2-benzoquinone) (2), has been investigated by broken-symmetry (BS) hybrid density functional (DFT) methods. BS DFT computations for the Ru complexes have elucidated that the closed-shell structure (2b) Ru(II)-Q complex is less stable than the open-shell structure (2bb) consisting of Ru(III) and semiquinone (SQ) radical fragments. These computations have also elucidated eight different electronic and spin structures of tetraradical intermediates that may be generated in the course of water splitting reaction. The Heisenberg spin Hamiltonian model for these species has been derived to elucidate six different effective exchange interactions (J) for four spin systems. Six J values have been determined using total energies of the eight (or seven) BS solutions for different spin configurations. The natural orbital analyses of these BS DFT solutions have also been performed in order to obtain natural orbitals and their occupation numbers, which are useful for the lucid understanding of the nature of chemical bonds of the Ru complexes. Implications of the computational results are discussed in relation to the proposed reaction mechanisms of water splitting reaction in artificial photosynthesis systems and the similarity between artificial and native water splitting systems.