Redox-coupled crystal structural changes in bovine heart cytochrome c oxidase

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.

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.

Energetic Mechanism of Cytochrome c-Cytochrome c Oxidase Electron Transfer Complex Formation under Turnover Conditions Revealed by Mutational Effects and Docking Simulation

Based on the mutational effects on the steady-state kinetics of the electron transfer reaction and our NMR analysis of the interaction site (Sakamoto, K., Kamiya, M., Imai, M., Shinzawa-Itoh, K., Uchida, T., Kawano, K., Yoshikawa, S., and Ishimori, K. (2011) Proc. Natl. Acad. Sci. U.S.A. 108, 12271-12276), we determined the structure of the electron transfer complex between cytochrome c (Cyt c) and cytochrome c oxidase (CcO) under turnover conditions and energetically characterized the interactions essential for complex formation. The complex structures predicted by the protein docking simulation were computationally selected and validated by the experimental kinetic data for mutant Cyt c in the electron transfer reaction to CcO. The interaction analysis using the selected Cyt c-CcO complex structure revealed the electrostatic and hydrophobic contributions of each amino acid residue to the free energy required for complex formation. Several charged residues showed large unfavorable (desolvation) electrostatic interactions that were almost cancelled out by large favorable (Columbic) electrostatic interactions but resulted in the destabilization of the complex. The residual destabilizing free energy is compensated by the van der Waals interactions mediated by hydrophobic amino acid residues to give the stabilized complex. Thus, hydrophobic interactions are the primary factors that promote complex formation between Cyt c and CcO under turnover conditions, whereas the change in the electrostatic destabilization free energy provides the variance of the binding free energy in the mutants. The distribution of favorable and unfavorable electrostatic interactions in the interaction site determines the orientation of the binding of Cyt c on CcO.

How a novel tyrosine-heme cross-link fine-tunes the structure and functions of heme proteins: a direct comparitive study of L29H/F43Y myoglobin

A heme-protein cross-link is a key post-translational modification (PTM) of heme proteins. Meanwhile, the structural and functional consequences of heme-protein cross-links are not fully understood, due to limited studies on a direct comparison of the same protein with and without the cross-link. A Tyr-heme cross-link with a C-O bond is a newly discovered PTM of heme proteins, and is spontaneously formed in F43Y myoglobin (Mb) between the Tyr hydroxyl group and the heme 4-vinyl group in vivo. In this study, we found that with an additional distal His29 introduced in the heme pocket, the double mutant L29H/F43Y Mb can form two distinct forms under different protein purification conditions, with and without a novel Tyr-heme cross-link. By solving the X-ray structures of both forms of L29H/F43Y Mb and performing spectroscopic studies, we made a direct structural and functional comparison in the same protein scaffold. It revealed that the Tyr-heme cross-link regulates the heme distal hydrogen-bonding network, and fine-tunes not only the spectroscopic and ligand binding properties, but also the protein reactivity. Moreover, the formation of the Tyr-heme cross-link in the double mutant L29H/F43Y Mb was investigated in vitro. This study addressed the key issue of how Tyr-heme cross-link fine-tunes the structure and functions of the heme protein, and provided a plausible mechanism for the formation of the newly discovered Tyr-heme cross-link.

SecA: a potential antimicrobial target

There is a consensus in the medical profession of the pressing need for novel antimicrobial agents due to issues related to drug resistance. In practice, solutions to this problem to a large degree lie with the identification of new and vital targets in bacteria and subsequently designing their inhibitors. We consider SecA a very promising antimicrobial target. In this review, we compile and analyze information available on SecA to show that inhibition of SecA has a multitude of consequences. Furthermore, we discuss issues critical to the design and evaluation of SecA inhibitors.

Calcium ions inhibit reduction of heme a in bovine cytochrome c oxidase

The effect of Ca(2+) on the rate of heme a reduction by dithionite and hexaammineruthenium (RuAm) was studied in the cyanide-complexed bovine cytochrome oxidase (CcO). The rate of heme a reduction is proportional to RuAm concentration below 300 μM with kv of 0.53×10(6) M(-1) s(-1). Ca(2+) inhibits the rate of heme a reduction by dithionite by ∼25%. As the reaction speeds up with increased concentrations of RuAm, the inhibition by Ca(2+) disappears. The inhibition of heme a reduction may contribute to recently described partial inhibition of CcO by Ca(2+) in the enzymatic assays. The inhibitory effect of Ca(2+) on heme a reduction indicates that ET through heme a may be coupled to proton movement in the exit part of the proton channel H.

ATP generation in a host cell in early-phase infection is increased by upregulation of cytochrome c oxidase activity via the p2 peptide from human immunodeficiency virus type 1 Gag

Background

Human immunodeficiency virus type 1 (HIV-1) must take advantage of its own proteins with two or more functions to successfully replicate. Although many attempts have been made to determine the function of viral proteins encoded in the HIV-1 genome, the role of the p2 peptide, a spacer between the capsid and the nucleocapsid in HIV-1 Gag in early-phase HIV infection still remains unclarified.

Results

In this study, we show that the p2 peptide enhances HIV-1 acute infection by increasing intracellular ATP production via the activation of mitochondrial cytochrome c oxidase (MT-CO) involved in the respiratory chain. We found that cell-permeable p2-peptide-treated cells were more effectively infected by HIV-1 than control cells. To characterize the effect of the p2 peptide on HIV-1 replication in MAGIC-5 cells, various HIV-1 cDNA products were measured by quantitative real-time PCR. The levels of the late (R/gag), 2-LTR circular (2-LTR), and integrated (Alu) forms of viral cDNAs increased in the presence of the p2 peptide. Interestingly, yeast two-hybrid analysis revealed a novel interaction between the p2 peptide and the mitochondrial intermembrane space domain (N²¹⁴–F²³⁵) of MT-CO subunit I (MT-CO1). Mutational analysis indicated that Gln⁶ in the p2 peptide is important for the interaction with MT-CO1. The p2 peptide activated MT-CO1 in vitro in a concentration-dependent manner, and fluorescence-microscopy analysis demonstrated that the p2 peptide had a significant effect on mitochondrial targeting. Furthermore, the analysis of HIV-1 lacking a functional p2 peptide demonstrated the inhibition of intracellular ATP production in MT-4 cells and monocyte-derived macrophages (MDMs) and a decrease in reverse transcription efficiency following infection of MT-4 cells and MDMs.

Conclusions

These findings provide evidence that the p2 peptide is a viral positive allosteric modulator of MT-CO and the increased intracellular ATP production after HIV infection in a p2-peptide-dependent manner is essential for efficient reverse transcription in early-phase HIV-1 infection.

Loop recognition and copper-mediated disulfide reduction underpin metal site assembly of CuA in human cytochrome oxidase

Maturation of cytochrome oxidases is a complex process requiring assembly of several subunits and adequate uptake of the metal cofactors. Two orthologous Sco proteins (Sco1 and Sco2) are essential for the correct assembly of the dicopper CuA site in the human oxidase, but their function is not fully understood. Here, we report an in vitro biochemical study that shows that Sco1 is a metallochaperone that selectively transfers Cu(I) ions based on loop recognition, whereas Sco2 is a copper-dependent thiol reductase of the cysteine ligands in the oxidase. Copper binding to Sco2 is essential to elicit its redox function and as a guardian of the reduced state of its own cysteine residues in the oxidizing environment of the mitochondrial intermembrane space (IMS). These results provide a detailed molecular mechanism for CuA assembly, suggesting that copper and redox homeostasis are intimately linked in the mitochondrion.

Roles of mitochondrial energy dissipation systems in plant development and acclimation to stress

BACKGROUND

Plants are sessile organisms that have the ability to integrate external cues into metabolic and developmental signals. The cues initiate specific signal cascades that can enhance the tolerance of plants to stress, and these mechanisms are crucial to the survival and fitness of plants. The adaption of plants to stresses is a complex process that involves decoding stress inputs as energy-deficiency signals. The process functions through vast metabolic and/or transcriptional reprogramming to re-establish the cellular energy balance. Members of the mitochondrial energy dissipation pathway (MEDP), alternative oxidases (AOXs) and uncoupling proteins (UCPs), act as energy mediators and might play crucial roles in the adaption of plants to stresses. However, their roles in plant growth and development have been relatively less explored.

SCOPE

This review summarizes current knowledge about the role of members of the MEDP in plant development as well as recent advances in identifying molecular components that regulate the expression of AOXs and UCPs. Highlighted in particular is a comparative analysis of the expression, regulation and stress responses between AOXs and UCPs when plants are exposed to stresses, and a possible signal cross-talk that orchestrates the MEDP, reactive oxygen species (ROS), calcium signalling and hormone signalling.

CONCLUSIONS

The MEDP might act as a cellular energy/metabolic mediator that integrates ROS signalling, energy signalling and hormone signalling with plant development and stress accumulation. However, the regulation of MEDP members is complex and occurs at transcriptional, translational, post-translational and metabolic levels. How this regulation is linked to actual fluxes through the AOX/UCP in vivo remains elusive.

Heterogeneous natural selection on oxidative phosphorylation genes among fishes with extreme high and low aerobic performance

Background

Oxidative phosphorylation (OXPHOS) is the primary source of ATP in eukaryotes and serves as a mechanistic link between variation in genotypes and energetic phenotypes. While several physiological and anatomical factors may lead to increased aerobic capacity, variation in OXPHOS proteins may influence OXPHOS efficiency and facilitate adaptation in organisms with varied energy demands. Although there is evidence that natural selection acts on OXPHOS genes, the focus has been on detection of directional (positive) selection on specific phylogenetic branches where traits that increase energetic demands appear to have evolved. We examined patterns of selection in a broader evolutionary context, i.e., on multiple lineages of fishes with extreme high and low aerobic performance.

Results

We found that patterns of natural selection on mitochondrial OXPHOS genes are complex among fishes with different swimming performance. Positive selection is not consistently associated with high performance taxa and appears to be strongest on lineages containing low performance taxa. In contrast, within high performance lineages, purifying (negative) selection appears to predominate.

Conclusions

We provide evidence that selection on OXPHOS varies in both form and intensity within and among lineages through evolutionary time. These results provide evidence for fluctuating selection on OXPHOS associated with divergence in aerobic performance. However, in contrast to previous studies, positive selection was strongest on low performance taxa suggesting that adaptation of OXPHOS involves many factors beyond enhancing ATP production in high performance taxa. The broader pattern indicates a complex interplay between organismal adaptations, ATP demand, and OXPHOS function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0453-7) contains supplementary material, which is available to authorized users.

Three toxic gases meet in the mitochondria

The rationale of the study was two-fold: (i) develop a functional synthetic model of the Cytochrome c oxidase (CcO) active site, (ii) use it as a convenient tool to understand or predict the outcome of the reaction of CcO with ligands (physiologically relevant gases and other ligands). At physiological pH and potential, the model catalyzes the 4-electron reduction of oxygen. This model was immobilized on self-assembled-monolayer (SAM) modified electrode. During catalytic oxygen reduction, electron delivery through SAMs is rate limiting, similar to the situation in CcO. This model contains all three redox-active components in CcO's active site, which are required to minimize the production of partially-reduced-oxygen-species (PROS): Fe-heme (“heme a3”) in a myoglobin-like model fitted with a proximal imidazole ligand, and a distal tris-imidazole Copper (“Cu_B”) complex, where one imidazole is cross-linked to a phenol (mimicking “Tyr244”). This functional CcO model demonstrates how CcO itself might tolerate the hormone NO (which diffuses through the mitochondria). It is proposed that Cu_B delivers superoxide to NO bound to Fe-heme forming peroxynitrite, then nitrate that diffuses away. Another toxic gas, H₂S, has exceptional biological effects: at ~80 ppm, H₂S induces a state similar to hibernation in mice, lowering the animal's temperature and slowing respiration. Using our functional CcO model, we have demonstrated that at the same concentration range H₂S can reversibly inhibit catalytic oxygen reduction. Such a reversible catalytic process on the model was also demonstrated with an organic compound, tetrazole (TZ). Following studies showed that TZ reversibly inhibits respiration in isolated mitochondria, and induces deactivation of platelets, a mitochondria-rich key component of blood coagulation. Hence, this program is a rare example illustrating the use of a functional model to understand and predict physiologically important reactions at the active site of CcO.

Synthetic heme/copper assemblies: toward an understanding of cytochrome c oxidase interactions with dioxygen and nitrogen oxides

Our long-time niche in synthetic biological inorganic chemistry has been to design ligands and generate coordination complexes of copper or iron ions or both, those reacting with dioxygen (O2) or nitrogen oxides (e.g., nitric oxide (NO(g)) and nitrite (NO2(-))) or both. As inspiration for this work, we turn to mitochondrial cytochrome c oxidase, which is responsible for dioxygen consumption and is also the predominant target for NO(g) and nitrite within mitochondria. In this Account, we highlight recent advances in studying synthetic heme/Cu complexes in two respects. First, there is the design, synthesis, and characterization of new O2 adducts whose further study will add insights into O2 reductive cleavage chemistry. Second, we describe how related heme/Cu constructs reduce nitrite ion to NO(g) or the reverse, oxidize NO(g) to nitrite. The reactions of nitrogen oxides occur as part of CcO's function, which is intimately tied to cellular O2 balance. We had first discovered that reduced heme/Cu compounds react with O2 giving μ-oxo heme-Fe(III)-O-Cu(II)(L) products; their properties are discussed. The O-atom is derived from dioxygen, and interrogations of these systems led to the construction and characterization of three distinctive classes of heme-peroxo complexes, two high-spin and one low-spin species. Recent investigations include a new approach to the synthesis of low-spin heme-peroxo-Cu complexes, employing a "naked" synthon, where the copper ligand denticity and geometric types can be varied. The result is a collection of such complexes; spectroscopic and structural features (by DFT calculations) are described. Some of these compounds are reactive toward reductants/protons effecting subsequent O-O cleavage. This points to how subtle improvements in ligand environment lead to a desired local structure and resulting optimized reactivity, as known to occur at enzyme active sites. The other sector of research is focused on heme/Cu assemblies mediating the redox interplay between nitrite and NO(g). In the nitrite reductase chemistry, the cupric center serves as a Lewis acid, while the heme is the redox active center providing the electron. The orientation of nitrite in approaching the ferrous heme center and N-atom binding are important. Also, detailed spectroscopic and kinetic studies of the NO(g) oxidase chemistry, in excellent agreement with theoretical calculations, reveal the intermediates and key mechanistic steps. Thus, we suggest that both chemical and biochemical heme/Cu-mediated nitrite reductase and NO(g) oxidase chemistry require N-atom binding to a ferrous heme along with cupric ion O-atom coordination, proceeding via a three-membered O-Fe-N chelate ring transition state. These important mechanistic features of heme/Cu systems interconverting NO(g) and nitrite are discussed for the first time.

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.

The subunit composition and function of mammalian cytochrome c oxidase

Cytochrome c oxidase (COX) from mammals and birds is composed of 13 subunits. The three catalytic subunits I-III are encoded by mitochondrial DNA, the ten nuclear-coded subunits (IV, Va, Vb, VIa, VIb, VIc, VIIa, VIIb, VIIc, VIII) by nuclear DNA. The nuclear-coded subunits are essentially involved in the regulation of oxygen consumption and proton translocation by COX, since their removal or modification changes the activity and their mutation causes mitochondrial diseases. Respiration, the basis for ATP synthesis in mitochondria, is differently regulated in organs and species by expression of tissue-, developmental-, and species-specific isoforms for COX subunits IV, VIa, VIb, VIIa, VIIb, and VIII, but the holoenzyme in mammals is always composed of 13 subunits. Various proteins and enzymes were shown, e.g., by co-immunoprecipitation, to bind to specific COX subunits and modify its activity, but these interactions are reversible, in contrast to the tightly bound 13 subunits. In addition, the formation of supercomplexes with other oxidative phosphorylation complexes has been shown to be largely variable. The regulatory complexity of COX is increased by protein phosphorylation. Up to now 18 phosphorylation sites have been identified under in vivo conditions in mammals. However, only for a few phosphorylation sites and four nuclear-coded subunits could a specific function be identified. Research on the signaling pathways leading to specific COX phosphorylations remains a great challenge for understanding the regulation of respiration and ATP synthesis in mammalian organisms. This article reviews the function of the individual COX subunits and their isoforms, as well as proteins and small molecules interacting and regulating the enzyme.

Respiratory conservation of energy with dioxygen: cytochrome C oxidase

Cytochrome c oxidase (CcO) is the terminal oxidase of cell respiration which reduces molecular oxygen (O₂) to H2O coupled with the proton pump. For elucidation of the mechanism of CcO, the three-dimensional location and chemical reactivity of each atom composing the functional sites have been extensively studied by various techniques, such as crystallography, vibrational and time-resolved electronic spectroscopy, since the X-ray structures (2.8 Å resolution) of bovine and bacterial CcO have been published in 1995.X-ray structures of bovine CcO in different oxidation and ligand binding states showed that the O₂reduction site, which is composed of Fe (heme a 3) and Cu (CuB), drives a non-sequential four-electron transfer for reduction of O₂to water without releasing any reactive oxygen species. These data provide the crucial structural basis to solve a long-standing problem, the mechanism of the O₂reduction.Time-resolved resonance Raman and charge translocation analyses revealed the mechanism for coupling between O₂reduction and the proton pump: O₂is received by the O₂reduction site where both metals are in the reduced state (R-intermediate), giving the O₂-bound form (A-intermediate). This is spontaneously converted to the P-intermediate, with the bound O₂fully reduced to 2 O²⁻. Hereafter the P-intermediate receives four electron equivalents from the second Fe site (heme a), one at a time, to form the three intermediates, F, O, and E to regenerate the R-intermediate. Each electron transfer step from heme a to the O₂reduction site is coupled with the proton pump.X-ray structural and mutational analyses of bovine CcO show three possible proton transfer pathways which can transfer pump protons (H) and chemical (water-forming) protons (K and D). The structure of the H-pathway of bovine CcO indicates that the driving force of the proton pump is the electrostatic repulsion between the protons on the H-pathway and positive charges of heme a, created upon oxidation to donate electrons to the O₂reduction site. On the other hand, mutational and time-resolved electrometric findings for the bacterial CcO strongly suggest that the D-pathway transfers both pump and chemical protons. However, the structure for the proton-gating system in the D-pathway has not been experimentally identified. The structural and functional diversities in CcO from various species suggest a basic proton pumping mechanism in which heme a pumps protons while heme a 3 reduces O₂as proposed in 1978.

Mitochondrial Morphology and Fundamental Parameters of the Mitochondrial Respiratory Chain Are Altered in Caenorhabditis elegans Strains Deficient in Mitochondrial Dynamics and Homeostasis Processes

Mitochondrial dysfunction has been linked to myriad human diseases and toxicant exposures, highlighting the need for assays capable of rapidly assessing mitochondrial health in vivo. Here, using the Seahorse XFe24 Analyzer and the pharmacological inhibitors dicyclohexylcarbodiimide and oligomycin (ATP-synthase inhibitors), carbonyl cyanide 4-(trifluoromethoxy) phenylhydrazone (mitochondrial uncoupler) and sodium azide (cytochrome c oxidase inhibitor), we measured the fundamental parameters of mitochondrial respiratory chain function: basal oxygen consumption, ATP-linked respiration, maximal respiratory capacity, spare respiratory capacity and proton leak in the model organism Caenhorhabditis elegans. Since mutations in mitochondrial homeostasis genes cause mitochondrial dysfunction and have been linked to human disease, we measured mitochondrial respiratory function in mitochondrial fission (drp-1)-, fusion (fzo-1)-, mitophagy (pdr-1, pink-1)-, and electron transport chain complex III (isp-1)-deficient C. elegans. All showed altered function, but the nature of the alterations varied between the tested strains. We report increased basal oxygen consumption in drp-1; reduced maximal respiration in drp-1, fzo-1, and isp-1; reduced spare respiratory capacity in drp-1 and fzo-1; reduced proton leak in fzo-1 and isp-1; and increased proton leak in pink-1 nematodes. As mitochondrial morphology can play a role in mitochondrial energetics, we also quantified the mitochondrial aspect ratio for each mutant strain using a novel method, and for the first time report increased aspect ratios in pdr-1- and pink-1-deficient nematodes.

X-ray structure of cyanide-bound bovine heart cytochrome c oxidase in the fully oxidized state at 2.0 Å resolution

The X-ray structure of cyanide-bound bovine heart cytochrome c oxidase in the fully oxidized state was determined at 2.0 Å resolution. The structure reveals that the peroxide that bridges the two metals in the fully oxidized state is replaced by a cyanide ion bound in a nearly symmetric end-on fashion without significantly changing the protein conformation outside the two metal sites.

Nitrogen Oxide Atom-Transfer Redox Chemistry; Mechanism of NO(g) to Nitrite Conversion Utilizing μ-oxo Heme-Fe(III)-O-Cu(II)(L) Constructs

While nitric oxide (NO, nitrogen monoxide) is a critically important signaling agent, its cellular concentrations must be tightly controlled, generally through its oxidative conversion to nitrite (NO2(-)) where it is held in reserve to be reconverted as needed. In part, this reaction is mediated by the binuclear heme a3/CuB active site of cytochrome c oxidase. In this report, the oxidation of NO(g) to nitrite is shown to occur efficiently in new synthetic μ-oxo heme-Fe(III)-O-Cu(II)(L) constructs (L being a tridentate or tetradentate pyridyl/alkylamino ligand), and spectroscopic and kinetic investigations provide detailed mechanistic insights. Two new X-ray structures of μ-oxo complexes have been determined and compared to literature analogs. All μ-oxo complexes react with 2 mol equiv NO(g) to give 1:1 mixtures of discrete [(L)Cu(II)(NO2(-))](+) plus ferrous heme-nitrosyl compounds; when the first NO(g) equiv reduces the heme center and itself is oxidized to nitrite, the second equiv of NO(g) traps the ferrous heme thus formed. For one μ-oxo heme-Fe(III)-O-Cu(II)(L) compound, the reaction with NO(g) reveals an intermediate species ("intermediate"), formally a bis-NO adduct, [(NO)(porphyrinate)Fe(II)-(NO2(-))-Cu(II)(L)](+) (λmax = 433 nm), confirmed by cryo-spray ionization mass spectrometry and EPR spectroscopy, along with the observation that cooling a 1:1 mixture of [(L)Cu(II)(NO2(-))](+) and heme-Fe(II)(NO) to -125 °C leads to association and generation of the key 433 nm UV-vis feature. Kinetic-thermodynamic parameters obtained from low-temperature stopped-flow measurements are in excellent agreement with DFT calculations carried out which describe the sequential addition of NO(g) to the μ-oxo complex.

Defining the role of tyrosine and rational tuning of oxidase activity by genetic incorporation of unnatural tyrosine analogs

While a conserved tyrosine (Tyr) is found in oxidases, the roles of phenol ring pKa and reduction potential in O2 reduction have not been defined despite many years of research on numerous oxidases and their models. These issues represent major challenges in our understanding of O2 reduction mechanism in bioenergetics. Through genetic incorporation of unnatural amino acid analogs of Tyr, with progressively decreasing pKa of the phenol ring and increasing reduction potential, in the active site of a functional model of oxidase in myoglobin, a linear dependence of both the O2 reduction activity and the fraction of H2O formation with the pKa of the phenol ring has been established. By using these unnatural amino acids as spectroscopic probe, we have provided conclusive evidence for the location of a Tyr radical generated during reaction with H2O2, by the distinctive hyperfine splitting patterns of the halogenated tyrosines and one of its deuterated derivatives incorporated at the 33 position of the protein. These results demonstrate for the first time that enhancing the proton donation ability of the Tyr enhances the oxidase activity, allowing the Tyr analogs to augment enzymatic activity beyond that of natural Tyr.

Significant Improvement of Oxidase Activity through the Genetic Incorporation of a Redox-active Unnatural Amino Acid

Incorporation of 3-methoxytyrosine boosts the oxidase activity of the myoglobin model of oxidase, stressing the importance of the redox potential tuning of tyrosine.

While nature employs various covalent and non-covalent strategies to modulate tyrosine (Y) redox potential and pK_a in order to optimize enzyme activities, such approaches have not been systematically applied for the design of functional metalloproteins. Through the genetic incorporation of 3-methoxytyrosine (OMeY) into myoglobin, we replicated important features of cytochrome c oxidase (CcO) in this small soluble protein, which exhibits selective O₂ reduction activity while generating a small amount of reactive oxygen species (ROS). These results demonstrate that the electron donating ability of a tyrosine residue in the active site is important for CcO function. Moreover, we elucidated the structural basis for the genetic incorporation of OMeY into proteins by solving the X-ray structure of OMeY specific aminoacyl-tRNA synthetase complexed with OMeY.

CO-dynamics in the active site of cytochrome c oxidase

The transfer of CO from heme a3 to the Cu(B) site in Cytochrome c oxidase (CcO) after photolysis is studied using molecular dynamics simulations using an explicitly reactive, parametrized potential energy surface based on density functional theory calculations. After photodissociation from the heme-Fe, the CO ligand rebinds to the Cu(B) site on the sub-picosecond time scale. Depending on the simulation protocol the characteristic time ranges from 260 fs to 380 fs which compares with an estimated 450 fs from experiment based on the analysis of the spectral changes as a function of time delay after the photodissociating pulse. Following photoexcitation ≈90% of the ligands are found to rebind to either the Cu(B) (major component, 85%) or the heme-Fe (minor component, 2%) whereas about 10% remain in an unbound state. The infrared spectra of unbound CO in the active site is broad and featureless and no appreciable shift relative to gas-phase CO is found, which is in contrast to the situation in myoglobin. These observations explain why experimentally, unbound CO in the binuclear site of CcO has not been found as yet.

Structure and reaction mechanism of a novel enone reductase

Recently, a novel gut-bacterial fatty acid metabolism, saturation of polyunsaturated fatty acid, that modifies fatty acid composition of the host and is expected to improve our health by altering lipid metabolism related to the onset of metabolic syndrome, was discovered in Lactobacillus plantarum AKU 1009a. Enzymes constituting the pathway catalyze sequential reactions of free fatty acids without CoA or acyl carrier protein. Among these enzymes, CLA-ER was identified as an enone reductase that can saturate the C=C bond in the 10-oxo-trans-11-octadecenoic acid (KetoB) to produce 10-oxo-octadecanoic acid (KetoC). This enzyme is the sole member of the NADH oxidase/flavin reductase family that has been identified to exert an enone reduction activity. Here, we report both the structure of holo CLA-ER with cofactor FMN and the KetoC-bound structure, which elucidate the structural basis of enone group recognition of free fatty acids and provide the unique catalytic mechanism as an enone reductase in the NADH oxidase/flavin reductase family. A 'cap' structure of CLA-ER underwent a large conformational change upon KetoC binding. The resulting binding site adopts a sandglass shape and is positively charged at one side, which is suitable to recognize a fatty acid molecule with enone group. Based on the crystal structures and enzymatic activities of several mutants, we identified C51, F126 and Y101 as the critical residues for the reaction and proposed an alternative electron transfer pathway of CLA-ER. These findings expand our understanding of the complexity of fatty acid metabolism.

DATABASE

The atomic coordinates have been deposited in the Protein Data Bank (PDB), www.pdb.org (PDB ID 4QLX, 4QLY).

Proton-coupled electron transfer and the role of water molecules in proton pumping by cytochrome c oxidase

Molecular oxygen acts as the terminal electron sink in the respiratory chains of aerobic organisms. Cytochrome c oxidase in the inner membrane of mitochondria and the plasma membrane of bacteria catalyzes the reduction of oxygen to water, and couples the free energy of the reaction to proton pumping across the membrane. The proton-pumping activity contributes to the proton electrochemical gradient, which drives the synthesis of ATP. Based on kinetic experiments on the O-O bond splitting transition of the catalytic cycle (A → P(R)), it has been proposed that the electron transfer to the binuclear iron-copper center of O2 reduction initiates the proton pump mechanism. This key electron transfer event is coupled to an internal proton transfer from a conserved glutamic acid to the proton-loading site of the pump. However, the proton may instead be transferred to the binuclear center to complete the oxygen reduction chemistry, which would constitute a short-circuit. Based on atomistic molecular dynamics simulations of cytochrome c oxidase in an explicit membrane-solvent environment, complemented by related free-energy calculations, we propose that this short-circuit is effectively prevented by a redox-state-dependent organization of water molecules within the protein structure that gates the proton transfer pathway.

Higd1a is a positive regulator of cytochrome c oxidase

Cytochrome c oxidase (CcO) is the only enzyme that uses oxygen to produce a proton gradient for ATP production during mitochondrial oxidative phosphorylation. Although CcO activity increases in response to hypoxia, the underlying regulatory mechanism remains elusive. By screening for hypoxia-inducible genes in cardiomyocytes, we identified hypoxia inducible domain family, member 1A (Higd1a) as a positive regulator of CcO. Recombinant Higd1a directly integrated into highly purified CcO and increased its activity. Resonance Raman analysis revealed that Higd1a caused structural changes around heme a, the active center that drives the proton pump. Using a mitochondria-targeted ATP biosensor, we showed that knockdown of endogenous Higd1a reduced oxygen consumption and subsequent mitochondrial ATP synthesis, leading to increased cell death in response to hypoxia; all of these phenotypes were rescued by exogenous Higd1a. These results suggest that Higd1a is a previously unidentified regulatory component of CcO, and represents a therapeutic target for diseases associated with reduced CcO activity.

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.

Oxygen as Acceptor

Like most bacteria, Escherichia coli has a flexible and branched respiratory chain that enables the prokaryote to live under a variety of environmental conditions, from highly aerobic to completely anaerobic. In general, the bacterial respiratory chain is composed of dehydrogenases, a quinone pool, and reductases. Substrate-specific dehydrogenases transfer reducing equivalents from various donor substrates (NADH, succinate, glycerophosphate, formate, hydrogen, pyruvate, and lactate) to a quinone pool (menaquinone, ubiquinone, and dimethylmenoquinone). Then electrons from reduced quinones (quinols) are transferred by terminal reductases to different electron acceptors. Under aerobic growth conditions, the terminal electron acceptor is molecular oxygen. A transfer of electrons from quinol to O₂ is served by two major oxidoreductases (oxidases), cytochrome bo₃ encoded by cyoABCDE and cytochrome bd encoded by cydABX. Terminal oxidases of aerobic respiratory chains of bacteria, which use O₂ as the final electron acceptor, can oxidize one of two alternative electron donors, either cytochrome c or quinol. This review compares the effects of different inhibitors on the respiratory activities of cytochrome bo₃ and cytochrome bd in E. coli. It also presents a discussion on the genetics and the prosthetic groups of cytochrome bo₃ and cytochrome bd. The E. coli membrane contains three types of quinones that all have an octaprenyl side chain (C₄₀). It has been proposed that the bo₃ oxidase can have two ubiquinone-binding sites with different affinities. "WHAT'S NEW" IN THE REVISED ARTICLE: The revised article comprises additional information about subunit composition of cytochrome bd and its role in bacterial resistance to nitrosative and oxidative stresses. Also, we present the novel data on the electrogenic function of appBCX-encoded cytochrome bd-II, a second bd-type oxidase that had been thought not to contribute to generation of a proton motive force in E. coli, although its spectral properties closely resemble those of cydABX-encoded cytochrome bd.

Exploring O2 diffusion in A-type cytochrome c oxidases: molecular dynamics simulations uncover two alternative channels towards the binuclear site

Cytochrome c oxidases (Ccoxs) are the terminal enzymes of the respiratory chain in mitochondria and most bacteria. These enzymes couple dioxygen (O₂) reduction to the generation of a transmembrane electrochemical proton gradient. Despite decades of research and the availability of a large amount of structural and biochemical data available for the A-type Ccox family, little is known about the channel(s) used by O₂ to travel from the solvent/membrane to the heme a₃-Cu_B binuclear center (BNC). Moreover, the identification of all possible O₂ channels as well as the atomic details of O₂ diffusion is essential for the understanding of the working mechanisms of the A-type Ccox. In this work, we determined the O₂ distribution within Ccox from Rhodobacter sphaeroides, in the fully reduced state, in order to identify and characterize all the putative O₂ channels leading towards the BNC. For that, we use an integrated strategy combining atomistic molecular dynamics (MD) simulations (with and without explicit O₂ molecules) and implicit ligand sampling (ILS) calculations. Based on the 3D free energy map for O₂ inside Ccox, three channels were identified, all starting in the membrane hydrophobic region and connecting the surface of the protein to the BNC. One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the BNC. Both alternative O₂ channels start in the membrane spanning region and terminate close to Y288_I. These channels are a combination of multiple transiently interconnected hydrophobic cavities, whose opening and closure is regulated by the thermal fluctuations of the lining residues. Furthermore, our results show that, in this Ccox, the most likely (energetically preferred) routes for O₂ to reach the BNC are the alternative channels, rather than the X-ray inferred pathway.

Cytochrome c oxidases (Ccoxs), the terminal enzymes of the respiratory electron transport chain in eukaryotes and many prokaryotes, are key enzymes in aerobic respiration. These proteins couple the reduction of molecular dioxygen to water with the creation of a transmembrane electrochemical proton gradient. Over the last decades, most of the Ccoxs research focused on the mechanisms and energetics of reduction and/or proton pumping, and little emphasis has been given to the pathways used by dioxygen to reach the binuclear center, where dioxygen reduction takes place. In particular, the existence and the characteristics of the channel(s) used by O₂ to travel from the solvent/membrane to the binuclear site are still unclear. In this work, we combine all-atom molecular dynamics simulations and implicit ligand sampling calculations in order to identify and characterize the O₂ delivery channels in the Ccox from Rhodobacter sphaeroides. Altogether, our results suggest that, in this Ccox, O₂ can diffuse via three well-defined channels that start in membrane region (where O₂ solubility is higher than in the water). One of these channels corresponds to the pathway inferred from the X-ray data available, whereas the other two are alternative routes for O₂ to reach the binuclear center.

Understanding and applying tyrosine biochemical diversity

This review highlights some of the recent advances made in our understanding of the diversity of tyrosine biochemistry and shows how this has inspired novel applications in numerous areas of molecular design and synthesis, including chemical biology and bioconjugation. The pathophysiological implications of tyrosine biochemistry will be presented from a molecular perspective and the opportunities for therapeutic intervention explored.

Protein-protein interfaces from cytochrome c oxidase I evolve faster than nonbinding surfaces, yet negative selection is the driving force

Respiratory complexes are encoded by two genomes (mitochondrial DNA [mtDNA] and nuclear DNA [nDNA]). Although the importance of intergenomic coadaptation is acknowledged, the forces and constraints shaping such coevolution are largely unknown. Previous works using cytochrome c oxidase (COX) as a model enzyme have led to the so-called “optimizing interaction” hypothesis. According to this view, mtDNA-encoded residues close to nDNA-encoded residues evolve faster than the rest of positions, favoring the optimization of protein–protein interfaces. Herein, using evolutionary data in combination with structural information of COX, we show that failing to discern the effects of interaction from other structural and functional effects can lead to deceptive conclusions such as the “optimizing hypothesis.” Once spurious factors have been accounted for, data analysis shows that mtDNA-encoded residues engaged in contacts are, in general, more constrained than their noncontact counterparts. Nevertheless, noncontact residues from the surface of COX I subunit are a remarkable exception, being subjected to an exceptionally high purifying selection that may be related to the maintenance of a suitable heme environment. We also report that mtDNA-encoded residues involved in contacts with other mtDNA-encoded subunits are more constrained than mtDNA-encoded residues interacting with nDNA-encoded polypeptides. This differential behavior cannot be explained on the basis of predicted thermodynamic stability, as interactions between mtDNA-encoded subunits contribute more weakly to the complex stability than those interactions between subunits encoded by different genomes. Therefore, the higher conservation observed among mtDNA-encoded residues involved in intragenome interactions is likely due to factors other than structural stability.

Reconstruction of absolute absorption spectrum of reduced heme a in cytochrome C oxidase from bovine heart

This paper presents a new experimental approach for determining the individual optical characteristics of reduced heme a in bovine heart cytochrome c oxidase starting from a small selective shift of the heme a absorption spectrum induced by calcium ions. The difference spectrum induced by Ca2+ corresponds actually to a first derivative (differential) of the heme a(2+) absolute absorption spectrum. Such an absolute spectrum was obtained for the mixed-valence cyanide complex of cytochrome oxidase (a(2+)a3(3+)-CN) and was subsequently used as a basis spectrum for further procession and modeling. The individual absorption spectrum of the reduced heme a in the Soret region was reconstructed as the integral of the difference spectrum induced by addition of Ca2+. The spectrum of heme a(2+) in the Soret region obtained in this way is characterized by a peak with a maximum at 447 nm and half-width of 17 nm and can be decomposed into two Gaussians with maxima at 442 and 451 nm and half-widths of ~10 nm (589 cm(-1)) corresponding to the perpendicularly oriented electronic π→π* transitions B0x and B0y in the porphyrin ring. The reconstructed spectrum in the Soret band differs significantly from the "classical" absorption spectrum of heme a(2+) originally described by Vanneste (Vanneste, W. H. (1966) Biochemistry, 65, 838-848). The differences indicate that the overall γ-band of heme a(2+) in cytochrome oxidase contains in addition to the B0x and B0y transitions extra components that are not sensitive to calcium ions, or, alternatively, that the Vanneste's spectrum of heme a(2+) contains significant contribution from heme a3(2+). The reconstructed absorption band of heme a(2+) in the α-band with maximum at 605 nm and half-width of 18 nm (850 cm(-1)) corresponds most likely to the individual Q0y transition of heme a, whereas the Q0x transition contributes only weakly to the spectrum.

Infrared and Raman spectroscopic investigation of the reaction mechanism of cytochrome c oxidase

Recent progress in studies on the proton-pumping and O₂reduction mechanisms of cytochrome c oxidase (CcO) elucidated by infrared (IR) and resonance Raman (rR) spectroscopy, is reviewed. CcO is the terminal enzyme of the respiratory chain and its O₂reduction reaction is coupled with H⁺ pumping activity across the inner mitochondrial membrane. The former is catalyzed by heme a3 and its mechanism has been determined using a rR technique, while the latter used the protein moiety and has been investigated with an IR technique. The number of H⁺ relative to e⁻ transferred in the reaction is 1:1, and their coupling is presumably performed by heme a and nearby residues. To perform this function, different parts of the protein need to cooperate with each other spontaneously and sequentially. It is the purpose of this article to describe the structural details on the coupling on the basis of the vibrational spectra of certain specified residues and chromophores involved in the reaction. Recent developments in time-resolved IR and Raman technology concomitant with protein manipulation methods have yielded profound insights into such structural changes. In particular, the new IR techniques that yielded the breakthrough are reviewed and assessed in detail. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

Characterizing the proton loading site in cytochrome c oxidase

Cytochrome c oxidase (CcO) uses the energy released by reduction of O2 to H2O to drive eight charges from the high pH to low pH side of the membrane, increasing the electrochemical gradient. Four electrons and protons are used for chemistry, while four more protons are pumped. Proton pumping requires that residues on a pathway change proton affinity through the reaction cycle to load and then release protons. The protonation states of all residues in CcO are determined in MultiConformational Continuum Electrostatics simulations with the protonation and redox states of heme a, a3, Cu(B), Y288, and E286 used to define the catalytic cycle. One proton is found to be loaded and released from residues identified as the proton loading site (PLS) on the P-side of the protein in each of the four CcO redox states. Thus, the same proton pumping mechanism can be used each time CcO is reduced. Calculations with structures of Rhodobacter sphaeroides, Paracoccus denitrificans, and bovine CcO derived by crystallography and molecular dynamics show the PLS functions similarly in different CcO species. The PLS is a cluster rather than a single residue, as different structures show 1-4 residues load and release protons. However, the proton affinity of the heme a3 propionic acids primarily determines the number of protons loaded into the PLS; if their proton affinity is too low, less than one proton is loaded.

The pore-lining regions in cytochrome c oxidases: A computational analysis of caveolin, cholesterol and transmembrane helix contributions to proton movement

Cytochrome c oxidase (CcO) is the terminal enzyme in the electron transfer chain. CcO catalyzes a four electron reduction of O2 to water at a catalytic site formed by high-spin heme (a3) and copper atoms (CuB). While it is recognized that proton movement is coupled to oxygen reduction, the proton channel(s) have not been well defined. Using computational methods developed to study protein topology, membrane channels and 3D packing arrangements within transmembrane (TM) helix arrays, we find that subunit-1 (COX-1), subunit-2 (COX-2) and subunit-3 (COX-3) contribute 139, 46 and 25 residues, respectively, to channel formation between the mitochondrial matrix and intermembrane space. Nine of 12 TM helices in COX-1, both helices in COX-2 and 5 of the 6 TM helices in COX-3 are pore-lining regions (possible channel formers). Heme a3 and the CuB sites (as well as the CuA center of COX-2) are located within the channel that includes TM-6, TM-7, TM-10 and TM-11 of COX-1 and are associated with multiple cholesterol and caveolin-binding (CB) motifs. Sequence analysis identifies five CB motifs within COX-1, two within COX-2 and four within COX-3; each caveolin containing a pore-lining helix C-terminal to a TM helix-turn-helix. Channel formation involves interaction between multiple pore-lining regions within protein subunits and/or dimers. PoreWalker analysis lends support to the D-channel model of proton translocation. Under physiological conditions, caveolins may introduce channel formers juxtaposed to those in COX-1, COX-2 and COX-3, which together with cholesterol may form channel(s) essential for proton translocation through the inner mitochondrial membrane.

Constant pH Molecular Dynamics in Explicit Solvent with Enveloping Distribution Sampling and Hamiltonian Exchange

We present a new computational approach for constant pH simulations in explicit solvent based on the combination of the enveloping distribution sampling (EDS) and Hamiltonian replica exchange (HREX) methods. Unlike constant pH methods based on variable and continuous charge models, our method is based on discrete protonation states. EDS generates a hybrid Hamiltonian of different protonation states. A smoothness parameter s is used to control the heights of energy barriers of the hybrid-state energy landscape. A small s value facilitates state transitions by lowering energy barriers. Replica exchange between EDS potentials with different s values allows us to readily obtain a thermodynamically accurate ensemble of multiple protonation states with frequent state transitions. The analysis is performed with an ensemble obtained from an EDS Hamiltonian without smoothing, s = ∞, which strictly follows the minimum energy surface of the end states. The accuracy and efficiency of this method is tested on aspartic acid, lysine, and glutamic acid, which have two protonation states, a histidine with three states, a four-residue peptide with four states, and snake cardiotoxin with eight states. The pKa values estimated with the EDS-HREX method agree well with the experimental pKa values. The mean absolute errors of small benchmark systems range from 0.03 to 0.17 pKa units, and those of three titratable groups of snake cardiotoxin range from 0.2 to 1.6 pKa units. This study demonstrates that EDS-HREX is a potent theoretical framework, which gives the correct description of multiple protonation states and good calculated pKa values.

The pathway of O₂to the active site in heme-copper oxidases

The route of O₂to and from the high-spin heme in heme-copper oxidases has generally been believed to emulate that of carbon monoxide (CO). Time-resolved and stationary infrared experiments in our laboratories of the fully reduced CO-bound enzymes, as well as transient optical absorption saturation kinetics studies as a function of CO pressure, have provided strong support for CO binding to CuB⁺ on the pathway to and from the high-spin heme. The presence of CO on CuB⁺ suggests that O₂binding may be compromised in CO flow-flash experiments. Time-resolved optical absorption studies show that the rate of O₂and NO binding in the bovine enzyme (1 × 10⁸M⁻¹s⁻¹) is unaffected by the presence of CO, which is consistent with the rapid dissociation (t½ = 1.5μs) of CO from CuB⁺. In contrast, in Thermus thermophilus (Tt) cytochrome ba3 the O₂and NO binding to heme a3 slows by an order of magnitude in the presence of CO (from 1 × 10⁹ to 1 × 10⁸M⁻¹s⁻¹), but is still considerably faster (~10μs at 1atm O₂) than the CO off-rate from CuB in the absence of O₂(milliseconds). These results show that traditional CO flow-flash experiments do not give accurate results for the physiological binding of O₂and NO in Tt ba3, namely, in the absence of CO. They also raise the question whether in CO flow-flash experiments on Tt ba3 the presence of CO on CuB⁺ impedes the binding of O₂to CuB⁺ or, if O₂does not bind to CuB⁺ prior to heme a3, whether the CuB⁺-CO complex sterically restricts access of O₂to the heme. Both possibilities are discussed, and we argue that O₂binds directly to heme a3 in Tt ba3, causing CO to dissociate from CuB⁺ in a concerted manner through steric and/or electronic effects. This would allow CuB⁺ to function as an electron donor during the fast (5μs) breaking of the OO bond. These results suggest that the binding of CO to CuB⁺ on the path to and from heme a3 may not be applicable to O₂and NO in all heme-copper oxidases. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

How hydrogen peroxide is metabolized by oxidized cytochrome c oxidase

In the absence of external electron donors, oxidized bovine cytochrome c oxidase (CcO) exhibits the ability to decompose excess H2O2. Depending on the concentration of peroxide, two mechanisms of degradation were identified. At submillimolar peroxide concentrations, decomposition proceeds with virtually no production of superoxide and oxygen. In contrast, in the millimolar H2O2 concentration range, CcO generates superoxide from peroxide. At submillimolar concentrations, the decomposition of H2O2 occurs at least at two sites. One is the catalytic heme a3-CuB center where H2O2 is reduced to water. During the interaction of the enzyme with H2O2, this center cycles back to oxidized CcO via the intermediate presence of two oxoferryl states. We show that at pH 8.0 two molecules of H2O2 react with the catalytic center accomplishing one cycle. In addition, the reactions at the heme a3-CuB center generate the surface-exposed lipid-based radical(s) that participates in the decomposition of peroxide. It is also found that the irreversible decline of the catalytic activity of the enzyme treated with submillimolar H2O2 concentrations results specifically from the decrease in the rate of electron transfer from heme a to the heme a3-CuB center during the reductive phase of the catalytic cycle. The rates of electron transfer from ferrocytochrome c to heme a and the kinetics of the oxidation of the fully reduced CcO with O2 were not affected in the peroxide-modified CcO.

IR signatures of the metal centres of bovine cytochrome c oxidase: assignments and redox-linkage

Assignments of IR bands of reduced minus oxidized IR difference spectra of bovine and related cytochrome c oxidases are reviewed and their linkages to specific metal centres are assessed. To aid this, redox-poised difference spectra in the presence of cyanide or carbon monoxide are presented. These ligands fix the redox states of either haem a3 alone or haem a3 and CuB respectively, while allowing redox cycling of the remaining centres.

Cardiac mitochondria and reactive oxygen species generation

Mitochondrial reactive oxygen species (ROS) have emerged as an important mechanism of disease and redox signaling in the cardiovascular system. Under basal or pathological conditions, electron leakage for ROS production is primarily mediated by the electron transport chain and the proton motive force consisting of a membrane potential (ΔΨ) and a proton gradient (ΔpH). Several factors controlling ROS production in the mitochondria include flavin mononucleotide and flavin mononucleotide-binding domain of complex I, ubisemiquinone and quinone-binding domain of complex I, flavin adenine nucleotide-binding moiety and quinone-binding pocket of complex II, and unstable semiquinone mediated by the Q cycle of complex III. In mitochondrial complex I, specific cysteinyl redox domains modulate ROS production from the flavin mononucleotide moiety and iron-sulfur clusters. In the cardiovascular system, mitochondrial ROS have been linked to mediating the physiological effects of metabolic dilation and preconditioning-like mitochondrial ATP-sensitive potassium channel activation. Furthermore, oxidative post-translational modification by glutathione in complex I and complex II has been shown to affect enzymatic catalysis, protein-protein interactions, and enzyme-mediated ROS production. Conditions associated with oxidative or nitrosative stress, such as myocardial ischemia and reperfusion, increase mitochondrial ROS production via oxidative injury of complexes I and II and superoxide anion radical-induced hydroxyl radical production by aconitase. Further insight into cellular mechanisms by which specific redox post-translational modifications regulate ROS production in the mitochondria will enrich our understanding of redox signal transduction and identify new therapeutic targets for cardiovascular diseases in which oxidative stress perturbs normal redox signaling.