Redox-driven membrane-bound proton pumps

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.

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

Cryo-EM structures of intermediates suggest an alternative catalytic reaction cycle for cytochrome c oxidase

Cytochrome c oxidases are among the most important and fundamental enzymes of life. Integrated into membranes they use four electrons from cytochrome c molecules to reduce molecular oxygen (dioxygen) to water. Their catalytic cycle has been considered to start with the oxidized form. Subsequent electron transfers lead to the E-state, the R-state (which binds oxygen), the P-state (with an already split dioxygen bond), the F-state and the O-state again. Here, we determined structures of up to 1.9 Å resolution of these intermediates by single particle cryo-EM. Our results suggest that in the O-state the active site contains a peroxide dianion and in the P-state possibly an intact dioxygen molecule, the F-state may contain a superoxide anion. Thus, the enzyme's catalytic cycle may have to be turned by 180 degrees.

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.

Probing the Proton-Loading Site of Cytochrome C Oxidase Using Time-Resolved Fourier Transform Infrared Spectroscopy

Crystal structure analyses at atomic resolution and FTIR spectroscopic studies of cytochrome c oxidase have yet not revealed protonation or deprotonation of key sites of proton transfer in a time-resolved mode. Here, a sensitive technique to detect protolytic transitions is employed. In this work, probing a proton-loading site of cytochrome c oxidase from Paracoccus denitrificans with time-resolved Fourier transform infrared spectroscopy is presented for the first time. For this purpose, variants with single-site mutations of N131V, D124N, and E278Q, the key residues in the D-channel, were studied. The reaction of mutated CcO enzymes with oxygen was monitored and analyzed. Seven infrared bands in the “fast” kinetic spectra were found based on the following three requirements: (1) they are present in the “fast” phases of N131V and D124N mutants, (2) they have reciprocal counterparts in the “slow” kinetic spectra in these mutants, and (3) they are absent in “fast” kinetic spectra of the E278Q mutant. Moreover, the double-difference spectra between the first two mutants and E278Q revealed more IR bands that may belong to the proton-loading site protolytic transitions. From these results, it is assumed that several polar residues and/or water molecule cluster(s) share a proton as a proton-loading site. This site can be propionate itself (holding only a fraction of H+), His403, and/or water cluster(s).

Comparison Between O and OH Intermediates of Cytochrome c Oxidase Studied by FTIR Spectroscopy

Cytochrome c oxidase is terminal enzyme in the respiratory chain of mitochondria and many aerobic bacteria. It catalyzes reduction of oxygen to water. During its catalysis, CcO proceeds through several quite stable intermediates (R, A, PR/M, O/OH, E/EH). This work is concentrated on the elucidation of the differences between structures of oxidized intermediates O and O H in different CcO variants and at different pH values. Oxidized intermediates of wild type and mutated CcO from Paracoccus denitrificans were studied by means of static and time-resolved Fourier-transform infrared spectroscopy in acidic and alkaline conditions in the infrared region 1800-1000 cm-1. No reasonable differences were found between all variants in these conditions, and in this spectral region. This finding means that the binuclear center of oxygen reduction keeps a very similar structure and holds the same ligands in the studied conditions. The further investigation in search of differences should be performed in the 4000-2000 cm-1 IR region where water ligands absorb.

Integral caa3-Cytochrome c Oxidase from Thermus thermophilus: Purification and Crystallization

Cytochrome c oxidase is a respiratory enzyme catalyzing the energy-conserving reduction of molecular oxygen to water-a fundamental biological process of cell respiration. The first crystal structures of the type A cytochrome c oxidases, bovine heart and Paracoccus denitrificans cytochrome c oxidases, were published in 1995 and contributed immensely to the understanding of the enzyme's mechanism of action. The senior author's research focus was directed toward understanding the structure and function of the type B cytochrome c oxidases, ba₃-oxidase and type A2 caa₃-oxidase, both from the extreme thermophilic bacterium Thermus thermophilus. While the ba₃-oxidase structure was published in 2000 and functional characterization is well-documented in the literature, we recently successfully solved the structure of the caa₃-nature made enzyme-substrate complex. This chapter is dedicated to the purification and crystallization process of caa₃-cytochrome c oxidase.

An update of the chemiosmotic theory as suggested by possible proton currents inside the coupling membrane

Understanding how biological systems convert and store energy is a primary purpose of basic research. However, despite Mitchell's chemiosmotic theory, we are far from the complete description of basic processes such as oxidative phosphorylation (OXPHOS) and photosynthesis. After more than half a century, the chemiosmotic theory may need updating, thanks to the latest structural data on respiratory chain complexes. In particular, up-to date technologies, such as those using fluorescence indicators following proton displacements, have shown that proton translocation is lateral rather than transversal with respect to the coupling membrane. Furthermore, the definition of the physical species involved in the transfer (proton, hydroxonium ion or proton currents) is still an unresolved issue, even though the latest acquisitions support the idea that protonic currents, difficult to measure, are involved. Moreover, F_oF₁-ATP synthase ubiquitous motor enzyme has the peculiarity (unlike most enzymes) of affecting the thermodynamic equilibrium of ATP synthesis. It seems that the concept of diffusion of the proton charge expressed more than two centuries ago by Theodor von Grotthuss is to be taken into consideration to resolve these issues. All these uncertainties remind us that also in biology it is necessary to consider the Heisenberg indeterminacy principle, which sets limits to analytical questions.

The structure of the oxidized state of cytochrome c oxidase - experiments and theory compared

Cytochrome c oxidase (CcO), the terminal enzyme in the respiratory chain, reduces molecular oxygen to water. Experimental data on the midpoint potentials of the heme iron/copper active site cofactors do not match the overall reaction energetics, and are also in conflict with the observed efficiency of energy conservation in CcO. Therefore it has been postulated that the ferric/cupric intermediate (the oxidized state) exists in two forms. One form, labelled O_H, is presumably involved during catalytic turnover, and should have a high Cu_B midpoint potential due to a metastable high energy structure. When no more electrons are supplied, the O_H state supposedly relaxes to the resting form, labelled O, with a lower energy and a lower midpoint potential. It has been suggested that there is a pure geometrical difference between the O_H and O states, obtained by moving a water molecule inside the active site. It is shown here that the difference between the two forms of the oxidized state must be of a more chemical nature. The reason is that all types of geometrically relaxed structures of the oxidized intermediate have similar energies, all with a high proton coupled reduction potential in accordance with the postulated O_H state. One hypothesized chemical modification of the O_H state is the transfer of an extra proton, possibly internal, into the active site. Such a protonated state has several properties that agree with experimental data on the relaxed oxidized state, including a decreased midpoint potential.

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.

Long-range proton-coupled electron transfer in biological energy conversion: towards mechanistic understanding of respiratory complex I

Biological energy conversion is driven by efficient enzymes that capture, store and transfer protons and electrons across large distances. Recent advances in structural biology have provided atomic-scale blueprints of these types of remarkable molecular machinery, which together with biochemical, biophysical and computational experiments allow us to derive detailed energy transduction mechanisms for the first time. Here, I present one of the most intricate and least understood types of biological energy conversion machinery, the respiratory complex I, and how its redox-driven proton-pump catalyses charge transfer across approximately 300 Å distances. After discussing the functional elements of complex I, a putative mechanistic model for its action-at-a-distance effect is presented, and functional parallels are drawn to other redox- and light-driven ion pumps.

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.

The proton pumping bo oxidase from Vitreoscilla

The cytochrome bo3 quinol oxidase from Vitreoscilla (vbo3) catalyses oxidation of ubiquinol and reduction of O2 to H2O. Data from earlier studies suggested that the free energy released in this reaction is used to pump sodium ions instead of protons across a membrane. Here, we have studied the functional properties of heterologously expressed vbo3 with a variety of methods. (i) Following oxygen consumption with a Clark-type electrode, we did not observe a measurable effect of Na+ on the oxidase activity of purified vbo3 solubilized in detergent or reconstituted in liposomes. (ii) Using fluorescent dyes, we find that vbo3 does not pump Na+ ions, but H+ across the membrane, and that H+-pumping is not influenced by the presence of Na+. (iii) Using an oxygen pulse method, it was found that 2 H+/e- are ejected from proteoliposomes, in agreement with the values found for the H+-pumping bo3 oxidase of Escherichia coli (ecbo3). This coincides with the interpretation that 1 H+/e- is pumped across the membrane and 1 H+/e- is released during quinol oxidation. (iv) When the electron transfer kinetics of vbo3 upon reaction with oxygen were followed in single turnover experiments, a similar sequence of reaction steps was observed as reported for the E. coli enzyme and none of these reactions was notably affected by the presence of Na+. Overall the data show that vbo3 is a proton pumping terminal oxidase, behaving similarly to the Escherichia coli bo3 quinol oxidase.

Fermentation Revisited: How Do Microorganisms Survive Under Energy-Limited Conditions?

During fermentation F_OF₁ hydrolyzes ATP, coupling proton transport to proton-motive force (pmf) generation. Despite that, pmf generated by ATP hydrolysis does not satisfy the energy budget of a fermenting cell. However, pmf can also be generated by extrusion of weak organic acids such as lactate and by hydrogen cycling catalyzed by hydrogenases (Hyds). Here we highlight recent advances in our understanding of how the transport of weak organic acids and enzymes contributes to pmf generation during fermentation. The potential impact of these processes on metabolism and energy conservation during microbial fermentation have been overlooked and they not only expand on Mitchell's chemiosmotic theory but also are of relevance to the fields of microbial biochemistry and human and animal health.

A Systematic DFT Approach for Studying Mechanisms of Redox Active Enzymes

When DFT has been applied to study mechanisms of redox processes a common procedure has been to study the results for many different functionals. For redox reactions involving the first row transition metals, this approach has given very different results for different functionals. The conclusion has been that DFT cannot be used for these reactions. In the meantime, results with strong predictability have been generated, most noteworthy for photosystem II, where all DFT predictions have been verified by experiments performed later. In order to obtain these predictive results using DFT, an alternative, systematic approach has been used, where the key differences between the results for different functionals can be rationalized by using a single parameter, rather than using the very large number of differences in the functionals.

pH-sensitive vibrational probe reveals a cytoplasmic protonated cluster in bacteriorhodopsin

Infrared spectroscopy has been used in the past to probe the dynamics of internal proton transfer reactions taking place during the functional mechanism of proteins but has remained mostly silent to protonation changes in the aqueous medium. Here, by selectively monitoring vibrational changes of buffer molecules with a temporal resolution of 6 µs, we have traced proton release and uptake events in the light-driven proton-pump bacteriorhodopsin and correlate these to other molecular processes within the protein. We demonstrate that two distinct chemical entities contribute to the temporal evolution and spectral shape of the continuum band, an unusually broad band extending from 2,300 to well below 1,700 cm^-1 The first contribution corresponds to deprotonation of the proton release complex (PRC), a complex in the extracellular domain of bacteriorhodopsin where an excess proton is shared by a cluster of internal water molecules and/or ionic E194/E204 carboxylic groups. We assign the second component of the continuum band to the proton uptake complex, a cluster with an excess proton reminiscent to the PRC but located in the cytoplasmic domain and possibly stabilized by D38. Our findings refine the current interpretation of the continuum band and call for a reevaluation of the last proton transfer steps in bacteriorhodopsin.

The obligate respiratory supercomplex from Actinobacteria

Actinobacteria are closely linked to human life as industrial producers of bioactive molecules and as human pathogens. Respiratory cytochrome bcc complex and cytochrome aa3 oxidase are key components of their aerobic energy metabolism. They form a supercomplex in the actinobacterial species Corynebacterium glutamicum. With comprehensive bioinformatics and phylogenetic analysis we show that genes for cyt bcc-aa3 supercomplex are characteristic for Actinobacteria (Actinobacteria and Acidimicrobiia, except the anaerobic orders Actinomycetales and Bifidobacteriales). An obligatory supercomplex is likely, due to the lack of genes encoding alternative electron transfer partners such as mono-heme cyt c. Instead, subunit QcrC of bcc complex, here classified as short di-heme cyt c, will provide the exclusive electron transfer link between the complexes as in C. glutamicum. Purified to high homogeneity, the C. glutamicum bcc-aa3 supercomplex contained all subunits and cofactors as analyzed by SDS-PAGE, BN-PAGE, absorption and EPR spectroscopy. Highly uniform supercomplex particles in electron microscopy analysis support a distinct structural composition. The supercomplex possesses a dimeric stoichiometry with a ratio of a-type, b-type and c-type hemes close to 1:1:1. Redox titrations revealed a low potential bcc complex (Em(ISP)=+160mV, Em(bL)=-291mV, Em(bH)=-163mV, Em(cc)=+100mV) fined-tuned for oxidation of menaquinol and a mixed potential aa3 oxidase (Em(CuA)=+150mV, Em(a/a3)=+143/+317mV) mediating between low and high redox potential to accomplish dioxygen reduction. The generated molecular model supports a stable assembled supercomplex with defined architecture which permits energetically efficient coupling of menaquinol oxidation and dioxygen reduction in one supramolecular entity.

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.

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 complexity, challenges and benefits of comparing two transporter classification systems in TCDB and Pfam

Transport systems comprise roughly 10% of all proteins in a cell, playing critical roles in many processes. Improving and expanding their classification is an important goal that can affect studies ranging from comparative genomics to potential drug target searches. It is not surprising that different classification systems for transport proteins have arisen, be it within a specialized database, focused on this functional class of proteins, or as part of a broader classification system for all proteins. Two such databases are the Transporter Classification Database (TCDB) and the Protein family (Pfam) database. As part of a long-term endeavor to improve consistency between the two classification systems, we have compared transporter annotations in the two databases to understand the rationale for differences and to improve both systems. Differences sometimes reflect the fact that one database has a particular transporter family while the other does not. Differing family definitions and hierarchical organizations were reconciled, resulting in recognition of 69 Pfam ‘Domains of Unknown Function’, which proved to be transport protein families to be renamed using TCDB annotations. Of over 400 potential new Pfam families identified from TCDB, 10% have already been added to Pfam, and TCDB has created 60 new entries based on Pfam data. This work, for the first time, reveals the benefits of comprehensive database comparisons and explains the differences between Pfam and TCDB.

How cytochrome c oxidase can pump four protons per oxygen molecule at high electrochemical gradient

Experiments have shown that the A-family cytochrome c oxidases pump four protons per oxygen molecule, also at a high electrochemical gradient. This has been considered a puzzle, since two of the reduction potentials involved, Cu(II) and Fe(III), were estimated from experiments to be too low to afford proton pumping at a high gradient. The present quantum mechanical study (using hybrid density functional theory) suggests a solution to this puzzle. First, the calculations show that the charge compensated Cu(II) potential for CuB is actually much higher than estimated from experiment, of the same order as the reduction potentials for the tyrosyl radical and the ferryl group, which are also involved in the catalytic cycle. The reason for the discrepancy between theory and experiment is the very large uncertainty in the experimental observations used to estimate the equilibrium potentials, mainly caused by the lack of methods for direct determination of reduced CuB. Second, the calculations show that a high energy metastable state, labeled EH, is involved during catalytic turnover. The EH state mixes the low reduction potential of Fe(III) in heme a3 with another, higher potential, here suggested to be that of the tyrosyl radical, resulting in enough exergonicity to allow proton pumping at a high gradient. In contrast, the corresponding metastable oxidized state, OH, is not significantly higher in energy than the resting state, O. Finally, to secure the involvement of the high energy EH state it is suggested that only one proton is taken up via the K-channel during catalytic turnover.

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

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

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.

Current advances in research of cytochrome c oxidase

The function of cytochrome c oxidase as a biomolecular nanomachine that transforms energy of redox reaction into protonmotive force across a biological membrane has been subject of intense research, debate, and controversy. The structure of the enzyme has been solved for several organisms; however details of its molecular mechanism of proton pumping still remain elusive. Particularly, the identity of the proton pumping site, the key element of the mechanism, is still open to dispute. The pumping mechanism has been for a long time one of the key unsolved issues of bioenergetics and biochemistry, but with the accelerating progress in this field many important details and principles have emerged. Current advances in cytochrome oxidase research are reviewed here, along with a brief discussion of the most complete proton pumping mechanism proposed to date, and a molecular basis for control of its efficiency.

The pKa Cooperative: a collaborative effort to advance structure-based calculations of pKa values and electrostatic effects in proteins

The pK(a) Cooperative (http://www.pkacoop.org) was organized to advance development of accurate and useful computational methods for structure-based calculation of pK(a) values and electrostatic energies in proteins. The Cooperative brings together laboratories with expertise and interest in theoretical, computational, and experimental studies of protein electrostatics. To improve structure-based energy calculations, it is necessary to better understand the physical character and molecular determinants of electrostatic effects. Thus, the Cooperative intends to foment experimental research into fundamental aspects of proteins that depend on electrostatic interactions. It will maintain a depository for experimental data useful for critical assessment of methods for structure-based electrostatics calculations. To help guide the development of computational methods, the Cooperative will organize blind prediction exercises. As a first step, computational laboratories were invited to reproduce an unpublished set of experimental pK(a) values of acidic and basic residues introduced in the interior of staphylococcal nuclease by site-directed mutagenesis. The pK(a) values of these groups are unique and challenging to simulate owing to the large magnitude of their shifts relative to normal pK(a) values in water. Many computational methods were tested in this first Blind Prediction Challenge and critical assessment exercise. A workshop was organized in the Telluride Science Research Center to objectively assess the performance of many computational methods tested on this one extensive data set. This volume of Proteins: Structure, Function, and Bioinformatics introduces the pK(a) Cooperative, presents reports submitted by participants in the Blind Prediction Challenge, and highlights some of the problems in structure-based calculations identified during this exercise.

Coupled electron and proton transfer reactions during the O→E transition in bovine cytochrome c oxidase

A combined DFT/electrostatic approach is employed to study the coupling of proton and electron transfer reactions in cytochrome c oxidase (CcO) and its proton pumping mechanism. The coupling of the chemical proton to the internal electron transfer within the binuclear center is examined for the O→E transition. The novel features of the His291 pumping model are proposed, which involve timely well-synchronized sequence of the proton-coupled electron transfer reactions. The obtained pK(a)s and E(m)s of the key ionizable and redox-active groups at the different stages of the O→E transition are consistent with available experimental data. The PT step from E242 to H291 is examined in detail for various redox states of the hemes and various conformations of E242 side-chain. Redox potential calculations of the successive steps in the reaction cycle during the O→E transition are able to explain a cascade of equilibria between the different intermediate states and electron redistribution between the metal centers during the course of the catalytic activity. All four electrometric phases are discussed in the light of the obtained results, providing a robust support for the His291 model of proton pumping in CcO.

Poly(aniline) nanowires in sol-gel coated ITO: a pH-responsive substrate for planar supported lipid bilayers

Facilitated ion transport across an artificial lipid bilayer coupled to a solid substrate is a function common to several types of bioelectronic devices based on supported membranes, including biomimetic fuel cells and ion channel biosensors. Described here is fabrication of a pH-sensitive transducer composed of a porous sol-gel layer derivatized with poly(aniline) (PANI) nanowires grown from an underlying planar indium-tin oxide (ITO) electrode. The upper sol-gel surface is hydrophilic, smooth, and compatible with deposition of a planar supported lipid bilayer (PSLB) formed via vesicle fusion. Conducting tip AFM was used to show that the PANI wires are connected to the ITO, which convert this electrode into a potentiometric pH sensor. The response to changes in the pH of the buffer contacting the PANI nanowire/sol-gel/ITO electrode is blocked by the very low ion permeability of the overlying fluid PSLB. The feasibility of using this assembly to monitor facilitated proton transport across the PSLB was demonstrated by doping the membrane with lipophilic ionophores that respond to a transmembrane pH gradient, which produced an apparent proton permeability several orders of magnitude greater than values measured for undoped lipid bilayers.

Stabilization of the peroxy intermediate in the oxygen splitting reaction of cytochrome cbb(3)

The proton-pumping cbb(3)-type cytochrome c oxidases catalyze cell respiration in many pathogenic bacteria. For reasons not yet understood, the apparent dioxygen (O(2)) affinity in these enzymes is very high relative to other members of the heme-copper oxidase (HCO) superfamily. Based on density functional theory (DFT) calculations on intermediates of the oxygen scission reaction in active-site models of cbb(3)- and aa(3)-type oxidases, we find that a transient peroxy intermediate (I(P), Fe[III]-OOH(-)) is ~6kcal/mol more stable in the former case, resulting in more efficient kinetic trapping of dioxygen and hence in a higher apparent oxygen affinity. The major molecular basis for this stabilization is a glutamate residue, polarizing the proximal histidine ligand of heme b(3) in the active site.

Redox Bohr effects and the role of heme a in the proton pump of bovine heart cytochrome c oxidase

Structural and functional observations are reviewed which provide evidence for a central role of redox Bohr effect linked to the low-spin heme a in the proton pump of bovine heart cytochrome c oxidase. Data on the membrane sidedness of Bohr protons linked to anaerobic oxido-reduction of the individual metal centers in the liposome reconstituted oxidase are analysed. Redox Bohr protons coupled to anaerobic oxido-reduction of heme a (and Cu(A)) and Cu(B) exhibit membrane vectoriality, i.e. protons are taken up from the inner space upon reduction of these centers and released in the outer space upon their oxidation. Redox Bohr protons coupled to anaerobic oxido-reduction of heme a(3) do not, on the contrary, exhibit vectorial nature: protons are exchanged only with the outer space. A model of the proton pump of the oxidase, in which redox Bohr protons linked to the low-spin heme a play a central role, is described. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

B3LYP study on reduction mechanisms from O2 to H2O at the catalytic sites of fully reduced and mixed-valence bovine cytochrome c oxidases

Reduction mechanisms of oxygen molecule to water molecules in the fully reduced (FR) and mixed-valence (MV) bovine cytochrome c oxidases (CcO) have been systematically examined based on the B3LYP calculations. The catalytic cycle using four electrons and four protons has been also shown consistently. The MV CcO catalyses reduction to produce one water molecule, while the FR CcO catalyses to produce two water molecules. One water molecule is added into vacant space between His240 and His290 in the catalytic site. This water molecule constructs the network of hydrogen bonds of Tyr244, farnesyl ethyl, and Thr316 that is a terminal residue of the K-pathway. It plays crucial roles for the proton transfer to the dioxygen to produce the water molecules in both MV and FR CcOs. Tyr244 functions as a relay of the proton transfer from the K-pathway to the added water molecule, not as donors of a proton and an electron to the dioxygen. The reduction mechanisms of MV and FR CcOs are strictly distinguished. In the FR CcO, the Cu atom at the Cu(B) site maintains the reduced state Cu(I) during the process of formation of first water molecule and plays an electron storage. At the final stage of formation of first water molecule, the Cu(I) atom releases an electron to Fe-O. During the process of formation of second water molecule, the Cu atom maintains the oxidized state Cu(II). In contrast with experimental proposals, the K-pathway functions for formation of first water molecule, while the D-pathway functions for second water molecule. The intermediates, P(M), P(R), F, and O, obtained in this work are compared with those proposed experimentally.

Microscopic pKa analysis of Glu286 in cytochrome c oxidase (Rhodobacter sphaeroides): toward a calibrated molecular model

As stringent tests for the molecular model and computational protocol, microscopic pK(a) calculations are performed for the key residue, Glu286, in cytochrome c oxidase (CcO) using a combined quantum mechanical/molecular mechanical (QM/MM) potential and a thermodynamic integration protocol. The impact of the number of water molecules in the hydrophobic cavity and protonation state of several key residues (e.g., His334, Cu(B)-bound water, and PRD(a3)) on the computed microscopic pK(a) values of Glu286 has been systematically examined. To help evaluate the systematic errors in the QM/MM-based protocol, microscopic pK(a) calculations have also been carried out for sites in a soluble protein (Asp70 in T4 lysozyme) and a better-characterized membrane protein (Asp85 in bacteriorhodopsin). Overall, the results show a significant degree of internal consistency and reproducibility that support the effectiveness of the computational framework. Although the number of water molecules in the hydrophobic cavity does not greatly influence the computed pK(a) of Glu286, the protonation states of several residues, some of which are rather far away, have more significant impacts. Adopting the standard protonation state for all titratable residues leaves a large net charge on the system and a significantly elevated pK(a) for Glu286, highlighting that any attempt to address the energetics of proton transfers in CcO at a microscopic level should carefully select the protonation state of residues, even those not in the immediate neighborhood of the active site. The calculations indirectly argue against the deprotonation of His334 for the proton pumping process, although further studies that explicitly compute its pK(a) are required for a more conclusive statement. Finally, the deprotonated Glu286 is found to be in a stable water-mediated connection with PRD(a3) for at least several nanoseconds when this presumed pumping site is protonated. This does not support the proposed role of Glu286 as a robust gating valve that prevents proton leakage, although a conclusive statement awaits a more elaborate characterization of the Glu286-PRD(a3) connectivity with free energy simulations and a protonated PRD(a3). The large sets of microscopic simulations performed here have provided useful guidance to the establishment of a meaningful molecular model and effective computational protocol for explicitly analyzing the proton transfer kinetics in CcO, which is required for answering key questions regarding the pumping function of this fascinating and complex system.

Role of copper ion in regulating ligand binding in a myoglobin-based cytochrome C oxidase model

Cytochrome c oxidase (CcO), the terminal enzyme in the mitochondrial respiratory chain, catalyzes the four-electron reduction of dioxygen to water in a binuclear center comprised of a high-spin heme (heme a(3)) and a copper atom (Cu(B)) coordinated by three histidine residues. As a minimum model for CcO, a mutant of sperm whale myoglobin, named Cu(B)Mb, has been engineered, in which a copper atom is held in the distal heme pocket by the native E7 histidine and two nonnative histidine residues. In this work, the role of the copper in regulating ligand binding in Cu(B)Mb was investigated. Resonance Raman studies show that the presence of copper in CO-bound Cu(B)Mb leads to a CcO-like distal heme pocket. Stopped-flow data show that, upon the initiation of the CO binding reaction, the ligand first binds to the Cu(+); it subsequently transfers from Cu(+) to Fe(2+) in an intramolecular process, similar to that reported for CcO. The high CO affinity toward Cu(+) and the slow intramolecular CO transfer rate between Cu(+) and Fe(2+) in the Cu(B)Mb/Cu(+) complex are analogous to those in Thermus thermophilus CcO (TtCcO) but distinct from those in bovine CcO (bCcO). Additional kinetic studies show that, upon photolysis of the NO-bound Cu(B)Mb/Cu(+) complex, the photolyzed ligand transiently binds to Cu(+) and subsequently rebinds to Fe(2+), accounting for the 100% geminate recombination yield, similar to that found in TtCcO. The data demonstrate that the Cu(B)Mb/Cu(+) complex reproduces essential structural and kinetic features of CcO and that the complex is more akin to TtCcO than to bCcO.

Critical structural role of R481 in cytochrome c oxidase from Rhodobacter sphaeroides

The R481 residue in cytochrome c oxidase from Rhodobacter sphaeroides forms hydrogen bonds with the propionate groups of both heme a and heme a(3). It has been postulated that R481 is the proton loading site in the proton exit pathway essential for proton translocation. A recent functional study showed that the mutations of R481 to His, Leu and Gln cause the reduction of the activity to approximately 5-18% of the native level, and the absence of proton pumping in R481Q but retention of approximately 40% efficiency in R481H and R481L (H.J. Lee, L. Ojemyr, A. Vakkasoglu, P. Brzezinski and R. B. Gennis, manuscript submitted). To decipher the molecular mechanism underlying the perturbed functionalities, we have used resonance Raman spectroscopy to examine the structural properties of the three mutants. The data show that the frequencies of the formyl CO stretching modes of both the heme a and a(3) in the mutants are characteristic of formyl groups exposed to an aqueous environment, indicating that the mutations disrupt the native H-bonding interaction between the formyl group of heme a and R52, as well as the hydrophobic environment surrounding the formyl group of heme a(3). In addition to the change in the environments of heme a and a(3), the Raman data show that the mutations induce a partial conversion of the heme a(3) from a high-spin to a low-spin state, suggesting that the mutations are associated with the rearrangement of the Cu(B)-heme a(3) binuclear center. The Raman results reported here demonstrate that R481 plays a critical role in supporting efficient proton pumping, by holding the heme groups in a proper environment.

Properties of Arg481 mutants of the aa3-type cytochrome c oxidase from Rhodobacter sphaeroides suggest that neither R481 nor the nearby D-propionate of heme a3 is likely to be the proton loading site of the proton pump

Cytochrome c oxidase utilizes the energy from electron transfer and reduction of oxygen to water and pumps protons across the membrane, generating a proton motive force. A large body of biochemical work has shown that all the pumped protons enter the enzyme through the D-channel, which is apparent in X-ray structures as a chain of water molecules connecting D132 at the cytoplasmic surface of the enzyme to E286, near the enzyme active site. The exit pathway utilized by pumped protons beyond this point and leading to the bacterial periplasm is not known. Also not known is the proton loading site (or sites) which undergoes changes in pKa in response to the chemistry at the enzyme active site and drives the proton pump mechanism. In this paper, we examine the role of R481, a highly conserved arginine that forms an ion pair with the D-propionate of heme a3. The R481H, R481N, R481Q, and R481L mutants were examined. The R481H mutant oxidase is approximately 18% active and pumps protons with approximately 40% of the stoichiometry of the wild type. The R481N, R481Q, and R481L mutants each retain only approximately 5% of the steady-state activity, and this is shown to be due to inhibition of steps in the reaction of O(2) with the reduced enzyme. Neither the R481N mutant nor the R481Q mutant oxidases pump protons, but remarkably, the R481L mutant does pump protons with the same efficiency as the R481H mutant. Since the proton pump is clearly operating in the R481L mutant, these results rule out an essential role in the proton pump mechanism for R481 or its hydrogen bond partner, the D-propionate of heme a3.

Kinetic gating of the proton pump in cytochrome c oxidase

Cytochrome c oxidase (CcO), the terminal enzyme of the respiratory chain, reduces oxygen to water and uses the released energy to pump protons across a membrane. Here, we use kinetic master equations to explore the energetic and kinetic control of proton pumping in CcO. We construct models consistent with thermodynamic principles, the structure of CcO, experimentally known proton affinities, and equilibrium constants of intermediate reactions. The resulting models are found to capture key properties of CcO, including the midpoint redox potentials of the metal centers and the electron transfer rates. We find that coarse-grained models with two proton sites and one electron site can pump one proton per electron against membrane potentials exceeding 100 mV. The high pumping efficiency of these models requires strong electrostatic couplings between the proton loading (pump) site and the electron site (heme a), and kinetic gating of the internal proton transfer. Gating is achieved by enhancing the rate of proton transfer from the conserved Glu-242 to the pump site on reduction of heme a, consistent with the predictions of the water-gated model of proton pumping. The model also accounts for the phenotype of D-channel mutations associated with loss of pumping but retained turnover. The fundamental mechanism identified here for the efficient conversion of chemical energy into an electrochemical potential should prove relevant also for other molecular machines and novel fuel-cell designs.

A combined picture from theory and experiments on water oxidation, oxygen reduction and proton pumping

In order to illustrate how theory and experiments can be combined, examples are taken from two projects that have been going on for a decade. The goal is to obtain the full mechanistic picture of water oxidation in photosystem II and proton pumping in cytochrome c oxidase. It is argued that for obtaining a complete description of these processes, both experiments and theoretical calculations are needed. It is obvious that there are aspects, which are out of reach for computations, but there are also key aspects that can not be obtained by experiments. This concerns very short-lived species but also, in the case of photosynthesis in particular, structural information that is presently out of reach. The main contributions from theory in the present cases, is for photosynthesis a mechanism for O-O bond formation including new and improved structures for all S-states, and for proton pumping a plausible and simple mechanism for proton gating. The examples also illustrate that sometimes rather qualitative experimental information can be of highest importance.