Exploring the proton pump and exit pathway for pumped protons in cytochrome ba3 from Thermus thermophilus

Diversity and evolution of nitric oxide reduction in bacteria and archaea

Nitrous oxide is a potent greenhouse gas whose production is catalyzed by nitric oxide reductase (NOR) members of the heme-copper oxidoreductase (HCO) enzyme superfamily. We identified several previously uncharacterized HCO families, four of which (eNOR, sNOR, gNOR, and nNOR) appear to perform NO reduction. These families have novel active-site structures and several have conserved proton channels, suggesting that they might be able to couple NO reduction to energy conservation. We isolated and biochemically characterized a member of the eNOR family from the bacterium Rhodothermus marinus and found that it performs NO reduction. These recently identified NORs exhibited broad phylogenetic and environmental distributions, greatly expanding the diversity of microbes in nature capable of NO reduction. Phylogenetic analyses further demonstrated that NORs evolved multiple times independently from oxygen reductases, supporting the view that complete denitrification evolved after aerobic respiration.

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.

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

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

The Unusual Homodimer of a Heme-Copper Terminal Oxidase Allows Itself to Utilize Two Electron Donors

The heme-copper oxidase superfamily comprises cytochrome c and ubiquinol oxidases. These enzymes catalyze the transfer of electrons from different electron donors onto molecular oxygen. A B-family cytochrome c oxidase from the hyperthermophilic bacterium Aquifex aeolicus was discovered previously to be able to use both cytochrome c and naphthoquinol as electron donors. Its molecular mechanism as well as the evolutionary significance are yet unknown. Here we solved its 3.4 Å resolution electron cryo-microscopic structure and discovered a novel dimeric structure mediated by subunit I (CoxA2) that would be essential for naphthoquinol binding and oxidation. The unique structural features in both proton and oxygen pathways suggest an evolutionary adaptation of this oxidase to its hyperthermophilic environment. Our results add a new conceptual understanding of structural variation of cytochrome c oxidases in different species.

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

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

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

A Water Dimer Shift Activates a Proton Pumping Pathway in the PR → F Transition of ba3 Cytochrome c Oxidase

Broken-symmetry density functional calculations have been performed on the [Fe_a3⁴⁺,Cu_B²⁺] state of the dinuclear center (DNC) for the P_R → F part of the catalytic cycle of ba₃ cytochrome c oxidase (CcO) from Thermus thermophilus (Tt), using the OLYP-D3-BJ functional. The calculations show that the movement of the H₂O molecules in the DNC affects the pK_a values of the residue side chains of Tyr237 and His376⁺, which are crucial for proton transfer/pumping in ba₃ CcO from Tt. The calculated lowest energy structure of the DNC in the [Fe_a3⁴⁺,Cu_B²⁺] state (state F) is of the form Fe_a3⁴⁺═O^2-···Cu_B²⁺, in which the H₂O ligand that resulted from protonation of the OH^- ligand in the P_R state is dissociated from the Cu_B²⁺ site. The calculated Fe_a3⁴⁺═O^2- distance in F (1.68 Å) is 0.03 Å longer than that in P_R (1.65 Å), which can explain the different Fe_a3⁴⁺═O^2- stretching modes in P (804 cm^-1) and F (785 cm^-1) identified by resonance Raman experiments. In this F state, the Cu_B²⁺···O^2- (ferryl-oxygen) distance is only around 2.4 Å. Hence, the subsequent O_H state [Fe_a3³⁺-OH^--Cu_B²⁺] with a μ-hydroxo bridge can be easily formed, as shown by our calculations.

Insights into proton translocation in cbb3 oxidase from MD simulations

Heme-copper oxidases are membrane protein complexes that catalyse the final step of the aerobic respiration, namely the reduction of oxygen to water. The energy released during catalysis is coupled to the active translocation of protons across the membrane, which contributes to the establishment of an electrochemical gradient that is used for ATP synthesis. The distinctive C-type (or cbb₃) cytochrome c oxidases, which are mostly present in proteobacteria, exhibit a number of unique structural and functional features, including high catalytic activity at low oxygen concentrations. At the moment, the functioning mechanism of C-type oxidases, in particular the proton transfer/pumping mechanism presumably via a single proton channel, is still poorly understood. In this work we used all-atom molecular dynamics simulations and continuum electrostatics calculations to obtain atomic-level insights into the hydration and dynamics of a cbb₃ oxidase. We provide the details of the water dynamics and proton transfer pathways for both the "chemical" and "pumped" protons, and show that formation of protonic connections is strongly affected by the protonation state of key residues, namely H243, E323 and H337.

The cellular membrane as a mediator for small molecule interaction with membrane proteins

The cellular membrane constitutes the first element that encounters a wide variety of molecular species to which a cell might be exposed. Hosting a large number of structurally and functionally diverse proteins associated with this key metabolic compartment, the membrane not only directly controls the traffic of various molecules in and out of the cell, it also participates in such diverse and important processes as signal transduction and chemical processing of incoming molecular species. In this article, we present a number of cases where details of interaction of small molecular species such as drugs with the membrane, which are often experimentally inaccessible, have been studied using advanced molecular simulation techniques. We have selected systems in which partitioning of the small molecule with the membrane constitutes a key step for its final biological function, often binding to and interacting with a protein associated with the membrane. These examples demonstrate that membrane partitioning is not only important for the overall distribution of drugs and other small molecules into different compartments of the body, it may also play a key role in determining the efficiency and the mode of interaction of the drug with its target protein. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

ns-μs Time-Resolved Step-Scan FTIR of ba₃ Oxidoreductase from Thermus thermophilus: Protonic Connectivity of w941-w946-w927

Time-resolved step-scan FTIR spectroscopy has been employed to probe the dynamics of the ba₃ oxidoreductase from Thermus thermophilus in the ns-μs time range and in the pH/pD 6–9 range. The data revealed a pH/pD sensitivity of the D372 residue and of the ring-A propionate of heme a₃. Based on the observed transient changes a model in which the protonic connectivity of w941-w946-927 to the D372 and the ring-A propionate of heme a₃ is described.

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

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

Multiscale simulations reveal key features of the proton-pumping mechanism in cytochrome c oxidase

Cytochrome c oxidase (CcO) reduces oxygen to water and uses the released free energy to pump protons across the membrane. We have used multiscale reactive molecular dynamics simulations to explicitly characterize (with free-energy profiles and calculated rates) the internal proton transport events that enable proton pumping during first steps of oxidation of the fully reduced enzyme. Our results show that proton transport from amino acid residue E286 to both the pump loading site (PLS) and to the binuclear center (BNC) are thermodynamically driven by electron transfer from heme a to the BNC, but that the former (i.e., pumping) is kinetically favored whereas the latter (i.e., transfer of the chemical proton) is rate-limiting. The calculated rates agree with experimental measurements. The backflow of the pumped proton from the PLS to E286 and from E286 to the inside of the membrane is prevented by large free-energy barriers for the backflow reactions. Proton transport from E286 to the PLS through the hydrophobic cavity and from D132 to E286 through the D-channel are found to be strongly coupled to dynamical hydration changes in the corresponding pathways and, importantly, vice versa.

Nanosecond ligand migration and functional protein relaxation in ba3 oxidoreductase: Structures of the B0, B1 and B2 intermediate states

Nanosecond time-resolved step-scan FTIR spectroscopy (nTRS (2) -FTIR) has been applied to literally probe the active site of the carbon monoxide (CO)-bound thermophilic ba3 heme-copper oxidoreductase as it executes its function. The nTRS (2) - snapshots of the photolysed heme a3 Fe-CO/CuB species captured a "transition state" whose side chains prevent the photolysed CO to enter the docking cavity. There are three sets of ba3 photoproduct bands of docked CO with different orientation exhibiting different kinetics. The trajectories of the "docked" CO at 2122, 2129 and 2137cm(-1) is referred to in the literature as B2, B1 and B0 intermediate states, respectively. The present data provided direct evidence for the role of water in controlling ligand orientation in an intracavity protein environment.

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

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

Investigating the Proton Donor in the NO Reductase from Paracoccus denitrificans

Variant nomenclature: the variants were made in the NorB subunit if not indicated by the superscript c, which are variants in the NorC subunit (e.g. E122A = exchange of Glu-122 in NorB for an Ala, E71cD; exchange of Glu-71 in NorC for an Asp). Bacterial NO reductases (NORs) are integral membrane proteins from the heme-copper oxidase superfamily. Most heme-copper oxidases are proton-pumping enzymes that reduce O2 as the last step in the respiratory chain. With electrons from cytochrome c, NO reductase (cNOR) from Paracoccus (P.) denitrificans reduces NO to N2O via the following reaction: 2NO+2e-+2H+→N2O+H2O. Although this reaction is as exergonic as O2-reduction, cNOR does not contribute to the electrochemical gradient over the membrane. This means that cNOR does not pump protons and that the protons needed for the reaction are taken from the periplasmic side of the membrane (since the electrons are donated from this side). We previously showed that the P. denitrificans cNOR uses a single defined proton pathway with residues Glu-58 and Lys-54 from the NorC subunit at the entrance. Here we further strengthened the evidence in support of this pathway. Our further aim was to define the continuation of the pathway and the immediate proton donor for the active site. To this end, we investigated the region around the calcium-binding site and both propionates of heme b3 by site directed mutagenesis. Changing single amino acids in these areas often had severe effects on cNOR function, with many variants having a perturbed active site, making detailed analysis of proton transfer properties difficult. Our data does however indicate that the calcium ligation sphere and the region around the heme b3 propionates are important for proton transfer and presumably contain the proton donor. The possible evolutionary link between the area for the immediate donor in cNOR and the proton loading site (PLS) for pumped protons in oxygen-reducing heme-copper oxidases is discussed.

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

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

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

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

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.

Detection of functional hydrogen-bonded water molecules with protonated/deprotonated key carboxyl side chains in the respiratory enzyme ba3-oxidoreductase

In this work we report FTIR data on the detection of functional hydrogen-bonded water molecules with protonated/deprotonated key carboxyl side chains in the respiratory enzyme ba₃-oxidoreductase.

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 transmembrane channel-like protein family and human papillomaviruses: Insights into epidermodysplasia verruciformis and progression to squamous cell carcinoma

Epidermodysplasia verruciformis (EV) is a rare genodermatosis characterized by increased sensitivity to infection by the β-subtype of human papillomaviruses (β-HPVs), causing persistent, tinea versicolor-like dermal lesions. In a majority of affected individuals, these macular lesions progress to invasive cutaneous squamous cell carcinoma (CSCC) in sun-exposed areas. While mutations in transmembrane channel-like 6 (TMC6 / EVER1) and 8 (TMC8 / EVER2) have been causally linked to EV, their molecular functions are unclear. It is likely that their protective effects involve regulation of the β-HPV life cycle, host keratinocyte apoptosis vs. survival balance and/or T-cell interaction with infected host cells.

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.

The K(C) channel in the cbb3-type respiratory oxygen reductase from Rhodobacter capsulatus is required for both chemical and pumped protons

The heme-copper superfamily of proton-pumping respiratory oxygen reductases are classified into three families (A, B, and C families) based on structural and phylogenetic analyses. Most studies have focused on the A family, which includes the eukaryotic mitochondrial cytochrome c oxidase as well as many bacterial homologues. Members of the C family, also called the cbb3-type oxygen reductases, are found only in prokaryotes and are of particular interest because of their presence in a number of human pathogens. All of the heme-copper oxygen reductases require proton-conducting channels to convey chemical protons to the active site for water formation and to convey pumped protons across the membrane. Previous work indicated that there is only one proton-conducting input channel (the K(C) channel) present in the cbb3-type oxygen reductases, which, if correct, must be utilized by both chemical protons and pumped protons. In this work, the effects of mutations in the K(C) channel of the cbb3-type oxygen reductase from Rhodobacter capsulatus were investigated by expressing the mutants in a strain lacking other respiratory oxygen reductases. Proton pumping was evaluated by using intact cells, and catalytic oxygen reductase activity was measured in isolated membranes. Two mutations, N346M and Y374F, severely reduced catalytic activity, presumably by blocking the chemical protons required at the active site. One mutation, T272A, resulted in a substantially lower proton-pumping stoichiometry but did not inhibit oxygen reductase activity. These are the first experimental data in support of the postulate that pumped protons are taken up from the bacterial cytoplasm through the K(C) channel.

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

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

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

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

The nitric-oxide reductase from Paracoccus denitrificans uses a single specific proton pathway

The NO reductase from Paracoccus denitrificans reduces NO to N2O (2NO + 2H(+) + 2e(-) → N2O + H2O) with electrons donated by periplasmic cytochrome c (cytochrome c-dependent NO reductase; cNOR). cNORs are members of the heme-copper oxidase superfamily of integral membrane proteins, comprising the O2-reducing, proton-pumping respiratory enzymes. In contrast, although NO reduction is as exergonic as O2 reduction, there are no protons pumped in cNOR, and in addition, protons needed for NO reduction are derived from the periplasmic solution (no contribution to the electrochemical gradient is made). cNOR thus only needs to transport protons from the periplasm into the active site without the requirement to control the timing of opening and closing (gating) of proton pathways as is needed in a proton pump. Based on the crystal structure of a closely related cNOR and molecular dynamics simulations, several proton transfer pathways were suggested, and in principle, these could all be functional. In this work, we show that residues in one of the suggested pathways (denoted pathway 1) are sensitive to site-directed mutation, whereas residues in the other proposed pathways (pathways 2 and 3) could be exchanged without severe effects on turnover activity with either NO or O2. We further show that electron transfer during single-turnover reduction of O2 is limited by proton transfer and can thus be used to study alterations in proton transfer rates. The exchange of residues along pathway 1 showed specific slowing of this proton-coupled electron transfer as well as changes in its pH dependence. Our results indicate that only pathway 1 is used to transfer protons in cNOR.

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.