Effect of Distal Interactions on O2 Binding to Heme

Heme-copper and Heme O2-derived synthetic (bioinorganic) chemistry toward an understanding of cytochrome c oxidase dioxygen chemistry

Cytochrome c oxidase (CcO), also widely known as mitochondrial electron-transport-chain complex IV, is a multi-subunit transmembrane protein responsible for catalyzing the last step of the electron transport chain, dioxygen reduction to water, which is essential to the establishment and maintenance of the membrane proton gradient that drives ATP synthesis. Although many intermediates in the CcO catalytic cycle have been spectroscopically and/or computationally authenticated, the specifics regarding the IP intermediate, hypothesized to be a heme-Cu (hydro)peroxo species whose O-O bond homolysis is supported by a hydrogen-bonding network of water molecules, are largely obscured by the fast kinetics of the A (FeIII-O2•-/CuI/Tyr) → PM (FeIV=O/CuII-OH/Tyr•) step. In this review, we have focused on the recent advancements in the design, development, and characterization of synthetic heme-peroxo‑copper model complexes, which can circumvent the abovementioned limitation, for the investigation of the formation of IP and its O-O cleavage chemistry. Novel findings regarding (a) proton and electron transfer (PT/ET) processes, together with their contributions to exogenous phenol induced O-O cleavage, (b) the stereo-electronic tunability of the secondary coordination sphere (especially hydrogen-bonding) on the geometric and spin state alteration of the heme-peroxo‑copper unit, and (c) a plausible mechanism for the Tyr-His cofactor biogenesis, are discussed in great detail. Additionally, since the ferric-superoxide and the ferryl-oxo (Compound II) species are critically involved in the CcO catalytic cycle, this review also highlights a few fundamental aspects of these heme-only (i.e., without copper) species, including the structural and reactivity influences of electron-donating trans-axial ligands and Lewis acid-promoted H-bonding.

How Nanoparticles Modify Adsorbed Proteins: Impact of Silica Nanoparticles on the Hemoglobin Active Site

The adsorption of proteins on surfaces has been studied for a long time, but the relationship between the structural and functional properties of the adsorbed protein and the adsorption mechanism remains unclear. Using hemoglobin adsorbed on silica nanoparticles, we have previously shown that hemoglobin's affinity towards oxygen increases with adsorption. Nevertheless, it was also shown that there were no significant changes in the quaternary and secondary structures. In order to understand the change in activity, we decided in this work to focus on the active sites of hemoglobin, the heme and its iron. After measuring adsorption isotherms of porcine hemoglobin on Ludox silica nanoparticles, we analyzed the structural modifications of adsorbed hemoglobin by X-ray absorption spectroscopy and circular dichroism spectra in the Soret region. It was found that upon adsorption, there were modifications in the heme pocket environment due to changes in the angles of the heme vinyl functions. These alterations can explain the greater affinity observed.

Insights into the effect of distal histidine and water hydrogen bonding on NO ligation to ferrous and ferric heme: a DFT study

The effect of distal histidine on ligation of NO to ferrous and ferric-heme, has been investigated with the high-level density functional theoretical (DFT) method. It has been predicted that the distal histidine significantly stabilizes the interaction of NO ferrous-heme (by −2.70 kcal mol⁻¹). Also, water hydrogen bonding is quite effective in strengthening the Fe–NO bond in ferrous heme. In contrast in ferric heme, due to the large distance between the H₂O and O(NO) and lack of hydrogen bonding, the distal histidine exhibits only a slight effect on the binding of NO to the ferric analogue. Concerning the bond nature of Fe^II–NO and Fe^III–NO in heme, a QTAIM analysis predicts a partially covalent and ionic bond nature in both systems.

The effect of distal histidine on ligation of NO to ferrous and ferric-heme, has been investigated with the high-level density functional theoretical (DFT) method.

Application of Lever's EL Parameter Scale toward Fe(II)/Fe(III) versus Pc(2-)/Pc(1-) Oxidation Process Crossover Point in Axially Coordinated Iron(II) Phthalocyanine Complexes

The electronic structures and, particularly, the nature of the HOMO in a series of PcFeL₂, PcFeL'L″, and [PcFeX₂]^2- complexes (Pc = phthalocyaninato(2-) ligand; L = NH₃, n-BuNH₂, imidazole (Im), pyridine (Py), PMe₃, PBu₃, t-BuNC, P(OBu)₃, and DMSO; L' = CO; L″ = NH₃ or n-BuNH₂; X = NCO^-, NCS^-, CN^-, imidazolate (Im^-), or 1,2,4-triazolate(Tz^-)) were probed by electrochemical, spectroelectrochemical, and chemical oxidation as well as theoretical (density functional theory, DFT) studies. In general, energies of the metal-centered occupied orbitals in various six-coordinate iron phthalocyanine complexes correlate well with Lever Electrochemical Parameter E_L and intercross the phthalocyanine-centered a_1u orbital in several compounds with moderate-to-strong π-accepting axial ligands. In these cases, an oxidation of the phthalocyanine macrocycle (Pc(2-)/Pc(1-)) rather than the central metal ion (Fe(II)/Fe(III)) was theoretically predicted and experimentally confirmed.

Accurate Prediction of the Excited States in the Fully Conjugated Porphyrin Tapes across the Full Spectral Range: A Story of the Interplay between π-π* and Intramolecular Charge-Transfer Transitions in Soft Chromophores

The ability of density functional theory (DFT) and time-dependent DFT (TDDFT) methods for the accurate prediction of the energies and oscillator strengths of the excited states in a series of fully conjugated meso-meso β-β β-β triple-linked porphyrin oligomers (porphyrin tapes 2-12) was probed in the gas phase and solution using several exchange-correlation functionals. It was demonstrated that the use of the hybrid B3LYP functional provides a good compromise for the accurate prediction of the localized π-π* and intramolecular charge-transfer transitions, thus allowing confident interpretation of the UV-vis-NIR spectra of porphyrin oligomers. The TDDFT-based sum-over-state (SOS) calculations for the porphyrin tape dimer 2 and trimer 3 as well as parent monomer 1 correctly predicted the signs and shapes of the magnetic circular dichroism (MCD) signals in the low-energy region of the spectra.

Accurate Prediction of Mössbauer Hyperfine Parameters in Bis-Axially Coordinated Iron(II) Phthalocyanines Using Density Functional Theory Calculations: A Story of a Single Orbital Revealed by Natural Bond Orbital Analysis

Density Functional Theory (DFT) calculations coupled with several exchange-correlation functionals were used for the prediction of Mössbauer hyperfine parameters of 36 bis-axially coordinated iron(II) phthalocyanine complexes with the general formulas PcFeL₂, PcFeL'L″, and [PcFeX₂]^2-, including four new compounds. Both gas-phase and PCM calculations using BPW91 and MN12L exchange-correlation functionals were found to accurately predict both Mössbauer quadrupole splittings and the correct trends in experimentally observed isomer shifts. In comparison, hybrid exchange-correlation functionals underestimated quadrupole splittings, while still accurately predicted isomer shifts. Out of ∼40 exchange-correlation functionals tested, only MN12L was found to correctly reproduce quadrupole splitting trends in the PcFeL₂ complexes coordinated with phosphorus-donor axial ligands (i.e., P(OnBu)₃ ≈ P(OEt)₃ < PMe₃ < P[(CH₂O)₂CH₂]-p-C₆H₄NO₂ < PEt₃ ≈ PnBu₃). Natural Bond Orbital (NBO) analysis was successfully used to explain the general trends in the observed quadrupole splitting for all compounds of interest. In particular, the general trends in the quadrupole splitting correlate well with the axial ligand dependent, NBO-predicted population of the 3d_z2 orbital of the Fe ion and are reflective of the hypothesis proposed by Ohya and co-workers ( Inorg. Chem., 1984, 23, 1303) on the adaptability of the phthalocyanine's π-system toward Fe-L_ax interactions. The first X-ray crystal structure of a PcFeL₂ complex with axial phosphine ligands is also reported.

Single gold nanoparticle-driven heme cofactor nanozyme as an unprecedented oxidase mimetic

The catalytic diversity of heme enzymes is a perpetuating pursuit for biomimetic chemistry, but heme nanozymes exhibit catalytic activity only reminiscent of peroxidases. Miraculously, the oxidase-like catalytic function of a heme cofactor is elicited with the help of gold nanoparticles (AuNPs) by maintaining heme with a low-valence state (ferrous) in a confined configuration.

Dioxygen Binding to all 3d, 4d, and 5d Transition Metals from Coupled-Cluster Theory

Understanding how transition metals bind and activate dioxygen (O₂ ) is limited by experimental and theoretical uncertainties, making accurate quantum mechanical descriptors of interest. Here we report coupled-cluster CCSD(T) energies with large basis sets and vibrational and relativistic corrections for 160 3d, 4d, and 5d metal-O₂ systems. We define four reaction energies (120 in total for the 30 metals) that quantify O-O activation and reveal linear relationships between metal-oxygen and O-O binding energies. The CCSD(T) data can be combined with thermochemical cycles to estimate chemisorption and physisorption energies for each metal from metal oxide embedding energies, in good correlation with atomization enthalpies (R² =0.75). Spin-geometry variations can break the linearities, of interest to circumventing the Sabatier principle. Pt, Pd, Co, and Fe form a distinct group with the weakest O₂ binding. R² up to 0.84 between surface adsorption energies and our energies for MO₂ systems indicate relevance also to real catalytic systems.

Energy decomposition analysis based on broken symmetry unrestricted density functional theory

In this paper, the generalized Kohn-Sham energy decomposition analysis (GKS-EDA) scheme is extended to molecular interactions in open shell singlet states, which is a challenge for many popular EDA methods due to the multireference character. Based on broken symmetry (BS) unrestricted density functional theory with a spin projection approximation, the extension scheme, named GKS-EDA(BS) in this paper, divides the total interaction energy into electrostatic, exchange-repulsion, polarization, correlation, and dispersion terms. Test examples include the pancake bond in the phenalenyl dimer, the ligand interactions in the Fe(ii)-porphyrin complexes, and the radical interactions in dehydrogenated guanine-cytosine base pairs and show that GKS-EDA(BS) is a practical EDA tool for open shell singlet systems.

DFT Fea3-O/O-O Vibrational Frequency Calculations over Catalytic Reaction Cycle States in the Dinuclear Center of Cytochrome c Oxidase

Density functional vibrational frequency calculations have been performed on eight geometry optimized cytochrome c oxidase (CcO) dinuclear center (DNC) reaction cycle intermediates and on the oxymyoglobin (oxyMb) active site. The calculated Fe-O and O-O stretching modes and their frequency shifts along the reaction cycle have been compared with the available resonance Raman (rR) measurements. The calculations support the proposal that in state A[Fea33+-O2-•···CuB+] of CcO, O2 binds with Fea32+ in a similar bent end-on geometry to that in oxyMb. The calculations show that the observed 20 cm-1 shift of the Fea3-O stretching mode from the PR to F state is caused by the protonation of the OH- ligand on CuB2+ (PR[Fea34+═O2-···HO--CuB2+] → F[Fea34+═O2-···H2O-CuB2+]), and that the H2O ligand is still on the CuB2+ site in the rR identified F[Fea34+═O2-···H2O-CuB2+] state. Further, the observed rR band at 356 cm-1 between states PR and F is likely an O-Fea3-porphyrin bending mode. The observed 450 cm-1 low Fea3-O frequency mode for the OH active oxidized state has been reproduced by our calculations on a nearly symmetrically bridged Fea33+-OH-CuB2+ structure with a relatively long Fea3-O distance near 2 Å. Based on Badger's rule, the calculated Fea3-O distances correlate well with the calculated νFe-O-2/3 (νFe-O is the Fea3-O stretching frequency) with correlation coefficient R = 0.973.

Effect of toxic ligands on O2 binding to heme and their toxicity mechanism

Heme, as the cofactor and active site of Hb, enables Hb to carry out the necessary function required for O2 management for life, that is, reversible O2 binding for transport. In this paper, the microscopic mechanism of heme-associated poisoning has been elucidated from the perspective of electronic interaction by performing first-principles calculations. The results show that the functional groups (-CHO, -COOH, -NO2, -NH2) and CN exhibit a stronger affinity for heme than O2 and are more likely to occupy the O2 binding site, which results in the loss of the ability of heme to carry O2. Moreover, the addition of functional groups, CO and CN to heme at the side site can cause a pronounced enhancement toward the O2 binding characteristics of heme, which prevents heme from releasing O2 to oxygen-consuming tissues as the blood circulates. The reversible O2 binding function of heme is disrupted by the presence of these toxic ligands in the heme binding pocket, which greatly affects O2 transport in the blood. The inability of tissues to obtain O2 leads to tissue hypoxia, which is the main cause of poisoning. Based on the energy, geometry and electronic properties, the hypoxia mechanism proposed by us coincides well with experiment, and the research has the potential to provide a theoretical reference for the relevant areas of bioscience.

Spin Interconversion of Heme-Peroxo-Copper Complexes Facilitated by Intramolecular Hydrogen-Bonding Interactions

Synthetic peroxo-bridged high-spin (HS) heme-(μ-η2:η1-O22-)-Cu(L) complexes incorporating (as part of the copper ligand) intramolecular hydrogen-bond (H-bond) capabilities and/or steric effects are herein demonstrated to affect the complex's electronic and geometric structure, notably impacting the spin state. An H-bonding interaction with the peroxo core favors a low-spin (LS) heme-(μ-η1:η1-O22-)-Cu(L) structure, resulting in a reversible temperature-dependent interconversion of spin state (5 coordinate HS to 6 coordinate LS). The LS state dominates at low temperatures, even in the absence of a strong trans-axial heme ligand. Lewis base addition inhibits the H-bond facilitated spin interconversion by competition for the H-bond donor, illustrating the precise H-bonding interaction required to induce spin-crossover (SCO). Resonance Raman spectroscopy (rR) shows that the H-bonding pendant interacts with the bridging peroxide ligand to stabilize the LS but not the HS state. The H-bond (to the Cu-bound O atom) acts to weaken the O-O bond and strengthen the Fe-O bond, exhibiting ν(M-O) and ν(O-O) values comparable to analogous known LS complexes with a strong donating trans-axial ligand, 1,5-dicyclohexylimidazole, (DCHIm)heme-(μ-η1:η1-O22-)-Cu(L). Variable-temperature (-90 to -130 °C) UV-vis and 2H NMR spectroscopies confirm the SCO process and implicate the involvement of solvent binding. Examining a case of solvent binding without SCO, thermodynamic parameters were obtained from a van't Hoff analysis, accounting for its contribution in SCO. Taken together, these data provide evidence for the H-bond group facilitating a core geometry change and allowing solvent to bind, stabilizing a LS state. The rR data, complemented by DFT analysis, reveal a stronger H-bonding interaction with the peroxo core in the LS compared to the HS complexes, which enthalpically favors the LS state. These insights enhance our fundamental understanding of secondary coordination sphere influences in metalloenzymes.

Resonant inelastic X-ray scattering determination of the electronic structure of oxyhemoglobin and its model complex

Hemoglobin and myoglobin are oxygen-binding proteins with S = 0 heme {FeO₂}⁸ active sites. The electronic structure of these sites has been the subject of much debate. This study utilizes Fe K-edge X-ray absorption spectroscopy (XAS) and 1s2p resonant inelastic X-ray scattering (RIXS) to study oxyhemoglobin and a related heme {FeO₂}⁸ model compound, [(pfp)Fe(1-MeIm)(O₂)] (pfp = meso-tetra(α,α,α,α-o-pivalamido-phenyl)porphyrin, or TpivPP, 1-MeIm = 1-methylimidazole) (pfpO₂), which was previously analyzed using L-edge XAS. The K-edge XAS and RIXS data of pfpO₂ and oxyhemoglobin are compared with the data for low-spin Fe^II and Fe^III [Fe(tpp)(Im)₂]^0/+ (tpp = tetra-phenyl porphyrin) compounds, which serve as heme references. The X-ray data show that pfpO₂ is similar to Fe^II, while oxyhemoglobin is qualitatively similar to Fe^III, but with significant quantitative differences. Density-functional theory (DFT) calculations show that the difference between pfpO₂ and oxyhemoglobin is due to a distal histidine H bond to O₂ and the less hydrophobic environment in the protein, which lead to more backbonding into the O₂ A valence bond configuration interaction multiplet model is used to analyze the RIXS data and show that pfpO₂ is dominantly Fe^II with 6-8% Fe^III character, while oxyhemoglobin has a very mixed wave function that has 50-77% Fe^III character and a partially polarized Fe-O₂ π-bond.

Synthetic Fe/Cu Complexes: Toward Understanding Heme-Copper Oxidase Structure and Function

Heme-copper oxidases (HCOs) are terminal enzymes on the mitochondrial or bacterial respiratory electron transport chain, which utilize a unique heterobinuclear active site to catalyze the 4H+/4e- reduction of dioxygen to water. This process involves a proton-coupled electron transfer (PCET) from a tyrosine (phenolic) residue and additional redox events coupled to transmembrane proton pumping and ATP synthesis. Given that HCOs are large, complex, membrane-bound enzymes, bioinspired synthetic model chemistry is a promising approach to better understand heme-Cu-mediated dioxygen reduction, including the details of proton and electron movements. This review encompasses important aspects of heme-O2 and copper-O2 (bio)chemistries as they relate to the design and interpretation of small molecule model systems and provides perspectives from fundamental coordination chemistry, which can be applied to the understanding of HCO activity. We focus on recent advancements from studies of heme-Cu models, evaluating experimental and computational results, which highlight important fundamental structure-function relationships. Finally, we provide an outlook for future potential contributions from synthetic inorganic chemistry and discuss their implications with relevance to biological O2-reduction.

Chemical Bond Energies of 3d Transition Metals Studied by Density Functional Theory

Despite their vast importance to inorganic chemistry, materials science, and catalysis, the accuracy of modeling the formation or cleavage of metal-ligand (M-L) bonds depends greatly on the chosen functional and the type of bond in a way that is not systematically understood. In order to approach a state of high-accuracy DFT for rational prediction of chemistry and catalysis, such system-dependencies need to be resolved. We studied 30 different density functionals applied to a "balanced data set" of 60 experimental diatomic M-L bond energies; this data set has no bias toward any d^q configuration, metal, bond type, or ligand as all of these occur to the same extent, and we can therefore identify accuracy bottlenecks. We show that the performance of a functional is very dependent on data set choice, and we dissect these effects into system type. In addition to the use of balanced data sets, we also argue that the precision (rather than just accuracy) of a functional is of interest, measured by standard deviations of the errors. There are distinct system dependencies both in the ligand and metal series: Hydrides are best described by a very large HF exchange percentage, possibly due to self-interaction error, whereas halides are best described by very small (0-10%) HF exchange fractions, and double-bond enforcing oxides and sulfides favor 10-25% HF exchange, as is also average for the full data set. Thus, average HF requirements hide major system-dependent requirements. For late transition metals Co-Zn, HF percentage of 0-10% is favored, whereas for the early transition metals Sc-Fe hybrid functionals with 20% HF exchange or higher are commonly favored. Accordingly, B3LYP is an excellent choice for early d-block but a poor choice for late transition metals. We conclude that DFT intrinsically underestimates the bond strengths of late vs early transition metals, correlating with increased effective nuclear charge. Thus, the revised RPBE, which reduces the overbinding tendency of PBE, is mainly an advantage for the early and mid transition metals and not very much for the late transition metals, i.e. there is a metal-dependent effect of the relative performance of RPBE vs PBE, which are widely used to study adsorption energetics on metal surfaces. Overall, the best performing functionals are PW6B95, the MN15 and MN15-L functionals, and the double hybrid B2PLYP.

Redox Potentials and Electronic States of Iron Porphyrin IX Adsorbed on Single Crystal Gold Electrode Surfaces

Metalloporphyrins are active sites in metalloproteins and synthetic catalysts. They have also been studied extensively by electrochemistry as well as being prominent targets in electrochemical scanning tunneling microscopy (STM). Previous studies of FePPIX adsorbed on graphite and alkylthiol modified Au electrodes showed a pair of reversible Fe(III/II)PPIX peaks at about -0.41 V (vs NHE) at high solution pH. We recently used iron protoporphyrin IX (FePPIX) as an intercalating probe for long-range electrochemical electron transfer through a G-quadruplex oligonucleotide (DNAzyme); this study disclosed two, rather than a single pair of voltammetric peaks with a new and dominating peak, shifted 200 mV positive relative to the ≈-0.4 V peak. Prompted by this unexpected observation, we report here a study of the voltammetry of FePPIX itself on single-crystal Au(111), (100), and (110) and polycrystalline Au electrode surfaces. In all cases the dominating pair of new Fe(III/II)PPIX redox peaks, shifted positively by more than 200 mV compared to those of previous studies appeared. This observation is supported by density functional theory (DFT) which shows that strong dispersion forces in the FePPIX/Au electronic interaction drive the midpoint potential toward positive values. The FePPIX spin states depend on interaction with the Au(111) interface, converting all the Fe(II)/(III)PPIX species into low-spin states. These results support electrochemical evidence for the nature of the electronic coupling between FePPIX and Au-surfaces, and the electronic states of adsorbate molecules, with a bearing also on recent reports of magnetic FePPIX/Au(111) interactions in ultrahigh vacuum (UHV).

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.

Ligand-Directed Reactivity in Dioxygen and Water Binding to cis-[Pd(NHC)2(η2-O2)]

Reaction of [Pd(IPr)₂] (IPr = 1,3-bis(2,6-diisopropylphenyl)imidazol-2-ylidene) and O₂ leads to the surprising discovery that at low temperature the initial reaction product is a highly labile peroxide complex cis-[Pd(IPr)₂(η²-O₂)]. At temperatures ≳ -40 °C, cis-[Pd(IPr)₂(η²-O₂)] adds a second O₂ to form trans-[Pd(IPr)₂(η¹-O₂)₂]. Squid magnetometry and EPR studies yield data that are consistent with a singlet diradical ground state with a thermally accessible triplet state for this unique bis-superoxide complex. In addition to reaction with O₂, cis-[Pd(IPr)₂(η²-O₂)] reacts at low temperature with H₂O in methanol/ether solution to form trans-[Pd(IPr)₂(OH)(OOH)]. The crystal structure of trans-[Pd(IPr)₂(OOH)(OH)] is reported. Neither reaction with O₂ nor reaction with H₂O occurs under comparable conditions for cis-[Pd(IMes)₂(η²-O₂)] (IMes = 1,3-bis(2,4,6-trimethylphenyl)imidazol-2-ylidene). The increased reactivity of cis-[Pd(IPr)₂(η²-O₂)] is attributed to the enthalpy of binding of O₂ to [Pd(IPr)₂] (-14.5 ± 1.0 kcal/mol) that is approximately one-half that of [Pd(IMes)₂] (-27.9 ± 1.5 kcal/mol). Computational studies identify the cause as interligand repulsion forcing a wider C-Pd-C angle and tilting of the NHC plane in cis-[Pd(IPr)₂(η²-O₂)]. Arene-arene interactions are more favorable and serve to further stabilize cis-[Pd(IMes)₂(η²-O₂)]. Inclusion of dispersion effects in DFT calculations leads to improved agreement between experimental and computational enthalpies of O₂ binding. A complete reaction diagram is constructed for formation of trans-[Pd(IPr)₂(η¹-O₂)₂] and leads to the conclusion that kinetic factors inhibit formation of trans-[Pd(IMes)₂(η¹-O₂)₂] at the low temperatures at which it is thermodynamically favored. Failure to detect the predicted T-shaped intermediate trans-[Pd(NHC)₂(η¹-O₂)] for either NHC = IMes or IPr is attributed to dynamic effects. A partial potential energy diagram for initial binding of O₂ is constructed. A range of low-energy pathways at different angles of approach are present and blur the distinction between pure "side-on" or "end-on" trajectories for oxygen binding.

The dioxygen adducts of iron and manganese porphyrins: electronic structure and binding energy

The electronic structures of adducts of O₂ and metal porphyrins were thoroughly investigated by highly accurate DMRG-CASPT2.

The secondary coordination sphere and axial ligand effects on oxygen reduction reaction by iron porphyrins: a DFT computational study

Oxygen reduction reaction (ORR) catalyzed by a bio-inspired iron porphyrin bearing a hanging carboxylic acid group over the porphyrin ring, and a tethered axial imidazole ligand was studied by DFT calculations. BP86 free energy calculations of the redox potentials and pK a's of reaction components involved in the proton coupled electron transfer (PCET) reactions of the ferric-hydroxo and -superoxo complexes were performed based on Born-Haber thermodynamic cycle in conjunction with a continuum solvation model. The comparison was made with iron porphyrins that lack either in the hanging acid group or axial ligand, suggesting that H-bond interaction between the carboxylic acid and iron-bound hydroxo, aquo, superoxo, and peroxo ligands (de)stabilizes the Fe-O bonding, resulting in the increase in the reduction potential of the ferric complexes. The axial ligand interaction with the imidazole raises the affinity of the iron-bound superoxo and peroxo ligands for proton. In addition, a low-spin end-on ferric-hydroperoxo intermediate, a key precursor for O-O cleavage, can be stabilized in the presence of axial ligation. Thus, selective and efficient ORR of iron porphyrin can be achieved with the aid of the secondary coordination sphere and axial ligand interactions.

The secondary coordination sphere controlled reactivity of a ferric-superoxo heme: unexpected conversion to a ferric hydroperoxo intermediate by reaction with a high-spin ferrous heme

A bio-inspired heme complex involving both a proton donor and an axial imidazole ligand reduces the activation energy for the formation of a ferric hydroperoxo intermediate. A high-spin ferrous heme is shown to be capable of reducing its superoxy species to generate a ferric hydroperoxo intermediate for the first time.

Theoretical Study of Spin Crossover in 30 Iron Complexes

Iron complexes are important spin crossover (SCO) systems with vital roles in oxidative metabolism and promising technological potential. The SCO tendency depends on the free energy balance of high- and low-spin states, which again depends on physical effects such as dispersion, relativistic effects, and vibrational entropy. This work studied 30 different iron SCO systems with experimentally known thermochemical data, using 12 different density functionals. Remarkably general entropy-enthalpy compensation across SCO systems was identified (R = 0.82, p = 0.002) that should be considered in rational SCO design. Iron(II) complexes displayed higher ΔH and ΔS values than iron(III) complexes and also less steep compensation effects. First-coordination sphere ΔS values computed from numerical frequencies reproduce most of the experimental entropy and should thus be included when modeling spin-state changes in inorganic chemistry (R = 0.52, p = 3.4 × 10(-3); standard error in TΔS ≈ 4.4 kJ/mol at 298 K vs 16 kJ/mol of total TΔS on average). Zero-point energies favored high-spin states by 9 kJ/mol on average. Interestingly, dispersion effects are surprisingly large for the SCO process (average: 9 kJ/mol, but up to 33 kJ/mol) and favor the more compact low-spin state. Relativistic effects favor low-spin by ∼9 kJ/mol on average, but up to 24 kJ/mol. B3LYP*, TPSSh, B2PLYP, and PW6B95 performed best for the typical calculation scheme that includes ZPE. However, if relativistic and dispersion effects are included, only B3LYP* remained accurate. On average, high-spin was favored by LYP by 11-15 kJ/mol relative to other correlation functionals, and by 4.2 kJ/mol per 1% HF exchange in hybrids. 13% HF exchange was optimal without dispersion, and 15% was optimal with all effects included for these systems.

Highly efficient As(v)/Sb(v) removal by magnetic sludge composite: synthesis, characterization, equilibrium, and mechanism studies

Magnetic sludge with the unique flexibility, low thermal conductance and plentiful adsorption sites was fabricated for highly efficient As/Sb removal.

Role of Spin States in Nitric Oxide Binding to Cobalt(II) and Manganese(II) Porphyrins. Is Tighter Binding Always Stronger?

Binding of nitric oxide (NO) to metalloporphyrins and heme groups is important in biochemistry while challenging to describe accurately by density functional theory (DFT) calculations. Here, the structural and thermochemical aspect of NO binding to Co(II) and Mn(II) porphyrins is investigated by DFT and DFT-D (dispersion-corrected) calculations, supported by reliable coupled-cluster methodology (CCSD(T)), and critically correlated with the experimental data. It is argued that whereas the bonding of NO to Co(II) porphyrin is a simple radical recombination, the bonding of NO to Mn(II) porphyrin is accompanied by a crossing of spin states. For this reason, the spin-state conversion energy contributes to the Mn-NO bond energy, and the paradigmatic correlation between bond length and bond energy is violated for the considered nitrosyl complexes: the Mn-NO bond is (structurally) shorter by ∼0.2 Å, albeit (energetically) weaker by a few kcal/mol, compared with the Co-NO bond. Moreover, none of the many tested DFT methods can reproduce the Co-NO and Mn-NO bond energies simultaneously, except for calculations with B3LYP*-D3, TPSSh-D3, and M06-D3 methods supplemented with the proposed spin-state energy correction (to compensate for an error on the calculated spin-state conversion energy). The results of this study are important to appreciate the role of spin-state changes in ligand binding properties of heme-related models. They also highlight the need for accurate calculations for correct interpretation of experimental data, including the qualitative structure-energy relationship.

Functionalization of nanostructured cerium oxide films with histidine

The surfaces of polycrystalline cerium oxide films were modified by histidine adsorption under vacuum and characterized by the synchrotron based techniques of core and valence level photoemission, resonant photoemission and near edge X-ray absorption spectroscopy, as well as atomic force microscopy. Histidine is strongly bound to the oxide surface in the anionic form through the deprotonated carboxylate group, and forms a disordered molecular adlayer. The imidazole ring and the amino side group do not form bonds with the substrate but are involved in the intermolecular hydrogen bonding which stabilizes the molecular adlayer. The surface reaction with histidine results in water desorption accompanied by oxide reduction, which is propagated into the bulk of the film. Previously studied, well-characterized surfaces are a guide to the chemistry of the present polycrystalline surface and histidine bonds via the carboxylate group in both cases. In contrast, bonding via the imidazole group occurs on the well-ordered surface but not in the present case. The morphology and structure of the cerium oxide are decisive factors which define the adsorption geometry of the histidine adlayer.

Halide binding and inhibition of laccase copper clusters: the role of reorganization energy

Laccase-like proteins are multicopper oxidases involved in several biological and industrial processes. Their application is commonly limited due to inhibition by fluoride and chloride, and as-isolated proteins are often substantially activated by heat, suggesting that multiple redox states can complicate characterization. Understanding these processes at the molecular level is thus desirable but theoretically unexplored. This paper reports systematic calculations of geometries, reorganization energies, and ionization energies for all partly oxidized states of the trinuclear copper clusters in realistic models with ∼200 atoms. Corrections for scalar-relativistic effects, dispersion, and thermal effects were estimated. Fluoride, chloride, hydroxide, or water was bound to the T2 copper site of the oxidized resting state, and the peroxo intermediate was also computed for reference. Antiferromagnetic coupling, assigned oxidation states, and general structures were consistent with known spectroscopic data. The computations show that (i) ligands bound to the T2 site substantially increase the reorganization energy of the second reduction of the resting state and reduce the redox potentials, providing a possible mechanism for inhibition; (ii) the reorganization energy is particularly large for F(-) but also high for Cl(-), consistent with the experimental tendency of inhibition; (iii) reduction leads to release of Cl(-) from the T2 site, suggesting a mechanism for heat/reduction activation of laccases by dissociation of inhibiting halides or hydroxide from T2.

Co-C dissociation of adenosylcobalamin (coenzyme B12): role of dispersion, induction effects, solvent polarity, and relativistic and thermal corrections

Quantum-chemical cluster modeling is challenged in the limit of large, soft systems by the effects of dispersion and solvent, and well as other physical interactions. Adenosylcobalamin (AdoCbl, coenzyme B12), as one of the most complex cofactors in life, constitutes such a challenge. The cleavage of its unique organometallic Co-C bond has inspired multiple studies of this cofactor. This paper reports the fully relaxed potential energy surface of Co-C cleavage of AdoCbl, including for the first time all side-chain interactions with the dissociating Ado group. Various methods and corrections for dispersion, relativistic effects, solvent polarity, basis set superposition error, and thermal and vibrational effects were investigated, totaling more than 550 single-point energies for the large model. The results show immense variability depending on method, including solvation, functional type, and dispersion, challenging the conceived accuracy of methods used for such systems. In particular, B3LYP-D3 seems to severely underestimate the Co-C bond strength, consistent with previous results, and BP86 remains accurate for cobalamins when dispersion interactions are accounted for.