Cytochrome c oxidase: structure and spectroscopy

Cytochrome c oxidase dysfunction in oxidative stress

Cytochrome c oxidase (CcO) is the terminal oxidase of the mitochondrial electron transport chain. This bigenomic enzyme in mammals contains 13 subunits of which the 3 catalytic subunits are encoded by the mitochondrial genes. The remaining 10 subunits with suspected roles in the regulation, and/or assembly, are coded by the nuclear genome. The enzyme contains two heme groups (heme a and a3) and two Cu(2+) centers (Cu(2+) A and Cu(2+) B) as catalytic centers and handles more than 90% of molecular O(2) respired by the mammalian cells and tissues. CcO is a highly regulated enzyme which is believed to be the pacesetter for mitochondrial oxidative metabolism and ATP synthesis. The structure and function of the enzyme are affected in a wide variety of diseases including cancer, neurodegenerative diseases, myocardial ischemia/reperfusion, bone and skeletal diseases, and diabetes. Despite handling a high O(2) load the role of CcO in the production of reactive oxygen species still remains a subject of debate. However, a volume of evidence suggests that CcO dysfunction is invariably associated with increased mitochondrial reactive oxygen species production and cellular toxicity. In this paper we review the literature on mechanisms of multimodal regulation of CcO activity by a wide spectrum of physiological and pathological factors. We also review an array of literature on the direct or indirect roles of CcO in reactive oxygen species production.

Spectroscopic and kinetic investigation of the fully reduced and mixed valence states of ba3-cytochrome c oxidase from Thermus thermophilus: a Fourier transform infrared (FTIR) and time-resolved step-scan FTIR study

The complete understanding of a molecular mechanism of action requires the thermodynamic and kinetic characterization of different states and intermediates. Cytochrome c oxidase reduces O(2) to H(2)O, a reaction coupled to proton translocation across the membrane. Therefore, it is necessary to undertake a thorough characterization of the reduced form of the enzyme and the determination of the electron transfer processes and pathways between the redox-active centers. In this study Fourier transform infrared (FTIR) and time-resolved step-scan FTIR spectroscopy have been applied to study the fully reduced and mixed valence states of cytochrome ba(3) from Thermus thermophilus. We used as probe carbon monoxide (CO) to characterize both thermodynamically and kinetically the cytochrome ba(3)-CO complex in the 5.25-10.10 pH/pD range and to study the reverse intramolecular electron transfer initiated by the photolysis of CO in the two-electron reduced form. The time-resolved step-scan FTIR data revealed no pH/pD dependence in both the decay of the transient Cu(B)(1+)-CO complex and rebinding to heme a(3) rates, suggesting that no structural change takes place in the vicinity of the binuclear center. Surprisingly, photodissociation of CO from the mixed valence form of the enzyme does not lead to reverse electron transfer from the reduced heme a(3) to the oxidized low-spin heme b, as observed in all the other aa(3) and bo(3) oxidases previously examined. The heme b-heme a(3) electron transfer is guaranteed, and therefore, there is no need for structural rearrangements and complex synchronized cooperativities. Comparison among the available structures of ba(3)- and aa(3)-cytochrome c oxidases identifies possible active pathways involved in the electron transfer processes and key structural elements that contribute to the different behavior observed in cytochrome ba(3).

Frequency matrix approach demonstrates high sequence quality in avian BARCODEs and highlights cryptic pseudogenes

The accuracy of DNA barcode databases is critical for research and practical applications. Here we apply a frequency matrix to assess sequencing errors in a very large set of avian BARCODEs. Using 11,000 sequences from 2,700 bird species, we show most avian cytochrome c oxidase I (COI) nucleotide and amino acid sequences vary within a narrow range. Except for third codon positions, nearly all (96%) sites were highly conserved or limited to two nucleotides or two amino acids. A large number of positions had very low frequency variants present in single individuals of a species; these were strongly concentrated at the ends of the barcode segment, consistent with sequencing error. In addition, a small fraction (0.1%) of BARCODEs had multiple very low frequency variants shared among individuals of a species; these were found to represent overlooked cryptic pseudogenes lacking stop codons. The calculated upper limit of sequencing error was 8 × 10(-5) errors/nucleotide, which was relatively high for direct Sanger sequencing of amplified DNA, but unlikely to compromise species identification. Our results confirm the high quality of the avian BARCODE database and demonstrate significant quality improvement in avian COI records deposited in GenBank over the past decade. This approach has potential application for genetic database quality control, discovery of cryptic pseudogenes, and studies of low-level genetic variation.

Interaction between overtraining and the interindividual variability may (not) trigger muscle oxidative stress and cardiomyocyte apoptosis in rats

Severe endurance training (overtraining) may cause underperformance related to muscle oxidative stress and cardiomyocyte alterations. Currently, such relationship has not been empirically established. In this study, Wistar rats (n = 19) underwent eight weeks of daily exercise sessions followed by three overtraining weeks in which the daily frequency of exercise sessions increased. After the 11th training week, eight rats exhibited a reduction of 38% in performance (nonfunctional overreaching group (NFOR)), whereas eleven rats exhibited an increase of 18% in performance (functional overreaching group (FOR)). The red gastrocnemius of NFOR presented significantly lower citrate synthase activity compared to FOR, but similar to that of the control. The activity of mitochondrial complex IV in NFOR was lower than that of the control and FOR. This impaired mitochondrial adaptation in NFOR was associated with increased antioxidant enzyme activities and increased lipid peroxidation (in muscle and plasma) relative to FOR and control. Cardiomyocyte apoptosis was higher in NFOR. Plasma creatine kinase levels were unchanged. We observed that some rats that presented evidence of muscle oxidative stress are also subject to cardiomyocyte apoptosis under endurance overtraining. Blood lipid peroxides may be a suitable biomarker for muscle oxidative stress that is unrelated to severe muscle damage.

Mobility of Xe atoms within the oxygen diffusion channel of cytochrome ba(3) oxidase

We use a form of "freeze-trap, kinetic crystallography" to explore the migration of Xe atoms away from the dinuclear heme a(3)/Cu(B) center in Thermus thermophilus cytochrome ba(3) oxidase. This enzyme is a member of the heme-copper oxidase superfamily and is thus crucial for dioxygen-dependent life. The mechanisms involved in the migration of oxygen, water, electrons, and protons into and/or out of the specialized channels of the heme-copper oxidases are generally not well understood. Pressurization of crystals with Xe gas previously revealed a O(2) diffusion channel in cytochrome ba(3) oxidase that is continuous, Y-shaped, 18-20 Å in length and comprised of hydrophobic residues, connecting the protein surface within the bilayer to the a(3)-Cu(B) center in the active site. To understand movement of gas molecules within the O(2) channel, we performed crystallographic analysis of 19 Xe laden crystals freeze-trapped in liquid nitrogen at selected times between 0 and 480 s while undergoing outgassing at room temperature. Variation in Xe crystallographic occupancy at five discrete sites as a function of time leads to a kinetic model revealing relative degrees of mobility of Xe atoms within the channel. Xe egress occurs primarily through the channel formed by the Xe1 → Xe5 → Xe3 → Xe4 sites, suggesting that ingress of O(2) is likely to occur by the reverse of this process. The channel itself appears not to undergo significant structural changes during Xe migration, thereby indicating a passive role in this important physiological function.

Bioenergetics at extreme temperature: Thermus thermophilus ba(3)- and caa(3)-type cytochrome c oxidases

Seven years into the completion of the genome sequencing projects of the thermophilic bacterium Thermus thermophilus strains HB8 and HB27, many questions remain on its bioenergetic mechanisms. A key fact that is occasionally overlooked is that oxygen has a very limited solubility in water at high temperatures. The HB8 strain is a facultative anaerobe whereas its relative HB27 is strictly aerobic. This has been attributed to the absence of nitrate respiration genes from the HB27 genome that are carried on a mobilizable but highly-unstable plasmid. In T. thermophilus, the nitrate respiration complements the primary aerobic respiration. It is widely known that many organisms encode multiple biochemically-redundant components of the respiratory complexes. In this minireview, the presence of the two cytochrome c oxidases (CcO) in T. thermophilus, the ba(3)- and caa(3)-types, is outlined along with functional considerations. We argue for the distinct evolutionary histories of these two CcO including their respective genetic and molecular organizations, with the caa(3)-oxidase subunits having been initially 'fused'. Coupled with sequence analysis, the ba(3)-oxidase crystal structure has provided evolutionary and functional information; for example, its subunit I is more closely related to archaeal sequences than bacterial and the substrate-enzyme interaction is hydrophobic as the elevated growth temperature weakens the electrostatic interactions common in mesophiles. Discussion on the role of cofactors in intra- and intermolecular electron transfer and proton pumping mechanism is also included.

Factors that control catalytic two- versus four-electron reduction of dioxygen by copper complexes

The selective two-electron reduction of O(2) by one-electron reductants such as decamethylferrocene (Fc*) and octamethylferrocene (Me(8)Fc) is efficiently catalyzed by a binuclear Cu(II) complex [Cu(II)(2)(LO)(OH)](2+) (D1) {LO is a binucleating ligand with copper-bridging phenolate moiety} in the presence of trifluoroacetic acid (HOTF) in acetone. The protonation of the hydroxide group of [Cu(II)(2)(LO)(OH)](2+) with HOTF to produce [Cu(II)(2)(LO)(OTF)](2+) (D1-OTF) makes it possible for this to be reduced by 2 equiv of Fc* via a two-step electron-transfer sequence. Reactions of the fully reduced complex [Cu(I)(2)(LO)](+) (D3) with O(2) in the presence of HOTF led to the low-temperature detection of the absorption spectra due to the peroxo complex [Cu(II)(2)(LO)(OO)] (D) and the protonated hydroperoxo complex [Cu(II)(2)(LO)(OOH)](2+) (D4). No further Fc* reduction of D4 occurs, and it is instead further protonated by HOTF to yield H(2)O(2) accompanied by regeneration of [Cu(II)(2)(LO)(OTF)](2+) (D1-OTF), thus completing the catalytic cycle for the two-electron reduction of O(2) by Fc*. Kinetic studies on the formation of Fc*(+) under catalytic conditions as well as for separate examination of the electron transfer from Fc* to D1-OTF reveal there are two important reaction pathways operating. One is a rate-determining second reduction of D1-OTF, thus electron transfer from Fc* to a mixed-valent intermediate [Cu(II)Cu(I)(LO)](2+) (D2), which leads to [Cu(I)(2)(LO)](+) that is coupled with O(2) binding to produce [Cu(II)(2)(LO)(OO)](+) (D). The other involves direct reaction of O(2) with the mixed-valent compound D2 followed by rapid Fc* reduction of a putative superoxo-dicopper(II) species thus formed, producing D.

Exploration of respiratory chain of Nocardia asteroides: purification of succinate quinone oxidoreductase

Nocardia asteroides is a pathogenic bacterium that causes severe pulmonary infections and plays a vital role in HIV development. Its electron transport chain containing cytochromes as electron carriers is still undiscovered. Information regarding cytochromes is important during drug synthesis based on cytochrome inhibitions. In this study we explored the electron transport of N. asteroides. Spectroscopic analysis of cytoplasm and membranes isolated from N. asteroides indicates the presence of soluble cytochrome-c, complex-II and the modified a(1)c(1) complex as the terminal oxidase. The molecular weight of the respiratory complex-II isolated and purified from the given bacterium was 103 kDa and was composed of three subunits, of 14, 26 and 63 kDa. Complex-II showed symmetrical α-absorption peaks at 561 nm in the reduced state. Spectral analysis revealed the presence of only one heme b molecule (14-kDa subunit) in complex-II, which was confirmed by heme staining. Heme b content was found to be 9.5 nmol/mg in complex-II. The electron transport chain of N. asteroides showed the presence of soluble cytochrome-c, cytochrome-a(1)c(1) and cytochrome-b.

Net proton uptake is preceded by multiple proton transfer steps upon electron injection into cytochrome c oxidase

Cytochrome c oxidase (COX), the last enzyme of the respiratory chain of aerobic organisms, catalyzes the reduction of molecular oxygen to water. It is a redox-linked proton pump, whose mechanism of proton pumping has been controversially discussed, and the coupling of proton and electron transfer is still not understood. Here, we investigated the kinetics of proton transfer reactions following the injection of a single electron into the fully oxidized enzyme and its transfer to the hemes using time-resolved absorption spectroscopy and pH indicator dyes. By comparison of proton uptake and release kinetics observed for solubilized COX and COX-containing liposomes, we conclude that the 1-μs electron injection into Cu(A), close to the positive membrane side (P-side) of the enzyme, already results in proton uptake from both the P-side and the N (negative)-side (1.5 H(+)/COX and 1 H(+)/COX, respectively). The subsequent 10-μs transfer of the electron to heme a is accompanied by the release of 1 proton from the P-side to the aqueous bulk phase, leaving ∼0.5 H(+)/COX at this side to electrostatically compensate the charge of the electron. With ∼200 μs, all but 0.4 H(+) at the N-side are released to the bulk phase, and the remaining proton is transferred toward the hemes to a so-called "pump site." Thus, this proton may already be taken up by the enzyme as early as during the first electron transfer to Cu(A). These results support the idea of a proton-collecting antenna, switched on by electron injection.

Neurometabolic mechanisms for memory enhancement and neuroprotection of methylene blue

This paper provides the first review of the memory-enhancing and neuroprotective metabolic mechanisms of action of methylene blue in vivo. These mechanisms have important implications as a new neurobiological approach to improve normal memory and to treat memory impairment and neurodegeneration associated with mitochondrial dysfunction. Methylene blue's action is unique because its neurobiological effects are not determined by regular drug-receptor interactions or drug-response paradigms. Methylene blue shows a hormetic dose-response, with opposite effects at low and high doses. At low doses, methylene blue is an electron cycler in the mitochondrial electron transport chain, with unparalleled antioxidant and cell respiration-enhancing properties that affect the function of the nervous system in a versatile manner. A major role of the respiratory enzyme cytochrome oxidase on the memory-enhancing effects of methylene blue is supported by available data. The memory-enhancing effects have been associated with improvement of memory consolidation in a network-specific and use-dependent fashion. In addition, low doses of methylene blue have also been used for neuroprotection against mitochondrial dysfunction in humans and experimental models of disease. The unique auto-oxidizing property of methylene blue and its pleiotropic effects on a number of tissue oxidases explain its potent neuroprotective effects at low doses. The evidence reviewed supports a mechanistic role of low-dose methylene blue as a promising and safe intervention for improving memory and for the treatment of acute and chronic conditions characterized by increased oxidative stress, neurodegeneration and memory impairment.

The superfamily of heme-copper oxygen reductases: types and evolutionary considerations

Heme-copper oxygen reductases (HCO) reduce O(2) to water being the last enzymatic complexes of most aerobic respiratory chains. These enzymes promote energy conservation coupling the catalytic reaction to charge separation and charge translocation across the prokaryotic cytoplasmatic or mitochondrial membrane. In this way they contribute to the establishment and maintenance of the transmembrane difference of electrochemical potential, which is vital for solute/nutrient cell import, synthesis of ATP and motility. The HCO enzymes most probably share with the nitric oxide reductases, NORs, a common ancestor. We have proposed the classification of HCOs into three different types, A, B and C; based on the constituents of their proton channels (Pereira, Santana and Teixeira (2001) Biochim Biophys Acta, 1505, 185-208). This classification was recently challenged by the suggestion of other different types of HCOs. Using an enlarged sampling we performed an exhaustive bioinformatic reanalysis of HCOs family. Our results strengthened our previously proposed classification and showed no need for the existence of more divisions. Now, we analyze the taxonomic distribution of HCOs and NORs and the congruence of their sequence trees with the 16S rRNA tree. We observed that HCOs are widely distributed in the two prokaryotic domains and that the different types of enzymes are not confined to a specific taxonomic group or environmental niche.

Distinguishing two groups of flavin reductases by analyzing the protonation state of an active site carboxylic acid

Flavin-containing reductases are involved in a wide variety of physiological reactions such as photosynthesis, nitric oxide synthesis, and detoxification of foreign compounds, including therapeutic drugs. Ferredoxin-NADP(H)-reductase (FNR) is the prototypical enzyme of this family. The fold of this protein is highly conserved and occurs as one domain of several multidomain enzymes such as the members of the diflavin reductase family. The enzymes of this family have emerged as fusion of a FNR and a flavodoxin. Although the active sites of these enzymes are very similar, different enzymes function in opposite directions, that is, some reduce oxidized nicotinamide adenine dinucleotide phosphate (NADP(+)) and some oxidize reduced nicotinamide adenine dinucleotide phosphate (NADPH). In this work, we analyze the protonation behavior of titratable residues of these enzymes through electrostatic calculations. We find that a highly conserved carboxylic acid in the active site shows a different titration behavior in different flavin reductases. This residue is deprotonated in flavin reductases present in plastids, but protonated in bacterial counterparts and in diflavin reductases. The protonation state of the carboxylic acid may also influence substrate binding. The physiological substrate for plastidic enzymes is NADP(+), but it is NADPH for the other mentioned reductases. In this article, we discuss the relevance of the environment of this residue for its protonation and its importance in catalysis. Our results allow to reinterpret and explain experimental data.

A bioinformatics classifier and database for heme-copper oxygen reductases

Background

Heme-copper oxygen reductases (HCOs) are the last enzymatic complexes of most aerobic respiratory chains, reducing dioxygen to water and translocating up to four protons across the inner mitochondrial membrane (eukaryotes) or cytoplasmatic membrane (prokaryotes). The number of completely sequenced genomes is expanding exponentially, and concomitantly, the number and taxonomic distribution of HCO sequences. These enzymes were initially classified into three different types being this classification recently challenged.

Methodology

We reanalyzed the classification scheme and developed a new bioinformatics classifier for the HCO and Nitric oxide reductases (NOR), which we benchmark against a manually derived gold standard sequence set. It is able to classify any given sequence of subunit I from HCO and NOR with a global recall and precision both of 99.8%. We use this tool to classify this protein family in 552 completely sequenced genomes.

Conclusions

We concluded that the new and broader data set supports three functional and evolutionary groups of HCOs. Homology between NORs and HCOs is shown and NORs closest relationship with C Type HCOs demonstrated. We established and made available a classification web tool and an integrated Heme-Copper Oxygen reductase and NOR protein database (www.evocell.org/hco).

The diagnosis of acute lower limb compartment syndrome: Applications of near infrared spectroscopy

Acute compartment syndrome of the lower limb is a significant problem in surgical practice, the successful management of which depends upon swift diagnosis and intervention.

Conventionally, diagnosis has been based upon clinical assessment; however, this can be unreliable and the potential for missed compartment syndrome remains. The supplementary use of compartment pressure monitoring has addressed some of these issues, but it remains an invasive technique, the exact role of which is still debated in the literature.

Near infrared spectroscopy (NIRS) is an emerging technique in medical practice which provides a non-invasive, continuous and real time measure of local tissue oxygenation. Early experimental work and subsequent clinical studies have demonstrated that NIRS provides an accurate means of detecting compartment syndrome, and that its sensitivity in some circumstances may exceed that of monitoring compartment pressures. Despite this promise, limitations of the technique, such as difficultly monitoring the deep posterior compartment of leg using current systems, and the relative expense of the equipment, have hindered broader adoption.

Electron transfer pathways in cytochrome c oxidase

Mixed quantum mechanical/molecular mechanics calculations were used to explore the electron pathway of the terminal electron transfer enzyme, cytochrome c oxidase. This enzyme catalyzes the reduction of molecular oxygen to water in a multiple step process. Density functional calculations on the three redox centers allowed for the characterization of the electron transfer mechanism, following the sequence Cu(A)→heme a→heme a(3). This process is largely affected by the presence of positive charges, confirming the possibility of a proton coupled electron transfer. An extensive mapping of all residues involved in the electron transfer, between the Cu(A) center (donor) and the O(2) reduction site heme a(3)-Cu(B) (receptor), was obtained by selectively activating/deactivating different quantum regions. The method employed, called QM/MM e-pathway, allowed the identification of key residues along the possible electron transfer paths, consistent with experimental data. In particular, the role of arginines 481 and 482 appears crucial in the Cu(A)→heme a and in the heme a→heme a(3) electron transfer processes. This article is part of a Special Issue entitled: Allosteric cooperativity in respiratory proteins.

Exploration of the cytochrome c oxidase pathway puzzle and examination of the origin of elusive mutational effects

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is a problem of great current interest. Despite promising mechanistic proposals, so far, a physically consistent model that would reproduce all the relevant barriers needed to create a working pump has not been presented. In addition, there are major problems in elucidating the origin of key mutational effects and in understanding the nature of the apparent pK(a) values associated with the pH dependencies of specific proton transfer (PT) reactions in CcO. This work takes a key step in resolving the above problems, by considering mutations, such as the Asn139Asp replacement, that blocks proton pumping without affecting PT to the catalytic site. We first introduce a formulation that makes it possible to relate the apparent pK(a) of Glu286 to different conformational states of this residue. We then use the new formulation along with the calculated pK(a) values of Glu286 at these different conformations to reproduce the experimentally observed apparent pK(a) of the residue. Next, we take the X-ray structures of the native and Asn139Asp mutant of the Paracoccus denitrificans CcO (N131D in this system) and reproduce for the first time the change in the primary PT pathways (and other key features) based on simulations that start with the observed structural changes. We also consider the competition between proton transport to the catalytic site and the pump site, as a function of the bulk pH, as well as the H/D isotope effect, and use this information to explore the relative height of the two barriers. The paper emphasizes the crucial role of energy-based considerations that include the PT process, and the delicate control of PT in CcO.

Mutagenesis of tyrosine residues within helix VII in subunit I of the cytochrome cbb₃ oxidase from Rhodobacter capsulatus

The cbb (3)-type oxidases are members of the heme-copper oxidase superfamily, distant by sequence comparisons, but sharing common functional characteristics. The cbb (3) oxidases are missing an active-site tyrosine residue that is absolutely conserved in all A and B-type heme-copper oxidases. This tyrosine is known to play a critical role in the catalytic mechanisms of A and B-type oxidases. The absence of this tyrosine in the cbb (3) oxidases raises the possibility that the cbb (3) oxidases utilize a different catalytic mechanism from that of the other members of the superfamily, or have this conserved residue in different helices. Recently sequence comparisons indicate that, a tyrosine residues that might be analogous to the active-site tyrosine in other oxidases are present in the cbb (3) oxidases but these tyrosines originates from a different transmembrane helix within the protein. In this research, three conserved tyrosine residues, Y294, Y308 and Y318, in helix VII were substituted for phenylalanine. Y318F mutant in the Rhodobacter capsulatus oxidase resulted in a fully assembled enzyme with nativelike structure and activity, but Y294F mutant is not assembled and have a catalytic activity. On the other hand, Y308F mutant is fully assembled enzyme with nativelike structure, but lacking catalytic activity. This result indicates that Y308 should be crucial in catalytic activity of the cbb (3) oxidase of R. capsulatus. These findings support the assumption that all of the heme-copper oxidases utilize the same catalytic mechanism and provide a residue originates from different places within the primary sequence for different members of the same superfamily.

Genes or culture: are mitochondrial genes associated with tool use in bottlenose dolphins (Tursiops sp.)?

Some bottlenose dolphins use marine sponges as foraging tools ('sponging'), which appears to be socially transmitted from mothers mainly to their female offspring. Yet, explanations alternative to social transmission have been proposed. Firstly, the propensity to engage in sponging might be due to differences in diving ability caused by variation of mitochondrial genes coding for proteins of the respiratory chain. Secondly, the cultural technique of sponging may have selected for changes in these same genes (or other autosomal ones) among its possessors. We tested whether sponging can be predicted by mitochondrial coding genes and whether these genes are under selection. In 29 spongers and 54 non-spongers from two study sites, the non-coding haplotype at the HVRI locus was a significant predictor of sponging, whereas the coding mitochondrial genes were not. There was no evidence of selection in the investigated genes. Our study shows that mitochondrial gene variation is unlikely to be a viable alternative to cultural transmission as a primary driver of tool use in dolphins.

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

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

Vibrational resonances and CuB displacement controlled by proton motion in cytochrome c oxidase

Cytochrome c oxidase (CcO), found in the inner mitochondrial membranes or in many bacteria, catalyzes the four-electron reduction of molecular oxygen to water. Four protons are pumped across the inner mitochondrial membrane through CcO. In this study, quantum mechanics/molecular mechanics and molecular dynamics calculations are used to probe the spectroscopic characteristics of the ferryl intermediates in the aa(3) CcO/O(2) reaction. These highly elaborate calculations, supported by several calculations on smaller model systems, demonstrate the sensitivity of vibrational frequencies on the Coulombic field of heme a(3) and their dependence on the distance of the adjacent Cu(B) to the heme a(3)-Fe atom. This distance seems to be associated with the protonation state of the heme a(3) propionate A, and we propose that it plays a crucial role on the mechanism of action of CcO. In detail, we link proton pumping activity in CcO enzyme (a) to a multiple (1:1:2) resonance among the frequencies of Fe(IV)=O bond stretching, the breathing mode of Histidine 411, and a bending mode of the His411-Fe(IV)=O species (aa(3) from Paracoccus denitrificans numbering) and (b) to Cu(B) displacement by electrostatic interactions toward the heme a(3) iron. We find that the vibrations of the His411-Fe(IV)=O unit become highly coupled depending on the protonation state of the heme a(3) ring A propionate/Asp399 pair, and we propose a mechanism for the resonance Raman enhancement of the bending mode delta(His411-Fe(IV)=O). Calculations on model systems demonstrate that the position of Cu(B) in relation to heme a(3) iron-oxo plays a crucial role in regulating that resonance. We also discuss the origin of the coupling between bending, delta(His411-Fe(IV)=O) and nu(Fe=O) stretching modes, and the role played by such vibrational coupling interactions or Cu(B) position in controlling functional properties of the enzyme, including electron/proton coupling as well as experimental spectra.

Multifrequency pulsed electron paramagnetic resonance on metalloproteins

Metalloproteins often contain metal centers that are paramagnetic in some functional state of the protein; hence electron paramagnetic resonance (EPR) spectroscopy can be a powerful tool for studying protein structure and function. Dipolar spectroscopy allows the determination of the dipole-dipole interactions between metal centers in protein complexes, revealing the structural arrangement of different paramagnetic centers at distances of up to 8 nm. Hyperfine spectroscopy can be used to measure the interaction between an unpaired electron spin and nuclear spins within a distance of 0.8 nm; it therefore permits the characterization of the local structure of the paramagnetic center's ligand sphere with very high precision. In this Account, we review our laboratory's recent applications of both dipolar and hyperfine pulsed EPR methods to metalloproteins. We used pulsed dipolar relaxation methods to investigate the complex of cytochrome c and cytochrome c oxidase, a noncovalent protein-protein complex involved in mitochondrial electron-transfer reactions. Hyperfine sublevel correlation spectroscopy (HYSCORE) was used to study the ligand sphere of iron-sulfur clusters in complex I of the mitochondrial respiratory chain and substrate binding to the molybdenum enzyme polysulfide reductase. These examples demonstrate the potential of the two techniques; however, they also highlight the difficulties of data interpretation when several paramagnetic species with overlapping spectra are present in the protein. In such cases, further approaches and data are very useful to enhance the information content. Relaxation filtered hyperfine spectroscopy (REFINE) can be used to separate the individual components of overlapping paramagnetic species on the basis of differences in their longitudinal relaxation rates; it is applicable to any kind of pulsed hyperfine or dipolar spectroscopy. Here, we show that the spectra of the iron-sulfur clusters in complex I can be separated by this method, allowing us to obtain hyperfine (and dipolar) information from the individual species. Furthermore, performing pulsed EPR experiments at different magnetic fields is another important tool to disentangle the spectral components in such complex systems. Despite the fact that high magnetic fields do not usually lead to better spectral separation for metal centers, they provide additional information about the relative orientation of different paramagnetic centers. Our high-field EPR studies on cytochrome c oxidase reveal essential information regarding the structural arrangement of the binuclear Cu(A) center with respect to both the manganese ion within the enzyme and the cytochrome in the protein-protein complex with cytochrome c.

Screening of genes related to ovary development in Chinese shrimp Fenneropenaeus chinensis by suppression subtractive hybridization

The ovary of triploid shrimp Fenneropenaeus chinensis was apparently impaired compared to that of the diploid shrimp at the same age. Therefore triploid shrimp ovary is possible to be taken as a model to understand the mechanism of ovary development of shrimp compared to that of the ovary of diploid shrimp at the same age. In the present study, a suppression subtractive hybridization (SSH) technique was applied to identify differentially expressed genes in the ovary between diploid and triploid shrimp. For the forward library (RNA from the ovary of triploid shrimp as the tester), 54 genes were identified. For the reverse library (RNA from the ovary of diploid shrimp as the tester), 16 genes were identified. The identified genes encoded proteins with multiple functions, including extracellular matrix components, cytoskeleton, cell growth and death, metabolism, genetic information processing, signal transduction/transport or immunity related proteins. Eleven differentially expressed genes were selected to be confirmed in the ovaries of triploid and diploid shrimp by semi-quantitative RT-PCR. Genes encoding spermatogonial stem-cell renewal factor, cytochrome c oxidase subunits I and II, clottable protein, antimicrobial peptide and transposase showed up-regulated expressions in the ovary of triploid shrimp. Genes encoding tubulin, cellular apoptosis susceptibility protein, farnesoic acid O-methyltransferase, thrombospondin and heat shock protein 90 genes showed higher expressions in the ovary of diploid shrimp. The differential expressions of the above genes are suggested to be related to the ovary development of shrimp. It will provide a new clue to uncover the molecular mechanisms underlying the ovarian development in penaeid shrimp.

Low temperature 65Cu NMR spectroscopy of the Cu+ site in azurin

(65)Cu central-transition NMR spectroscopy of the blue copper protein azurin in the reduced Cu(I) state, conducted at 18.8 T and 10 K, gave a strongly second order quadrupole perturbed spectrum, which yielded a (65)Cu quadrupole coupling constant of +/-71.2 +/- 1 MHz, corresponding to an electric field gradient of +/-1.49 atomic units at the copper site, and an asymmetry parameter of approximately 0.2. Quantum chemical calculations employing second order Møller-Plesset perturbation theory and large basis sets successfully reproduced these experimental results. Sensitivity and relaxation times were quite favorable, suggesting that NMR may be a useful probe of the electronic state of copper sites in proteins.

Identification of Escherichia coli HemG as a novel, menadione-dependent flavodoxin with protoporphyrinogen oxidase activity

Protoporphyrinogen oxidase (PPO, EC 1.3.3.4) catalyzes the six-electron oxidation of protoporphyrinogen IX to the fully conjugated protoporphyrin IX. Eukaryotes and Gram-positive bacteria possess an oxygen-dependent, FAD-containing enzyme for this step, while the majority of Gram-negative bacteria lack this oxygen-dependent PPO. In Escherichia coli, PPO activity is known to be linked to respiration and the quinone pool. In E. coli SASX38, the knockout of hemG causes a loss of measurable PPO activity. HemG is a small soluble protein typical of long chain flavodoxins. Herein, purified recombinant HemG was shown to be capable of a menadione-dependent conversion of protoporphyrinogen IX to protoporphyrin IX. Electrochemical analysis of HemG revealed similarities to other flavodoxins. Interestingly, HemG, a member of a class of the long chain flavodoxin family that is unique to the gamma-proteobacteria, possesses a 22-residue sequence that, when transferred into E. coli flavodoxin A, produces a chimera that will complement an E. coli hemG mutant, indicating that this region confers PPO activity to the flavodoxin. These findings reveal a previously unidentified class of PPO enzymes that do not utilize oxygen as an electron acceptor, thereby allowing gamma-proteobacteria to synthesize heme in both aerobic and anaerobic environments.

Dielectric relaxation of cytochrome c oxidase: Comparison of the microscopic and continuum models

We have studied a charge-insertion process that models the deprotonation of a histidine side chain in the active site of cytochrome c oxidase (CcO) using both the continuum electrostatic calculations and the microscopic simulations. The group of interest is a ligand to Cu(B) center of CcO, which has been previously suggested to play the role of the proton pumping element in the enzyme; the group is located near a large internal water cavity in the protein. Using the nonpolarizable Amber-99 force field in molecular dynamics (MD) simulations, we have calculated the nuclear part of the reaction-field energy of charging of the His group and combined it with the electronic part, which we estimated in terms of the electronic continuum (EC) model, to obtain the total reaction-field energy of charging. The total free energy obtained in this MDEC approach was then compared with that calculated using pure continuum electrostatic model with variable dielectric parameters. The dielectric constant for the "dry" protein and that of the internal water cavity of CcO were determined as those parameters that provide best agreement between the continuum and microscopic MDEC model. The nuclear (MD) polarization alone (without electronic part) of a dry protein was found to correspond to an unphysically low dielectric constant of only about 1.3, whereas the inclusion of electronic polarizability increases the protein dielectric constant to 2.6-2.8. A detailed analysis is presented as to how the protein structure should be selected for the continuum calculations, as well as which probe and atomic radii should be used for cavity definition. The dielectric constant of the internal water cavity was found to be 80 or even higher using "standard" parameters of water probe radius, 1.4 A, and protein atomic radii from the MD force field for cavity description; such high values are ascribed to the fact that the standard procedure produces unphysically small cavities. Using x-ray data for internal water in CcO, we have explored optimization of the parameters and the algorithm of cavity description. For Amber radii, the optimal probe size was found to be 1.25 A; the dielectric of water cavity in this case is in the range of 10-16. The most satisfactory cavity description, however, was achieved with ProtOr atomic radii, while keeping the probe radius to be standard 1.4 A. In this case, the value of cavity dielectric constant was found to be in the range of 3-6. The obtained results are discussed in the context of recent calculations and experimental measurements of dielectric properties of proteins.

Molecular dynamics simulation of water in cytochrome c oxidase reveals two water exit pathways and the mechanism of transport

We have examined the network of connected internal cavities in cytochrome c oxidase along which water produced at the catalytic center is removed from the enzyme. Using combination of structural analysis, molecular dynamics simulations, and free energy calculations we have identified two exit pathways that connect the Mg2+ ion cavity to the outside of the enzyme. Each pathway has a well-defined bottleneck, which determines the overall rate of water traffic along the exit pathway, and a specific cooperative mechanism of passing it. One of the pathways is going via Arg438/439 (in bovine numbering) toward the CuA center, approaching closely its His204B ligand and Lys171B residue; and the other is going toward Asp364 and Thr294. Comparison of the pathways among different aa3-type enzymes shows that they are well conserved. Possible connections of the finding to redox-coupled proton pumping mechanism are discussed. We propose specific mutations near the bottlenecks of the exit pathways that can test some of our hypotheses.

High resolution crystal structure of Paracoccus denitrificans cytochrome c oxidase: new insights into the active site and the proton transfer pathways

The structure of the two-subunit cytochrome c oxidase from Paracoccus denitrificans has been refined using X-ray cryodata to 2.25 A resolution in order to gain further insights into its mechanism of action. The refined structural model shows a number of new features including many additional solvent and detergent molecules. The electron density bridging the heme a(3) iron and Cu(B) of the active site is fitted best by a peroxo-group or a chloride ion. Two waters or OH(-) groups do not fit, one water (or OH(-)) does not provide sufficient electron density. The analysis of crystals of cytochrome c oxidase isolated in the presence of bromide instead of chloride appears to exclude chloride as the bridging ligand. In the D-pathway a hydrogen bonded chain of six water molecules connects Asn131 and Glu278, but the access for protons to this water chain is blocked by Asn113, Asn131 and Asn199. The K-pathway contains two firmly bound water molecules, an additional water chain seems to form its entrance. Above the hemes a cluster of 13 water molecules is observed which potentially form multiple exit pathways for pumped protons. The hydrogen bond pattern excludes that the Cu(B) ligand His326 is present in the imidazolate form.

Heme-copper oxidases and their electron donors in cyanobacterial respiratory electron transport

Cyanobacteria are the paradigmatic organisms of oxygenic (plant-type) photosynthesis and aerobic respiration. Since there is still an amazing lack of knowledge on the role and mechanism of their respiratory electron transport, we have critically analyzed all fully or partially sequenced genomes for heme-copper oxidases and their (putative) electron donors cytochrome c(6), plastocyanin, and cytochrome c(M). Well-known structure-function relationships of the two branches of heme-copper oxidases, namely cytochrome c (aa(3)-type) oxidase (COX) and quinol (bo-type) oxidase (QOX), formed the base for a critical inspection of genes and ORFs found in cyanobacterial genomes. It is demonstrated that at least one operon encoding subunits I-III of COX is found in all cyanobacteria, whereas many non-N(2)-fixing species lack QOX. Sequence analysis suggests that both cyanobacterial terminal oxidases should be capable of both the four-electron reduction of dioxygen and proton pumping. All diazotrophic organisms have at least one operon that encodes QOX. In addition, the highly refined specialization in heterocyst forming Nostocales is reflected by the presence of two paralogs encoding COX. The majority of cyanobacterial genomes contain one gene or ORF for plastocyanin and cytochrome c(M), whereas 1-4 paralogs for cytochrome c(6) were found. These findings are discussed with respect to published data about the role of respiration in wild-type and mutated cyanobacterial strains in normal metabolism, stress adaptation, and nitrogen fixation. A model of the branched electron-transport pathways downstream of plastoquinol in cyanobacteria is presented.

Acidic lipids, H(+)-ATPases, and mechanism of oxidative phosphorylation. Physico-chemical ideas 30 years after P. Mitchell's Nobel Prize award

Peter D. Mitchell, who was awarded the Nobel Prize in Chemistry 30 years ago, in 1978, formulated the chemiosmotic theory of oxidative phosphorylation. This review initially analyzes the major aspects of this theory, its unresolved problems, and its modifications. A new physico-chemical mechanism of energy transformation and coupling of oxidation and phosphorylation is then suggested based on recent concepts regarding proteins, including ATPases that work as molecular motors, and acidic lipids that act as hydrogen ion (H(+)) carriers. According to this proposed mechanism, the chemical energy of a redox substrate is transformed into nonequilibrium states of electron-transporting chain (ETC) coupling proteins. This leads to nonequilibrium pumping of H(+) into the membrane. An acidic lipid, cardiolipin, binds with this H(+) and carries it to the ATP-synthase along the membrane surface. This transport generates gradients of surface tension or electric field along the membrane surface. Hydrodynamic effects on a nanolevel lead to rotation of ATP-synthase and finally to the release of ATP into aqueous solution. This model also explains the generation of a transmembrane protonmotive force that is used for regulation of transmembrane transport, but is not necessary for the coupling of electron transport and ATP synthesis.

Toward a chemical mechanism of proton pumping by the B-type cytochrome c oxidases: application of density functional theory to cytochrome ba3 of Thermus thermophilus

A mechanism for proton pumping by the B-type cytochrome c oxidases is presented in which one proton is pumped in conjunction with the weakly exergonic, two-electron reduction of Fe-bound O 2 to the Fe-Cu bridging peroxodianion and three protons are pumped in conjunction with the highly exergonic, two-electron reduction of Fe(III)- (-)O-O (-)-Cu(II) to form water and the active oxidized enzyme, Fe(III)- (-)OH,Cu(II). The scheme is based on the active-site structure of cytochrome ba 3 from Thermus thermophilus, which is considered to be both necessary and sufficient for coupled O 2 reduction and proton pumping when appropriate gates are in place (not included in the model). Fourteen detailed structures obtained from density functional theory (DFT) geometry optimization are presented that are reasonably thought to occur during the four-electron reduction of O 2. Each proton-pumping step takes place when a proton resides on the imidazole ring of I-His376 and the large active-site cluster has a net charge of +1 due to an uncompensated, positive charge formally associated with Cu B. Four types of DFT were applied to determine the energy of each intermediate, and standard thermochemical approaches were used to obtain the reaction free energies for each step in the catalytic cycle. This application of DFT generally conforms with previously suggested criteria for a valid model (Siegbahn, P. E. M.; Blomberg, M. A. R. Chem. Rev. 2000, 100, 421-437) and shows how the chemistry of O 2 reduction in the heme a 3 -Cu B dinuclear center can be harnessed to generate an electrochemical proton gradient across the lipid bilayer.

Age-related alterations in oxidatively damaged proteins of mouse skeletal muscle mitochondrial electron transport chain complexes

Age-associated mitochondrial dysfunction is a major source of reactive oxygen species (ROS) and oxidative modification to proteins. Mitochondrial electron transport chain (ETC) complexes I and III are the sites of ROS production and we hypothesize that proteins of the ETC complexes are primary targets of ROS-mediated modification which impairs their structure and function. The pectoralis, primarily an aerobic red muscle, and quadriceps, primarily an anaerobic white muscle, have different rates of respiration and oxygen-carrying capacity, and hence, different rates of ROS production. This raises the question of whether these muscles exhibit different levels of oxidative protein modification. Our studies reveal that the pectoralis shows a dramatic age-related decline in almost all complex activities that correlates with increased oxidative modification. Similar complex proteins were modified in the quadriceps, at a significantly lower level with less change in enzyme and ETC coupling function. We postulate that mitochondrial ROS causes damage to specific ETC subunits which increases with age and leads to further mitochondrial dysfunction. We conclude that physiological characteristics of the pectoralis vs quadriceps may play a role in age-associated rate of mitochondrial dysfunction and in the decline in tissue function.

Expressed sequence tags of Trichinella spiralis muscle stage larvae

In order to obtain greater insight into the relevant genomic expression patterns of Trichinella spiralis, 992 expressed sequence tags (ESTs) were collected from a cDNA library of T. spiralis muscle stage larvae and assembled into 60 clusters and 385 singletons. Of them, 445 (44.7%) ESTs were annotated to their homologous genes, and small fractions were matched to known genes of nematodes. The annotated ESTs were classified into 25 eukaryotic orthologous groups (KOG). Cytochrome C oxidase (34 clones) was found to be most frequent species.

Thermodynamic redox behavior of the heme centers in A-type heme-copper oxygen reductases: comparison between the two subfamilies

The study of the thermodynamic redox behavior of the hemes from two members of the A family of heme-copper oxygen reductases, Paracoccus denitrificans aa3 (A1 subfamily) and Rhodothermus marinus caa3 (A2 subfamily) enzymes, is presented. At different pH values, midpoint reduction potentials and interaction potentials were obtained in the framework of a pairwise model for two interacting redox centers. In both enzymes, the hemes have different reduction potentials. For the A1-type enzyme, it was shown that heme a has a pH-dependent midpoint reduction potential, whereas that of heme a3 is pH independent. For the A2-type enzyme the opposite was observed. The midpoint reduction potential of heme c from subunit II of the caa3 enzyme was determined by fitting the data with a single-electron Nernst curve, and it was shown to be pH dependent. The results presented here for these A-type enzymes are compared with those previously obtained for representative members of the B and C families.

Electrostatic basis for the unidirectionality of the primary proton transfer in cytochrome c oxidase

Gaining detailed understanding of the energetics of the proton-pumping process in cytochrome c oxidase (CcO) is one of the challenges of modern biophysics. Despite promising mechanistic proposals, most works have not related the activation barriers of the different assumed steps to the protein structure, and there has not been a physically consistent model that reproduced the barriers needed to create a working pump. This work reevaluates the activation barriers for the primary proton transfer (PT) steps by calculations that reflect all relevant free energy contributions, including the electrostatic energies of the generated charges, the energies of water insertion, and large structural rearrangements of the donor and acceptor. The calculations have reproduced barriers that account for the directionality and sequence of events in the primary PT in CcO. It has also been found that the PT from Glu-286 (E) to the propionate of heme a(3) (Prd) provides a gate for an initial back leakage from the high pH side of the membrane. Interestingly, the rotation of E that brings it closer to Prd appears to provide a way for blocking competing pathways in the primary PT. Our study elucidates and quantifies the nature of the control of the directionality in the primary PT in CcO and provides instructive insight into the role of the water molecules in biological PT, showing that "bridges" of several water molecules in hydrophobic regions present a problem (rather than a solution) that is minimized in the primary PT.

Looking for the minimum common denominator in haem-copper oxygen reductases: towards a unified catalytic mechanism

Haem-copper oxygen reductases are transmembrane protein complexes that reduce dioxygen to water and pump protons across the mitochondrial or periplasmatic membrane, contributing to the transmembrane difference of electrochemical potential. Seven years ago we proposed a classification of these enzymes into three different families (A, B and C), based on the amino acid residues of their proton channels and amino acid sequence comparison, later supported by the so far identified characteristics of the catalytic centre of members from each family. The three families have in common the same general structural fold of the catalytic subunit, which contains the same or analogous prosthetic groups, and proton channels. These observations raise the hypothesis that the mechanisms for dioxygen reduction, proton pumping and the coupling of the two processes may be the same for all these enzymes. Under this hypothesis, they should be performed and controlled by the same or equivalent elements/events, and the identification of retained elements in all families will reveal their importance and may prompt the definition of the enzyme operating mode. Thus, we believe that the search for a minimum common denominator has a crucial importance, and in this article we highlight what is already established for the haem-copper oxygen reductases and emphasize the main questions still unanswered in a comprehensive basis.