Cytochrome c oxidase: Intermediates of the catalytic cycle and their energy-coupled interconversion

[Pathogenetic basis of optic nerve atrophy in methanol poisoning]

Optic nerve atrophy is a pathomorphological consequence of diseases of the peripheral neuron of the visual pathway, manifested as atrophy of nerve fibers of varying severity. The toxic effect of methanol is mainly associated with formic acid and formaldehyde, which suppress the cytochrome system, inhibit oxidative phosphorylation, and thereby cause a deficiency of adenosine triphosphoric acid, to which brain and retinal tissues are especially susceptible. When formiate accumulates, tissue respiration is disrupted, leading to pronounced tissue hypoxia. As a result of such methanol metabolism, metabolic acidosis occurs. Tissue hypoxia develops in the first few hours as a result of the action of formic acid on the respiratory enzyme chain at the cytochrome oxidase level. Hypoxia and, as a consequence, a decrease in energy supply lead to a disruption of biological oxidation and the development of apoptosis in the optic nerve fibers. Understanding the process of optic nerve atrophy development at the pathogenetic level in methyl alcohol intoxication will help make a correct early diagnosis and prescribe timely treatment.

The catalytic reaction of cytochrome c oxidase probed by in situ gas titrations and FTIR difference spectroscopy

Cytochrome c oxidase (CcO) is a transmembrane heme‑copper metalloenzyme that catalyzes the reduction of O2 to H2O at the reducing end of the respiratory electron transport chain. To understand this reaction, we followed the conversion of CcO from Rhodobacter sphaeroides between several active-ready and carbon monoxide-inhibited states via attenuated total reflection Fourier-transform infrared (ATR FTIR) difference spectroscopy. Utilizing a novel gas titration setup, we prepared the mixed-valence, CO-inhibited R2CO state as well as the fully-reduced R4 and R4CO states and induced the "active ready" oxidized state OH. These experiments are performed in the dark yielding FTIR difference spectra exclusively triggered by exposure to O2, the natural substrate of CcO. Our data demonstrate that the presence of CO at heme a3 does not impair the catalytic oxidation of CcO when the cycle starts from the fully-reduced states. Interestingly, when starting from the R2CO state, the release of the CO ligand upon purging with inert gas yield a product that is indistinguishable from photolysis-induced states. The observed changes at heme a3 in the catalytic binuclear center (BNC) result from the loss of CO and are unrelated to electronic excitation upon illumination. Based on our experiments, we re-evaluate the assignment of marker bands that appear in time-resolved photolysis and perfusion-induced experiments on CcO.

Individual heme a and heme a3 contributions to the Soret absorption spectrum of the reduced bovine cytochrome c oxidase

Bovine cytochrome c oxidase (CcO) contains two hemes, a and a₃, chemically identical but differing in coordination and spin state. The Soret absorption band of reduced aa₃-type cytochrome c oxidase consists of overlapping bands of the hemes a²⁺ and a₃²⁺. It shows a peak at ∼444 nm and a distinct shoulder at ∼425 nm. However, attribution of individual spectral lineshapes to hemes a²⁺ and a₃²⁺ in the Soret is controversial. In the present work, we characterized spectral contributions of hemes a²⁺ and a₃²⁺ using two approaches. First, we reconstructed bovine CcO heme a²⁺ spectrum using a selective Ca²⁺-induced spectral shift of the heme a²⁺. Second, we investigated photobleaching of the reduced Thermus thermophilus ba₃- and bovine aa₃-oxidases in the Soret induced by femtosecond laser pulses in the Q-band. The resolved spectra show splitting of the electronic B_0x-, B_0y-transitions of both reduced hemes. The heme a²⁺ spectrum is shifted to the red relative to heme a₃²⁺ spectrum. The ∼425 nm shoulder is mostly attributed to heme a₃²⁺.

Interaction of Amphipathic Peptide from Influenza Virus M1 Protein with Mitochondrial Cytochrome Oxidase

The Bile Acid Binding Site (BABS) of cytochrome oxidase (CcO) binds numerous amphipathic ligands. To determine which of the BABS-lining residues are critical for interaction, we used the peptide P4 and its derivatives A1-A4. P4 is composed of two flexibly bound modified α-helices from the M1 protein of the influenza virus, each containing a cholesterol-recognizing CRAC motif. The effect of the peptides on the activity of CcO was studied in solution and in membranes. The secondary structure of the peptides was examined by molecular dynamics, circular dichroism spectroscopy, and testing the ability to form membrane pores. P4 was found to suppress the oxidase but not the peroxidase activity of solubilized CcO. The K_i(app) is linearly dependent on the dodecyl-maltoside (DM) concentration, indicating that DM and P4 compete in a 1:1 ratio. The true K_i is 3 μM. The deoxycholate-induced increase in K_i(app) points to a competition between P4 and deoxycholate. A1 and A4 inhibit solubilized CcO with K_i(app)~20 μM at 1 mM DM. A2 and A3 hardly inhibit CcO either in solution or in membranes. The mitochondrial membrane-bound CcO retains sensitivity to P4 and A4 but acquires resistance to A1. We associate the inhibitory effect of P4 with its binding to BABS and dysfunction of the proton channel K. Trp residue is critical for inhibition. The resistance of the membrane-bound enzyme to inhibition may be due to the disordered secondary structure of the inhibitory peptide.

Interactions of reactive sulfur species with metalloproteins

Reactive sulfur species (RSS) entail a diverse family of sulfur derivatives that have emerged as important effector molecules in H₂S-mediated biological events. RSS (including H₂S) can exert their biological roles via widespread interactions with metalloproteins. Metalloproteins are essential components along the metabolic route of oxygen in the body, from the transport and storage of O₂, through cellular respiration, to the maintenance of redox homeostasis by elimination of reactive oxygen species (ROS). Moreover, heme peroxidases contribute to immune defense by killing pathogens using oxygen-derived H₂O₂ as a precursor for stronger oxidants. Coordination and redox reactions with metal centers are primary means of RSS to alter fundamental cellular functions. In addition to RSS-mediated metalloprotein functions, the reduction of high-valent metal centers by RSS results in radical formation and opens the way for subsequent per- and polysulfide formation, which may have implications in cellular protection against oxidative stress and in redox signaling. Furthermore, recent findings pointed out the potential role of RSS as substrates for mitochondrial energy production and their cytoprotective capacity, with the involvement of metalloproteins. The current review summarizes the interactions of RSS with protein metal centers and their biological implications with special emphasis on mechanistic aspects, sulfide-mediated signaling, and pathophysiological consequences. A deeper understanding of the biological actions of reactive sulfur species on a molecular level is primordial in H₂S-related drug development and the advancement of redox medicine.

Kinetic Mathematical Modeling of Oxidative Phosphorylation in Cardiomyocyte Mitochondria

Oxidative phosphorylation (OXPHOS) is an oxygen-dependent process that consumes catabolized nutrients to produce adenosine triphosphate (ATP) to drive energy-dependent biological processes such as excitation-contraction coupling in cardiomyocytes. In addition to in vivo and in vitro experiments, in silico models are valuable for investigating the underlying mechanisms of OXPHOS and predicting its consequences in both physiological and pathological conditions. Here, we compare several prominent kinetic models of OXPHOS in cardiomyocytes. We examine how their mathematical expressions were derived, how their parameters were obtained, the conditions of their experimental counterparts, and the predictions they generated. We aim to explore the general landscape of energy production mechanisms in cardiomyocytes for future in silico models.

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.

Mitochondria: intracellular sentinels of infections

Structure and integrity of the mitochondrial network play important roles in many cellular processes. Loss of integrity can lead to the activation of a variety of signalling pathways and affect the cell’s response to infections. The activation of such mitochondria-mediated cellular responses has implications for infection recognition, signal transduction and pathogen control. Although we have a basic understanding of mitochondrial factors such as mitochondrial DNA or RNA that may be involved in processes like pro-inflammatory signalling, the diverse roles of mitochondria in host defence remain unclear. Here we will first summarise the functions of mitochondria in the host cell and provide an overview of the major known mitochondrial stress responses. We will then present recent studies that have contributed to the understanding of the role of mitochondria in infectious diseases and highlight a number of recently investigated models of bacterial and viral infections.

Direct Interaction of Mitochondrial Cytochrome c Oxidase with Thyroid Hormones: Evidence for Two Binding Sites

Thyroid hormones regulate tissue metabolism to establish an energy balance in the cell, in particular, by affecting oxidative phosphorylation. Their long-term impact is mainly associated with changes in gene expression, while the short-term effects may differ in their mechanisms. Our work was devoted to studying the short-term effects of hormones T2, T3 and T4 on mitochondrial cytochrome c oxidase (CcO) mediated by direct contact with the enzyme. The data obtained indicate the existence of two separate sites of CcO interaction with thyroid hormones, differing in their location, affinity and specificity to hormone binding. First, we show that T3 and T4 but not T2 inhibit the oxidase activity of CcO in solution and on membrane preparations with K_i ≈ 100–200 μM. In solution, T3 and T4 compete in a 1:1 ratio with the detergent dodecyl-maltoside to bind to the enzyme. The peroxidase and catalase partial activities of CcO are not sensitive to hormones, but electron transfer from heme a to the oxidized binuclear center is affected. We believe that T3 and T4 could be ligands of the bile acid-binding site found in the 3D structure of CcO by Ferguson-Miller’s group, and hormone-induced inhibition is associated with dysfunction of the K-proton channel. A possible role of this interaction in the physiological regulation of the enzyme is discussed. Second, we find that T2, T3, and T4 inhibit superoxide generation by oxidized CcO in the presence of excess H₂O₂. Inhibition is characterized by K_i values of 0.3–5 μM and apparently affects the formation of O₂^●− at the protein surface. The second binding site for thyroid hormones presumably coincides with the point of tight T2 binding on the Va subunit described in the literature.

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.

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.

A Resonance Raman Marker Band Characterizes the Slow and Fast Form of Cytochrome c Oxidase

Cytochrome c oxidase (CcO) in its as-isolated form is known to exist in a slow and fast form, which differ drastically in their ability to bind oxygen and other ligands. While preparation methods have been established that yield either the fast or the slow form of the protein, the underlying structural differences have not been identified yet. Here, we have performed surface enhanced resonance Raman (SERR) spectroscopy of CcO immobilized on electrodes in both forms. SERR spectra obtained in resonance with the heme a₃ metal-to-ligand charge transfer (MLCT) transition at 650 nm displayed a sharp vibrational band at 748 or 750 cm^-1 when the protein was in its slow or fast form, respectively. DFT calculations identified the band as a mode of the His-419 ligand that is sensitive to the oxygen ligand and the protonation state of Tyr-288 within the binuclear complex. Potential-dependent SERR spectroscopy showed a redox-induced change of this band around 525 mV versus Ag/AgCl exclusively for the fast form, which coincides with the redox potential of the Tyr-O/Tyr-O^- transition. Our data points to a peroxide ligand in the resting state of CcO for both forms. The observed frequencies and redox sensitivities of the Raman marker band suggest that a radical Tyr-288 is present in the fast form and a protonated Tyr-288 in the slow form.

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.

Acute toxicity of cyanide in aerobic respiration: Theoretical and experimental support for murburn explanation

The inefficiency of cyanide/HCN (CN) binding with heme proteins (under physiological regimes) is demonstrated with an assessment of thermodynamics, kinetics, and inhibition constants. The acute onset of toxicity and CN's mg/Kg LD50 (μM lethal concentration) suggests that the classical hemeFe binding-based inhibition rationale is untenable to account for the toxicity of CN. In vitro mechanistic probing of CN-mediated inhibition of hemeFe reductionist systems was explored as a murburn model for mitochondrial oxidative phosphorylation (mOxPhos). The effect of CN in haloperoxidase catalyzed chlorine moiety transfer to small organics was considered as an analogous probe for phosphate group transfer in mOxPhos. Similarly, inclusion of CN in peroxidase-catalase mediated one-electron oxidation of small organics was used to explore electron transfer outcomes in mOxPhos, leading to water formation. The free energy correlations from a Hammett study and IC50/Hill slopes analyses and comparison with ligands ( CO/ H 2 S/ N 3 - ) $\left( {\text{CO}}/{{{{\text{H}}_{2}}\text{S}}/{\text{N}_{3}^{\text{-}}}\;}\; \right)$ provide insights into the involvement of diffusible radicals and proton-equilibriums, explaining analogous outcomes in mOxPhos chemistry. Further, we demonstrate that superoxide (diffusible reactive oxygen species, DROS) enables in vitro ATP synthesis from ADP+phosphate, and show that this reaction is inhibited by CN. Therefore, practically instantaneous CN ion-radical interactions with DROS in matrix catalytically disrupt mOxPhos, explaining the acute lethal effect of CN.

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.

One small step for cytochrome P450 in its catalytic cycle, one giant leap for enzymology

The intermediates operating in the cytochrome P450 catalytic cycle have been investigated for more than half a century, fascinating many enzymologists. Each intermediate has its unique role to carry out diverse oxidations. Natural time course of the catalytic cycle is quite fast, hence, not all of the reactive intermediates could be isolated during physiological catalysis. Different high-valent iron intermediates have been proposed as primary oxidants: the candidates are compound 0 (Cpd 0, [FeOOH][Formula: see text]P450) and compound I (Cpd I, Fe(IV)[Formula: see text]O por[Formula: see text]P450). Among them, the role of Cpd I in hydroxylation is fairly well understood due the discovery of the peroxide shunt. This review endeavors to put the outstanding research efforts conducted to isolate and characterize the intermediates together. In addition to spectral features of each intermediate in the catalytic cycle, the oxidizing powers of Cpd 0 and Cpd I will be discussed along with most recent scientific findings.

Mitochondrial electron transport chain, ROS generation and uncoupling (Review)

The mammalian mitochondrial electron transport chain (ETC) includes complexes I-IV, as well as the electron transporters ubiquinone and cytochrome c. There are two electron transport pathways in the ETC: Complex I/III/IV, with NADH as the substrate and complex II/III/IV, with succinic acid as the substrate. The electron flow is coupled with the generation of a proton gradient across the inner membrane and the energy accumulated in the proton gradient is used by complex V (ATP synthase) to produce ATP. The first part of this review briefly introduces the structure and function of complexes I-IV and ATP synthase, including the specific electron transfer process in each complex. Some electrons are directly transferred to O₂ to generate reactive oxygen species (ROS) in the ETC. The second part of this review discusses the sites of ROS generation in each ETC complex, including sites I_F and I_Q in complex I, site II_F in complex II and site III_Qo in complex III, and the physiological and pathological regulation of ROS. As signaling molecules, ROS play an important role in cell proliferation, hypoxia adaptation and cell fate determination, but excessive ROS can cause irreversible cell damage and even cell death. The occurrence and development of a number of diseases are closely related to ROS overproduction. Finally, proton leak and uncoupling proteins (UCP_S) are discussed. Proton leak consists of basal proton leak and induced proton leak. Induced proton leak is precisely regulated and induced by UCPs. A total of five UCPs (UCP1-5) have been identified in mammalian cells. UCP1 mainly plays a role in the maintenance of body temperature in a cold environment through non-shivering thermogenesis. The core role of UCP2-5 is to reduce oxidative stress under certain conditions, therefore exerting cytoprotective effects. All diseases involving oxidative stress are associated with UCPs.

Investigating the Many Roles of Internal Water in Cytochrome c Oxidase

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain. As part of its catalytic cycle, C cO transfers protons to its Fe-Cu binuclear center (BNC) to reduce oxygen, and in addition, it pumps protons across the mitochondrial inner, or bacterial, membrane where it is located. It is believed that this proton transport is facilitated by a network of water chains inside the enzyme. Here we present an analysis of the hydration of C cO, including the BNC region, using a semi-empirical hydration program, Dowser++, recently developed in our group. Using high-resolution X-ray data, we show that Dowser++ predictions match very accurately the water molecules seen in the D- and K-channels of C cO, as well as in the vicinity of its BNC. Moreover, Dowser++ predicts many more internal water molecules than is typically seen in the experiment. However, no significant hydration of the catalytic cavity in C cO described recently in the literature is observed. As Dowser++ itself does not account for structural changes of the protein, this result supports the earlier assessment that the proposed wetting transition in the catalytic cavity can only either be due to structural rearrangements of BNC, possibly induced by the charges during the catalytic cycle, or occur transiently, in concert with the proton transfer. Molecular dynamics simulations were performed to investigate the global dynamic nature of Dowser++ waters in C cO, and the results suggest a consistent explanation as to why some predicted water molecules would be missing in the experimental structures. Furthermore, in light of the significant protein hydration predicted by Dowser++, the dielectric constant of the hydrated cavities in C cO was also investigated using the Fröhlich-Kirkwood model; the results indicate that in the cavities where water is packed sufficiently densely the dielectric constant can approach values comparable even to that of bulk water.

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.

Effect of Ca2+ on the redox potential of heme a in cytochrome c oxidase

Subunit I of cytochrome c oxidase (CcO) from mitochondria and many bacteria contains a cation binding site (CBS) located at the outer positively charged aqueous phase not far from heme a. Binding of Ca²⁺ with the CBS in bovine CcO inhibits activity of the enzyme 2-3 -fold [Vygodina, T., Kirichenko, A. & Konstantinov A.A. (2013) Direct Regulation of Cytochrome c Oxidase by Calcium Ions, PLoS One.8 e74436]. Here we show that binding of Ca²⁺ at CBS of bovine CcO shifts E_m of heme a to the positive by 15-20 mV. Na⁺ ions that bind to the same site and compete with Ca²⁺ do not affect E_m of heme a and also prevent and reverse the effect of Ca²⁺. No effect of Ca²⁺ or EGTA is observed on E_m of heme a with the wild type bacterial oxidases from R.sphaeroides or P.denitrificans that contain tightly-bound calcium at the site. In the D477A mutant CcO from P. denitrificans that binds Ca²⁺ reversibly like the mitochondrial CcO, calcium shifts redox titration curve of heme a to the positive by ∼35-50 mV that is in good agreement with the results of electrostatic calculations; however, as shown earlier, it does not inhibit CcO activity of the mutant enzyme. Therefore the data do not support the proposal that the inhibitory effect of Ca²⁺ on CcO activity may be explained by the Ca²⁺-induced shift of E_m of heme a. Rather, Ca²⁺ retards electron transfer by inhibition of charge dislocation in the exit part of the proton channel H in mammalian CcO, that is absent in the bacterial oxidases.

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.

Mitochondrial Respiratory Chain Complexes

Mitochondria are the power stations of the eukaryotic cell, using the energy released by the oxidation of glucose and other sugars to produce ATP. Electrons are transferred from NADH, produced in the citric acid cycle in the mitochondrial matrix, to oxygen by a series of large protein complexes in the inner mitochondrial membrane, which create a transmembrane electrochemical gradient by pumping protons across the membrane. The flow of protons back into the matrix via a proton channel in the ATP synthase leads to conformational changes in the nucleotide binding pockets and the formation of ATP. The three proton pumping complexes of the electron transfer chain are NADH-ubiquinone oxidoreductase or complex I, ubiquinone-cytochrome c oxidoreductase or complex III, and cytochrome c oxidase or complex IV. Succinate dehydrogenase or complex II does not pump protons, but contributes reduced ubiquinone. The structures of complex II, III and IV were determined by x-ray crystallography several decades ago, but complex I and ATP synthase have only recently started to reveal their secrets by advances in x-ray crystallography and cryo-electron microscopy. The complexes I, III and IV occur to a certain extent as supercomplexes in the membrane, the so-called respirasomes. Several hypotheses exist about their function. Recent cryo-electron microscopy structures show the architecture of the respirasome with near-atomic detail. ATP synthase occurs as dimers in the inner mitochondrial membrane, which by their curvature are responsible for the folding of the membrane into cristae and thus for the huge increase in available surface that makes mitochondria the efficient energy plants of the eukaryotic cell.

Splitting of the O-O bond at the heme-copper catalytic site of respiratory oxidases

Heme-copper oxidases catalyze the four-electron reduction of O₂ to H₂O at a catalytic site that is composed of a heme group, a copper ion (Cu_B), and a tyrosine residue. Results from earlier experimental studies have shown that the O-O bond is cleaved simultaneously with electron transfer from a low-spin heme (heme a/b), forming a ferryl state (P_R ; Fe⁴⁺=O^2-, Cu_B²⁺-OH^-). We show that with the Thermus thermophilus ba₃ oxidase, at low temperature (10°C, pH 7), electron transfer from the low-spin heme b to the catalytic site is faster by a factor of ~10 (τ ≅ 11 μs) than the formation of the P_R ferryl (τ ≅110 μs), which indicates that O₂ is reduced before the splitting of the O-O bond. Application of density functional theory indicates that the electron acceptor at the catalytic site is a high-energy peroxy state [Fe³⁺-O^--O^-(H⁺)], which is formed before the P_R ferryl. The rates of heme b oxidation and P_R ferryl formation were more similar at pH 10, indicating that the formation of the high-energy peroxy state involves proton transfer within the catalytic site, consistent with theory. The combined experimental and theoretical data suggest a general mechanism for O₂ reduction by heme-copper oxidases.

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.

Critical Aspects of Heme-Peroxo-Cu Complex Structure and Nature of Proton Source Dictate Metal-O(peroxo) Breakage versus Reductive O-O Cleavage Chemistry

The 4H+/4e- reduction of O2 to water, a key fuel-cell reaction also carried out in biology by oxidase enzymes, includes the critical O-O bond reductive cleavage step. Mechanistic investigations on active-site model compounds, which are synthesized by rational design to incorporate systematic variations, can focus on and resolve answers to fundamental questions, including protonation and/or H-bonding aspects, which accompany electron transfer. Here, we describe the nature and comparative reactivity of two low-spin heme-peroxo-Cu complexes, LS-4DCHIm, [(DCHIm)F8FeIII-(O22-)-CuII(DCHIm)4]+, and LS-3DCHIm, [(DCHIm)F8FeIII-(O22-)-CuII(DCHIm)3]+ (F8 = tetrakis(2,6-difluorophenyl)-porphyrinate; DCHIm = 1,5-dicyclohexylimidazole), toward different proton (4-nitrophenol and [DMF·H+](CF3SO3-)) (DMF = dimethyl-formamide) or electron (decamethylferrocene (Fc*)) sources. Spectroscopic reactivity studies show that differences in structure and electronic properties of LS-3DCHIm and LS-4DCHIm lead to significant differences in behavior. LS-3DCHIm is resistant to reduction, is unreactive toward weakly acidic 4-NO2-phenol, and stronger acids cleave the metal-O bonds, releasing H2O2. By contrast, LS-4DCHIm forms an adduct with 4-NO2-phenol, which includes an H-bond to the peroxo O-atom distal to Fe (resonance Raman (rR) spectroscopy and DFT). With addition of Fc* (2 equiv overall required), O-O reductive cleavage occurs, giving water, Fe(III), and Cu(II) products; however, a kinetic study reveals a one-electron rate-determining process, ket = 1.6 M-1 s-1 (-90 °C). The intermediacy of a high-valent [(DCHIm)F8FeIV═O] species is thus implied, and separate experiments show that one-electron reduction-protonation of [(DCHIm)F8FeIV═O] occurs faster (ket2 = 5.0 M-1 s-1), consistent with the overall postulated mechanism. The importance of the H-bonding interaction as a prerequisite for reductive cleavage is highlighted.

Modeling the detailed kinetics of mitochondrial cytochrome c oxidase: Catalytic mechanism and nitric oxide inhibition

Cytochrome c oxidase (CcO) catalyzes the exothermic reduction of O₂ to H₂O by using electrons from cytochrome c, and hence plays a crucial role in ATP production. Although details on the enzyme structure and redox centers involved in O₂ reduction have been known, there still remains a considerable ambiguity on its mechanism of action, e.g., the number of sequential electrons donated to O₂ in each catalytic step, the sites of protonation and proton pumping, and nitric oxide (NO) inhibition mechanism. In this work, we developed a thermodynamically constrained mechanistic mathematical model for the catalytic action of CcO based on available kinetic data. The model considers a minimal number of redox centers on CcO and couples electron transfer and proton pumping driven by proton motive force (PMF), and accounts for the inhibitory effects of NO on the reaction kinetics. The model is able to fit well all the available kinetic data under diverse experimental conditions with a physiologically realistic unique parameter set. The model predictions show that: 1) the apparent K_m of O₂ varies considerably and increases from fully reduced to fully oxidized cytochrome c depending on pH and the energy state of mitochondria, and 2) the intermediate enzyme states depend on pH and cytochrome c redox fraction and play a central role in coupling mitochondrial respiration to PMF. The developed CcO model can easily be integrated into existing mitochondrial bioenergetics models to understand the role of the enzyme in controlling oxidative phosphorylation in normal and disease conditions.

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.

Structural Changes and Proton Transfer in Cytochrome c Oxidase

In cytochrome c oxidase electron transfer from cytochrome c to O2 is linked to transmembrane proton pumping, which contributes to maintaining a proton electrochemical gradient across the membrane. The mechanism by which cytochrome c oxidase couples the exergonic electron transfer to the endergonic proton translocation is not known, but it presumably involves local structural changes that control the alternating proton access to the two sides of the membrane. Such redox-induced structural changes have been observed in X-ray crystallographic studies at residues 423-425 (in the R. sphaeroides oxidase), located near heme a. The aim of the present study is to investigate the functional effects of these structural changes on reaction steps associated with proton pumping. Residue Ser425 was modified using site-directed mutagenesis and time-resolved spectroscopy was used to investigate coupled electron-proton transfer upon reaction of the oxidase with O2. The data indicate that the structural change at position 425 propagates to the D proton pathway, which suggests a link between redox changes at heme a and modulation of intramolecular proton-transfer rates.

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

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.

Sulfide inhibition of and metabolism by cytochrome c oxidase

Hydrogen sulfide (H2S), a classic cytochrome c oxidase inhibitor, is also an in vitro oxidase substrate and an in vivo candidate hormonal ('gasotransmitter') species affecting sleep and hibernation. H2S, nitric oxide (NO) and carbon monoxide (CO) share some common features. All are low-molecular-mass physiological effectors and also oxidase inhibitors, capable of binding more than one enzyme site, and each is an oxidizable 'substrate'. The oxidase oxidizes CO to CO2, NO to nitrite and sulfide to probable persulfide species. Mitochondrial cytochrome c oxidase in an aerobic steady state with ascorbate and cytochrome c is rapidly inhibited by sulfide in a biphasic manner. At least two successive inhibited species are involved, probably partially reduced. The oxidized enzyme, in the absence of turnover, occurs in at least two forms: the 'pulsed' and 'resting' states. The pulsed form reacts aerobically with sulfide to form two intermediates, 'P' and 'F', otherwise involved in the reaction of oxygen with reduced enzyme. Sulfide can directly reduce the oxygen-reactive a3CuB binuclear centre in the pulsed state. The resting enzyme does not undergo such a step, but only a very slow one-electron reduction of the electron-transferring haem a. In final reactivation phases, both the steady-state inhibition of catalysis and the accumulation of P and F states are reversed by slow sulfide oxidation. A model for this complex reaction pattern is presented.

Cation binding site of cytochrome c oxidase: progress report

Cytochrome c oxidase from bovine heart binds Ca(2+) reversibly at a specific Cation Binding Site located near the outer face of the mitochondrial membrane. Ca(2+) shifts the absorption spectrum of heme a, which allowed earlier the determination of the kinetic and equilibrium characteristics of the binding, and, as shown recently, the binding of calcium to the site inhibits cytochrome oxidase activity at low turnover rates of the enzyme [Vygodina, Т., Kirichenko, A., Konstantinov, A.A (2013). Direct Regulation of Cytochrome c Oxidase by Calcium Ions. PloS ONE 8, e74436]. This paper summarizes further progress in the studies of the Cation Binding Site in this group presenting the results to be reported at 18th EBEC Meeting in Lisbon, 2014. The paper revises specificity of the bovine oxidase Cation Binding Site for different cations, describes dependence of the Ca(2+)-induced inhibition on turnover rate of the enzyme and reports very high affinity binding of calcium with the "slow" form of cytochrome oxidase. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference. Guest Editors: Manuela Pereira and Miguel Teixeira.

Density functional study for the bridged dinuclear center based on a high-resolution X-ray crystal structure of ba3 cytochrome c oxidase from Thermus thermophilus

Strong electron density for a peroxide type dioxygen species bridging the Fea3 and CuB dinuclear center (DNC) was observed in the high-resolution (1.8 Å) X-ray crystal structures (PDB entries 3S8G and 3S8F) of ba3 cytochrome c oxidase (CcO) from Thermus thermophilus. The crystals represent the as-isolated X-ray photoreduced CcO structures. The bridging peroxide was proposed to arise from the recombination of two radiation-produced HO(•) radicals formed either very near to or even in the space between the two metals of the DNC. It is unclear whether this peroxide species is in the O2(2-), O2(•)(-), HO2(-), or the H2O2 form and what is the detailed electronic structure and binding geometry including the DNC. In order to answer what form of this dioxygen species was observed in the DNC of the 1.8 Å X-ray CcO crystal structure (3S8G), we have applied broken-symmetry density functional theory (BS-DFT) geometric and energetic calculations (using OLYP potential) on large DNC cluster models with different Fea3-CuB oxidation and spin states and with O2(2-), O2(•)(-), HO2(-), or H2O2 in the bridging position. By comparing the DFT optimized geometries with the X-ray crystal structure (3S8G), we propose that the bridging peroxide is HO2(-). The X-ray crystal structure is likely to represent the superposition of the Fea3(2+)-(HO2(-))-CuB(+) DNC's in different states (Fe(2+) in low spin (LS), intermediate spin (IS), or high spin (HS)) with the majority species having the proton of the HO2(-) residing on the oxygen atom (O1) which is closer to the Fea3(2+) site in the Fea3(2+)-(HO-O)(-)-CuB(+) conformation. Our calculations show that the side chain of Tyr237 is likely trapped in the deprotonated Tyr237(-) anion form in the 3S8G X-ray crystal structure.

Functions of the hydrophilic channels in protonmotive cytochrome c oxidase

The structures and functions of hydrophilic channels in electron-transferring membrane proteins are discussed. A distinction is made between proton channels that can conduct protons and dielectric channels that are non-conducting but can dielectrically polarize in response to the introduction of charge changes in buried functional centres. Functions of the K, D and H channels found in A1-type cytochrome c oxidases are reviewed in relation to these ideas. Possible control of function by dielectric channels and their evolutionary relation to proton channels is explored.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

Efficient and selective isotopic labeling of hemes to facilitate the study of multiheme proteins

Specific isotopic labeling of hemes provides a unique opportunity to characterize the structure and function of heme-proteins. Unfortunately, current methods do not allow efficient labeling in high yields of multiheme cytochromes c, which are of great biotechnological interest. Here, a method for production of recombinant multiheme cytochromes c in Escherichia coli with isotopically labeled hemes is reported. A small tetraheme cytochrome of 12 kDa from Shewanella oneidensis MR-1 was used to demonstrate the method, achieving a production of 4 mg pure protein per liter. This method achieves, in a single step, efficient expression and incorporation of hemes isotopically labeled in specific atom positions adequate for spectroscopic characterization of these complex heme proteins. It is, furthermore, of general application to heme proteins, opening new possibilities for the characterization of this important class of proteins.