Linking chemical electron-proton transfer to proton pumping in cytochrome c oxidase: broken-symmetry DFT exploration of intermediates along the catalytic reaction pathway of the iron-copper dinuclear complex

A Bioinspired Nonheme FeIII-(O22-)-CuII Complex with an St = 1 Ground State

Cytochrome c oxidase (CcO) is a heme copper oxidase (HCO) that catalyzes the natural reduction of oxygen to water. A profound understanding of some of the elementary steps leading to the intricate 4e-/4H+ reduction of O2 is presently lacking. A total spin St = 1 FeIII-(O22-)-CuII (IP) intermediate is proposed to reduce the overpotentials associated with the reductive O-O bond rupture by allowing electron transfer from a tyrosine moiety without the necessity of any spin-surface crossing. Direct evidence of the involvement of IP in the CcO catalytic cycle is, however, missing. A number of heme copper peroxido complexes have been prepared as synthetic models of IP, but all of them possess the catalytically nonrelevant St = 0 ground state resulting from antiferromagnetic coupling between the S = 1/2 FeIII and CuII centers. In a complete nonheme approach, we now report the spectroscopic characterization and reactivity of the FeIII-(O22-)-CuII intermediates 1 and 2, which differ only by a single -CH3 versus -H substituent on the central amine of the tridentate ligands binding to copper. Complex 1 with an end-on peroxido core and ferromagnetically (St = 1) coupled FeIII and CuII centers performs H-bonding-mediated O-O bond cleavage in the presence of phenol to generate oxoiron(IV) and exchange-coupled copper(II) and PhO• moieties. In contrast, the μ-η2:η1 peroxido complex 2, with a St = 0 ground state, is unreactive toward phenol. Thus, the implications for spin topology contributions to O-O bond cleavage, as proposed for the heme FeIII-(O22-)-CuII intermediate in CcO, can be extended to nonheme chemistry.

Reaction pathways, proton transfer, and proton pumping in ba3 class cytochrome c oxidase: perspectives from DFT quantum chemistry and molecular dynamics

After drawing comparisons between the reaction pathways of cytochrome c oxidase (CcO, Complex 4) and the preceding complex cytochrome bc1 (Complex 3), both being proton pumping complexes along the electron transport chain, we provide an analysis of the reaction pathways in bacterial ba3 class CcO, comparing spectroscopic results and kinetics observations with results from DFT calculations. For an important arc of the catalytic cycle in CcO, we can trace the energy pathways for the chemical protons and show how these pathways drive proton pumping of the vectorial protons. We then explore the proton loading network above the Fe heme a3-CuB catalytic center, showing how protons are loaded in and then released by combining DFT-based reaction energies with molecular dynamics simulations over states of that cycle. We also propose some additional reaction pathways for the chemical and vector protons based on our recent work with spectroscopic support.

New insights into the proton pumping mechanism of ba3 cytochrome c oxidase: the functions of key residues and water

As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.

Added Complexity!-Mechanistic Aspects of Heterobimetallic Complexes for Application in Homogeneous Catalysis

Inspired by multimetallic assemblies and their role in enzyme catalysis, chemists have developed a plethora of heterobimetallic complexes for application in homogeneous catalysis. Starting with small heterobimetallic complexes with σ-donating and π-accepting ligands, such as N-heterocyclic carbene and carbonyl ligands, more and more complex systems have been developed over the past two decades. These systems can show a significant increase in catalytic activity compared with their monometallic counterparts. This increase can be attributed to new reaction pathways enabled by the presence of a second metal center in the active catalyst. This review focuses on mechanistic aspects of heterobimetallic complexes in homogeneous catalysis. Depending on the type of interaction of the second metal with the substrates, heterobimetallic complexes can be subdivided into four classes. Each of these classes is illustrated with multiple examples, showcasing the versatility of both, the types of interactions possible, and the reactions accessible.

Second-Sphere Hydrogen-Bond Donors and Acceptors Affect the Rate and Selectivity of Electrochemical Oxygen Reduction by Iron Porphyrins Differently

The factors that control the rate and selectivity of 4e^-/4H⁺ O₂ reduction are important for efficient energy transformation as well as for understanding the terminal step of respiration in aerobic organisms. Inspired by the design of naturally occurring enzymes which are efficient catalysts for O₂ and H₂O₂ reduction, several artificial systems have been generated where different second-sphere residues have been installed to enhance the rate and efficiency of the 4e^-/4H⁺ O₂ reduction. These include hydrogen-bonding residues like amines, carboxylates, ethers, amides, phenols, etc. In some cases, improvements in the catalysis were recorded, whereas in some cases improvements were marginal or nonexistent. In this work, we use an iron porphyrin complex with pendant 1,10-phenanthroline residues which show a pH-dependent variation of the rate of the electrochemical O₂ reduction reaction (ORR) over 2 orders of magnitude. In-situ surface-enhanced resonance Raman spectroscopy reveals the presence of different intermediates at different pH's reflecting different rate-determining steps at different pH's. These data in conjunction with density functional theory calculations reveal that when the distal 1,10-phenanthroline is neutral it acts as a hydrogen-bond acceptor which stabilizes H₂O (product) binding to the active Fe^II state and retards the reaction. However, when the 1,10-phenanthroline is protonated, it acts as a hydrogen-bond donor which enhances O₂ reduction by stabilizing Fe^III-O₂^.- and Fe^III-OOH intermediates and activating the O-O bond for cleavage. On the basis of these data, general guidelines for controlling the different possible rate-determining steps in the complex multistep 4e^-/4H⁺ ORR are developed and a bioinspired principle-based design of an efficient electrochemical ORR is presented.

Second Sphere Effects on Oxygen Reduction and Peroxide Activation by Mononuclear Iron Porphyrins and Related Systems

Activation and reduction of O₂ and H₂O₂ by synthetic and biosynthetic iron porphyrin models have proved to be a versatile platform for evaluating second-sphere effects deemed important in naturally occurring heme active sites. Advances in synthetic techniques have made it possible to install different functional groups around the porphyrin ligand, recreating artificial analogues of the proximal and distal sites encountered in the heme proteins. Using judicious choices of these substituents, several of the elegant second-sphere effects that are proposed to be important in the reactivity of key heme proteins have been evaluated under controlled environments, adding fundamental insight into the roles played by these weak interactions in nature. This review presents a detailed description of these efforts and how these have not only demystified these second-sphere effects but also how the knowledge obtained resulted in functional mimics of these heme enzymes.

Mössbauer Property Calculations on Fea33+∙∙∙H2O∙∙∙CuB2+ Dinuclear Center Models of the Resting Oxidized as-Isolated Cytochrome c Oxidase

Mössbauer isomer shift and quadrupole splitting properties have been calculated using the OLYP‐D3(BJ) density functional method on previously obtained (W.‐G. Han Du, et al., Inorg Chem. 2020, 59, 8906–8915) geometry optimized Fe_a3³⁺−H₂O−Cu_B²⁺ dinuclear center (DNC) clusters of the resting oxidized (O state) “as‐isolated” cytochrome c oxidase (CcO). The calculated results are highly consistent with the available experimental observations. The calculations have also shown that the structural heterogeneities of the O state DNCs implicated by the Mössbauer experiments are likely consequences of various factors, particularly the variable positions of the central H₂O molecule between the Fe_a3³⁺ and Cu_B²⁺ sites in different DNCs, whether or not this central H₂O molecule has H‐bonding interaction with another H₂O molecule, the different spin states having similar energies for the Fe_a3³⁺ sites, and whether the Fe_a3³⁺ and Cu_B²⁺ sites are ferromagnetically or antiferromagnetically spin‐coupled.

Activation of O2 and NO in heme-copper oxidases - mechanistic insights from computational modelling

Heme-copper oxidases are transmembrane enzymes involved in aerobic and anaerobic respiration. The largest subgroup contains the cytochrome c oxidases (CcO), which reduce molecular oxygen to water. A significant part of the free energy released in this exergonic process is conserved as an electrochemical gradient across the membrane, via two processes, electrogenic chemistry and proton pumping. A deviant subgroup is the cytochrome c dependent NO reductases (cNOR), which reduce nitric oxide to nitrous oxide and water. This is also an exergonic reaction, but in this case none of the released free energy is conserved. Computational studies applying hybrid density functional theory to cluster models of the bimetallic active sites in the heme-copper oxidases are reviewed. To obtain a reliable description of the reaction mechanisms, energy profiles of the entire catalytic cycles, including the reduction steps have to be constructed. This requires a careful combination of computational results with certain experimental data. Computational studies have elucidated mechanistic details of the chemical parts of the reactions, involving cleavage and formation of covalent bonds, which have not been obtainable from pure experimental investigations. Important insights regarding the mechanisms of energy conservation have also been gained. The computational studies show that the reduction potentials of the active site cofactors in the CcOs are large enough to afford electrogenic chemistry and proton pumping, i.e. efficient energy conservation. These results solve a conflict between different types of experimental data. A mechanism for the proton pumping, involving a specific and crucial role for the active site tyrosine, conserved in all CcOs, is suggested. For the cNORs, the calculations show that the low reduction potentials of the active site cofactors are optimized for fast elimination of the toxic NO molecules. At the same time, the low reduction potentials lead to endergonic reduction steps with high barriers. To prevent even higher barriers, which would lead to a too slow reaction, when the electrochemical gradient across the membrane is present, the chemistry must occur in a non-electrogenic manner. This explains why there is no energy conservation in cNOR.

Coupled transport of electrons and protons in a bacterial cytochrome c oxidase-DFT calculated properties compared to structures and spectroscopies

After a general introduction to the features and mechanisms of cytochrome c oxidases (CcOs) in mitochondria and aerobic bacteria, we present DFT calculated physical and spectroscopic properties for the catalytic reaction cycle compared with experimental observations in bacterial ba3 type CcO, also with comparisons/contrasts to aa3 type CcOs. The Dinuclear Complex (DNC) is the active catalytic reaction center, containing a heme a3 Fe center and a near lying Cu center (called CuB) where by successive reduction and protonation, molecular O2 is transformed to two H2O molecules, and protons are pumped from an inner region across the membrane to an outer region by transit through the CcO integral membrane protein. Structures, energies and vibrational frequencies for Fe-O and O-O modes are calculated by DFT over the catalytic cycle. The calculated DFT frequencies in the DNC of CcO are compared with measured frequencies from Resonance Raman spectroscopy to clarify the composition, geometry, and electronic structures of different intermediates through the reaction cycle, and to trace reaction pathways. X-ray structures of the resting oxidized state are analyzed with reference to the known experimental reaction chemistry and using DFT calculated structures in fitting observed electron density maps. Our calculations lead to a new proposed reaction pathway for coupling the PR → F → OH (ferryl-oxo → ferric-hydroxo) pathway to proton pumping by a water shift mechanism. Through this arc of the catalytic cycle, major shifts in pKa's of the special tyrosine and a histidine near the upper water pool activate proton transfer. Additional mechanisms for proton pumping are explored, and the role of the CuB+ (cuprous state) in controlling access to the dinuclear reaction site is proposed.

A Water Molecule Residing in the Fea33+···CuB2+ Dinuclear Center of the Resting Oxidized as-Isolated Cytochrome c Oxidase: A Density Functional Study

Although the dinuclear center (DNC) of the resting oxidized "as-isolated" cytochrome c oxidase (CcO) is not a catalytically active state, its detailed structure, especially the nature of the bridging species between the Fea33+ and CuB2+ metal sites, is still both relevant and unsolved. Recent crystallographic work has shown an extended electron density for a peroxide type dioxygen species (O1-O2) bridging the Fea3 and CuB centers. In this paper, our density functional theory (DFT) calculations show that the observed peroxide type electron density between the two metal centers is most likely a mistaken analysis due to overlap of the electron density of a water molecule located at different positions between apparent O1 and O2 sites in DNCs of different CcO molecules with almost the same energy. Because the diffraction pattern and the resulting electron density map represent the effective long-range order averaged over many molecules and unit cells in the X-ray structure, this averaging can lead to an apparent observed superposition of different water positions between the Fea33+ and CuB2+ metal sites.

Proton-Electron Transfer to the Active Site Is Essential for the Reaction Mechanism of Soluble Δ9-Desaturase

A full understanding of the catalytic action of non-heme iron (NHFe) and non-heme diiron (NHFe₂) enzymes is still beyond the grasp of contemporary computational and experimental techniques. Many of these enzymes exhibit fascinating chemo-, regio-, and stereoselectivity, in spite of employing highly reactive intermediates which are necessary for activations of most stable chemical bonds. Herein, we study in detail one intriguing representative of the NHFe₂ family of enzymes: soluble Δ⁹ desaturase (Δ⁹D), which desaturates rather than performing the thermodynamically favorable hydroxylation of substrate. Its catalytic mechanism has been explored in great detail by using QM(DFT)/MM and multireference wave function methods. Starting from the spectroscopically observed 1,2-μ-peroxo diferric P intermediate, the proton-electron uptake by P is the favored mechanism for catalytic activation, since it allows a significant reduction of the barrier of the initial (and rate-determining) H-atom abstraction from the stearoyl substrate as compared to the "proton-only activated" pathway. Also, we ruled out that a Q-like intermediate (high-valent diamond-core bis-μ-oxo-[Fe^IV]₂ unit) is involved in the reaction mechanism. Our mechanistic picture is consistent with the experimental data available for Δ⁹D and satisfies fairly stringent conditions required by Nature: the chemo-, stereo-, and regioselectivity of the desaturation of stearic acid. Finally, the mechanisms evaluated are placed into a broader context of NHFe₂ chemistry, provided by an amino acid sequence analysis through the families of the NHFe₂ enzymes. Our study thus represents an important contribution toward understanding the catalytic action of the NHFe₂ enzymes and may inspire further work in NHFe₍₂₎ biomimetic chemistry.

Heme-Cu Binucleating Ligand Supports Heme/O2 and FeII-CuI/O2 Reactivity Providing High- and Low-Spin FeIII-Peroxo-CuII Complexes

The focus of this study is in the description of synthetic heme/copper/O₂ chemistry employing a heme-containing binucleating ligand which provides a tridentate chelate for copper ion binding. The addition of O₂ (-80 °C, tetrahydrofuran (THF) solvent) to the reduced heme compound (P^ImH)Fe^II (1), gives the oxy-heme adduct, formally a heme-superoxide complex Fe^III-(O₂^•-) (2) (resonance Raman spectroscopy (rR): ν_O-O, 1171 cm^-1 (Δ¹⁸O₂, -61 cm^-1); ν_Fe-O, 575 cm^-1 (Δ¹⁸O₂, -24 cm^-1)). Simple warming of 2 to room temperature regenerates reduced complex 1; this reaction is reversible, as followed by UV-vis spectroscopy. Complex 2 is electron paramagnetic resonance (EPR)-silent and exhibits upfield-shifted pyrrole resonances (δ 9.12 ppm) in ²H NMR spectroscopy, indicative of a six-coordinate low-spin heme. The coordination of the tethered imidazolyl arm to the heme-superoxide complex as an axial base ligand is suggested. We also report the new fully reduced heme-copper complex [(P^ImH)Fe^IICu^I]⁺ (3), where the copper ion is bound to the tethered tridentate portion of P^ImH. This reacts with O₂ to give a distinctive low-temperature-stable, high-spin (S = 2, overall) peroxo-bridged complex [(P^ImH)Fe^III-(O₂^2-)-Cu^II]⁺ (3a): λ_max, 420 (Soret), 545, 565 nm; δ_pyrr, 93 ppm; ν_O-O, 799 cm^-1 (Δ¹⁸O₂, -48 cm^-1); ν_Fe-O, 524 cm^-1 (Δ¹⁸O₂, -23 cm^-1). To 3a, the addition of dicyclohexylimidazole (DCHIm), which serves as a heme axial base, leads to low-spin (S = 0 overall) species complex [(DCHIm)(P^ImH)Fe^III-(O₂^2-)-Cu^II]⁺ (3b): λ_max, 425 (Soret), 538 nm; δ_pyrr, 10.2 ppm; ν_O-O, 817 cm^-1 (Δ¹⁸O₂, -55 cm^-1); ν_Fe-O, 610 cm^-1 (Δ¹⁸O₂, -26 cm^-1). These investigations into the characterization of the O₂-adducts from (P^ImH)Fe^II (1) with/without additional copper chelation advance our understanding of the dioxygen reactivity of heme-only and heme/Cu-ligand heterobinuclear system, thus potentially relevant to O₂ reduction in heme-copper oxidases or fuel-cell chemistry.

Importance of the iron-sulfur component and of the siroheme modification in the resting state of sulfite reductase

The active site of sulfite reductase (SiR) consists of an unusual siroheme-Fe₄S₄ assembly coupled via a cysteinate sulfur, and serves for multi-electron reduction reactions. Clear explanations have not been demonstrated for the reasons behind the choice of siroheme (vs. other types of heme) or for the single-atom coupling to an Fe₄S₄ center (as opposed to simple adjacency or to coupling via chains consisting of more than one atom). Possible explanations for these choices have previously been invoked, relating to the control of the spin state of the substrate-binding (siro)heme iron, modulation of the trans effect of the (Fe₄S₄-bound) cysteinate, or modulation of the redox potential. Reported here is a density functional theory (DFT) investigation of the structural interplay (in terms of geometry, molecular orbitals and magnetic interactions) between the siroheme and the Fe₄S₄ center as well as the importance of the covalent modifications within siroheme compared to the more common heme b, aiming to verify the role of the siroheme modification and of the Fe₄S₄ cluster at the SiR active site, with focus on previously-formulated hypotheses (geometrical/sterics, spin state, redox and electron-transfer control). A calibration of various DFT methods/variants for the correct description of ground state spin multiplicity is performed using a set of problematic cases of bioinorganic Fe centers; out of 11 functionals tested, M06-L and B3LYP offer the best results - though none of them correctly predict the spin state for all test cases. Upon examination of the relative energies of spin states, reduction potentials, energy decomposition (electrostatic, exchange-repulsion, orbital relaxation, correlation and dispersion interactions) and Mayer bond indices in SiR models, the following main roles of the siroheme and cubane are identified: (1) the cubane cofactor decreases the reduction potential of the siroheme and stabilizes the siroheme-cysteine bond interaction, and (2) the siroheme removes the quasi-degeneracy between the intermediate and high-spin states found in ferrous systems by preserving the latter as ground state; the higher-spin preference and the increased accessibility of multiple spin states are likely to be important in selective binding of the substrate and of the subsequent reaction intermediates, and in efficient changes in redox states throughout the catalytic cycle.

Geometric and Electronic Structure Contributions to O-O Cleavage and the Resultant Intermediate Generated in Heme-Copper Oxidases

This study investigates the mechanism of O-O bond cleavage in heme-copper oxidase (HCO) enzymes, combining experimental and computational insights from enzyme intermediates and synthetic models. It is determined that HCOs undergo a proton-initiated O-O cleavage mechanism where a single water molecule in the active site enables proton transfer (PT) from the cross-linked tyrosine to a peroxo ligand bridging the heme Fe^III and Cu^II, and multiple H-bonding interactions lower the tyrosine p K_a. Due to sterics within the active site, the proton must either transfer initially to the O(Fe) (a high-energy intermediate), or from another residue over a ∼10 Å distance to reach the O(Cu) atom directly. While the distance between the H⁺ donor (Tyr) and acceptor (O(Cu)) results in a barrier to PT, this separation is critical for the low barrier to O-O cleavage as it enhances backbonding from Fe into the O₂^2- σ* orbital. Thus, PT from Tyr precedes O-O elongation and is rate-limiting, consistent with available kinetic data. The electron transfers from tyrosinate after the barrier via a superexchange pathway provided by the cross-link, generating intermediate P_M. P_M is evaluated using available experimental data. The geometric structure contains an Fe^IV═O that is H-bonded to the Cu^II-OH. The electronic structure is a singlet, where the Fe^IV and Cu^II are antiferromagnetically coupled through the H-bond between the oxo(Fe) and hydroxo(Cu) ligands, while the Cu^II and Tyr^• are ferromagnetically coupled due their delocalization into orthogonal magnetic orbitals on the cross-linked His residue. These findings provide critical insights into the mechanism of efficient O₂ reduction in HCOs, and the nature of the P_M intermediate that couples this reaction to proton pumping.

Influence of intramolecular secondary sphere hydrogen-bonding interactions on cytochrome c oxidase inspired low-spin heme-peroxo-copper complexes

Dioxygen reduction by heme-copper oxidases is a critical biochemical process, wherein hydrogen bonding is hypothesized to participate in the critical step involving the active-site reductive cleavage of the O-O bond. Sixteen novel synthetic heme-(μ-O2 2-)-Cu(XTMPA) complexes, whose design is inspired by the cytochrome c oxidase active site structure, were generated in an attempt to form the first intramolecular H-bonded complexes. Derivatives of the "parent" ligand (XTMPA, TMPA = (tris((2-pyridyl)methyl)amine)) possessing one or two amine pendants preferentially form an H-bond with the copper-bound O-atom of the peroxide bridge. This is evidenced by a characteristic blue shift in the ligand-to-metal charge transfer (LMCT) bands observed in UV-vis spectroscopy (consistent with lowering of the peroxo π* relative to the iron orbitals) and a weakening of the O-O bond determined by resonance Raman spectroscopy (rR), with support from Density Functional Theory (DFT) calculations. Remarkably, with the TMPA-based infrastructure (versus similar heme-peroxo-copper complexes with different copper ligands), the typically undetected Cu-O stretch for these complexes was observed via rR, affording critical insights into the nature of the O-O peroxo core for the complexes studied. While amido functionalities have been shown to have greater H-bonding capabilities than their amino counterparts, in these heme-peroxo-copper complexes amido substituents distort the local geometry such that H-bonding with the peroxo core only imparts a weak electronic effect; optimal H-bonding interactions are observed by employing two amino groups on the copper ligand. The amino-substituted systems presented in this work reveal a key orientational anisotropy in H-bonding to the peroxo core for activating the O-O bond, offering critical insights into effective O-O cleavage chemistry. These findings indirectly support computational and protein structural studies suggesting the presence of an interstitial H-bonding water molecule in the CcO active site, which is critical for the desired reactivity. The results are evaluated with appropriate controls and discussed with respect to potential O2-reduction capabilities.

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.

A Water Dimer Shift Activates a Proton Pumping Pathway in the PR → F Transition of ba3 Cytochrome c Oxidase

Broken-symmetry density functional calculations have been performed on the [Fe_a3⁴⁺,Cu_B²⁺] state of the dinuclear center (DNC) for the P_R → F part of the catalytic cycle of ba₃ cytochrome c oxidase (CcO) from Thermus thermophilus (Tt), using the OLYP-D3-BJ functional. The calculations show that the movement of the H₂O molecules in the DNC affects the pK_a values of the residue side chains of Tyr237 and His376⁺, which are crucial for proton transfer/pumping in ba₃ CcO from Tt. The calculated lowest energy structure of the DNC in the [Fe_a3⁴⁺,Cu_B²⁺] state (state F) is of the form Fe_a3⁴⁺═O^2-···Cu_B²⁺, in which the H₂O ligand that resulted from protonation of the OH^- ligand in the P_R state is dissociated from the Cu_B²⁺ site. The calculated Fe_a3⁴⁺═O^2- distance in F (1.68 Å) is 0.03 Å longer than that in P_R (1.65 Å), which can explain the different Fe_a3⁴⁺═O^2- stretching modes in P (804 cm^-1) and F (785 cm^-1) identified by resonance Raman experiments. In this F state, the Cu_B²⁺···O^2- (ferryl-oxygen) distance is only around 2.4 Å. Hence, the subsequent O_H state [Fe_a3³⁺-OH^--Cu_B²⁺] with a μ-hydroxo bridge can be easily formed, as shown by our calculations.

Phenol-Induced O-O Bond Cleavage in a Low-Spin Heme-Peroxo-Copper Complex: Implications for O2 Reduction in Heme-Copper Oxidases

This study evaluates the reaction of a biomimetic heme-peroxo-copper complex, {[(DCHIm)(F8)FeIII]-(O22-)-[CuII(AN)]}+ (1), with a phenolic substrate, involving a net H-atom abstraction to cleave the bridging peroxo O-O bond that produces FeIV═O, CuII-OH, and phenoxyl radical moieties, analogous to the chemistry carried out in heme-copper oxidases (HCOs). A 3D potential energy surface generated for this reaction reveals two possible reaction pathways: one involves nearly complete proton transfer (PT) from the phenol to the peroxo ligand before the barrier; the other involves O-O homolysis, where the phenol remains H-bonding to the peroxo OCu in the transition state (TS) and transfers the H+ after the barrier. In both mechanisms, electron transfer (ET) from phenol occurs after the PT (and after the barrier); therefore, only the interaction with the H+ is involved in lowering the O-O cleavage barrier. The relative barriers depend on covalency (which governs ET from Fe), and therefore vary with DFT functional. However, as these mechanisms differ by the amount of PT at the TS, kinetic isotope experiments were conducted to determine which mechanism is active. It is found that the phenolic proton exhibits a secondary kinetic isotope effect, consistent with the calculations for the H-bonded O-O homolysis mechanism. The consequences of these findings are discussed in relation to O-O cleavage in HCOs, supporting a model in which a peroxo intermediate serves as the active H+ acceptor, and both the H+ and e- required for O-O cleavage derive from the cross-linked Tyr residue present at the active site.

Why copper is preferred over iron for oxygen activation and reduction in haem-copper oxidases

Haem-copper oxidase (HCO) catalyses the natural reduction of oxygen to water using a haem-copper centre. Despite decades of research on HCOs, the role of non-haem metal and the reason for nature's choice of copper over other metals such as iron remains unclear. Here, we use a biosynthetic model of HCO in myoglobin that selectively binds different non-haem metals to demonstrate 30-fold and 11-fold enhancements in the oxidase activity of Cu- and Fe-bound HCO mimics, respectively, as compared with Zn-bound mimics. Detailed electrochemical, kinetic and vibrational spectroscopic studies, in tandem with theoretical density functional theory calculations, demonstrate that the non-haem metal not only donates electrons to oxygen but also activates it for efficient O-O bond cleavage. Furthermore, the higher redox potential of copper and the enhanced weakening of the O-O bond from the higher electron density in the d orbital of copper are central to its higher oxidase activity over iron. This work resolves a long-standing question in bioenergetics, and renders a chemical-biological basis for the design of future oxygen-reduction catalysts.

Insights into proton translocation in cbb3 oxidase from MD simulations

Heme-copper oxidases are membrane protein complexes that catalyse the final step of the aerobic respiration, namely the reduction of oxygen to water. The energy released during catalysis is coupled to the active translocation of protons across the membrane, which contributes to the establishment of an electrochemical gradient that is used for ATP synthesis. The distinctive C-type (or cbb₃) cytochrome c oxidases, which are mostly present in proteobacteria, exhibit a number of unique structural and functional features, including high catalytic activity at low oxygen concentrations. At the moment, the functioning mechanism of C-type oxidases, in particular the proton transfer/pumping mechanism presumably via a single proton channel, is still poorly understood. In this work we used all-atom molecular dynamics simulations and continuum electrostatics calculations to obtain atomic-level insights into the hydration and dynamics of a cbb₃ oxidase. We provide the details of the water dynamics and proton transfer pathways for both the "chemical" and "pumped" protons, and show that formation of protonic connections is strongly affected by the protonation state of key residues, namely H243, E323 and H337.

Water exit pathways and proton pumping mechanism in B-type cytochrome c oxidase from molecular dynamics simulations

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. While proton uptake channels as well as water exit channels have been identified for A-type CcOs, the means by which water and protons exit B-type CcOs remain unclear. In this work, we investigate potential mechanisms for proton transport above the dinuclear center (DNC) in ba3-type CcO of Thermus thermophilus. Using long-time scale, all-atom molecular dynamics (MD) simulations for several relevant protonation states, we identify a potential mechanism for proton transport that involves propionate A of the active site heme a3 and residues Asp372, His376 and Glu126(II), with residue His376 acting as the proton-loading site. The proposed proton transport process involves a rotation of residue His376 and is in line with experimental findings. We also demonstrate how the strength of the salt bridge between residues Arg225 and Asp287 depends on the protonation state and that this salt bridge is unlikely to act as a simple electrostatic gate that prevents proton backflow. We identify two water exit pathways that connect the water pool above the DNC to the outer P-side of the membrane, which can potentially also act as proton exit transport pathways. Importantly, these water exit pathways can be blocked by narrowing the entrance channel between residues Gln151(II) and Arg449/Arg450 or by obstructing the entrance through a conformational change of residue Tyr136, respectively, both of which seem to be affected by protonation of residue His376.

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.

Data for molecular dynamics simulations of B-type cytochrome c oxidase with the Amber force field

Cytochrome c oxidase (CcO) is a vital enzyme that catalyzes the reduction of molecular oxygen to water and pumps protons across mitochondrial and bacterial membranes. This article presents parameters for the cofactors of ba₃-type CcO that are compatible with the all-atom Amber ff12SB and ff14SB force fields. Specifically, parameters were developed for the Cu_A pair, heme b, and the dinuclear center that consists of heme a₃ and Cu_B bridged by a hydroperoxo group. The data includes geometries in XYZ coordinate format for cluster models that were employed to compute proton transfer energies and derive bond parameters and point charges for the force field using density functional theory. Also included are the final parameter files that can be employed with the Amber leap program to generate input files for molecular dynamics simulations with the Amber software package. Based on the high resolution (1.8 Å) X-ray crystal structure of the ba₃-type CcO from Thermus thermophilus (Protein Data Bank ID number PDB: 3S8F), we built a model that is embedded in a POPC lipid bilayer membrane and solvated with TIP3P water molecules and counterions. We provide PDB data files of the initial model and the equilibrated model that can be used for further studies.

Dinickel complexes with anthyridine-based ligands

Two new dinickel complexes with anthyridine-based ligands, 5-phenyl-2,8-bis(2-pyridinyl)-1,9,10-anthyridine () and 5-phenyl-2,8-bis(6'-bipyridinyl)-1,9,10-anthyridine (), are reported. Complexation of Ni(OAc)2 with and in trifluoroacetic acid provided the corresponding dinickel complexes [{Ni2()(H2O)6(CF3COO)2}(CF3COO)2] () and [Ni2()(CF3COO)4(H2O)] (), respectively. Both complexes were characterized by spectroscopic methods and further confirmed by X-ray crystallography. Structural analyses of and revealed the NiNi distances in the complexes to be 5.4086(6) and 5.0138(7) Å, respectively. The catalytic activities of complexes and toward the reduction of carboxylic acids were evaluated. It appears that complex shows a good catalytic activity toward the reduction of carboxylates into the corresponding alcohols by diphenylsilane.

A broken-symmetry density functional study of structures, energies, and protonation states along the catalytic O-O bond cleavage pathway in ba3 cytochrome c oxidase from Thermus thermophilus

Broken-symmetry density functional calculations have been performed on the [Fea3, CuB] dinuclear center (DNC) of ba3 cytochrome c oxidase from Thermus thermophilus in the states of [Fea3(3+)-(HO2)(-)-CuB(2+), Tyr237(-)] and [Fea3(4+)[double bond, length as m-dash]O(2-), OH(-)-CuB(2+), Tyr237˙], using both PW91-D3 and OLYP-D3 functionals. Tyr237 is a special tyrosine cross-linked to His233, a ligand of CuB. The calculations have shown that the DNC in these states strongly favors the protonation of His376, which is above propionate-A, but not of the carboxylate group of propionate-A. The energies of the structures obtained by constrained geometry optimizations along the O-O bond cleavage pathway between [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] and [Fea3(4+)[double bond, length as m-dash]O(2-)HO(-)-CuB(2+), Tyr237˙] have also been calculated. The transition of [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] → [Fea3(4+)[double bond, length as m-dash]O(2-)HO(-)-CuB(2+), Tyr237˙] shows a very small barrier, which is less than 3.0/2.0 kcal mol(-1) in PW91-D3/OLYP-D3 calculations. The protonation state of His376 does not affect this O-O cleavage barrier. The rate limiting step of the transition from state A (in which O2 binds to Fea3(2+)) to state PM ([Fea3(4+)[double bond, length as m-dash]O(2-), OH(-)-CuB(2+), Tyr237˙], where the O-O bond is cleaved) in the catalytic cycle is, therefore, the proton transfer originating from Tyr237 to O-O to form the hydroperoxo [Fea3(3+)-(O-OH)(-)-CuB(2+), Tyr237(-)] state. The importance of His376 in proton uptake and the function of propionate-A/neutral-Asp372 as a gate to prevent the proton from back-flowing to the DNC are also shown.

Broken Symmetry DFT Calculations/Analysis for Oxidized and Reduced Dinuclear Center in Cytochrome c Oxidase: Relating Structures, Protonation States, Energies, and Mössbauer Properties in ba3 Thermus thermophilus

The Fea3(3+)···CuB(2+) dinuclear center (DNC) structure of the as-isolated oxidized ba3 cytochrome c oxidase (CcO) from Thermus thermophilus (Tt) is still not fully understood. When the proteins are initially crystallized in the oxidized state, they typically become radiolyticly reduced through X-ray irradiation. Several X-ray crystal structures of reduced ba3 CcO from Tt are available. However, depending on whether the crystals were prepared in a lipidic cubic phase environment or in detergent micelles, and whether the CcO's were chemically or radiolyticly reduced, the X-ray diffraction analysis of the crystals showed different Fea3(2+)···CuB(+) DNC structures. On the other hand, Mössbauer spectroscopic experiments on reduced and oxidized ba3 CcOs from Tt (Zimmermann et al., Proc. Natl. Acad. Sci. USA 1988, 85, 5779-5783) revealed multiple (57)Fea3(2+) and (57)Fea3(3+) components. Moreover, one of the (57)Fea3(3+) components observed at 4.2 K transformed from a proposed "low-spin" state to a different high-spin species when the temperature was increased above 190 K, whereas the other high-spin (57)Fea3(3+) component remained unchanged. In the current Article, in order to understand the heterogeneities of the DNC in both Mössbauer spectra and X-ray crystal structures, the spin crossover of one of the (57)Fea3(3+) components, and how the coordination and spin states of the Fea3(3+/2+) and Cu(2+/1+) sites relate to the heterogeneity of the DNC structures, we have applied density functional OLYP calculations to the DNC clusters established based on the different X-ray crystal structures of ba3 CcO from Tt. As a result, specific oxidized and reduced DNC structures related to the observed Mössbauer spectra and to spectral changes with temperature have been proposed. Our calculations also show that, in certain intermediate states, the His233 and His283 ligand side chains may dissociate from the CuB(+) site, and they may become potential proton loading sites during the catalytic cycle.

Proton transfer in the K-channel analog of B-type Cytochrome c oxidase from Thermus thermophilus

A key enzyme in aerobic metabolism is cytochrome c oxidase (CcO), which catalyzes the reduction of molecular oxygen to water in the mitochondrial and bacterial membranes. Substrate electrons and protons are taken up from different sides of the membrane and protons are pumped across the membrane, thereby generating an electrochemical gradient. The well-studied A-type CcO uses two different entry channels for protons: the D-channel for all pumped and two consumed protons, and the K-channel for the other two consumed protons. In contrast, the B-type CcO uses only a single proton input channel for all consumed and pumped protons. It has the same location as the A-type K-channel (and thus is named the K-channel analog) without sharing any significant sequence homology. In this study, we performed molecular-dynamics simulations and electrostatic calculations to characterize the K-channel analog in terms of its energetic requirements and functionalities. The function of Glu-15B as a proton sink at the channel entrance is demonstrated by its rotational movement out of the channel when it is deprotonated and by its high pKA value when it points inside the channel. Tyr-244 in the middle of the channel is identified as the valve that ensures unidirectional proton transfer, as it moves inside the hydrogen-bond gap of the K-channel analog only while being deprotonated. The electrostatic energy landscape was calculated for all proton-transfer steps in the K-channel analog, which functions via proton-hole transfer. Overall, the K-channel analog has a very stable geometry without large energy barriers.

Broken-Symmetry DFT Computations for the Reaction Pathway of IspH, an Iron-Sulfur Enzyme in Pathogenic Bacteria

The recently discovered methylerythritol phosphate (MEP) pathway provides new targets for the development of antibacterial and antimalarial drugs. In the final step of the MEP pathway, the [4Fe-4S] IspH protein catalyzes the 2e(-)/2H(+) reductive dehydroxylation of (E)-4-hydroxy-3-methyl-but-2-enyl diphosphate (HMBPP) to afford the isoprenoid precursors isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP). Recent experiments have attempted to elucidate the IspH catalytic mechanism to drive inhibitor development. Two competing mechanisms have recently emerged, differentiated by their proposed HMBPP binding modes upon 1e(-) reduction of the [4Fe-4S] cluster: (1) a Birch reduction mechanism, in which HMBPP remains bound to the [4Fe-4S] cluster through its terminal C4-OH group (ROH-bound) until the -OH is cleaved as water; and (2) an organometallic mechanism, in which the C4-OH group rotates away from the [4Fe-4S] cluster, allowing the HMBPP olefin group to form a metallacycle complex with the apical iron (η(2)-bound). We perform broken-symmetry density functional theory computations to assess the energies and reduction potentials associated with the ROH- and η(2)-bound states implicated by these competing mechanisms. Reduction potentials obtained for ROH-bound states are more negative (-1.4 to -1.0 V) than what is typically expected of [4Fe-4S] ferredoxin proteins. Instead, we find that η(2)-bound states are lower in energy than ROH-bound states when the [4Fe-4S] cluster is 1e(-) reduced. Furthermore, η(2)-bound states can already be generated in the oxidized state, yielding reduction potentials of ca. -700 mV when electron addition occurs after rotation of the HMBPP C4-OH group. We demonstrate that such η(2)-bound states are kinetically accessible both when the IspH [4Fe-4S] cluster is oxidized and 1e(-) reduced. The energetically preferred pathway gives 1e(-) reduction of the cluster after substrate conformational change, generating the 1e(-) reduced intermediate proposed in the organometallic mechanism.

Use of Broken-Symmetry Density Functional Theory To Characterize the IspH Oxidized State: Implications for IspH Mechanism and Inhibition

With current therapies becoming less efficacious due to increased drug resistance, new inhibitors of both bacterial and malarial targets are desperately needed. The recently discovered methylerythritol phosphate (MEP) pathway for isoprenoid synthesis provides novel targets for the development of such drugs. Particular attention has focused on the IspH protein, the final enzyme in the MEP pathway, which uses its [4Fe-4S] cluster to catalyze the formation of the isoprenoid precursors IPP and DMAPP from HMBPP. IspH catalysis is achieved via a 2e-/2H+ reductive dehydroxylation of HMBPP; the mechanism by which catalysis is achieved, however, is highly controversial. The work presented herein provides the first step in assessing different routes to catalysis by using computational methods. By performing broken-symmetry density functional theory (BS-DFT) calculations that employ both the conductor-like screening solvation model (DFT/COSMO) and a finite-difference Poisson-Boltzmann self-consistent reaction field methodology (DFT/SCRF), we evaluate geometries, energies, and Mössbauer signatures of the different protonation states that may exist in the oxidized state of the IspH catalytic cycle. From DFT/SCRF computations performed on the oxidized state, we find a state where the substrate, HMBPP, coordinates the apical iron in the [4Fe-4S] cluster as an alcohol group (ROH) to be one of two, isoenergetic, lowest-energy states. In this state, the HMBPP pyrophosphate moiety and an adjacent glutamate residue (E126) are both fully deprotonated, making the active site highly anionic. Our findings that this low-energy state also matches the experimental geometry of the active site and that its computed isomer shifts agree with experiment validate the use of the DFT/SCRF method to assess relative energies along the IspH reaction pathway. Additional studies of IspH catalytic intermediates are currently being pursued.