The mechanism for proton pumping in cytochrome c oxidase from an electrostatic and quantum chemical perspective

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.

Coupled chemical reactions: Effects of electric field, diffusion, and boundary control

Chemical reactions involve the movement of charges, and this paper presents a mathematical model for describing chemical reactions in electrolytes. The model is developed using an energy variational method that aligns with classical thermodynamics principles. It encompasses both electrostatics and chemical reactions within consistently defined energetic and dissipative functionals. Furthermore, the energy variation method is extended to account for open systems that involve the input and output of charge and mass. Such open systems have the capability to convert one form of input energy into another form of output energy. In particular, a two-domain model is developed to study a reaction system with self-regulation and internal switching, which plays a vital role in the electron transport chain of mitochondria responsible for ATP generation-a crucial process for sustaining life. Simulations are conducted to explore the influence of electric potential on reaction rates and switching dynamics within the two-domain system. It shows that the electric potential inhibits the oxidation reaction while accelerating the reduction reaction.

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

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.

Beating Heart of Cytochrome c Oxidase: The Shared Proton of Heme a3 Propionates

Cytochrome c oxidase (CcO) pumps protons from the N-side to the P-side and consumes electrons from the P-side of the mitochondrial membrane driven by energy gained from reduction of dioxygen to water. ATP synthesis uses the resulting proton gradient and electrostatic potential difference. Since the distance a proton travels through CcO is too large for a one-step transfer process, proton-loading sites (PLS) that can carry protons transiently are necessary. One specific pump-active PLS couples to the redox reaction, thus energizing the proton to move across the membrane against electric potential and proton gradient. The PLS should also prevent proton backflow. Therefore, the propionates of the two redox-active hemes in CcO were suggested as PLS candidates although, according to CcO crystal structures, none of the four propionates can be protonated on account of strong H-bonds. Here, we show that modeling the local structure around heme a₃ propionates enhances significantly their capability of carrying a proton jointly. This was not possible for the propionates of heme a. The modeled structures are stable in molecular dynamics simulations (MDS) and are energetically similar to the crystal structure. Precise electrostatic energy computations of MDS data are used to estimate the pK_A values of all titratable residues in CcO. For the modeled structures, the heme a₃ propionates have pK_A values high enough to host a proton transiently but not too high to fix the proton permanently. The change in pK_A throughout the redox reaction is sufficient to push the proton to the P-side of the membrane and to provide the protons with the necessary amount of energy for ATP synthesis.

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

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

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.

Modulation of the electron-proton coupling at cytochrome a by the ligation of the oxidized catalytic center in bovine cytochrome c oxidase

Cytochrome a was suggested as the key redox center in the proton pumping process of bovine cytochrome c oxidase (CcO). Recent studies showed that both the structure of heme a and its immediate vicinity are sensitive to the ligation and the redox state of the distant catalytic center composed of iron of cytochrome a₃ (Fe_a3) and copper (Cu_B). Here, the influence of the ligation at the oxidized Fe_a3³⁺-Cu_B²⁺ center on the electron-proton coupling at heme a was examined in the wide pH range (6.5-11). The strength of the coupling was evaluated by the determination of pH dependence of the midpoint potential of heme a (E_m(a)) for the cyanide (the low-spin Fe_a3³⁺) and the formate-ligated CcO (the high-spin Fe_a3³⁺). The measurements were performed under experimental conditions when other three redox centers of CcO are oxidized. Two slightly differing linear pH dependencies of E_m(a) were found for the CN- and the formate-ligated CcO with slopes of -13 mV/pH unit and -23 mV/pH unit, respectively. These linear dependencies indicate only a weak and unspecific electron-proton coupling at cytochrome a in both forms of CcO. The lack of the strong electron-proton coupling at the physiological pH values is also substantiated by the UV-Vis absorption and electron-paramagnetic resonance spectroscopy investigations of the cyanide-ligated oxidized CcO. It is shown that the ligand exchange at Fe_a³⁺ between His-Fe_a³⁺-His and His-Fe_a³⁺-OH^- occurs only at pH above 9.5 with the estimated pK >11.0.

Explaining leak states in the proton pump of heme-copper oxidases observed in single-molecule experiments

Heme-copper oxidases couple the exergonic oxygen reduction with the endergonic proton translocation. Redox-linked structural changes have been localized in deeply buried regions of the protein, near the low-potential heme. How these movements can modulate distant gating events along the intramolecular proton path, where the entry (exit) of pumped proton occurs, is a major concern for the proton pump models. Generally, these models associate, more or less directly, all translocation events with redox transitions. Although they can account for many phenomenological aspects of the pump, evidences from single-molecules experiments about leak states of the pump represent a formidable challenge. Disconnecting the redox-linked pK_a shifts of the proton loading site from the external barriers, we obtain a simple stochastic mechanism which behaves similarly to the real enzyme, able to reverse the flow of the proton transfer.

Structural changes at the surface of cytochrome c oxidase alter the proton-pumping stoichiometry

Data from earlier studies showed that minor structural changes at the surface of cytochrome c oxidase, in one of the proton-input pathways (the D pathway), result in dramatically decreased activity and a lower proton-pumping stoichiometry. To further investigate how changes around the D pathway orifice influence functionality of the enzyme, here we modified the nearby C-terminal loop of subunit I of the Rhodobacter sphaeroides cytochrome c oxidase. Removal of 16 residues from this flexible surface loop resulted in a decrease in the proton-pumping stoichiometry to <50% of that of the wild-type enzyme. Replacement of the protonatable residue Glu552, part of the same loop, by an Ala, resulted in a similar decrease in the proton-pumping stoichiometry without loss of the O₂-reduction activity or changes in the proton-uptake kinetics. The data show that minor structural changes at the orifice of the D pathway, at a distance of ~40 Å from the proton gate of cytochrome c oxidase, may alter the proton-pumping stoichiometry of the enzyme.

Kinetics and Efficiency of Energy-Transducing Enzymes

Complexes I to IV, with the exception of Complex II, are redox-driven proton pumps that convert redox energy of oxygen reduction to proton gradient across the mitochondrial or bacterial membrane; in turn, the created electrochemical gradient drives the adenosine triphosphate synthesis in the cells by utilizing complex V of the chain. Here we address a general question of the efficiency of such enzymes, considering them as molecular machines that couple endergonic and exergonic reactions and converting one form of free energy into another. One well-known example of the efficiency is given by Carnot's theorem for heat engines. Here we extend the concept to respiratory enzymes and specifically focus on the proton pumping by Complex I of the respiratory chain, nicotinamide adenine dinucleotide dehydrogenase. To discuss the efficiency issues, we develop a model of enzyme kinetics, which generalizes the Michaelis-Menten model. Our model includes several substrates and products and, in general, can be considered as Generalized Michaelis-Menten Kinetic model. The model might be useful for describing complex enzyme kinetics, regardless of the efficiency issues that are addressed in this paper.

Tracing the Pathways of Waters and Protons in Photosystem II and Cytochrome c Oxidase

Photosystem II (PSII) uses water as the terminal electron donor, producing oxygen in the Mn4CaO5 oxygen evolving complex (OEC), while cytochrome c oxidase (CcO) reduces O2 to water in its heme–Cu binuclear center (BNC). Each protein is oriented in the membrane to add to the proton gradient. The OEC, which releases protons, is located near the P-side (positive, at low-pH) of the membrane. In contrast, the BNC is in the middle of CcO, so the protons needed for O2 reduction must be transferred from the N-side (negative, at high pH). In addition, CcO pumps protons from N- to P-side, coupled to the O2 reduction chemistry, to store additional energy. Thus, proton transfers are directly coupled to the OEC and BNC redox chemistry, as well as needed for CcO proton pumping. The simulations that study the changes in proton affinity of the redox active sites and the surrounding protein at different states of the reaction cycle, as well as the changes in hydration that modulate proton transfer paths, are described.

The yeast mitochondrial proteins Rcf1 and Rcf2 support the enzymology of the cytochrome c oxidase complex and generation of the proton motive force

The yeast mitochondrial proteins Rcf1 and Rcf2 are associated with a subpopulation of the cytochrome bc ₁-cytochrome c oxidase supercomplex and have been proposed to play a role in the assembly and/or modulation of the activity of the cytochrome c oxidase (complex IV, CIV). Yeast mutants deficient in either Rcf1 or Rcf2 proteins can use aerobic respiration-based metabolism for growth, but the absence of both proteins results in a strong growth defect. In this study, using assorted biochemical and biophysical analyses of Rcf1/Rcf2 single and double null-mutant yeast cells and mitochondria, we further explored how Rcf1 and Rcf2 support aerobic respiration and growth. We show that the absence of Rcf1 physically reduces the levels of CIV and diminishes the ability of the CIV that is present to maintain a normal mitochondrial proton motive force (PMF). Although the absence of Rcf2 did not noticeably affect the physical content of CIV, the PMF generated by CIV was also lower than normal. Our results indicate that the detrimental effects of the absence of Rcf1 and Rcf2 proteins on the CIV complex are distinct in terms of CIV assembly/accumulation and additive in terms of the ability of CIV to generate PMF. Thus, the combined absence of Rcf1 and Rcf2 alters both CIV physiology and assembly. We conclude that the slow aerobic growth of the Rcf1/Rcf2 double null mutant results from diminished generation of mitochondrial PMF by CIV and limits the level of CIV activity required for maintenance of the PMF and growth under aerobic conditions.

Cytochrome c oxidase structures suggest a four-state stochastic pump mechanism

A simple stochastic model for a cytochrome c oxidase proton pump.

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.

Redox-Driven Proton Pumps of the Respiratory Chain

In aerobic cells, the proton gradient that drives ATP synthesis is created by three different proton pumps-membrane enzymes of the respiratory electron transport chain known as complex I, III, and IV. Despite the striking dissimilarity of structures and apparent differences in molecular mechanisms of proton pumping, all three enzymes have much in common and employ the same universal physical principles of converting redox energy to proton pumping. In this study, we describe a simple mathematical model that illustrates the general principles of redox-driven proton pumps and discuss their implementation in complex I, III, and IV of the respiratory chain.

Insights into functions of the H channel of cytochrome c oxidase from atomistic molecular dynamics simulations

Proton pumping A-type cytochrome c oxidase (CcO) terminates the respiratory chains of mitochondria and many bacteria. Three possible proton transfer pathways (D, K, and H channels) have been identified based on structural, functional, and mutational data. Whereas the D channel provides the route for all pumped protons in bacterial A-type CcOs, studies of bovine mitochondrial CcO have led to suggestions that its H channel instead provides this route. Here, we have studied H-channel function by performing atomistic molecular dynamics simulations on the entire, as well as core, structure of bovine CcO in a lipid-solvent environment. The majority of residues in the H channel do not undergo large conformational fluctuations. Its upper and middle regions have adequate hydration and H-bonding residues to form potential proton-conducting channels, and Asp51 exhibits conformational fluctuations that have been observed crystallographically. In contrast, throughout the simulations, we do not observe transient water networks that could support proton transfer from the N phase toward heme a via neutral His413, regardless of a labile H bond between Ser382 and the hydroxyethylfarnesyl group of heme a In fact, the region around His413 only became sufficiently hydrated when His413 was fixed in its protonated imidazolium state, but its calculated pK_a is too low for this to provide the means to create a proton transfer pathway. Our simulations show that the electric dipole moment of residues around heme a changes with the redox state, hence suggesting that the H channel could play a more general role as a dielectric well.

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.

The mechanism for oxygen reduction in cytochrome c dependent nitric oxide reductase (cNOR) as obtained from a combination of theoretical and experimental results

Bacterial NO-reductases (NOR) belong to the heme-copper oxidase (HCuO) superfamily, in which most members are O₂-reducing, proton-pumping enzymes. This study is one in a series aiming to elucidate the reaction mechanisms of the HCuOs, including the mechanisms for cellular energy conservation. One approach towards this goal is to compare the mechanisms for the different types of HCuOs, cytochrome c oxidase (CcO) and NOR, reducing the two substrates O₂ and NO. Specifically in this study, we describe the mechanism for oxygen reduction in cytochrome c dependent NOR (cNOR). Hybrid density functional calculations were performed on large cluster models of the cNOR binuclear active site. Our results are used, together with published experimental information, to construct a free energy profile for the entire catalytic cycle. Although the overall reaction is quite exergonic, we show that during the reduction of molecular oxygen in cNOR, two of the reduction steps are endergonic with high barriers for proton uptake, which is in contrast to oxygen reduction in CcO, where all reduction steps are exergonic. This difference between the two enzymes is suggested to be important for their differing capabilities for energy conservation. An additional result from this study is that at least three of the four reduction steps are initiated by proton transfer to the active site, which is in contrast to CcO, where electrons always arrive before the protons to the active site. The roles of the non-heme metal ion and the redox-active tyrosine in the active site are also discussed.

Mitochondrial cytochrome c oxidase: catalysis, coupling and controversies

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

Effects of membrane curvature and pH on proton pumping activity of single cytochrome bo3 enzymes

The molecular mechanism of proton pumping by heme-copper oxidases (HCO) has intrigued the scientific community since it was first proposed. We have recently reported a novel technology that enables the continuous characterisation of proton transport activity of a HCO and ubiquinol oxidase from Escherichia coli, cytochrome bo₃, for hundreds of seconds on the single enzyme level (Li et al. J Am Chem Soc 137 (2015) 16055-16063). Here, we have extended these studies by additional experiments and analyses of the proton transfer rate as a function of proteoliposome size and pH at the N- and P-side of single HCOs. Proton transport activity of cytochrome bo₃ was found to decrease with increased curvature of the membrane. Furthermore, proton uptake at the N-side (proton entrance) was insensitive to pH between pH6.4-8.4, while proton release at the P-side had an optimum pH of ~7.4, suggesting that the pH optimum is related to proton release from the proton exit site. Our previous single-enzyme experiments identified rare, long-lived conformation states of cytochrome bo₃ where protons leak back under turn-over conditions. Here, we analyzed and found that ~23% of cytochrome bo₃ proteoliposomes show ΔpH half-lives below 50s after stopping turnover, while only ~5% of the proteoliposomes containing a non-pumping mutant, E286C cytochrome bo₃ exhibit such fast decays. These single-enzyme results confirm our model in which HCO exhibit heterogeneous pumping rates and can adopt rare leak states in which protons are able to rapidly flow back.

A Protonated Water Cluster as a Transient Proton-Loading Site in Cytochrome c Oxidase

Cytochrome c oxidase (CcO) is a redox-driven proton pump that powers aerobic respiratory chains. We show here by multi-scale molecular simulations that a protonated water cluster near the active site is likely to serve as the transient proton-loading site (PLS) that stores a proton during the pumping process. The pKa of this water cluster is sensitive to the redox states of the enzyme, showing distinct similarities to other energy converting proton pumps.

Monte Carlo Simulations of Glu-242 in Cytochrome c Oxidase

Monte Carlo (MC) simulations of conformational changes and protonation of Glu-242, a key residue that shuttles protons in cytochrome c oxidase (CcO), are reported. Previous studies suggest that this residue may play a role of the valve of the enzyme proton pump. Here we examine how sensitive the results of simulations are to the computational method used. We applied both molecular mechanic (MM) and hybrid quantum mechanic:molecular mechanic (QM:MM) methods and find that the results are qualitatively different. The results indicate that the mechanism for proton gating in CcO is still an open issue.

Mechanism of Oxygen Reduction in Cytochrome c Oxidase and the Role of the Active Site Tyrosine

Cytochrome c oxidase, the terminal enzyme in the respiratory chain, reduces molecular oxygen to water and stores the released energy through electrogenic chemistry and proton pumping across the membrane. Apart from the heme-copper binuclear center, there is a conserved tyrosine residue in the active site (BNC). The tyrosine delivers both an electron and a proton during the O-O bond cleavage step, forming a tyrosyl radical. The catalytic cycle then occurs in four reduction steps, each taking up one proton for the chemistry (water formation) and one proton to be pumped. It is here suggested that in three of the reduction steps the chemical proton enters the center of the BNC, leaving the tyrosine unprotonated with radical character. The reproprotonation of the tyrosine occurs first in the final reduction step before binding the next oxygen molecule. It is also suggested that this reduction mechanism and the presence of the tyrosine are essential for the proton pumping. Density functional theory calculations on large cluster models of the active site show that only the intermediates with the proton in the center of the BNC and with an unprotonated tyrosyl radical have a high electron affinity of similar size as the electron donor, which is essential for the ability to take up two protons per electron and thus for the proton pumping. This type of reduction mechanism is also the only one that gives a free energy profile in accordance with experimental observations for the amount of proton pumping in the working enzyme.

Single Enzyme Experiments Reveal a Long-Lifetime Proton Leak State in a Heme-Copper Oxidase

Heme-copper oxidases (HCOs) are key enzymes in prokaryotes and eukaryotes for energy production during aerobic respiration. They catalyze the reduction of the terminal electron acceptor, oxygen, and utilize the Gibbs free energy to transport protons across a membrane to generate a proton (ΔpH) and electrochemical gradient termed proton motive force (PMF), which provides the driving force for the adenosine triphosphate (ATP) synthesis. Excessive PMF is known to limit the turnover of HCOs, but the molecular mechanism of this regulatory feedback remains relatively unexplored. Here we present a single-enzyme study that reveals that cytochrome bo₃ from Escherichia coli, an HCO closely homologous to Complex IV in human mitochondria, can enter a rare, long-lifetime leak state during which proton flow is reversed. The probability of entering the leak state is increased at higher ΔpH. By rapidly dissipating the PMF, we propose that this leak state may enable cytochrome bo₃, and possibly other HCOs, to maintain a suitable ΔpH under extreme redox conditions.

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.

Protonation of the binuclear active site in cytochrome c oxidase decreases the reduction potential of CuB

One of the remaining mysteries regarding the respiratory enzyme cytochrome c oxidase is how proton pumping can occur in all reduction steps in spite of the low reduction potentials observed in equilibrium titration experiments for two of the active site cofactors, CuB(II) and Fea3(III). It has been speculated that, at least the copper cofactor can acquire two different states, one metastable activated state occurring during enzyme turnover, and one relaxed state with lower energy, reached only when the supply of electrons stops. The activated state should have a transiently increased CuB(II) reduction potential, allowing proton pumping. The relaxed state should have a lower reduction potential, as measured in the titration experiments. However, the structures of these two states are not known. Quantum mechanical calculations show that the proton coupled reduction potential for CuB is inherently high in the active site as it appears after reaction with oxygen, which explains the observed proton pumping. It is suggested here that, when the flow of electrons ceases, a relaxed resting state is formed by the uptake of one extra proton, on top of the charge compensating protons delivered in each reduction step. The extra proton in the active site decreases the proton coupled reduction potential for CuB by almost half a volt, leading to agreement with titration experiments. Furthermore, the structure for the resting state with an extra proton is found to have a hydroxo-bridge between CuB(II) and Fea3(III), yielding a magnetic coupling that can explain the experimentally observed EPR silence.

The role of conformational selection in the molecular recognition of the wild type and mutants XPA67-80 peptides by ERCC1

Molecular recognition is a fundamental step in the coordination of biomolecular pathways. Understanding how recognition and binding occur between highly flexible protein domains is a complex task. The conformational selection theory provides an elegant rationalization of the recognition mechanism, especially valid in cases when unstructured protein regions are involved. The recognition of a poorly structured peptide, namely XPA67-80 , by its target receptor ERCC1, falls in this challenging study category. The microsecond molecular dynamics (MD) simulations, discussed in this work, show that the conformational propensity of the wild type XPA67-80 peptide in solution supports conformational selection as the key mechanism driving its molecular recognition by ERCC1. Moreover, all the mutations of the XPA67-80 peptide studied here cause a significant increase of its conformational disorder, relative to the wild type. Comparison to experimental data suggests that the loss of the recognized structural motifs at the microscopic time scale can contribute to the critical decrease in binding observed for one of the mutants, further substantiating the key role of conformational selection in recognition. Ultimately, because of the high sequence identity and analogy in binding, it is conceivable that the conclusions of this study on the XPA67-80 peptide also apply to the ERCC1-binding domain of the XPA protein.

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.

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.

Microscopic basis for kinetic gating in Cytochrome c oxidase: insights from QM/MM analysis

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps.

Understanding the mechanism of vectorial proton pumping in biomolecules requires establishing the microscopic basis for the regulation of both thermodynamic and kinetic features of the relevant proton transfer steps. For the proton pump cytochrome c oxidase, while the regulation of thermodynamic driving force for key proton transfers has been discussed in great detail, the microscopic basis for the control of proton transfer kinetics has been poorly understood. Here we carry out extensive QM/MM free energy simulations to probe the kinetics of relevant proton transfer steps and analyze the effects of local structure and hydration level. We show that protonation of the proton loading site (PLS, taken to be a propionate of heme a₃) requires a concerted process in which a key glutamic acid (Glu286H) delivers the proton to the PLS while being reprotonated by an excess proton coming from the D-channel. The concerted nature of the mechanism is a crucial feature that enables the loading of the PLS before the cavity containing Glu286 is better hydrated to lower its pK_a to experimentally measured range; the charged rather than dipolar nature of the process also ensures a tight coupling with heme a reduction, as emphasized by Siegbahn and Blomberg. In addition, we find that rotational flexibility of the PLS allows its protonation before that of the binuclear center (the site where oxygen gets reduced to water). Together with our recent study (P. Goyal, et al., Proc. Natl. Acad. Sci. U. S. A., 2013, 110, 18886–18891) that focused on the modulation of Glu286 pK_a, the current work suggests a mechanism that builds in a natural sequence for the protonation of the PLS prior to that of the binuclear center. This provides microscopic support to the kinetic constraints revealed by kinetic network analysis as essential elements that ensure an efficient vectorial proton transport in cytochrome c oxidase.

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

After a summary of the problem of coupling electron and proton transfer to proton pumping in cytochrome c oxidase, we present the results of our earlier and recent density functional theory calculations for the dinuclear Fe-a3-CuB reaction center in this enzyme. A specific catalytic reaction wheel diagram is constructed from the calculations, based on the structures and relative energies of the intermediate states of the reaction cycle. A larger family of tautomers/protonation states is generated compared to our earlier work, and a new lowest-energy pathway is proposed. The entire reaction cycle is calculated for the new smaller model (about 185-190 atoms), and two selected arcs of the wheel are chosen for calculations using a larger model (about 205 atoms). We compare the structural and redox energetics and protonation calculations with available experimental data. The reaction cycle map that we have built is positioned for further improvement and testing against experiment.

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.

Conformational determinants for the recruitment of ERCC1 by XPA in the nucleotide excision repair (NER) Pathway: structure and dynamics of the XPA binding motif

XPA is an essential protein in the nucleotide excision repair (NER) pathway, in charge of recruiting the ERCC1-XPF endonuclease complex to the DNA damage site. The only currently available structural insight into the binding of XPA to ERCC1 derives from the solution NMR structure of a complex between the ERCC1 central fragment and a 14-residue peptide, corresponding to the highly conserved binding motif of the XPA N-terminus, XPA₆₇₋₈₀. The extensive all-atom molecular-dynamics simulation study of the XPA₆₇₋₈₀ peptide both bound to the ERCC1 central fragment and free in solution presented here completes the profile of the structural determinants responsible for the ERCC1/XPA₆₇₋₈₀ complex stability. In addition to the wild-type, this study also looks at specific XPA₆₇₋₈₀ mutants in complex with the ERCC1 central domain and thus contributes to defining the conformational determinants for binding, as well as all of the essential structural elements necessary for the rational design of an XPA-based, ERCC1-specific inhibitor.

Proton pumping in cytochrome c oxidase: energetic requirements and the role of two proton channels

Cytochrome c oxidase is a superfamily of membrane bound enzymes catalyzing the exergonic reduction of molecular oxygen to water, producing an electrochemical gradient across the membrane. The gradient is formed both by the electrogenic chemistry, taking electrons and protons from opposite sides of the membrane, and by proton pumping across the entire membrane. In the most efficient subfamily, the A-family of oxidases, one proton is pumped in each reduction step, which is surprising considering the fact that two of the reduction steps most likely are only weakly exergonic. Based on a combination of quantum chemical calculations and experimental information, it is here shown that from both a thermodynamic and a kinetic point of view, it should be possible to pump one proton per electron also with such an uneven distribution of the free energy release over the reduction steps, at least up to half the maximum gradient. A previously suggested pumping mechanism is developed further to suggest a reason for the use of two proton transfer channels in the A-family. Since the rate of proton transfer to the binuclear center through the D-channel is redox dependent, it might become too slow for the steps with low exergonicity. Therefore, a second channel, the K-channel, where the rate is redox-independent is needed. A redox-dependent leakage possibility is also suggested, which might be important for efficient energy conservation at a high gradient. A mechanism for the variation in proton pumping stoichiometry over the different subfamilies of cytochrome oxidase is also suggested. This article is part of a Special Issue entitled: 18th European Bioenergetic Conference.

Mutations in the D-channel of cytochrome c oxidase causes leakage of the proton pump

It has experimentally been found that certain mutations close to the entry point of the proton transfer channel in cytochrome c oxidase stop proton translocation but not the oxygen reduction chemistry. This effect is termed uncoupling. Since the mutations are 20Å away from the catalytic center, this is very surprising. A new explanation for this phenomenon is suggested here, involving a local effect at the entry point of the proton channel, rather than the long range effects suggested earlier.

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.

Redox-state sensing by hydrogen bonds in the CuA center of cytochrome c oxidase

Cytochrome c oxidases (CcO) couple electron transfer to active proton translocation through a gated mechanism that minimizes energy losses by preventing protons from flowing backwards or leaking. Such a complex mechanism requires that information about the redox and protonation states of the different centers be transmitted between different parts of the oxidase. Here we report a network of residues located around the electron entry point of CcO, the CuA site in subunit II, that experience collective pH equilibria around neutral pH. This network starts at the occluded side of the CuA site and extends to the interface between subunits I and II of the CcO, where the proton exit is located and through which electrons flow into subunit I. One of the residues in this network is directly involved in a hydrogen bond to one of the CuA ligands, whose strength is highly sensitive to the redox state of the metal center. We propose that this interaction mediates the transmission of redox changes from ET centers to other functional regions of the oxidase, and possibly also in other similar machineries, as part of their gating and regulatory mechanisms.

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.

Role of aspartate 132 at the orifice of a proton pathway in cytochrome c oxidase

Proton transfer across biological membranes underpins central processes in biological systems, such as energy conservation and transport of ions and molecules. In the membrane proteins involved in these processes, proton transfer takes place through specific pathways connecting the two sides of the membrane via control elements within the protein. It is commonly believed that acidic residues are required near the orifice of such proton pathways to facilitate proton uptake. In cytochrome c oxidase, one such pathway starts near a conserved Asp-132 residue. Results from earlier studies have shown that replacement of Asp-132 by, e.g., Asn, slows proton uptake by a factor of ∼5,000. Here, we show that proton uptake at full speed (∼10(4) s(-1)) can be restored in the Asp-132-Asn oxidase upon introduction of a second structural modification further inside the pathway (Asn-139-Thr) without compensating for the loss of the negative charge. This proton-uptake rate was insensitive to Zn(2+) addition, which in the wild-type cytochrome c oxidase slows the reaction, indicating that Asp-132 is required for Zn(2+) binding. Furthermore, in the absence of Asp-132 and with Thr at position 139, at high pH (>9), proton uptake was significantly accelerated. Thus, the data indicate that Asp-132 is not strictly required for maintaining rapid proton uptake. Furthermore, despite the rapid proton uptake in the Asn-139-Thr/Asp-132-Asn mutant cytochrome c oxidase, proton pumping was impaired, which indicates that the segment around these residues is functionally linked to pumping.

A Kinetic Model of Proton Transport in a Multi-Redox Centre Protein: Cytochrome c Oxidase

We use chemical reaction kinetics to explore the stepwise electron and proton transfer reactions of cytochrome c oxidase (C cO) from R. sphaeroides. Proton transport coupled to electron transport (ET) is investigated in terms of a sequence of protonation-dependent second-order redox reactions. Thereby, we assume fixed rather than shifting dissociation constants of the redox sites. Proton transport can thus be simulated particularly when separate proton uptake and release sites are assumed rather than the same proton pump site for every ET step. In order to test these assumptions, we make use of a model system introduced earlier, which allows us to study direct ET of redox enzymes by electrochemistry. A four-electron transfer model of C cO had been developed before, according to which electrons are transferred from the electrode to CuA. Thereafter, electrons are transferred along the sequence heme a, heme a3 and CuB. In the present investigation, we consider protonation equilibria of the oxidised and reduced species for each of the four centres. Moreover, we add oxygen/H2O as the terminal (fifth) redox couple including protonation of reduced oxygen to water. Finally we arrive at a kinetic model comprising five protonation-dependent redox couples. The results from the simulations are compared with experimental data obtained in the absence and presence of oxygen. As a result, we can show that proton transport can be modelled in terms of protonation-dependent redox kinetics.