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

Time-resolved serial crystallography to track the dynamics of carbon monoxide in the active site of cytochrome c oxidase

Cytochrome c oxidase (CcO) is part of the respiratory chain and contributes to the electrochemical membrane gradient in mitochondria as well as in many bacteria, as it uses the energy released in the reduction of oxygen to pump protons across an energy-transducing biological membrane. Here, we use time-resolved serial femtosecond crystallography to study the structural response of the active site upon flash photolysis of carbon monoxide (CO) from the reduced heme a3 of ba3-type CcO. In contrast with the aa3-type enzyme, our data show how CO is stabilized on CuB through interactions with a transiently ordered water molecule. These results offer a structural explanation for the extended lifetime of the CuB-CO complex in ba3-type CcO and, by extension, the extremely high oxygen affinity of the enzyme.

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.

A simple goniometer-compatible flow cell for serial synchrotron X-ray crystallography

Serial femtosecond crystallography was initially developed for room-temperature X-ray diffraction studies of macromolecules at X-ray free electron lasers. When combined with tools that initiate biological reactions within microcrystals, time-resolved serial crystallography allows the study of structural changes that occur during an enzyme catalytic reaction. Serial synchrotron X-ray crystallography (SSX), which extends serial crystallography methods to synchrotron radiation sources, is expanding the scientific community using serial diffraction methods. This report presents a simple flow cell that can be used to deliver microcrystals across an X-ray beam during SSX studies. This device consists of an X-ray transparent glass capillary mounted on a goniometer-compatible 3D-printed support and is connected to a syringe pump via light-weight tubing. This flow cell is easily mounted and aligned, and it is disposable so can be rapidly replaced when blocked. This system was demonstrated by collecting SSX data at MAX IV Laboratory from microcrystals of the integral membrane protein cytochrome c oxidase from Thermus thermophilus, from which an X-ray structure was determined to 2.12 Å resolution. This simple SSX platform may help to lower entry barriers for non-expert users of SSX.

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

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

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.

Oxygen Activation and Energy Conservation by Cytochrome c Oxidase

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

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

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

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

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

Serial femtosecond crystallography structure of cytochrome c oxidase at room temperature

Cytochrome c oxidase catalyses the reduction of molecular oxygen to water while the energy released in this process is used to pump protons across a biological membrane. Although an extremely well-studied biological system, the molecular mechanism of proton pumping by cytochrome c oxidase is still not understood. Here we report a method to produce large quantities of highly diffracting microcrystals of ba₃-type cytochrome c oxidase from Thermus thermophilus suitable for serial femtosecond crystallography. The room-temperature structure of cytochrome c oxidase is solved to 2.3 Å resolution from data collected at an X-ray Free Electron Laser. We find overall agreement with earlier X-ray structures solved from diffraction data collected at cryogenic temperature. Previous structures solved from synchrotron radiation data, however, have shown conflicting results regarding the identity of the active-site ligand. Our room-temperature structure, which is free from the effects of radiation damage, reveals that a single-oxygen species in the form of a water molecule or hydroxide ion is bound in the active site. Structural differences between the ba₃-type and aa₃-type cytochrome c oxidases around the proton-loading site are also described.

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.

Reaction of S. cerevisiae mitochondria with ligands: Kinetics of CO and O2 binding to flavohemoglobin and cytochrome c oxidase

Kinetic methods used to investigate electron and proton transfer within cytochrome c oxidase (CytcO) are often based on the use of light to dissociate small ligands, such as CO, thereby initiating the reaction. Studies of intact mitochondria using these methods require identification of proteins that may bind CO and determination of the ligand-binding kinetics. In the present study we have investigated the kinetics of CO-ligand binding to S. cerevisiae mitochondria and cellular extracts. The data indicate that CO binds to two proteins, CytcO and a (yeast) flavohemoglobin (yHb). The latter has been shown previously to reside in both the cell cytosol and the mitochondrial matrix. Here, we found that yHb resides also in the intermembrane space and binds CO in its reduced state. As observed previously, we found that the yHb population in the mitochondrial matrix binds CO, but only after removal of the inner membrane. The mitochondrial yHb (in both the intermembrane space and the matrix) recombines with CO with τ≅270ms, which is significantly slower than observed with the cytosolic yHb (main component τ≅1.3ms). The data indicate that the yHb populations in the different cell compartments differ in structure.

Rapid Electron Transfer within the III-IV Supercomplex in Corynebacterium glutamicum

Complex III in C. glutamicum has an unusual di-heme cyt. c₁ and it co-purifies with complex IV in a supercomplex. Here, we investigated the kinetics of electron transfer within this supercomplex and in the cyt. aa₃ alone (cyt. bc₁ was removed genetically). In the reaction of the reduced cyt. aa₃ with O₂, we identified the same sequence of events as with other A-type oxidases. However, even though this reaction is associated with proton uptake, no pH dependence was observed in the kinetics. For the cyt. bc₁-cyt. aa₃ supercomplex, we observed that electrons from the c-hemes were transferred to Cu_A with time constants 0.1–1 ms. The b-hemes were oxidized with a time constant of 6.5 ms, indicating that this electron transfer is rate-limiting for the overall quinol oxidation/O₂ reduction activity (~210 e⁻/s). Furthermore, electron transfer from externally added cyt. c to cyt. aa₃ was significantly faster upon removal of cyt. bc₁ from the supercomplex, suggesting that one of the c-hemes occupies a position near Cu_A. In conclusion, isolation of the III-IV-supercomplex allowed us to investigate the kinetics of electron transfer from the b-hemes, via the di-heme cyt. c₁ and heme a to the heme a₃-Cu_B catalytic site of cyt. aa₃.

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.

Isolation of yeast complex IV in native lipid nanodiscs

We used the amphipathic styrene maleic acid (SMA) co-polymer to extract cytochrome c oxidase (CytcO) in its native lipid environment from S. cerevisiae mitochondria. Native nanodiscs containing one CytcO per disc were purified using affinity chromatography. The longest cross-sections of the native nanodiscs were 11nm×14nm. Based on this size we estimated that each CytcO was surrounded by ~100 phospholipids. The native nanodiscs contained the same major phospholipids as those found in the mitochondrial inner membrane. Even though CytcO forms a supercomplex with cytochrome bc₁ in the mitochondrial membrane, cyt. bc₁ was not found in the native nanodiscs. Yet, the loosely-bound Respiratory SuperComplex factors were found to associate with the isolated CytcO. The native nanodiscs displayed an O₂-reduction activity of ~130 electrons CytcO^-1s^-1 and the kinetics of the reaction of the fully reduced CytcO with O₂ was essentially the same as that observed with CytcO in mitochondrial membranes. The kinetics of CO-ligand binding to the CytcO catalytic site was similar in the native nanodiscs and the mitochondrial membranes. We also found that excess SMA reversibly inhibited the catalytic activity of the mitochondrial CytcO, presumably by interfering with cyt. c binding. These data point to the importance of removing excess SMA after extraction of the membrane protein. Taken together, our data shows the high potential of using SMA-extracted CytcO for functional and structural studies.

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.

Regulatory role of the respiratory supercomplex factors in Saccharomyces cerevisiae

The respiratory supercomplex factors (Rcf) 1 and 2 mediate supramolecular interactions between mitochondrial complexes III (ubiquinol-cytochrome c reductase; cyt. bc1) and IV (cytochrome c oxidase; CytcO). In addition, removal of these polypeptides results in decreased activity of CytcO, but not of cyt. bc1 In the present study, we have investigated the kinetics of ligand binding, the single-turnover reaction of CytcO with O2, and the linked cyt. bc1-CytcO quinol oxidation-oxygen-reduction activities in mitochondria in which Rcf1 or Rcf2 were removed genetically (strains rcf1Δ and rcf2Δ, respectively). The data show that in the rcf1Δ and rcf2Δ strains, in a significant fraction of the population, ligand binding occurs over a time scale that is ∼100-fold faster (τ ≅ 100 μs) than observed with the wild-type mitochondria (τ ≅ 10 ms), indicating structural changes. This effect is specific to removal of Rcf and not dissociation of the cyt. bc1-CytcO supercomplex. Furthermore, in the rcf1Δ and rcf2Δ strains, the single-turnover reaction of CytcO with O2 was incomplete. This observation indicates that the lower activity of CytcO is caused by a fraction of inactive CytcO rather than decreased CytcO activity of the entire population. Furthermore, the data suggest that the Rcf1 polypeptide mediates formation of an electron-transfer bridge from cyt. bc1 to CytcO via a tightly bound cyt. c We discuss the significance of the proposed regulatory mechanism of Rcf1 and Rcf2 in the context of supramolecular interactions between cyt. bc1 and CytcO.

Analyzing the electrogenicity of cytochrome c oxidase

Measurements of voltage changes in response to charge separation within membrane proteins can offer fundamental information on spectroscopically "invisible" steps. For example, results from studies of voltage changes associated with electron and proton transfer in cytochrome c oxidase could, in principle, be used to discriminate between different theoretical models describing the molecular mechanism of proton pumping. Earlier analyses of data from these measurements have been based on macroscopic considerations that may not allow for exploring the actual molecular mechanisms. Here, we have used a coarse-grained model describing the relation between observed voltage changes and specific charge-transfer reactions, which includes an explicit description of the membrane, the electrolytes, and the electrodes. The results from these calculations offer mechanistic insights at the molecular level. Our main conclusion is that previously assumed mechanistic evidence that was based on electrogenic measurements is not unique. However, the ability of our calculations to obtain reliable voltage changes means that we have a tool that can be used to describe a wide range of electrogenic charge transfers in channels and transporters, by combining voltage measurements with other experiments and simulations to analyze new mechanistic proposals.

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.