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

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

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

Structural insights into functional properties of the oxidized form of cytochrome c oxidase

Cytochrome c oxidase (CcO) is an essential enzyme in mitochondrial and bacterial respiration. It catalyzes the four-electron reduction of molecular oxygen to water and harnesses the chemical energy to translocate four protons across biological membranes. The turnover of the CcO reaction involves an oxidative phase, in which the reduced enzyme (R) is oxidized to the metastable OH state, and a reductive phase, in which OH is reduced back to the R state. During each phase, two protons are translocated across the membrane. However, if OH is allowed to relax to the resting oxidized state (O), a redox equivalent to OH, its subsequent reduction to R is incapable of driving proton translocation. Here, with resonance Raman spectroscopy and serial femtosecond X-ray crystallography (SFX), we show that the heme a3 iron and CuB in the active site of the O state, like those in the OH state, are coordinated by a hydroxide ion and a water molecule, respectively. However, Y244, critical for the oxygen reduction chemistry, is in the neutral protonated form, which distinguishes O from OH, where Y244 is in the deprotonated tyrosinate form. These structural characteristics of O provide insights into the proton translocation mechanism of CcO.

Charge Regulation in a Rieske Proton Pump Pinpoints Zero, One, and Two Proton-Coupled Electron Transfer

The degree to which redox-driven proton pumps regulate net charge during electron transfer (ΔZ_ET) remains undetermined due to difficulties in measuring the net charge of solvated proteins. Values of ΔZ_ET can reflect reorganization energies or redox potentials associated with ET and can be used to distinguish ET from proton(s)-coupled electron transfer (PCET). Here, we synthesized protein "charge ladders" of a Rieske [2Fe-2S] subunit from Thermus thermophilus (truncTtRp) and made 120 electrostatic measurements of ΔZ_ET across pH. Across pH 5-10, truncTtRp is suspected of transitioning from ET to PCET, and then to two proton-coupled ET (2PCET). Upon reduction, we found that truncTtRp became more negative at pH 6.0 by one unit (ΔZ_ET = -1.01 ± 0.14), consistent with single ET; was isoelectric at pH 8.8 (ΔZ_ET = -0.01 ± 0.45), consistent with PCET; and became more positive at pH 10.6 (ΔZ_ET = +1.37 ± 0.60), consistent with 2PCET. These ΔZ_ET values are attributed to protonation of H154 and H134. Across pH, redox potentials of TtRp (measured previously) correlated with protonation energies of H154 and H134 and ΔZ_ET for truncTtRp, supporting a discrete proton pumping mechanism for Rieske proteins at the Fe-coordinating histidines.

Structural basis for functional properties of cytochrome c oxidase

Cytochrome c oxidase (CcO) is an essential enzyme in mitochondrial and bacterial respiration. It catalyzes the four-electron reduction of molecular oxygen to water and harnesses the chemical energy to translocate four protons across biological membranes, thereby establishing the proton gradient required for ATP synthesis1. The full turnover of the CcO reaction involves an oxidative phase, in which the reduced enzyme (R) is oxidized by molecular oxygen to the metastable oxidized OH state, and a reductive phase, in which OH is reduced back to the R state. During each of the two phases, two protons are translocated across the membranes2. However, if OH is allowed to relax to the resting oxidized state (O), a redox equivalent to OH, its subsequent reduction to R is incapable of driving proton translocation2,3. How the O state structurally differs from OH remains an enigma in modern bioenergetics. Here, with resonance Raman spectroscopy and serial femtosecond X-ray crystallography (SFX)4, we show that the heme a3 iron and CuB in the active site of the O state, like those in the OH state5,6, are coordinated by a hydroxide ion and a water molecule, respectively. However, Y244, a residue covalently linked to one of the three CuB ligands and critical for the oxygen reduction chemistry, is in the neutral protonated form, which distinguishes O from OH, where Y244 is in the deprotonated tyrosinate form. These structural characteristics of O provide new insights into the proton translocation mechanism of CcO.

Electric fields control water-gated proton transfer in cytochrome c oxidase

Aerobic life is powered by membrane-bound enzymes that catalyze the transfer of electrons to oxygen and protons across a biological membrane. Cytochrome c oxidase (CcO) functions as a terminal electron acceptor in mitochondrial and bacterial respiratory chains, driving cellular respiration and transducing the free energy from O₂ reduction into proton pumping. Here we show that CcO creates orientated electric fields around a nonpolar cavity next to the active site, establishing a molecular switch that directs the protons along distinct pathways. By combining large-scale quantum chemical density functional theory (DFT) calculations with hybrid quantum mechanics/molecular mechanics (QM/MM) simulations and atomistic molecular dynamics (MD) explorations, we find that reduction of the electron donor, heme a, leads to dissociation of an arginine (Arg438)-heme a₃ D-propionate ion-pair. This ion-pair dissociation creates a strong electric field of up to 1 V Å^-1 along a water-mediated proton array leading to a transient proton loading site (PLS) near the active site. Protonation of the PLS triggers the reduction of the active site, which in turn aligns the electric field vectors along a second, "chemical," proton pathway. We find a linear energy relationship of the proton transfer barrier with the electric field strength that explains the effectivity of the gating process. Our mechanism shows distinct similarities to principles also found in other energy-converting enzymes, suggesting that orientated electric fields generally control enzyme catalysis.

Second coordination sphere effects on the mechanistic pathways for dioxygen activation by a ferritin: involvement of a Tyr radical and the identification of a cation binding site

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp65, which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center.

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.

Theoretical Studies of Acetyl-CoA Synthase Catalytic Mechanism

DFT calculations were performed for the A-cluster from the enzyme Acetyl-CoA synthase (ACS). The acid constants (pKa), reduction potentials, and pH-dependent reduction potential for the A-cluster with different oxidation states and ligands were calculated. Good agreement of the reduction potentials, dependent on pH in the experiment, was obtained. On the basis of the calculations, a mechanism for the methylation reaction involving two–electron reduction and protonation on the proximal nickel atom of the reduced A-cluster is proposed.

Multidimensional redox potential/pKa coupling in multicopper oxidases from molecular dynamics: implications for the proton transfer mechanism

Bilirubin oxidases (BOD) are metalloenzymes that catalyze the conversion of O₂ and bilirubin to biliverdin and water in the metabolism of chlorophyll and porphyrin. In this work we have used the CpHMD method to analyze the effects of the different oxidation states on the BOD trinuclear cluster (TNC). Our results demonstrate that there is a link between the different oxidation states of copper ions and the protonation capacity of nearby titratable residues. Each configuration affects pK_a differently, creating proton gradients within the enzyme that act in an extremely orderly manner. This order is closely linked to the catalytic mechanism and leads us to the conclusion of the entry of the O₂ molecule and its reduction in water molecules is associated with the probability of the release of protons from nearby acid groups. With this information, we deduce that under the initial reaction conditions the acidic side chains of nearby residues can be protonated; this allows the enzyme to reduce the activation energy of the reaction by coupling the proton transfer to oxidation state changes in the metallic center.

Protonation of the oxo-bridged heme/copper assemblies: Modeling the oxidized state of the cytochrome c oxidase active site

In this study on model compounds for the resting oxidized state of the iron‑copper binuclear center in cytochrome c oxidase (CcO), we describe the synthesis of a new μ-oxo-heme/Cu complex, [(TPP)Fe^III-O-Cu^II(tmpa)][B(C₆F₅)₄] (2) {TPP: tetraphenyl porphyrinate(2-); TMPA: tris(2-pyridylmethylamine)}, as well as two protonation events for three μ-oxo-heme/Cu complexes with varying peripheral substituents on the heme site. The addition of increasing amounts of strong acid to these μ-oxo-heme/Cu systems successively led to the generation of the corresponding μ-hydroxo, μ-aquo, and the dissociated complexes. The heme/Cu assemblies bridged through a water ligand are reported here for the first time and the ¹H NMR and ¹⁹F NMR spectral properties are consistent with antiferromagnetically coupled high-spin iron(III) and copper(II) centers. By titration using a series of protonated amines, the pK_a values for the corresponding μ-hydroxo-heme/Cu species (i.e., the first protonation event) have been reported and compared with the pK_a ranges previously estimated for related systems. These synthetic systems may represent structural models for the oxidized Fe^III-X-Cu^II resting state, or turnover intermediates and can be employed to clarify the nature of proton/electron transfer events in CcO. SYNOPSIS: The resting oxidized state of the cytochrome c oxidase active site contains an Fe_a3-OH_x-Cu_B moiety. Here, we investigated two successive protonation events, for a series of μ-oxo-heme/Cu assemblies and reported the pK_a values for the first protonation event. The μ-aquo-heme/Cu complexes described here are the first examples of such systems.

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.