Proton transfer in ba3 cytochrome c oxidase from Thermus thermophilus

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

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

Membrane proteome of the thermoalkaliphile Caldalkalibacillus thermarum TA2.A1

Proteomics has greatly advanced the understanding of the cellular biochemistry of microorganisms. The thermoalkaliphile Caldalkalibacillus thermarum TA2.A1 is an organism of interest for studies into how alkaliphiles adapt to their extreme lifestyles, as it can grow from pH 7.5 to pH 11. Within most classes of microbes, the membrane-bound electron transport chain (ETC) enables a great degree of adaptability and is a key part of metabolic adaptation. Knowing what membrane proteins are generally expressed is crucial as a benchmark for further studies. Unfortunately, membrane proteins are the category of proteins hardest to detect using conventional cellular proteomics protocols. In part, this is due to the hydrophobicity of membrane proteins as well as their general lower absolute abundance, which hinders detection. Here, we performed a combination of whole cell lysate proteomics and proteomics of membrane extracts solubilised with either SDS or FOS-choline-12 at various temperatures. The combined methods led to the detection of 158 membrane proteins containing at least a single transmembrane helix (TMH). Within this data set we revealed a full oxidative phosphorylation pathway as well as an alternative NADH dehydrogenase type II (Ndh-2) and a microaerophilic cytochrome oxidase ba3. We also observed C. thermarum TA2.A1 expressing transporters for ectoine and glycine betaine, compounds that are known osmolytes that may assist in maintaining a near neutral internal pH when the external pH is highly alkaline.

A 2.2 Å cryoEM structure of a quinol-dependent NO Reductase shows close similarity to respiratory oxidases

Quinol-dependent nitric oxide reductases (qNORs) are considered members of the respiratory heme-copper oxidase superfamily, are unique to bacteria, and are commonly found in pathogenic bacteria where they play a role in combating the host immune response. qNORs are also essential enzymes in the denitrification pathway, catalysing the reduction of nitric oxide to nitrous oxide. Here, we determine a 2.2 Å cryoEM structure of qNOR from Alcaligenes xylosoxidans, an opportunistic pathogen and a denitrifying bacterium of importance in the nitrogen cycle. This high-resolution structure provides insight into electron, substrate, and proton pathways, and provides evidence that the quinol binding site not only contains the conserved His and Asp residues but also possesses a critical Arg (Arg720) observed in cytochrome bo3, a respiratory quinol oxidase.

Interaction of Terminal Oxidases with Amphipathic Molecules

The review focuses on recent advances regarding the effects of natural and artificial amphipathic compounds on terminal oxidases. Terminal oxidases are fascinating biomolecular devices which couple the oxidation of respiratory substrates with generation of a proton motive force used by the cell for ATP production and other needs. The role of endogenous lipids in the enzyme structure and function is highlighted. The main regularities of the interaction between the most popular detergents and terminal oxidases of various types are described. A hypothesis about the physiological regulation of mitochondrial-type enzymes by lipid-soluble ligands is considered.

Bioenergetics and Reactive Nitrogen Species in Bacteria

The production of reactive nitrogen species (RNS) by the innate immune system is part of the host’s defense against invading pathogenic bacteria. In this review, we summarize recent studies on the molecular basis of the effects of nitric oxide and peroxynitrite on microbial respiration and energy conservation. We discuss possible molecular mechanisms underlying RNS resistance in bacteria mediated by unique respiratory oxygen reductases, the mycobacterial bcc-aa₃ supercomplex, and bd-type cytochromes. A complete picture of the impact of RNS on microbial bioenergetics is not yet available. However, this research area is developing very rapidly, and the knowledge gained should help us develop new methods of treating infectious diseases.

Bioinorganic Chemistry on Electrodes: Methods to Functional Modeling

One of the major goals of bioinorganic chemistry has been to mimic the function of elegant metalloenzymes. Such functional modeling has been difficult to attain in solution, in particular, for reactions that require multiple protons and multiple electrons (nH⁺/ne^-). Using a combination of heterogeneous electrochemistry, electrode and molecule design one may control both electron transfer (ET) and proton transfer (PT) of these nH⁺/ne^- reactions. Such control can allow functional modeling of hydrogenases (H⁺ + e^- → 1/2 H₂), cytochrome c oxidase (O₂ + 4 e^- + 4 H⁺ → 2 H₂O), monooxygenases (RR'CH₂ + O₂ + 2 e^- + 2 H⁺ → RR'CHOH + H₂O) and dioxygenases (S + O₂ → SO₂; S = organic substrate) in aqueous medium and at room temperatures. In addition, these heterogeneous constructs allow probing unnatural bioinspired reactions and estimation of the inner- and outer-sphere reorganization energy of small molecules and proteins.

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.

Impact of Hydrogen Sulfide on Mitochondrial and Bacterial Bioenergetics

This review focuses on the effects of hydrogen sulfide (H2S) on the unique bioenergetic molecular machines in mitochondria and bacteria-the protein complexes of electron transport chains and associated enzymes. H2S, along with nitric oxide and carbon monoxide, belongs to the class of endogenous gaseous signaling molecules. This compound plays critical roles in physiology and pathophysiology. Enzymes implicated in H2S metabolism and physiological actions are promising targets for novel pharmaceutical agents. The biological effects of H2S are biphasic, changing from cytoprotection to cytotoxicity through increasing the compound concentration. In mammals, H2S enhances the activity of FoF1-ATP (adenosine triphosphate) synthase and lactate dehydrogenase via their S-sulfhydration, thereby stimulating mitochondrial electron transport. H2S serves as an electron donor for the mitochondrial respiratory chain via sulfide quinone oxidoreductase and cytochrome c oxidase at low H2S levels. The latter enzyme is inhibited by high H2S concentrations, resulting in the reversible inhibition of electron transport and ATP production in mitochondria. In the branched respiratory chain of Escherichia coli, H2S inhibits the bo3 terminal oxidase but does not affect the alternative bd-type oxidases. Thus, in E. coli and presumably other bacteria, cytochrome bd permits respiration and cell growth in H2S-rich environments. A complete picture of the impact of H2S on bioenergetics is lacking, but this field is fast-moving, and active ongoing research on this topic will likely shed light on additional, yet unknown biological effects.

Proton Pumping and Non-Pumping Terminal Respiratory Oxidases: Active Sites Intermediates of These Molecular Machines and Their Derivatives

Terminal respiratory oxidases are highly efficient molecular machines. These most important bioenergetic membrane enzymes transform the energy of chemical bonds released during the transfer of electrons along the respiratory chains of eukaryotes and prokaryotes from cytochromes or quinols to molecular oxygen into a transmembrane proton gradient. They participate in regulatory cascades and physiological anti-stress reactions in multicellular organisms. They also allow microorganisms to adapt to low-oxygen conditions, survive in chemically aggressive environments and acquire antibiotic resistance. To date, three-dimensional structures with atomic resolution of members of all major groups of terminal respiratory oxidases, heme-copper oxidases, and bd-type cytochromes, have been obtained. These groups of enzymes have different origins and a wide range of functional significance in cells. At the same time, all of them are united by a catalytic reaction of four-electron reduction in oxygen into water which proceeds without the formation and release of potentially dangerous ROS from active sites. The review analyzes recent structural and functional studies of oxygen reduction intermediates in the active sites of terminal respiratory oxidases, the features of catalytic cycles, and the properties of the active sites of these enzymes.

Rejigging Electron and Proton Transfer to Transition between Dioxygenase, Monooxygenase, Peroxygenase, and Oxygen Reduction Activity: Insights from Bioinspired Constructs of Heme Enzymes

Nature has employed heme proteins to execute a diverse set of vital life processes. Years of research have been devoted to understanding the factors which bias these heme enzymes, with all having a heme cofactor, toward distinct catalytic activity. Among them, axial ligation, distal super structure, and substrate binding pockets are few very vividly recognized ones. Detailed mechanistic investigation of these heme enzymes suggested that several of these enzymes, while functionally divergent, use similar intermediates. Furthermore, the formation and decay of these intermediates depend on proton and electron transfer processes in the enzyme active site. Over the past decade, work in this group, using in situ surface enhanced resonance Raman spectroscopy of synthetic and biosynthetic analogues of heme enzymes, a general idea of how proton and electron transfer rates relate to the lifetime of different O2 derived intermediates has been developed. These findings suggest that the enzymatic activities of all these heme enzymes can be integrated into one general cycle which can be branched out to different catalytic pathways by regulating the lifetime and population of each of these intermediates. This regulation can further be achieved by tuning the electron and proton transfer steps. By strategically populating one of these intermediates during oxygen reduction, one can navigate through different catalytic processes to a desired direction by altering proton and electron transfer steps.

ROS Defense Systems and Terminal Oxidases in Bacteria

Reactive oxygen species (ROS) comprise the superoxide anion (O₂^•−), hydrogen peroxide (H₂O₂), hydroxyl radical (^•OH), and singlet oxygen (¹O₂). ROS can damage a variety of macromolecules, including DNA, RNA, proteins, and lipids, and compromise cell viability. To prevent or reduce ROS-induced oxidative stress, bacteria utilize different ROS defense mechanisms, of which ROS scavenging enzymes, such as superoxide dismutases, catalases, and peroxidases, are the best characterized. Recently, evidence has been accumulating that some of the terminal oxidases in bacterial respiratory chains may also play a protective role against ROS. The present review covers this role of terminal oxidases in light of recent findings.

Specific inhibition of proton pumping by the T315V mutation in the K channel of cytochrome ba3 from Thermus thermophilus

Cytochrome ba₃ from Thermus thermophilus belongs to the B family of heme-copper oxidases and pumps protons across the membrane with an as yet unknown mechanism. The K channel of the A family heme-copper oxidases provides delivery of a substrate proton from the internal water phase to the binuclear heme-copper center (BNC) during the reductive phase of the catalytic cycle, while the D channel is responsible for transferring both substrate and pumped protons. By contrast, in the B family oxidases there is no D-channel and the structural equivalent of the K channel seems to be responsible for the transfer of both categories of protons. Here we have studied the effect of the T315V substitution in the K channel on the kinetics of membrane potential generation coupled to the oxidative half-reaction of the catalytic cycle of cytochrome ba₃. The results suggest that the mutated enzyme does not pump protons during the reaction of the fully reduced form with molecular oxygen in a single turnover. Specific inhibition of proton pumping in the T315V mutant appears to be a consequence of inability to provide rapid (τ ~ 100 μs) reprotonation of the internal transient proton donor(s) of the K channel. In contrast to the A family, the K channel of the B-type oxidases is necessary for the electrogenic transfer of both pumped and substrate protons during the oxidative half-reaction of the catalytic cycle.

The Redox-Active Tyrosine Is Essential for Proton Pumping in Cytochrome c Oxidase

Cellular respiration involves electron transport via a number of enzyme complexes to the terminal Cytochrome c oxidase (CcO), in which molecular oxygen is reduced to water. The free energy released in the reduction process is used to establish a transmembrane electrochemical gradient, via two processes, both corresponding to charge transport across the membrane in which the enzymes are embedded. First, the reduction chemistry occurring in the active site of CcO is electrogenic, which means that the electrons and protons are delivered from opposite sides of the membrane. Second, the exergonic chemistry is coupled to translocation of protons across the entire membrane, referred to as proton pumping. In the largest subfamily of the CcO enzymes, the A-family, one proton is pumped for every electron needed for the chemistry, making the energy conservation particularly efficient. In the present study, hybrid density functional calculations are performed on a model of the A-family CcOs. The calculations show that the redox-active tyrosine, conserved in all types of CcOs, plays an essential role for the energy conservation. Based on the calculations a reaction mechanism is suggested involving a tyrosyl radical (possibly mixed with tyrosinate character) in all reduction steps. The result is that the free energy released in each reduction step is large enough to allow proton pumping in all reduction steps without prohibitively high barriers when the gradient is present. Furthermore, the unprotonated tyrosine provides a mechanism for coupling the uptake of two protons per electron in every reduction step, i.e. for a secure proton pumping.

Nitric Oxide Does Not Inhibit but Is Metabolized by the Cytochrome bcc-aa3 Supercomplex

Nitric oxide (NO) is a well-known active site ligand and inhibitor of respiratory terminal oxidases. Here, we investigated the interaction of NO with a purified chimeric bcc-aa₃ supercomplex composed of Mycobacterium tuberculosis cytochrome bcc and Mycobacterium smegmatisaa₃-type terminal oxidase. Strikingly, we found that the enzyme in turnover with O₂ and reductants is resistant to inhibition by the ligand, being able to metabolize NO at 25 °C with an apparent turnover number as high as ≈303 mol NO (mol enzyme)⁻¹ min⁻¹ at 30 µM NO. The rate of NO consumption proved to be proportional to that of O₂ consumption, with 2.65 ± 0.19 molecules of NO being consumed per O₂ molecule by the mycobacterial bcc-aa₃. The enzyme was found to metabolize the ligand even under anaerobic reducing conditions with a turnover number of 2.8 ± 0.5 mol NO (mol enzyme)⁻¹ min⁻¹ at 25 °C and 8.4 µM NO. These results suggest a protective role of mycobacterial bcc-aa₃ supercomplexes against NO stress.

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.

Identifying the proton loading site cluster in the ba3 cytochrome c oxidase that loads and traps protons

Cytochrome c Oxidase (CcO) is the terminal electron acceptor in aerobic respiratory chain, reducing O₂ to water. The released free energy is stored by pumping protons through the protein, maintaining the transmembrane electrochemical gradient. Protons are held transiently in a proton loading site (PLS) that binds and releases protons driven by the electron transfer reaction cycle. Multi-Conformation Continuum Electrostatics (MCCE) was applied to crystal structures and Molecular Dynamics snapshots of the B-type Thermus thermophilus CcO. Six residues are identified as the PLS, binding and releasing protons as the charges on heme b and the binuclear center are changed: the heme a₃ propionic acids, Asp287, Asp372, His376 and Glu126B. The unloaded state has one proton and the loaded state two protons on these six residues. Different input structures, modifying the PLS conformation, show different proton distributions and result in different proton pumping behaviors. One loaded and one unloaded protonation states have the loaded/unloaded states close in energy so the PLS binds and releases a proton through the reaction cycle. The alternative proton distributions have state energies too far apart to be shifted by the electron transfers so are locked in loaded or unloaded states. Here the protein can use active states to load and unload protons, but has nearby trapped states, which stabilize PLS protonation state, providing new ideas about the CcO proton pumping mechanism. The distance between the PLS residues Asp287 and His376 correlates with the energy difference between loaded and unloaded states.

Proton transfer in uncoupled variants of cytochrome c oxidase

Cytochrome c oxidase is a membrane-bound redox-driven proton pump that harbors two proton-transfer pathways, D and K, which are used at different stages of the reaction cycle. Here, we address the question if a D pathway with a modified energy landscape for proton transfer could take over the role of the K pathway when the latter is blocked by a mutation. Our data indicate that structural alterations near the entrance of the D pathway modulate energy barriers that influence proton transfer to the proton-loading site. The data also suggest that during reduction of the catalytic site, its protonation has to occur via the K pathway and that this proton transfer to the catalytic site cannot take place through the D pathway.

The mechanism for oxygen reduction in the C family cbb3 cytochrome c oxidases - Implications for the proton pumping stoichiometry

Cytochrome c oxidases (CcOs) couple the exergonic reduction of molecular oxygen to proton pumping across the membrane in which they are embedded, thereby conserving a significant part of the free energy. The A family CcOs are known to pump four protons per oxygen molecule, while there is no consensus regarding the proton pumping stoichiometry for the C family cbb₃ oxidases. Hybrid density functional theory is used here to investigate the catalytic mechanism for oxygen reduction in cbb₃ oxidases. A surprising result is that the barrier for O O bond cleavage at the mixed valence reduction level seems to be too high compared to the overall reaction rate of the enzyme. It is therefore suggested that the O O bond is cleaved only after the first proton coupled reduction step, and that this reduction step most likely is not coupled to proton pumping. Furthermore, since the cbb₃ oxidases have only one proton channel leading to the active site, it is proposed that the activated E_H intermediate, suggested to be responsible for proton pumping in one of the reduction steps in the A family, cannot be involved in the catalytic cycle for cbb₃, which results in the lack of proton pumping also in the E to R reduction step. In summary, the calculations indicate that only two protons are pumped per oxygen molecule in cbb₃ oxidases. However, more experimental information on this divergent enzyme is needed, e.g. whether the flow of electrons resembles that in the other more well-studied CcO families.

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.

Active Site Midpoint Potentials in Different Cytochrome c Oxidase Families: A Computational Comparison

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain, reducing molecular oxygen to water. The binuclear active site in C cO comprises a high-spin heme associated with a Cu_B complex and a redox active tyrosine. The electron transport in the respiratory chain is driven by increasing midpoint potentials of the involved cofactors, resulting in a release of free energy, which is stored by coupling the electron transfer to proton translocation across a membrane, building up an electrochemical gradient. In this context, the midpoint potentials of the active site cofactors in the C cOs are of special interest, since they determine the driving forces for the individual oxygen reduction steps and thereby affect the efficiency of the proton pumping. It has been difficult to obtain useful information on some of these midpoint potentials from experiments. However, since each of the reduction steps in the catalytic cycle of oxygen reduction to water corresponds to the formation of an O-H bond, they can be calculated with a reasonably high accuracy using quantum chemical methods. From the calculated O-H bond strengths, the proton-coupled midpoint potentials of the active site cofactors can be estimated. Using models representing the different families of C cO's (A, B, and C), the calculations give midpoint potentials that should be relevant during catalytic turnover. The calculations also suggest possible explanations for why some experimentally measured potentials deviate significantly from the calculated ones, i.e., for Cu_B in all oxidase families, and for heme b₃ in the C family.

Control of transmembrane charge transfer in cytochrome c oxidase by the membrane potential

The respiratory chain in mitochondria is composed of membrane-bound proteins that couple electron transfer to proton translocation across the inner membrane. These charge-transfer reactions are regulated by the proton electrochemical gradient that is generated and maintained by the transmembrane charge transfer. Here, we investigate this feedback mechanism in cytochrome c oxidase in intact inner mitochondrial membranes upon generation of an electrochemical potential by hydrolysis of ATP. The data indicate that a reaction step that involves proton uptake to the catalytic site and presumably proton translocation is impaired by the potential, but electron transfer is not affected. These results define the order of electron and proton-transfer reactions and suggest that the proton pump is regulated by the transmembrane electrochemical gradient through control of internal proton transfer rather than by control of electron transfer.

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.

Oxygen Activation and Energy Conservation by Cytochrome c Oxidase

This review focuses on the type A cytochrome c oxidases (CcO), which are found in all mitochondria and also in several aerobic bacteria. CcO catalyzes the respiratory reduction of dioxygen (O₂) to water by an intriguing mechanism, the details of which are fairly well understood today as a result of research for over four decades. Perhaps even more intriguingly, the membrane-bound CcO couples the O₂ reduction chemistry to translocation of protons across the membrane, thus contributing to generation of the electrochemical proton gradient that is used to drive the synthesis of ATP as catalyzed by the rotary ATP synthase in the same membrane. After reviewing the structure of the core subunits of CcO, the active site, and the transfer paths of electrons, protons, oxygen, and water, we describe the states of the catalytic cycle and point out the few remaining uncertainties. Finally, we discuss the mechanism of proton translocation and the controversies in that area that still prevail.

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.

Cavity hydration dynamics in cytochrome c oxidase and functional implications

Cytochrome c oxidase (CcO) is a transmembrane protein that uses the free energy of O₂ reduction to generate the proton concentration gradient across the membrane. The regulation of competitive proton transfer pathways has been established to be essential to the vectorial transport efficiency of CcO, yet the underlying mechanism at the molecular level remains lacking. Recent studies have highlighted the potential importance of hydration-level change in an internal cavity that connects the proton entrance channel, the site of O₂ reduction, and the putative proton exit route. In this work, we use atomistic molecular dynamics simulations to investigate the energetics and timescales associated with the volume fluctuation and hydration-level change in this central cavity. Extensive unrestrained molecular dynamics simulations (accumulatively [Formula: see text]4 [Formula: see text]s) and free energy computations for different chemical states of CcO support a model in which the volume and hydration level of the cavity are regulated by the protonation state of a propionate group of heme a₃ and, to a lesser degree, the redox state of heme a and protonation state of Glu286. Markov-state model analysis of [Formula: see text]2-[Formula: see text]s trajectories suggests that hydration-level change occurs on the timescale of 100-200 ns before the proton-loading site is protonated. The computed energetic and kinetic features for the cavity wetting transition suggest that reversible hydration-level change of the cavity can indeed be a key factor that regulates the branching of proton transfer events and therefore contributes to the vectorial efficiency of proton transport.

Time-resolved generation of membrane potential by ba3 cytochrome c oxidase from Thermus thermophilus coupled to single electron injection into the O and OH states

Two electrogenic phases with characteristic times of ~14μs and ~290μs are resolved in the kinetics of membrane potential generation coupled to single-electron reduction of the oxidized "relaxed" O state of ba₃ oxidase from T. thermophilus (O→E transition). The rapid phase reflects electron redistribution between Cu_A and heme b. The slow phase includes electron redistribution from both Cu_A and heme b to heme a₃, and electrogenic proton transfer coupled to reduction of heme a₃. The distance of proton translocation corresponds to uptake of a proton from the inner water phase into the binuclear center where heme a₃ is reduced, but there is no proton pumping and no reduction of Cu_B. Single-electron reduction of the oxidized "unrelaxed" state (O_H→E_H transition) is accompanied by electrogenic reduction of the heme b/heme a₃ pair by Cu_A in a "fast" phase (~22μs) and transfer of protons in "middle" and "slow" electrogenic phases (~0.185ms and ~0.78ms) coupled to electron redistribution from the heme b/heme a₃ pair to the Cu_B site. The "middle" and "slow" electrogenic phases seem to be associated with transfer of protons to the proton-loading site (PLS) of the proton pump, but when all injected electrons reach Cu_B the electronic charge appears to be compensated by back-leakage of the protons from the PLS into the binuclear site. Thus proton pumping occurs only to the extent of ~0.1 H⁺/e^-, probably due to the formed membrane potential in the experiment.

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.

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

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

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

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

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

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

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

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

A method aimed at assessing the functional consequences of the supramolecular organization of the respiratory electron transfer chain by time-resolved studies

A steadily increasing number of physiological, biochemical, and structural studies have provided a growing support to the notion that the respiratory electron transfer chain may contain supra-molecular edifices made of the assembly of some, if not all, of its individual links. This structure, usually referred to as the solid state model-in comparison to the liquid state model in which the electron transfer reactions between the membrane bound enzymes are diffusion controlled-is seen as conferring specific kinetic properties to the chain and thus as being highly relevant from a functional point of view. Although the assumption that structural changes are mirrored by functional adjustment is undoubtedly legitimate, experimental evidences supporting it remain scarce. Here we review a recent methodological development aimed at tackling the functional relevance of the supramolecular organization of the respiratory electron transfer chain in intact cells.

Moving protons and electrons in biomimetic systems

An enormous variety of biological redox reactions are accompanied by changes in proton content at enzyme active sites, in their associated cofactors, in substrates and/or products, and between protein interfaces. Understanding this breadth of reactivity is an ongoing chemical challenge. A great many workers have developed and investigated biomimetic model complexes to build new ways of thinking about the mechanistic underpinnings of such complex biological proton-coupled electron transfer (PCET) reactions. Of particular importance are those model reactions that involve transfer of one proton (H(+)) and one electron (e(-)), which is equivalent to transfer of a hydrogen atom (H(•)). In this Current Topic, we review key concepts in PCET reactivity and describe important advances in biomimetic PCET chemistry, with a special emphasis on research that has enhanced efforts to understand biological PCET reactions.

Mutation of a single residue in the ba3 oxidase specifically impairs protonation of the pump site

The ba3-type cytochrome c oxidase from Thermus thermophilus is a membrane-bound protein complex that couples electron transfer to O2 to proton translocation across the membrane. To elucidate the mechanism of the redox-driven proton pumping, we investigated the kinetics of electron and proton transfer in a structural variant of the ba3 oxidase where a putative "pump site" was modified by replacement of Asp372 by Ile. In this structural variant, proton pumping was uncoupled from internal electron transfer and O2 reduction. The results from our studies show that proton uptake to the pump site (time constant ∼65 μs in the wild-type cytochrome c oxidase) was impaired in the Asp372Ile variant. Furthermore, a reaction step that in the wild-type cytochrome c oxidase is linked to simultaneous proton uptake and release with a time constant of ∼1.2 ms was slowed to ∼8.4 ms, and in Asp372Ile was only associated with proton uptake to the catalytic site. These data identify reaction steps that are associated with protonation and deprotonation of the pump site, and point to the area around Asp372 as the location of this site in the ba3 cytochrome c oxidase.

A structural and functional perspective on the evolution of the heme-copper oxidases

The heme-copper oxidases (HCOs) catalyze the reduction of O2 to water, and couple the free energy to proton pumping across the membrane. HCOs are divided into three sub-classes, A, B and C, whose order of emergence in evolution has been controversial. Here we have analyzed recent structural and functional data on HCOs and their homologues, the nitric oxide reductases (NORs). We suggest that the C-type oxidases are ancient enzymes that emerged from the NORs. In contrast, the A-type oxidases are the most advanced from both structural and functional viewpoints, which we interpret as evidence for having evolved later.

Evidence for distinct electron transfer processes in terminal oxidases from different origin by means of protein film voltammetry

Cytochrome aa₃ from Paracoccus denitrificans and cytochrome ba₃ from Thermus thermophilus, two distinct members of the heme–copper oxidase superfamily, were immobilized on electrodes modified with gold nanoparticles. This procedure allowed us to achieve direct electron transfer between the enzyme and the gold nanoparticles and to obtain evidence for different electrocatalytic properties of the two enzymes. The pH dependence and thermostability reveal that the enzymes are highly adapted to their native environments. These results suggest that evolution resulted in different solutions to the common problem of electron transfer to oxygen.

FTIR studies of metal ligands, networks of hydrogen bonds, and water molecules near the active site Mn₄CaO₅ cluster in Photosystem II

The photosynthetic conversion of water to molecular oxygen is catalyzed by the Mn₄CaO₅ cluster in Photosystem II and provides nearly our entire supply of atmospheric oxygen. The Mn₄CaO₅ cluster accumulates oxidizing equivalents in response to light-driven photochemical events within Photosystem II and then oxidizes two molecules of water to oxygen. The Mn₄CaO₅ cluster converts water to oxygen much more efficiently than any synthetic catalyst because its protein environment carefully controls the cluster's reactivity at each step in its catalytic cycle. This control is achieved by precise choreography of the proton and electron transfer reactions associated with water oxidation and by careful management of substrate (water) access and proton egress. This review describes the FTIR studies undertaken over the past two decades to identify the amino acid residues that are responsible for this control and to determine the role of each. In particular, this review describes the FTIR studies undertaken to characterize the influence of the cluster's metal ligands on its activity, to delineate the proton egress pathways that link the Mn₄CaO₅ cluster with the thylakoid lumen, and to characterize the influence of specific residues on the water molecules that serve as substrate or as participants in the networks of hydrogen bonds that make up the water access and proton egress pathways. This information will improve our understanding of water oxidation by the Mn₄CaO₅ catalyst in Photosystem II and will provide insight into the design of new generations of synthetic catalysts that convert sunlight into useful forms of storable energy. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

The causes of reduced proton-pumping efficiency in type B and C respiratory heme-copper oxidases, and in some mutated variants of type A

The heme-copper oxidases may be divided into three categories, A, B, and C, which include cytochrome c and quinol-oxidising enzymes. All three types are known to be proton pumps and are found in prokaryotes, whereas eukaryotes only contain A-type cytochrome c oxidase in their inner mitochondrial membrane. However, the bacterial B- and C-type enzymes have often been reported to pump protons with an H(+)/e(-) ratio of only one half of the unit stoichiometry in the A-type enzyme. We will show here that these observations are likely to be the result of difficulties with the measuring technique together with a higher sensitivity of the B- and C-type enzymes to the protonmotive force that opposes pumping. We find that under optimal conditions the H(+)/e(-) ratio is close to unity in all the three heme-copper oxidase subfamilies. A higher tendency for proton leak in the B- and C-type enzymes may result from less efficient gating of a proton pump mechanism that we suggest evolved before the so-called D-channel of proton transfer. There is also a discrepancy between results using whole bacterial cells vs. phospholipid vesicles inlaid with oxidase with respect to the observed proton pumping after modification of the D-channel residue asparagine-139 (Rhodobacter sphaeroides numbering) to aspartate in A-type cytochrome c oxidase. This discrepancy might also be explained by a higher sensitivity of proton pumping to protonmotive force in the mutated variant. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

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.

Single mutations that redirect internal proton transfer in the ba3 oxidase from Thermus thermophilus

The ba3-type cytochrome c oxidase from Thermus thermophilus is a membrane-bound proton pump. Results from earlier studies have shown that with the aa3-type oxidases proton uptake to the catalytic site and "pump site" occurs simultaneously. However, with ba3 oxidase the pump site is loaded before proton transfer to the catalytic site because the proton transfer to the latter is slower than that with the aa3 oxidases. In addition, the timing of formation and decay of catalytic intermediates is different in the two types of oxidases. In the present study, we have investigated two mutant ba3 CytcOs in which residues of the proton pathway leading to the catalytic site as well as the pump site were exchanged, Thr312Val and Tyr244Phe. Even though ba3 CytcO uses only a single proton pathway for transfer of the substrate and "pumped" protons, the amino-acid residue substitutions had distinctly different effects on the kinetics of proton transfer to the catalytic site and the pump site. The results indicate that the rates of these reactions can be modified independently by replacement of single residues within the proton pathway. Furthermore, the data suggest that the Thr312Val and Tyr244Phe mutations interfere with a structural rearrangement in the proton pathway that is rate limiting for proton transfer to the catalytic site.

Redundancy and modularity in membrane-associated dissimilatory nitrate reduction in Bacillus

The genomes of two phenotypically denitrifying type strains of the genus Bacillus were sequenced and the pathways for dissimilatory nitrate reduction were reconstructed. Results suggest that denitrification proceeds in the periplasmic space and in an analogous fashion as in Gram-negative organisms, yet with the participation of proteins that tend to be membrane-bound or membrane-associated. A considerable degree of functional redundancy was observed with marked differences between B. azotoformans LMG 9581(T) and B. bataviensis LMG 21833(T). In addition to the already characterized menaquinol/cyt c-dependent nitric oxide reductase (Suharti et al., 2001, 2004) of which the encoding genes could be identified now, evidence for another novel nitric oxide reductase (NOR) was found. Also, our analyses confirm earlier findings on branched electron transfer with both menaquinol and cytochrome c as reductants. Quite unexpectedly, both bacilli have the disposal of two parallel pathways for nitrite reduction enabling a life style as a denitrifier and as an ammonifying bacterium.