Active site intermediates in the reduction of O(2) by cytochrome oxidase, and their derivatives

Radical in the Peroxide-Produced F-Type Ferryl Form of Bovine Cytochrome c Oxidase

The reduction of O₂ in respiratory cytochrome c oxidases (CcO) is associated with the generation of the transmembrane proton gradient by two mechanisms. In one of them, the proton pumping, two different types of the ferryl intermediates of the catalytic heme a₃-Cu_B center P and F forms, participate. Equivalent ferryl states can be also formed by the reaction of the oxidized CcO (O) with H₂O₂. Interestingly, in acidic solutions a single molecule of H₂O₂ can generate from the O an additional F-type ferryl form (F^•) that should contain, in contrast to the catalytic F intermediate, a free radical at the heme a₃-Cu_B center. In this work, the formation and the endogenous decay of both the ferryl iron of heme a₃ and the radical in F^• intermediate were examined by the combination of four experimental approaches, isothermal titration calorimetry, electron paramagnetic resonance, and electronic absorption spectroscopy together with the reduction of this form by the defined number of electrons. The results are consistent with the generation of radicals in F^• form. However, the radical at the catalytic center is more rapidly quenched than the accompanying ferryl state of heme a₃, very likely by the intrinsic oxidation of the enzyme itself.

Temperature-dependent structural transition following X-ray-induced metal center reduction in oxidized cytochrome c oxidase

Cytochrome c oxidase (CcO) is the terminal enzyme in the electron transfer chain in the inner membrane of mitochondria. It contains four metal redox centers, two of which, Cu_B and heme a₃, form the binuclear center (BNC), where dioxygen is reduced to water. Crystal structures of CcO in various forms have been reported, from which ligand-binding states of the BNC and conformations of the protein matrix surrounding it have been deduced to elucidate the mechanism by which the oxygen reduction chemistry is coupled to proton translocation. However, metal centers in proteins can be susceptible to X-ray-induced radiation damage, raising questions about the reliability of conclusions drawn from these studies. Here, we used microspectroscopy-coupled X-ray crystallography to interrogate how the structural integrity of bovine CcO in the fully oxidized state (O) is modulated by synchrotron radiation. Spectroscopic data showed that, upon X-ray exposure, O was converted to a hybrid O∗ state where all the four metal centers were reduced, but the protein matrix was trapped in the genuine O conformation and the ligands in the BNC remained intact. Annealing the O∗ crystal above the glass transition temperature induced relaxation of the O∗ structure to a new R∗ structure, wherein the protein matrix converted to the fully reduced R conformation with the exception of helix X, which partly remained in the O conformation because of incomplete dissociation of the ligands from the BNC. We conclude from these data that reevaluation of reported CcO structures obtained with synchrotron light sources is merited.

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.

Critical roles of the CuB site in efficient proton pumping as revealed by crystal structures of mammalian cytochrome c oxidase catalytic intermediates

Mammalian cytochrome c oxidase (CcO) reduces O₂ to water in a bimetallic site including Fe_a3 and Cu_B giving intermediate molecules, termed A-, P-, F-, O-, E-, and R-forms. From the P-form on, each reaction step is driven by single-electron donations from cytochrome c coupled with the pumping of a single proton through the H-pathway, a proton-conducting pathway composed of a hydrogen-bond network and a water channel. The proton-gradient formed is utilized for ATP production by F-ATPase. For elucidation of the proton pumping mechanism, crystal structural determination of these intermediate forms is necessary. Here we report X-ray crystallographic analysis at ∼1.8 Å resolution of fully reduced CcO crystals treated with O₂ for three different time periods. Our disentanglement of intermediate forms from crystals that were composed of multiple forms determined that these three crystallographic data sets contained ∼45% of the O-form structure, ∼45% of the E-form structure, and ∼20% of an oxymyoglobin-type structure consistent with the A-form, respectively. The O- and E-forms exhibit an unusually long Cu_B²⁺-OH^- distance and Cu_B¹⁺-H₂O structure keeping Fe_a3³⁺-OH^- state, respectively, suggesting that the O- and E-forms have high electron affinities that cause the O→E and E→R transitions to be essentially irreversible and thus enable tightly coupled proton pumping. The water channel of the H-pathway is closed in the O- and E-forms and partially open in the R-form. These structures, together with those of the recently reported P- and F-forms, indicate that closure of the H-pathway water channel avoids back-leaking of protons for facilitating the effective proton pumping.

Axial Ligand Coordination Tuning of the Electrocatalytic Activity of Iron Porphyrin Electrografted onto Carbon Nanotubes for the Oxygen Reduction Reaction

The oxygen reduction reaction (ORR) is essential in many life processes and energy conversion systems. It is desirable to design transition metal molecular catalysts inspired by enzymatic oxygen activation/reduction processes as an alternative to noble-metal-Pt-based ORR electrocatalysts, especially in view point of fuel cell commercialization. We have fabricated bio-inspired molecular catalysts electrografted onto multiwalled carbon nanotubes (MWCNTs) in which 5,10,15,20-tetra(pentafluorophenyl) iron porphyrin (iron porphyrin FeF₂₀ TPP) is coordinated with covalently electrografted axial ligands varying from thiophene to imidazole on the MWCNTs' surface. The catalysts' electrocatalytic activity varied with the axial coordination environment (i. e., S-thiophene, N-imidazole, and O-carboxylate); the imidazole-coordinated catalyst MWCNTs-Im-FeF₂₀ TPP exhibited the highest ORR activity among the prepared catalysts. When MWCNT-Im-FeF₂₀ TPP was loaded onto the cathode of a zinc-air battery, an open-cell voltage (OCV) of 1.35 V and a maximum power density (P_max ) of 110 mW cm^-2 were achieved; this was higher than those of MWCNTs-Thi-FeF₂₀ TPP (OCV=1.30 V, P_max =100 mW cm^-2 ) and MWCNTs-Ox-FeF₂₀ TPP (OCV=1.28 V, P_max =86 mW cm^-2 ) and comparable with a commercial Pt/C catalyst (OCV=1.45 V, P_max =120 mW cm^-2 ) under similar experimental conditions. This study provides a time-saving method to prepare covalently immobilized molecular electrocatalysts on carbon-based materials with structure-performance correlation that is also applicable to the design of other electrografted catalysts for energy conversion.

Peroxide stimulated transition between the ferryl intermediates of bovine cytochrome c oxidase

During catalysis of cytochrome c oxidases (CcO) several ferryl intermediates of the catalytic heme a₃-Cu_B center are observed. In the P_M ferryl state, produced by the reaction of two-electron reduced CcO with O₂, the ferryl iron of heme a₃ and a free radical are present at the catalytic center. The radical reduction stimulates the transition of the P_M into another ferryl F state. Similar ferryl states can be also generated from the oxidized CcO (O) in the reaction with H₂O₂. The P_M, the product of the reaction of the O with one molecule of peroxide, is transformed into the F state by the second molecule of H₂O₂. However, the chemical nature of this transition has not been unambiguously elucidated yet. Here, we examined the redox state of the peroxide-produced P_M and F states by the one-electron reduction. The F form and interestingly also the major fraction of the P_M sample, likely another P-type ferryl form (P_R), were found to be the one oxidizing equivalent above the O state. However, the both P-type forms are transformed into the F state by additional molecule of H₂O₂. It is suggested that the P_R-to-F transition is due to the binding of H₂O₂ to Cu_B triggering a structural change together with the uptake of H⁺ at the catalytic center. In the P_M-to-F conversion, these two events are complemented with the annihilation of radical by the intrinsic oxidation of the enzyme.

Repurposing metalloproteins as mimics of natural metalloenzymes for small-molecule activation

Artificial metalloenzymes (ArMs) consist of an unnatural metal or cofactor embedded in a protein scaffold, and are an excellent platform for applying the concepts of protein engineering to catalysis. In this Focused Review, we describe the application of ArMs as simple, tunable artificial models of the active sites of complex natural metalloenzymes for small-molecule activation. In this sense, ArMs expand the strategies of synthetic model chemistry to protein-based supporting ligands with potential for participation from the second coordination sphere. We focus specifically on ArMs that are structural, spectroscopic, and functional models of enzymes for activation of small molecules like CO, CO₂, O₂, N₂, and NO, as well as production/consumption of H₂. These ArMs give insight into the identities and roles of metalloenzyme structural features within and near the cofactor. We give examples of ArM work relevant to hydrogenases, acetyl-coenzyme A synthase, superoxide dismutase, heme oxygenases, nitric oxide reductase, methyl-coenzyme M reductase, copper-O₂ enzymes, and nitrogenases.

Molecular understanding of heteronuclear active sites in heme-copper oxidases, nitric oxide reductases, and sulfite reductases through biomimetic modelling

Heme-copper oxidases (HCO), nitric oxide reductases (NOR), and sulfite reductases (SiR) catalyze the multi-electron and multi-proton reductions of O2, NO, and SO32-, respectively. Each of these reactions is important to drive cellular energy production through respiratory metabolism and HCO, NOR, and SiR evolved to contain heteronuclear active sites containing heme/copper, heme/nonheme iron, and heme-[4Fe-4S] centers, respectively. The complexity of the structures and reactions of these native enzymes, along with their large sizes and/or membrane associations, make it challenging to fully understand the crucial structural features responsible for the catalytic properties of these active sites. In this review, we summarize progress that has been made to better understand these heteronuclear metalloenzymes at the molecular level though study of the native enzymes along with insights gained from biomimetic models comprising either small molecules or proteins. Further understanding the reaction selectivity of these enzymes is discussed through comparisons of their similar heteronuclear active sites, and we offer outlook for further investigations.

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.

The redox-coupled proton-channel opening in cytochrome c oxidase

The interplay of cytochrome c oxidase's cofactor electrostatics, long-range conformational changes, H-bond rearrangement, and water dynamics enables transient proton-channel activation.

The Reactions of O2 and NO with Mixed-Valence ba3 Cytochrome c Oxidase from Thermus thermophilus

Earlier CO flow-flash experiments on the fully reduced Thermus thermophilus ba₃ (Tt ba₃) cytochrome oxidase revealed that O₂ binding was slowed down by a factor of 10 in the presence of CO (Szundi et al., 2010, PNAS 107, 21010-21015). The goal of the current study is to explore whether the long apparent lifetime (∼50 ms) of the Cu_B⁺-CO complex generated upon photolysis of the CO-bound mixed-valence Tt ba₃ (Koutsoupakis et al., 2019, Acc. Chem. Res. 52, 1380-1390) affects O₂ and NO binding and the ability of Cu_B to act as an electron donor during O-O bond splitting. The CO recombination, NO binding, and the reaction of mixed-valence Tt ba₃ with O₂ were investigated by time-resolved optical absorption spectroscopy using the CO flow-flash approach and photolabile O₂ and NO carriers. No electron backflow was detected after photolysis of the mixed-valence CO-bound Tt ba₃. The rate of O₂ and NO binding was two times slower than in the fully reduced enzyme in the presence of CO and 20 times slower than in the absence of CO. The purported long-lived Cu_B⁺-CO complex did not prevent O-O bond splitting and the resulting P_M formation, which was significantly faster (5-10 times) than in the bovine heart enzyme. We propose that O₂ binding to heme a₃ in Tt ba₃ causes CO to dissociate from Cu_B⁺ in a concerted manner through steric and/or electronic effects, thus allowing Cu_B⁺ to act as an electron donor in the mixed-valence enzyme. The significantly faster O₂ binding and O-O bond cleavage in Tt ba₃ compared to analogous steps in the aa₃ oxidases could reflect evolutionary adaptation of the enzyme to the microaerobic conditions of the T. thermophilus HB8 species.

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.

Silver Ions as a Tool for Understanding Different Aspects of Copper Metabolism

In humans, copper is an important micronutrient because it is a cofactor of ubiquitous and brain-specific cuproenzymes, as well as a secondary messenger. Failure of the mechanisms supporting copper balance leads to the development of neurodegenerative, oncological, and other severe disorders, whose treatment requires a detailed understanding of copper metabolism. In the body, bioavailable copper exists in two stable oxidation states, Cu(I) and Cu(II), both of which are highly toxic. The toxicity of copper ions is usually overcome by coordinating them with a wide range of ligands. These include the active cuproenzyme centers, copper-binding protein motifs to ensure the safe delivery of copper to its physiological location, and participants in the Cu(I) ↔ Cu(II) redox cycle, in which cellular copper is stored. The use of modern experimental approaches has allowed the overall picture of copper turnover in the cells and the organism to be clarified. However, many aspects of this process remain poorly understood. Some of them can be found out using abiogenic silver ions (Ag(I)), which are isoelectronic to Cu(I). This review covers the physicochemical principles of the ability of Ag(I) to substitute for copper ions in transport proteins and cuproenzyme active sites, the effectiveness of using Ag(I) to study copper routes in the cells and the body, and the limitations associated with Ag(I) remaining stable in only one oxidation state. The use of Ag(I) to restrict copper transport to tumors and the consequences of large-scale use of silver nanoparticles for human health are also discussed.

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.

Bacterial denitrifying nitric oxide reductases and aerobic respiratory terminal oxidases use similar delivery pathways for their molecular substrates

The superfamily of heme‑copper oxidoreductases (HCOs) include both NO and O₂ reductases. Nitric oxide reductases (NORs) are bacterial membrane enzymes that catalyze an intermediate step of denitrification by reducing nitric oxide (NO) to nitrous oxide (N₂O). They are structurally similar to heme‑copper oxygen reductases (HCOs), which reduce O₂ to water. The experimentally observed apparent bimolecular rate constant of NO delivery to the deeply buried catalytic site of NORs was previously reported to approach the diffusion-controlled limit (10⁸-10⁹ M^-1 s^-1). Using the crystal structure of cytochrome-c dependent NOR (cNOR) from Pseudomonas aeruginosa, we employed several protocols of molecular dynamics (MD) simulation, which include flooding simulations of NO molecules, implicit ligand sampling and umbrella sampling simulations, to elucidate how NO in solution accesses the catalytic site of this cNOR. The results show that NO partitions into the membrane, enters the enzyme from the lipid bilayer and diffuses to the catalytic site via a hydrophobic tunnel that is resolved in the crystal structures. This is similar to what has been found for O₂ diffusion through the closely related O₂ reductases. The apparent second order rate constant approximated using the simulation data is ~5 × 10⁸ M^-1 s^-1, which is optimized by the dynamics of the amino acid side chains lining in the tunnel. It is concluded that both NO and O₂ reductases utilize well defined hydrophobic tunnels to assure that substrate diffusion to the buried catalytic sites is not rate limiting under physiological conditions.

The electron distribution in the "activated" state of cytochrome c oxidase

Cytochrome c oxidase catalyzes reduction of O₂ to H₂O at a catalytic site that is composed of a copper ion and heme group. The reaction is linked to translocation of four protons across the membrane for each O₂ reduced to water. The free energy associated with electron transfer to the catalytic site is unequal for the four electron-transfer events. Most notably, the free energy associated with reduction of the catalytic site in the oxidized cytochrome c oxidase (state O) is not sufficient for proton pumping across the energized membrane. Yet, this electron transfer is mechanistically linked to proton pumping. To resolve this apparent discrepancy, a high-energy oxidized state (denoted O_H) was postulated and suggested to be populated only during catalytic turnover. The difference between states O and O_H was suggested to be manifested in an elevated midpoint potential of Cu_B in the latter. This proposal predicts that one-electron reduction of cytochrome c oxidase after its oxidation would yield re-reduction of essentially only Cu_B. Here, we investigated this process and found ~5% and ~6% reduction of heme a₃ and Cu_B, respectively, i.e. the apparent redox potentials for heme a₃ and Cu_B are lower than that of heme a.

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.

Investigation of dioxygen activation by copper(ii)–iminate/aminate complexes

Cu^II amidate/iminate complexes activate dioxygen by a ligated to Cu^II, –HCN– moiety.

The reductive phase of Rhodobacter sphaeroides cytochrome c oxidase disentangled by CO ligation

Cytochrome c oxidase (CcO) is a membrane protein of the respiratory chain that catalytically reduces molecular oxygen (O2) to water while translocating protons across the membrane. The enzyme hosts two copper and two heme iron moieties (heme a/heme a3). The atomic details of the sequential steps that go along with this redox-driven proton translocation are a matter of debate. Particularly for the reductive phase of CcO that precedes oxygen binding experimental data are scarce. Here, we use CcO under anaerobic conditions where carbon monoxide (CO) is bound to heme a3 which in tandem with CuB forms the binuclear center (BNC). Fourier-transform infrared (FTIR) absorption spectroscopy is combined with electro-chemistry to probe different redox and protonation states populated by variation of the external electrostatic potential. With this approach, the redox behavior of heme a and the BNC could be separated and the corresponding redox potentials were determined. We also infer the protonation of one of the propionate side chains of heme a3 to correlate with the oxidation of heme a. Experimental changes in the local electric field surrounding CO bound to heme a3 are determined by their vibrational Stark effect and agree well with electrostatic computations. The comparison of experimental and computational results indicates that changes of the heme a3/CuB redox state are coupled to proton transfer towards heme a3. The latter supports the role of the heme a3 propionate D as proton loading site.

Mitochondrial cytochrome c oxidase: catalysis, coupling and controversies

Mitochondrial cytochrome c oxidase is a member of a diverse superfamily of haem-copper oxidases. Its mechanism of oxygen reduction is reviewed in terms of the cycle of catalytic intermediates and their likely chemical structures. This reaction cycle is coupled to the translocation of protons across the inner mitochondrial membrane in which it is located. The likely mechanism by which this occurs, derived in significant part from studies of bacterial homologues, is presented. These mechanisms of catalysis and coupling, together with current alternative proposals of underlying mechanisms, are critically reviewed.

Factors Determining the Rate and Selectivity of 4e-/4H+ Electrocatalytic Reduction of Dioxygen by Iron Porphyrin Complexes

Reactivity as well as selectivity are crucial in the activation and electrocatalytic reduction of molecular oxygen. Recent developments in the understanding of the mechanism of electrocatalytic O₂ reduction by iron porphyrin complexes in situ using surface enhanced resonance Raman spectroscopy coupled to rotating disc electrochemistry (SERRS-RDE) in conjunction with H/D isotope effects on electrocatalytic current reveals that the rate of O₂ reduction, ∼10⁴ to 10⁵ M^-1 s^-1 for simple iron porphyrins, is limited by the rate of O-O bond cleavage of an intermediate ferric peroxide species (Fe^III-OOH). SERRS-RDE probes the system in operando when it is under steady state such that any intermediate species that has a greater rate of formation relative to its rate of decay, including the rate determining species, would accumulate and can be identified. This technique is particularly well suited to investigate iron porphyrin electrocatalysts as the intense symmetric ligand vibrations allow determination of the oxidation and spin states of the bound iron with high fidelity. The rate of O₂ reduction could be tuned up by 3 orders of magnitude by incorporating residues in the catalyst design that can exert "push" or "pull" effects, that is, axial phenolate and thiolate ligands and distal arginine residues. Similarly the rate of O-O bond cleavage can be enhanced by several orders of magnitude upon incorporating a distal Cu site and installing the active site in a hydrophobic protein environment in synthetic models and biosynthetic protein scaffolds. The selectivity, however, is solely determined by the site of protonation of a ferric peroxide (Fe^III-OOH) intermediate and can be governed by installing preorganized second sphere residues in the distal pocket. The 4e^-/4H⁺ reduction of O₂ entails protonation of the distal oxygen of the Fe^III-OOH species, while 2e^-/2H⁺ reduction requires the proximal oxygen to be protonated. Mechanistic investigations of CO₂ reduction by iron porphyrins reveal that the rate-determining step is the C-O bond cleavage of a Fe^II-COOH species analogous to the O-O bond cleavage step of a Fe^III-OOH species in O₂ reduction. The selectivity, resulting in either CO or HCOOH, is determined by the site of protonation of this species. These similarities suggests that the chemical principles governing the rate and selectivity of reduction of small molecules like O₂, CO₂, NO_x, and SO_x may be quite similar in nature.

Electron flow through biological molecules: does hole hopping protect proteins from oxidative damage?

Biological electron transfers often occur between metal-containing cofactors that are separated by very large molecular distances. Employing photosensitizer-modified iron and copper proteins, we have shown that single-step electron tunneling can occur on nanosecond to microsecond timescales at distances between 15 and 20 Å. We also have shown that charge transport can occur over even longer distances by hole hopping (multistep tunneling) through intervening tyrosines and tryptophans. In this perspective, we advance the hypothesis that such hole hopping through Tyr/Trp chains could protect oxygenase, dioxygenase, and peroxidase enzymes from oxidative damage. In support of this view, by examining the structures of P450 (CYP102A) and 2OG-Fe (TauD) enzymes, we have identified candidate Tyr/Trp chains that could transfer holes from uncoupled high-potential intermediates to reductants in contact with protein surface sites.

Kinetics of cytochrome c oxidase from R. sphaeroides initiated by direct electron transfer followed by tr-SEIRAS

Time-resolved surface-enhanced IR-absorption spectroscopy (tr-SEIRAS) has been performed on cytochrome c oxidase from Rhodobacter sphaeroides. The enzyme was converted electrochemically into the fully reduced state. Thereafter, in the presence of oxygen, the potential was switched to open circuit potential (OCP). Under these conditions, the enzyme is free to undergo enzymatic oxidation in the absence of an external electric field. Tr-SEIRAS was performed using the step-scan technique, triggered by periodic potential pulses switching between - 800mV and OCP. Single bands were resolved in a broad band in the amide I region using phase sensitive detection. Amplitudes of these bands were analyzed as a function of time. Time constants in the ms time scale were considered in terms of conformational changes of the protein secondary structures associated with the enzymatic turnover of the protein.

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.

All the O2 Consumed by Thermus thermophilus Cytochrome ba3 Is Delivered to the Active Site through a Long, Open Hydrophobic Tunnel with Entrances within the Lipid Bilayer

Cytochrome ba3 is a proton-pumping heme-copper oxygen reductase from the extreme thermophile Thermus thermophilus. Despite the fact that the enzyme's active site is buried deep within the protein, the apparent second order rate constant for the initial binding of O2 to the active-site heme has been experimentally found to be 10(9) M(-1) s(-1) at 298 K, at or near the diffusion limit, and 2 orders of magnitude faster than for O2 binding to myoglobin. To provide quantitative and microscopic descriptions of the O2 delivery pathway and mechanism in cytochrome ba3, extensive molecular dynamics simulations of the enzyme in its membrane-embedded form have been performed, including different protocols of explicit ligand sampling (flooding) simulations with O2, implicit ligand sampling analysis, and in silico mutagenesis. The results show that O2 diffuses to the active site exclusively via a Y-shaped hydrophobic tunnel with two 25-Å long membrane-accessible branches that coincide with the pathway previously suggested by the crystallographically identified xenon binding sites. The two entrances of the bifurcated tunnel of cytochrome ba3 are located within the lipid bilayer, where O2 is preferentially partitioned from the aqueous phase. The largest barrier to O2 migration within the tunnel is estimated to be only 1.5 kcal/mol, allowing O2 to reach the enzyme active site virtually impeded by one-dimensional diffusion once it reaches a tunnel entrance at the protein surface. Unlike other O2-utilizing proteins, the tunnel is "open" with no transient barriers observed due to protein dynamics. This unique low-barrier passage through the protein ensures that O2 transit through the protein is never rate-limiting.

The role of the K-channel and the active-site tyrosine in the catalytic mechanism of cytochrome c oxidase

The active site of cytochrome c oxidase (CcO) comprises an oxygen-binding heme, a nearby copper ion (CuB), and a tyrosine residue that is covalently linked to one of the histidine ligands of CuB. Two proton-conducting pathways are observed in CcO, namely the D- and the K-channels, which are used to transfer protons either to the active site of oxygen reduction (substrate protons) or for pumping. Proton transfer through the D-channel is very fast, and its role in efficient transfer of both substrate and pumped protons is well established. However, it has not been fully clear why a separate K-channel is required, apparently for the supply of substrate protons only. In this work, we have analysed the available experimental and computational data, based on which we provide new perspectives on the role of the K-channel. Our analysis suggests that proton transfer in the K-channel may be gated by the protonation state of the active-site tyrosine (Tyr244) and that the neutral radical form of this residue has a more general role in the CcO mechanism than thought previously. This article is part of a Special Issue entitled 'EBEC 2016: 19th European Bioenergetics Conference, Riva del Garda, Italy, July 2-6, 2016', edited by Prof. Paolo Bernardi.

A Designed Metalloenzyme Achieving the Catalytic Rate of a Native Enzyme

Terminal oxidases catalyze four-electron reduction of oxygen to water, and the energy harvested is utilized to drive the synthesis of adenosine triphosphate. While much effort has been made to design a catalyst mimicking the function of terminal oxidases, most biomimetic catalysts have much lower activity than native oxidases. Herein we report a designed oxidase in myoglobin with an O₂ reduction rate (52 s^–1) comparable to that of a native cytochrome (cyt) cbb₃ oxidase (50 s^–1) under identical conditions. We achieved this goal by engineering more favorable electrostatic interactions between a functional oxidase model designed in sperm whale myoglobin and its native redox partner, cyt b₅, resulting in a 400-fold electron transfer (ET) rate enhancement. Achieving high activity equivalent to that of native enzymes in a designed metalloenzyme offers deeper insight into the roles of tunable processes such as ET in oxidase activity and enzymatic function and may extend into applications such as more efficient oxygen reduction reaction catalysts for biofuel cells.

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.

Role of subunit III and its lipids in the molecular mechanism of cytochrome c oxidase

The terminal respiratory enzyme cytochrome c oxidase (CcO) reduces molecular oxygen to water, and pumps protons across the inner mitochondrial membrane, or the plasma membrane of bacteria. A two-subunit CcO harbors all the elements necessary for oxygen reduction and proton pumping. However, it rapidly undergoes turnover-induced irreversible damage, which is effectively prevented by the presence of subunit III and its tightly bound lipids. We have performed classical atomistic molecular dynamics (MD) simulations on a three-subunit CcO, which show the formation of water wires between the polar head groups of lipid molecules bound to subunit III and the proton uptake site Asp91 (Bos taurus enzyme numbering). Continuum electrostatic calculations suggest that these lipids directly influence the proton affinity of Asp91 by 1-2pK units. We surmise that lipids bound to subunit III influence the rate of proton uptake through the D-pathway, and therefore play a key role in preventing turnover-induced inactivation. Atomistic MD simulations show that subunit III is rapidly hydrated in the absence of internally bound lipids, which is likely to affect the rate of O2 diffusion into the active-site. The role of subunit III with its indigenous lipids in the molecular mechanism of CcO is discussed.

Defining the role of tyrosine and rational tuning of oxidase activity by genetic incorporation of unnatural tyrosine analogs

While a conserved tyrosine (Tyr) is found in oxidases, the roles of phenol ring pKa and reduction potential in O2 reduction have not been defined despite many years of research on numerous oxidases and their models. These issues represent major challenges in our understanding of O2 reduction mechanism in bioenergetics. Through genetic incorporation of unnatural amino acid analogs of Tyr, with progressively decreasing pKa of the phenol ring and increasing reduction potential, in the active site of a functional model of oxidase in myoglobin, a linear dependence of both the O2 reduction activity and the fraction of H2O formation with the pKa of the phenol ring has been established. By using these unnatural amino acids as spectroscopic probe, we have provided conclusive evidence for the location of a Tyr radical generated during reaction with H2O2, by the distinctive hyperfine splitting patterns of the halogenated tyrosines and one of its deuterated derivatives incorporated at the 33 position of the protein. These results demonstrate for the first time that enhancing the proton donation ability of the Tyr enhances the oxidase activity, allowing the Tyr analogs to augment enzymatic activity beyond that of natural Tyr.

Significant Improvement of Oxidase Activity through the Genetic Incorporation of a Redox-active Unnatural Amino Acid

Incorporation of 3-methoxytyrosine boosts the oxidase activity of the myoglobin model of oxidase, stressing the importance of the redox potential tuning of tyrosine.

While nature employs various covalent and non-covalent strategies to modulate tyrosine (Y) redox potential and pK_a in order to optimize enzyme activities, such approaches have not been systematically applied for the design of functional metalloproteins. Through the genetic incorporation of 3-methoxytyrosine (OMeY) into myoglobin, we replicated important features of cytochrome c oxidase (CcO) in this small soluble protein, which exhibits selective O₂ reduction activity while generating a small amount of reactive oxygen species (ROS). These results demonstrate that the electron donating ability of a tyrosine residue in the active site is important for CcO function. Moreover, we elucidated the structural basis for the genetic incorporation of OMeY into proteins by solving the X-ray structure of OMeY specific aminoacyl-tRNA synthetase complexed with OMeY.

Silica nanoparticles for the oriented encapsulation of membrane proteins into artificial bilayer lipid membranes

An artificial bilayer lipid membrane system is presented, featuring the oriented encapsulation of membrane proteins in a functionally active form. Nickel nitrilo-triacetic acid-functionalized silica nanoparticles, of a diameter of around 25 nm, are used to attach the proteins via a genetically engineered histidine tag in a uniform orientation. Subsequently, the proteins are reconstituted within a phospholipid bilayer, formed around the particles by in situ dialysis to form so-called proteo-lipobeads (PLBs). With a final size of about 50 nm, the PLBs can be employed for UV/vis spectroscopy studies, particularly of multiredox center proteins, because the effects of light scattering are negligible. As a proof of concept, we use cytochrome c oxidase (CcO) from P. denitrificans with the his tag genetically engineered to subunit I. In this orientation, the P side of CcO is directed to the outside and hence electron transfer can be initiated by reduced cytochrome c (cc). UV/vis measurements are used in order to determine the occupancy by CcO molecules encapsulated in the lipid bilayer as well as the kinetics of electron transfer between CcO and cc. The kinetic data are analyzed in terms of the Michaelis-Menten kinetics showing that the turnover rate of CcO is significantly decreased compared to that of solubilized protein, whereas the binding characteristics are improved. The data demonstrate the suitability of PLBs for functional cell-free bioassays of membrane proteins.

Proton-coupled electron transfer and the role of water molecules in proton pumping by cytochrome c oxidase

Molecular oxygen acts as the terminal electron sink in the respiratory chains of aerobic organisms. Cytochrome c oxidase in the inner membrane of mitochondria and the plasma membrane of bacteria catalyzes the reduction of oxygen to water, and couples the free energy of the reaction to proton pumping across the membrane. The proton-pumping activity contributes to the proton electrochemical gradient, which drives the synthesis of ATP. Based on kinetic experiments on the O-O bond splitting transition of the catalytic cycle (A → P(R)), it has been proposed that the electron transfer to the binuclear iron-copper center of O2 reduction initiates the proton pump mechanism. This key electron transfer event is coupled to an internal proton transfer from a conserved glutamic acid to the proton-loading site of the pump. However, the proton may instead be transferred to the binuclear center to complete the oxygen reduction chemistry, which would constitute a short-circuit. Based on atomistic molecular dynamics simulations of cytochrome c oxidase in an explicit membrane-solvent environment, complemented by related free-energy calculations, we propose that this short-circuit is effectively prevented by a redox-state-dependent organization of water molecules within the protein structure that gates the proton transfer pathway.

Prediction of high- and low-affinity quinol-analogue-binding sites in the aa3 and bo3 terminal oxidases from Bacillus subtilis and Escherichia coli1

Haem-copper oxidases are the terminal enzymes in both prokaryotic and eukaryotic respiratory chains. They catalyse the reduction of dioxygen to water and convert redox energy into a transmembrane electrochemical proton gradient during their catalytic activity. Haem-copper oxidases show substantial structure similarity, but spectroscopic and biochemical analyses indicate that these enzymes contain diverse prosthetic groups and use different substrates (i.e. cytochrome c or quinol). Owing to difficulties in membrane protein crystallization, there are no definitive structural data about the quinol oxidase physiological substrate-binding site(s). In the present paper, we propose an atomic structure model for the menaquinol:O2 oxidoreductase of Bacillus subtilis (QOx.aa3). Furthermore, a multistep computational approach is used to predict residues involved in the menaquinol/menaquinone binding within B. subtilis QOx.aa3 as well as those involved in quinol/quinone binding within Escherichia coli QOx.bo3. Two specific sequence motifs, R70GGXDX4RXQX3PX3FX[D/N/E/Q]X2HYNE97 and G159GSPX2GWX2Y169 (B. subtilis numbering), were highlighted within QOx from Bacillales. Specific residues within the first and the second sequence motif participate in the high- and low-affinity substrate-binding sites respectively. Using comparative analysis, two analogous motifs, R71GFXDX4RXQX8[Y/F]XPPHHYDQ101 and G163EFX3GWX2Y173 (E. coli numbering) were proposed to be involved in Enterobacteriales/Rhodobacterales/Rhodospirillales QOx high- and low-affinity quinol-derivative-binding sites. Results and models are discussed in the context of the literature.

Time-resolved surface-enhanced IR-absorption spectroscopy of direct electron transfer to cytochrome c oxidase from R. sphaeroides

Time-resolved surface-enhanced IR-absorption spectroscopy triggered by electrochemical modulation has been performed on cytochrome c oxidase from Rhodobacter sphaeroides. Single bands isolated from a broad band in the amide I region using phase-sensitive detection were attributed to different redox centers. Their absorbances changing on the millisecond timescale could be fitted to a model based on protonation-dependent chemical reaction kinetics established previously. Substantial conformational changes of secondary structures coupled to redox transitions were revealed.

How hydrogen peroxide is metabolized by oxidized cytochrome c oxidase

In the absence of external electron donors, oxidized bovine cytochrome c oxidase (CcO) exhibits the ability to decompose excess H2O2. Depending on the concentration of peroxide, two mechanisms of degradation were identified. At submillimolar peroxide concentrations, decomposition proceeds with virtually no production of superoxide and oxygen. In contrast, in the millimolar H2O2 concentration range, CcO generates superoxide from peroxide. At submillimolar concentrations, the decomposition of H2O2 occurs at least at two sites. One is the catalytic heme a3-CuB center where H2O2 is reduced to water. During the interaction of the enzyme with H2O2, this center cycles back to oxidized CcO via the intermediate presence of two oxoferryl states. We show that at pH 8.0 two molecules of H2O2 react with the catalytic center accomplishing one cycle. In addition, the reactions at the heme a3-CuB center generate the surface-exposed lipid-based radical(s) that participates in the decomposition of peroxide. It is also found that the irreversible decline of the catalytic activity of the enzyme treated with submillimolar H2O2 concentrations results specifically from the decrease in the rate of electron transfer from heme a to the heme a3-CuB center during the reductive phase of the catalytic cycle. The rates of electron transfer from ferrocytochrome c to heme a and the kinetics of the oxidation of the fully reduced CcO with O2 were not affected in the peroxide-modified CcO.